Understanding the underlying physical mechanism of the impressive rate accelerations achieved by enzymes has been a central question in biology. Transition-state (TS) theory has provided the fundamental framework since the rate of a reaction is dictated by the free-energy difference between the ground state and TS (Truhlar, 2015). A wealth of research has led to the general notion that enzymatic rate acceleration is due to the enzymes’ much higher affinity to the TS relative to its substrates. This idea has been supported experimentally by the high affinities measured for transition-state analogs (TSA) relative to substrates, leading to the design of TSA as high-affine enzyme inhibitors (Lienhard, 1973; Schramm, 2007). Since the TS represents a maximum in the free-energy landscape, its structure cannot be determined directly due to its low probability and ephemeral nature, TSA-bound species have been extensively used as structural proxies.

The classic TS theory and structural visualization of enzyme/TSA complexes lead to a current view of quite unique structures on the dividing surface at the maximum in the free-energy landscape. In contrast, it is now generally accepted that proteins are large ensembles of conformations, a concept rooted in pioneering work by Frauenfelder and coworkers for protein function, in analogy to the ensemble/funnel concept of protein folding (Frauenfelder and Wolynes, 1985). Consequently, enzyme–substrate (ES) complexes – experimentally accessible as minima in the free-energy landscape – are composed of a vast ensemble of molecular configurations due to their macromolecular nature. Hence, catalysis by pathing through exclusive TS’s would pose a huge entropic bottleneck/energy barrier for enzyme-catalyzed reactions.

Motivated by this consideration, together with the marginal outcomes in current enzyme design relative to the catalytic power of naturally evolved enzymes, we set out to investigate the atomistic nature of TSs of the chemical step in an enzymatic cycle and its implication for rate acceleration by combining quantum-mechanics/molecular-mechanics (QM/MM) calculations and experiments. We take advantage of the groundbreaking work on TS theory by QM simulation developed in the field (Kamerlin et al., 2013; Masgrau and Truhlar, 2015; Warshel and Bora, 2016; Zinovjev and Tuñón, 2017), well summarized in several reviews (Cui, 2016; Senn and Thiel, 2009; van der Kamp and Mulholland, 2013). We chose a phosphoryl-transfer (P-transfer) reaction due to its impressive enzymatic rate acceleration, as the uncatalyzed reactions have extremely high-energy barriers (Kerns et al., 2015; Lassila et al., 2011). The chemical rationale for these high barriers and reasons as to why P-transfer reactions are ubiquitous in living organisms and are involved in almost all essential biological processes is elegantly discussed in a fundamental paper by Westheimer in 1987: ‘Why nature chose phosphates’ (Westheimer, 1987). The central role of enzyme-catalyzed P-transfer reactions in biology (i.e., genetic code, energy storage, signaling, and compartmentalization) has led to extensive discussions and controversies of how enzymes catalyze these vital reactions so efficiently (Allen and Dunaway-Mariano, 2016; Kamerlin et al., 2013; Kerns et al., 2015; Lassila et al., 2011; Pabis et al., 2016; Westheimer, 1987). Here, we combine QM/MM calculations with experimental temperature- and pH-dependent kinetic studies and X-ray crystallography of adenylate kinase (Adk) to shed light on the question of TS structures for the chemical step in its catalytic cycle. While QM/MM methods have been extensively applied to study enzymatic reactions, yielding accurate and detailed descriptions of the molecular events during catalysis (Acevedo and Jorgensen, 2010; Cheng et al., 2005; Dinner et al., 2001; Hahn et al., 2015; Hirvonen et al., 2020; Karplus and Kuriyan, 2005; Lai and Cui, 2020a; López-Canut et al., 2011; Mones et al., 2013; Palermo et al., 2015; Rosta et al., 2014; Roston et al., 2018; Roston and Cui, 2016; Schwartz and Schramm, 2009; Senn and Thiel, 2009; Turjanski et al., 2009), we test here the reaction free-energy profiles (FEPs) from simulations subsequently by designed experiments. We deliver a revised notion of TS stabilization by the enzyme through our discovery of a structurally broad transition-state ensemble (TSE).

Adk is an extensively studied phosphotransferase that catalyzes the reversible, roughly isoenergetic conversion of two ADP molecules into ATP and AMP (Figure 1A, B), thereby maintaining the cellular concentrations of these nucleotides (Dzeja and Terzic, 2009). It is an essential enzyme found in every cell and organism. Extensive structural, dynamic, and kinetic studies led to a comprehensive picture of the overall Adk reaction mechanism and its underlying energy landscape (Figure 1A; Beckstein et al., 2009; Berry et al., 1994; Berry et al., 2006; Henzler-Wildman et al., 2007; Kerns et al., 2015; Müller et al., 1996; Müller and Schulz, 1992; Wolf-Watz et al., 2004). The conformational change of the AMP- and ATP-lids closing and opening is crucial in the enzymatic cycle: Lid-closing positions the two substrates and the active-site residues for efficient chemistry (i.e., P-transfer) and prohibits the alternative, energetically favorable reaction of phosphoryl hydrolysis (Kerns et al., 2015). Lid-opening, essential for product release, and not P-transfer is the rate-limiting step in the catalytic cycle and Mg²⁺ accelerates both steps (Kerns et al., 2015). Here, we investigate the actual P-transfer step computationally, followed by experimental testing of our computational results. This chemical step would take about 7000 years without the enzyme (Stockbridge and Wolfenden, 2009) compared to >5000 s⁻¹ with the enzyme (Kerns et al., 2015).

Figure 1 with 6 supplements see all

Download asset Open asset

QM/MM calculations of the P-transfer step in Adk

In order to shed light on the mechanism by which Adk from Aquifex aeolicus catalyzes the phosphoryl-transfer reaction by more than 12 orders of magnitude (Kerns et al., 2015), we performed QM/MM simulations starting with two ADP molecules and Mg²⁺ bound in the active site. The starting structures for the simulations were prepared using the X-ray structure of Adk in complex with Ap5A (P1,P5-Di(adenosine-5′) pentaphosphate) and coordinated to Zn²⁺ (2RGX; Henzler-Wildman et al., 2007). ADP–ADP coordinates were built using Ap5A as a template and Zn²⁺ was replaced by Mg²⁺. For the QM/MM simulations, the QM region was defined as the diphosphate moiety of both ADP molecules, the Mg²⁺ ion, plus the four coordinating water molecules. The rest of the system was described at molecular-mechanics level using the AMBER ff99sb force field (Hornak et al., 2006) and was solvated with TIP3P water molecules (Figure 1—figure supplement 1).

The equilibrated starting structures agree well with the X-ray structures of the enzyme bound to Mg²⁺/ADP (4CF7; Kerns et al., 2015; Figure 1—figure supplement 2). Steered molecular dynamics simulations were run in both the forward (ADP/ADP to ATP/AMP) and reverse direction (ATP/AMP to ADP/ADP) with Mg²⁺ present in the active site (Figure 1—figure supplement 3). Since it is unknown whether the fully charged or monoprotonated nucleotide (on one β-ADP oxygen) state is the more reactive configuration, we performed QM/MM simulations for both cases. FEPs of the P-transfer were determined using Multiple Steered Molecular Dynamics and Jarzynski’s Relationship (Crespo et al., 2005; Jarzynski, 1997; Ramírez et al., 2014; Figure 1C, Figure 1—figure supplements 3–5). The results reveal a much smaller free energy of activation for the fully charged nucleotide state (Δ_fG^‡ of 13 ± 0.9 kcal/mol) relative to monoprotonated state (Δ_fG^‡ of 23 ± 0.9 kcal/mol) (Table 1); hence, we conclude that this is the most reactive enzyme configuration.

Table 1

	With Mg2+	Without Mg2+	With Mg2+, ADP monoprotonated
ΔG*	−6 (1.7)	+4 (2.5)	+6 (1.9)
ΔfG†	13 (0.9)	34 (1.6)	23 (0.9)
ΔbG†	20 (0.8)	30 (0.9)	18 (0.9)
ξ(TS)	−0.5 to 0.7 ‡ 0.0 to 0.2	−0.2 to 0 ‡ −0.1	−0.3 to −0.1 ‡ Not calculated
ζ(TS)	165 to 190	165 to 180	170 to 180

