QM/MM molecular dynamics study of the galactopyranose → galactofuranose reaction catalysed by Trypanosoma cruzi UDP-galactopyranose mutase

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5′-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2′-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2′-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.

Molecular Dynamics Simulations of Substrate Release from Trypanosoma cruzi UDP-Galactopyranose Mutase

The enzyme UDP-galactopyranose mutase (UGM) represents a promising drug target for the treatment of infections with Trypanosoma cruzi. We have computed the Potential of Mean Force for the release of UDP-galactopyranose from UGM, using Umbrella Sampling simulations. The simulations revealed the conformational changes that both substrate and enzyme undergo during the process. It was determined that the galactopyranose portion of the substrate is highly mobile and that the opening/closing of the active site occurs in stages. Previously uncharacterized interactions with highly conserved residues were also identified. These findings provide new pieces of information that contribute to the rational design of drugs against T. cruzi.

Steric Control of the Rate-Limiting Step of UDP-Galactopyranose Mutase

Galactose is an abundant monosaccharide found exclusively in mammals as galactopyranose (Gal p), the six-membered ring form of this sugar. In contrast, galactose appears in many pathogenic microorganisms as the five-membered ring form, galactofuranose (Gal f). Gal f biosynthesis begins with the conversion of UDP-Gal p to UDP-Gal f catalyzed by the flavoenzyme UDP-galactopyranose mutase (UGM). Because UGM is essential for the survival and proliferation of several pathogens, there is interest in understanding the catalytic mechanism to aid inhibitor development. Herein, we have used kinetic measurements and molecular dynamics simulations to explore the features of UGM that control the rate-limiting step (RLS). We show that UGM from the pathogenic fungus Aspergillus fumigatus also catalyzes the isomerization of UDP-arabinopyranose (UDP-Ara p), which differs from UDP-Gal p by lacking a -CH₂-OH substituent at the C5 position of the hexose ring. Unexpectedly, the RLS changed from a chemical step for the natural substrate to product release with UDP-Ara p. This result implicated residues that contact the -CH₂-OH of UDP-Gal p in controlling the mechanistic path. The mutation of one of these residues, Trp315, to Ala changed the RLS of the natural substrate to product release, similar to the wild-type enzyme with UDP-Ara p. Molecular dynamics simulations suggest that steric complementarity in the Michaelis complex is responsible for this distinct behavior. These results provide new insight into the UGM mechanism and, more generally, how steric factors in the enzyme active site control the free energy barriers along the reaction path.

Multiple functionalities of reduced flavin in the non-redox reaction catalyzed by UDP-galactopyranose mutase

Flavin cofactors are widely used by enzymes to catalyze a broad range of chemical reactions. Traditionally, flavins in enzymes are regarded as redox centers, which enable enzymes to catalyze the oxidation or reduction of substrates. However, a new class of flavoenzyme has emerged over the past quarter century in which the flavin functions as a catalytic center in a non-redox reaction. Here we introduce the unifying concept of flavin hot spots to understand and categorize the mechanisms and reactivities of both traditional and noncanonical flavoenzymes. The major hot spots of reactivity include the N5, C4a, and C4O atoms of the isoalloxazine, and the 2' hydroxyl of the ribityl chain. The role of hot spots in traditional flavoenzymes, such as monooxygenases, is briefly reviewed. A more detailed description of flavin hot spots in noncanonical flavoenzymes is provided, with a focus on UDP-galactopyranose mutase, where the N5 functions as a nucleophile that attacks the anomeric carbon atom of the substrate. Recent results from mechanistic enzymology, kinetic crystallography, and computational chemistry provide a complete picture of the chemical mechanism of UDP-galactopyranose mutase.

In Crystallo Capture of a Covalent Intermediate in the UDP-Galactopyranose Mutase Reaction

UDP-galactopyranose mutase (UGM) plays an essential role in galactofuranose biosynthesis in pathogens by catalyzing the conversion of UDP-galactopyranose to UDP-galactofuranose. Here we report the first crystal structure of a covalent intermediate in the UGM reaction. The 2.3 Å resolution structure reveals UDP bound in the active site and galactopyranose linked to the FAD through a covalent bond between the anomeric C of galactopyranose and N5 of the FAD. The structure confirms the role of the flavin as nucleophile and supports the hypothesis that the proton destined for O5 of galactofuranose is shuttled from N5 of the FAD via O4 of the FAD.

Insights into the mechanism of oxidation of dihydroorotate to orotate catalysed by human class 2 dihydroorotate dehydrogenase: a QM/MM free energy study

The dihydroorotate dehydrogenase (DHOD) enzyme catalyzes the unique redox reaction in the de novo pyrimidine biosynthesis pathway. In this reaction, the oxidation of dihydroorotate (DHO) to orotate (OA) and reduction of the flavin mononucleotide (FMN) cofactor is catalysed by DHOD. The class 2 DHOD, to which the human enzyme belongs, was experimentally shown to follow a stepwise mechanism but the data did not allow the determination of the order of bond-breaking in a stepwise oxidation of DHO. The goal of this study is to understand the reaction mechanism at the molecular level of class 2 DHOD, which may aid in the design of inhibitors that selectively impact the activity of only certain members of the enzyme family. In this paper, the catalytic mechanism of oxidation of DHO to OA in human DHOD was studied using a hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) approach and Molecular Dynamics (MD) simulations. The free energy barriers calculated reveal that the mechanism in human DHOD occurs via a stepwise reaction pathway. In the first step, a proton is abstracted from the C5 of DHO to the deprotonated Ser215 side chain. Whereas, in the second step, the transfer of the hydride or hydride equivalent from the C6 of DHO to the N5 of FMN, where free energy barrier calculated by the DFT/MM level is 10.84 kcal mol(-1). Finally, a residual decomposition analysis was carried out in order to elucidate the influence of the catalytic region residues during DHO oxidation.