<figure><img src="image_1.png" /></figure>

HAL
open science

# Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies.

Chitra Narayanan, David N Bernard, Nicolas Doucet

► To cite this version:

Chitra Narayanan, David N Bernard, Nicolas Doucet. Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies. Catalysts, 2017, Special Issue "Enzymes Catalysts", 6 (6), pp.81. 10.3390/catal6060081. pasteur-01534404

HAL Id: pasteur-01534404

https://riip.hal.science/pasteur-01534404

Submitted on 7 Jun 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies

Chitra Narayanan¹,†, David N. Bernard¹,†, and Nicolas Doucet¹,²,³,*

¹INRS—Institut Armand-Frappier, Université du Québec, 531 Boul. des Prairies, Laval, QC H7V 1B7, Canada (C.N.)

²PROTEO, the Québec Network for Research on Protein Function, Engineering, and Applications, 1045 Avenue de la Médecine, Université Laval, Québec, QC G1V 0A6, Canada

³GRASP, the Groupe de Recherche Axé sur la Structure des Protéines, 3649 Promenade Sir William Osler, McGill University, Montréal, QC H3G 0B1, Canada

## Abstract

It is now common knowledge that enzymes are mobile entities relying on complex atomic-scale dynamics and coordinated conformational events for proper ligand recognition and catalysis. However, the exact role of protein dynamics in enzyme function remains either poorly understood or difficult to interpret. This mini-review intends to reconcile biophysical observations and biological significance by first describing a number of common experimental and computational methodologies employed to characterize atomic-scale residue motions on various timescales in enzymes, and second by illustrating how the knowledge of these motions can be used to describe the functional behavior of enzymes and even act upon it. Two biologically relevant examples will be highlighted, namely the HIV-1 protease and DNA polymerase β enzyme systems.

## Keywords

enzyme catalysis; conformational dynamics; nuclear magnetic resonance; molecular dynamics simulations; HIV-1 protease; DNA polymerase β

## 1. Introduction

Proteins and enzymes are essential components of living cells. Among a variety of other functions, they act as hormones in cell signaling, protein transporters, antibodies in host defense, etc. Concomitantly, enzymes are some of the most efficient catalysts known to date, enhancing biochemical reaction rates up to 19 orders of magnitude relative to uncatalyzed

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).

**Author Contributions:** All authors made substantial contributions to the conception, design, and writing of this manuscript. All authors have read and approved the final version of the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.

<footnote>*Correspondence: nicolas.doucet@iaf.inrs.ca; Tel.: +450-687-5010 (ext. 4212).</footnote>

<footnote>†These authors contributed equally to this work.</footnote>

<footnote>Academic Editor: David D. Boehr</footnote>
reactions [1,2]. However, factors contributing to these large rate enhancements in enzyme-catalyzed reactions remain largely uncharacterized. The original structure-function paradigm, popular for many decades, portrayed enzymes as rigid structures with shapes that facilitate substrates, ligand and/or cofactor binding. This original model evolved over time, and theories were developed that led to the now-accepted induced fit, conformational selection, and transition-state stabilization models to explain the behavior of protein-ligand recognition and catalysis in the molecular function of enzymes [3]. Increasing evidence suggests that proteins sample a variety of distinct conformations (or sub-states) enabled by concerted atomic-scale dynamical fluctuations occurring over a wide range of timescales and acting on the primary, secondary, tertiary, and quaternary organization of their molecular structure [4]. Conformational transitions between highly and rarely populated states have been shown to play important roles in substrate recognition, binding, and product release, among others [5–8]. Advances in experimental and computational methodologies continue to offer new insights into enzyme conformational motions over functionally relevant timescales [9–12]. Interestingly, experimental [13,14], computational [15,16], and sequence-based [17,18] approaches have also revealed functional networks of concerted residue motions distant from the active site in selected enzyme systems. Correlations between the timescale of conformational fluctuations and that of catalytic turnover have been well established in a variety of enzyme systems, including cyclophilin A [14], RNase A [5,19,20], triosephosphate isomerase [21], and HIV-1 protease [22], among others. Further, the rate of conformational exchange has been shown to coincide with the rate-limiting step, such as product release, in some of these systems [19,20]. Conformational exchange between sub-states allows enzymes to sample higher energy conformations with structural and dynamical properties important for function such as ligand binding and allosteric regulation [11,16,23]. In addition to the millisecond conformational exchange, dynamics on faster timescales modulate the chemical environment through rearrangements in the active site, thus affecting enzyme function [24]. Taken together, a view is emerging whereby conformational fluctuations occurring over a range of timescales can affect enzyme function through conformational sampling along preferred pathways.

Over the past few years, a number of controversial statements relating ‘catalysis’ with ‘dynamics’ have been published in the literature, often giving rise to heated debates between experimentalists and theoreticians [25,26]. As recently outlined, these debates are often semantic in nature and can be traced back to actual definitions (and research field perceptions) of what ‘protein dynamics’ represents, in addition to which atomic-scale events are being observed and/or over which timescales they occur during enzyme-catalyzed reactions [25]. In the current report, we consider motions affecting enzyme function in a broad sense, i.e., “any motion that, if impeded, would reduce the ability of the enzyme to function by the mechanism that it has evolved to execute” [27]. As such, we take into account motions on a wide range of timescales, from the fast, local picosecond-nanosecond (ps-ns) fluctuations of individual residues to the slower microsecond-millisecond-second (µs-ms-s) global conformational exchange experienced by enzymes. Conserved motional events occurring on a broad range of timescales could be involved in many and/or a selected subset of events occurring along the catalytic cycle of an enzyme, including ligand binding, recognition, substrate discrimination, structural rearrangements leading to the transition
state, product release, etc. [25]. To clarify and prevent confusion, we never imply that
'motions contributing to catalysis' exclusively describe the fast femtosecond atomic motions
involved in transition-state chemistry.

In the current account, we draw attention to a number of experimental and computational
methodologies that have recently improved our understanding of catalytic and functional
properties in two enzyme systems that closely rely on conformational dynamics for proper
biological function. We do not pretend to cover the overwhelming number of theoretical
approaches, methodologies, or enzymatic systems that previously illustrated the role of
protein dynamics in enzyme function, but instead focus on recent reports where the
combination of experimental NMR and computational techniques have emphasized the role
of conformational exchange in inhibitor binding (HIV-1 protease) and fidelity in DNA repair
through the selection of correct nucleotides (DNA polymerase β).

## 2. Selected Experimental and Computational Approaches for Sampling Conformational Motions in Proteins

A variety of approaches, including (but not limited to) NMR spectroscopy, X-ray diffraction, single molecule FRET, and computational simulations, have been used to probe conformational dynamics in proteins over nanosecond to millisecond timescales, offering the means to extract motions potentially relevant to biological function [4,28]. One of the most commonly used experiments to characterize local and global conformational events experienced by enzymes on the timescale of their catalytic rate kcat—typically occurring on the order of milliseconds—is the NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiment [29–33]. This methodology garnered considerable attention in recent years since motions experienced by enzymes on this physiologically relevant timescale can improve our understanding of catalysis and allostery, as the flexibility events uncovered can be analyzed in light of catalytic events. The CPMG method enables the detection of atomic-scale millisecond dynamics of enzymes in solution, either in the free form or bound to biologically significant ligands/analogues that mimic intermediate states along the chemical reaction [13,34,35]. On this time frame, molecular flexibility may involve the rearrangement of single residues and/or entire secondary structures moving in a coordinated fashion (for reviews see [9,36–38]) (Figure 1). In its typical form, the method involves the characterization of an equilibrium exchange process experienced by the N–H amide bond of amino acid residues, performed by recording the transverse relaxation rate constant (R2) as a function of vCPMG, the frequency with which refocusing pulses are applied in a CPMG pulse train [30,39]. These refocusing pulses are effectively invisible to amino acid residues that do not undergo conformational exchange in the protein, but they add an additional contribution (Rx) to R2 for exchanging residues—those experiencing motions on this millisecond time frame—and thus causing a decrease in NMR resonance intensity. Plotting this decrease in signal intensity as a function of vCPMG not only allows the qualitative identification of residues experiencing conformational exchange in the enzyme, but fitting the curves to relaxation dispersion equations also allows the quantitative extraction of structural information related to the excited (less populated) sub-states (differences in chemical shifts,
$\Delta\omega$), in addition to their population dynamics in solution (conformational exchange rates, $k_{ex}$, and equilibrium populations, $p_A$ and $p_B$) (see [39] and references therein).

