Enzyme Dynamics During Catalysis Measured by NMR Spectroscopy

Conformational dynamics of adenylate kinase in crystals

Adenylate kinase is a ubiquitous enzyme in living systems and undergoes dramatic conformational changes during its catalytic cycle. For these reasons, it is widely studied by genetic, biochemical, and biophysical methods, both experimental and theoretical. We have determined the basic crystal structures of three differently liganded states of adenylate kinase from Methanotorrus igneus, a hyperthermophilic organism whose adenylate kinase is a homotrimeric oligomer. The multiple copies of each protomer in the asymmetric unit of the crystal provide a unique opportunity to study the variation in the structure and were further analyzed using advanced crystallographic refinement methods and analysis tools to reveal conformational heterogeneity and, thus, implied dynamic behaviors in the catalytic cycle.

Biomolecular NMR spectroscopy in the era of artificial intelligence

Biomolecular nuclear magnetic resonance (NMR) spectroscopy and artificial intelligence (AI) have a burgeoning synergy. Deep learning-based structural predictors have forever changed structural biology, yet these tools currently face limitations in accurately characterizing protein dynamics, allostery, and conformational heterogeneity. We begin by highlighting the unique abilities of biomolecular NMR spectroscopy to complement AI-based structural predictions toward addressing these knowledge gaps. We then highlight the direct integration of deep learning approaches into biomolecular NMR methods. AI-based tools can dramatically improve the acquisition and analysis of NMR spectra, enhancing the accuracy and reliability of NMR measurements, thus streamlining experimental processes. Additionally, deep learning enables the development of novel types of NMR experiments that were previously unattainable, expanding the scope and potential of biomolecular NMR spectroscopy. Ultimately, a combination of AI and NMR promises to further revolutionize structural biology on several levels, advance our understanding of complex biomolecular systems, and accelerate drug discovery efforts.

New experimental evidence for pervasive dynamics in proteins

There is ample computational, but only sparse experimental data suggesting that pico-ns motions with 1 Å amplitude are pervasive in proteins in solution. Such motions, if present in reality, must deeply affect protein function and protein entropy. Several NMR relaxation experiments have provided insights into motions of proteins in solution, but they primarily report on azimuthal angle variations of vectors of covalently-linked atoms. As such, these measurements are not sensitive to distance fluctuations, and cannot but under-represent the dynamical properties of proteins. Here we analyze a novel NMR relaxation experiment to measure amide proton transverse relaxation rates in uniformly 15 N labeled proteins, and present results for protein domain GB1 at 283 and 303 K. These relaxation rates depend on fluctuations of dipolar interactions between 1 HN and many nearby protons on both the backbone and sidechains. Importantly, they also report on fluctuations in the distances between these protons. We obtained a large mismatch between rates computed from the crystal structure of GB1 and the experimental rates. But when the relaxation rates were calculated from a 200 ns molecular dynamics trajectory using a novel program suite, we obtained a substantial improvement in the correspondence of experimental and theoretical rates. As such, this work provides novel experimental evidence of widespread motions in proteins. Since the improvements are substantial, but not sufficient, this approach may also present a new benchmark to help improve the theoretical forcefields underlying the molecular dynamics calculations.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Significant Loop Motions in the SsoPTP Protein Tyrosine Phosphatase Allow for Dual General Acid Functionality

Conformational dynamics are important factors in the function of enzymes, including protein tyrosine phosphatases (PTPs). Crystal structures of PTPs first revealed the motion of a protein loop bearing a conserved catalytic aspartic acid, and subsequent nuclear magnetic resonance and computational analyses have shown the presence of motions, involved in catalysis and allostery, within and beyond the active site. The tyrosine phosphatase from the thermophilic and acidophilic Sulfolobus solfataricus (SsoPTP) displays motions of its acid loop together with dynamics of its phosphoryl-binding P-loop and the Q-loop, the first instance of such motions in a PTP. All three loops share the same exchange rate, implying their motions are coupled. Further evidence of conformational flexibility comes from mutagenesis, kinetics, and isotope effect data showing that E40 can function as an alternate general acid to protonate the leaving group when the conserved acid, D69, is mutated to asparagine. SsoPTP is not the first PTP to exhibit an alternate general acid (after VHZ and TkPTP), but E40 does not correspond to the sequence or structural location of the alternate general acids in those precedents. A high-resolution X-ray structure with the transition state analogue vanadate clarifies the role of the active site arginine R102, which varied in structures of substrates bound to a catalytically inactive mutant. The coordinated motions of all three functional loops in SsoPTP, together with the function of an alternate general acid, suggest that catalytically competent conformations are present in solution that have not yet been observed in crystal structures.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

Biogenesis, quality control, and structural dynamics of proteins as explored in living cells via site-directed photocrosslinking

