Prechemistry barriers and checkpoints do not contribute to fidelity and catalysis as long as they are not rate limiting

History of advances in enzyme kinetic methods: From minutes to milliseconds

The last review on transient-state kinetic methods in The Enzymes was published three decades ago (Johnson, K.A., 1992. The Enzymes, XX, 1-61). In that review the foundations were laid out for the logic behind the design and interpretation of experiments. In the intervening years the instrumentation has improved mainly in providing better sample economy and shorter dead times. More significantly, in 1992 we were just introducing methods for fitting data based on numerical integration of rate equations, but the software was slow and difficult to use. Today, advances in numerical integration methods for data fitting have led to fast and dynamic software, making it easy to fit data without simplifying approximations. This approach overcomes multiple disadvantages of traditional data fitting based on equations derived by analytical integration of rate equations, requiring simplifying approximations. Mechanism-based fitting using computer simulation resolves mechanisms by accounting for the concentration dependence of the rates and amplitudes of the reaction to find a set of intrinsic rate constants that reproduce the experimental data, including details about how the experiment was performed in modeling the data. Rather than discuss how to design and interpret rapid-quench and stopped-flow experiments individually, we now focus on how to fit them simultaneously so that the quench-flow data anchor the interpretation of fluorescence signals. Computer simulation streamlines the fitting of multiple experiments globally to yield a single unifying model to account for all available data.

High fidelity DNA strand-separation is the major specificity determinant in DNA methyltransferase CcrM's catalytic mechanism

Strand-separation is emerging as a novel DNA recognition mechanism but the underlying mechanisms and quantitative contribution of strand-separation to fidelity remain obscure. The bacterial DNA adenine methyltransferase, CcrM, recognizes 5'GANTC'3 sequences through a DNA strand-separation mechanism with unusually high selectivity. To explore this novel recognition mechanism, we incorporated Pyrrolo-dC into cognate and noncognate DNA to monitor the kinetics of strand-separation and used tryptophan fluorescence to follow protein conformational changes. Both signals are biphasic and global fitting showed that the faster phase of DNA strand-separation was coincident with the protein conformational transition. Non-cognate sequences did not display strand-separation and methylation was reduced > 300-fold, providing evidence that strand-separation is a major determinant of selectivity. Analysis of an R350A mutant showed that the enzyme conformational step can occur without strand-separation, so the two events are uncoupled. A stabilizing role for the methyl-donor (SAM) is proposed; the cofactor interacts with a critical loop which is inserted between the DNA strands, thereby stabilizing the strand-separated conformation. The results presented here are broadly applicable to the study of other N6-adenine methyltransferases that contain the structural features implicated in strand-separation, which are found widely dispersed across many bacterial phyla, including human and animal pathogens, and some Eukaryotes.

Design and interpretation of experiments to establish enzyme pathway and define the role of conformational changes in enzyme specificity

We describe the experimental methods and analysis to define the role of enzyme conformational changes in specificity based on published studies using DNA polymerases as an ideal model system. Rather than give details of how to perform transient-state and single-turnover kinetic experiments, we focus on the rationale of the experimental design and interpretation. We show how initial experiments to measure k_cat and k_cat/K_m can accurately quantify specificity but do not define its underlying mechanistic basis. We describe methods to fluorescently label enzymes to monitor conformational changes and to correlate fluorescence signals with rapid-chemical-quench flow assays to define the steps in the pathway. Measurements of the rate of product release and of the kinetics of the reverse reaction complete the kinetic and thermodynamic description of the full reaction pathway. This analysis showed that the substrate-induced change in enzyme structure from an open to a closed state was much faster than rate-limiting chemical bond formation. However, because the reverse of the conformational change was much slower than chemistry, specificity is governed solely by the product of the binding constant for the initial weak substrate binding and the rate constant for the conformational change (k_cat/K_m=K₁k₂) so that the specificity constant does not include k_cat. The enzyme conformational change leads to a closed complex in which the substrate is bound tightly and is committed to the forward reaction. In contrast, an incorrect substrate is bound weakly, and the rate of chemistry is slow, so the mismatch is released from the enzyme rapidly. Thus, the substrate-induced-fit is the major determinant of specificity. The methods outlined here should be applicable to other enzyme systems.