With regards to the timescale of protein motions, NMR is not limited to the µs-ms range. To this day, the model-free formalism developed by Lipari and Szabo in the early 1980s remains one of the preferred approaches to determine the atomic-scale dynamics of a protein in solution [41,42]. Using spin relaxation NMR data, the method primarily provides information on residue dynamics occurring on the ps-ns timescale, which are most often attributed to single residue side-chain rotations or local residue motions. The method measures spin relaxation rates for each residue amide N-H pair to extract order parameters ($S^2$), conformational exchange parameters ($R_{ex}$), and the correlation time ($\tau_c$), offering an interpretation of the internal amplitude of the N-H bond vector in solution [43–45]. The order parameter $S^2$ provides an easily interpretable measure of the rigidity of each protein residue, where the $S^2$ of a completely unrestricted N-H bond would theoretically correspond to a value of 0, while that of a fully rigid residue would be 1. For a number of residues that do not fit to simple calculation models, an $R_{ex}$ exchange parameter can be estimated to achieve an adequate fit to the model-free analysis [45].

While the aforementioned NMR relaxation analyses provide detailed information on the molecular flexibility experienced by enzymes in solution, other experimental methodologies offer additional details pertaining to potentially significant biological motions. Single-molecule Förster Resonance Energy Transfer (single-molecule FRET) [4,46] and site-directed spin labeling Double Electron-Electron Resonance (DEER) spectroscopy [47]—also called pulsed Electron Paramagnetic Resonance (pulsed-EPR)—are both used to determine the distance between two labels, ranging up to a few tens of Å. In DEER, a pair of cysteine residues is chemically modified with electronic spin labels, namely the nitroxide radical of a methylthiosulfonate moiety [22,48]. The magnitude of the dipolar coupling between the two spins can be extracted from the measured DEER echo curves, and this magnitude is inversely proportional to the cube of the distance between the two spins [49], thus producing a distance profile. Populations of the different states of the protein can then be evaluated by deconvoluting the distance profile. Much like with protein crystallography, a disadvantage of the technique is the requirement of cryogenic temperatures [47], which does not provide an accurate description of biological conditions. On the other hand, single-molecule FRET can be used in solution and at physiologically relevant temperatures, much like NMR. This technique also requires labeling of two distinct sites on a protein, this time using two different fluorophores. Since these fluorophores are selected so that the absorption spectrum of the first (the acceptor) matches the emission spectrum of the second (the donor), the fluorescence energy can be transmitted from the donor to the acceptor if they are in close proximity to each other. Thus, the emission wavelength of the donor will disappear and the emission wavelength of the acceptor will be detected as a function of distance. By measuring the FRET efficiency, the inter-fluorophore distance can be calculated. The single-molecule version of this experiment can thus allow one to witness whole domain reorientations in a protein. This technique also has the advantage of being performed in real-time [4,46].

While most experimental methodologies yield quantitative information on flexibility events, they often fail to provide a clear visual representation of the actual atomic-scale motions
experienced by an enzyme in solution [44]. In contrast, computational molecular dynamics (MD) simulations trace the position and interaction of atoms over time based on classical mechanics principles, to provide an atomic-scale time evolution of conformational properties in proteins and other biomolecules. Until recently, MD simulations typically sampled over hundreds of nanoseconds, while protein motions span a much wider range of timescales [4], with many catalytically relevant motions occurring on the same timescale as the catalytic constant $k_{cat}$, i.e., typically milliseconds or slower [19,25,50,51]. This presents an existing limitation of computational MD simulations. Advances in software and hardware capabilities have facilitated access to longer timescales, with several recent studies reporting ms timescale simulations [52–54]. While a single ms trajectory has been achieved using Anton [55], a specialized supercomputer for protein simulations, access to the ms timescale with traditional, unbiased MD simulations still remains limited.

A variety of enhanced sampling techniques such as accelerated MD (aMD) [56,57], umbrella sampling [58], steered MD [59], and metadynamics [60–62], among others, have been developed to improve sampling of the conformational landscape of proteins. In aMD, the potential energy landscape is modified by raising energy minima that are below a threshold, therefore minimizing the energy barriers separating states, and facilitating access to conformational states not easily accessible through traditional MD [57]. This approach has been used to characterize dynamics on the µs-ms timescale in a variety of protein systems [63]. A combination of normal and accelerated MD simulations was recently used to characterize the role of conformational dynamics in stabilizing the transition state in the enzyme cyclophilin A [64]. Umbrella sampling uses a biasing potential to characterize the higher energy states along a reaction coordinate [65]. Metadynamics involves the use of a history-dependent biasing potential that acts on collective variables corresponding to select degrees of freedom [60]. A combination of NMR and metadynamics was recently used to characterize the conformational ensemble and the free energy landscape of the highly dynamic helix 1 of the prion protein [66]. While these approaches have clear benefits from enhanced sampling, the requirement of the knowledge of start and end structures (umbrella sampling), or the selection of suitable reaction coordinates (metadynamics), among others, limit their application [67].

Other advanced sampling methods, such as the Markov State Model (MSM, also known as the kinetic network model) [68–70] and Quasi Anharmonic Analysis (QAA) [11,71], facilitate the identification of conformational sub-states, equilibrium populations of sub-states, and transition probabilities between these sub-states. MSMs are constructed from trajectories obtained from MD simulations. In these models, the conformational space is discretized into multiple states, and a network connecting these states is created from the MD simulations. A transition matrix corresponding to the probability of transitions between different states is constructed to characterize the kinetics of transition between states. MSMs have been applied to sample a wide range of timescales associated with protein folding, functional dynamics, ligand binding and for characterizing intrinsically disordered proteins [72–74]. QAA uses higher order statistics to describe positional fluctuations and the coupling of these fluctuations in different regions of the protein [71]. It allows the characterization of conformational sub-states along a reaction pathway and the identification of sub-states with structural and dynamical properties important for function [75]. This
approach was used to identify the hierarchical organization of conformational sub-states in
ubiquitin and the T4 lysozyme, and the conformational sub-states associated with the cis/
trans isomerization catalyzed by cyclophilin A were also characterized [75]. This technique
was recently used to characterize conformational sub-states associated with the
interconversion between reactant and product states in wild-type and mutant forms of RNase
A [76]. Recent years have also seen advances in the integration of experimental and
computational approaches for characterizing the dynamics of proteins. NMR chemical shifts
and residual dipolar couplings have been used as replica-averaged structural restraints to
characterize the conformational fluctuations in RNase A and the hen lysozyme, respectively
[77–79].

3. Case Studies

**3.1. HIV Protease**

The protease of the type 1 human immunodeficiency virus (HIV-1 PR) is an aspartyl
proteinase that cleaves the polyprotein product of the *gag-pro-pol* gene and the Gag protein
itself, both of which are vital steps in the maturation of the virus [80]. The three most
prevalent subtypes of the virus are those whose proteases have been the most studied:
subtypes B, C, and CRF_01 A/E, which are the prevalent HIV subtypes in North America
and Europe, Africa, and Asia, respectively. Even though all three subtypes have a distinct
wild-type (WT) protease, subtype B is the most widely characterized [47]. To date, nine
Food and Drug Administration (FDA)-approved inhibitors targeting this protein have been
used to treat HIV/acquired immune deficiency syndrome (AIDS)-infected patients.
However, multidrug resistance to protease inhibitors has sparked renewed interests in
uncovering distinct inhibitors targeting HIV-1 PR, and thus understanding drug resistance
acquisition mechanisms is of utmost importance [81].