Protein biogenesis and quality control are essential to maintaining a functional pool of proteins and involve numerous protein factors that dynamically and transiently interact with each other and with the substrate proteins in living cells. Conventional methods are hardly effective for studying dynamic, transient, and weak protein-protein interactions that occur in cells. Herein, we review how the site-directed photocrosslinking approach, which relies on the genetic incorporation of a photoreactive unnatural amino acid into a protein of interest at selected individual amino acid residue positions and the covalent trapping of the interacting proteins upon ultraviolent irradiation, has become a highly efficient way to explore the aspects of protein contacts in living cells. For example, in the past decade, this approach has allowed the profiling of the in vivo substrate proteins of chaperones or proteases under both physiologically optimal and stressful (e.g., acidic) conditions, mapping residues located at protein interfaces, identifying new protein factors involved in the biogenesis of membrane proteins, trapping transiently formed protein complexes, and snapshotting different structural states of a protein. We anticipate that the site-directed photocrosslinking approach will play a fundamental role in dissecting the detailed mechanisms of protein biogenesis, quality control, and dynamics in the future.

Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase

The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.

Surface engineering of polyester-degrading enzymes to improve efficiency and tune specificity

Certain members of the carboxylesterase superfamily can act at the interface between water and water-insoluble substrates. However, nonnatural bulky polyesters usually are not efficiently hydrolyzed. In the recent years, the potential of enzyme engineering to improve hydrolysis of synthetic polyesters has been demonstrated. Regions on the enzyme surface have been modified by using site-directed mutagenesis in order to tune sorption processes through increased hydrophobicity of the enzyme surface. Such modifications can involve specific amino acid substitutions, addition of binding modules, or truncation of entire domains improving sorption properties and/or dynamics of the enzyme. In this review, we provide a comprehensive overview on different strategies developed in the recent years for enzyme surface engineering to improve the activity of polyester-hydrolyzing enzymes.

The Exact Nuclear Overhauser Enhancement: Recent Advances

Although often depicted as rigid structures, proteins are highly dynamic systems, whose motions are essential to their functions. Despite this, it is difficult to investigate protein dynamics due to the rapid timescale at which they sample their conformational space, leading most NMR-determined structures to represent only an averaged snapshot of the dynamic picture. While NMR relaxation measurements can help to determine local dynamics, it is difficult to detect translational or concerted motion, and only recently have significant advances been made to make it possible to acquire a more holistic representation of the dynamics and structural landscapes of proteins. Here, we briefly revisit our most recent progress in the theory and use of exact nuclear Overhauser enhancements (eNOEs) for the calculation of structural ensembles that describe their conformational space. New developments are primarily targeted at increasing the number and improving the quality of extracted eNOE distance restraints, such that the multi-state structure calculation can be applied to proteins of higher molecular weights. We then review the implications of the exact NOE to the protein dynamics and function of cyclophilin A and the WW domain of Pin1, and finally discuss our current research and future directions.

Protein NMR Studies of Substrate Binding to Human Blood Group A and B Glycosyltransferases

Donor and acceptor substrate binding to human blood group A and B glycosyltransferases (GTA, GTB) has been studied by a variety of protein NMR experiments. Prior crystallographic studies had shown these enzymes to adopt an open conformation in the absence of substrates. Binding either of the donor substrate UDP-Gal or of UDP induces a semiclosed conformation. In the presence of both donor and acceptor substrates, the enzymes shift towards a closed conformation with ordering of an internal loop and the C-terminal residues, which then completely cover the donor-binding pocket. Chemical-shift titrations of uniformly ² H,¹⁵ N-labeled GTA or GTB with UDP affected about 20 % of all crosspeaks in ¹ H,¹⁵ N TROSY-HSQC spectra, reflecting substantial plasticity of the enzymes. On the other hand, it is this conformational flexibility that impedes NH backbone assignments. Chemical-shift-perturbation experiments with δ1-[¹³ C]methyl-Ile-labeled samples revealed two Ile residues-Ile123 at the bottom of the UDP binding pocket, and Ile192 as part of the internal loop-that were significantly disturbed upon stepwise addition of UDP and H-disaccharide, also revealing long-range perturbations. Finally, methyl TROSY-based relaxation dispersion experiments do not reveal micro- to millisecond timescale motions. Although this study reveals substantial conformational plasticity of GTA and GTB, the matter of how binding of substrates shifts the enzymes into catalytically competent states remains enigmatic.

Perspective: Defining and quantifying the role of dynamics in enzyme catalysis

Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme's conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.

Networks of Dynamic Allostery Regulate Enzyme Function

Many protein systems rely on coupled dynamic networks to allosterically regulate function. However, the broad conformational space sampled by non-coherently dynamic systems has precluded detailed analysis of their communication mechanisms. Here, we have developed a methodology that combines the high sensitivity afforded by nuclear magnetic resonance relaxation techniques and single-site multiple mutations, termed RASSMM, to identify two allosterically coupled dynamic networks within the non-coherently dynamic enzyme cyclophilin A. Using this methodology, we discovered two key hotspot residues, Val6 and Val29, that communicate through these networks, the mutation of which altered active-site dynamics, modulating enzymatic turnover of multiple substrates. Finally, we utilized molecular dynamics simulations to identify the mechanism by which one of these hotspots is coupled to the larger dynamic networks. These studies confirm a link between enzyme dynamics and the catalytic cycle of cyclophilin A and demonstrate how dynamic allostery may be engineered to tune enzyme function.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Determination of the Full Catalytic Cycle among Multiple Cyclophilin Family Members and Limitations on the Application of CPMG-RD in Reversible Catalytic Systems