Single-molecule Taq DNA polymerase dynamics

Taq DNA polymerase functions at elevated temperatures with fast conformational dynamics-regimes previously inaccessible to mechanistic, single-molecule studies. Here, single-walled carbon nanotube transistors recorded the motions of Taq molecules processing matched or mismatched template-deoxynucleotide triphosphate pairs from 22° to 85°C. By using four enzyme orientations, the whole-enzyme closures of nucleotide incorporations were distinguished from more rapid, 20-μs closures of Taq's fingers domain testing complementarity and orientation. On average, one transient closure was observed for every nucleotide binding event; even complementary substrate pairs averaged five transient closures between each catalytic incorporation at 72°C. The rate and duration of the transient closures and the catalytic events had almost no temperature dependence, leaving all of Taq's temperature sensitivity to its rate-determining open state.

Conformational dynamics during misincorporation and mismatch extension defined using a DNA polymerase with a fluorescent artificial amino acid

High-fidelity DNA polymerases select the correct nucleotide over the structurally similar incorrect nucleotides with extremely high specificity while maintaining fast rates of incorporation. Previous analysis revealed the conformational dynamics and complete kinetic pathway governing correct nucleotide incorporation using a high-fidelity DNA polymerase variant containing a fluorescent unnatural amino acid. Here we extend this analysis to investigate the kinetics of nucleotide misincorporation and mismatch extension. We report the specificity constants for all possible misincorporations and characterize the conformational dynamics of the enzyme during misincorporation and mismatch extension. We present free energy profiles based on the kinetic measurements and discuss the effect of different steps on specificity. During mismatch incorporation and subsequent extension with the correct nucleotide, the rates of the conformational change and chemistry are both greatly reduced. The nucleotide dissociation rate, however, increases to exceed the rate of chemistry. To investigate the structural basis for discrimination against mismatched nucleotides, we performed all atom molecular dynamics simulations on complexes with either the correct or mismatched nucleotide bound at the polymerase active site. The simulations suggest that the closed form of the enzyme with a mismatch bound is greatly destabilized due to weaker interactions with active site residues, nonideal base pairing, and a large increase in the distance from the 3'-OH group of the primer strand to the α-phosphate of the incoming nucleotide, explaining the reduced rates of misincorporation. The observed kinetic and structural mechanisms governing nucleotide misincorporation reveal the general principles likely applicable to other high-fidelity DNA polymerases.

Optimized incorporation of an unnatural fluorescent amino acid affords measurement of conformational dynamics governing high-fidelity DNA replication

DNA polymerase from bacteriophage T7 undergoes large, substrate-induced conformational changes that are thought to account for high replication fidelity, but prior studies were adversely affected by mutations required to construct a Cys-lite variant needed for site-specific fluorescence labeling. Here we have optimized the direct incorporation of a fluorescent un-natural amino acid, (7-hydroxy-4-coumarin-yl)-ethylglycine, using orthogonal amber suppression machinery in Escherichia coli MS methods verify that the unnatural amino acid is only incorporated at one position with minimal background. We show that the single fluorophore provides a signal to detect nucleotide-induced conformational changes through equilibrium and stopped-flow kinetic measurements of correct nucleotide binding and incorporation. Pre-steady-state chemical quench methods show that the kinetics and fidelity of DNA replication catalyzed by the labeled enzyme are largely unaffected by the unnatural amino acid. These advances enable rigorous analysis to establish the kinetic and mechanistic basis for high-fidelity DNA replication.

Conformational dynamics during high-fidelity DNA replication and translocation defined using a DNA polymerase with a fluorescent artificial amino acid

