Next Article in Journal
Co-Assembled Supported Catalysts: Synthesis of Nano-Structured Supported Catalysts with Hierarchic Pores through Combined Flow and Radiation Induced Co-Assembled Nano-Reactors
Next Article in Special Issue
Solvent-Free Lipase-Catalyzed Synthesis of Technical-Grade Sugar Esters and Evaluation of Their Physicochemical and Bioactive Properties
Previous Article in Journal
Synthesis and Characterization of Ag-Ag2O/TiO2@polypyrrole Heterojunction for Enhanced Photocatalytic Degradation of Methylene Blue
Previous Article in Special Issue
Preparation of Cross-Linked Glucoamylase Aggregates Immobilization by Using Dextrin and Xanthan Gum as Protecting Agents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

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

by
Chitra Narayanan
1,†,
David N. Bernard
1,† and
Nicolas Doucet
1,2,3,*
1
INRS—Institut Armand-Frappier, Université du Québec, 531 Boul. des Prairies, Laval, QC H7V 1B7, Canada
2
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
3
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
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Catalysts 2016, 6(6), 81; https://doi.org/10.3390/catal6060081
Submission received: 5 April 2016 / Revised: 11 May 2016 / Accepted: 20 May 2016 / Published: 27 May 2016
(This article belongs to the Special Issue Enzyme Catalysis)

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.