Cyclophilins catalyze cis ↔ trans isomerization of peptidyl-prolyl bonds, influencing protein folding along with a breadth of other biological functions such as signal transduction. Here, we have determined the microscopic rate constants defining the full enzymatic cycle for three human cyclophilins and a more distantly related thermophilic bacterial cyclophilin when catalyzing interconversion of a biologically representative peptide substrate. The cyclophilins studied here exhibit variability in on-enzyme interconversion as well as an up to 2-fold range in rates of substrate binding and release. However, among the human cyclophilins, the microscopic rate constants appear to have been tuned to maintain remarkably similar isomerization rates without a concurrent conservation of apparent binding affinities. While the structures and active site compositions of the human cyclophilins studied here are highly conserved, we find that the enzymes exhibit significant variability in microsecond to millisecond time scale mobility, suggesting a role for the inherent conformational fluctuations that exist within the cyclophilin family as being functionally relevant in regulating substrate interactions. We have additionally modeled the relaxation dispersion profile given by the commonly employed Carr-Purcell-Meiboom-Gill relaxation dispersion (CPMG-RD) experiment when applied to a reversible enzymatic system such as cyclophilin isomerization and identified a significant limitation in the applicability of this approach for monitoring on-enzyme turnover. Specifically, we show both computationally and experimentally that the CPMG-RD experiment is sensitive to noncatalyzed substrate binding and release in reversible systems even at saturating substrate concentrations unless the on-enzyme interconversion rate is much faster than the substrate release rate.

Structure and Dynamics of GeoCyp: A Thermophilic Cyclophilin with a Novel Substrate Binding Mechanism That Functions Efficiently at Low Temperatures

Thermophilic proteins have found extensive use in research and industrial applications because of their high stability and functionality at elevated temperatures while simultaneously providing valuable insight into our understanding of protein folding, stability, dynamics, and function. Cyclophilins, constituting a ubiquitously expressed family of peptidyl-prolyl isomerases with a range of biological functions and disease associations, have been utilized both for conferring stress tolerances and in exploring the link between conformational dynamics and enzymatic function. To date, however, no active thermophilic cyclophilin has been fully biophysically characterized. Here, we determine the structure of a thermophilic cyclophilin (GeoCyp) from Geobacillus kaustophilus, characterize its dynamic motions over several time scales using an array of methodologies that include chemical shift-based methods and relaxation experiments over a range of temperatures, and measure catalytic activity over a range of temperatures to compare its structure, dynamics, and function to those of a mesophilic counterpart, human cyclophilin A (CypA). Unlike those of most thermophile/mesophile pairs, GeoCyp catalysis is not substantially impaired at low temperatures as compared to that of CypA, retaining ~70% of the activity of its mesophilic counterpart. Examination of substrate-bound ensembles reveals a mechanism by which the two cyclophilins may have adapted to their environments through altering dynamic loop motions and a critical residue that acts as a clamp to regulate substrate binding differentially in CypA and GeoCyp. Fast time scale (pico- to nanosecond) dynamics are largely conserved between the two proteins, in accordance with the high degree of structural similarity, although differences do exist in their temperature dependencies. Slower (microsecond) time scale motions are likewise localized to similar regions in the two proteins with some variability in their magnitudes yet do not exhibit significant temperature dependencies in either enzyme.

Enzyme dynamics from NMR spectroscopy

Biological activities of enzymes, including regulation or coordination of mechanistic stages preceding or following the chemical step, may depend upon kinetic or equilibrium changes in protein conformations. Exchange of more open or flexible conformational states with more closed or constrained states can influence inhibition, allosteric regulation, substrate recognition, formation of the Michaelis complex, side reactions, and product release. NMR spectroscopy has long been applied to the study of conformational dynamic processes in enzymes because these phenomena can be characterized over multiple time scales with atomic site resolution. Laboratory-frame spin-relaxation measurements, sensitive to reorientational motions on picosecond–nanosecond time scales, and rotating-frame relaxation-dispersion measurements, sensitive to chemical exchange processes on microsecond–millisecond time scales, provide information on both conformational distributions and kinetics. This Account reviews NMR spin relaxation studies of the enzymes ribonuclease HI from mesophilic (Escherichia coli) and thermophilic (Thermus thermophilus) bacteria, E. coli AlkB, and Saccharomyces cerevisiae triosephosphate isomerase to illustrate the contributions of conformational flexibility and dynamics to diverse steps in enzyme mechanism.