We address the role of enzyme conformational dynamics in specificity for a high-fidelity DNA polymerase responsible for genome replication. We present the complete characterization of the conformational dynamics during the correct nucleotide incorporation forward and reverse reactions using stopped-flow and rapid-quench methods with a T7 DNA polymerase variant containing a fluorescent unnatural amino acid, (7-hydroxy-4-coumarin-yl) ethylglycine, which provides a signal for enzyme conformational changes. We show that the forward conformational change (>6000 s⁻¹) is much faster than chemistry (300 s⁻¹) while the enzyme opening to allow release of bound nucleotide (1.7 s⁻¹) is much slower than chemistry. These parameters show that the conformational change selects a correct nucleotide for incorporation through an induced-fit mechanism. We also measured conformational changes occurring after chemistry and during pyrophosphorolysis, providing new insights into processive polymerization. Pyrophosphorolysis occurs via a conformational selection mechanism as the pyrophosphate binds to a rare pretranslocation state of the enzyme–DNA complex. Global data fitting was achieved by including experiments in the forward and reverse directions to correlate conformational changes with chemical reaction steps. This analysis provided well-constrained values for nine rate constants to establish a complete free-energy profile including the rates of DNA translocation during processive synthesis. Translocation does not follow Brownian ratchet or power stroke models invoking nucleotide binding as the driving force. Rather, translocation is rapid and thermodynamically favorable after enzyme opening and pyrophosphate release, and it appears to limit the rate of processive synthesis at 4 °C.

Enzyme dynamics: Looking beyond a single structure

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

Exploring the Catalytic Mechanism of Cas9 Using Information Inferred from Endonuclease VII

Elucidating the nature of the gene editing mechanism of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) is an important task in view of the role of this breakthrough to the advancement of human medicine. In particular, it is crucial to understand the catalytic mechanism of Cas9 (one of the CRISPR associated proteins) and its role in confirming accurate editing. Thus, we focus in this work on an attempt to analyze the catalytic mechanism of Cas9. Considering the absence of detailed structural information on the active form of Cas9, we use an empirical valence bond (EVB) which is calibrated on the closely related mechanism of T4 endonuclease VII. The calibrated EVB is then used in studying the reaction of Cas9, while trying several structural models. It is found that the catalytic activation requires a large conformational change, where K848 or other positively charged group moves from a relatively large distance toward the scissile phosphate. This conformational change leads to the change in position of the Mg2+ ion and to a major reduction in the activation barrier for the catalytic reaction. Our finding provides an important clue on the nature of the catalytic activation of CAS9 and thus should help in elucidating a key aspect of the gene editing process. For example, the approach used here should be effective in exploring the nature of off target activation and its relationship to the energetics of the unwinding process. This strategy may offer ways to improve the selectivity of Cas9.

New standards for collecting and fitting steady state kinetic data

The Michaelis-Menten equation is usually expressed in terms of k _cat and K _m values: v = k _cat[S]/(K _m + [S]). However, it is impossible to interpret K _m in the absence of additional information, while the ratio of k _cat/K _m provides a measure of enzyme specificity and is proportional to enzyme efficiency and proficiency. Moreover, k _cat/K _m provides a lower limit on the second order rate constant for substrate binding. For these reasons it is better to redefine the Michaelis-Menten equation in terms of k _cat and k _cat/K _m values: v = k _SP[S]/(1 + k _SP[S]/k _cat), where the specificity constant, k _SP = k _cat/K _m. In this short review, the rationale for this assertion is explained and it is shown that more accurate measurements of k _cat/K _m can be derived directly using the modified form of the Michaelis-Menten equation rather than calculated from the ratio of k _cat and K _m values measured separately. Even greater accuracy is achieved with fitting the raw data directly by numerical integration of the rate equations rather than using analytically derived equations. The importance of fitting to derive k _cat and k _cat/K _m is illustrated by considering the role of conformational changes in enzyme specificity where k _cat and k _cat/K _m can reflect different steps in the pathway. This highlights the pitfalls in attempting to interpret K _m, which is best understood as the ratio of k _cat divided by k _cat/K _m.

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.

A new general method for simultaneous fitting of temperature and concentration dependence of reaction rates yields kinetic and thermodynamic parameters for HIV reverse transcriptase specificity

