Theoretical Insights into the Reaction and Inhibition Mechanism of Metal-Independent Retaining Glycosyltransferase Responsible for Mycothiol Biosynthesis

Unlocking Precision Docking for Metalloproteins

Metalloproteins play a fundamental role in molecular biology, contributing to various biological processes. However, the discovery of high-affinity ligands targeting metalloproteins has been delayed due, in part, to a lack of suitable tools and data. Molecular docking, a widely used technique for virtual screening of small-molecule ligand interactions with proteins, often faces challenges when applied to metalloproteins due to the particular nature of the ligand metal bond. To address these limitations associated with docking metalloproteins, we introduce a knowledge-driven docking approach known as "metalloprotein bias docking" (MBD), which extends the AutoDock Bias technique. We assembled a comprehensive data set of metalloprotein-ligand complexes from 15 different metalloprotein families, encompassing Ca, Co, Fe, Mg, Mn, and Zn metal ions. Subsequently, we conducted a performance analysis of our MBD method and compared it to the conventional docking (CD) program AutoDock4, applied to various metalloprotein targets within our data set. Our results demonstrate that MBD outperforms CD, significantly enhancing accuracy, selectivity, and precision in ligand pose prediction. Additionally, we observed a positive correlation between our predicted ligand free energies and the corresponding experimental values. These findings underscore the potential of MBD as a valuable tool for the effective exploration of metalloprotein-ligand interactions.

Revealing the conformational dynamics of UDP-GlcNAc recognition by O-GlcNAc transferase via Markov state model

The O-linked N-acetylglucosamine (O-GlcNAc) glycosylation is a critical post-translational modification and closely linked to various physiological and pathological conditions. The O-GlcNAc transferase (OGT) functions as the only glycosyltransferase of O-GlcNAc glycosylation by transferring GlcNAc from UDP-GlcNAc to serine or threonine residues on protein substrates. The interaction mode of UDP-GlcNAc against OGT has been preliminarily revealed by the crystal structures, yet an atomic-level comprehension for the conformational dynamics of the recognition process remains elusive. Here, we construct the Markov state model based on extensive all-atom molecular dynamics (MD) simulations with an aggregated simulation time of ∼9 μs, and reveal that the UDP-GlcNAc recognition process by OGT encompasses four key metastable states, occurring within an estimated timescale of ∼10 μs. During UDP-GlcNAc recognition process, we find the pyrophosphate moiety (P2O52-) initially anchors to the active pocket via salt bridge and hydrogen bonds, facilitating subsequent binding of the uridine and GlcNAc moieties. Furthermore, the functional roles of K842 involved in the salt bridge with P2O52- were evaluated through extra mutant MD simulations. Overall, our study provides valuable insights into the UDP-GlcNAc recognition mechanism by OGT, which could further aid in mechanistic studies of O-GlcNAc glycosylation and drug development targeting on OGT.

Best Practices on QM/MM Simulations of Biological Systems

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

How inverting β-1,4-galactosyltransferase-1 can quench a high charge of the by-product UDP3- in catalysis: a QM/MM study of enzymatic reaction with native and UDP-5' thio galactose substrates

