Understanding the R882H mutation effects of DNA methyltransferase DNMT3A: a combination of molecular dynamics simulations and QM/MM calculations

Insights into the Inhibitory Mechanisms of the Covalent Drugs for DNMT3A

The perturbations of DNA methyltransferase 3 alpha (DNMT3A) may cause uncontrolled gene expression, resulting in cancers and tumors. The DNMT inhibitors Azacytidine (AZA) and Zebularine (ZEB) inhibit the DNMT family with no specificities, and consequently would bring side effects during the treatment. Therefore, it is vital to understand the inhibitory mechanisms in DNMT3A to inform the new inhibitor design for DNMTs. Herein, we carried out molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) simulations to investigate the inhibitory mechanisms of the AZA and ZEB. The results were compared to the methyl transfer of cytosine. We showed how the AZA might stop the methyl transfer process, whereas the ZEB might be stuck in a methyl-transferred intermediate (IM3). The IM3 state then fails the elimination due to the unique protein dynamics that result in missing the catalytic water chain. Our results brought atomic-level insights into the mechanisms of the two drugs in DNMT3A, which could benefit the new generation of drug design for the DNMTs.

Effectiveness of ruminal xylanase with an extra proline-rich C-terminus on lignocellulosic biomass degradation

The efficient degradation of plant polysaccharides in agricultural waste requires xylanases with high catalytic activity. In this study, the C-terminal proline-rich GH10 xylanase XynA from sheep rumen was investigated using product analysis, structural characterization, truncated and site-directed mutagenesis, molecular dynamics simulation, and application evaluation, revealing that the proline-rich C-terminus contributes to the interaction at the substrate-binding pocket to reduce the binding free energy. Compared to the C-terminally truncated enzyme XynA-Tr, XynA has a more favorable conformation for proton transfer and affinity attack, facilitating the degradation of oligomeric and beechwood xylan without altering the hydrolysis pattern. Moreover, both the reduced sugar yield and weight loss of the pretreated wheat bran, corn cob, and corn stalk hydrolyzed by XynA for 12 h increased by more than 30 %. These findings are important to better understand the relationship between enzyme activities and their terminal regions and suggest candidate materials for lignocellulosic biomass utilization.

Unveiling the Methyl Transfer Mechanisms in the Epigenetic Machinery DNMT3A-3L: A Comprehensive Study Integrating Assembly Dynamics with Catalytic Reactions

In epigenetic mechanisms, DNA methyltransferase 3 alpha (DNMT3A) acts as an initiator for DNA methylation and prevents the downstream genes from expressing. Perturbations of DNMT3A functions may cause uncontrolled gene expression, resulting in pathogenic consequences such as cancers. It is, therefore, vitally important to understand the catalytic process of DNMT3A in its biological macromolecule assembly, viz., heterotetramer: (DNMT3A-3 L)^dimer. In this study, we utilized molecular dynamics (MD) simulations, Markov State Models (MSM), and quantum mechanics/molecular mechanics simulations (QM/MM) to investigate the de novo methyl transfer process. We identified the dynamics of the key residues relevant to the insertion of the target cytosine (dC) into the catalytic domain of DNMT3A, and the detailed potential energy surface of the seven-step reaction referring to methyl transfer. Our calculated potential energy barrier (22.51 kcal/mol) approximates the former experimental data (23.12 kcal/mol). The conformational change of the 5-methyl-cytosine (5mC) intermediate was found necessary in forming a four-water chain for the elimination step, which is unique to the other DNMTs. The biological assembly facilitates the creation of such a water chain, and the elimination occurs in an asynchronized mechanism in the two catalytic pockets. We anticipate the findings can enable a better understanding of the general mechanisms of the de novo methyl transfer for fulfilling the key enzymatic functions in epigenetics. And the unique elimination of DNMT3A might ignite novel methods for designing anti-cancer and tumor inhibitors of DNMTs.

Theoretical Studies of Mutual Effects between 6-m-r Hemiketalization and 26-m-r Lactonization in Pimaricin Thioesterase