HIV-1 PR is only active in its homodimeric form. The 99-amino-acid protomers assemble in a symmetric configuration, forming a tunnel. The protomers each provide one of the catalytic aspartate residues (Asp25, located in the tunnel), while access to the catalytic site is allowed through the opening and closing of a pair of β-hairpins, termed “flaps” (residues 43–58) (colored red in Figure 2). The protomers interact with each other through their N- and C-termini, as well as at the tip of the flaps [82]. A previous report briefly discusses the atomic-scale dynamics experienced by WT HIV-1 PR, describing how movements of the flaps originating from the isomerization of Gly51 are transmitted to the active site [83]. This vital conformational transition, identified as early as 1995, involves the 180° reorientation of the Gly51 backbone N–H vector, exchanging between β-turn types I and II, a process taking place over a 10 µs timescale, as determined from R<sub>ex</sub> values extracted using the Lipari-Szabo model-free formalism mentioned above [84]. This unique residue reorganization represents the sole significant structural difference between otherwise symmetrical protomers in the crystal structures. Most commonly, Gly51 is found in the *L*-conformation in one protomer and in the *D*-conformation in the opposing protomer [22]. Model-free analyses also confirmed that the flap region, along with the hinge loop (residues 34–42), are the primary regions in which sub-nanosecond timescale residue dynamics are localized [85].
NMR relaxation and computational MD simulations demonstrated that flaps can adopt a variety or conformations as the enzyme proceeds through its catalytic cycle [47,86,87]. These flaps alternate between a “closed” conformation, a “semi-open” conformation, a “wide open” conformation, and a “curled/tucked” conformation (Figure 2). The flaps close over the active site to maintain the substrate inside the catalytic pocket and to provide a gating mechanism for substrate binding, but they also interact with each other in the absence of substrate or inhibitor [85]. A combination of DEER measurements, NMR relaxation experiments, and MD simulations provided information on Asp25 and flap dynamics in each protomer during catalytic turnover. The µs-ms timescale motions of Asp25 and Gly27 were found to correlate with those of the flaps [22]. It was also observed that the more flexible L-Gly51 must be located on the same protomer as the protonated general acid Asp25-COOH for efficient catalysis to occur, while the more rigid D-Gly51 must be located on the protomer bearing the unprotonated general base Asp25-COO-. The rate-limiting step of catalytic turnover was shown to be linked to isomerization of the flap tips towards the type I and type II β-turn configurations [22], further illustrating the functional importance of conformational exchange in this enzyme system.

The kinetic network model, constructed by discretizing the conformational space of atomistic simulations into a series of related sub-states, was used to determine the transitions between the semi-open, closed, and open conformations of HIV-1 PR [88]. This approach facilitates the characterization of interconverting rates between kinetically relevant sub-states, and can predict kinetic quantities on slower (millisecond) timescales [72,89]. Deng *et al.* combined replica exchange molecular dynamics simulations with transition path theory to characterize the diversity and temperature-dependence of kinetic pathways between the three conformational states of HIV-1 PR. They showed that transitions between semi-open and closed states occur on a fast timescale (~33 ns), while transitions between semi-open and open states were infrequent (~375 ns), consistent with NMR observations showing motions on the fast sub-nanosecond and slower microsecond timescales [88].

Not all studies agree on which conformation is the most likely populated form of the apo WT protease in solution. Studies using the DEER method concluded that the semi-open conformation is the most populated state, independent of the subtype [22,47], while NMR residual dipolar coupling (RDC) measurements suggested that the closed conformation was predominant [90]. However, all studies agree that the populations of these conformers change upon inhibitor binding, with the most potent inhibitors favoring the “closed” conformation [47,91–93]. In all common subtypes of WT HIV-1 PR (subtypes B, C, and CRF_01 A/E), the residence time of inhibitors in the active-site pocket, and thus the inhibitory strength, is correlated with the degree to which said inhibitor induces flap closure. Consequently, many drug resistance-inducing mutations perturb the dynamic behavior of the flaps relative to the WT enzyme [47]. Different inhibitors cause varying population shifts among the four possible flap conformations. According to DEER measurements, while inhibitors indinavir, nelfinavir, and atazanavir trigger very small differences in populations in all three subtypes—the semi-open conformation being the most populated—all other six FDA-approved inhibitors strongly shift the population equilibrium towards the closed configuration. Interestingly, the same six inhibitors also make the curled/tucked conformer much more likely in subtype CRF_01 A/E [47]. While smaller in scope than those of the
flaps, conformational changes in the hydrophobic core predicted by MD simulations were also confirmed to impact protease activity [94]. In the WT HIV-1 PR, the two halves of the hydrophobic core slide over each other during these rearrangements, and cross-linking these two halves—thus impeding hydrophobic core rearrangements—causes a 150-fold reduction in catalytic efficiency [94,95]. Also, mutations distal to the active site in the hydrophobic core are hypothesized to cause drug resistance by altering the flexibility of the protein [94].

Previous reports have highlighted that many distinct mutations can arise in patients developing drug resistance to HIV treatment [96]. Some of these mutations are located in the active-site pocket and directly impair inhibitor binding (primary mutations), while others are more distal and indirectly induce resistance (secondary mutations) (Figure 3). None of these distal mutations cause significant changes in either $k_{cat}$ or $K_M$ in HIV-1 PR, suggesting that conformational events experienced by the enzyme play a significant role in resistance acquisition, which was confirmed by NMR, MD simulations, and DEER spectroscopy. Mutations in 45 of the 99 residues in HIV-1 PR have been linked to drug resistance *in vitro* [97], and a quarter of the residues have been linked to drug resistance in patients, arising through evolutionary selective pressure driven by protease inhibitor therapy. However, mutations at only seven sites can be linked to direct active-site modifications (Figure 3a) [92,97,98]. These primary mutations often decrease ligand binding affinity, both for inhibitors and natural substrates. Secondary mutations often arise to restore activity towards substrates, sometimes as a result of conformational behavior. For instance, incorporation of I15V, E35D, R41K, and R57K was demonstrated by DEER to increase backbone flexibility in WT HIV-1 PR [48]. These mutations, which are all believed to act as function-restoring secondary replacements, result in the increased population of the curled/tucked flap conformation, which does not normally occur in apo subtype B. Other secondary mutations can favor drug resistance. For example, 13 out of the 19 hydrophobic core residues are associated with drug resistance acquisition, despite the absence of direct contacts with inhibitors [94]. Importantly, the accumulation of mutations in HIV-1 PR does not necessarily lead to increased drug resistance [48,93,99]. Indeed, alterations in inhibitor-protein interactions caused by secondary mutations correlate only partly with $K_i$ variations. The explanation most often put forward is that these secondary mutations alter the dynamic ensemble of the protein, which propagates through the entire structure (Figure 3b), sometimes even affecting the dynamically constrained catalytic Asp25 residue [90,92,93,98]. Thus, if a secondary mutation does not alter protein dynamics, it will have no effect on inhibitor binding. Nonetheless, increased drug resistance is the usual outcome of mutational accumulation. For instance, a DEER study comparing a pediatric AIDS patient-derived HIV-1 PR before and after protease inhibitor treatment found that the accumulation of secondary mutations shifts the populations of the conformers to favor the open-like conformations, when it originally favored the closed conformation [100]. These authors emphasize the importance of the L63P replacement on the behavior of the protein in the pretreatment sample, arguing that it likely affects hydrophobic core packing, and thus favoring the closed conformation in the apo protein.

These results illustrate how knowledge of the conformational transitions in HIV-1 PR may help develop new drug strategies against AIDS, further exemplifying the importance of concerted and organized molecular flexibility for the proper catalytic function in this
enzyme. Targeting conformational transitions has already provided promising treatment alternatives for new HIV-1 PR drug leads. The use of MD to screen for new competitive or allosteric inhibitors that target MD-identified, transiently open surface cavities was shown to be effective against mutants of HIV-1 PR [92,101,102]. In one case, a designed allosteric inhibitor was demonstrated to be effective *in vitro* using enzymatic essays, even against drug-resistant variants [102]. Interestingly, the allosteric inhibitor strategy used in these cases was designed to lock a conformation in place, thus preventing the catalytic cycle from talking place.

