Catalytic Dynamics Research Articles

Perturbative diffraction methods resolve a conformational switch that facilitates a two-step enzymatic mechanism

Enzymes catalyze biochemical reactions through precise positioning of substrates, cofactors, and amino acids to modulate the transition-state free energy. However, the role of conformational dynamics remains poorly understood due to poor experimental access. This shortcoming is evident with Escherichia coli dihydrofolate reductase (DHFR), a model system for the role of protein dynamics in catalysis, for which it is unknown how the enzyme regulates the different active site environments required to facilitate proton and hydride transfer. Here, we describe ligand-, temperature-, and electric-field-based perturbations during X-ray diffraction experiments to map the conformational dynamics of the Michaelis complex of DHFR. We resolve coupled global and local motions and find that these motions are engaged by the protonated substrate to promote efficient catalysis. This result suggests a fundamental design principle for multistep enzymes in which pre-existing dynamics enable intermediates to drive rapid electrostatic reorganization to facilitate subsequent chemical steps.

HDX-MS Reveals Substrate-Dependent, Localized EX1 Conformational Dynamics in the Retaining GT-B Glycosyltransferase, MshA.

Glycosyltransferases (GTs) are well-characterized with respect to static 3D structures and molecular dynamics simulations, but there is a lack of reports on in-solution dynamics on time scales relevant to turnover. Here, backbone amide hydrogen/deuterium exchange followed by mass spectrometry (HDX-MS) was used to investigate the in-solution dynamics of the model retaining GT MshA from Corynebacterium glutamicum (CgMshA). CgMshA has a GT-B fold and catalyzes the transfer of N-acetyl-glucosamine (GlcNAc) from UDP-GlcNAc to l-myo-inositol-1-phosphate in the first step in mycothiol biosynthesis. HDX-MS results identify several key regions of conformational changes in response to UDP-GlcNAc binding, including residues 159-198 in the N-terminal domain and residues 323-354 in the C-terminal domain. These regions also exhibited substrate-dependent EX1 exchange kinetics consistent with conformational tension on the milliseconds to seconds time scale. A potential source of this conformational change is the flexible β4/α5 loop in the C-terminal domain, which sits at the interface of the two domains and likely interacts with the GlcNAc ring of UDP-GlcNAc. In contrast to UDP-GlcNAc, the UDP-CgMshA product complex exhibited severe decreases in deuterium incorporation, suggesting a less dynamic conformation. The HDX-MS results are complemented by solvent viscosity effects of 1.8-2.3 on the CgMshA kcat value, which are consistent with product release as a rate-determining step and possibly a direct role for protein dynamics in catalysis. The identification of in-solution dynamics that are sensitive to substrate binding allows for the proposal of a more detailed mechanism in CgMshA including conformation tension between the donor sugar and the flexible C-terminal domain β4/α5 loop providing sufficient conformational sampling for substrate-assisted catalysis to occur.

Comparison of the Role of Protein Dynamics in Catalysis by Dihydrofolate Reductase from E. coli and H. sapiens.

Dihydrofolate reductase (DHFR) is a well-studied, clinically relevant enzyme known for being highly dynamic over the course of its catalytic cycle. However, the role dynamic motions play in the explicit hydride transfer from the nicotinamide cofactor to the dihydrofolate substrate remains unclear because reaction initiation and direct spectroscopic examination on the appropriate time scale for such femtosecond to picosecond motions is challenging. Here, we employ pre-steady-state kinetics to observe the hydride transfer as directly as possible in two different species of DHFR: Escherichia coli and Homo sapiens. While the hydride transfer has been well-characterized in DHFR from E. coli, improvements in time resolution now allow for sub-millisecond dead times for stopped-flow spectroscopy, which reveals that the maximum rate is indeed faster than previously recorded. The rate in the human enzyme, previously only estimated, is also able to be directly observed using cutting-edge stopped-flow instrumentation. In addition to the pH dependence of the hydride transfer rates for both enzymes, we examine the primary H/D kinetic isotope effect to reveal a temperature dependence in the human enzyme that is absent from the E. coli counterpart. This dependence, which appears above a temperature of 15 °C is a shared feature among other hydride transfer enzymes and is also consistent with computational work suggesting the presence of a fast promoting-vibration that provides donor-acceptor compression on the time scale of catalysis to facilitate the chemistry step.