The catalysis of inverting glycosyltransferases consists of several biophysical and biochemical processes during which the transfer of a sugar residue from the purine phosphate donor substrate to an acceptor substrate occurs with stereo-inversion of the anomeric C1 center at a product. During catalysis a highly charged phosphate by-product (UDP3-) is formed and a mechanism of how the enzyme stabilizes it back to the UDP2- form is not known. Using methods of molecular modeling (hybrid DFT-QM/MM calculations) we proposed and validated a catalytic mechanism of bovine inverting β-1,4-galactosyltransferase-1 (β4Gal-T1) with native (UDP-galactose) and thio donor substrates (UDP-5' thio galactose). We focused on three aspects of the mechanism not yet investigated: (i) the formation of an oxocarbenium ion intermediate, which was only found for the retaining glycosyltransferases for the time being; (ii) the mechanism of stabilization of a highly charged phosphate by-product (UDP3-) back to its standard in vivo form (UDP2-); (iii) explanation for why in experimental measurements the rate of catalysis with the thio donor substrate is only 8% of the rate of that with the natural substrate. To understand the differences in the interaction patterns between the complexes enzyme : UDP-Gal and enzyme : UDP-5S-Gal, fragmented molecular orbital (FMO) decomposition energy analysis was carried out at the DFT level.

Dissecting the Structural and Chemical Determinants of the "Open-to-Closed" Motion in the Mannosyltransferase PimA from Mycobacteria

The phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential peripheral membrane glycosyltransferase that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides (PIMs), key structural elements and virulence factors of Mycobacterium tuberculosis. PimA undergoes functionally important conformational changes, including (i) α-helix-to-β-strand and β-strand-to-α-helix transitions and (ii) an "open-to-closed" motion between the two Rossmann-fold domains, a conformational change that is necessary to generate a catalytically competent active site. In previous work, we established that GDP-Man and GDP stabilize the enzyme and facilitate the switch to a more compact active state. To determine the structural contribution of the mannose ring in such an activation mechanism, we analyzed a series of chemical derivatives, including mannose phosphate (Man-P) and mannose pyrophosphate-ribose (Man-PP-RIB), and additional GDP derivatives, such as pyrophosphate ribose (PP-RIB) and GMP, by the combined use of X-ray crystallography, limited proteolysis, circular dichroism, isothermal titration calorimetry, and small angle X-ray scattering methods. Although the β-phosphate is present, we found that the mannose ring, covalently attached to neither phosphate (Man-P) nor PP-RIB (Man-PP-RIB), does promote the switch to the active compact form of the enzyme. Therefore, the nucleotide moiety of GDP-Man, and not the sugar ring, facilitates the "open-to-closed" motion, with the β-phosphate group providing the high-affinity binding to PimA. Altogether, the experimental data contribute to a better understanding of the structural determinants involved in the "open-to-closed" motion not only observed in PimA but also visualized and/or predicted in other glycosyltransfeases. In addition, the experimental data might prove to be useful for the discovery and/or development of PimA and/or glycosyltransferase inhibitors.

Leloir Glycosyltransferases in Applied Biocatalysis: A Multidisciplinary Approach

Enzymes are nature’s catalyst of choice for the highly selective and efficient coupling of carbohydrates. Enzymatic sugar coupling is a competitive technology for industrial glycosylation reactions, since chemical synthetic routes require extensive use of laborious protection group manipulations and often lack regio- and stereoselectivity. The application of Leloir glycosyltransferases has received considerable attention in recent years and offers excellent control over the reactivity and selectivity of glycosylation reactions with unprotected carbohydrates, paving the way for previously inaccessible synthetic routes. The development of nucleotide recycling cascades has allowed for the efficient production and reuse of nucleotide sugar donors in robust one-pot multi-enzyme glycosylation cascades. In this way, large glycans and glycoconjugates with complex stereochemistry can be constructed. With recent advances, LeLoir glycosyltransferases are close to being applied industrially in multi-enzyme, programmable cascade glycosylations.

Digalactosyldiacylglycerol Synthase Gene MtDGD1 Plays an Essential Role in Nodule Development and Nitrogen Fixation

Little is known about the genes participating in digalactosyldiacylglycerol (DGDG) synthesis during nodule symbiosis. Here, we identified full-length MtDGD1, a synthase of DGDG, and characterized its effect on symbiotic nitrogen fixation in Medicago truncatula. Immunofluorescence and immunoelectron microscopy showed that MtDGD1 was located on the symbiosome membranes in the infected cells. β-Glucuronidase histochemical staining revealed that MtDGD1 was highly expressed in the infection zone of young nodules as well as in the whole mature nodules. Compared with the control, MtDGD1-RNA interference transgenic plants exhibited significant decreases in nodule number, symbiotic nitrogen fixation activity, and DGDG abundance in the nodules, as well as abnormal nodule and symbiosome development. Overexpression of MtDGD1 resulted in enhancement of nodule number and nitrogen fixation activity. In response to phosphorus starvation, the MtDGD1 expression level was substantially upregulated and the abundance of nonphospholipid DGDG was significantly increased in the roots and nodules, accompanied by corresponding decreases in the abundance of phospholipids such as phosphatidylcholine, phosphatidylethanolamine, and phosphatidylinositol. Overall, our results indicate that DGD1 contributes to effective nodule organogenesis and nitrogen fixation by affecting the synthesis and content of DGDG during symbiosis.

Boronic-Acid-Catalyzed Regioselective and 1,2- cis-Stereoselective Glycosylation of Unprotected Sugar Acceptors via SNi-Type Mechanism

Regio- and 1,2- cis-stereoselective chemical glycosylation of unprotected glycosyl acceptors has been in great demand for the efficient synthesis of natural glycosides. However, simultaneously regulating these selectivities has been a longstanding problem in synthetic organic chemistry. In nature, glycosyl transferases catalyze regioselective 1,2- cis-glycosylations via the S_Ni mechanism, yet no useful chemical glycosylations based on this mechanism have been developed. In this paper, we report a highly regio- and 1,2- cis-stereoselective S_Ni-type glycosylation of 1,2-anhydro donors and unprotected sugar acceptors using p-nitrophenylboronic acid (10e) as a catalyst in the presence of water under mild conditions. Highly controlled regio- and 1,2- cis-stereoselectivities were achieved via the combination of boron-mediated carbohydrate recognition and the S_Ni-type mechanism. Mechanistic studies using the KIEs and DFT calculations were consistent with a highly dissociative concerted S_Ni mechanism. This glycosylation method was applied successfully to the direct glycosylation of unprotected natural glycosides and the efficient synthesis of a complex oligosaccharide with minimal protecting groups.

Chiral Phosphoric Acid Catalyzed Enantioselective Ring Expansion Reaction of 1,3-Dithiane Derivatives: Case Study of the Nature of Ion-Pairing Interaction

Chiral counterion controlled asymmetric catalysis via an ion-pairing interaction has attracted immense attention in recent years. Despite a number of successful studies, the mechanistic elucidation of the stereocontrolling element in the ion-pairing interaction is rarely conducted and hence its nature is still far from being well understood. Herein we report an in-depth mechanistic case study of a newly developed enantioselective ring expansion reaction of 1,3-dithiane derivatives catalyzed by chiral phosphoric acid (CPA). An unprecedented enantioselective 1,2-sulfur rearrangement/stereospecific nucleophilic addition sequence was proven to be the stereoselective pathway. More importantly, by thorough investigation of the intrinsic nature of the stereospecific nucleophilic addition to the cationic thionium intermediate, we discovered that the key interaction in this process is the nonclassical C-H···O hydrogen bonds formed between the conjugate base of the CPA catalyst and the cationic intermediate. These C-H···O hydrogen bonds not only bind the catalyst to the substrates to form energetically favored states throughout the overall processes but also firmly maintain the relative positions of these fragments as the "fixed" contact ion pair to sustain the chiral information generated at the initial sulfur rearrangement step. This mechanistic case study provides a very clear understanding of the nature of the ion-pairing interaction in organocatalysis. The conclusion encourages the further development of the research field with the focus to design new organocatalysts and cultivate novel organocatalytic transformations.