## 3.2. DNA Polymerase β

DNA polymerase β (Pol β) is a 39 kDa monomeric enzyme involved in the base excision repair of damaged DNA. Pol β prefers short-gapped or single-base-gapped DNA substrates [103] and is responsible for the repair of around 20,000 DNA lesions per cell per day [104]. Mutations of some residues in Pol β have been shown to increase the error rate of the enzyme; further, overexpression of Pol β has been implicated in nearly 30% of human cancers [105,106]. Incorporation of correct nucleotides is an essential step in the catalytic mechanism of the enzyme. The structure of Pol β consists of two domains: an 8 kDa lyase domain involved in the removal of deoxyribose 5'-phosphate, and a 31 kDa polymerase domain which performs template-directed DNA synthesis [107]. The polymerase domain is further composed of nucleoside 5'-triphosphate selection, nucleotide transferase, and DNA-binding subdomains, corresponding to the fingers, palm, and thumb regions, respectively. The catalytic mechanism of Pol β, which requires two Mg²⁺ ions as cofactors, is shown in Figure 4.

Crystallographic studies have shown that Pol β undergoes large conformational changes upon substrate binding [110–115] and suggest an induced fit mechanism for the binding of correct nucleotides [111]. The apo form of the enzyme is characterized by an extended conformation of the lyase domain (Figure 5a), which adopts a compact open conformation upon DNA binding (Figure 5b). Binding of the nucleotide 5'-triphosphate (dNTP) promotes closing of the N-helix in the finger subdomain (Figure 5c). NMR studies were recently performed to characterize the millisecond timescale conformational exchange experienced by the apo and DNA-bound binary complex forms of Pol β [8]. The authors of that study showed that conformational exchange in the apo form occurs primarily in the lyase domain and the base of the thumb in the polymerase domain, and corresponds to residues that are either in direct contact or in close proximity to the substrate in the binary complex. In contrast, the number of residues displaying conformational exchange in the binary complex with DNA was significantly reduced. These residues are primarily located in the polymerase domain, suggesting that DNA binding limits the motions in the lyase domain (Figure 6). Chemical shift variations were located in the lyase domain, associated with large conformational changes experienced by Pol β upon DNA binding (Figure 5a,b). Comparison of the conformational exchange relaxation rates with chemical shift variations between the apo and binary forms indicates that the enzyme in the apo form samples conformations adopted by the DNA-bound binary form, thus contradicting the induced fit mechanism proposed by the aforementioned crystallography experiments, at least for the formation of the binary complex (step 1 in Figure 4). Interestingly, residues displaying conformational
exchange in the apo and binary complex forms correspond to mutational sites known to affect Pol β activity and identified in colon cancer tumors [116]. Taken together, these results hint at the potential role of conformational motions in the catalytic function of Pol β.

Conformational changes in both Pol β and substrates have been suggested to play a critical role in the selection and incorporation of correct nucleotides [117]. Experimental [109,111–113,118] and computational [119,120] studies suggest an induced fit mechanism for the selection of correct nucleotides (step 2 in Figure 4), which form a Watson-Crick base pair with the DNA substrate. The mechanism involved in the selection of correct dNTP binding and discrimination is a hotly debated topic [109]. Crystal structures with open (Figure 5d), closed (Figure 5e), and intermediate conformations have been obtained with mismatched or incorrect dNTPs, such as dA that cannot form a correct base pair with dG of the DNA substrate. Previous fluorescence studies [118,121,122] showed movements in Pol β upon incorporation of dNTPs and cofactors, while more recent studies [108] indicate the rapid movement (ms timescale) of the fingers domain upon binding to correct dNTPs. Incorporation of incorrect dNTPs did not induce these conformational changes.

Moscato *et al.* characterized the effects of Pol β binding to correct and incorrect dNTPs on the millisecond timescale using NMR relaxation dispersion experiments [109]. Binding of correct dNTP to Pol β led to large chemical shift variations localized in the fingers subdomain, similar to previous observations [111,114]. No residue displayed conformational exchange, suggesting a stable conformation of the ternary complex on the millisecond timescale, similar to that observed for the binary complex. In contrast, chemical shift changes upon addition of the incorrect dNTP were much smaller than those observed for the correct dNTP. Comparison of the direction of chemical shift variations during titrations, performed using chemical shift perturbation [123], revealed that the conformation of the mismatching complex is not on-pathway to the closed conformation observed upon binding of the correct nucleotide. Further, binding of the mismatched dNTP led to enhanced millisecond dynamics of 10 residues with conformational exchange rates varying between 500 and 2700 s⁻¹. These results support an induced fit mechanism for the ternary complex (step 2 in Figure 4), where binding of the correct nucleotide promotes a shift from the open binary complex form to a closed ternary complex conformation.

## 4. Concluding Remarks

The intrinsic dynamical properties of enzymes have been shown to impact the catalytic function in a variety of well-characterized systems. Advances in experimental and computational approaches continue to offer insights into the important role of conformational fluctuations in function on timescales relevant for enzyme catalysis and other biological events. Conformational exchange on the catalytically relevant timescale facilitates sampling of excited-state conformations or sub-states, aiding in various steps along the reaction pathway, such as substrate recognition, binding, and product release. In this review, we depicted the important effects of conformational dynamics on correct substrate selection and substrate/inhibitor recognition and binding in two enzyme systems. These observations also emphasize the role of dynamic fluctuations in tuning the
conformational landscape for efficient catalysis, and present a potential approach for using
dynamics to modulate function and for designing better inhibitors in drug design.

## Acknowledgments

This work was supported by the National Institute of General Medical Sciences (NIGMS) of the US National Institutes of Health (NIH) under award number R01GM105978 (to N.D.) and a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant under award number RGPIN-2016-05557 (to N.D.). D.N.B. is the recipient of an NSERC Postgraduate Scholarships-Doctoral (PGS D) Award and N.D. holds a Fonds de Recherche Québec-Santé (FRQS) Research Scholar Junior 2 Career Award.

## References

1. Garcia-Viloca M, Gao J, Karplus M, Truhlar DG. How enzymes work: Analysis by modern rate theory and computer simulations. Science. 2004; 303:186-195. [PubMed: 14716003]

2. Wolfenden R. Degrees of difficulty of water-consuming reactions in the absence of enzymes. Chem Rev. 2006; 106:3379-3396. [PubMed: 16895333]

3. Hammes GG, Chang YC, Oas TG. Conformational selection or induced fit: A flux description of reaction mechanism. Proc Natl Acad Sci USA. 2009; 106:13737-13741. [PubMed: 19666553]

4. Henzler-Wildman K, Kern D. Dynamic personalities of proteins. Nature. 2007; 450:964-972. [PubMed: 18075575]

5. Watt ED, Shimada H, Kovrigin EL, Loria JP. The mechanism of rate-limiting motions in enzyme function. Proc Natl Acad Sci USA. 2007; 104:11981-11986. [PubMed: 17615241]

6. Venkitakrishnan RP, Zaborowski E, McElheny D, Benkovic SJ, Dyson HJ, Wright PE. Conformational changes in the active site loops of dihydrofolate reductase during the catalytic cycle. Biochemistry. 2004; 43:16046-16055. [PubMed: 15609999]

7. Boehr DD. Promiscuity in protein-RNA interactions: Conformational ensembles facilitate molecular recognition in the spliceosome. Bioessays. 2012; 34:174-180. [PubMed: 22144099]

8. Berlow RB, Swain M, Dalal S, Sweasy JB, Loria JP. Substrate-dependent millisecond domain motions in DNA polymerase β. J Mol Biol. 2012; 419:171-182. [PubMed: 22446382]

9. Kay LE. NMR studies of protein structure and dynamics. J Magn Reson. 2005; 173:193-207. [PubMed: 15780912]