Spin relaxation measurements and molecular dynamics (MD) simulations of the bacterial ribonuclease H enzymes show that the handle region, one of three loop regions that interact with substrates, interconverts between two conformations. Comparison of these conformations with the structure of the complex between Homo sapiens ribonuclease H and a DNA:RNA substrate suggests that the more closed state is inhibitory to binding. The large population of the closed conformation in T. thermophilus ribonuclease H contributes to the increased Michaelis constant compared with the E. coli enzyme. NMR spin relaxation and fluorescence spectroscopy have characterized a conformational transition in AlkB between an open state, in which the side chains of methionine residues in the active site are disordered, and a closed state, in which these residues are ordered. The open state is highly populated in the AlkB/Zn(II) complex, and the closed state is highly populated in the AlkB/Zn(II)/2OG/substrate complex, in which 2OG is the 2-oxoglutarate cosubstrate and the substrate is a methylated DNA oligonucleotide. The equilibrium is shifted to approximately equal populations of the two conformations in the AlkB/Zn(II)/2OG complex. The conformational shift induced by 2OG ensures that 2OG binds to AlkB/Zn(II) prior to the substrate. In addition, the opening rate of the closed conformation limits premature release of substrate, preventing generation of toxic side products by reaction with water. Closure of active site loop 6 in triosephosphate isomerase is critical for forming the Michaelis complex, but reopening of the loop after the reaction is (partially) rate limiting. NMR spin relaxation and MD simulations of triosephosphate isomerase in complex with glycerol 3-phosphate demonstrate that closure of loop 6 is a highly correlated rigid-body motion. The MD simulations also indicate that motions of Gly173 in the most flexible region of loop 6 contribute to opening of the active site loop for product release. Considered together, these three enzyme systems illustrate the power of NMR spin relaxation investigations in providing global insights into the role of conformational dynamic processes in the mechanisms of enzymes from initial activation to final product release.

Exploring the multiscale signaling behavior of phototropin1 from Chlamydomonas reinhardtii using a full-residue space kinetic Monte Carlo molecular dynamics technique

Devising analysis tools for elucidating the regulatory mechanism of complex enzymes has been a challenging task for many decades. It generally requires the determination of the structural-dynamical information of protein solvent systems far from equilibrium over multiple length and time scales, which is still difficult both theoretically and experimentally. To cope with the problem, we introduce a full-residue space multiscale simulation method based on a combination of the kinetic Monte Carlo and molecular dynamics techniques, in which the rates of the rate-determining processes are evaluated from a biomolecular forcefield on the fly during the simulation run by taking into account the full space of residues. To demonstrate its reliability and efficiency, we explore the light-induced functional behavior of the full-length phototropin1 from Chlamydomonas reinhardtii (Cr-phot1) and its various subdomains. Our results demonstrate that in the dark state the light oxygen voltage-2-Jα (LOV2-Jα) photoswitch inhibits the enzymatic activity of the kinase, whereas the LOV1-Jα photoswitch controls the dimerization with the LOV2 domain. This leads to the repulsion of the LOV1-LOV2 linker out of the interface region between both LOV domains, which results in a positively charged surface suitable for cell-membrane interaction. By contrast, in the light state, we observe that the distance between both LOV domains is increased and the LOV1-LOV2 linker forms a helix-turn-helix (HTH) motif, which enables gene control through nucleotide binding. Finally, we find that the kinase is activated through the disruption of the Jα-helix from the LOV2 domain, which is followed by a stretching of the activation loop (A-loop) and broadening of the catalytic cleft of the kinase.

Towards a true protein movie: a perspective on the potential impact of the ensemble-based structure determination using exact NOEs

Confined by the Boltzmann distribution of the energies of the states, a multitude of structural states are inherent to biomolecules. For a detailed understanding of a protein's function, its entire structural landscape at atomic resolution and insight into the interconversion between all the structural states (i.e. dynamics) are required. Whereas dedicated trickery with NMR relaxation provides aspects of local dynamics, and 3D structure determination by NMR is well established, only recently have several attempts been made to formulate a more comprehensive description of the dynamics and the structural landscape of a protein. Here, a perspective is given on the use of exact NOEs (eNOEs) for the elucidation of structural ensembles of a protein describing the covered conformational space.

Alanine scan of core positions in ubiquitin reveals links between dynamics, stability, and function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.

Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: function follows fluctuations

Proteins exist as dynamic ensembles of molecules, implying that protein amino acid sequences evolved to code for both the ground state structure as well as the entire energy landscape of excited states. Accumulating theoretical and experimental evidence suggests that enzymes use such conformational fluctuations to facilitate allosteric processes important for substrate binding and possibly catalysis. This phenomenon can be clearly demonstrated in Escherichia coli adenylate kinase, where experimentally observed local unfolding of the LID subdomain, as opposed to a more commonly postulated rigid-body opening motion, is related to substrate binding. Because "entropy promoting" glycine mutations designed to increase specifically the local unfolding of the LID domain also affect substrate binding, changes in the excited energy landscape effectively tune the function of this enzyme without changing the ground state structure or the catalytic site. Thus, additional thermodynamic information, above and beyond the single folded structure of an enzyme-substrate complex, is likely required for a full and quantitative understanding of how enzymes work.

Heteronuclear Adiabatic Relaxation Dispersion (HARD) for quantitative analysis of conformational dynamics in proteins

