Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Exploring Product Release from Yeast Cytosine Deaminase with Metadynamics

The yeast cytosine deaminase (yCD) enzyme/5-fluorocytosine prodrug system is a promising candidate for targeted chemotherapeutics. After conversion of the prodrug into the toxic chemotherapeutic drug, 5-fluorouracil (5-FU), the slow product release from the enzyme limits the overall catalytic efficiency of the enzyme/prodrug system. Here, we present a computational study of the product release of the anticancer drug, 5-FU, from yCD using metadynamics. We present a comparison of the 5-FU drug to the natural enzyme product, uracil. We use volume-based metadynamics to compute the free energy landscape for product release and show a modest binding affinity for the product to the enzyme, consistent with experiments. Next, we use infrequent metadynamics to estimate the unbiased release rate from Kramers time-dependent rate theory and find a favorable comparison to experiment with a slower rate of product release for the 5-FU system. Our work demonstrates how adaptive sampling methods can be used to study the protein-ligand unbinding process for engineering enzyme/prodrug systems and gives insights into the molecular mechanism of product release for the yCD/5-FU system.

New Role of Water in Transketolase Catalysis

Transketolase catalyzes the interconversion of keto and aldo sugars. Its coenzyme is thiamine diphosphate. The binding of keto sugar with thiamine diphosphate is possible only after C2 deprotonation of its thiazole ring. It is believed that deprotonation occurs due to the direct transfer of a proton to the amino group of its aminopyrimidine ring. Using mass spectrometry, it is shown that a water molecule is directly involved in the deprotonation process. After the binding of thiamine diphosphate with transketolase and its subsequent cleavage, a thiamine diphosphate molecule is formed with a mass increased by one oxygen molecule. After fragmentation, a thiamine diphosphate molecule is formed with a mass reduced by one and two hydrogen atoms, that is, HO and H₂O are split off. Based on these data, it is assumed that after the formation of holotransketolase, water is covalently bound to thiamine diphosphate, and carbanion is formed as a result of its elimination. This may be a common mechanism for other thiamine enzymes. The participation of a water molecule in the catalysis of the one-substrate transketolase reaction and a possible reason for the effect of the acceptor substrate on the affinity of the donor substrate for active sites are also shown.

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.

Insights into the Dual Shuttle Catalytic Mechanism of Guanine Deaminase

Guanine deaminases (GD) are essential enzymes that help in regulating the nucleobase pool. Since the deamination reaction can result in the accumulation of mutagenic bases that can lead to genomic instability, these enzymes are tightly regulated and are nonpromiscuous. Here, we delineate the basis of their substrate fidelity via entailing the reaction mechanism of deamination by employing density functional theory (DFT) calculations on NE0047, a GD from Nitrosomonas europaea. The results show that, unlike pyrimidine deaminases, which require a single glutamic acid as a proton shuttle, GDs involve two amino acids, E79 and E143 (numbering in NE0047), which control its reactivity. The hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have shown that the first Zn-bound proton transfer to the N3 atom of the substrate is mediated by the E79 residue, and the second proton is transferred to the amine nitrogen of substrate via E143. Moreover, cluster models reveal that the crystallographic water molecules near the active site control the reactivity. A comparison with human GD reveals that the proposed catalytic mechanism is generic, and the knowledge generated here can be effectively applied to design selective inhibitors.

Study of Photochemical Cytosine to Uracil Transition via Ultrafast Photo-Cross-Linking Using Vinylcarbazole Derivatives in Duplex DNA

Gene therapies, including genome editing, RNAi, anti-sense technology and chemical DNA editing are becoming major methods for the treatment of genetic disorders. Techniques like CRISPR-Cas9, zinc finger nuclease (ZFN) and transcription activator-like effector-based nuclease (TALEN) are a few such enzymatic techniques. Most enzymatic genome editing techniques have their disadvantages. Thus, non-enzymatic and non-invasive technologies for nucleic acid editing has been reported in this study which might possess some advantages over the older methods of DNA manipulation. 3-cyanovinyl carbazole (^CNVK) based nucleic acid editing takes advantage of photo-cross-linking between a target pyrimidine and the ^CNVK to afford deamination of cytosine and convert it to uracil. This method previously required the use of high temperatures but, in this study, it has been optimized to take place at physiological conditions. Different counter bases (inosine, guanine and cytosine) complementary to the target cytosine were used, along with derivatives of ^CNVK (^NH2VK and ^OHVK) to afford the deamination at physiological conditions.

Theoretical study of mechanisms for the hydrolytic deamination of cytosine via steered molecular dynamic simulations

Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations. Two pathways were considered to explain the mechanism of this process, where the water molecule attacks the C–N bond to give an intermediate (an amino–hydroxy–oxo structure in the A-path, and a hydroxy–oxo in the B-path) as the determinant step of reaction. Stationary structures along both energy profiles were analyzed at molecular dynamics level, obtaining states with higher free energies than those from electronic calculations in the gas phase and in solution described as a continuous medium. From the results obtained, the more complex A-pathway, with five steps, was kinetically the most favorable (with an endergonic reaction energy of 7.41 kcal mol⁻¹, a high barrier of 67.53 kcal mol⁻¹, and a small velocity constant k₂ = 1.80 × 10⁻³⁷ s⁻¹), concluding that the uracil base can participate in a spontaneous genetic mutation since the uracil–ammonia complex has a long lifetime of 6.10 × 10²⁷ s. This process turns out exergonic and faster when carried out in gas phase simulation or electronic calculation with a continuous medium, due to the disappearance of explicit water molecules that can compete with the assistant molecule.

Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations.

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.

Effect of nucleobase change on cytosine deamination through DNA photo-cross-linking reaction via 3-cyanovinylcarbazole nucleoside

Photo-chemical deamination of cytosine using 3-cyanovinylcarbazole nucleoside (^CNVK) mediated photo-cross-linking is a technique for site-directed mutagenesis. Using this technique in vivo requires the elimination of a high-temperature incubation step; instead, incubation should be carried out under physiological conditions. To improve the reactivity of ^CNVK mediated photo-cross-link induced deamination of cytosine under physiological conditions, an evaluation of base pairing in cytosine was carried out with respect to its deamination. Guanine was replaced with 4 different counter bases (inosine, 2-aminopurine, 5-nitroindole, and nebularine), showing distinct hydrogen bonding patterns with target cytosine, which was incorporated at the -1 position with respect to ^CNVK in the ^CNVK-modified photo-responsive oligodeoxyribonucleotides to ascertain the role of hydrogen bonding in deamination under physiological conditions. Among the counter bases, inosine showed the highest acceleration towards the photo-induced deamination reaction.

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.