1. ξ(TS) is the range of the reaction coordinate in the TSE (in Å); ζ(TS) is the improper dihedral angle of the transferring phosphate in the TS. The estimated errors of the free energies are in parenthesis and are computed as described in Materials and methods.

2. *

   overall reaction free energy.

3. †

   activation free energy of the forward reaction; activation free energy of the backward reaction.

4. ‡

   Values were analyzed by a visual analysis of the multiple steered molecular dynamics (MSMD) trajectories, no committor distribution were calculated.

When the calculations were repeated in the absence of Mg²⁺ (with the nucleotide monoprotonated, since the fully charged nucleotide state prohibited the reaction), a large increase in the free-energy activation barrier was observed relative to the Mg²⁺-bound system (Figure 1C, Figure 1—figure supplement 5, Δ_fG^‡ of 34 ± 1.6 kcal/mol), in agreement with expectations from experiments (Kerns et al., 2015). We note that only a lower limit for the overall acceleration by Mg²⁺ (>10⁵-fold) could be estimated from published results (Kerns et al., 2015), since the P-transfer was too fast to be measured experimentally in the presence of Mg²⁺.

During revision, we tested the accuracy and robustness of our results from our original QM/MM level of theory by employing higher-level pure DFT(PBE) free-energy calculations similar to the approaches used by Ganguly et al., 2020; Figure 1—figure supplement 6. The results are in excellent agreement with the FEPs obtained before (Figure 1C, E; Figure 2A, C) thus further supporting the computational results.

Figure 2

Download asset Open asset

The TSE – transferring phosphoryl group delocalized

The nature of the TS of enzyme-catalyzed P-transfer reaction with respect to its associative and dissociative character has been of significant interest and heated debate (Kamerlin et al., 2013; Kamerlin and Wilkie, 2007; Lassila et al., 2011; Roston and Cui, 2016). The definition as well as theoretical and experimental approaches to distinguish between them are rooted in elegant and fundamental work on nonenzymatic P-transfer reactions (Duarte et al., 2015; Hengge, 2002; Hou et al., 2012; Hou and Cui, 2012; Kamerlin et al., 2013; Kerns et al., 2015; Kirby and Nome, 2015; Stockbridge and Wolfenden, 2009). For Adk, a concerted mechanism is observed with no intermediates (Figures 1 and 2). Instead of using the classification of associative versus dissociative character (Lassila et al., 2011), we rather apply the clear definition by Cui and coworkers of a tight versus loose TS for concerted P-transfer reactions as they are unambiguously related to the bond-order in the TS (Lai and Cui, 2020a; Lai and Cui, 2020b; Roston et al., 2018; Roston and Cui, 2016). For Adk, a contraction of the active site is observed while passing through the TSE, as captured by the decrease in the distance between acceptor and leaving group during the transition (Figure 1D, E). In the absence of Mg²⁺, the acceptor and donor need to be within 4 Å for the reaction to occur (Figure 1D) whereas in the presence of Mg²⁺ an acceptor–donor distance of 4.5 Å is sufficient (Figure 1E). The character of the TSE can best be seen from the widely used Moore-O’Ferrall–Jencks diagram (Jencks, 1972; O’Ferrall, 1970) that plots the distance from the transferring phosphate to the oxygen of the donor and acceptor, respectively (Figure 2). The TSE changes from tight/synchronous without Mg²⁺ (meaning high bond-order for both bonds between transferring phosphate to the leaving group and the attacking nucleophile) into loose with Mg²⁺ (lower bond-order for these two bonds; Figures 2 and 3; Lai and Cui, 2020a; Lai and Cui, 2020b; Roston et al., 2018; Roston and Cui, 2016).

Importantly, from these two-dimensional (Figure 2) and the one-dimensional (Figure 1C) free-energy plots, we noticed a novel striking feature: a large ensemble of conformations in the TS region with vast differences in the position of the transferring phosphate but equal values in free energy for the enzyme with Mg²⁺ (light blue dots in Figure 2, yellow area in Figure 1C). In other words, the enzyme seems to operate with a wide TSE. In contrast, without Mg²⁺ fewer conformations seem to comprise the TSE. To dive more into this unexpected computational result and better visualize the TSE, we show zoom-ins of representative TSE snapshots and their superposition (Figure 3, Figure 3—figure supplement 1). Notably, the TSE conformations are distributed along the reaction coordinate. The difference in position of the transferring phosphate along the reaction coordinate within the TSE (about 1 Å with Mg²⁺; Figure 1c, Figure 2 and Figure 3b) implies that the TSE contains many highly asymmetric conformations, meaning that the transferring phosphate can be much closer to the leaving oxygen than the attacking oxygen and vice versa (Figure 3D). This asymmetry is logically tied with nonplanar configurations of the transferring phosphate (Figure 3D and Table 1).

Figure 3 with 9 supplements see all

Download asset Open asset

Comparison of the TSE snapshots for the fully active, Mg²⁺-bound enzyme obtained from QM/MM with an X-ray structure of the same enzyme in complex with the TSA (ADP–Mg²⁺–AlF_4—AMP) (Kerns et al., 2015) serves as an initial experimental validation of our simulations (Figure 3C, F, Figure 3—figure supplement 2). At the same time, this comparison highlights the power of the QM/MM simulations to investigate the catalytic mechanism, as the TSA-bound X-ray structure seems to imply quite a unique TS conformation, in sharp contrast to the broad TSE discovered in our simulations.

A detailed analysis of the pathways with and without Mg²⁺ reveals well-known features from other P-transfer enzymes, such as coordinated changes of the Mulliken charge populations coupled with proton transfers (Figure 3—figure supplements 3–5). Mg²⁺ achieves its catalytic effect on the chemical step playing several roles. It enables the donor and acceptor phosphates to come close enough to react, thus achieving a reactive conformation. In the absence of Mg²⁺, a proton must bridge the terminal phosphates to prevent the strong electrostatic repulsion between the two nucleotides. In addition, Mg²⁺ stabilizes the charge of the transferring phosphate and thereby lowering the enthalpic barrier. Most interestingly, the presence of a cation seems to result in a wide TSE. Structures of the TSE show that Mg²⁺ is optimally coordinated to the transferring phosphates and water molecules (Figure 3G).

Besides the multifaceted role of Mg²⁺, our simulations provide insights into the function of the fully conserved arginine residues in the active site for lowering the activation barrier for this chemical reaction (Figure 3—figure supplements 6–8). R85 appears to arrange the beta-phosphates of the ADPs in the proper position for the reaction to be started. R150 and R161 are involved in anchoring ATP and AMP residues in the backward reaction. These last two residues and R85 interact with the Pα of AMP, stabilizing its negative charge. R36 is near to the alpha-phosphate of AD(M)P along the entire catalytic reaction, stabilizing both the TS and products. Overall, interactions between these arginine side chains and several backbone amides and the Mg²⁺ with the phosphates of the substrate make up a well-organized, asymmetric active site enabling efficient reversible P-transfer with a wide TSE (Figure 3G).

Umbrella sampling and committor analysis buttress wide TSE for Adk with Mg²⁺

To verify our computational finding of a wide TSE in the fully assembled enzyme with Mg²⁺ in contrast to a narrower TSE without a divalent metal, we performed QM/MM umbrella sampling and a committor analysis as described in methods. The umbrella sampling FEPs (Figure 4a–d) corroborate the key difference in the range of the TSE. Noteworthy, the activation free energies are smaller for the umbrella sampling when compared to the multiple steered molecular dynamics (MSMD). The largest difference being that for the reaction in the absence of Mg (Figure 1C). Higher barriers when comparing the results using Jarzynski’s to those obtained with umbrella sampling is not totally unexpected. It is well known that exponential averaging tends to overestimate free-energy barriers, and particularly when they are too high (Park et al., 2003).

Figure 4

Download asset Open asset

Notably, the calculated difference in activation barriers with and without Mg from these umbrella simulations of about 9.5 kcal/mol is in good agreement with the experimentally determined ones of ≥11 kcal/mol (P-transfer with Mg is >500 s⁻¹ (Kerns et al., 2015) compared to 7.5 × 10⁻⁴ s⁻¹ measured here).

