Mediation of donor-acceptor distance in an enzymatic methyl transfer reaction

Methyl transfer in psilocybin biosynthesis

Psilocybin, the natural hallucinogen produced by Psilocybe ("magic") mushrooms, holds great promise for the treatment of depression and several other mental health conditions. The final step in the psilocybin biosynthetic pathway, dimethylation of the tryptophan-derived intermediate norbaeocystin, is catalysed by PsiM. Here we present atomic resolution (0.9 Å) crystal structures of PsiM trapped at various stages of its reaction cycle, providing detailed insight into the SAM-dependent methylation mechanism. Structural and phylogenetic analyses suggest that PsiM derives from epitranscriptomic N6-methyladenosine writers of the METTL16 family, which is further supported by the observation that bound substrates physicochemically mimic RNA. Inherent limitations of the ancestral monomethyltransferase scaffold hamper the efficiency of psilocybin assembly and leave PsiM incapable of catalysing trimethylation to aeruginascin. The results of our study will support bioengineering efforts aiming to create novel variants of psilocybin with improved therapeutic properties.

High-throughput computational investigation of protein electrostatics and cavity for SAM-dependent methyltransferases

S-adenosyl methionine (SAM)-dependent methyl transferases (MTases) are a ubiquitous class of enzymes catalyzing dozens of essential life processes. Despite targeting a large space of substrates with diverse intrinsic reactivity, SAM MTases have similar catalytic efficiency. While understanding of MTase mechanism has grown tremendously through the integration of structural characterization, kinetic assays, and multiscale simulations, it remains elusive how these enzymes have evolved to fit the diverse chemical needs of their respective substrates. In this work, we performed a high-throughput molecular modeling analysis of 91 SAM MTases to better understand how their properties (i.e., electric field [EF] strength and active site volumes) help achieve similar catalytic efficiency toward substrates of different reactivity. We found that EF strengths have largely adjusted to make the target atom a better methyl acceptor. For MTases that target RNA/DNA and histone proteins, our results suggest that EF strength accommodates formal hybridization state and variation in cavity volume trends with diversity of substrate classes. Metal ions in SAM MTases contribute negatively to EF strength for methyl donation and enzyme scaffolds tend to offset these contributions.

Fragment-Based Calculations of Enzymatic Thermochemistry Require Dielectric Boundary Conditions

Electronic structure calculations on enzymes require hundreds of atoms to obtain converged results, but fragment-based approximations offer a cost-effective solution. We present calculations on enzyme models containing 500-600 atoms using the many-body expansion, comparing to benchmarks in which the entire enzyme-substrate complex is described at the same level of density functional theory. When the amino acid fragments contain ionic side chains, the many-body expansion oscillates under vacuum boundary conditions but rapid convergence is restored using low-dielectric boundary conditions. This implies that full-system calculations in the gas phase are inappropriate benchmarks for assessing errors in fragment-based approximations. A three-body protocol retains sub-kilocalorie per mole fidelity with respect to a supersystem calculation, as does a two-body calculation combined with a full-system correction at a low-cost level of theory. These protocols pave the way for application of high-level quantum chemistry to large systems via rigorous, ab initio treatment of many-body polarization.

Ectopic expression of MmSERT, a mouse serotonin transporter gene, regulates salt tolerance and ABA sensitivity in apple and Arabidopsis

5-hydroxytryptamine (5-HT) is ubiquitously present in animals and plants, playing a vital regulatory role. SERT, a conserved serotonin reuptake transporter in animals, regulates intracellular and extracellular concentrations of 5-HT. Few studies have reported 5-HT transporters in plants. Hence, we cloned MmSERT, a serotonin reuptake transporter, from Mus musculus. Ectopic expression of MmSERT into apple calli, apple roots and Arabidopsis. Because 5-HT plays a momentous role in plant stress tolerance, we used MmSERT transgenic materials for stress treatment. We found that MmSERT transgenic materials, including apple calli, apple roots and Arabidopsis, exhibited a stronger salt tolerance phenotype. The reactive oxygen species (ROS) produced were significantly lower in MmSERT transgenic materials compared with controls under salt stress. Meanwhile, MmSERT induced the expression of SOS1, SOS3, NHX1, LEA5 and LTP1 in response to salt stress. 5-HT is the precursor of melatonin, which regulates plant growth under adversity and effectively scavenges ROS. Detection of MmSERT transgenic apple calli and Arabidopsis revealed higher melatonin levels than controls. Besides, MmSERT decreased the sensitivity of apple calli and Arabidopsis to abscisic acid (ABA). In summary, these results demonstrated that MmSERT plays a vital role in plant stress resistances, which perhaps serves as a reference for the application of transgenic technology to improve crops in the future.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification

While QM/MM studies of enzymatic reactions are widely used in computational chemistry, the results of such studies are subject to numerous sources of uncertainty, and the effect of different choices by the simulation scientist that are required when setting up QM/MM calculations is often unclear. In particular, the selection of the QM region is crucial for obtaining accurate and reliable results. Simply including amino acids by their distance to the active site is mostly not sufficient as necessary residues are missing or unimportant residues are included without evidence. Here, we take a first step toward quantifying uncertainties in QM/MM calculations by assessing the sensitivity of QM/MM reaction energies with respect to variations of the MM point charges. We show that such a point charge variation analysis (PCVA) can be employed to judge the accuracy of QM/MM reaction energies obtained with a selected QM region and devise a protocol to systematically construct QM regions that minimize this uncertainty. We apply such a PCVA to the example of catechol O-methyltransferase and demonstrate that it provides a simple and reliable approach for the construction of the QM region. Our PCVA-based scheme is computationally efficient and requires only calculations for a system with a minimal QM region. Our work highlights the promise of applying methods of uncertainty quantification in computational chemistry.

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.

Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins

Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.

2'-O methylation of RNA cap in SARS-CoV-2 captured by serial crystallography

The genome of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) coronavirus has a capping modification at the 5'-untranslated region (UTR) to prevent its degradation by host nucleases. These modifications are performed by the Nsp10/14 and Nsp10/16 heterodimers using S-adenosylmethionine as the methyl donor. Nsp10/16 heterodimer is responsible for the methylation at the ribose 2'-O position of the first nucleotide. To investigate the conformational changes of the complex during 2'-O methyltransferase activity, we used a fixed-target serial synchrotron crystallography method at room temperature. We determined crystal structures of Nsp10/16 with substrates and products that revealed the states before and after methylation, occurring within the crystals during the experiments. Here we report the crystal structure of Nsp10/16 in complex with Cap-1 analog (m7GpppAm2'-O). Inhibition of Nsp16 activity may reduce viral proliferation, making this protein an attractive drug target.

Nano-second protein dynamics of key residue at Position 38 in catechol-O-methyltransferase system: a time-resolved fluorescence study

Human catechol-O-methyltransferase, a key enzyme related to neurotransmitter metabolism, catalyses a methyl transfer from S-adenosylmethionine to catechol. Although extensive studies aim to understand the enzyme mechanisms, the connection of protein dynamics and enzyme catalysis is still not clear. Here, W38in (Trp143Phe) and W38in/Y68A (Trp143Phe with Tyr68Ala) mutants were carried out to study the relationship of dynamics and catalysis in nano-second timescale using time-resolved fluorescence lifetimes and Stokes shifts in various solvents. The comprehensive data implied the mutant W38in/Y68A with lower activity is more rigid than the ‘WT’−W38in, suggesting the importance of flexibility at residue 38 to maintain the optimal catalysis.

Transition-State Vibrational Analysis and Isotope Effects for COMT-Catalyzed Methyl Transfer

Isotopic partition-function ratios (IPFRs) computed for transition structures (TSs) of the methyl-transfer reaction catalyzed by catechol O-methyltransferase and modeled by hybrid QM/MM methods are analyzed. The ability of smaller Hessians to reproduce trends in α-3H3 and 14Cα IPFRs as obtained using the much larger subset QM/MM Hessians from which they are extracted is investigated critically. A 6-atom-extracted Hessian reproduces perfectly the α-T3 IPFR values from the full-subset Hessians of all the TSs but not the α-14CIPFRs. Average AM1/OPLS-AA harmonic frequencies and mean-square amplitudes are presented for the 12 normal modes of the α-CH3 moiety within the active site of several enzymic transition structures, together with QM/MM potential energy scans along each of these modes to assess the degree of anharmonicity. A novel investigation of ponderal effects upon IPFRs suggests that the value for α-14C tends toward a limiting minimum whereas that for α-T3 tends toward a limiting maximum as the mass of the rest of the system increases. The transition vector is dominated by motions of atoms within the donor and acceptor moieties and is very well described as a simple combination of Walden-inversion "umbrella" bending and asymmetric stretching of the SCα and CαO bonds. The contribution of atoms of the protein residues Met40, Tyr68, and Asp141 to the transition vector is extremely small. Average valence force constants for the COMT TS show significant differences from early BEBOVIB estimates which were used in support of the compression hypothesis for catalysis. There is no correlation between TS IPFRs and the nonbonded distances for close contacts between the S atom of SAM and Tyr68 or between any of the H atoms of the transferring methyl group and either Met40 or Asp141.

Hydrogen deuterium exchange defines catalytically linked regions of protein flexibility in the catechol O-methyltransferase reaction