Recent studies have demonstrated the dominant role of induced fit in enzyme specificity of HIV reverse transcriptase and many other enzymes. However, relevant thermodynamic parameters are lacking, and equilibrium thermodynamic methods are of no avail because the key parameters can only be determined by kinetic measurement. By modifying KinTek Explorer software, we present a new general method for globally fitting data collected over a range of substrate concentrations and temperatures and apply it to HIV reverse transcriptase. Fluorescence stopped-flow methods were used to record the kinetics of enzyme conformational changes that monitor nucleotide binding and incorporation. The nucleotide concentration dependence was measured at temperatures ranging from 5 to 37 °C, and the raw data were fit globally to derive a single set of rate constants at 37 °C and a set of activation enthalpy terms to account for the kinetics at all other temperatures. This comprehensive analysis afforded thermodynamic parameters for nucleotide binding (K_d , ΔG, ΔH, and ΔS at 37 °C) and kinetic parameters for enzyme conformational changes and chemistry (rate constants and activation enthalpy). Comparisons between wild-type enzyme and a mutant resistant to nucleoside analogs used to treat HIV infections reveal that the ground state binding is weaker and the activation enthalpy for the conformational change step is significantly larger for the mutant. Further studies to explore the structural underpinnings of the observed thermodynamics and kinetics of the conformational change step may help to design better analogs to treat HIV infections and other diseases. Our new method is generally applicable to enzyme and chemical kinetics.

Induced fit of the peptidyl-transferase center of the ribosome and conformational freedom of the esterified amino acids

The catalytic site of most enzymes can efficiently handle only one substrate. In contrast, the ribosome is capable of polymerizing at a similar rate at least 20 different kinds of amino acids from aminoacyl-tRNA carriers while using just one catalytic site, the peptidyl-transferase center (PTC). An induced-fit mechanism has been uncovered in the PTC, but a possible connection between this mechanism and the uniform handling of the substrates has not been investigated. We present an analysis of published ribosome structures supporting the hypothesis that the induced fit eliminates unreactive rotamers predominantly populated for some A-site aminoacyl esters before induction. We show that this hypothesis is fully consistent with the wealth of kinetic data obtained with these substrates. Our analysis reveals that induction constrains the amino acids into a reactive conformation in a side-chain independent manner. It allows us to highlight the rationale of the PTC structural organization, which confers to the ribosome the very unusual ability to handle large as well as small substrates.

Exploring the mechanism of DNA polymerases by analyzing the effect of mutations of active site acidic groups in Polymerase β

Elucidating the catalytic mechanism of DNA polymerase is crucial for a progress in the understanding of the control of replication fidelity. This work tries to advance the mechanistic understanding by analyzing the observed effect of mutations of the acidic groups in the active site of Polymerase β as well as the pH effect on the rate constant. The analysis involves both empirical valence bond (EVB) free energy calculations and considerations of the observed pH dependence of the reaction. The combined analysis indicates that the proton transfer (PT) from the nucleophilic O3' has two possible pathways, one to D256 and the second to the bulk. We concluded based on calculations and the experimental pH profile that the most likely path for the wild-type (WT) and the D256E and D256A mutants is a PT to the bulk, although the WT may also use a PT to Asp 256. Our analysis highlights the need for very extensive sampling in the calculations of the activation barrier and also clearly shows that ab initio QM/MM calculations that do not involve extensive sampling are unlikely to give a clear quantitative picture of the reaction mechanism. Proteins 2016; 84:1644-1657. © 2016 Wiley Periodicals, Inc.

Enzyme Selectivity of HIV Reverse Transcriptase: Conformations, Ligands, and Free Energy Partition

Atomically detailed simulations of HIV RT are performed to investigate the contributions of the conformational transition to the overall rate and specificity of enzyme catalysis. A number of different scenarios are considered within Milestoning theory to provide a more complete picture of the process of opening and closing the enzyme. We consider the open to closed transition in the absence of and with the correct and incorrect substrates. We also consider the free energy profile and the kinetics of the conformational change after the chemistry step in which a new base was added to the DNA, but the DNA was not yet displaced. We partition the free energy along the reaction coordinate and analyze the importance of different protein domains. Strikingly, significant influence on the free energy profile is detected for amino acids far from the active site. The overall long-range impact is about 50 percent of the total. We also illustrate that the overall rate is not necessarily determined by the highest free energy barrier along the reaction path (with respect to the free enzyme and substrate) and that the specificity is not necessarily determined by the same reaction step that determines the rate.

The Closing Mechanism of DNA Polymerase I at Atomic Resolution