A consequent comparative commitment analysis that calculates the probability of reaching either reactants or products displays a very shallow change along the reaction coordinate in the presence of Mg²⁺, further corroborating a wide TSE for the P-transfer step in the fully active enzyme, in stark contrast to a narrow and steep change in the committors indicative of a narrower TSE without Mg²⁺ (Figure 4d). The TSE is defined in function of the committor distribution (Figure 3—figure supplement 9), showing the TSE is at −0.1 Å of the reaction coordinate for the reaction without Mg²⁺. In contrast, the TSE is wider for the reaction with Mg²⁺ (reaction coordinate values in the range from 0.0 to 0.2 Å; Figure 3—figure supplement 9 and Table 1).

Experimental characterization of the activation parameters of the P-transfer step

We felt the need to additionally experimentally test our major new finding from the QM/MM simulations: a delocalized TSE for the fully active enzyme. This feature would result in a lowering of the entropic barrier for the chemical step, thereby contributing to the enzyme-catalyzed rate enhancement. To avoid the issues of accuracy of QM/MM simulations that are well documented (Acevedo and Jorgensen, 2010; Elstner, 2007; Gaus et al., 2014; Roston et al., 2018), our system has the advantage of comparing the exact same reaction coordinate with the only difference being the presence/absence of a single atom, Mg²⁺, thereby studying differences rather than absolute values. Since the reaction with Mg²⁺ (fully active enzyme) revealed a more pronounced delocalization than the reaction in its absence, the divalent cation is predicted to lower the activation entropy. To test this prediction, we experimentally determined the enthalpic and entropic contributions to the chemical reaction barrier in Adk by measuring its temperature dependence in the presence and absence of Ca²⁺. The experimental trick of replacing Mg²⁺ with Ca²⁺ was used since the chemical step with Mg²⁺ is too fast to be experimentally measured, and previous studies showed that Ca²⁺ is an appropriate mimic to selectively probe the chemical step (Kerns et al., 2015).

The resulting Eyring plots for the two experiments (Figure 5A) deliver the enthalpic (ΔH^‡ of 15.3 ± 0.5 and 16.7 ± 1 kcal/mol) and entropic (ΔS^‡ = −5.7 ± 1.6 and −19.1 ± 3.1 cal/mol/K) contributions to the chemical reaction barrier in the presence and absence of Ca²⁺, respectively. The presence of Ca²⁺ reduces both barriers considerably. The decrease in the enthalpic barrier by Ca²⁺, which is even more pronounced with the optimal Mg²⁺ ion, is not new and has been well established in the literature for P-transfer reactions, including the Adk chemical reaction (Kerns et al., 2015; Pérez-Gallegos et al., 2017; Rosta et al., 2014; Yang et al., 2012). The striking new result, central to this work, is the large decrease in the entropic barrier in the presence of a divalent cation, thereby strongly supporting our key finding of a broader TSE from the QM/MM simulations.

Figure 5

Download asset Open asset

We were able to further test our QM/MM simulations with a second designed set of experiments by measuring the pH dependence of the P-transfer step. We find that the reaction rate increases with higher pH in the presence of Ca²⁺, whereas without a metal the opposite trend is observed (Figure 5B). These experimental results match the simulation results: the fully charged nucleotides state was the most reactive in the presence of the metal, whereas the monoprotonated state had a much higher free-energy barrier. In contrast, without metal only the monoprotonated nucleotides were reactive. Notably, at low pH (pH 5), the P-transfer step becomes rate limiting even with Mg²⁺, allowing for an estimate of the rate enhancement of Mg²⁺ versus Ca²⁺ of 24-fold.

Finally, we measured the forward and backward reaction rate constants of several arginine mutants that were identified as important for the P-transfer from the QM/MM simulations: R36K, R85K, R124K, R150K, and R161K (Table 2). The large decrease in the catalytic rate for each single Arg to Lys mutation further highlights the notion of a highly choreographed active site for a well-coordinated P-transfer (Figure 3G), in which arginine residues that are not even directly coordinated to the transferring phosphate play an equally important role.

Table 2

AAdk mutant	Rate constant (s−1)2 ADP → ATP + AMP	Rate constant (s−1)ATP + AMP → 2 ADP
R124K	4.6 ± 0.5	1.6 ± 0.4
R150K	1.3 ± 0.3	1.1 ± 0.3
R161K	0.4 ± 0.1	0.5 ± 0.2
R85K	0.3 ± 0.1	1.6 ± 0.2
R36K	0.8 ± 0.2	14 ± 2

Adk is an impressively efficient enzyme, being able to accelerate the P-transfer reaction over a billion times. Previous work showed the complexity of the Adk free-energy landscape for this two-substrate, two-product reaction (Figure 1A), providing direct experimental evidence for the widely accepted concept of a multidimension free-energy landscape with multiple intermediates and transitions states for the full enzymatic cycle, already described and demonstrated by Benkovic and Hammes-Schiffer (Benkovic et al., 2008). Full lid-closure results in an excellent preorientation of the donor and acceptor groups of the two bound nucleotides and the aiding active-site residues in the enzyme/substrate complex, in agreement with the general view of necessary pre-organization in enzyme catalysis. Lid-opening had been identified as the rate-limiting step in the enzymatic cycle (Kerns et al., 2015), and the corresponding TSE for this conformational change (Stiller et al., 2019) and the high-resolution structure of next intermediate, the EP complex before product dissociation, recently been described (Stiller et al., 2022).

For a comprehensive understanding of the complete free-energy landscape of Adk catalysis, we here investigate the key step, catalyzing the actual chemical step of P-transfer. While ES and EP complexes such as Adk bound to ADP/Mg²⁺ can be structurally characterized by traditional experimental methods as they represent a minimum in the free-energy landscape, the reaction path from substrate to product involving breaking and forming covalent bonds (chemical step) and traversing the crucial TS can only be ‘visualized’ by quantum mechanics-based molecular simulations. The power of such simulations to examine the chemical steps in enzymes, including P-transfer reactions, has been extensively documented (Jin et al., 2017; Lai and Cui, 2020b; Mokrushina et al., 2020; Pérez-Gallegos et al., 2015; Pérez-Gallegos et al., 2017; Roston and Cui, 2016; Valiev et al., 2007). The focus in the literature has been on the comparison of the enzyme-catalyzed and uncatalyzed reaction with respect to the bond character at the transitions state (i.e., associate versus dissociative or tight versus loose) (Admiraal and Herschlag, 1995; Hengge, 2002; Kamerlin et al., 2013; Lassila et al., 2011). QM calculations for uncatalyzed reactions have been well documented in the literature showing narrow and symmetric TSE (Klähn et al., 2006: Wang et al., 2015).

Here, we uncover a central new result that provokes a modified TS theory. Enzymes, due to their macromolecular nature, provide a fundamentally different, advantageous way to catalyze these chemical reactions compared to the uncatalyzed reaction by employing a broad TSE. Many different molecular configurations can be accommodated in the TS region with comparable energies via collective motions, spanning a 1-Å range along the reaction coordinate for the transferring phosphate that travels a total distance of only 2 Å during the P-transfer. This features resemblance the now well-established conformational sampling of proteins in ground states such as ES complexes. Our findings explain why enzymes do not face an entropic bottleneck for catalysis. Furthermore, as enzyme active sites are asymmetric in contrast to the symmetric nature of the solvent for uncatalyzed reactions, we find that the TSE comprises also highly asymmetric conformations. Our findings help to resolve the controversy about the nature of the TS in enzyme-catalyzed P-transfer reactions between theory and experiments (Lassila et al., 2011). The complex nature of the active site of enzymes, in contrast to simple solvent, results in different mechanisms in the enzyme-catalyzed reaction.

We note that a previous QM/MM minimum energy path calculations for a different Adk enzyme, E. coli Adk, using a semiempirical method had proposed a different mechanism with a stable metaphosphate intermediate, reporting an even lower energy for this metaphosphate intermediate than the ES and EP complexes (Shibanuma et al., 2020). In complex systems, such as enzymes, it is possible to observe artificial local minima when using minimum energy path searching strategies due to inadequate sampling (Mendieta-Moreno, 2015; Quesne et al., 2016). From our experimental NMR and X-ray data, we know that such a stable metaphosphate does not exist in Adk-catalyzed reactions, highlighting the importance of experimental verification of simulations as performed here, and the use of extensive sampling with proper thermodynamic treatment.