Pimaricin is a small polyene macrolide antibiotic and has been broadly used as an antimycotic and antiprotozoal agent in both humans and foods. As a thioesterase in type-I polyketide synthase, pimTE controls the 26-m-r macrolide main chain release in pimaricin biosynthesis. In this work, we sought to determine whether the 6-m-r hemiketal formation was linked to pimTE-catalyzed 26-m-r lactonization. Compared to non-hemiketal TEs, pimTE is characterized by an aspartic acid residue (D179) accessible to the U-turn motif in the acyl-enzyme intermediate. Both the covalent docking and molecular dynamics simulations demonstrate that the reactive conformations for macrocyclic lactonization are drastically promoted by the 6-m-r hemiketal. Moreover, the small-model quantum mechanistic calculations suggest that protic residues can significantly accelerate the 6-m-r hemiketal cyclization. In addition, the post-hemiketal molecular dynamic simulations demonstrate that hydrogen-bonding networks surrounding the substrate U-turn of the hairpin-shaped conformation changes significantly when the 6-m-r hemiketal is formed. In particular, the R-hemiketal intermediate is not only catalyzed by the D179 residue, but also twists the hairpin structure to the 26-m-r lactonizing pre-reaction state. By contrast, the S-hemiketal formation is unlikely catalyzed by D179, which twists the hairpin in an opposite direction. Our results propose that pimTE could be a bi-functional enzyme, which can synergistically catalyze tandem 6-m-r and 26-m-r formations during the main-chain release of pimaricin biosynthesis.

Evolutionary coupling-inspired engineering of alcohol dehydrogenase reveals the influence of distant sites on its catalytic efficiency for stereospecific synthesis of chiral alcohols

Alcohol dehydrogenase (ADH) has attracted much attention due to its ability to catalyze the synthesis of important chiral alcohol pharmaceutical intermediates with high stereoselectivity. ADH protein engineering efforts have generally focused on reshaping the substrate-binding pocket. However, distant sites outside the pocket may also affect its activity, although the underlying molecular mechanism remains unclear. The current study aimed to apply evolutionary coupling-inspired engineering to the ADH CpRCR and to identify potential mutation sites. Through conservative analysis, phylogenic analysis and residues distribution analysis, the co-evolution hotspots Leu34 and Leu137 were confirmed to be highly evolved under the pressure of natural selection and to be possibly related to the catalytic function of the protein. Hence, Leu34 and Leu137, far away from the active center, were selected for mutation. The generated CpRCR-L34A and CpRCR-L137V variants showed high stereoselectivity and 1.24-7.81 fold increase in k cat /K m value compared with that of the wild type, when reacted with 8 aromatic ketones or β-ketoesters. Corresponding computational study implied that L34 and L137 may extend allosteric fluctuation in the protein structure from the distal mutational site to the active site. Moreover, the L34 and L137 mutations modified the pre-reaction state in multiple ways, in terms of position of the hydride with respect to the target carbonyl. These findings provide insights into the catalytic mechanism of the enzyme and facilitate its regulation from the perspective of the site interaction network.

Both Configuration and QM Region Size Matter: Zinc Stability in QM/MM Models of DNA Methyltransferase

Quantum-mechanical/molecular-mechanical (QM/MM) methods are essential to the study of metalloproteins, but the relative importance of sampling and degree of QM treatment in achieving quantitative predictions is poorly understood. We study the relative magnitude of configurational and QM-region sensitivity of energetic and electronic properties in a representative Zn²⁺ metal binding site of a DNA methyltransferase. To quantify property variations, we analyze snapshots extracted from 250 ns of molecular dynamics simulation. To understand the degree of QM-region sensitivity, we perform analysis using QM regions ranging from a minimal 49-atom region consisting only of the Zn²⁺ metal and its four coordinating Cys residues up to a 628-atom QM region that includes residues within 12 Å of the metal center. Over the configurations sampled, we observe that illustrative properties (e.g., rigid Zn²⁺ removal energy) exhibit large fluctuations that are well captured with even minimal QM regions. Nevertheless, for both energetic and electronic properties, we observe a slow approach to asymptotic limits with similarly large changes in absolute values that converge only with larger (ca. 300-atom) QM region sizes. For the smaller QM regions, the electronic description of Zn²⁺ binding is incomplete: the metal binds too tightly and is too stabilized by the strong electrostatic potential of MM point charges, and the Zn-S bond covalency is overestimated. Overall, this work suggests that efficient sampling with QM/MM in small QM regions is an effective method to explore the influence of enzyme structure on target properties. At the same time, accurate descriptions of electronic and energetic properties require a larger QM region than the minimal metal-coordinating residues in order to converge treatment of both metal-local bonding and the overall electrostatic environment.