Different QM/MM Approaches To Elucidate Enzymatic Reactions: Case Study on ppGalNAcT2

Engineering the Reaction Pathway of a Non-heme Iron Oxygenase Using Ancestral Sequence Reconstruction

Non-heme iron (FeII), α-ketoglutarate (α-KG)-dependent oxygenases are a family of enzymes that catalyze an array of transformations that cascade forward after the formation of radical intermediates. Achieving control over the reaction pathway is highly valuable and a necessary step toward broadening the applications of these biocatalysts. Numerous approaches have been used to engineer the reaction pathway of FeII/α-KG-dependent enzymes, including site-directed mutagenesis, DNA shuffling, and site-saturation mutagenesis, among others. Herein, we showcase a novel ancestral sequence reconstruction (ASR)-guided strategy in which evolutionary information is used to pinpoint the residues critical for controlling different reaction pathways. Following this, a combinatorial site-directed mutagenesis approach was used to quickly evaluate the importance of each residue. These results were validated using a DNA shuffling strategy and through quantum mechanical/molecular mechanical (QM/MM) simulations. Using this approach, we identified a set of active site residues together with a key hydrogen bond between the substrate and an active site residue, which are crucial for dictating the dominant reaction pathway. Ultimately, we successfully converted both extant and ancestral enzymes that perform benzylic hydroxylation into variants that can catalyze an oxidative ring-expansion reaction, showcasing the potential of utilizing ASR to accelerate the reaction pathway engineering within enzyme families that share common structural and mechanistic features.

Glycosylation Modulates the Structure and Functions of Collagen: A Review

Collagens are fundamental constituents of the extracellular matrix and are the most abundant proteins in mammals. Collagens belong to the family of fibrous or fiber-forming proteins that self-assemble into fibrils that define their mechanical properties and biological functions. Up to now, 28 members of the collagen superfamily have been recognized. Collagen biosynthesis occurs in the endoplasmic reticulum, where specific post-translational modification-glycosylation-is also carried out. The glycosylation of collagens is very specific and adds β-d-galactopyranose and β-d-Glcp-(1→2)-d-Galp disaccharide through β-O-linkage to hydroxylysine. Several glycosyltransferases, namely COLGALT1, COLGALT2, LH3, and PGGHG glucosidase, were associated the with glycosylation of collagens, and recently, the crystal structure of LH3 has been solved. Although not fully understood, it is clear that the glycosylation of collagens influences collagen secretion and the alignment of collagen fibrils. A growing body of evidence also associates the glycosylation of collagen with its functions and various human diseases. Recent progress in understanding collagen glycosylation allows for the exploitation of its therapeutic potential and the discovery of new agents. This review will discuss the relevant contributions to understanding the glycosylation of collagens. Then, glycosyltransferases involved in collagen glycosylation, their structure, and catalytic mechanism will be surveyed. Furthermore, the involvement of glycosylation in collagen functions and collagen glycosylation-related diseases will be discussed.

Modeling Ensembles of Enzyme Reaction Pathways with Hi-MSM Reveals the Importance of Accounting for Pathway Diversity

Conventional quantum mechanical-molecular mechanics (QM/MM) simulation approaches for modeling enzyme reactions often assume that there is one dominant reaction pathway and that this pathway can be sampled starting from an X-ray structure of the enzyme. These assumptions reduce computational cost; however, their validity has not been extensively tested. This is due in part to the lack of a rigorous formalism for integrating disparate pathway information from dynamical QM/MM calculations. Here, we present a way to model ensembles of reaction pathways efficiently using a divide-and-conquer strategy through Hierarchical Markov State Modeling (Hi-MSM). This approach allows information on multiple, distinct pathways to be incorporated into a chemical kinetic model, and it allows us to test these two assumptions. Applying Hi-MSM to the reaction carried out by dihydrofolate reductase (DHFR) we find (i) there are multiple, distinct pathways significantly contributing to the overall flux of the reaction that the conventional approach does not identify and (ii) that the conventional approach does not identify the dominant reaction pathway. Thus, both assumptions underpinning the conventional approach are violated. Since DHFR is a relatively small enzyme, and configuration space scales exponentially with protein size, accounting for multiple reaction pathways is likely to be necessary for most enzymes.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

The catalytic reaction mechanism of tyrosylprotein sulfotransferase-1

Tyrosine sulfation alters the biological activity of many proteins involved in different physiological and pathophysiological conditions, such as non-specific immune reaction, response to inflammation and ischemia, targeting of leukocytes and stem cells, or the formation of cancer metastases. Tyrosine sulfation is catalyzed by the enzymes tyrosylprotein sulfotransferases (TPST). In this study, we used QM/MM Car-Parrinello metadynamics simulations together with QM/MM potential energy calculations to investigate the catalytic mechanism of isoform TPST-1. The structural changes along the reaction coordinate are analyzed and discussed. Furthermore, both the methods supported the S_N2 type of catalytic mechanism. The reaction barrier obtained from CPMD metadynamics was 12.8 kcal mol^-1, and the potential energy scan led to reaction barriers of 11.6 kcal mol^-1 and 13.7 kcal mol^-1 with the B3LYP and OPBE functional, respectively. The comparison of the two methods (metadynamics and potential energy scan) may be helpful for future mechanistic studies. The insight into the reaction mechanism of TPST-1 might help with the rational design of transition-state TPST inhibitors.

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.

Selectins-The Two Dr. Jekyll and Mr. Hyde Faces of Adhesion Molecules-A Review

Selectins belong to a group of adhesion molecules that fulfill an essential role in immune and inflammatory responses and tissue healing. Selectins are glycoproteins that decode the information carried by glycan structures, and non-covalent interactions of selectins with these glycan structures mediate biological processes. The sialylated and fucosylated tetrasaccharide sLe^x is an essential glycan recognized by selectins. Several glycosyltransferases are responsible for the biosynthesis of the sLe^x tetrasaccharide. Selectins are involved in a sequence of interactions of circulated leukocytes with endothelial cells in the blood called the adhesion cascade. Recently, it has become evident that cancer cells utilize a similar adhesion cascade to promote metastases. However, like Dr. Jekyll and Mr. Hyde’s two faces, selectins also contribute to tissue destruction during some infections and inflammatory diseases. The most prominent function of selectins is associated with the initial stage of the leukocyte adhesion cascade, in which selectin binding enables tethering and rolling. The first adhesive event occurs through specific non-covalent interactions between selectins and their ligands, with glycans functioning as an interface between leukocytes or cancer cells and the endothelium. Targeting these interactions remains a principal strategy aimed at developing new therapies for the treatment of immune and inflammatory disorders and cancer. In this review, we will survey the significant contributions to and the current status of the understanding of the structure of selectins and the role of selectins in various biological processes. The potential of selectins and their ligands as therapeutic targets in chronic and acute inflammatory diseases and cancer will also be discussed. We will emphasize the structural characteristic of selectins and the catalytic mechanisms of glycosyltransferases involved in the biosynthesis of glycan recognition determinants. Furthermore, recent achievements in the synthesis of selectin inhibitors will be reviewed with a focus on the various strategies used for the development of glycosyltransferase inhibitors, including substrate analog inhibitors and transition state analog inhibitors, which are based on knowledge of the catalytic mechanism.