Human catechol O-methyltransferase (COMT) has emerged as a model for understanding enzyme-catalyzed methyl transfer from S-adenosylmethionine (AdoMet) to small-molecule catecholate acceptors. Mutation of a single residue (tyrosine 68) behind the methyl-bearing sulfonium of AdoMet was previously shown to impair COMT activity by interfering with methyl donor-acceptor compaction within the activated ground state of the wild type enzyme [J. Zhang, H. J. Kulik, T. J. Martinez, J. P. Klinman, Proc. Natl. Acad. Sci. U.S.A. 112, 7954-7959 (2015)]. This predicts the involvement of spatially defined protein dynamical effects that further tune the donor/acceptor distance and geometry as well as the electrostatics of the reactants. Here, we present a hydrogen/deuterium exchange (HDX)-mass spectrometric study of wild type and mutant COMT, comparing temperature dependences of HDX against corresponding kinetic and cofactor binding parameters. The data show that the impaired Tyr68Ala mutant displays similar breaks in Arrhenius plots of both kinetic and HDX properties that are absent in the wild type enzyme. The spatial resolution of HDX below a break point of 15-20 °C indicates changes in flexibility across ∼40% of the protein structure that is confined primarily to the periphery of the AdoMet binding site. Above 20 °C, Tyr68Ala behaves more like WT in HDX, but its rate and enthalpic barrier remain significantly altered. The impairment of catalysis by Tyr68Ala can be understood in the context of a mutationally induced alteration in protein motions that becomes manifest along and perpendicular to the primary group transfer coordinate.

Moving Through Barriers in Science and Life

This first serious attempt at an autobiographical accounting has forced me to sit still long enough to compile my thoughts about a long personal and scientific journey. I especially hope that my trajectory will be of interest and perhaps beneficial to much younger women who are just getting started in their careers. To paraphrase from Virginia Woolf's writings in A Room of One's Own at the beginning of the 20th century, "for most of history Anonymous was a Woman." However, Ms. Woolf is also quoted as saying "nothing has really happened until it has been described," a harbinger of the enormous historical changes that were about to be enacted and recorded by women in the sciences and other disciplines. The progress in my chosen field of study-the chemical basis of enzyme action-has also been remarkable, from the first description of an enzyme's 3D structure to a growing and deep understanding of the origins of enzyme catalysis.

Examining the Origin of Catalytic Power of Catechol O-Methyltransferase

For decades, there has been debate regarding the origin of the catalytic power of enzymes. In this work, we use the approach of computational chemistry to study the enzyme catechol O-methyltransferase (COMT) and reveal that the two current views on the catalytic mechanism of enzymes, the rate-promoting vibrations and the electric field, may both be viewed as part of the chemical step catalyzed by COMT. However, we show that the rate-promoting vibrations cause the electrostatic effect. This work provides insight into the catalytic mechanism of COMT and resolves a longstanding controversy regarding this enzyme's mechanism.

Quantum Mechanical Description of Electrostatics Provides a Unified Picture of Catalytic Action Across Methyltransferases

Methyl transferases (MTases) are a well-studied class of enzymes for which competing enzymatic enhancement mechanisms have been suggested, ranging from structural methyl group CH···X hydrogen bonds (HBs) to electrostatic- and charge-transfer-driven stabilization of the transition state (TS). We identified all Class I MTases for which reasonable resolution (<2.0 Å) crystal structures could be used to form catalytically competent ternary complexes for multiscale (i.e., quantum-mechanical/molecular-mechanical or QM/MM) simulation of the S_N2 methyl transfer reaction coordinate. The four Class I MTases studied have both distinct functions (e.g., protein repair or biosynthesis) and substrate nucleophiles (i.e., C, N, or O). While CH···X HBs stabilize all reactant complexes, no universal TS stabilization role is found for these interactions in MTases. A consistent picture is instead obtained through analysis of charge transfer and electrostatics, wherein much of cofactor-substrate charge separation is maintained in the TS region, and electrostatic potential is correlated with substrate nucleophilicity (i.e., intrinsic reactivity).

Equatorial Active Site Compaction and Electrostatic Reorganization in Catechol-O-methyltransferase