10. Vallurupalli P, Hansen DF, Kay LE. Structures of invisible, excited protein states by relaxation dispersion NMR spectroscopy. Proc Natl Acad Sci USA. 2008; 105:11766-11771. [PubMed: 18701719]

11. Ramanathan A, Savol A, Burger V, Chennubhotla CS, Agarwal PK. Protein conformational populations and functionally relevant substates. Acc Chem Res. 2014; 47:149-156. [PubMed: 23988159]

12. Vanwart AT, Eargle J, Luthey-Schulten Z, Amaro RE. Exploring residue component contributions to dynamical network models of allostery. J Chem Theory Comput. 2012; 8:2949-2961. [PubMed: 23139645]

13. Gagné D, Charest LA, Morin S, Kovrigin EL, Doucet N. Conservation of flexible residue clusters among structural and functional enzyme homologues. J Biol Chem. 2012; 287:44289-44300. [PubMed: 23135272]

14. Eisenmesser EZ, Millet O, Labeikovsky W, Korzhnev DM, Wolf-Watz M, Bosco DA, Skalicky JJ, Kay LE, Kern D. Intrinsic dynamics of an enzyme underlies catalysis. Nature. 2005; 438:117-121. [PubMed: 16267559]

15. Agarwal PK, Geist A, Gorin A. Protein dynamics and enzymatic catalysis: Investigating the peptidyl-prolyl cis-trans isomerization activity of cyclophilin A. Biochemistry. 2004; 43:10605-10618. [PubMed: 15311922]

16. Goodey NM, Benkovic SJ. Allosteric regulation and catalysis emerge via a common route. Nat Chem Biol. 2008; 4:474-482. [PubMed: 18641628]
17. Suel GM, Lockless SW, Wall MA, Ranganathan R. Evolutionarily conserved networks of residues mediate allosteric communication in proteins. Nat Struct Biol. 2003; 10:59–69. [PubMed: 12483203]

18. Reynolds KA, McLaughlin RN, Ranganathan R. Hot spots for allosteric regulation on protein surfaces. Cell. 2011; 147:1564–1575. [PubMed: 22196731]

19. Cole R, Loria JP. Evidence for flexibility in the function of ribonuclease A. Biochemistry. 2002; 41:6072–6081. [PubMed: 11994002]

20. Gagné D, Doucet N. Structural and functional importance of local and global conformational fluctuations in the RNase A superfamily. FEBS J. 2013; 280:5596–5607. [PubMed: 23763751]

21. Pareek V, Samanta M, Joshi NV, Balaram H, Murthy MR, Balaram P. Connecting active site loop conformations and catalysis in triosephosphate isomerase: Insights from a rare variation at residue 96 in the plasmodial enzyme. Chembiochem. 2016 in press.

22. Torbeev VY, Raghuraman H, Hamelberg D, Tonelli M, Westler WM, Perozo E, Kent SB. Protein conformational dynamics in the mechanism of HIV-1 protease catalysis. Proc Natl Acad Sci USA. 2011; 108:20982–20987. [PubMed: 22158985]

23. Kar G, Keskin O, Gursoy A, Nussinov R. Allostery and population shift in drug discovery. Curr Opin Pharmacol. 2010; 10:715–722. [PubMed: 20884293]

24. Agarwal PK. Enzymes: An integrated view of structure, dynamics and function. Microb Cell Fact. 2006; 5:2. [PubMed: 16409630]

25. Kohen A. Role of dynamics in enzyme catalysis: Substantial versus semantic controversies. Acc Chem Res. 2015; 48:466–473. [PubMed: 25539442]

26. Kamerlin SC, Warshel A. At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis? Proteins. 2010; 78:1339–1375. [PubMed: 20099310]

27. Nagel ZD, Klinman JP. A 21st century revisionist's view at a turning point in enzymology. Nat Chem Biol. 2009; 5:543–550. [PubMed: 19620995]

28. Hanoian P, Liu CT, Hammes-Schiffer S, Benkovic S. Perspectives on electrostatics and conformational motions in enzyme catalysis. Acc Chem Res. 2015; 48:482–489. [PubMed: 25565178]

29. Grutsch S, Bruschweiler S, Tollinger M. NMR methods to study dynamic allostery. PLoS Comput Biol. 2016; 12:e1004620. [PubMed: 26964042]

30. Palmer AG 3rd. Chemical exchange in biomacromolecules: Past, present, and future. J Magn Reson. 2014; 241:3–17. [PubMed: 24656076]

31. Kleckner IR, Foster MP. An introduction to NMR-based approaches for measuring protein dynamics. Biochim Biophys Acta. 2011; 1814:942–968. [PubMed: 21059410]

32. Loria JP, Rance M, Palmer AG 3rd. A trosy cpmg sequence for characterizing chemical exchange in large proteins. J Biomol NMR. 1999; 15:151–155. [PubMed: 10605088]

33. Loria JP, Rance M, Palmer AG. A relaxation-compensated carr-purcell-meiboom-gill sequence for characterizing chemical exchange by NMR spectroscopy. J Am Chem Soc. 1999; 121:2331–2332.

34. Oyen D, Fenwick RB, Stanfield RL, Dyson HJ, Wright PE. Cofactor-mediated conformational dynamics promote product release from Escherichia Coli dihydrofolate reductase via an allosteric pathway. J Am Chem Soc. 2015; 137:9459–9468. [PubMed: 26147643]

35. Bhabha G, Biel JT, Fraser JS. Keep on moving: Discovering and perturbing the conformational dynamics of enzymes. Acc Chem Res. 2015; 48:423–430. [PubMed: 25539415]

36. Doucet N. Can enzyme engineering benefit from the modulation of protein motions? Lessons learned from NMR relaxation dispersion experiments. Protein Pept Lett. 2011; 18:336–343. [PubMed: 21222641]

37. Palmer AG 3rd. NMR probes of molecular dynamics: Overview and comparison with other techniques. Annu Rev Biophys Biomol Struct. 2001; 30:129–155. [PubMed: 11340055]

38. Palmer AG 3rd. NMR characterization of the dynamics of biomacromolecules. Chem Rev. 2004; 104:3623–3640. [PubMed: 15303831]

39. Manley G, Loria JP. NMR insights into protein allostery. Arch Biochem Biophys. 2012; 519:223–231. [PubMed: 22198279]
40. Luk LY, Loveridge EJ, Allemann RK. Protein motions and dynamic effects in enzyme catalysis. Phys Chem Chem Phys. 2015; 17:30817–30827. [PubMed: 25854702]

41. Lipari G, Szabo A. Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity. J Am Chem Soc. 1982; 104:4546–4559.

42. Lipari G, Szabo A. Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 2. Analysis of experimental results. J Am Chem Soc. 1982; 104:4559–4570.

43. Morin S. A practical guide to protein dynamics from 15N spin relaxation in solution. Prog Nucl Magn Reson Spectrosc. 2011; 59:245–262. [PubMed: 21920220]

44. Fisette O, Lagüe P, Gagné S, Morin S. Synergistic applications of MD and NMR for the study of biological systems. J Biomed Biotechnol. 2012; 2012:254208. [PubMed: 22319241]

45. Mandel AM, Akke M, Palmer AG 3rd. Backbone dynamics of Escherichia Coli Ribonuclease HI: Correlations with structure and function in an active enzyme. J Mol Biol. 1995; 246:144–163. [PubMed: 7531772]

46. Stephenson JD, Kenyon JC, Symmons MF, Lever AM. Characterizing 3D RNA structure by single molecule FRET. Methods. 2016 in press.

47. Huang X, de Vera IM, Veloro AM, Blackburn ME, Kear JL, Carter JD, Rocca JR, Simmerling C, Dunn BM, Fanucci GE. Inhibitor-induced conformational shifts and ligand-exchange dynamics for HIV-1 protease measured by pulsed epr and NMR spectroscopy. J Phys Chem B. 2012; 116:14235–14244. [PubMed: 23167829]

