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

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.

Improving the Accuracy of Quantum Mechanics/Molecular Mechanics (QM/MM) Models with Polarized Fragment Charges

This paper introduces an economical approach for improving the accuracy and convergence of quantum mechanics/molecular mechanics (QM/MM) models. The approach is tested on a series of neutral and charged amino acids embedded in a 160-water cluster, where their intramolecular proton transfer energies (neutral amino acid → zwitterionic amino acid) were previously obtained at the ωB97X-D/6-31G(d) level of theory. When the charges on the MM atoms were replaced with those obtained at the same QM level of theory used to treat the QM atoms, this significantly improved the accuracy and convergence of the QM/MM models. In particular, the QM/MM model converged to within 1.4 kcal mol^-1 of directly calculated DFT energies for smaller (by as many as 20 waters) QM regions. The use of atomic charges obtained from the natural population analysis yielded the most significant improvement, while other charge schemes such as Mulliken, electrostatic potential, or CM5 led to poorer outcomes. It is further demonstrated that the QM atomic charges can be accurately estimated in a highly efficient manner using an iterative fragmentation approach based on the moving-domain QM/MM method. Similar observations were made when the approach was used to predict the barrier of an S_N2 reaction. Thus, the use of QM-quality atomic charges on MM atoms represents a simple and easy-to-implement strategy for improving the accuracy of QM/MM models.

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.