More recently, another group performed QM/MM simulations for E. coli Adk using the semiempirical AM1/d-phoT method (Nam et al., 2007) and umbrella sampling, and measured and computed rate reduction in the chemical step by similar Arg mutants as reported here (Ojeda-May et al., 2021; Dulko-Smith et al., 2023). In contrast to the paper described above (Shibanuma et al., 2020), their results are in agreement with our data. First, the experimentally measured rate reductions in Ojeda-May et al., 2021; Dulko-Smith et al., 2023 are in full agreement with our results. Second, our activation barriers from umbrella sampling are qualitatively similar to their values, with the obvious differences in the enzyme species (we study the thermophilic Aquifex Adk), and their choice of analyzing the QM/MM simulations in the opposite direction ATP/AMP to 2 ADP molecules. The latter lead them to singly protonate one oxygen on ATP, since this ‘back reaction’ is faster for the monoprotonated state, in agreement with our simulations. This nice agreement in both experiments and computation by our two groups strengthens the validity of reported results.

Intriguingly, the Nam/Wolf-Watz team (Ojeda-May et al., 2021; Dulko-Smith et al., 2023) and our group focus on completely different interpretations of the data! Nam/Wolf-Watz propose an interconnection, synchronization, between the chemical step and lid-opening. This language can be misleading since it can imply that the chemical step directly effects the opening, which of course cannot be the case. There is no ‘memory’ of the chemical step, the energy from the closing/opening step is dissipated within picoseconds. Lid-opening/closing is an energetically fully independent step from the P-transfer step, see our Figure 1A. We agree with the obvious result/interpretation by the authors that the Arg residues in the active site are important for both the P-transfer and the consequent lid-opening, but there is no synchronization between these independent steps. In contrast, we discover here a fundamental concept for rate enhancement by an optimal enzyme, the reduction in the activation entropy by a wide TSE. New experiments were triggered by our finding that then delivered experimental validation of this concept.

A wide TSE agrees with the definition of the TS being a surface on the potential energy surface instead of a saddle as elegantly discussed in a theoretical paper by Nussinov and colleagues (Ma et al., 2000). Our results deliver concrete evidence for this fundamental concept from QM/MM calculations and commitment analysis on Adk that were further buttressed by the consequently designed experiments. Although we have only quantitatively demonstrated this TSE concept here for Adk-catalyzed P-transfer, we hypothesize such mechanism to be more general for enzyme catalysis as it is rooted in the macroscopic nature of proteins. A conformationally spread-out TSE is likely one reason why enzymes are much bigger than a few active-site residues: On the one hand residues remote from the site of the chemical action allow for conformational flexibility without the protein unfolding, on the other hand they allow for increasing the probability of highly choreographed motions along the reaction coordinate (Hammes-Schiffer and Benkovic, 2006; Klinman and Kohen, 2013; Saen-Oon et al., 2008; Schramm and Schwartz, 2018) while reducing excursions into unproductive conformations (Otten et al., 2020). The latter mechanism is directly demonstrated here in the TSE, as the configurations seen in the simulations are along the P-transfer reaction coordinate. We note that we focus here on equilibrium effects and an atomistic description of the TSE, this is not related to the fundamental work on dynamical effects manifested in transmission coefficients due to barrier recrossing (Antoniou and Schwartz, 2016; Zinovjev and Tuñón, 2017). The latter has been quantitatively demonstrated to play a minor role.

A question often being asked is why would nature evolve a conformational step to be rate limiting, as the lid-opening step for Adk? We propose that this is the ‘price’ the enzyme pays for stabilizing the ES complex to enable efficient rate acceleration of the chemical step, and for suppressing detrimental alternate reactions such as phosphate hydrolysis. Release of product requires disassembly of such closed, low-energy structure. These principles have now been directly revealed on an atomic scale for Adk: extensive electrostatic interactions of five Arg and one Lys with the phosphoryl groups to efficiently accelerate the P-transfer step, which then need to be broken for lid-opening to allow product release. The effect of different Arg mutations on the lid-opening has been recently nicely demonstrated (Ojeda-May et al., 2021; Dulko-Smith et al., 2023).

As long as the rate of this conformational transition (k_open) and hence k_cat is not limiting for organismal fitness, there is no evolutionary pressure to further enhance it.

The importance of TSE including multiple transition pathways for protein folding (Frauenfelder et al., 2006; Royer, 2008), and more recently for protein conformational transitions within the folded state (Pontiggia et al., 2015; Stiller et al., 2019; Tsai and Nussinov, 2014) has been well documented. TSE in enzyme catalysis of chemical steps where covalent bonds are broken and formed as demonstrated here further promotes a unifying concept of protein folding and function, as the same principal concept of TSE is seen for these different types of conformational changes. Consequently, protein folding and function is embedded in a unifying energy landscape which is simply altered by ligands and solvent.

All the computational data produced in the current study is available in the following repository: ZENODO: https://doi.org/10.5281/zenodo.14647770.

1. 1. Jara G
   2. Marti MA
   3. Kern D
   4. Pontiggia F
   5. Otten R
   6. Agafonov R

   (2025) Zenodo

   Wide Transition-State Ensemble as Key Component for Enzyme Catalysis.

   https://doi.org/10.5281/zenodo.14647770

1. 1. Acevedo O
   2. Jorgensen WL

   (2010) Advances in quantum and molecular mechanical (QM/MM) simulations for organic and enzymatic reactions

   Accounts of Chemical Research 43:142–151.

   https://doi.org/10.1021/ar900171c

   • PubMed
   • Google Scholar
2. 1. Admiraal SJ
   2. Herschlag D

   (1995) Mapping the transition state for ATP hydrolysis: implications for enzymatic catalysis

   Chemistry & Biology 2:729–739.

   https://doi.org/10.1016/1074-5521(95)90101-9

   • Google Scholar
3. 1. Allen KN
   2. Dunaway-Mariano D

   (2016) Catalytic scaffolds for phosphoryl group transfer

   Current Opinion in Structural Biology 41:172–179.

   https://doi.org/10.1016/j.sbi.2016.07.017

   • PubMed
   • Google Scholar
4. 1. Allnér O
   2. Nilsson L
   3. Villa A

   (2012) Magnesium ion-water coordination and exchange in biomolecular simulations

   Journal of Chemical Theory and Computation 8:1493–1502.

   https://doi.org/10.1021/ct3000734

   • PubMed
   • Google Scholar
5. 1. Antoniou D
   2. Schwartz SD

   (2016) Phase space bottlenecks in enzymatic reactions

   The Journal of Physical Chemistry. B 120:433–439.

   https://doi.org/10.1021/acs.jpcb.5b11157

   • PubMed
   • Google Scholar
6. 1. Beckstein O
   2. Denning EJ
   3. Perilla JR
   4. Woolf TB

   (2009) Zipping and unzipping of adenylate kinase: atomistic insights into the ensemble of open<-->closed transitions

   Journal of Molecular Biology 394:160–176.

   https://doi.org/10.1016/j.jmb.2009.09.009

   • PubMed
   • Google Scholar
7. 1. Benkovic SJ
   2. Hammes GG
   3. Hammes-Schiffer S

   (2008) Free-energy landscape of enzyme catalysis

   Biochemistry 47:3317–3321.

   https://doi.org/10.1021/bi800049z

   • PubMed
   • Google Scholar
8. 1. Berendsen HJC
   2. Postma JPM
   3. van Gunsteren WF
   4. DiNola A
   5. Haak JR

   (1984) Molecular dynamics with coupling to an external bath

   The Journal of Chemical Physics 81:3684–3690.

   https://doi.org/10.1063/1.448118

   • Google Scholar
9. 1. Berry MB
   2. Meador B
   3. Bilderback T
   4. Liang P
   5. Glaser M
   6. Phillips GN

   (1994) The closed conformation of a highly flexible protein: the structure of E. coli adenylate kinase with bound AMP and AMPPNP

   Proteins 19:183–198.

   https://doi.org/10.1002/prot.340190304

   • PubMed
   • Google Scholar
10. 1. Berry MB
    2. Bae E
    3. Bilderback TR
    4. Glaser M
    5. Phillips GN

    (2006) Crystal structure of ADP/AMP complex of Escherichia coli adenylate kinase

    Proteins 62:555–556.

    https://doi.org/10.1002/prot.20699

    • PubMed
    • Google Scholar