48. Huang X, Britto MD, Kear-Scott JL, Boone CD, Rocca JR, Simmerling C, McKenna R, Bieri M, Gooley PR, Dunn BM, et al. The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics. J Biol Chem. 2014; 289:17203–17214. [PubMed: 24742668]

49. Jeschke G, Polyhach Y. Distance measurements on spin-labelled biomacromolecules by pulsed electron paramagnetic resonance. Phys Chem Chem Phys. 2007; 9:1895–1910. [PubMed: 17431518]

50. Boehr DD, McElheny D, Dyson HJ, Wright PE. Millisecond timescale fluctuations in dihydrofolate reductase are exquisitely sensitive to the bound ligands. Proc Natl Acad Sci USA. 2010; 107:1373–1378. [PubMed: 20080605]

51. Boehr DD, Liu X, Yang X. Targeting structural dynamics of the RNA-dependent RNA polymerase for anti-viral strategies. Curr Opin Virol. 2014; 9:194–200. [PubMed: 25224392]

52. Voelz VA, Bowman GR, Beauchamp K, Pande VS. Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1–39). J Am Chem Soc. 2010; 132:1526–1528. [PubMed: 20070076]

53. Bowman GR, Voelz VA, Pande VS. Atomistic folding simulations of the five-helix bundle protein λ6–85. J Am Chem Soc. 2011; 133:664–667. [PubMed: 21174461]

54. Voelz VA, Jager M, Yao S, Chen Y, Zhu L, Waldauer SA, Bowman GR, Friedrichs M, Bakajin O, Lapidus LJ; et al. Slow unfolded-state structuring in Acyl-coA binding protein folding revealed by simulation and experiment. J Am Chem Soc. 2012; 134:12565–12577. [PubMed: 22747188]

55. Shaw DE, Maragakis P, Lindorff-Larsen K, Piana S, Dror RO, Eastwood MP, Bank JA, Jumper JM, Salmon JK, Shan Y; et al. Atomic-level characterization of the structural dynamics of proteins. Science. 2010; 330:341–346. [PubMed: 20947758]

56. Bucher D, Grant BJ, Markwick PR, McCammon JA. Accessing a hidden conformation of the maltose binding protein using accelerated molecular dynamics. PLoS Comput Biol. 2011; 7:e1002034. [PubMed: 21533070]

57. Hamelberg D, Mongan J, McCammon JA. Accelerated molecular dynamics: A promising and efficient simulation method for biomolecules. J Chem Phys. 2004; 120:11919–11929. [PubMed: 15268227]

58. Wojtas-Niziurski W, Meng Y, Roux B, Berneche S. Self-learning adaptive umbrella sampling method for the determination of free energy landscapes in multiple dimensions. J Chem Theory Comput. 2013; 9:1885–1895. [PubMed: 23814508]
59. Da LT, Wang D, Huang X. Dynamics of pyrophosphate ion release and its coupled trigger loop motion from closed to open state in RNA polymerase II. J Am Chem Soc. 2012; 134:2399-2406. [PubMed: 22206270]

60. Tagle AM. Percutaneous encoscopic gastrostomy. Rev Gastroenterol Peru. 1998; 18:56-61. [PubMed: 12271336]

61. Sutto L, Gervasio FL. Effects of oncogenic mutations on the conformational free-energy landscape of EGFR kinase. Proc Natl Acad Sci USA. 2013; 110:10616-10621. [PubMed: 23754386]

62. Abrams C, Bussi G. Enhanced sampling in molecular dynamics using metadynamics, replica-exchange, and temperature-acceleration. Entropy. 2014; 16:163-199.

63. Doshi U, Hamelberg D. Towards fast, rigorous and efficient conformational sampling of biomolecules: Advances in accelerated molecular dynamics. Biochim Biophys Acta. 2015; 1850:878-888. [PubMed: 25153688]

64. Doshi U, McGowan LC, Ladani ST, Hamelberg D. Resolving the complex role of enzyme conformational dynamics in catalytic function. Proc Natl Acad Sci USA. 2012; 109:5699-5704. [PubMed: 22451902]

65. Torrie GM, Valleau JP. Nonphysical sampling distributions in monte carlo free-energy estimation: Umbrella sampling. J Comput Phys. 1977; 23:187-199.

66. Camilloni C, Schaal D, Schweimer K, Schwarzinger S, De Simone A. Energy landscape of the prion protein helix 1 probed by metadynamics and NMR. Biophys J. 2012; 102:158-167. [PubMed: 22225810]

67. Andersen OJ, Grouleff J, Needham P, Walker RC, Jensen F. Toward an enhanced sampling molecular dynamics method for studying ligand-induced conformational changes in proteins. J Phys Chem B. 2015; 119:14594-14603. [PubMed: 26482713]

68. Singhal N, Snow CD, Pande VS. Using path sampling to build better markovian state models: Predicting the folding rate and mechanism of a tryptophan zipper β hairpin. J Chem Phys. 2004; 121:415-425. [PubMed: 15260562]

69. Bowman GR, Beauchamp KA, Boxer G, Pande VS. Progress and challenges in the automated construction of Markov state models for full protein systems. J Chem Phys. 2009; 131:124101. [PubMed: 19791846]

70. Noé F, Fischer S. Transition networks for modeling the kinetics of conformational change in macromolecules. Curr Opin Struct Biol. 2008; 18:154-162. [PubMed: 18378442]

71. Ramanathan A, Agarwal PK. Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis. PLoS Biol. 2011; 9:e1001193. [PubMed: 22087074]

72. Shukla D, Hernandez CX, Weber JK, Pande VS. Markov state models provide insights into dynamic modulation of protein function. Acc Chem Res. 2015; 48:414-422. [PubMed: 25625937]

73. Chodera JD, Noe F. Markov state models of biomolecular conformational dynamics. Curr Opin Struct Biol. 2014; 25:135-144. [PubMed: 24836551]

74. Morcos F, Chatterjee S, McClendon CL, Brenner PR, Lopez-Rendon R, Zintsmaster J, Ercsey-Ravasz M, Sweet CR, Jacobson MP, Peng JW; et al. Modeling conformational ensembles of slow functional motions in Pin1-WW. PLoS Comput Biol. 2010; 6:e1001015. [PubMed: 21152000]

75. Ramanathan A, Savol AJ, Langmead CJ, Agarwal PK, Chennubhotla CS. Discovering conformational sub-states relevant to protein function. PLoS ONE. 2011; 6:e15827. [PubMed: 21297978]

76. Gagné D, French RL, Narayanan C, Simonovic M, Agarwal PK, Doucet N. Perturbation of the conformational dynamics of an active-site loop alters enzyme activity. Structure. 2015; 23:2256-2266. [PubMed: 26655472]

77. Camilloni C, Robustelli P, De Simone A, Cavalli A, Vendruscolo M. Characterization of the conformational equilibrium between the two major substates of RNase A using NMR chemical shifts. J Am Chem Soc. 2012; 134:3968-3971. [PubMed: 22320129]

78. De Simone A, Montalvao RW, Dobson CM, Vendruscolo M. Characterization of the interdomain motions in hen lysozyme using residual dipolar couplings as replica-averaged structural restraints in molecular dynamics simulations. Biochemistry. 2013; 52:6480-6486. [PubMed: 23941501]
79. Camilloni C, Cavalli A, Vendruscolo M. Assessment of the use of NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations to characterize the dynamics of proteins. J Phys Chem B. 2013; 117:1838–1843. [PubMed: 23327201]

80. Sundquist WI, Krausslich HG. HIV-1 assembly, budding, and maturation. Cold Spring Harb Perspect Med. 2012; 2:a006924. [PubMed: 22762019]

81. Ghosh AK, Osswald HL, Prato G. Recent progress in the development of HIV-1 protease inhibitors for the treatment of HIV/AIDS. J Med Chem. 2016 in press.

82. Spinelli S, Liu QZ, Alzari PM, Hirel PH, Poljak RJ. The three-dimensional structure of the aspartyl protease from the HIV-1 isolate BRU. Biochimie. 1991; 73:1391–1396. [PubMed: 1799632]