DNA polymerases must quickly and accurately distinguish between similar nucleic acids to form Watson-Crick base pairs and avoid DNA replication errors. Deoxynucleoside triphosphate (dNTP) binding to the DNA polymerase active site induces a large conformational change that is difficult to characterize experimentally on an atomic level. Here, we report an X-ray crystal structure of DNA polymerase I bound to DNA in the open conformation with a dNTP present in the active site. We use this structure to computationally simulate the open to closed transition of DNA polymerase in the presence of a Watson-Crick base pair. Our microsecond simulations allowed us to characterize the key steps involved in active site assembly, and propose the sequence of events involved in the prechemistry steps of DNA polymerase catalysis. They also reveal new features of the polymerase mechanism, such as a conserved histidine as a potential proton acceptor from the primer 3'-hydroxyl.

Noncognate DNA damage prevents the formation of the active conformation of the Y-family DNA polymerases DinB and DNA polymerase κ

Y-family DNA polymerases are specialized to copy damaged DNA, and are associated with increased mutagenesis, owing to their low fidelity. It is believed that the mechanism of nucleotide selection by Y-family DNA polymerases involves conformational changes preceding nucleotidyl transfer, but there is limited experimental evidence for such structural changes. In particular, nucleotide-induced conformational changes in bacterial or eukaryotic Y-family DNA polymerases have, to date, not been extensively characterized. Using hydrogen-deuterium exchange mass spectrometry, we demonstrate here that the Escherichia coli Y-family DNA polymerase DinB and its human ortholog DNA polymerase κ undergo a conserved nucleotide-induced conformational change in the presence of undamaged DNA and the correct incoming nucleotide. Notably, this holds true for damaged DNA containing N(2) -furfuryl-deoxyguanosine, which is efficiently copied by these two polymerases, but not for damaged DNA containing the major groove modification O(6) -methyl-deoxyguanosine, which is a poor substrate. Our observations suggest that DinB and DNA polymerase κ utilize a common mechanism for nucleotide selection involving a conserved prechemical conformational transition promoted by the correct nucleotide and only preferred DNA substrates.

DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation

By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.

Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization

We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.

Conformational dynamics of Thermus aquaticus DNA polymerase I during catalysis

Despite the fact that DNA polymerases have been investigated for many years and are commonly used as tools in a number of molecular biology assays, many details of the kinetic mechanism they use to catalyze DNA synthesis remain unclear. Structural and kinetic studies have characterized a rapid, pre-catalytic open-to-close conformational change of the Finger domain during nucleotide binding for many DNA polymerases including Thermus aquaticus DNA polymerase I (Taq Pol), a thermostable enzyme commonly used for DNA amplification in PCR. However, little has been performed to characterize the motions of other structural domains of Taq Pol or any other DNA polymerase during catalysis. Here, we used stopped-flow Förster resonance energy transfer to investigate the conformational dynamics of all five structural domains of the full-length Taq Pol relative to the DNA substrate during nucleotide binding and incorporation. Our study provides evidence for a rapid conformational change step induced by dNTP binding and a subsequent global conformational transition involving all domains of Taq Pol during catalysis. Additionally, our study shows that the rate of the global transition was greatly increased with the truncated form of Taq Pol lacking the N-terminal domain. Finally, we utilized a mutant of Taq Pol containing a de novo disulfide bond to demonstrate that limiting protein conformational flexibility greatly reduced the polymerization activity of Taq Pol.

Structure-based energetics of mRNA decoding on the ribosome

The origin of high fidelity in bacterial protein synthesis on the ribosome remains a fundamental unsolved problem despite available three-dimensional structures of different stages of the translation process. However, these structures open up the possibility of directly computing the energetics of tRNA selection that is required for an authentic understanding of fidelity in decoding. Here, we report extensive computer simulations that allow us to quantitatively calculate tRNA discrimination and uncover the energetics underlying accuracy in code translation. We show that the tRNA-mRNA interaction energetics varies drastically along the path from initial selection to peptide bond formation. While the selection process is obviously controlled by kinetics, the underlying thermodynamics explains the origin of the high degree of accuracy. The existence of both low- and high-selectivity states provides an efficient mechanism for initial selection and proofreading that does not require codon-dependent long-range structural signaling within the ribosome. It is instead the distinctly unequal population of the high-selectivity states for cognate and noncognate substrates that is the key discriminatory factor. The simulations reveal the essential roles played both by the 30S subunit conformational switch and by the common tRNA modification at position 37 in amplifying the accuracy.