11. 1. Bussi G
    2. Donadio D
    3. Parrinello M

    (2007) Canonical sampling through velocity rescaling

    The Journal of Chemical Physics 126:014101.

    https://doi.org/10.1063/1.2408420

    • PubMed
    • Google Scholar
12. 1. Cai ZL
    2. Lopez P
    3. Reimers JR
    4. Cui Q
    5. Elstner M

    (2007) Application of the computationally efficient self-consistent-charge density-functional tight-binding method to magnesium-containing molecules

    The Journal of Physical Chemistry. A 111:5743–5750.

    https://doi.org/10.1021/jp071701m

    • PubMed
    • Google Scholar
13. 1. Case DA
    2. Cheatham TE
    3. Darden T
    4. Gohlke H
    5. Luo R
    6. Merz KM
    7. Onufriev A
    8. Simmerling C
    9. Wang B
    10. Woods RJ

    (2005) The Amber biomolecular simulation programs

    Journal of Computational Chemistry 26:1668–1688.

    https://doi.org/10.1002/jcc.20290

    • PubMed
    • Google Scholar
14. 1. Cheng Y
    2. Zhang Y
    3. McCammon JA

    (2005) How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study

    Journal of the American Chemical Society 127:1553–1562.

    https://doi.org/10.1021/ja0464084

    • PubMed
    • Google Scholar
15. 1. Crespo A
    2. Martí MA
    3. Estrin DA
    4. Roitberg AE

    (2005) Multiple-steering QM-MM calculation of the free energy profile in chorismate mutase

    Journal of the American Chemical Society 127:6940–6941.

    https://doi.org/10.1021/ja0452830

    • PubMed
    • Google Scholar
16. 1. Cui Q

    (2016) Perspective: Quantum mechanical methods in biochemistry and biophysics

    The Journal of Chemical Physics 145:140901.

    https://doi.org/10.1063/1.4964410

    • PubMed
    • Google Scholar
17. 1. Defelipe LA
    2. Lanzarotti E
    3. Gauto D
    4. Marti MA
    5. Turjanski AG

    (2015) Protein topology determines cysteine oxidation fate: the case of sulfenyl amide formation among protein families

    PLOS Computational Biology 11:e1004051.

    https://doi.org/10.1371/journal.pcbi.1004051

    • PubMed
    • Google Scholar
18. 1. de M Seabra G
    2. Walker RC
    3. Elstner M
    4. Case DA
    5. Roitberg AE

    (2007) Implementation of the SCC-DFTB method for hybrid QM/MM simulations within the amber molecular dynamics package

    The Journal of Physical Chemistry. A 111:5655–5664.

    https://doi.org/10.1021/jp070071l

    • PubMed
    • Google Scholar
19. 1. Dinner AR
    2. Blackburn GM
    3. Karplus M

    (2001) Uracil-DNA glycosylase acts by substrate autocatalysis

    Nature 413:752–755.

    https://doi.org/10.1038/35099587

    • PubMed
    • Google Scholar
20. 1. Duarte F
    2. Åqvist J
    3. Williams NH
    4. Kamerlin SCL

    (2015) Resolving apparent conflicts between theoretical and experimental models of phosphate monoester hydrolysis

    Journal of the American Chemical Society 137:1081–1093.

    https://doi.org/10.1021/ja5082712

    • PubMed
    • Google Scholar
21. 1. Dulko-Smith B
    2. Ojeda-May P
    3. Ådén J
    4. Wolf-Watz M
    5. Nam K

    (2023) Mechanistic basis for a connection between the catalytic step and slow opening dynamics of adenylate kinase

    Journal of Chemical Information and Modeling 63:1556–1569.

    https://doi.org/10.1021/acs.jcim.2c01629

    • PubMed
    • Google Scholar
22. 1. Dzeja P
    2. Terzic A

    (2009) Adenylate kinase and AMP signaling networks: metabolic monitoring, signal communication and body energy sensing

    International Journal of Molecular Sciences 10:1729–1772.

    https://doi.org/10.3390/ijms10041729

    • PubMed
    • Google Scholar
23. 1. Elstner M

    (2007) SCC-DFTB: what is the proper degree of self-consistency?

    The Journal of Physical Chemistry. A 111:5614–5621.

    https://doi.org/10.1021/jp071338j

    • PubMed
    • Google Scholar
24. 1. Frauenfelder H
    2. Wolynes PG

    (1985) Rate theories and puzzles of hemeprotein kinetics

    Science 229:337–345.

    https://doi.org/10.1126/science.4012322

    • PubMed
    • Google Scholar
25. 1. Frauenfelder H
    2. Fenimore PW
    3. Chen G
    4. McMahon BH

    (2006) Protein folding is slaved to solvent motions

    PNAS 103:15469–15472.

    https://doi.org/10.1073/pnas.0607168103

    • PubMed
    • Google Scholar
26. 1. Ganguly A
    2. Weissman BP
    3. Giese TJ
    4. Li NS
    5. Hoshika S
    6. Rao S
    7. Benner SA
    8. Piccirilli JA
    9. York DM

    (2020) Confluence of theory and experiment reveals the catalytic mechanism of the Varkud satellite ribozyme

    Nature Chemistry 12:193–201.

    https://doi.org/10.1038/s41557-019-0391-x

    • PubMed
    • Google Scholar
27. 1. Gaus M
    2. Cui Q
    3. Elstner M

    (2014) Density functional tight binding: application to organic and biological molecules

    WIREs Computational Molecular Science 4:49–61.

    https://doi.org/10.1002/wcms.1156

    • Google Scholar
28. Software

    1. Grossfield A

    (2010) WHAM: an implementation of the weighted histogram analysis method, version 2.0.10.1

    WordPress.

    http://membrane.urmc.rochester.edu/wordpress/?page_id=126
29. 1. Hahn DK
    2. Tusell JR
    3. Sprang SR
    4. Chu X

    (2015) Catalytic mechanism of mammalian adenylyl cyclase: a computational investigation

    Biochemistry 54:6252–6262.

    https://doi.org/10.1021/acs.biochem.5b00655

    • PubMed
    • Google Scholar
30. 1. Hammes-Schiffer S
    2. Benkovic SJ

    (2006) Relating protein motion to catalysis

    Annual Review of Biochemistry 75:519–541.

    https://doi.org/10.1146/annurev.biochem.75.103004.142800

    • PubMed
    • Google Scholar
31. 1. Hengge AC

    (2002) Isotope effects in the study of phosphoryl and sulfuryl transfer reactions

    Accounts of Chemical Research 35:105–112.

    https://doi.org/10.1021/ar000143q

    • PubMed
    • Google Scholar
32. 1. Henzler-Wildman KA
    2. Thai V
    3. Lei M
    4. Ott M
    5. Wolf-Watz M
    6. Fenn T
    7. Pozharski E
    8. Wilson MA
    9. Petsko GA
    10. Karplus M
    11. Hübner CG
    12. Kern D

    (2007) Intrinsic motions along an enzymatic reaction trajectory

    Nature 450:838–844.

    https://doi.org/10.1038/nature06410

    • PubMed
    • Google Scholar
33. 1. Hirvonen VHA
    2. Mulholland AJ
    3. Spencer J
    4. van der Kamp MW

    (2020) Small changes in hydration determine cephalosporinase activity of OXA-48 β-Lactamases

    ACS Catalysis 10:6188–6196.

    https://doi.org/10.1021/acscatal.0c00596

    • Google Scholar
34. 1. Hornak V
    2. Abel R
    3. Okur A
    4. Strockbine B
    5. Roitberg A
    6. Simmerling C

    (2006) Comparison of multiple Amber force fields and development of improved protein backbone parameters

    Proteins 65:712–725.

    https://doi.org/10.1002/prot.21123

    • PubMed
    • Google Scholar
35. 1. Hou G
    2. Cui Q

    (2012) QM/MM analysis suggests that Alkaline Phosphatase (AP) and nucleotide pyrophosphatase/phosphodiesterase slightly tighten the transition state for phosphate diester hydrolysis relative to solution: implication for catalytic promiscuity in the AP superfamily

    Journal of the American Chemical Society 134:229–246.

    https://doi.org/10.1021/ja205226d

    • PubMed
    • Google Scholar