Catechol-O-methyltransferase (COMT) is a model S-adenosyl-l-methionine (SAM) dependent methyl transferase, which catalyzes the methylation of catecholamine neurotransmitters such as dopamine in the primary pathway of neurotransmitter deactivation in animals. Despite extensive study, there is no consensus view of the physical basis of catalysis in COMT. Further progress requires experimental data that directly probes active site geometry, protein dynamics and electrostatics, ideally in a range of positions along the reaction coordinate. Here we establish that sinefungin, a fungal-derived inhibitor of SAM-dependent enzymes that possess transition state-like charge on the transferring group, can be used as a transition state analog of COMT when combined with a catechol. X-ray crystal structures and NMR backbone assignments of the ternary complexes of the soluble form of human COMT containing dinitrocatechol, Mg2+ and SAM or sinefungin were determined. Comparison and further analysis with the aid of density functional theory calculations and molecular dynamics simulations provides evidence for active site "compaction", which is driven by electrostatic stabilization between the transferring methyl group and "equatorial" active site residues that are orthogonal to the donor-acceptor (pseudo reaction) coordinate. We propose that upon catecholamine binding and subsequent proton transfer to Lys 144, the enzyme becomes geometrically preorganized, with little further movement along the donor-acceptor coordinate required for methyl transfer. Catalysis is then largely facilitated through stabilization of the developing charge on the transferring methyl group via "equatorial" H-bonding and electrostatic interactions orthogonal to the donor-acceptor coordinate.

Revealing quantum mechanical effects in enzyme catalysis with large-scale electronic structure simulation

Enzymes have evolved to facilitate challenging reactions at ambient conditions with specificity seldom matched by other catalysts. Computational modeling provides valuable insight into catalytic mechanism, and the large size of enzymes mandates multi-scale, quantum mechanical-molecular mechanical (QM/MM) simulations. Although QM/MM plays an essential role in balancing simulation cost to enable sampling with full QM treatment needed to understand electronic structure in enzyme active sites, the relative importance of these two strategies for understanding enzyme mechanism is not well known. We explore challenges in QM/MM for studying the reactivity and stability of three diverse enzymes: i) Mg2+-dependent catechol O-methyltransferase (COMT), ii) radical enzyme choline trimethylamine lyase (CutC), and iii) DNA methyltransferase (DNMT1), which has structural Zn2+ binding sites. In COMT, strong non-covalent interactions lead to long range coupling of electronic structure properties across the active site, but the more isolated nature of the metallocofactor in DNMT1 leads to faster convergence of some properties. We quantify these effects in COMT by computing covariance matrices of by-residue electronic structure properties during dynamics and along the reaction coordinate. In CutC, we observe spontaneous bond cleavage following initiation events, highlighting the importance of sampling and dynamics. We use electronic structure analysis to quantify the relative importance of CHO and OHO non-covalent interactions in imparting reactivity. These three diverse cases enable us to provide some general recommendations regarding QM/MM simulation of enzymes.

Crystallographic and Computational Characterization of Methyl Tetrel Bonding in S-Adenosylmethionine-Dependent Methyltransferases

Tetrel bonds represent a category of non-bonding interaction wherein an electronegative atom donates a lone pair of electrons into the sigma antibonding orbital of an atom in the carbon group of the periodic table. Prior computational studies have implicated tetrel bonding in the stabilization of a preliminary state that precedes the transition state in S_N2 reactions, including methyl transfer. Notably, the angles between the tetrel bond donor and acceptor atoms coincide with the prerequisite geometry for the S_N2 reaction. Prompted by these findings, we surveyed crystal structures of methyltransferases in the Protein Data Bank and discovered multiple instances of carbon tetrel bonding between the methyl group of the substrate S-adenosylmethionine (AdoMet) and electronegative atoms of small molecule inhibitors, ions, and solvent molecules. The majority of these interactions involve oxygen atoms as the Lewis base, with the exception of one structure in which a chlorine atom of an inhibitor functions as the electron donor. Quantum mechanical analyses of a representative subset of the methyltransferase structures from the survey revealed that the calculated interaction energies and spectral properties are consistent with the values for bona fide carbon tetrel bonds. The discovery of methyl tetrel bonding offers new insights into the mechanism underlying the S_N2 reaction catalyzed by AdoMet-dependent methyltransferases. These findings highlight the potential of exploiting these interactions in developing new methyltransferase inhibitors.

Large-scale QM/MM free energy simulations of enzyme catalysis reveal the influence of charge transfer

Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations provide key insights into enzyme structure-function relationships. Numerous studies have demonstrated that large QM regions are needed to systematically converge ground state, zero temperature properties with electrostatic embedding QM/MM. However, it is not well known if ab initio QM/MM free energy simulations have this same dependence, in part due to the hundreds of thousands of energy evaluations required for free energy estimations that in turn limit QM region size. Here, we leverage recent advances in electronic structure efficiency and accuracy to carry out range-separated hybrid density functional theory free energy simulations in a representative methyltransferase. By studying 200 ps of ab initio QM/MM dynamics for each of five QM regions from minimal (64 atoms) to one-sixth of the protein (544 atoms), we identify critical differences between large and small QM region QM/MM in charge transfer between substrates and active site residues as well as in geometric structure and dynamics that coincide with differences in predicted free energy barriers. Distinct geometric and electronic structure features in the largest QM region indicate that important aspects of enzymatic rate enhancement in methyltransferases are identified with large-scale electronic structure.

