Stabilization of a Transition-State Analogue at the Active Site of Yeast Cytosine Deaminase: Importance of Proton Transfers

Combined Theoretical, Bioinformatic, and Biochemical Analyses of RNA Editing by Adenine Base Editors

Adenine base editors (ABEs) have been subjected to multiple rounds of mutagenesis with the goal of optimizing their function as efficient and precise genome editing agents. Despite an ever-expanding data set of ABE mutants and their corresponding DNA or RNA-editing activity, the molecular mechanisms defining these changes remain to be elucidated. In this study, we provide a systematic interpretation of the nature of these mutations using an entropy-based classification model that relies on evolutionary data from extant protein sequences. Using this model in conjunction with experimental analyses, we identify two previously reported mutations that form an epistatic pair in the RNA-editing functional landscape of ABEs. Molecular dynamics simulations reveal the atomistic details of how these two mutations affect substrate-binding and catalytic activity, via both individual and cooperative effects, hence providing insights into the mechanisms through which these two mutations are epistatically coupled.

Product release mechanism and the complete enzyme catalysis cycle in yeast cytosine deaminase (yCD): A computational study

Yeast cytosine deaminase (yCD) is critical in gene-directed enzyme prodrug therapy as it catalyzes the hydrolytic deamination of cytosine. The product (uracil) release process is considered as rate-limiting in the whole enzymatic catalysis and includes the cleavage of the uracil-metal bond and the delivery of free uracil out of the reactive site. Herein extensive combined random acceleration molecular dynamics (RAMD) and molecular dynamics (MD) simulations coupled with the umbrella sampling technique have been performed to study the product transport mechanism. Five channels have been identified, and the thermodynamic and dynamic characterizations for the two most favorable channels have been determined and analyzed. The free energy barrier for the most beneficial pathway is about 13kcal/mol and mainly results from the cleavage of hydrogen bonds between the ligand uracil and surrounding residues Asn51, Glu64, and Asp155. The conjugated rings of Phe114 and Trp152 play gating and guiding roles in the product delivery via π⋯π van der Waals interactions with the product. Finally, the full cycle of the enzymatic catalysis has been determined, making the whole process computationally more precise.

Mechanisms for the deamination reaction of 8-oxoguanine catalyzed by 8-oxoguanine deaminase: A combined QM/MM molecular dynamics study

The deamination reaction of 8-oxoguanine (8-oxoG) catalyzed by 8-oxoguanine deaminase (8-oxoGD) plays a critically important role in the DNA repair activity for oxidative damage. In order to elucidate the complete enzymatic catalysis mechanism at the stages of 8-oxoguanine binding, departure of 2-hydroxy-1H-purine-6,8(7H,9H)-dione from the active site, and formation of 8-oxoxanthine, extensive combined QM(PM3)/MM molecular dynamics simulations have been performed. Computations show that the rate-limiting step corresponds to the nucleophilic attack from zinc-coordinate hydroxide group to free 8-oxoguanine. Through conformational analyses, we demonstrate that Trp115, Trp123 and Leu119 connect to O8@8-oxoguanine with hydrogen bonds, and we suggest that mutations of tryptophan (115 and 123) to histidine or phenylalanine and mutation of leucine (119) to alanine could potentially lead to a mutant with enhanced activity. On this ground, a proton transfer mechanism for the formation of 8-oxoxanthine was further discussed. Both Glu218 and water molecule could be used as proton shuttles, and water molecule plays a major role in proton transfer in substrate. On the other hand, comparative simulations on the deamination of guanine and isocytosine reveal that, for the helping of hydrogen bonds between O8@8-oxoguanine and enzyme, O8@8-oxoguanine is the fastest to be deaminated among the three substrates which are also supported by the experimental kinetic constants.

QM/MM MD and Free Energy Simulation Study of Methyl Transfer Processes Catalyzed by PKMTs and PRMTs

Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here, we review the results of quantum mechanics/molecular mechanics molecular dynamics and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.

QM/MM MD and free energy simulation study of methyl transfer processes catalyzed by PKMTs and PRMTs

Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here we review the results of quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.

WHICH IS THE PROTON-SHUTTLE IN ISOXANTHOPTERIN DEAMINASE? QM/MM MD UNDERSTANDING

The deamination process of isoxanthopterin catalyzed by isoxanthopterin deaminase was determined using the combined QM(PM3)/MM molecular dynamics simulations. In this paper, the updated PM3 parameters were employed for zinc ions and the initial model was built up based on the crystal structure. Proton transfer and following steps have been investigated in two paths: Asp336 and His285 serve as the proton shuttle, respectively. Our simulations showed that His285 is more effective than Aap336 in proton transfer for deamination of isoxanthopterin. As hydrogen bonds between the substrate and surrounding residues play a key role in nucleophilic attack, we suggested mutating Thr195 to glutamic acid, which could enhance the hydrogen bonds and help isoxanthopterin get close to the active site. The simulations which change the substrate to pterin 6-carboxylate also performed for comparison. Our results provide reference for understanding of the mechanism of deaminase and for enhancing the deamination rate of isoxanthopterin deaminase.