83. Pastor N, Amero C. Information flow and protein dynamics: The interplay between nuclear magnetic resonance spectroscopy and molecular dynamics simulations. Front Plant Sci. 2015; 6:306. [PubMed: 25999971]

84. Nicholson LK, Yamazaki T, Torchia DA, Grzesiek S, Bax A, Stahl SJ, Kaufman JD, Wingfield PT, Lam PY, Jadhav PK; et al. Flexibility and function in HIV-1 protease. Nat Struct Biol. 1995; 2:274–280. [PubMed: 7796263]

85. Ishima R, Louis JM. A diverse view of protein dynamics from NMR studies of HIV-1 protease flaps. Proteins. 2008; 70:1408–1415. [PubMed: 17894346]

86. Xia J, Deng NJ, Levy RM. NMR relaxation in proteins with fast internal motions and slow conformational exchange: Model-free framework and Markov state simulations. J Phys Chem B. 2013; 117:6625–6634. [PubMed: 23638941]

87. Ishima R, Torchia DA, Lynch SM, Gronenborn AM, Louis JM. Solution structure of the mature HIV-1 protease monomer: Insight into the tertiary fold and stability of a precursor. J Biol Chem. 2003; 278:43311–43319. [PubMed: 12933791]

88. Deng NJ, Zheng W, Gallicchio E, Levy RM. Insights into the dynamics of HIV-1 protease: A kinetic network model constructed from atomistic simulations. J Am Chem Soc. 2011; 133:9387–9394. [PubMed: 21561098]

89. Shukla D, Meng Y, Roux B, Pande VS. Activation pathway of SRC kinase reveals intermediate states as targets for drug design. Nat Commun. 2014; 5:3397. [PubMed: 24584478]

90. Roche J, Louis JM, Bax A. Conformation of inhibitor-free HIV-1 protease derived from nmr spectroscopy in a weakly oriented solution. Chembiochem. 2015; 16:214–218. [PubMed: 25470009]

91. Yamazaki T, Hinck AP, Wang YX, Nicholson LK, Torchia DA, Wingfield P, Stahl SJ, Kaufman JD, Chang CH, Domaille PJ; et al. Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy. Protein Sci. 1996; 5:495–506. [PubMed: 8868486]

92. Leonis G, Steinbrecher T, Papadopoulos MG. A contribution to the drug resistance mechanism of darunavir, amprenavir, indinavir, and saquinavir complexes with HIV-1 protease due to flap mutation I50V: A systematic MM-PBSA and thermodynamic integration study. J Chem Inf Model. 2013; 53:2141–2153. [PubMed: 23834142]

93. Gupta A, Jamal S, Goyal S, Jain R, Wahi D, Grover A. Structural studies on molecular mechanisms of nelfinavir resistance caused by non-active site mutation V77I in HIV-1 protease. BMC Bioinform. 2015; 16:S10.

94. Mittal S, Cai Y, Nalam MN, Bolon DN, Schiffer CA. Hydrophobic core flexibility modulates enzyme activity in HIV-1 protease. J Am Chem Soc. 2012; 134:4163–4168. [PubMed: 22295904]

95. Goldfarb NE, Ohanessian M, Biswas S, McGee TD Jr, Mahon BP, Ostrov DA, Garcia J, Tang Y, McKenna R, Roitberg A; et al. Defective hydrophobic sliding mechanism and active site expansion in HIV-1 protease drug resistant variant Gly48Thr/Leu89Met: Mechanisms for the loss of saquinavir binding potency. Biochemistry. 2015; 54:422–433. [PubMed: 25513833]

96. Lee J, Goodey NM. Catalytic contributions from remote regions of enzyme structure. Chem Rev. 2011; 111:7595–7624. [PubMed: 21923192]

97. Wu TD, Schiffer CA, Gonzales MJ, Taylor J, Kantor R, Chou S, Israelski D, Zolopa AR, Fessel WJ, Shafer RW. Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments. J Virol. 2003; 77:4836–4847. [PubMed: 12663790]
98. Ragland DA, Nalivaika EA, Nalam MN, Prachanronarong KL, Cao H, Bandaranayake RM, Cai Y, Kurt-Yilmaz N, Schiffer CA. Drug resistance conferred by mutations outside the active site through alterations in the dynamic and structural ensemble of HIV-1 protease. J Am Chem Soc. 2014; 136:11956–11963. [PubMed: 25091085]

99. Tzoupis H, Leonis G, Avramopoulos A, Mavromoustakos T, Papadopoulos MG. Systematic molecular dynamics, MM-PBSA, and ab initio approaches to the saquinavir resistance mechanism in HIV-1 PR due to 11 double and multiple mutations. J Phys Chem B. 2014; 118:9538–9552. [PubMed: 25036111]

100. Carter JD, Gonzales EG, Huang X, Smith AN, de Vera IM, D'Amore PW, Rocca JR, Goodenow MM, Dunn BM, Fanucci GE. Effects of pre and post therapy drug-pressure selected mutations on HIV-1 protease conformational sampling. FEBS Lett. 2014; 588:3123–3128. [PubMed: 24983495]

101. Kunze J, Todoroff N, Schneider P, Rodrigues T, Geppert T, Reisen F, Schreuder H, Saas J, Hessler G, Baringhaus KH; et al. Targeting dynamic pockets of HIV-1 protease by structure-based computational screening for allosteric inhibitors. J Chem Inf Model. 2014; 54:987–991. [PubMed: 24528206]

102. Ung PM, Dunbar JB Jr, Gestwicki JE, Carlson HA. An allosteric modulator of HIV-1 protease shows equipotent inhibition of wild-type and drug-resistant proteases. J Med Chem. 2014; 57:6468–6478. [PubMed: 25062388]

103. Beard WA, Wilson SH. Structure and mechanism of DNA polymerase β. Chem Rev. 2006; 106:361–382. [PubMed: 16464010]

104. Barnes DE, Lindahl T. Repair and genetic consequences of endogenous DNA base damage in mammalian cells. Annu Rev Genet. 2004; 38:445–476. [PubMed: 15568983]

105. Starcevic D, Dalal S, Sweasy JB. Is there a link between DNA polymerase β and cancer? Cell Cycle. 2004; 3:998–1001. [PubMed: 15280658]

106. Loeb LA, Monnat RJ Jr. DNA polymerases and human disease. Nat Rev Genet. 2008; 9:594–604. [PubMed: 18626473]

107. Beard WA, Shock DD, Yang XP, DeLauder SF, Wilson SH. Loss of DNA polymerase β stacking interactions with templating purines, but not pyrimidines, alters catalytic efficiency and fidelity. J Biol Chem. 2002; 277:8235–8242. [PubMed: 11756435]

108. Towle-Weicksel JB, Dalal S, Sohl CD, Doublie S, Anderson KS, Sweasy JB. Fluorescence resonance energy transfer studies of DNA polymerase β: The critical role of fingers domain movements and a novel non-covalent step during nucleotide selection. J Biol Chem. 2014; 289:16541–16550. [PubMed: 24764311]

109. Moscato B, Swain M, Loria JP. Induced fit in the selection of correct versus incorrect nucleotides by DNA polymerase β. Biochemistry. 2016; 55:382–395. [PubMed: 26678253]

110. Sawaya MR, Pelletier H, Kumar A, Wilson SH, Kraut J. Crystal structure of rat DNA polymerase β: Evidence for a common polymerase mechanism. Science. 1994; 264:1930–1935. [PubMed: 7516581]

111. Sawaya MR, Prasad R, Wilson SH, Kraut J, Pelletier H. Crystal structures of human DNA polymerase β complexed with gapped and nicked DNA: Evidence for an induced fit mechanism. Biochemistry. 1997; 36:11205–11215. [PubMed: 9287163]

112. Batra VK, Beard WA, Shock DD, Krahn JM, Pedersen LC, Wilson SH. Magnesium-induced assembly of a complete DNA polymerase catalytic complex. Structure. 2006; 14:757–766. [PubMed: 16615916]