Kinetic Characterization and Computational Modeling of Escherichia coli Heptosyltransferase II: Exploring the Role of Protein Dynamics in Catalysis for GT-B Glycosyltransferase.

Glycosyltransferase (GT) enzymes promote the formation of glycosidic bonds between a sugar molecule and a diversity of substrates. Heptosyltransferase II (HepII) is a GT involved in the lipopolysaccharide (LPS) biosynthetic pathway that transfers the seven-carbon sugar (l-glycero-d-manno-heptose, Hep) onto a lipid-anchored glycopolymer (heptosylated Kdo2-Lipid A, Hep-Kdo2-Lipid A, or HLA). LPS plays a key role in Gram-negative bacterial sepsis, biofilm formation, and host colonization, and as such, LPS biosynthetic enzymes are targets for novel antimicrobial therapeutics. Three heptosyltransferases are involved in the inner-core LPS biosynthesis, with Escherichia coli HepII being the last to be quantitatively characterized in vivo. HepII shares modest sequence similarity with heptosyltransferase I (HepI) while maintaining a high degree of structural homology. Here, we report the first kinetic and biophysical characterization of HepII and demonstrate the properties of HepII that are shared with HepI, including sugar donor promiscuity and sugar acceptor-induced secondary structural changes, which results in significant thermal stabilization. HepII also has an increased catalytic efficiency and a significantly tighter binding affinity for both of its substrates compared to HepI. A structural model of the HepII ternary complex, refined by molecular dynamics simulations, was developed to probe the potentially important substrate-protein contacts. Ligand binding-induced changes in Trp fluorescence in HepII enabled the determination of substrate dissociation constants. Combined, these efforts meaningfully enhance our understanding of the heptosyltransferase family of enzymes and will aid in future efforts to design novel, potent, and specific inhibitors for this family of enzymes.

Examining enzyme active site dynamics using linear and nonlinear infrared spectroscopy: substrate orientation, cofactor composition, and solvent effects

The active sites of enzymes undergo motions spanning a broad range of timescales but the role of dynamics in catalysis is still unknown. We apply linear and nonlinear infrared spectroscopic techniques, combined with site-specific vibrational labeling, to gain insights into the ultrafast dynamics of the active site as a function of orientation of the substrate, composition of the active site, and the presence of cosolvents in the solution. We are examining these phenomena in three model enzymes bound to vibrationally labeled substrate analog probes.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects.

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

Crystal structure of the Thr316Ala mutant of a yeast JAMM deubiquitinase: implication of active-site loop dynamics in catalysis.