Role of glutamate 64 in the activation of the prodrug 5-fluorocytosine by yeast cytosine deaminase

Yeast cytosine deaminase (yCD) catalyzes the hydrolytic deamination of cytosine to uracil as well as the deamination of the prodrug 5-fluorocytosine (5FC) to the anticancer drug 5-fluorouracil. In this study, the role of Glu64 in the activation of the prodrug 5FC was investigated by site-directed mutagenesis, biochemical, nuclear magnetic resonance (NMR), and computational studies. Steady-state kinetics studies showed that the mutation of Glu64 causes a dramatic decrease in k(cat) and a dramatic increase in K(m), indicating Glu64 is important for both binding and catalysis in the activation of 5FC. (19)F NMR experiments showed that binding of the inhibitor 5-fluoro-1H-pyrimidin-2-one (5FPy) to the wild-type yCD causes an upfield shift, indicating that the bound inhibitor is in the hydrated form, mimicking the transition state or the tetrahedral intermediate in the activation of 5FC. However, binding of 5FPy to the E64A mutant enzyme causes a downfield shift, indicating that the bound 5FPy remains in an unhydrated form in the complex with the mutant enzyme. (1)H and (15)N NMR analysis revealed trans-hydrogen bond D/H isotope effects on the hydrogen of the amide of Glu64, indicating that the carboxylate of Glu64 forms two hydrogen bonds with the hydrated 5FPy. ONIOM calculations showed that the wild-type yCD complex with the hydrated form of the inhibitor 1H-pyrimidin-2-one is more stable than the initial binding complex, and in contrast, with the E64A mutant enzyme, the hydrated inhibitor is no longer favored and the conversion has a higher activation energy, as well. The hydrated inhibitor is stabilized in the wild-type yCD by two hydrogen bonds between it and the carboxylate of Glu64 as revealed by (1)H and (15)N NMR analysis. To explore the functional role of Glu64 in catalysis, we investigated the deamination of cytosine catalyzed by the E64A mutant by ONIOM calculations. The results showed that without the assistance of Glu64, both proton transfers before and after the formation of the tetrahedral reaction intermediate become partially rate-limiting steps. The results of the experimental and computational studies together indicate that Glu64 plays a critical role in both the binding and the chemical transformation in the conversion of the prodrug 5FC to the anticancer drug 5-fluorouracil.

QM/MM free energy simulations of salicylic acid methyltransferase: effects of stabilization of TS-like structures on substrate specificity

Salicylic acid methyltransferases (SAMTs) synthesize methyl salicylate (MeSA) using salicylate as the substrate. MeSA synthesized in plants may function as an airborne signal to activate the expression of defense-related genes and could also be a critical mobile signaling molecule that travels from the site of plant infection to establish systemic immunity in the induction of disease resistance. Here the results of QM/MM free energy simulations for the methyl transfer process in Clarkia breweri SAMT (CbSAMT) are reported to determine the origin of the substrate specificity of SAMTs. The free energy barrier for the methyl transfer from S-adenosyl-L-methionine (AdoMet) to 4-hydroxybenzoate in CbSAMT is found to be about 5 kcal/mol higher than that from AdoMet to salicylate, consistent with the experimental observations. It is suggested that the relatively high efficiency for the methylation of salicylate compared to 4-hydroxybenzoate is due, at least in part, to the reason that a part of the stabilization of the transition state (TS) configuration is already reflected in the reactant complex, presumably, through the binding. The results seem to indicate that the creation of the substrate complex (e.g., through mutagenesis and substrate modifications) with its structure closely resembling TS might be fruitful for improving the catalytic efficiency for some enzymes. The results show that the computer simulations may provide important insights into the origin of the substrate specificity for the SABATH family and could be used to help experimental efforts in generating engineered enzymes with altered substrate specificity.

QM/MM X-ray refinement of zinc metalloenzymes

Zinc metalloenzymes play an important role in biology. However, due to the limitation of molecular force field energy restraints used in X-ray refinement at medium or low resolutions, the precise geometry of the zinc coordination environment can be difficult to distinguish from ambiguous electron density maps. Due to the difficulties involved in defining accurate force fields for metal ions, the QM/MM (quantum-mechanical/molecular-mechanical) method provides an attractive and more general alternative for the study and refinement of metalloprotein active sites. Herein we present three examples that indicate that QM/MM based refinement yields a superior description of the crystal structure based on R and R(free) values and on the inspection of the zinc coordination environment. It is concluded that QM/MM refinement is an useful general tool for the improvement of the metal coordination sphere in metalloenzyme active sites.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.