NMR relaxation methods probe biomolecular motions over a wide range of timescales. In particular, the rotating frame spin-lock R(1ρ) and Carr-Purcell-Meiboom-Gill (CPMG) R(2) experiments are commonly used to characterize μs to ms dynamics, which play a critical role in enzyme folding and catalysis. In an effort to complement these approaches, we introduced the Heteronuclear Adiabatic Relaxation Dispersion (HARD) method, where dispersion in rotating frame relaxation rate constants (longitudinal R(1ρ) and transverse R(2ρ)) is created by modulating the shape and duration of adiabatic full passage (AFP) pulses. Previously, we showed the ability of the HARD method to detect chemical exchange dynamics in the fast exchange regime (k(ex)∼10(4)-10(5) s(-1)). In this article, we show the sensitivity of the HARD method to slower exchange processes by measuring R(1ρ) and R(2ρ) relaxation rates for two soluble proteins (ubiquitin and 10C RNA ligase). One advantage of the HARD method is its nominal dependence on the applied radio frequency field, which can be leveraged to modulate the dispersion in the relaxation rate constants. In addition, we also include product operator simulations to define the dynamic range of adiabatic R(1ρ) and R(2ρ) that is valid under all exchange regimes. We conclude from both experimental observations and simulations that this method is complementary to CPMG-based and rotating frame spin-lock R(1ρ) experiments to probe conformational exchange dynamics for biomolecules. Finally, this approach is germane to several NMR-active nuclei, where relaxation rates are frequency-offset independent.

The role of long-range intermolecular interactions in discovery of new drugs

INTRODUCTION

Long-range intermolecular interactions (interactions at distances between 100 and 1000 Å) play an important role in the interaction between drugs and therapeutic targets, and design techniques based on this concept could significantly improve and accelerate new drug discovery. Understanding these long-range intermolecular interactions will also help further our understanding of the molecular mechanisms and the underlying basic biological processes.

AREAS COVERED

This article looks at the physical bases of long-range intermolecular interactions in biological systems with a brief review of the literature data to support this concept. The article also gives some examples of techniques used in drug discovery that were based on the long-range intermolecular interaction concept.

EXPERT OPINION

The electron-ion interaction potential (EIIP) and average quasivalence number (AQVN) concepts shed new light on the role of long-range intermolecular interactions in biological systems. Further research of physicochemical mechanisms underlying long-range interactions between biological molecules is necessary for a better understanding of the basic biological processes. The addition of the computer-aided design techniques based on the EIIP/AQVN concept to the research and development will lead not only to a significant reduction in cost but also to an acceleration in the development of new drugs.

The feasibility of parameterizing four-state equilibria using relaxation dispersion measurements

Coupled equilibria play important roles in controlling information flow in biochemical systems, including allosteric molecules and multidomain proteins. In the simplest case, two equilibria are coupled to produce four interconverting states. In this study, we assessed the feasibility of determining the degree of coupling between two equilibria in a four-state system via relaxation dispersion measurements. A major bottleneck in this effort is the lack of efficient approaches to data analysis. To this end, we designed a strategy to efficiently evaluate the smoothness of the target function surface (TFS). Using this approach, we found that the TFS is very rough when fitting benchmark CPMG data to all adjustable variables of the four-state equilibria. After constraining a portion of the adjustable variables, which can often be achieved through independent biochemical manipulation of the system, the smoothness of TFS improves dramatically, although it is still insufficient to pinpoint the solution. The four-state equilibria can be finally solved with further incorporation of independent chemical shift information that is readily available. We also used Monte Carlo simulations to evaluate how well each adjustable parameter can be determined in a large kinetic and thermodynamic parameter space and how much improvement can be achieved in defining the parameters through additional measurements. The results show that in favorable conditions the combination of relaxation dispersion and biochemical manipulation allow the four-state equilibrium to be parameterized, and thus coupling strength between two processes to be determined.

Exploring sparsely populated states of macromolecules by diamagnetic and paramagnetic NMR relaxation

Sparsely populated states of macromolecules, characterized by short lifetimes and high free-energies relative to the predominant ground state, often play a key role in many biological, chemical, and biophysical processes. In this review, we briefly summarize various new developments in NMR spectroscopy that permit these heretofore invisible, sparsely populated states to be detected, characterized, and in some instances visualized. Relaxation dispersion spectroscopy yields detailed kinetic information on processes involving species characterized by distinct chemical shifts with lifetimes in the ∼50 μs-10 ms range and populations as low as 0.5%. In the fast exchange regime (time scale less than ∼250-500 μs), the footprint of sparsely populated states can be observed on paramagnetic relaxation enhancement profiles measured on the resonances of the major species, thereby yielding structural information that is directly related to paramagnetic center-nuclei distances from which it is possible, under suitable circumstances, to compute a structure or ensemble of structures for the minor species. Finally, differential transverse relaxation measurements can be used to detect lifetime broadening effects that directly reflect the unidirectional rates for the conversion of NMR-visible into high-molecular weight NMR-invisible species. Examples of these various approaches are presented.

Domain 3 of NS5A protein from the hepatitis C virus has intrinsic alpha-helical propensity and is a substrate of cyclophilin A