Catechol-O-Methyltransferase and UDP-Glucuronosyltransferases in the Metabolism of Baicalein in Different Species

BACKGROUND

Baicalein is the major bioactive flavonoid in some herb medicines and dietary plants; however, the detailed metabolism pathway of its major metabolite oroxylin A-7-O-β-D-glucuronide in human was not clear. It was important to illustrate the major metabolic enzymes that participate in its elimination for the clinic use of baicalein.

OBJECTIVES

We first revealed a two-step metabolism profile for baicalein and illustrated the combination of catechol-O-methyltransferase (COMT) and uridine diphosphate-glucuronosyltransferases (UGTs) in drug metabolism, further evaluated its bioactivity variation during drug metabolism.

METHODS

The metabolism profiles were systematically characterized in different human biology preparations; after then, the anti-inflammatory activities of metabolites were evaluated in LPS-induced RAW264.7 cell.

RESULTS

The first-step metabolite of baicalein was isolated and identified as oroxylin A; soluble-bound COMT (S-COMT) was the major enzyme responsible for its biotransformation. Specially, position 108 mutation of S-COMT significantly decreases the elimination. Meantime, oroxylin A was rapidly metabolized by UGTs, UGT1A1, -1A3, -1A6, -1A7, -1A8, -1A9, and -1A10 which were involved in the glucuronidation. Considerable species differences were observed with 1060-fold K _m (3.05 ± 1.86-3234 ± 475 μM) and 330-fold CL_int (5.93-1973 μL/min/mg) variations for baicalein metabolism. Finally, the middle metabolite oroxylin A exhibited a potent anti-inflammatory activity with the IC₅₀ value of 28 μM.

CONCLUSION

The detailed kinetic parameters indicated that COMT provide convenience for the next glucuronidation; monkey would be a preferred animal model for the preclinical investigation of baicalein. Importantly, oroxylin A should be reconsidered in evaluating baicalein efficacy against inflammatory diseases.

Insights on the Origin of Catalysis on Glycine N-Methyltransferase from Computational Modeling

The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor-acceptor distance (DAD): Is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The computed activation free energies and KIEs are in good agreement with experimental data, but the mutations do not meaningfully affect the DAD: Compression cannot explain the experimental variations on KIEs. On the contrary, electrostatic properties in the active site correlate with the catalytic activity of wild type and mutants. The plasticity of the enzyme moderates the effects of the mutations, explaining the rather small degree of variation in KIEs and reactivities.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

Origins of Enzyme Catalysis: Experimental Findings for C-H Activation, New Models, and Their Relevance to Prevailing Theoretical Constructs

The physical basis for enzymatic rate accelerations is a subject of great fundamental interest and of direct relevance to areas that include the de novo design of green catalysts and the pursuit of new drug regimens. Extensive investigations of C-H activating systems have provided considerable insight into the relationship between an enzyme's overall structure and the catalytic chemistry at its active site. This Perspective highlights recent experimental data for two members of distinct, yet iconic C-H activation enzyme classes, lipoxygenases and prokaryotic alcohol dehydrogenases. The data necessitate a reformulation of the dominant textbook definition of biological catalysis. A multidimensional model emerges that incorporates a range of protein motions that can be parsed into a combination of global stochastic conformational thermal fluctuations and local donor-acceptor distance sampling. These motions are needed to achieve a high degree of precision with regard to internuclear distances, geometries, and charges within the active site. The available model also suggests a physical framework for understanding the empirical enthalpic barrier in enzyme-catalyzed processes. We conclude by addressing the often conflicting interface between computational and experimental chemists, emphasizing the need for computation to predict experimental results in advance of their measurement.

Integration of Tmc1/2 into the mechanotransduction complex in zebrafish hair cells is regulated by Transmembrane O-methyltransferase (Tomt)

Transmembrane O-methyltransferase (TOMT/LRTOMT) is responsible for non-syndromic deafness DFNB63. However, the specific defects that lead to hearing loss have not been described. Using a zebrafish model of DFNB63, we show that the auditory and vestibular phenotypes are due to a lack of mechanotransduction (MET) in Tomt-deficient hair cells. GFP-tagged Tomt is enriched in the Golgi of hair cells, suggesting that Tomt might regulate the trafficking of other MET components to the hair bundle. We found that Tmc1/2 proteins are specifically excluded from the hair bundle in tomt mutants, whereas other MET complex proteins can still localize to the bundle. Furthermore, mouse TOMT and TMC1 can directly interact in HEK 293 cells, and this interaction is modulated by His183 in TOMT. Thus, we propose a model of MET complex assembly where Tomt and the Tmcs interact within the secretory pathway to traffic Tmc proteins to the hair bundle.

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

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

Kinetic Isotope Effects and Transition State Structure for Human Phenylethanolamine N-Methyltransferase