Graphical Abstract

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 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,6,7,8]. Advances in experimental and computational methodologies continue to offer new insights into enzyme conformational motions over functionally relevant timescales [9,10,11,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,30,31,32,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,37,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 νCPMG, 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 (Rex) 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 νCPMG 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, Δω), in addition to their population dynamics in solution (conformational exchange rates, kex, and equilibrium populations, pA and pB) (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 (S2), conformational exchange parameters (Rex), and the correlation time (τc), offering an interpretation of the internal amplitude of the N–H bond vector in solution [43,44,45]. The order parameter S2 provides an easily interpretable measure of the rigidity of each protein residue, where the S2 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 Rex 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 kcat, 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,53,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,61,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,69,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,73,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,78,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 Rex 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 of 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,92,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 kcat or KM 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 Ki 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 pre-treatment 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 taking 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 Mg2+ ions as cofactors, is shown in Figure 4.
Crystallographic studies have shown that Pol β undergoes large conformational changes upon substrate binding [110,111,112,113,114,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,112,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−1. 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.

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.

References

  1. Garcia-Viloca, M.; Gao, J.; Karplus, M.; Truhlar, D.G. How enzymes work: Analysis by modern rate theory and computer simulations. Science 2004, 303, 186–195. [Google Scholar] [CrossRef] [PubMed]
  2. Wolfenden, R. Degrees of difficulty of water-consuming reactions in the absence of enzymes. Chem. Rev. 2006, 106, 3379–3396. [Google Scholar] [CrossRef] [PubMed]
  3. Hammes, G.G.; Chang, Y.C.; Oas, T.G. Conformational selection or induced fit: A flux description of reaction mechanism. Proc. Natl. Acad. Sci. USA 2009, 106, 13737–13741. [Google Scholar] [CrossRef] [PubMed]
  4. Henzler-Wildman, K.; Kern, D. Dynamic personalities of proteins. Nature 2007, 450, 964–972. [Google Scholar] [CrossRef] [PubMed]
  5. Watt, E.D.; Shimada, H.; Kovrigin, E.L.; Loria, J.P. The mechanism of rate-limiting motions in enzyme function. Proc. Natl. Acad. Sci. USA 2007, 104, 11981–11986. [Google Scholar] [CrossRef] [PubMed]
  6. Venkitakrishnan, R.P.; Zaborowski, E.; McElheny, D.; Benkovic, S.J.; Dyson, H.J.; Wright, P.E. Conformational changes in the active site loops of dihydrofolate reductase during the catalytic cycle. Biochemistry 2004, 43, 16046–16055. [Google Scholar] [CrossRef] [PubMed]
  7. Boehr, D.D. Promiscuity in protein-RNA interactions: Conformational ensembles facilitate molecular recognition in the spliceosome. Bioessays 2012, 34, 174–180. [Google Scholar] [CrossRef] [PubMed]
  8. Berlow, R.B.; Swain, M.; Dalal, S.; Sweasy, J.B.; Loria, J.P. Substrate-dependent millisecond domain motions in DNA polymerase β. J. Mol. Biol. 2012, 419, 171–182. [Google Scholar] [CrossRef] [PubMed]
  9. Kay, L.E. NMR studies of protein structure and dynamics. J. Magn. Reson. 2005, 173, 193–207. [Google Scholar] [CrossRef] [PubMed]
  10. Vallurupalli, P.; Hansen, D.F.; Kay, L.E. Structures of invisible, excited protein states by relaxation dispersion NMR spectroscopy. Proc. Natl. Acad. Sci. USA 2008, 105, 11766–11771. [Google Scholar] [CrossRef] [PubMed]
  11. Ramanathan, A.; Savol, A.; Burger, V.; Chennubhotla, C.S.; Agarwal, P.K. Protein conformational populations and functionally relevant substates. Acc. Chem. Res. 2014, 47, 149–156. [Google Scholar] [CrossRef] [PubMed]
  12. Vanwart, A.T.; Eargle, J.; Luthey-Schulten, Z.; Amaro, R.E. Exploring residue component contributions to dynamical network models of allostery. J. Chem. Theory Comput. 2012, 8, 2949–2961. [Google Scholar] [CrossRef] [PubMed]
  13. Gagné, D.; Charest, L.A.; Morin, S.; Kovrigin, E.L.; Doucet, N. Conservation of flexible residue clusters among structural and functional enzyme homologues. J. Biol. Chem. 2012, 287, 44289–44300. [Google Scholar] [CrossRef] [PubMed]
  14. Eisenmesser, E.Z.; Millet, O.; Labeikovsky, W.; Korzhnev, D.M.; Wolf-Watz, M.; Bosco, D.A.; Skalicky, J.J.; Kay, L.E.; Kern, D. Intrinsic dynamics of an enzyme underlies catalysis. Nature 2005, 438, 117–121. [Google Scholar] [CrossRef] [PubMed]
  15. Agarwal, P.K.; 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. [Google Scholar] [CrossRef] [PubMed]
  16. Goodey, N.M.; Benkovic, S.J. Allosteric regulation and catalysis emerge via a common route. Nat. Chem. Biol. 2008, 4, 474–482. [Google Scholar] [CrossRef] [PubMed]
  17. Suel, G.M.; Lockless, S.W.; Wall, M.A.; Ranganathan, R. Evolutionarily conserved networks of residues mediate allosteric communication in proteins. Nat. Struct. Biol. 2003, 10, 59–69. [Google Scholar] [CrossRef] [PubMed]
  18. Reynolds, K.A.; McLaughlin, R.N.; Ranganathan, R. Hot spots for allosteric regulation on protein surfaces. Cell 2011, 147, 1564–1575. [Google Scholar] [CrossRef] [PubMed]
  19. Cole, R.; Loria, J.P. Evidence for flexibility in the function of ribonuclease A. Biochemistry 2002, 41, 6072–6081. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  21. Pareek, V.; Samanta, M.; Joshi, N.V.; Balaram, H.; Murthy, M.R.; 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. [Google Scholar] [CrossRef] [PubMed]
  22. Torbeev, V.Y.; Raghuraman, H.; Hamelberg, D.; Tonelli, M.; Westler, W.M.; Perozo, E.; Kent, S.B. Protein conformational dynamics in the mechanism of HIV-1 protease catalysis. Proc. Natl. Acad. Sci. USA 2011, 108, 20982–20987. [Google Scholar] [CrossRef] [PubMed]
  23. Kar, G.; Keskin, O.; Gursoy, A.; Nussinov, R. Allostery and population shift in drug discovery. Curr. Opin. Pharmacol. 2010, 10, 715–722. [Google Scholar] [CrossRef] [PubMed]
  24. Agarwal, P.K. Enzymes: An integrated view of structure, dynamics and function. Microb. Cell Fact. 2006, 5, 2. [Google Scholar] [CrossRef] [PubMed]
  25. Kohen, A. Role of dynamics in enzyme catalysis: Substantial versus semantic controversies. Acc. Chem. Res. 2015, 48, 466–473. [Google Scholar] [CrossRef] [PubMed]
  26. Kamerlin, S.C.; Warshel, A. At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis? Proteins 2010, 78, 1339–1375. [Google Scholar] [CrossRef] [PubMed]
  27. Nagel, Z.D.; Klinman, J.P. A 21st century revisionist's view at a turning point in enzymology. Nat. Chem. Biol. 2009, 5, 543–550. [Google Scholar] [CrossRef] [PubMed]
  28. Hanoian, P.; Liu, C.T.; Hammes-Schiffer, S.; Benkovic, S. Perspectives on electrostatics and conformational motions in enzyme catalysis. Acc. Chem. Res. 2015, 48, 482–489. [Google Scholar] [CrossRef] [PubMed]
  29. Grutsch, S.; Bruschweiler, S.; Tollinger, M. NMR methods to study dynamic allostery. PLoS Comput. Biol. 2016, 12, e1004620. [Google Scholar] [CrossRef] [PubMed]
  30. Palmer, A.G., 3rd. Chemical exchange in biomacromolecules: Past, present, and future. J. Magn. Reson. 2014, 241, 3–17. [Google Scholar] [CrossRef] [PubMed]
  31. Kleckner, I.R.; Foster, M.P. An introduction to NMR-based approaches for measuring protein dynamics. Biochim. Biophys. Acta 2011, 1814, 942–968. [Google Scholar] [CrossRef] [PubMed]
  32. Loria, J.P.; Rance, M.; Palmer, A.G., 3rd. A trosy cpmg sequence for characterizing chemical exchange in large proteins. J. Biomol. NMR 1999, 15, 151–155. [Google Scholar] [CrossRef] [PubMed]
  33. Loria, J.P.; Rance, M.; Palmer, A.G. A relaxation-compensated carr-purcell-meiboom-gill sequence for characterizing chemical exchange by NMR spectroscopy. J. Am. Chem. Soc. 1999, 121, 2331–2332. [Google Scholar] [CrossRef]
  34. Oyen, D.; Fenwick, R.B.; Stanfield, R.L.; Dyson, H.J.; Wright, P.E. Cofactor-mediated conformational dynamics promote product release from Escherichia Coli dihydrofolate reductase via an allosteric pathway. J. Am. Chem. Soc. 2015, 137, 9459–9468. [Google Scholar] [CrossRef] [PubMed]
  35. Bhabha, G.; Biel, J.T.; Fraser, J.S. Keep on moving: Discovering and perturbing the conformational dynamics of enzymes. Acc. Chem. Res. 2015, 48, 423–430. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  37. Palmer, A.G., 3rd. NMR probes of molecular dynamics: Overview and comparison with other techniques. Annu. Rev. Biophys. Biomol. Struct. 2001, 30, 129–155. [Google Scholar] [PubMed]
  38. Palmer, A.G., 3rd. NMR characterization of the dynamics of biomacromolecules. Chem. Rev. 2004, 104, 3623–3640. [Google Scholar] [CrossRef] [PubMed]
  39. Manley, G.; Loria, J.P. NMR insights into protein allostery. Arch. Biochem. Biophys. 2012, 519, 223–231. [Google Scholar] [CrossRef] [PubMed]
  40. Luk, L.Y.; Loveridge, E.J.; Allemann, R.K. Protein motions and dynamic effects in enzyme catalysis. Phys. Chem. Chem. Phys. 2015, 17, 30817–30827. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  43. Morin, S. A practical guide to protein dynamics from 15N spin relaxation in solution. Prog. Nucl. Magn. Reson. Spectrosc. 2011, 59, 245–262. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  45. Mandel, A.M.; Akke, M.; Palmer, A.G., 3rd. Backbone dynamics of Escherichia Coli Ribonuclease HI: Correlations with structure and function in an active enzyme. J. Mol. Biol. 1995, 246, 144–163. [Google Scholar] [CrossRef] [PubMed]
  46. Stephenson, J.D.; Kenyon, J.C.; Symmons, M.F.; Lever, A.M. Characterizing 3D RNA structure by single molecule FRET. Methods 2016, in press. [Google Scholar] [CrossRef] [PubMed]
  47. Huang, X.; de Vera, I.M.; Veloro, A.M.; Blackburn, M.E.; Kear, J.L.; Carter, J.D.; Rocca, J.R.; Simmerling, C.; Dunn, B.M.; Fanucci, G.E. 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. [Google Scholar] [CrossRef] [PubMed]
  48. Huang, X.; Britto, M.D.; Kear-Scott, J.L.; Boone, C.D.; Rocca, J.R.; Simmerling, C.; McKenna, R.; Bieri, M.; Gooley, P.R.; Dunn, B.M.; et al. The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics. J. Biol. Chem. 2014, 289, 17203–17214. [Google Scholar] [CrossRef] [PubMed]
  49. Jeschke, G.; Polyhach, Y. Distance measurements on spin-labelled biomacromolecules by pulsed electron paramagnetic resonance. Phys. Chem. Chem. Phys. 2007, 9, 1895–1910. [Google Scholar] [CrossRef] [PubMed]
  50. Boehr, D.D.; McElheny, D.; Dyson, H.J.; Wright, P.E. Millisecond timescale fluctuations in dihydrofolate reductase are exquisitely sensitive to the bound ligands. Proc. Natl. Acad. Sci. USA 2010, 107, 1373–1378. [Google Scholar] [CrossRef] [PubMed]
  51. Boehr, D.D.; Liu, X.; Yang, X. Targeting structural dynamics of the RNA-dependent RNA polymerase for anti-viral strategies. Curr. Opin. Virol. 2014, 9, 194–200. [Google Scholar] [CrossRef] [PubMed]
  52. Voelz, V.A.; Bowman, G.R.; Beauchamp, K.; Pande, V.S. Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1–39). J. Am. Chem. Soc. 2010, 132, 1526–1528. [Google Scholar] [CrossRef] [PubMed]
  53. Bowman, G.R.; Voelz, V.A.; Pande, V.S. Atomistic folding simulations of the five-helix bundle protein λ6–85. J. Am. Chem. Soc. 2011, 133, 664–667. [Google Scholar] [CrossRef] [PubMed]
  54. Voelz, V.A.; Jager, M.; Yao, S.; Chen, Y.; Zhu, L.; Waldauer, S.A.; Bowman, G.R.; Friedrichs, M.; Bakajin, O.; Lapidus, L.J.; 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. [Google Scholar] [CrossRef] [PubMed]
  55. Shaw, D.E.; Maragakis, P.; Lindorff-Larsen, K.; Piana, S.; Dror, R.O.; Eastwood, M.P.; Bank, J.A.; Jumper, J.M.; Salmon, J.K.; Shan, Y.; et al. Atomic-level characterization of the structural dynamics of proteins. Science 2010, 330, 341–346. [Google Scholar] [CrossRef] [PubMed]
  56. Bucher, D.; Grant, B.J.; Markwick, P.R.; McCammon, J.A. Accessing a hidden conformation of the maltose binding protein using accelerated molecular dynamics. PLoS Comput. Biol. 2011, 7, e1002034. [Google Scholar] [CrossRef] [PubMed]
  57. Hamelberg, D.; Mongan, J.; McCammon, J.A. Accelerated molecular dynamics: A promising and efficient simulation method for biomolecules. J. Chem. Phys. 2004, 120, 11919–11929. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  59. Da, L.T.; 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. [Google Scholar] [CrossRef] [PubMed]
  60. Tagle, A.M. Percutaneous encoscopic gastrostomy. Rev. Gastroenterol. Peru 1998, 18, 56–61. [Google Scholar]
  61. Sutto, L.; Gervasio, F.L. Effects of oncogenic mutations on the conformational free-energy landscape of EGFR kinase. Proc. Natl. Acad. Sci. USA 2013, 110, 10616–10621. [Google Scholar] [CrossRef] [PubMed]
  62. Abrams, C.; Bussi, G. Enhanced sampling in molecular dynamics using metadynamics, replica-exchange, and temperature-acceleration. Entropy 2014, 16, 163–199. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  64. Doshi, U.; McGowan, L.C.; Ladani, S.T.; Hamelberg, D. Resolving the complex role of enzyme conformational dynamics in catalytic function. Proc. Natl. Acad. Sci. USA 2012, 109, 5699–5704. [Google Scholar] [CrossRef] [PubMed]
  65. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in monte carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  67. Andersen, O.J.; Grouleff, J.; Needham, P.; Walker, R.C.; 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. [Google Scholar] [CrossRef] [PubMed]
  68. Singhal, N.; Snow, C.D.; Pande, V.S. 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. [Google Scholar] [CrossRef] [PubMed]
  69. Bowman, G.R.; Beauchamp, K.A.; Boxer, G.; Pande, V.S. Progress and challenges in the automated construction of Markov state models for full protein systems. J. Chem. Phys. 2009, 131, 124101. [Google Scholar] [CrossRef] [PubMed]
  70. Noé, F.; Fischer, S. Transition networks for modeling the kinetics of conformational change in macromolecules. Curr. Opin. Struct. Biol. 2008, 18, 154–162. [Google Scholar] [CrossRef] [PubMed]
  71. Ramanathan, A.; Agarwal, P.K. Evolutionarily conserved linkage between enzyme fold, flexibility, and catalysis. PLoS Biol. 2011, 9, e1001193. [Google Scholar] [CrossRef] [PubMed]
  72. Shukla, D.; Hernandez, C.X.; Weber, J.K.; Pande, V.S. Markov state models provide insights into dynamic modulation of protein function. Acc. Chem. Res. 2015, 48, 414–422. [Google Scholar] [CrossRef] [PubMed]
  73. Chodera, J.D.; Noe, F. Markov state models of biomolecular conformational dynamics. Curr. Opin. Struct. Biol. 2014, 25, 135–144. [Google Scholar] [CrossRef] [PubMed]
  74. Morcos, F.; Chatterjee, S.; McClendon, C.L.; Brenner, P.R.; Lopez-Rendon, R.; Zintsmaster, J.; Ercsey-Ravasz, M.; Sweet, C.R.; Jacobson, M.P.; Peng, J.W.; et al. Modeling conformational ensembles of slow functional motions in Pin1-WW. PLoS Comput. Biol. 2010, 6, e1001015. [Google Scholar] [CrossRef] [PubMed]
  75. Ramanathan, A.; Savol, A.J.; Langmead, C.J.; Agarwal, P.K.; Chennubhotla, C.S. Discovering conformational sub-states relevant to protein function. PLoS ONE 2011, 6, e15827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Gagné, D.; French, R.L.; Narayanan, C.; Simonovic, M.; Agarwal, P.K.; Doucet, N. Perturbation of the conformational eynamics of an active-site loop alters enzyme activity. Structure 2015, 23, 2256–2266. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  78. De Simone, A.; Montalvao, R.W.; Dobson, C.M.; 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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  80. Sundquist, W.I.; Krausslich, H.G. HIV-1 assembly, budding, and maturation. Cold Spring Harb. Perspect. Med. 2012, 2, a006924. [Google Scholar] [CrossRef] [PubMed]
  81. Ghosh, A.K.; Osswald, H.L.; Prato, G. Recent progress in the development of HIV-1 protease inhibitors for the treatment of HIV/AIDS. J. Med. Chem. 2016, in press. [Google Scholar] [CrossRef] [PubMed]
  82. Spinelli, S.; Liu, Q.Z.; Alzari, P.M.; Hirel, P.H.; Poljak, R.J. The three-dimensional structure of the aspartyl protease from the HIV-1 isolate BRU. Biochimie 1991, 73, 1391–1396. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [PubMed]
  84. Nicholson, L.K.; Yamazaki, T.; Torchia, D.A.; Grzesiek, S.; Bax, A.; Stahl, S.J.; Kaufman, J.D.; Wingfield, P.T.; Lam, P.Y.; Jadhav, P.K.; et al. Flexibility and function in HIV-1 protease. Nat. Struct. Biol. 1995, 2, 274–280. [Google Scholar] [CrossRef] [PubMed]
  85. Ishima, R.; Louis, J.M. A diverse view of protein dynamics from NMR studies of HIV-1 protease flaps. Proteins 2008, 70, 1408–1415. [Google Scholar] [CrossRef] [PubMed]
  86. Xia, J.; Deng, N.J.; Levy, R.M. 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. [Google Scholar] [CrossRef] [PubMed]
  87. Ishima, R.; Torchia, D.A.; Lynch, S.M.; Gronenborn, A.M.; Louis, J.M. 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. [Google Scholar] [CrossRef] [PubMed]
  88. Deng, N.J.; Zheng, W.; Gallicchio, E.; Levy, R.M. Insights into the dynamics of HIV-1 protease: A kinetic network model constructed from atomistic simulations. J. Am. Chem. Soc. 2011, 133, 9387–9394. [Google Scholar] [CrossRef] [PubMed]
  89. Shukla, D.; Meng, Y.; Roux, B.; Pande, V.S. Activation pathway of SRC kinase reveals intermediate states as targets for drug design. Nat. Commun. 2014, 5, 3397. [Google Scholar] [CrossRef] [PubMed]
  90. Roche, J.; Louis, J.M.; Bax, A. Conformation of inhibitor-free HIV-1 protease derived from nmr spectroscopy in a weakly oriented solution. Chembiochem 2015, 16, 214–218. [Google Scholar] [CrossRef] [PubMed]
  91. Yamazaki, T.; Hinck, A.P.; Wang, Y.X.; Nicholson, L.K.; Torchia, D.A.; Wingfield, P.; Stahl, S.J.; Kaufman, J.D.; Chang, C.H.; Domaille, P.J.; 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. [Google Scholar] [CrossRef] [PubMed]
  92. Leonis, G.; Steinbrecher, T.; Papadopoulos, M.G. 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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
  94. Mittal, S.; Cai, Y.; Nalam, M.N.; Bolon, D.N.; Schiffer, C.A. Hydrophobic core flexibility modulates enzyme activity in HIV-1 protease. J. Am. Chem. Soc. 2012, 134, 4163–4168. [Google Scholar] [CrossRef] [PubMed]
  95. Goldfarb, N.E.; Ohanessian, M.; Biswas, S.; McGee, T.D., Jr.; Mahon, B.P.; Ostrov, D.A.; 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. [Google Scholar] [CrossRef] [PubMed]
  96. Lee, J.; Goodey, N.M. Catalytic contributions from remote regions of enzyme structure. Chem. Rev. 2011, 111, 7595–7624. [Google Scholar] [CrossRef] [PubMed]
  97. Wu, T.D.; Schiffer, C.A.; Gonzales, M.J.; Taylor, J.; Kantor, R.; Chou, S.; Israelski, D.; Zolopa, A.R.; Fessel, W.J.; Shafer, R.W. Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments. J. Virol. 2003, 77, 4836–4847. [Google Scholar] [CrossRef] [PubMed]
  98. Ragland, D.A.; Nalivaika, E.A.; Nalam, M.N.; Prachanronarong, K.L.; Cao, H.; Bandaranayake, R.M.; Cai, Y.; Kurt-Yilmaz, N.; Schiffer, C.A. 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. [Google Scholar] [CrossRef] [PubMed]
  99. Tzoupis, H.; Leonis, G.; Avramopoulos, A.; Mavromoustakos, T.; Papadopoulos, M.G. 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. [Google Scholar] [CrossRef] [PubMed]
  100. Carter, J.D.; Gonzales, E.G.; Huang, X.; Smith, A.N.; de Vera, I.M.; D'Amore, P.W.; Rocca, J.R.; Goodenow, M.M.; Dunn, B.M.; Fanucci, G.E. Effects of pre and post therapy drug-pressure selected mutations on HIV-1 protease conformational sampling. FEBS Lett. 2014, 588, 3123–3128. [Google Scholar] [CrossRef] [PubMed]
  101. Kunze, J.; Todoroff, N.; Schneider, P.; Rodrigues, T.; Geppert, T.; Reisen, F.; Schreuder, H.; Saas, J.; Hessler, G.; Baringhaus, K.H.; 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. [Google Scholar] [CrossRef] [PubMed]
  102. Ung, P.M.; Dunbar, J.B., Jr.; Gestwicki, J.E.; Carlson, H.A. An allosteric modulator of HIV-1 protease shows equipotent inhibition of wild-type and drug-resistant proteases. J. Med. Chem. 2014, 57, 6468–6478. [Google Scholar] [CrossRef] [PubMed]
  103. Beard, W.A.; Wilson, S.H. Structure and mechanism of DNA polymerase β. Chem. Rev. 2006, 106, 361–382. [Google Scholar] [CrossRef] [PubMed]
  104. Barnes, D.E.; Lindahl, T. Repair and genetic consequences of endogenous DNA base damage in mammalian cells. Annu. Rev. Genet. 2004, 38, 445–476. [Google Scholar] [CrossRef] [PubMed]
  105. Starcevic, D.; Dalal, S.; Sweasy, J.B. Is there a link between DNA polymerase β and cancer? Cell Cycle 2004, 3, 998–1001. [Google Scholar] [CrossRef] [PubMed]
  106. Loeb, L.A.; Monnat, R.J., Jr. DNA polymerases and human disease. Nat. Rev. Genet. 2008, 9, 594–604. [Google Scholar] [CrossRef] [PubMed]
  107. Beard, W.A.; Shock, D.D.; Yang, X.P.; DeLauder, S.F.; Wilson, S.H. Loss of DNA polymerase β stacking interactions with templating purines, but not pyrimidines, alters catalytic efficiency and fidelity. J. Biol. Chem. 2002, 277, 8235–8242. [Google Scholar] [CrossRef] [PubMed]
  108. Towle-Weicksel, J.B.; Dalal, S.; Sohl, C.D.; Doublie, S.; Anderson, K.S.; Sweasy, J.B. 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. [Google Scholar] [CrossRef] [PubMed]
  109. Moscato, B.; Swain, M.; Loria, J.P. Induced fit in the selection of correct versus incorrect nucleotides by DNA polymerase β. Biochemistry 2016, 55, 382–395. [Google Scholar] [CrossRef] [PubMed]
  110. Sawaya, M.R.; Pelletier, H.; Kumar, A.; Wilson, S.H.; Kraut, J. Crystal structure of rat DNA polymerase β: Evidence for a common polymerase mechanism. Science 1994, 264, 1930–1935. [Google Scholar] [CrossRef] [PubMed]
  111. Sawaya, M.R.; Prasad, R.; Wilson, S.H.; 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. [Google Scholar] [CrossRef] [PubMed]
  112. Batra, V.K.; Beard, W.A.; Shock, D.D.; Krahn, J.M.; Pedersen, L.C.; Wilson, S.H. Magnesium-induced assembly of a complete DNA polymerase catalytic complex. Structure 2006, 14, 757–766. [Google Scholar] [CrossRef] [PubMed]
  113. Batra, V.K.; Beard, W.A.; Shock, D.D.; Pedersen, L.C.; Wilson, S.H. Structures of DNA polymerase β with active-site mismatches suggest a transient abasic site intermediate during misincorporation. Mol. Cell. 2008, 30, 315–324. [Google Scholar] [CrossRef] [PubMed]
  114. Freudenthal, B.D.; Beard, W.A.; Wilson, S.H. Structures of dNTP intermediate states during DNA polymerase active site assembly. Structure 2012, 20, 1829–1837. [Google Scholar] [CrossRef] [PubMed]
  115. Freudenthal, B.D.; Beard, W.A.; Shock, D.D.; Wilson, S.H. Observing a DNA polymerase choose right from wrong. Cell 2013, 154, 157–168. [Google Scholar] [CrossRef] [PubMed]
  116. Donigan, K.A.; Sun, K.W.; Nemec, A.A.; Murphy, D.L.; Cong, X.; Northrup, V.; Zelterman, D.; Sweasy, J.B. Human POLB gene is mutated in high percentage of colorectal tumors. J. Biol. Chem. 2012, 287, 23830–23839. [Google Scholar] [CrossRef] [PubMed]
  117. Beard, W.A.; Shock, D.D.; Batra, V.K.; Prasad, R.; Wilson, S.H. Substrate-induced DNA polymerase β activation. J. Biol. Chem. 2014, 289, 31411–31422. [Google Scholar] [CrossRef] [PubMed]
  118. Bakhtina, M.; Lee, S.; Wang, Y.; Dunlap, C.; Lamarche, B.; Tsai, M.D. 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. [Google Scholar] [CrossRef] [PubMed]
  119. Arora, K.; Beard, W.A.; Wilson, S.H.; Schlick, T. Mismatch-induced conformational distortions in polymerase β support an induced-fit mechanism for fidelity. Biochemistry 2005, 44, 13328–13341. [Google Scholar] [CrossRef] [PubMed]
  120. Kirmizialtin, S.; Nguyen, V.; Johnson, K.A.; Elber, R. How conformational dynamics of DNA polymerase select correct substrates: Experiments and simulations. Structure 2012, 20, 618–627. [Google Scholar] [CrossRef] [PubMed]
  121. Zhong, X.; Patel, S.S.; Werneburg, B.G.; Tsai, M.D. DNA polymerase β: Multiple conformational changes in the mechanism of catalysis. Biochemistry 1997, 36, 11891–11900. [Google Scholar] [CrossRef] [PubMed]
  122. Balbo, P.B.; Wang, E.C.; Tsai, M.D. Kinetic mechanism of active site assembly and chemical catalysis of DNA polymerase β. Biochemistry 2011, 50, 9865–9875. [Google Scholar] [CrossRef] [PubMed]
  123. Williamson, M.P. Using chemical shift perturbation to characterise ligand binding. Prog. Nucl. Magn. 2013, 73, 1–16. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Conformational exchange experienced by the Met20 loop in E. coli DHFR, as probed by 15N-CPMG NMR relaxation dispersion experiments. The CPMG experiment is particularly well suited to extract exchange rates (kex) 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 (kcat). This method can also provide quantitative information on low-populated “invisible” excited sub-states and their populations in solution (pA and pB). 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 1. Conformational exchange experienced by the Met20 loop in E. coli DHFR, as probed by 15N-CPMG NMR relaxation dispersion experiments. The CPMG experiment is particularly well suited to extract exchange rates (kex) 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 (kcat). This method can also provide quantitative information on low-populated “invisible” excited sub-states and their populations in solution (pA and pB). 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).