Nonstructural protein 5A (NS5A) is essential for hepatitis C virus (HCV) replication and constitutes an attractive target for antiviral drug development. Although structural data for its in-plane membrane anchor and domain D1 are available, the structure of domains 2 (D2) and 3 (D3) remain poorly defined. We report here a comparative molecular characterization of the NS5A-D3 domains of the HCV JFH-1 (genotype 2a) and Con1 (genotype 1b) strains. Combining gel filtration, CD, and NMR spectroscopy analyses, we show that NS5A-D3 is natively unfolded. However, NS5A-D3 domains from both JFH-1 and Con1 strains exhibit a propensity to partially fold into an α-helix. NMR analysis identifies two putative α-helices, for which a molecular model could be obtained. The amphipathic nature of the first helix and its conservation in all genotypes suggest that it might correspond to a molecular recognition element and, as such, promote the interaction with relevant biological partner(s). Because mutations conferring resistance to cyclophilin inhibitors have been mapped into NS5A-D3, we also investigated the functional interaction between NS5A-D3 and cyclophilin A (CypA). CypA indeed interacts with NS5A-D3, and this interaction is completely abolished by cyclosporin A. NMR heteronuclear exchange experiments demonstrate that CypA has in vitro peptidyl-prolyl cis/trans-isomerase activity toward some, but not all, of the peptidyl-prolyl bonds in NS5A-D3. These studies lead to novel insights into the structural features of NS5A-D3 and its relationships with CypA.

Molecular evolution of protein conformational changes revealed by a network of evolutionarily coupled residues

An improved understanding of protein conformational changes has broad implications for elucidating the mechanisms of various biological processes and for the design of protein engineering experiments. Understanding rearrangements of residue interactions is a key component in the challenge of describing structural transitions. Evolutionary properties of protein sequences and structures are extensively studied; however, evolution of protein motions, especially with respect to interaction rearrangements, has yet to be explored. Here, we investigated the relationship between sequence evolution and protein conformational changes and discovered that structural transitions are encoded in amino acid sequences as coevolving residue pairs. Furthermore, we found that highly coevolving residues are clustered in the flexible regions of proteins and facilitate structural transitions by forming and disrupting their interactions cooperatively. Our results provide insight into the evolution of protein conformational changes and help to identify residues important for structural transitions.

Millisecond dynamics in glutaredoxin during catalytic turnover is dependent on substrate binding and absent in the resting states

Conformational dynamics is important for enzyme function. Which motions of enzymes determine catalytic efficiency and whether the same motions are important for all enzymes, however, are not well understood. Here we address conformational dynamics in glutaredoxin during catalytic turnover with a combination of NMR magnetization transfer, R(2) relaxation dispersion, and ligand titration experiments. Glutaredoxins catalyze a glutathione exchange reaction, forming a stable glutathinoylated enzyme intermediate. The equilibrium between the reduced state and the glutathionylated state was biochemically tuned to exchange on the millisecond time scale. The conformational changes of the protein backbone during catalysis were followed by (15)N nuclear spin relaxation dispersion experiments. A conformational transition that is well described by a two-state process with an exchange rate corresponding to the glutathione exchange rate was observed for 23 residues. Binding of reduced glutathione resulted in competitive inhibition of the reduced enzyme having kinetics similar to that of the reaction. This observation couples the motions observed during catalysis directly to substrate binding. Backbone motions on the time scale of catalytic turnover were not observed for the enzyme in the resting states, implying that alternative conformers do not accumulate to significant concentrations. These results infer that the turnover rate in glutaredoxin is governed by formation of a productive enzyme-substrate encounter complex, and that catalysis proceeds by an induced fit mechanism rather than by conformer selection driven by intrinsic conformational dynamics.

Two-dimensional near-ultraviolet spectroscopy of aromatic residues in amyloid fibrils: a first principles study

We report a first principles study of two dimensional electronic spectroscopy of aromatic side chain transitions in the 32-residue β-amyloid (Aβ(9-40)) fibrils in the near ultraviolet (250-300 nm). An efficient exciton Hamiltonian with electrostatic fluctuations (EHEF) algorithm is used to compute the electronic excitations in the presence of environmental fluctuations. The through-space inter- and intra-molecular interactions are calculated with high level quantum mechanics (QM) approaches, and interfaced with molecular mechanics (MM) simulations. Distinct two dimensional near ultraviolet (2DNUV) spectroscopic signatures are identified for different aromatic transitions, and the couplings between them. 2DNUV signals associated with the transition couplings are shown to be very sensitive to the change of residue-residue interactions induced by residue mutations. Our simulations suggest that 2DNUV spectra could provide a useful local probe for the structure and kinetics of fibrils.

Intrinsically disordered PEP-19 confers unique dynamic properties to apo and calcium calmodulin