Phenylethanolamine N-methyltransferase (PNMT) catalyzes the S-adenosyl-l-methionine (SAM)-dependent conversion of norepinephrine to epinephrine. Epinephrine has been associated with critical processes in humans including the control of respiration and blood pressure. Additionally, PNMT activity has been suggested to play a role in hypertension and Alzheimer's disease. In the current study, labeled SAM substrates were used to measure primary methyl-14C and 36S and secondary methyl-3H, 5'-3H, and 5'-14C intrinsic kinetic isotope effects for human PNMT. The transition state of human PNMT was modeled by matching kinetic isotope effects predicted via quantum chemical calculations to intrinsic values. The model provides information on the geometry and electrostatics of the human PNMT transition state structure and indicates that human PNMT catalyzes the formation of epinephrine through an early SN2 transition state in which methyl transfer is rate-limiting.

13C ENDOR Spectroscopy of Lipoxygenase-Substrate Complexes Reveals the Structural Basis for C-H Activation by Tunneling

In enzymatic C-H activation by hydrogen tunneling, reduced barrier width is important for efficient hydrogen wave function overlap during catalysis. For native enzymes displaying nonadiabatic tunneling, the dominant reactive hydrogen donor-acceptor distance (DAD) is typically ca. 2.7 Å, considerably shorter than normal van der Waals distances. Without a ground state substrate-bound structure for the prototypical nonadiabatic tunneling system, soybean lipoxygenase (SLO), it has remained unclear whether the requisite close tunneling distance occurs through an unusual ground state active site arrangement or by thermally sampling conformational substates. Herein, we introduce Mn2+ as a spin-probe surrogate for the SLO Fe ion; X-ray diffraction shows Mn-SLO is structurally faithful to the native enzyme. 13C ENDOR then reveals the locations of 13C10 and reactive 13C11 of linoleic acid relative to the metal; 1H ENDOR and molecular dynamics simulations of the fully solvated SLO model using ENDOR-derived restraints give additional metrical information. The resulting three-dimensional representation of the SLO active site ground state contains a reactive (a) conformer with hydrogen DAD of ∼3.1 Å, approximately van der Waals contact, plus an inactive (b) conformer with even longer DAD, establishing that stochastic conformational sampling is required to achieve reactive tunneling geometries. Tunneling-impaired SLO variants show increased DADs and variations in substrate positioning and rigidity, confirming previous kinetic and theoretical predictions of such behavior. Overall, this investigation highlights the (i) predictive power of nonadiabatic quantum treatments of proton-coupled electron transfer in SLO and (ii) sensitivity of ENDOR probes to test, detect, and corroborate kinetically predicted trends in active site reactivity and to reveal unexpected features of active site architecture.

Drug Metabolism in the Liver

Metabolism is a biotransformation process, where endogenous and exogenous compounds are converted to more polar products to facilitate their elimination from the body. The process of metabolism is divided into 3 phases. Phase I metabolism involves functionalization reactions. Phase II drug metabolism is a conjugation reaction. Phase III refers to transporter-mediated elimination of drug and/or metabolites from body normally via liver, gut, kidney, or lung. This review presents basic information on drug-metabolizing enzymes and potential factors that might affect the metabolic capacities of the enzyme or alter drug response or drug-mediated toxicities.