AMSH, an endosome-associated deubiquitinase (DUB) with a high specificity for Lys63-linked polyubiquitin chains, plays an important role in endosomal-lysosomal sorting and down-regulation of cell-surface receptors. AMSH belongs to the JAMM family of DUBs that contain two insertion segments, Ins-1 and Ins-2, in the catalytic domain relative to the JAMM core found in the archaebacterial AfJAMM. Structural analyses of the AMSH homologs human AMSH-LP and fission yeast Sst2 reveal a flap-like structure formed by Ins-2 near the active site that appears to open and close during its catalytic cycle. A conserved phenylalanine residue of the flap interacts with a conserved aspartate residue of the Ins-1 β-turn to form a closed `lid' over the active site in the substrate-bound state. Analyses of these two residues (Phe403 and Asp315) in Sst2 showed that their interaction plays an important role in controlling the flexibility of Ins-2. The Lys63-linked diubiquitin substrate-bound form of Sst2 showed that the conserved phenylalanine also interacts with Thr316 of Ins-1, which is substituted by tyrosine in other AMSH orthologs. Although Thr316 makes no direct interaction with the substrate, its mutation to alanine resulted in a significant loss of activity. In order to understand the contribution of Thr316 to catalysis, the crystal structure of this mutant was determined. In spite of the effect of the mutation on catalytic activity, the structure of the Sst2 Thr316Ala mutant did not reveal significant changes in either the overall structure or the active-site arrangement relative to the wild type. The Phe403-Thr316 van der Waals interaction is impaired by the Thr316Ala mutation, abrogating the adoption of the closed active-site conformation required for catalysis. Since van der Waals interactions with phenylalanine are conserved across substrate-bound forms of AMSH-LP and Sst2, these interactions may be critical for loop immobilization and the positioning of the isopeptide bond of Lys63-linked polyubiquitin-chain substrates.

Hot-carrier dynamics in catalysis

Abstract

Vicent Moliner

Vicent Moliner

Modulating Enzyme Activity by Altering Protein Dynamics with Solvent.

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized; however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ∼3 fold at 20% (v/v) IPA concentration. Stopped-flow kinetic measurements find that the pH-independent khydride rate decreases by 2.2 fold. X-ray crystallographic enzyme structures show no noticeable differences, while computational studies indicate that the transition state and electrostatic effects were identical for water and mixed solvent conditions; quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzyme's ability to access its functionally relevant conformational substates, explaining the decreased khydride rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

Protein Mass Effects on Formate Dehydrogenase.

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase.

Aromatic amine dehydrogenase (AADH) and related enzymes are at the heart of debates on the roles of quantum tunneling and protein dynamics in catalysis. The reaction of tryptamine in AADH involves significant quantum tunneling in the rate-limiting proton transfer step, shown by large H/D primary kinetic isotope effects (KIEs), with unusual temperature dependence. We apply correlated ab initio combined quantum mechanics/molecular mechanics (QM/MM) methods, at levels up to local coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), to calculate accurate potential energy surfaces for this reaction, which are necessary for quantitative analysis of tunneling contributions and reaction dynamics. Different levels of QM/MM treatment are tested. Multiple pathways are calculated with fully flexible transition state optimization by the climbing-image nudged elastic band method at the density functional QM/MM level. The average LCCSD(T) potential energy barriers to proton transfer are 16.7 and 14.0 kcal/mol for proton transfer to the two carboxylate atoms of the catalytic base, Asp128β. The results show that two similar, but distinct pathways are energetically accessible. These two pathways have different barriers, exothermicity and curvature, and should be considered in analyses of the temperature dependence of reaction and KIEs in AADH and other enzymes. These results provide a benchmark for this prototypical enzyme reaction and will be useful for developing empirical models, and analyzing experimental data, to distinguish between different conceptual models of enzyme catalysis.

Structural dynamics of Casein Kinase I (CKI) from malarial parasite Plasmodium falciparum (Isolate 3D7): Insights from theoretical modelling and molecular simulations

The protein kinases (PKs), belonging to serine/threonine kinase (STKs), are important drug targets for a wide spectrum of diseases in human. Among protein kinases, the Casein Kinases (CKs) are vastly expanded in various organisms, where, the malarial parasite Plasmodium falciparum possesses a single member i.e., PfCKI, which can phosphorylate various proteins in parasite extracts in vitro condition. But, the structure-function relationship of PfCKI and dynamics of ATP binding is yet to be understood. Henceforth, an attempt was made to study the dynamics, stability, and ATP binding mechanisms of PfCKI through computational modelling, docking, molecular dynamics (MD) simulations, and MM/PBSA binding free energy estimation. Bi-lobed catalytic domain of PfCKI shares a high degree of secondary structure topology with CKI domains of rice, human, and mouse indicating co-evolution of these kinases. Molecular docking study revealed that ATP binds to the active site where the glycine-rich ATP-binding motif (G16-X-G18-X-X-G21) along with few conserved residues plays a crucial role maintaining stability of the complex. Structural superposition of PfCKI with close structural homologs depicted that the location and length of important loops are different, indicating the dynamic properties of these loops among CKIs, which is consistent with principal component analysis (PCA). PCA displayed that the overall global motion of ATP-bound form is comparatively higher than that of apo form. The present study provides insights into the structural features of PfCKI, which could contribute towards further understanding of related protein structures, dynamics of catalysis and phosphorylation mechanism in these important STKs from malarial parasite in near future.

Structural and Kinetic Studies of Formate Dehydrogenase from Candida boidinii.

The structure of formate dehydrogenase from Candida boidinii (CbFDH) is of both academic and practical interests. First, this enzyme represents a unique model system for studies on the role of protein dynamics in catalysis, but so far these studies have been limited by the availability of structural information. Second, CbFDH and its mutants can be used in various industrial applications (e.g., CO2 fixation or nicotinamide recycling systems), and the lack of structural information has been a limiting factor in commercial development. Here, we report the crystallization and structural determination of both holo- and apo-CbFDH. The free-energy barrier for the catalyzed reaction was computed and indicates that this structure indeed represents a catalytically competent form of the enzyme. Complementing kinetic examinations demonstrate that the recombinant CbFDH has a well-organized reactive state. Finally, a fortuitous observation has been made: the apoenzyme crystal was obtained under cocrystallization conditions with a saturating concentration of both the cofactor (NAD(+)) and inhibitor (azide), which has a nanomolar dissociation constant. It was found that the fraction of the apoenzyme present in the solution is less than 1.7 × 10(-7) (i.e., the solution is 99.9999% holoenzyme). This is an extreme case where the crystal structure represents an insignificant fraction of the enzyme in solution, and a mechanism rationalizing this phenomenon is presented.

Ligand-Dependent Conformational Dynamics of Dihydrofolate Reductase.

Enzymes are known to change among several conformational states during turnover. The role of such dynamic structural changes in catalysis is not fully understood. The influence of dynamics in catalysis can be inferred, but not proven, by comparison of equilibrium structures of protein variants and protein-ligand complexes. A more direct way to establish connections between protein dynamics and the catalytic cycle is to probe the kinetics of specific protein motions in comparison to progress along the reaction coordinate. We have examined the enzyme model system dihydrofolate reductase (DHFR) from Escherichia coli with tryptophan fluorescence-probed temperature-jump spectroscopy. We aimed to observe the kinetics of the ligand binding and ligand-induced conformational changes of three DHFR complexes to establish the relationship among these catalytic steps. Surprisingly, in all three complexes, the observed kinetics do not match a simple sequential two-step process. Through analysis of the relationship between ligand concentration and observed rate, we conclude that the observed kinetics correspond to the ligand binding step of the reaction and a noncoupled enzyme conformational change. The kinetics of the conformational change vary with the ligand's identity and presence but do not appear to be directly related to progress along the reaction coordinate. These results emphasize the need for kinetic studies of DHFR with highly specific spectroscopic probes to determine which dynamic events are coupled to the catalytic cycle and which are not.

Dihydrofolate Reductase as a Model for Enzyme Catalysis

Dihydrofolate reductase (DHFR) has emerged as a model enzyme to investigate the role of dynamics in catalysis. This review begins with a description of the timing of the protonation of the dihydrofolate substrate and the subsequent hydride transfer reaction. This is followed by a summary of the results obtained from bioinformatics, theoretical and experimental studies of coevolving DHFR residues. Some of these residues appear to be dynamically coupled to each other and to the reaction coordinate. Importantly, the coupling appears to expand across the entire enzyme as many of the coevolved amino acids are remote from the active site. The importance of introducing mutations during evolution in the correct evolutionary order is also discussed. These studies involved a combination of steady state and pre-steady state rate measurements, kinetic isotope effects, NMR spectroscopy and computational investigations of the wild type DHFR from E. coli and several mutants. The cumulative results implicate a role of motions across the entire protein in DHFR catalysis, which may serve as a guide in studies of other enzymatic systems. Keywords: Dihydrofolate reductase, enzyme dynamics, enzyme evolution, global network, hydride transfer, kinetic isotope effect.

Role of dynamics in enzyme catalysis: substantial versus semantic controversies.

ConspectusThe role of the enzyme’s dynamic motions in catalysis is at the center of heated contemporary debates among both theoreticians and experimentalists. Resolving these apparent disputes is of both intellectual and practical importance: incorporation of enzyme dynamics could be critical for any calculation of enzymatic function and may have profound implications for structure-based drug design and the design of biomimetic catalysts.Analysis of the literature suggests that while part of the dispute may reflect substantial differences between theoretical approaches, much of the debate is semantic. For example, the term “protein dynamics” is often used by some researchers when addressing motions that are in thermal equilibrium with their environment, while other researchers only use this term for nonequilibrium events. The last cases are those in which thermal energy is “stored” in a specific protein mode and “used” for catalysis before it can dissipate to its environment (i.e., “nonstatistical dynamics”). This terminology issue aside, a debate has arisen among theoreticians around the roles of nonstatistical vs statistical dynamics in catalysis. However, the author knows of no experimental findings available today that examined this question in enzyme catalyzed reactions.Another source of perhaps nonsubstantial argument might stem from the varying time scales of enzymatic motions, which range from seconds to femtoseconds. Motions at different time scales play different roles in the many events along the catalytic cascade (reactant binding, reprotonation of reactants, structural rearrangement toward the transition state, product release, etc.). In several cases, when various experimental tools have been used to probe catalytic events at differing time scales, illusory contradictions seem to have emerged. In this Account, recent attempts to sort the merits of those questions are discussed along with possible future directions.A possible summary of current studies could be that enzyme, substrate, and solvent dynamics contribute to enzyme catalyzed reactions in several ways: first via mutual “induced-fit” shifting of their conformational ensemble upon binding; then via thermal search of the conformational space toward the reaction’s transition-state (TS) and the rare event of the barrier crossing toward products, which is likely to be on faster time scales then the first and following events; and finally via the dynamics associated with products release, which are rate-limiting for many enzymatic reactions. From a chemical perspective, close to the TS, enzymatic systems seem to stiffen, restricting motions orthogonal to the chemical coordinate and enabling dynamics along the reaction coordinate to occur selectively. Studies of how enzymes evolved to support those efficient dynamics at various time scales are still in their infancy, and further experiments and calculations are needed to reveal these phenomena in both enzymes and uncatalyzed reactions.

Stochastic gate dynamics regulate the catalytic activity of ubiquitination enzymes.

Initiation of the DNA damage and innate immune responses is dependent upon the flow of chemical information through coupled protein-protein interaction networks and driven by the synthesis and recognition of Lys 63 linked polyubiquitin (polyUb) chains on adaptor proteins. The central chemical step in Lys 63-linked protein ubiquitination involves the reaction of a specific lysine on a target protein with Ub that is covalently attached as a thioester conjugate to the Ub conjugating enzyme (E2) Ubc13. The active site cysteine of Ubc13, and E2 enzymes in general, is buttressed by a flexible loop. The role of loop dynamics in catalysis was investigated by mutating the central and hinge residues to glycine. The loop dynamics were experimentally characterized through measurement of enzyme kinetics, main chain NMR relaxation, X-ray crystallographic studies, and in vivo studies in yeast. The experimental data were complemented by analysis of MD simulations of the dynamics and kinetics for the loop motion. The results show that fast pico- to nanosecond time scale active site loop fluctuations play a crucial role in regulating the catalytic activity of Ubc13 by functioning as a stochastic active site gate, which is characterized by precisely balanced rates of opening and closing. In vivo functional complementation assays in yeast demonstrate that defects within this regulatory mechanism can have profound biological consequences, given that Ubc13 is the only E2 dedicated to synthesizing Lys 63-linked polyUb chains.

Mechanistic studies of different thymidylate synthases and the role of protein dynamics in catalysis (110.1)

Thymidylate synthases (TSases) catalyze the last committed step in the de novo biosynthesis of thymidylate, one of the four DNA bases (T). Classical TSases catalyze a complex cascade, and we developed methods to study the physical nature of its mechanism, and to test the role of enzyme dynamics in catalysis. In a related study, a newly discovered family of TSases was examined. We used a combination of isotopic labeling, mutations, quench-flow, and X-ray crystallography (Fig. 1), and found an unusual chemical mechanism for these flavin-dependent thymidylate synthases (FDTS). Figure 1. An active site view of crystal structures of TmFDTS in complex with FAD, dUMP, and the co-factor CH2H4folate.

Effects of Cavities at the Nicotinamide Binding Site of Liver Alcohol Dehydrogenase on Structure, Dynamics and Catalysis

A role for protein dynamics in enzymatic catalysis of hydrogen transfer has received substantial scientific support, but the connections between protein structure and catalysis remain to be established. Valine residues 203 and 207 are at the binding site for the nicotinamide ring of the coenzyme in liver alcohol dehydrogenase and have been suggested to facilitate catalysis with “protein-promoting vibrations” (PPV). We find that the V207A substitution has small effects on steady-state kinetic constants and the rate of hydrogen transfer; the introduced cavity is empty and is tolerated with minimal effects on structure (determined at 1.2 Å for the complex with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol). Thus, no evidence is found to support a role for Val-207 in the dynamics of catalysis. The protein structures and ligand geometries (including donor–acceptor distances) in the V203A enzyme complexed with NAD+ and 2,3,4,5,6-pentafluorobenzyl alcohol or 2,2,2-trifluoroethanol (determined at 1.1 Å) are very similar to those for the wild-type enzyme, except that the introduced cavity accommodates a new water molecule that contacts the nicotinamide ring. The structures of the V203A enzyme complexes suggest, in contrast to previous studies, that the diminished tunneling and decreased rate of hydride transfer (16-fold, relative to that of the wild-type enzyme) are not due to differences in ground-state ligand geometries. The V203A substitution may alter the PPV and the reorganization energy for hydrogen transfer, but the protein scaffold and equilibrium thermal motions within the Michaelis complex may be more significant for enzyme catalysis.