PEP-19 (Purkinje cell protein 4) is an intrinsically disordered protein with an IQ calmodulin (CaM) binding motif. Expression of PEP-19 was recently shown to protect cells from apoptosis and cell death due to Ca(2+) overload. Our initial studies showed that PEP-19 causes novel and dramatic increases in the rates of association of Ca(2+) with and dissociation of Ca(2+) from the C-domain of CaM. The goal of this work was to study interactions between the C-domain of CaM (C-CaM) and PEP-19 by solution nuclear magnetic resonance (NMR) to identify mechanisms by which PEP-19 regulates binding of Ca(2+) to CaM. Our results show that PEP-19 causes a greater structural change in apo C-CaM than in Ca(2+)-C-CaM, and that the first Ca(2+) binds preferentially to site IV in the presence of PEP-19 with exchange characteristics that are consistent with a decrease in Ca(2+) binding cooperativity. Relatively weak binding of PEP-19 has distinct effects on chemical and conformational exchange on the microsecond to millisecond time scale. In apo C-CaM, PEP-19 binding causes a redistribution of residues that experience conformational exchange, leading to an increase in the number of residues around Ca(2+) binding site IV that undergo conformational exchange on the microsecond to millisecond time scale. This appears to be caused by an allosteric effect because these residues are not localized to the PEP-19 binding site. In contrast, PEP-19 increases the number of residues that exhibit conformational exchange in Ca(2+)-C-CaM. These residues are primarily localized to the PEP-19 binding site but also include Asp93 in site III. These results provide working models for the role of protein dynamics in the regulation of binding of Ca(2+) to CaM by PEP-19.

An introduction to NMR-based approaches for measuring protein dynamics

Proteins are inherently flexible at ambient temperature. At equilibrium, they are characterized by a set of conformations that undergo continuous exchange within a hierarchy of spatial and temporal scales ranging from nanometers to micrometers and femtoseconds to hours. Dynamic properties of proteins are essential for describing the structural bases of their biological functions including catalysis, binding, regulation and cellular structure. Nuclear magnetic resonance (NMR) spectroscopy represents a powerful technique for measuring these essential features of proteins. Here we provide an introduction to NMR-based approaches for studying protein dynamics, highlighting eight distinct methods with recent examples, contextualized within a common experimental and analytical framework. The selected methods are (1) Real-time NMR, (2) Exchange spectroscopy, (3) Lineshape analysis, (4) CPMG relaxation dispersion, (5) Rotating frame relaxation dispersion, (6) Nuclear spin relaxation, (7) Residual dipolar coupling, (8) Paramagnetic relaxation enhancement. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein's functional dynamics. Normal mode analysis can provide such variables in cases where functional motions are dictated by the molecule's architecture. Principal component analysis is a more expensive way of inferring the essential variables from the protein's motions, one that requires a long molecular dynamics simulation. Finally, the two popular models for the allosteric switching mechanism, ‘preexisting equilibrium’ and ‘induced fit’, are interpreted within the energy landscape paradigm as extreme points of a continuum of transition mechanisms. Some experimental evidence illustrating each of these two models, as well as intermediate mechanisms, is presented and discussed.

Structural and biochemical characterization of the human cyclophilin family of peptidyl-prolyl isomerases

Peptidyl-prolyl isomerases catalyze the conversion between cis and trans isomers of proline. The cyclophilin family of peptidyl-prolyl isomerases is well known for being the target of the immunosuppressive drug cyclosporin, used to combat organ transplant rejection. There is great interest in both the substrate specificity of these enzymes and the design of isoform-selective ligands for them. However, the dearth of available data for individual family members inhibits attempts to design drug specificity; additionally, in order to define physiological functions for the cyclophilins, definitive isoform characterization is required. In the current study, enzymatic activity was assayed for 15 of the 17 human cyclophilin isomerase domains, and binding to the cyclosporin scaffold was tested. In order to rationalize the observed isoform diversity, the high-resolution crystallographic structures of seven cyclophilin domains were determined. These models, combined with seven previously solved cyclophilin isoforms, provide the basis for a family-wide structure:function analysis. Detailed structural analysis of the human cyclophilin isomerase explains why cyclophilin activity against short peptides is correlated with an ability to ligate cyclosporin and why certain isoforms are not competent for either activity. In addition, we find that regions of the isomerase domain outside the proline-binding surface impart isoform specificity for both in vivo substrates and drug design. We hypothesize that there is a well-defined molecular surface corresponding to the substrate-binding S2 position that is a site of diversity in the cyclophilin family. Computational simulations of substrate binding in this region support our observations. Our data indicate that unique isoform determinants exist that may be exploited for development of selective ligands and suggest that the currently available small-molecule and peptide-based ligands for this class of enzyme are insufficient for isoform specificity.

Ultraviolet spectroscopy of protein backbone transitions in aqueous solution: combined QM and MM simulations

A generalized approach combining quantum mechanics (QM) and molecular mechanics (MM) calculations is developed to simulate the n --> pi* and pi --> pi* backbone transitions of proteins in aqueous solution. These transitions, which occur in the ultraviolet (UV) at 180-220 nm, provide a sensitive probe for secondary structures. The excitation Hamiltonian is constructed using high-level electronic structure calculations of N-methylacetamide (NMA). Its electrostatic fluctuations are modeled using a new algorithm, EHEF, which combines a molecular dynamics (MD) trajectory obtained with a MM forcefield and electronic structures of sampled MD snapshots calculated by QM. The lineshapes and excitation splittings induced by the electrostatic environment in the experimental UV linear absorption (LA) and circular dichroism (CD) spectra of several proteins in aqueous solution are reproduced by our calculations. The distinct CD features of alpha-helix and beta-sheet protein structures are observed in the simulations and can be assigned to different backbone geometries. The fine structure of the UV spectra is accurately characterized and enables us to identify signatures of secondary structures.