36. 1. Hou G
    2. Zhu X
    3. Elstner M
    4. Cui Q

    (2012) A modified QM/MM Hamiltonian with the Self-Consistent-Charge Density-Functional-Tight-Binding Theory for highly charged QM regions

    Journal of Chemical Theory and Computation 8:4293–4304.

    https://doi.org/10.1021/ct300649f

    • PubMed
    • Google Scholar
37. 1. Jarzynski C

    (1997) Nonequilibrium equality for free energy differences

    Physical Review Letters 78:2690–2693.

    https://doi.org/10.1103/PhysRevLett.78.2690

    • Google Scholar
38. 1. Jencks WP

    (1972) General acid-base catalysis of complex reactions in water

    Chemical Reviews 72:705–718.

    https://doi.org/10.1021/cr60280a004

    • Google Scholar
39. 1. Jin Y
    2. Richards NG
    3. Waltho JP
    4. Blackburn GM

    (2017) Metal fluorides as analogues for studies on phosphoryl transfer enzymes

    Angewandte Chemie 56:4110–4128.

    https://doi.org/10.1002/anie.201606474

    • PubMed
    • Google Scholar
40. 1. Jorgensen WL
    2. Chandrasekhar J
    3. Madura JD
    4. Impey RW
    5. Klein ML

    (1983) Comparison of simple potential functions for simulating liquid water

    The Journal of Chemical Physics 79:926–935.

    https://doi.org/10.1063/1.445869

    • Google Scholar
41. 1. Kamerlin SCL
    2. Wilkie J

    (2007) The role of metal ions in phosphate ester hydrolysis

    Organic & Biomolecular Chemistry 5:2098–2108.

    https://doi.org/10.1039/b701274h

    • PubMed
    • Google Scholar
42. 1. Kamerlin SCL
    2. Sharma PK
    3. Prasad RB
    4. Warshel A

    (2013) Why nature really chose phosphate

    Quarterly Reviews of Biophysics 46:1–132.

    https://doi.org/10.1017/S0033583512000157

    • PubMed
    • Google Scholar
43. 1. Karplus M
    2. Kuriyan J

    (2005) Molecular dynamics and protein function

    PNAS 102:6679–6685.

    https://doi.org/10.1073/pnas.0408930102

    • PubMed
    • Google Scholar
44. 1. Kerns SJ
    2. Agafonov RV
    3. Cho YJJ
    4. Pontiggia F
    5. Otten R
    6. Pachov DV
    7. Kutter S
    8. Phung LA
    9. Murphy PN
    10. Thai V
    11. Alber T
    12. Hagan MF
    13. Kern D

    (2015) The energy landscape of adenylate kinase during catalysis

    Nature Structural & Molecular Biology 22:124–131.

    https://doi.org/10.1038/nsmb.2941

    • PubMed
    • Google Scholar
45. 1. Kirby AJ
    2. Nome F

    (2015) Fundamentals of phosphate transfer

    Accounts of Chemical Research 48:1806–1814.

    https://doi.org/10.1021/acs.accounts.5b00072

    • PubMed
    • Google Scholar
46. 1. Klähn M
    2. Rosta E
    3. Warshel A

    (2006) On the mechanism of hydrolysis of phosphate monoesters dianions in solutions and proteins

    Journal of the American Chemical Society 128:15310–15323.

    https://doi.org/10.1021/ja065470t

    • PubMed
    • Google Scholar
47. 1. Klinman JP
    2. Kohen A

    (2013) Hydrogen tunneling links protein dynamics to enzyme catalysis

    Annual Review of Biochemistry 82:471–496.

    https://doi.org/10.1146/annurev-biochem-051710-133623

    • PubMed
    • Google Scholar
48. 1. Kumar S
    2. Rosenberg JM
    3. Bouzida D
    4. Swendsen RH
    5. Kollman PA

    (1992) THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. THE method

    Journal of Computational Chemistry 13:1011–1021.

    https://doi.org/10.1002/jcc.540130812

    • Google Scholar
49. 1. Lai R
    2. Cui Q

    (2020a) Differences in the nature of the phosphoryl transfer transition state in protein phosphatase 1 and alkaline phosphatase: insights from qm cluster models

    The Journal of Physical Chemistry. B 124:9371–9384.

    https://doi.org/10.1021/acs.jpcb.0c07863

    • PubMed
    • Google Scholar
50. 1. Lai R
    2. Cui Q

    (2020b) What does the brønsted slope measure in the phosphoryl transfer transition state?

    ACS Catalysis 10:13932–13945.

    https://doi.org/10.1021/acscatal.0c03764

    • PubMed
    • Google Scholar
51. 1. Lassila JK
    2. Zalatan JG
    3. Herschlag D

    (2011) Biological phosphoryl-transfer reactions: understanding mechanism and catalysis

    Annual Review of Biochemistry 80:669–702.

    https://doi.org/10.1146/annurev-biochem-060409-092741

    • PubMed
    • Google Scholar
52. 1. Lienhard GE

    (1973) Enzymatic catalysis and transition-state theory

    Science 180:149–154.

    https://doi.org/10.1126/science.180.4082.149

    • PubMed
    • Google Scholar
53. 1. López-Canut V
    2. Roca M
    3. Bertrán J
    4. Moliner V
    5. Tuñón I

    (2011) Promiscuity in alkaline phosphatase superfamily: unraveling evolution through molecular simulations

    Journal of the American Chemical Society 133:12050–12062.

    https://doi.org/10.1021/ja2017575

    • PubMed
    • Google Scholar
54. 1. Ma B
    2. Kumar S
    3. Tsai CJ
    4. Hu Z
    5. Nussinov R

    (2000) Transition-state ensemble in enzyme catalysis: possibility, reality, or necessity?

    Journal of Theoretical Biology 203:383–397.

    https://doi.org/10.1006/jtbi.2000.1097

    • PubMed
    • Google Scholar
55. 1. Masgrau L
    2. Truhlar DG

    (2015) The importance of ensemble averaging in enzyme kinetics

    Accounts of Chemical Research 48:431–438.

    https://doi.org/10.1021/ar500319e

    • PubMed
    • Google Scholar
56. 1. Meagher KL
    2. Redman LT
    3. Carlson HA

    (2003) Development of polyphosphate parameters for use with the AMBER force field

    Journal of Computational Chemistry 24:1016–1025.

    https://doi.org/10.1002/jcc.10262

    • PubMed
    • Google Scholar
57. 1. Mendieta-Moreno JI

    (2015) A practical quantum mechanics molecular mechanics method for the dynamical study of reactions in biomolecules

    In Advances in Protein Chemistry and Structural Biology 100:67–68.

    https://doi.org/10.1016/bs.apcsb.2015.06.003

    • Google Scholar
58. 1. Mokrushina YA
    2. Golovin AV
    3. Smirnov IV
    4. Chatziefthimiou SD
    5. Stepanova AV
    6. Bobik TV
    7. Zalevsky AO
    8. Zlobin AS
    9. Konovalov KA
    10. Terekhov SS
    11. Stepanov AV
    12. Pipiya SO
    13. Shamborant OG
    14. Round E
    15. Belogurov AA Jr
    16. Bourenkov G
    17. Makarov AA
    18. Wilmanns M
    19. Xie J
    20. Blackburn GM
    21. Gabibov AG
    22. Lerner RA

    (2020) Multiscale computation delivers organophosphorus reactivity and stereoselectivity to immunoglobulin scavengers

    PNAS 117:22841–22848.

    https://doi.org/10.1073/pnas.2010317117

    • PubMed
    • Google Scholar
59. 1. Mones L
    2. Tang WJ
    3. Florián J

    (2013) Empirical valence bond simulations of the chemical mechanism of ATP to cAMP conversion by anthrax edema factor

    Biochemistry 52:2672–2682.

    https://doi.org/10.1021/bi400088y

    • PubMed
    • Google Scholar
60. 1. Müller CW
    2. Schulz GE

    (1992) Structure of the complex between adenylate kinase from Escherichia coli and the inhibitor Ap5A refined at 1.9 A resolution: a model for A catalytic transition state

    Journal of Molecular Biology 224:159–177.

    https://doi.org/10.1016/0022-2836(92)90582-5

    • PubMed
    • Google Scholar
61. 1. Müller CW
    2. Schlauderer GJ
    3. Reinstein J
    4. Schulz GE

    (1996) Adenylate kinase motions during catalysis: an energetic counterweight balancing substrate binding

    Structure 4:147–156.

    https://doi.org/10.1016/s0969-2126(96)00018-4

    • PubMed
    • Google Scholar
62. 1. Nam K
    2. Cui Q
    3. Gao J
    4. York DM

    (2007) Specific reaction parametrization of the AM1/d hamiltonian for phosphoryl transfer reactions: H, O, and P Atoms

    Journal of Chemical Theory and Computation 3:486–504.

    https://doi.org/10.1021/ct6002466

    • PubMed
    • Google Scholar
63. 1. Nitsche MA
    2. Ferreria M
    3. Mocskos EE
    4. González Lebrero MC

    (2014) GPU accelerated implementation of density functional theory for hybrid QM/MM simulations

    Journal of Chemical Theory and Computation 10:959–967.

    https://doi.org/10.1021/ct400308n

    • PubMed
    • Google Scholar
64. 1. O’Ferrall RAM

    (1970) Relationships between E2 and E1cB mechanisms of β-elimination

    J. Chem. Soc. B pp. 274–277.

    https://doi.org/10.1039/J29700000274

    • Google Scholar
65. 1. Ojeda-May P
    2. Mushtaq AU
    3. Rogne P
    4. Verma A
    5. Ovchinnikov V
    6. Grundström C
    7. Dulko-Smith B
    8. Sauer UH
    9. Wolf-Watz M
    10. Nam K

    (2021) Dynamic connection between enzymatic catalysis and collective protein motions

    Biochemistry 60:2246–2258.

    https://doi.org/10.1021/acs.biochem.1c00221

    • PubMed
    • Google Scholar
66. 1. Otten R
    2. Pádua RAP
    3. Bunzel HA
    4. Nguyen V
    5. Pitsawong W
    6. Patterson M
    7. Sui S
    8. Perry SL
    9. Cohen AE
    10. Hilvert D
    11. Kern D

    (2020) How directed evolution reshapes the energy landscape in an enzyme to boost catalysis

    Science 370:1442–1446.

    https://doi.org/10.1126/science.abd3623

    • PubMed
    • Google Scholar
67. 1. Ozer G
    2. Valeev EF
    3. Quirk S
    4. Hernandez R

    (2010) Adaptive steered molecular dynamics of the long-distance unfolding of neuropeptide Y

    Journal of Chemical Theory and Computation 6:3026–3038.

    https://doi.org/10.1021/ct100320g

    • PubMed
    • Google Scholar
68. 1. Ozer G
    2. Quirk S
    3. Hernandez R

    (2012) Thermodynamics of decaalanine stretching in water obtained by adaptive steered molecular dynamics simulations

    Journal of Chemical Theory and Computation 8:4837–4844.

    https://doi.org/10.1021/ct300709u

    • PubMed
    • Google Scholar
69. 1. Pabis A
    2. Duarte F
    3. Kamerlin SCL

    (2016) Promiscuity in the enzymatic catalysis of phosphate and sulfate transfer

    Biochemistry 55:3061–3081.

    https://doi.org/10.1021/acs.biochem.6b00297

    • PubMed
    • Google Scholar
70. 1. Palermo G
    2. Cavalli A
    3. Klein ML
    4. Alfonso-Prieto M
    5. Dal Peraro M
    6. De Vivo M

    (2015) Catalytic metal ions and enzymatic processing of DNA and RNA

    Accounts of Chemical Research 48:220–228.

    https://doi.org/10.1021/ar500314j

    • PubMed
    • Google Scholar
71. 1. Park S
    2. Khalili-Araghi F
    3. Tajkhorshid E
    4. Schulten K

    (2003) Free energy calculation from steered molecular dynamics simulations using Jarzynski’s equality

    The Journal of Chemical Physics 119:3559–3566.

    https://doi.org/10.1063/1.1590311

    • Google Scholar
72. 1. Perdew JP
    2. Burke K
    3. Ernzerhof M

    (1996) Generalized gradient approximation made simple

    Physical Review Letters 77:3865–3868.

    https://doi.org/10.1103/PhysRevLett.77.3865

    • PubMed
    • Google Scholar
73. 1. Pérez-Gallegos A
    2. Garcia-Viloca M
    3. González-Lafont À
    4. Lluch JM

    (2015) SP20 phosphorylation reaction catalyzed by protein kinase a: qm/mm calculations based on recently determined crystallographic structures

    ACS Catalysis 5:4897–4912.

    https://doi.org/10.1021/acscatal.5b01064

    • Google Scholar
74. 1. Pérez-Gallegos A
    2. Garcia-Viloca M
    3. González-Lafont À
    4. Lluch JM

    (2017) Understanding how cAMP-dependent protein kinase can catalyze phosphoryl transfer in the presence of Ca 2+ and Sr 2+: A QM/MM study

    Physical Chemistry Chemical Physics 19:10377–10394.

    https://doi.org/10.1039/C7CP00666G

    • Google Scholar
75. 1. Pontiggia F
    2. Pachov DV
    3. Clarkson MW
    4. Villali J
    5. Hagan MF
    6. Pande VS
    7. Kern D

    (2015) Free energy landscape of activation in a signalling protein at atomic resolution

    Nature Communications 6:7284.

    https://doi.org/10.1038/ncomms8284

    • PubMed
    • Google Scholar
76. 1. Quesne MG
    2. Borowski T
    3. de Visser SP

    (2016) Quantum mechanics/molecular mechanics modeling of enzymatic processes: caveats and breakthroughs

    Chemistry 22:2562–2581.

    https://doi.org/10.1002/chem.201503802

    • PubMed
    • Google Scholar
77. 1. Ramírez CL
    2. Zeida A
    3. Jara GE
    4. Roitberg AE
    5. Martí MA

    (2014) Improving efficiency in smd simulations through a hybrid differential relaxation algorithm

    Journal of Chemical Theory and Computation 10:4609–4617.

    https://doi.org/10.1021/ct500672d

    • PubMed
    • Google Scholar
78. 1. Rosta E
    2. Yang W
    3. Hummer G

    (2014) Calcium inhibition of ribonuclease H1 two-metal ion catalysis

    Journal of the American Chemical Society 136:3137–3144.

    https://doi.org/10.1021/ja411408x

    • PubMed
    • Google Scholar
79. 1. Roston D
    2. Cui Q

    (2016) Substrate and transition state binding in alkaline phosphatase analyzed by computation of oxygen isotope effects

    Journal of the American Chemical Society 138:11946–11957.

    https://doi.org/10.1021/jacs.6b07347

    • PubMed
    • Google Scholar
80. 1. Roston D
    2. Lu X
    3. Fang D
    4. Demapan D
    5. Cui Q

    (2018) Analysis of phosphoryl-transfer enzymes with qm/mm free energy simulations

    Methods in Enzymology 607:53–90.

    https://doi.org/10.1016/bs.mie.2018.05.005

    • PubMed
    • Google Scholar
81. 1. Royer CA

    (2008) The nature of the transition state ensemble and the mechanisms of protein folding: a review

    Archives of Biochemistry and Biophysics 469:34–45.

    https://doi.org/10.1016/j.abb.2007.08.022

    • PubMed
    • Google Scholar
82. 1. Ryckaert JP
    2. Ciccotti G
    3. Berendsen HJC

    (1977) Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes

    Journal of Computational Physics 23:327–341.

    https://doi.org/10.1016/0021-9991(77)90098-5

    • Google Scholar
83. 1. Saen-Oon S
    2. Ghanem M
    3. Schramm VL
    4. Schwartz SD

    (2008) Remote mutations and active site dynamics correlate with catalytic properties of purine nucleoside phosphorylase

    Biophysical Journal 94:4078–4088.

    https://doi.org/10.1529/biophysj.107.121913

    • PubMed
    • Google Scholar
84. 1. Schramm VL

    (2007) Enzymatic transition state theory and transition state analogue design

    The Journal of Biological Chemistry 282:28297–28300.

    https://doi.org/10.1074/jbc.R700018200

    • PubMed
    • Google Scholar
85. 1. Schramm VL
    2. Schwartz SD

    (2018) Promoting vibrations and the function of enzymes

    Emerging Theoretical and Experimental Convergence. Biochemistry 57:3299–3308.

    https://doi.org/10.1021/acs.biochem.8b00201

    • Google Scholar
86. 1. Schwartz SD
    2. Schramm VL

    (2009) Enzymatic transition states and dynamic motion in barrier crossing

    Nature Chemical Biology 5:551–558.

    https://doi.org/10.1038/nchembio.202

    • PubMed
    • Google Scholar
87. 1. Senn HM
    2. Thiel W

    (2009) QM/MM methods for biomolecular systems

    Angewandte Chemie 48:1198–1229.

    https://doi.org/10.1002/anie.200802019

    • PubMed
    • Google Scholar
88. 1. Shibanuma Y
    2. Nemoto N
    3. Yamamoto N
    4. Sampei GI
    5. Kawai G

    (2020) Crystal structure of adenylate kinase from an extremophilic archaeon Aeropyrum pernix with ATP and AMP

    Journal of Biochemistry 168:223–229.

    https://doi.org/10.1093/jb/mvaa043

    • PubMed
    • Google Scholar
89. 1. Stiller JB
    2. Kerns SJ
    3. Hoemberger M
    4. Cho YJ
    5. Otten R
    6. Hagan MF
    7. Kern D

    (2019) Probing the transition state in enzyme catalysis by high-pressure nmr dynamics

    Nature Catalysis 2:726–734.

    https://doi.org/10.1038/s41929-019-0307-6

    • PubMed
    • Google Scholar
90. 1. Stiller JB
    2. Otten R
    3. Häussinger D
    4. Rieder PS
    5. Theobald DL
    6. Kern D

    (2022) Structure determination of high-energy states in a dynamic protein ensemble

    Nature 603:528–535.

    https://doi.org/10.1038/s41586-022-04468-9

    • PubMed
    • Google Scholar
91. 1. Stockbridge RB
    2. Wolfenden R

    (2009) The intrinsic reactivity of ATP and the catalytic proficiencies of kinases acting on glucose, N-acetylgalactosamine, and homoserine: a thermodynamic analysis

    The Journal of Biological Chemistry 284:22747–22757.

    https://doi.org/10.1074/jbc.M109.017806

    • PubMed
    • Google Scholar
92. 1. Truhlar DG

    (2015) Transition state theory for enzyme kinetics

    Archives of Biochemistry and Biophysics 582:10–17.

    https://doi.org/10.1016/j.abb.2015.05.004

    • PubMed
    • Google Scholar
93. 1. Tsai CJ
    2. Nussinov R

    (2014) A unified view of “how allostery works”

    PLOS Computational Biology 10:e1003394.

    https://doi.org/10.1371/journal.pcbi.1003394

    • PubMed
    • Google Scholar
94. 1. Turjanski AG
    2. Hummer G
    3. Gutkind JS

    (2009) How mitogen-activated protein kinases recognize and phosphorylate their targets: A QM/MM study

    Journal of the American Chemical Society 131:6141–6148.

    https://doi.org/10.1021/ja8071995

    • PubMed
    • Google Scholar
95. 1. Valiev M
    2. Yang J
    3. Adams JA
    4. Taylor SS
    5. Weare JH

    (2007) Phosphorylation reaction in cAPK protein kinase-free energy quantum mechanical/molecular mechanics simulations

    The Journal of Physical Chemistry. B 111:13455–13464.

    https://doi.org/10.1021/jp074853q

    • PubMed
    • Google Scholar
96. 1. van der Kamp MW
    2. Mulholland AJ

    (2013) Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

    Biochemistry 52:2708–2728.

    https://doi.org/10.1021/bi400215w

    • PubMed
    • Google Scholar
97. 1. Walker RC
    2. Crowley MF
    3. Case DA

    (2008) The implementation of a fast and accurate QM/MM potential method in Amber

    Journal of Computational Chemistry 29:1019–1031.

    https://doi.org/10.1002/jcc.20857

    • PubMed
    • Google Scholar
98. 1. Wang C
    2. Huang W
    3. Liao JL

    (2015) QM/MM investigation of ATP hydrolysis in aqueous solution

    The Journal of Physical Chemistry. B 119:3720–3726.

    https://doi.org/10.1021/jp512960e

    • PubMed
    • Google Scholar
99. 1. Warshel A
    2. Bora RP

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

    The Journal of Chemical Physics 144:180901.

    https://doi.org/10.1063/1.4947037

    • PubMed
    • Google Scholar
100. 1. Westheimer FH

     (1987) Why nature chose phosphates

     Science 235:1173–1178.

     https://doi.org/10.1126/science.2434996

     • PubMed
     • Google Scholar
101. 1. Wolf-Watz M
     2. Thai V
     3. Henzler-Wildman K
     4. Hadjipavlou G
     5. Eisenmesser EZ
     6. Kern D

     (2004) Linkage between dynamics and catalysis in a thermophilic-mesophilic enzyme pair

     Nature Structural & Molecular Biology 11:945–949.

     https://doi.org/10.1038/nsmb821

     • PubMed
     • Google Scholar
102. 1. Yang Y
     2. Yu H
     3. York D
     4. Elstner M
     5. Cui Q

     (2008) Description of phosphate hydrolysis reactions with the Self-Consistent-Charge Density-Functional-Tight-Binding (SCC-DFTB) theory. 1. Parameterization

     Journal of Chemical Theory and Computation 4:2067–2084.

     https://doi.org/10.1021/ct800330d

     • PubMed
     • Google Scholar
103. 1. Yang L
     2. Liao RZ
     3. Ding WJ
     4. Liu K
     5. Yu JG
     6. Liu RZ

     (2012) Why calcium inhibits magnesium-dependent enzyme phosphoserine phosphatase? A theoretical study

     Theoretical Chemistry Accounts 131:1–12.

     https://doi.org/10.1007/s00214-012-1275-y

     • Google Scholar
104. 1. Zinovjev K
     2. Tuñón I

     (2017) Quantifying the limits of transition state theory in enzymatic catalysis

     PNAS 114:12390–12395.

     https://doi.org/10.1073/pnas.1710820114

     • PubMed
     • Google Scholar

Author details

1. Gabriel E Jara

   Departamento de Química Inorgánica, Analítica y Química-Física (INQUIMAE-CONICET), Universidad de Buenos Aires, Buenos Aires, Argentina

   Present address

   Brazilian Biosciences National Laboratory (LNBio), Brazilian Center for Research in Energy and Materials (CNPEM), Campinas, Brazil

   Contribution

   Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-5831-1392

2. Francesco Pontiggia

   Howard Hughes Medical Institute, Department of Biochemistry, Brandeis University, Waltham, United States

   Present address

   Psivant Therapeutics, Salem, United States

   Contribution

   Data curation, Formal analysis, Investigation, Visualization

   Competing interests

   No competing interests declared

3. Renee Otten

   Howard Hughes Medical Institute, Department of Biochemistry, Brandeis University, Waltham, United States

   Present address

   Treeline Biosciences, Watertown, United States

   Contribution

   Data curation, Formal analysis, Investigation, Visualization

   Competing interests

   No competing interests declared

4. Roman V Agafonov

   Howard Hughes Medical Institute, Department of Biochemistry, Brandeis University, Waltham, United States

   Present address

   C4 Therapeutics, Watertown, United States

   Contribution

   Data curation, Formal analysis, Investigation, Visualization

   Competing interests

   No competing interests declared

5. Marcelo A Martí

   Departamento de Química Biológica (IQUIBICEN-CONICET), Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina

   Contribution

   Conceptualization, Resources, Supervision, Investigation, Writing – original draft, Project administration, Writing – review and editing

   For correspondence

   marti.marcelo@gmail.com

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7911-9340

6. Dorothee Kern

   Howard Hughes Medical Institute, Department of Biochemistry, Brandeis University, Waltham, United States

   Present address

   Department of Integrative Structural & Computational Biology and HHMI, The Scripps Research Institute, La Jolla, United States

   Contribution

   Conceptualization, Resources, Supervision, Funding acquisition, Investigation, Methodology, Writing – original draft, Project administration, Writing – review and editing

   For correspondence

   dkern@brandeis.edu

   Competing interests

   D.K. is co-founder of Relay Therapeutics and MOMA Therapeutics

   "This ORCID iD identifies the author of this article:" 0000-0002-7631-8328

• 710

  views

• 39

  downloads

• 0

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.