How Large Should the QM Region Be in QM/MM Calculations? The Case of Catechol O-Methyltransferase

Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a prototypical methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500-600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (approximately 11-16 residues or 200-300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes.

Structural and functional aspects of decorsin and its analog as recognized by integrin αIIbβ3

Decorsin is an antagonist of platelet glycoprotein integrin αIIbβ3 on platelets; the protein is 39 amino acids long with three disulfide bridges in its tertiary structure. To demonstrate decorsin's mechanism of action, we applied the computational virtual technique and platelet aggregation inhibition assay, which showed that the flanking amino-acid residues of the Arg-Gly-Asp (RGD) motif play an important role in platelet aggregation. The computational simulations revealed that the RGD motif mainly contributes to the stability of the complex when decorsion interacts with integrin αIIbβ3. However, the C-terminal residues, such as 34A→W and 35D→R, was also found to possibly play a key role in their binding structures. Moreover, we produced a decorsin analog (A34W plus D35R decorsin), in which the 34A (alanine) and 35D (aspartic acid) residues were respectively substituted by W (tryptophan) and R (arginine). This isoform was then recombinantly expressed in Escherichia coli. Intriguingly, this mutant type showed higher anti-platelet aggregation activity than the wildtype. Our study may further contribute to finding decorsin mutants with higher anti-platelet aggregation activity.

Design of Potent and Druglike Nonphenolic Inhibitors for Catechol O-Methyltransferase Derived from a Fragment Screening Approach Targeting the S-Adenosyl-l-methionine Pocket

A fragment screening approach designed to target specifically the S-adenosyl-l-methionine pocket of catechol O-methyl transferase allowed the identification of structurally related fragments of high ligand efficiency and with activity on the described orthogonal assays. By use of a reliable enzymatic assay together with X-ray crystallography as guidance, a series of fragment modifications revealed an SAR and, after several expansions, potent lead compounds could be obtained. For the first time nonphenolic and small low nanomolar potent, SAM competitive COMT inhibitors are reported. These compounds represent a novel series of potent COMT inhibitors that might be further optimized to new drugs useful for the treatment of Parkinson's disease, as adjuncts in levodopa based therapy, or for the treatment of schizophrenia.

Functional AdoMet Isosteres Resistant to Classical AdoMet Degradation Pathways

S-adenosyl-l-methionine (AdoMet) is an essential enzyme cosubstrate in fundamental biology with an expanding range of biocatalytic and therapeutic applications. We report the design, synthesis, and evaluation of stable, functional AdoMet isosteres that are resistant to the primary contributors to AdoMet degradation (depurination, intramolecular cyclization, and sulfonium epimerization). Corresponding biochemical and structural studies demonstrate the AdoMet surrogates to serve as competent enzyme cosubstrates and to bind a prototypical class I model methyltransferase (DnrK) in a manner nearly identical to AdoMet. Given this conservation in function and molecular recognition, the isosteres presented are anticipated to serve as useful surrogates in other AdoMet-dependent processes and may also be resistant to, and/or potentially even inhibit, other therapeutically relevant AdoMet-dependent metabolic transformations (such as the validated drug target AdoMet decarboxylase). This work also highlights the ability of the prototypical class I model methyltransferase DnrK to accept non-native surrogate acceptors as an enabling feature of a new high-throughput methyltransferase assay.

Exploring the Dependence of QM/MM Calculations of Enzyme Catalysis on the Size of the QM Region

Although QM/MM calculations are the primary current tool for modeling enzymatic reactions, the reliability of such calculations can be limited by the size of the QM region. Thus, we examine in this work the dependence of QM/MM calculations on the size of the QM region, using the reaction of catechol-O-methyl transferase (COMT) as a test case. Our study focuses on the effect of adding residues to the QM region on the activation free energy, obtained with extensive QM/MM sampling. It is found that the sensitivity of the activation barrier to the size of the QM is rather limited, while the dependence of the reaction free energy is somewhat larger. Of course, the results depend on the inclusion of the first solvation shell in the QM regions. For example, the inclusion of the Mg²⁺ ion can change the activation barrier due to charge transfer effects. However, such effects can easily be included in semiempirical approaches by proper parametrization. Overall, we establish that QM/MM calculations of activation barriers of enzymatic reactions are not highly sensitive to the size of the QM region, beyond the immediate region that describes the reacting atoms.

Computational Investigation of the Interplay of Substrate Positioning and Reactivity in Catechol O-Methyltransferase

Catechol O-methyltransferase (COMT) is a SAM- and Mg2+-dependent methyltransferase that regulates neurotransmitters through methylation. Simulations and experiments have identified divergent catecholamine substrate orientations in the COMT active site: molecular dynamics simulations have favored a monodentate coordination of catecholate substrates to the active site Mg2+, and crystal structures instead preserve bidentate coordination along with short (2.65 Å) methyl donor-acceptor distances. We carry out longer dynamics (up to 350 ns) to quantify interconversion between bidentate and monodentate binding poses. We provide a systematic determination of the relative free energy of the monodentate and bidentate structures in order to identify whether structural differences alter the nature of the methyl transfer mechanism and source of enzymatic rate enhancement. We demonstrate that the bidentate and monodentate binding modes are close in energy but separated by a 7 kcal/mol free energy barrier. Analysis of interactions in the two binding modes reveals that the driving force for monodentate catecholate orientations in classical molecular dynamics simulations is derived from stronger electrostatic stabilization afforded by alternate Mg2+ coordination with strongly charged active site carboxylates. Mixed semi-empirical-classical (SQM/MM) substrate C-O distances (2.7 Å) for the bidentate case are in excellent agreement with COMT X-ray crystal structures, as long as charge transfer between the substrates, Mg2+, and surrounding ligands is permitted. SQM/MM free energy barriers for methyl transfer from bidentate and monodentate catecholate configurations are comparable at around 21-22 kcal/mol, in good agreement with experiment (18-19 kcal/mol). Overall, the work suggests that both binding poses are viable for methyl transfer, and accurate descriptions of charge transfer and electrostatics are needed to provide balanced relative barriers when multiple binding poses are accessible, for example in other transferases.

Convergent Mechanistic Features between the Structurally Diverse N- and O-Methyltransferases: Glycine N-Methyltransferase and Catechol O-Methyltransferase

Although an enormous and still growing number of biologically diverse methyltransferases have been reported and identified, a comprehensive understanding of the enzymatic methyl transfer mechanism is still lacking. Glycine N-methyltransferase (GNMT), a member of the family that acts on small metabolites as the substrate, catalyzes methyl transfer from S-adenosyl-l-methionine (AdoMet) to glycine to form S-adenosyl-l-homocysteine and sarcosine. We report primary carbon (¹²C/¹⁴C) and secondary (¹H₃/³H₃) kinetic isotope effects at the transferred methyl group, together with ¹H₃/³H₃ binding isotope effects for wild-type GNMT and a series of Tyr21 mutants. The data implicate a compaction effect in the methyl transfer step that is conferred by the protein structure. Furthermore, a remarkable similarity of properties is observed between GNMT and catechol O-methyltransferase, despite significant differences between these enzymes with regard to their active site structures and catalyzed reactions. We attribute these results to a catalytically relevant reduction in the methyl donor–acceptor distance that is dependent on a tyrosine side chain positioned behind the methyl-bearing sulfur of AdoMet.

Enhancing Paradynamics for QM/MM Sampling of Enzymatic Reactions

Despite the enormous increase in computer power, it is still extremely challenging to obtain computationally converging sampling of ab initio QM/MM (QM(ai)/MM) free energy surfaces in condensed phases. The sampling problem can be significantly reduced by the use of the reference potential paradynamics (PD) approach, but even this approach still requires major computer time in studies of enzymatic reactions. To further reduce the sampling problem we developed here a new PD version where we use an empirical valence bond reference potential that has a minimum rather than a maximum at the transition state region of the target potential (this is accomplished conveniently by shifting the EVB of the product state). Hence, we can map the TS region in a more efficient way. Here, we introduce and validate the inverted EVB PD approach. The validation involves the study of the S(N)2 step of the reaction catalyzed by haloakene dehalogenase (DhlA) and the GTP hydrolysis in the RasGAP system. In addition, we have also studied the corresponding reaction in water for each of the systems described here and the reaction involving trimethylsulfonium and dimethylamine in solution. The results are encouraging and the new strategy appears to provide a powerful way of evaluating QM(ai)/MM activation free energies.

Influence of Equatorial CH⋅⋅⋅O Interactions on Secondary Kinetic Isotope Effects for Methyl Transfer

DFT calculations for methyl cation complexed within a constrained cage of water molecules permit the controlled manipulation of the "axial" donor/acceptor distance and the "equatorial" distance to hydrogen-bond acceptors. The kinetic isotope effect k(CH3)/k(CT3) for methyl transfer within a cage with a short axial distance becomes less inverse for shorter equatorial C⋅⋅⋅O distances: a decrease of 0.5 Å results in a 3 % increase at 298 K. Kinetic isotope effects in AdoMet-dependent methyltransferases may be m∧odulated by CH⋅⋅⋅O hydrogen bonding, and factors other than axial compression may contribute, at least partially, to recently reported isotope-effect variations for catechol-O-methyltransferase and its mutant structures.

Transition state for the NSD2-catalyzed methylation of histone H3 lysine 36

Nuclear receptor SET domain containing protein 2 (NSD2) catalyzes the methylation of histone H3 lysine 36 (H3K36). It is a determinant in Wolf-Hirschhorn syndrome and is overexpressed in human multiple myeloma. Despite the relevance of NSD2 to cancer, there are no potent, selective inhibitors of this enzyme reported. Here, a combination of kinetic isotope effect measurements and quantum chemical modeling was used to provide subangstrom details of the transition state structure for NSD2 enzymatic activity. Kinetic isotope effects were measured for the methylation of isolated HeLa cell nucleosomes by NSD2. NSD2 preferentially catalyzes the dimethylation of H3K36 along with a reduced preference for H3K36 monomethylation. Primary Me-(14)C and (36)S and secondary Me-(3)H3, Me-(2)H3, 5'-(14)C, and 5'-(3)H2 kinetic isotope effects were measured for the methylation of H3K36 using specifically labeled S-adenosyl-l-methionine. The intrinsic kinetic isotope effects were used as boundary constraints for quantum mechanical calculations for the NSD2 transition state. The experimental and calculated kinetic isotope effects are consistent with an SN2 chemical mechanism with methyl transfer as the first irreversible chemical step in the reaction mechanism. The transition state is a late, asymmetric nucleophilic displacement with bond separation from the leaving group at (2.53 Å) and bond making to the attacking nucleophile (2.10 Å) advanced at the transition state. The transition state structure can be represented in a molecular electrostatic potential map to guide the design of inhibitors that mimic the transition state geometry and charge.