Allosteric regulation of protease activity by small molecules

Proteases regulate a plethora of biological processes. Because they irreversibly cleave peptide bonds, the activity of proteases is strictly controlled. While there are many ways to regulate protease activity, an emergent mechanism is the modulation of protease function by small molecules acting at allosteric sites. This mode of regulation holds the potential to allow for the specific and temporal control of a given biological process using small molecules. These compounds also serve as useful tools for studying protein dynamics and function. This review highlights recent advances in identifying and characterizing natural and synthetic small molecule allosteric regulators of proteases and discusses their utility in studies of protease function, drug discovery and protein engineering.

At the dawn of the 21st century: Is dynamics the missing link for understanding enzyme catalysis?

Enzymes play a key role in almost all biological processes, accelerating a variety of metabolic reactions as well as controlling energy transduction, the transcription, and translation of genetic information, and signaling. They possess the remarkable capacity to accelerate reactions by many orders of magnitude compared to their uncatalyzed counterparts, making feasible crucial processes that would otherwise not occur on biologically relevant timescales. Thus, there is broad interest in understanding the catalytic power of enzymes on a molecular level. Several proposals have been put forward to try to explain this phenomenon, and one that has rapidly gained momentum in recent years is the idea that enzyme dynamics somehow contributes to catalysis. This review examines the dynamical proposal in a critical way, considering basically all reasonable definitions, including (but not limited to) such proposed effects as "coupling between conformational and chemical motions," "landscape searches" and "entropy funnels." It is shown that none of these proposed effects have been experimentally demonstrated to contribute to catalysis, nor are they supported by consistent theoretical studies. On the other hand, it is clarified that careful simulation studies have excluded most (if not all) dynamical proposals. This review places significant emphasis on clarifying the role of logical definitions of different catalytic proposals, and on the need for a clear formulation in terms of the assumed potential surface and reaction coordinate. Finally, it is pointed out that electrostatic preorganization actually accounts for the observed catalytic effects of enzymes, through the corresponding changes in the activation free energies.

Single-molecule force spectroscopy approach to enzyme catalysis

Enzyme catalysis has been traditionally studied using a diverse set of techniques such as bulk biochemistry, x-ray crystallography, and NMR. Recently, single-molecule force spectroscopy by atomic force microscopy has been used as a new tool to study the catalytic properties of an enzyme. In this approach, a mechanical force ranging up to hundreds of piconewtons is applied to the substrate of an enzymatic reaction, altering the conformational energy of the substrate-enzyme interactions during catalysis. From these measurements, the force dependence of an enzymatic reaction can be determined. The force dependence provides valuable new information about the dynamics of enzyme catalysis with sub-angstrom resolution, a feat unmatched by any other current technique. To date, single-molecule force spectroscopy has been applied to gain insight into the reduction of disulfide bonds by different enzymes of the thioredoxin family. This minireview aims to present a perspective on this new approach to study enzyme catalysis and to summarize the results that have already been obtained from it. Finally, the specific requirements that must be fulfilled to apply this new methodology to any other enzyme will be discussed.

Multiple roles of mobile active center loops in the E1 component of the Escherichia coli pyruvate dehydrogenase complex - Linkage of protein dynamics to catalysis

The region encompassing residues 401-413 on the E1 component of the pyruvate dehydrogenase multienzyme complex from Escherichia coli comprises a loop (the inner loop) which was not seen in the X-ray structure in the presence of thiamin diphosphate, the required cofactor for the enzyme. This loop is seen in the presence of a stable analogue of the pre-decarboxylation intermediate, the covalent adduct between the substrate analogue methyl acetylphosphonate and thiamin diphosphate, C2α-phosphonolactylthiamin diphosphate. It has been shown that the residue H407 and several other residues on this loop are required to reduce the mobility of the loop so electron density corresponding to it can be seen once the pre-decarboxylation intermediate is formed. Concomitantly, the loop encompassing residues 541-557 (the outer loop) appears to work in tandem with the inner loop and there is a hydrogen bond between the two loops ensuring their correlated motion. The inner loop was shown to: a) sequester the active center from carboligase side reactions; b) assist the interaction between the E1 and the E2 components, thereby affecting the overall reaction rate of the entire multienzyme complex; c) control substrate access to the active center. Using viscosity effects on kinetics it was shown that formation of the pre-decarboxylation intermediate is specifically affected by loop movement. A cysteine-less variant was created for the E1 component, onto which cysteines were substituted at selected loop positions. Introducing an electron spin resonance spin label and an (19)F NMR label onto these engineered cysteines, the loop mobility was examined: a) both methods suggested that in the absence of ligand, the loop exists in two conformations; b) line-shape analysis of the NMR signal at different temperatures, enabled estimation of the rate constant for loop movement, and this rate constant was found to be of the same order of magnitude as the turnover number for the enzyme under the same conditions. Furthermore, this analysis gave important insights into rate-limiting thermal loop dynamics. Overall, the results suggest that the dynamic properties correlate with catalytic events on the E1 component of the pyruvate dehydrogenase complex.