Catalysts 06 00081 g001
Figure 2. 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 2. 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.
Catalysts 06 00081 g002
Figure 3. 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 3. 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.
Catalysts 06 00081 g003
Figure 4. Schematic representation of the catalytic mechanism of Pol β. Eo and Ec correspond to the open and closed conformations of the enzyme. Sequential binding of the DNA (DNAn) and dNTP leads to the formation of the binary (step 1) and ternary complexes (step 2), respectively. Subsequent conformational rearrangement of the enzyme (Ec) aligns the active site for catalysis (step 3). The chemical step of catalysis is followed by the release of pyrophosphate (PPi) and the DNA product (DNAn+1) (step 4). Figure adapted from [108,109].
Figure 4. Schematic representation of the catalytic mechanism of Pol β. Eo and Ec correspond to the open and closed conformations of the enzyme. Sequential binding of the DNA (DNAn) and dNTP leads to the formation of the binary (step 1) and ternary complexes (step 2), respectively. Subsequent conformational rearrangement of the enzyme (Ec) aligns the active site for catalysis (step 3). The chemical step of catalysis is followed by the release of pyrophosphate (PPi) and the DNA product (DNAn+1) (step 4). Figure adapted from [108,109].
Catalysts 06 00081 g004
Figure 5. 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 (PDB 3C2M) 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 ce, illustrating subtle yet functionally significant conformational exchange in this enzyme. DNA substrates in all the binary and ternary forms are shown in grey (be), the dNTPs in ce are shown as cyan sticks, and the Mg2+ ions are represented as magenta spheres in c.
Figure 5. 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 (PDB 3C2M) 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 ce, illustrating subtle yet functionally significant conformational exchange in this enzyme. DNA substrates in all the binary and ternary forms are shown in grey (be), the dNTPs in ce are shown as cyan sticks, and the Mg2+ ions are represented as magenta spheres in c.
Catalysts 06 00081 g005
Figure 6. Conformational exchange experienced by Pol β. Residues displaying elevated R2 and 15N-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 kex ~1400 s−1 are shown as smaller spheres whiles residues with kex ~4000 s−1 are shown as larger red spheres. Few residues in the other polymerase domains show elevated R2 values; (b) Only residue Glu21 displays 15N-CPMG relaxation dispersion in the binary form. Other residues displaying elevated R2 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].
Figure 6. Conformational exchange experienced by Pol β. Residues displaying elevated R2 and 15N-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 kex ~1400 s−1 are shown as smaller spheres whiles residues with kex ~4000 s−1 are shown as larger red spheres. Few residues in the other polymerase domains show elevated R2 values; (b) Only residue Glu21 displays 15N-CPMG relaxation dispersion in the binary form. Other residues displaying elevated R2 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].
Catalysts 06 00081 g006

Share and Cite

MDPI and ACS Style

Narayanan, C.; Bernard, D.N.; Doucet, N. Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies. Catalysts 2016, 6, 81. https://doi.org/10.3390/catal6060081

AMA Style

Narayanan C, Bernard DN, Doucet N. Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies. Catalysts. 2016; 6(6):81. https://doi.org/10.3390/catal6060081

Chicago/Turabian Style

Narayanan, Chitra, David N. Bernard, and Nicolas Doucet. 2016. "Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies" Catalysts 6, no. 6: 81. https://doi.org/10.3390/catal6060081

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop