Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective

The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.

Factors That Determine the Variation of Equilibrium and Kinetic Properties of QM/MM Enzyme Simulations: QM Region, Conformation, and Boundary Condition

To analyze the impact of various technical details on the results of quantum mechanical (QM)/molecular mechanical (MM) enzyme simulations, including the QM region size, catechol-O-methyltransferase (COMT) is studied as a model system using an approximate QM/MM method (DFTB3/CHARMM). The results show that key equilibrium and kinetic properties for methyl transfer in COMT exhibit limited variations with respect to the size of the QM region, which ranges from ∼100 to ∼500 atoms in this study. With extensive sampling, local and global structural characteristics of the enzyme are largely conserved across the studied QM regions, while the nature of the transition state (e.g., secondary kinetic isotope effect) and reaction exergonicity are largely maintained. Deviations in the free energy profile with different QM region sizes are similar in magnitude to those observed with changes in other simulation protocols, such as different initial enzyme conformations and boundary conditions. Electronic structural properties, such as the covariance matrix of residual charge fluctuations, appear to exhibit rather long-range correlations, especially when the peptide backbone is included in the QM region; this observation holds when a range-separated DFT approach is used as the QM region, suggesting that delocalization error is unlikely the origin. Overall, the analyses suggest that multiple simulation details determine the results of QM/MM enzyme simulations with comparable contributions.

Monitoring protein conformational changes using fluorescent nanoantennas

Understanding the relationship between protein structural dynamics and function is crucial for both basic research and biotechnology. However, methods for studying the fast dynamics of structural changes are limited. Here, we introduce fluorescent nanoantennas as a spectroscopic technique to sense and report protein conformational changes through noncovalent dye-protein interactions. Using experiments and molecular simulations, we detect and characterize five distinct conformational states of intestinal alkaline phosphatase, including the transient enzyme-substrate complex. We also explored the universality of the nanoantenna strategy with another model protein, Protein G and its interaction with antibodies, and demonstrated a rapid screening strategy to identify efficient nanoantennas. These versatile nanoantennas can be used with diverse dyes to monitor small and large conformational changes, suggesting that they could be used to characterize diverse protein movements or in high-throughput screening applications.

Machine Learning for Chemical Reactions

Machine learning (ML) techniques applied to chemical reactions have a long history. The present contribution discusses applications ranging from small molecule reaction dynamics to computational platforms for reaction planning. ML-based techniques can be particularly relevant for problems involving both computation and experiments. For one, Bayesian inference is a powerful approach to develop models consistent with knowledge from experiments. Second, ML-based methods can also be used to handle problems that are formally intractable using conventional approaches, such as exhaustive characterization of state-to-state information in reactive collisions. Finally, the explicit simulation of reactive networks as they occur in combustion has become possible using machine-learned neural network potentials. This review provides an overview of the questions that can and have been addressed using machine learning techniques, and an outlook discusses challenges in this diverse and stimulating field. It is concluded that ML applied to chemistry problems as practiced and conceived today has the potential to transform the way with which the field approaches problems involving chemical reactions, in both research and academic teaching.

Structures of an engineered phospholipase D with specificity for secondary alcohol transphosphatidylation: insights into plasticity of substrate binding and activation

Phospholipase D (PLD) is an enzyme useful for the enzymatic modification of phospholipids. In the presence of primary alcohols, the enzyme catalyses transphosphatidylation of the head group of phospholipid substrates to synthesise a modified phospholipid product. However, the enzyme is specific for primary alcohols and thus the limitation of the molecular size of the acceptor compounds has restricted the type of phospholipid species that can be synthesised. An engineered variant of PLD from Streptomyces antibioticus termed TNYR SaPLD was developed capable of synthesising 1-phosphatidylinositol with positional specificity of up to 98%. To gain a better understanding of the substrate binding features of the TNYR SaPLD, crystal structures have been determined for the free enzyme and its complexes with phosphate, phosphatidic acid and 1-inositol phosphate. Comparisons of these structures with the wild-type SaPLD show a larger binding site able to accommodate a bulkier secondary alcohol substrate as well as changes to the position of a flexible surface loop proposed to be involved in substrate recognition. The complex of the active TNYR SaPLD with 1-inositol phosphate reveals a covalent intermediate adduct with the ligand bound to H442 rather than to H168, the proposed nucleophile in the wild-type enzyme. This structural feature suggests that the enzyme exhibits plasticity of the catalytic mechanism different from what has been reported to date for PLDs. These structural studies provide insights into the underlying mechanism that governs the recognition of myo-inositol by TNYR SaPLD, and an important foundation for further studies of the catalytic mechanism.

Alkaline Phosphatases: in Silico Study on the Catalytic Effect of Conserved Active Site Residues Using Human Placental Alkaline Phosphatase (PLAP) As a Model Protein

The metalloenzymes from the alkaline phosphatase (AP) superfamily catalyze the hydrolysis and transphosphorylation of phosphate monoesters. The role of several amino acids highly conserved in the active site of this family of enzymes was examined, using human placental AP (PLAP) as a model protein. By employing an active-site model based on the X-ray crystal structure of PLAP, mutations of several key residues were modeled by quantum mechanical methods in order to determine their impact on the catalytic activity. Kinetic and thermodynamic estimations were achieved for each reaction step of the catalytic mechanism by characterization of the intermediates and transition states on the reaction pathway, and the effects of mutations on the activation barriers were analyzed. A good accordance was observed between the present computational results and experimental measurements reported in the literature.

What Does the Brønsted Slope Measure in the Phosphoryl Transfer Transition State?

The structural and energetic features of phosphate and phosphonate hydrolysis in Protein Phosphatase-1 (PP1) and water are studied using quantum mechanical (QM) cluster models. The calculations are able to reproduce observed kinetic isotope effects and capture several key trends in the experimental Brønsted plots: the β l g values are rather different for phosphate and phosphonate ester hydrolysis in solution but are similar in PP1. Detailed analyses of structure, charge distribution and bond order of computed transition states support the general conclusion from experimental study that phosphoryl transfer transition states are different for the two classes of substrates in solution but similar in PP1. On the other hand, the microscopic models also highlight notable differences between the phosphate and phosphonate transition states, which are manifested in not only structure but also kinetic isotope effects. Overall, we find that while β l g / β E Q , l g generally correlates with the partial charge on leaving group oxygen and the fractional bond order of the breaking P- O l g bond, the precise mapping between β l g / β E Q , l g and P- O l g bond order in the transition state is difficult due largely to the cross talk between breaking and forming P-O bonds. Therefore, further supporting previous analyses of limitations of free energy relations, our results suggest that while free energy relation is a valuable tool for probing the nature of transition state, a quantitative mapping of β l g and β l g / β E Q , l g values to structure or charge in the transition state should be conducted with great care.

Differences in the Nature of the Phosphoryl Transfer Transition State in Protein Phosphatase 1 and Alkaline Phosphatase: Insights from QM Cluster Models

Quantum mechanical (QM) cluster models are used to probe effects on the catalytic properties of protein phosphatase 1 (PP1) and alkaline phosphatase (AP) due to metal ions and active site residues. The calculations suggest that the phosphoryl transfer transition states in PP1 are synchronous in nature with a significant degree of P-Olg cleavage, while those in AP are tighter with a modest degree of P-Olg cleavage and a range of P-Onuc formation. Similar to observations made in our recent work, a significant degree of cross talk between the forming and breaking P-O bonds complicates the interpretation of the Brønsted relation, especially in regard to AP for which the computed βlg/βEQ,lg value does not correlate with the degree of P-Olg cleavage regardless of the metal ions in the active site. By comparison, the correlation between βlg/βEQ,lg and the P-Olg bond order is more applicable to PP1, which generally exhibits less variation in the transition state than AP. Results for computational models with swapped metal ions between PP1 and AP suggest that the metal ions modulate both the nature of the transition state and the degrees of sensitivity of the transition state to the leaving group. In the reactant state, the degree of the scissile bond polarization is also different in the two enzymes, although this difference appears to be largely determined by the active site residues rather than the metal ions. Therefore, both the identity of the metal ion and the positioning of polar or charged residues in the active site contribute to the distinct catalytic characteristics of these enzymes. Several discrepancies observed between the QM cluster results and the available experimental data highlight the need for further QM/MM method developments for the quantitative analysis of metalloenzymes that contain open-shell transition metal ions.

Single-step Replacement of an Unreactive C-H Bond by a C-S Bond Using Polysulfide as the Direct Sulfur Source in Anaerobic Ergothioneine Biosynthesis

Ergothioneine, a natural longevity vitamin and antioxidant, is a thiol-histidine derivative. Recently, two types of biosynthetic pathways were reported. In the aerobic ergothioneine biosynthesis, a non-heme iron enzyme incorporates a sulfoxide to an sp² C-H bond in trimethyl-histidine (hercynine) through oxidation reactions. In contrast, in the anaerobic ergothioneine biosynthetic pathway in a green sulfur bacterium, Chlorobium limicola, a rhodanese domain containing protein (EanB) directly replaces this unreactive hercynine C-H bond with a C-S bond. Herein, we demonstrate that polysulfide (HSS_nSR) is the direct sulfur-source in EanB-catalysis. After identifying EanB's substrates, X-ray crystallography of several intermediate states along with mass spectrometry results provide additional mechanistic details for this reaction. Further, quantum mechanics/molecular mechanics (QM/MM) calculations reveal that protonation of N_π of hercynine by Tyr353 with the assistance of Thr414 is a key activation step for the hercynine sp² C-H bond in this trans-sulfuration reaction.

One- and Two-Proton Transfer Mechanisms Coexist in One Active Site

Acibenzolar-S-methyl (ASM) is one of the most successfully commercialized plant activators of the systemic acquired resistance (SAR). However, its activation (hydrolysis) mechanism catalyzed by the salicylic acid binding protein 2 (SABP2) remains elusive. The fundamental catalytic mechanism of the SABP2-catalyzed hydrolysis of the ASM had been investigated by extensive computational and experimental studies, including QM/MM simulations, charge transfer analysis, small-molecule synthesis, and biochemical assays. Here we report that the promiscuous SABP2 shows different catalytic mechanisms toward different substrates. To catalyze the ASM hydrolysis, the SABP2 uses a two-proton transfer mechanism, and the key intermediate is stabilized by the charge transfer effect; to catalyze the ethyl 1,2,3-benzothiadiazole-7-carboxylate (BTM, an ASM analogue) hydrolysis, the SABP2 applies the one-proton transfer mechanism, and the classic tetrahedral intermediate is stabilized by the electrostatic effect. The HPLC analyses of the SABP2 esterase activities toward the ASM and the BTM show comparable results with our computaional results, suggesting that the obtained computational mechanism insights are reasonable. The obtained mechanism is not only an important supplement to the theory of enzymes' catalytic promiscuity, but it also contributes a possible strategy for the design of next generation plant SAR activators.

Insights into Catalytic Gas-Phase Hydrolysis of Organophosphate Chemical Warfare Agents by MOF-Supported Bimetallic Metal-Oxo Clusters

Zirconium-based metal-organic frameworks (Zr-MOFs) have been reported to be efficient catalysts for the hydrolysis of organophosphate chemical warfare agents (CWAs) in buffered solutions. However, for the gas-phase reaction, which is more relevant to the situation in a battlefield gas mask application, the kinetics of Zr-MOF catalysts may be severely hindered by strong product inhibition. To improve the catalytic performance, we computationally screened a series of synthetically accessible Zr-MOF-supported bimetallic metal-oxo clusters in which the metal-oxygen-metal active motif is preserved, aiming to find catalysts that have lower binding affinities to the hydrolysis product. For the promising catalyst Al₂O₂(OH)₂@NU-1000 identified from the screening using density functional theory, we mapped out the full reaction pathway of gas-phase dimethyl p-nitrophenolphosphate (DMNP) hydrolysis and analyzed the free energy profile as well as the turnover frequency (TOF). We found that the catalytic mechanism on the new catalyst is slightly different from the one on NU-1000, which also led to a different TOF-limiting step. Additional factors that can affect the overall catalytic performance in practical application, such as the amount of ambient moisture and the existence of acid gases that may poison the catalyst, have also been evaluated.

Extensive free-energy simulations identify water as the base in nucleotide addition by DNA polymerase

Transphosphorylation of nucleotide triphosphates is the central reaction in DNA replication by DNA polymerase as well as many other biological processes. Despite its importance, the microscopic chemical mechanism of transphosphorylation of nucleotide triphosphates is, in most cases, unknown. Here we use extensive simulations of DNA polymerase η to test mechanistic hypotheses. We systematically survey the reactive space by calculating 2D free-energy surfaces for 10 different plausible mechanisms that have been proposed. We supplement these free-energy surfaces with calculations of pKa for a number of potentially acidic protons in different states relevant to the catalytic cycle. We find that among all of the conditions that we test, the smallest activation barrier occurs for a reaction where a Mg2+-coordinated water deprotonates the nucleophilic 3'-OH, and this deprotonation is concerted with the phosphoryl transfer. The presence of a third Mg2+ in the active site lowers the activation barrier for the water-as-base mechanism, as does protonation of the pyrophosphate leaving group, which is consistent with general acid catalysis. The results demonstrate the value of simulations, when used in conjunction with experimental data, to help establish a microscopic chemical mechanism in a complex environment.

Mechanistic insights into the phosphoryl transfer reaction in cyclin-dependent kinase 2: A QM/MM study

Cyclin-dependent kinase 2 (CDK2) is an important member of the CDK family exerting its most important function in the regulation of the cell cycle. It catalyzes the transfer of the gamma phosphate group from an ATP (adenosine triphosphate) molecule to a Serine/Threonine residue of a peptide substrate. Due to the importance of this enzyme, and protein kinases in general, a detailed understanding of the reaction mechanism is desired. Thus, in this work the phosphoryl transfer reaction catalyzed by CDK2 was revisited and studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our results suggest that the base-assisted mechanism is preferred over the substrate-assisted pathway when one Mg²⁺ is present in the active site, in agreement with a previous theoretical study. The base-assisted mechanism resulted to be dissociative, with a potential energy barrier of 14.3 kcal/mol, very close to the experimental derived value. An interesting feature of the mechanism is the proton transfer from Lys129 to the phosphoryl group at the second transition state, event that could be helping in neutralizing the charge on the phosphoryl group upon the absence of a second Mg²⁺ ion. Furthermore, important insights into the mechanisms in terms of bond order and charge analysis were provided. These descriptors helped to characterize the synchronicity of bond forming and breaking events, and to characterize charge transfer effects. Local interactions at the active site are key to modulate the charge distribution on the phosphoryl group and therefore alter its reactivity.

The catalytic mechanism of S-acyltransferases: acylation is triggered on by a loose transition state and deacylation is turned off by a tight transition state

The reaction of S-acyltransferase is characterized by a loose transition state.

Two symmetric arginine residues play distinct roles in Thermus thermophilus Argonaute DNA guide strand-mediated DNA target cleavage

Bacterium Thermus thermophilus Argonaute (Ago; TtAgo) is a prokaryotic Ago (pAgo) that acts as the host defense against the uptake and propagation of foreign DNA by catalyzing the DNA cleavage reaction. The TtAgo active site consists of a plugged-in glutamate finger with two arginine residues (R545 and R486) located symmetrically around it. An interesting challenge is to understand how they can collaboratively facilitate enzymatic catalysis. In Kluyveromyces polysporus Ago, a eukaryotic Ago, the evolutionarily symmetrical residues are arginine and histidine, both of which function to stabilize the plugged-in catalytic tetrad conformation. Surprisingly, our simulation results indicated that, in TtAgo, only R545 is involved in the cleavage reaction by serving as a critical structural anchor to stabilize the catalytic tetrad Asp-Glu-Asp-Asp that is completed by the insertion of the glutamate finger, whereas R486 is not involved in target cleavage. The TtAgo-mediated target DNA cleavage occurs in a substrate-assisted mechanism, in which the pro-Rp (Rp, a tetrahedral phosphorus center with "R-type" chirality) oxygen of scissile phosphate acts as a general base to activate the nucleophilic water. Our unexpected theoretical findings on distinct roles played by R545 and R486 in TtAgo catalysis have been validated by single-point site-mutagenesis experiments, wherein the target cleavage is abolished for all mutants of R545. In sharp contrast, the cleavage activity is maintained for all mutants of R486. Our work provides mechanistic insights on the catalytic specificity of Ago proteins and could facilitate the design of new gene-editing tools in the long term.

Structural and mechanistic basis for preferential deadenylation of U6 snRNA by Usb1

Post-transcriptional modification of snRNA is central to spliceosome function. Usb1 is an exoribonuclease that shortens the oligo-uridine tail of U6 snRNA, resulting in a terminal 2',3' cyclic phosphate group in most eukaryotes, including humans. Loss of function mutations in human Usb1 cause the rare disorder poikiloderma with neutropenia (PN), and result in U6 snRNAs with elongated 3' ends that are aberrantly adenylated. Here, we show that human Usb1 removes 3' adenosines with 20-fold greater efficiency than uridines, which explains the presence of adenylated U6 snRNAs in cells lacking Usb1. We determined three high-resolution co-crystal structures of Usb1: wild-type Usb1 bound to the substrate analog adenosine 5'-monophosphate, and an inactive mutant bound to RNAs with a 3' terminal adenosine and uridine. These structures, along with QM/MM MD simulations of the catalytic mechanism, illuminate the molecular basis for preferential deadenylation of U6 snRNA. The extent of Usb1 processing is influenced by the secondary structure of U6 snRNA.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

Challenges and advances in the computational modeling of biological phosphate hydrolysis

Phosphate ester hydrolysis is fundamental to many life processes, and has been the topic of substantial experimental and computational research effort. However, even the simplest of phosphate esters can be hydrolyzed through multiple possible pathways that can be difficult to distinguish between, either experimentally, or computationally. Therefore, the mechanisms of both the enzymatic and non-enzymatic reactions have been historically controversial. In the present contribution, we highlight a number of technical issues involved in reliably modeling these computationally challenging reactions, as well as proposing potential solutions. We also showcase examples of our own work in this area, discussing both the non-enzymatic reaction in aqueous solution, as well as insights obtained from the computational modeling of organophosphate hydrolysis and catalytic promiscuity amongst enzymes that catalyze phosphoryl transfer.

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S_N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg²⁺ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

Building a More Predictive Protein Force Field: A Systematic and Reproducible Route to AMBER-FB15

The increasing availability of high-quality experimental data and first-principles calculations creates opportunities for developing more accurate empirical force fields for simulation of proteins. We developed the AMBER-FB15 protein force field by building a high-quality quantum chemical data set consisting of comprehensive potential energy scans and employing the ForceBalance software package for parameter optimization. The optimized potential surface allows for more significant thermodynamic fluctuations away from local minima. In validation studies where simulation results are compared to experimental measurements, AMBER-FB15 in combination with the updated TIP3P-FB water model predicts equilibrium properties with equivalent accuracy, and temperature dependent properties with significantly improved accuracy, in comparison with published models. We also discuss the effect of changing the protein force field and water model on the simulation results.

Perspective: Quantum mechanical methods in biochemistry and biophysics

In this perspective article, I discuss several research topics relevant to quantum mechanical (QM) methods in biophysical and biochemical applications. Due to the immense complexity of biological problems, the key is to develop methods that are able to strike the proper balance of computational efficiency and accuracy for the problem of interest. Therefore, in addition to the development of novel ab initio and density functional theory based QM methods for the study of reactive events that involve complex motifs such as transition metal clusters in metalloenzymes, it is equally important to develop inexpensive QM methods and advanced classical or quantal force fields to describe different physicochemical properties of biomolecules and their behaviors in complex environments. Maintaining a solid connection of these more approximate methods with rigorous QM methods is essential to their transferability and robustness. Comparison to diverse experimental observables helps validate computational models and mechanistic hypotheses as well as driving further development of computational methodologies.

Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin

The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis. Comparison of the observed hydrolysis pathways and the corresponding free energy profiles leads to detailed models for the mechanism of ATP hydrolysis in the pre-powerstroke state and proposes factors that regulate the hydrolysis activity in different conformational states. In the pre-powerstroke state, the scissile P_γ-O_3β bond breaks early in the reaction. Proton transfer from the lytic water to the γ-phosphate through active site residues is an important part of the kinetic bottleneck; several hydrolysis pathways that feature distinct proton transfer routes are found to have similar free energy barriers, suggesting a significant degree of plasticity in the hydrolysis mechanism. Comparison of hydrolysis in the pre-powerstroke state and the closed post-rigor model suggests that optimization of residues beyond the active site for electrostatic stabilization and preorganization is likely important to enzyme design.

Substrate and Transition State Binding in Alkaline Phosphatase Analyzed by Computation of Oxygen Isotope Effects

Enzymes are powerful catalysts, and a thorough understanding of the sources of their catalytic power will facilitate many medical and industrial applications. Here we have studied the catalytic mechanism of alkaline phosphatase (AP), which is one of the most catalytically proficient enzymes known. We have used quantum mechanics calculations and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to model a variety of isotope effects relevant to the reaction of AP. We have calculated equilibrium isotope effects (EIEs), binding isotope effects (BIEs), and kinetic isotope effects (KIEs) for a range of phosphate mono- and diester substrates. The results agree well with experimental values, but the model for the reaction's transition state (TS) differs from the original interpretation of those experiments. Our model indicates that isotope effects on binding make important contributions to measured KIEs on V/K, which complicated interpretation of the measured values. Our results provide a detailed interpretation of the measured isotope effects and make predictions that can test the proposed model. The model indicates that the substrate is deformed in the ground state (GS) of the reaction and partially resembles the TS. The highly preorganized active site preferentially binds conformations that resemble the TS and not the GS, which induces the substrate to adapt to the enzyme, rather than the other way around-as with classic "induced fit" models. The preferential stabilization of the TS over the GS is what lowers the barrier to the chemical step.

The Competing Mechanisms of Phosphate Monoester Dianion Hydrolysis

Despite the numerous experimental and theoretical studies on phosphate monoester hydrolysis, significant questions remain concerning the mechanistic details of these biologically critical reactions. In the present work we construct a linear free energy relationship for phosphate monoester hydrolysis to explore the effect of modulating leaving group pK_a on the competition between solvent- and substrate-assisted pathways for the hydrolysis of these compounds. Through detailed comparative electronic-structure studies of methyl phosphate and a series of substituted aryl phosphate monoesters, we demonstrate that the preferred mechanism is dependent on the nature of the leaving group. For good leaving groups, a strong preference is observed for a more dissociative solvent-assisted pathway. However, the energy difference between the two pathways gradually reduces as the leaving group pK_a increases and creates mechanistic ambiguity for reactions involving relatively poor alkoxy leaving groups. Our calculations show that the transition-state structures vary smoothly across the range of pK_as studied and that the pathways remain discrete mechanistic alternatives. Therefore, while not impossible, a biological catalyst would have to surmount a significantly higher activation barrier to facilitate a substrate-assisted pathway than for the solvent-assisted pathway when phosphate is bonded to good leaving groups. For poor leaving groups, this intrinsic preference disappears.

QM/MM Analysis of Transition States and Transition State Analogues in Metalloenzymes

Enzymology is approaching an era where many problems can benefit from computational studies. While ample challenges remain in quantitatively predicting behavior for many enzyme systems, the insights that often come from computations are an important asset for the enzymology community. Here we provide a primer for enzymologists on the types of calculations that are most useful for mechanistic problems in enzymology. In particular, we emphasize the integration of models that range from small active-site motifs to fully solvated enzyme systems for cross-validation and dissection of specific contributions from the enzyme environment. We then use a case study of the enzyme alkaline phosphatase to illustrate specific application of the methods. The case study involves examination of the binding modes of putative transition state analogues (tungstate and vanadate) to the enzyme. The computations predict covalent binding of these ions to the enzymatic nucleophile and that they adopt the trigonal bipyramidal geometry of the expected transition state. By comparing these structures with transition states found through free energy simulations, we assess the degree to which the transition state analogues mimic the true transition states. Technical issues worth treating with care as well as several remaining challenges to quantitative analysis of metalloenzymes are also highlighted during the discussion.

Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.