113. Batra VK, Beard WA, Shock DD, Pedersen LC, Wilson SH. Structures of DNA polymerase β with active-site mismatches suggest a transient abasic site intermediate during misincorporation. Mol Cell. 2008; 30:315–324. [PubMed: 18471977]

114. Freudenthal BD, Beard WA, Wilson SH. Structures of dNTP intermediate states during DNA polymerase active site assembly. Structure. 2012; 20:1829–1837. [PubMed: 22959623]

115. Freudenthal BD, Beard WA, Shock DD, Wilson SH. Observing a DNA polymerase choose right from wrong. Cell. 2013; 154:157–168. [PubMed: 23827680]
116. Donigan KA, Sun KW, Nemec AA, Murphy DL, Cong X, Northrup V, Zelterman D, Sweasy JB. Human POLB gene is mutated in high percentage of colorectal tumors. J Biol Chem. 2012; 287:23830-23839. [PubMed: 22577134]

117. Beard WA, Shock DD, Batra VK, Prasad R, Wilson SH. Substrate-induced DNA polymerase β activation. J Biol Chem. 2014; 289:31411-31422. [PubMed: 25261471]

118. Bakhtina M, Lee S, Wang Y, Dunlap C, Lamarche B, Tsai MD. Use of viscogens, dNTPαS, and rhodium(III) as probes in stopped-flow experiments to obtain new evidence for the mechanism of catalysis by DNA polymerase β. Biochemistry. 2005; 44:5177-5187. [PubMed: 15794655]

119. Arora K, Beard WA, Wilson SH, Schlick T. Mismatch-induced conformational distortions in polymerase β support an induced-fit mechanism for fidelity. Biochemistry. 2005; 44:13328-13341. [PubMed: 16201758]

120. Kirmizialtin S, Nguyen V, Johnson KA, Elber R. How conformational dynamics of DNA polymerase select correct substrates: Experiments and simulations. Structure. 2012; 20:618-627. [PubMed: 22483109]

121. Zhong X, Patel SS, Werneburg BG, Tsai MD. DNA polymerase β: Multiple conformational changes in the mechanism of catalysis. Biochemistry. 1997; 36:11891-11900. [PubMed: 9305982]

122. Balbo PB, Wang EC, Tsai MD. Kinetic mechanism of active site assembly and chemical catalysis of DNA polymerase β. Biochemistry. 2011; 50:9865-9875. [PubMed: 22010960]

123. Williamson MP. Using chemical shift perturbation to characterise ligand binding. Prog Nucl Magn. 2013; 73:1-16.
<figure><img src="image_3.png" /><figcaption>**Figure 1.**</figcaption></figure>

Conformational exchange experienced by the Met20 loop in *E. coli* DHFR, as probed by ¹⁵N-CPMG NMR relaxation dispersion experiments. The CPMG experiment is particularly well suited to extract exchange rates (*k*<sub>ex</sub>) between two or more conformations in solution on the timescale of catalysis (ms) in many enzyme systems, offering a measure of comparison between conformational exchange experienced by the enzyme and its catalytic rate (*k*<sub>cat</sub>). This method can also provide quantitative information on low-populated “invisible” excited sub-states and theft populations in solution (*p*<sub>A</sub> and *p*<sub>B</sub>). The cartoon representation illustrates a simplified view of the closed (green) and occluded (magenta) conformations sampled by the Met20 loop in *E. coli* DHFR as it catalyzes hydride transfer, an atomic-scale movement essential to bacterial DHFR function and correlated with substrate/cofactor recognition and turnover in this enzyme (reviewed in [28,35,40]). Depicted PDB structures are 1RX2 (green) and 1RX7 (magenta).
<figure><img src="image_4.png" /><figcaption>**Figure 2.**</figcaption></figure>

Conformational states experienced by HIV-1 PR as it proceeds through its catalytic cycle. Crystal structures of HIV-1 PR in the semi-open (PDB 1HHR), closed (PDB 1HVR) and wide open (PDB 1TW7) conformations. The catalytic Asp25 residues from each protomer are depicted as cyan sticks and the “flaps” are depicted in red. The hinge loop, corresponding to residues 34–42, is shown in blue. The curled/tucked conformation of HIV-1 PR is not depicted, since it has never been crystallographically resolved.
<figure><img src="image_5.png" /><figcaption>**Figure 3.**</figcaption></figure>

Location of drug resistance and conformational mutations in HIV-1 PR. Structural mapping of various mutations on the crystal structure of apo HIV-1 PR (1HHR). The catalytic Asp25 residue is depicted as cyan sticks on each enzyme protomer. All spheres represent a residue documented to undergo mutation in HIV-1 PR [97]. (a) Location of drug resistance-inducing mutations. Magenta spheres represent mutations that impair inhibitor binding through direct active-site contacts (primary mutations) and gray spheres represent distal indirect mutations (secondary mutations); (b) Location of mutations documented to affect the conformational sampling of HIV-1 PR, depicted as orange spheres. Gray spheres are other drug resistance-inducing mutations. All mutations were compiled from references listed in the text.
<figure><img src="image_6.png" /><figcaption>**Figure 4.**</figcaption></figure>

Schematic representation of the catalytic mechanism of Pol β. $E_o$ and $E_c$ correspond to the open and closed conformations of the enzyme. Sequential binding of the DNA ($DNA_n$) and dNTP leads to the formation of the binary (step 1) and ternary complexes (step 2), respectively. Subsequent conformational rearrangement of the enzyme ($E_c$) aligns the active site for catalysis (step 3). The chemical step of catalysis is followed by the release of pyrophosphate ($PP_i$) and the DNA product ($DNA_{n+1}$) (step 4). Figure adapted from [108,109].
a

<figure><img src="image_7.png" /></figure>

b

<figure><img src="image_8.png" /></figure>

c

<figure><img src="image_9.png" /></figure>

d

<figure><img src="image_10.png" /></figure>

e

<figure><img src="image_11.png" /><figcaption>**Figure 5.**</figcaption></figure>

Structure of free and ligand-bound states of Pol β. The lyase domain is shown in red while the thumb, palm, and fingers subdomains of the polymerase domain are shown in blue, green and yellow, respectively. Pol β in the (a) apo (PDB 1BPD) form is characterized by an open conformation; (b) binary complex (PDB 1BPX) formation is accompanied by a 40 Å conformational shift of the lyase domain to form a compact “open” conformation; (c) matched dNTP-bound ternary (PDB 1BPY) complex shows a 10 Å shift of the N-helix of the fingers domain relative to the binary form (in orange), leading to the formation of a “closed” conformation. Binding of mismatched dNTPs has been shown to adopt (d) closed (PDB3C2M) and (e) open (PDB 4F5P) conformations of the N-helix. For comparison, the N-helix from the binary form is shown in orange in the ternary complexes in c-e, illustrating subtle yet functionally significant conformational exchange in this enzyme. DNA substrates in all the binary and ternary forms are shown in grey (b-e), the dNTPs in c-e are shown as cyan sticks, and the Mg²⁺ ions are represented as magenta spheres in c.
<figure><img src="image_12.png" /><figcaption>**Figure 6.**</figcaption></figure>

Conformational exchange experienced by Pol β. Residues displaying elevated R₂ and ¹⁵N-CPMG relaxation dispersion curves are shown as blue and red spheres, respectively, corresponding to Cα atoms. (a) The apo form shows millisecond dynamics predominantly in the lyase domain (residues 1–90). Residues with *k*<sub>ex</sub> ~1400 s⁻¹ are shown as smaller spheres while residues with *k*<sub>ex</sub> ~4000 s⁻¹ are shown as larger red spheres. Few residues in the other polymerase domains show elevated R₂ values; (b) Only residue Glu21 displays ¹⁵N-CPMG relaxation dispersion in the binary form. Other residues displaying elevated R₂ values are located in the polymerase domain of Pol β. DNA in the binary form is shown in black. The relaxation data used for preparing this figure were taken from reference [8].