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.

How E-, L-, and P-Selectins Bind to sLex and PSGL-1: A Quantification of Critical Residue Interactions.

Selectins and their ability to interact with specific ligands are a cornerstone in cell communication. Over the last three decades, a considerable wealth of experimental and molecular modeling insights into their structure and modus operandi were gathered. Nonetheless, explaining the role of individual selectin residues on a quantitative level remained elusive, despite its importance in understanding the structure-function relationship in these molecules and designing their inhibitors. This work explores essential interactions of selectin-ligand binding, employing a multiscale approach that combines molecular dynamics, quantum-chemical calculations, and residue interaction network models. Such an approach successfully reproduces most of the experimental findings. It proves to be helpful, with the potential for becoming an established tool for quantitative predictions of residue contribution to the binding of biomolecular complexes. The results empower us to quantify the importance of particular residues and functional groups in the protein-ligand interface and to pinpoint differences in molecular recognition by the three selectins. We show that mutations in the E-, L-, and P-selectins, e.g., different residues in positions 46, 85, 97, and 107, present a crucial difference in how the ligand is engaged. We assess the role of sulfation of tyrosine residues in PSGL-1 and suggest that TyrSO3- in position 51 interacting with Arg85 in P-selectin is a significant factor in the increased affinity of P-selectin to PSGL-1 compared to E- and L-selectins. We propose an original pharmacophore targeting five essential PSGL-binding sites based on the analysis of the selectin···PSGL-1 interactions.

Pan-selectin inhibitors as potential therapeutics for COVID-19 treatment: in silico screening study.

Coronavirus disease 2019 (COVID-19) has spread rapidly throughout the globe. The spectrum of disease is broad but among hospitalized patients with COVID-19, respiratory failure from acute respiratory distress syndrome is the leading cause of mortality. There is an urgent need for an effective treatment. The current focus has been developing novel therapeutics, including antivirals, protease inhibitors, vaccines and targeting the overactive cytokine response with anti-cytokine therapy. The overproduction of early response proinflammatory cytokines results in what has been described as a "cytokine storm" is leading eventually to death when the cells fail to terminate the inflammatory response. Accumulating evidence shows that inflammatory cytokines induce selectin ligands that play a crucial role in the pathogenesis of inflammatory diseases by mediating leukocyte migration from the blood into the tissue. Thus, the selectins and selectin ligands represent a promising therapeutic target for the treatment of COVID-19. In this paper, potential pan-selectin inhibitors were identified employing a virtual screening using a docking procedure. For this purpose, the Asinex and ZINC databases of ligands, including approved drugs, biogenic compounds and glycomimetics, altogether 923,602 compounds, were screened against the P-, L- and E-selectin. At first, the experimentally confirmed inhibitors were docked into all three selectins' carbohydrate recognition domains to assess the suitability of the screening procedure. Finally, based on the evaluation of ligands binding, we propose 10 purchasable pan-selectin inhibitors to develop COVID-19 therapeutics.

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 SN2 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.

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 sLex is an essential glycan recognized by selectins. Several glycosyltransferases are responsible for the biosynthesis of the sLex 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.

How Mycobacterium tuberculosis Galactofuranosyl Transferase†2 (GlfT2) Generates Alternating ß-(1-6) and ß-(1-5) Linkages: A QM/MM Molecular Dynamics Study of the Chemical Steps.

Mycobacterium tuberculosis features a unique cell wall that protects the bacterium from the external environment. Disruption of the cell wall assembly is a promising direction for novel anti-tuberculotic drugs. A key component of the cell wall is galactan, a polysaccharide chain composed of galactofuranose (Galf) units connected by alternating ß-(1-5) and ß-(1-6) linkages. The majority of the galactan chain is biosynthesized by a bifunctional enzyme-galactofuranosyl transferase 2 (GlfT2). GlfT2 catalyzes two reactions: the formation of ß-(1-5) and ß-(1-6) linkages. It was suggested that the enzyme acts through a processive mechanism until it adds 30-35 Galf units in a single active site. We applied a QM/MM string method coupled with ab initio molecular dynamics simulations to study the two reactions catalyzed by GlfT2. We showed that both reactions proceed very similarly and feature similar transition-state structures. We also present novel information about the ring puckering behavior of the five-membered furanose ring during the glycosyltransferase reaction and a calculated transition-state structure with galactose in a furanose form that may be used as a guide for the rational design of very specific and extremely potent inhibitors, that is, transition-state analogues, for GlfT2. Due to the absence of a furanose form of galactose in humans, transition-state-analogous inhibitors represent an attractive scaffold for the development of novel antibacterial drugs.

Three-dimensional homology model of GlcNAc-TV glycosyltransferase.

The enzyme UDP-N-acetylglucosamine: α-d-mannoside ß-1-6 N-acetylglucosaminyltransferase V (GnT-V) catalyzes the transfer of GlcNAc from the UDP-GlcNAc donor to the α-1-6-linked mannose of the trimannosyl core structure of glycoproteins to produce the ß-1-6-linked branching of N-linked oligosaccharides. ß-1-6-GlcNAc-branched N-glycans are associated with cancer growth and metastasis. Therefore, the inhibition of GnT-V represents a key target for anti-cancer drug development. However, the development of potent and specific inhibitors of GnT-V is hampered by the lack of information on the three-dimensional structure of the enzyme and on the binding characteristics of its substrates. Here we present the first 3D structure of GnT-V as a result of homology modeling. Various alignment methods, docking the donor and acceptor substrates, and molecular dynamics simulation were used to construct seven homology models of GnT-V and characterize the binding of its substrates. The best homology model is consistent with available experimental data. The three-dimensional model, the structure of the enzyme catalytic site and binding information obtained for the donor and acceptor can be useful in studies of the catalytic mechanism and design of inhibitors of GnT-V.

Stepwise catalytic mechanism via short-lived intermediate inferred from combined QM/MM MERP and PES calculations on retaining glycosyltransferase ppGalNAcT2.

The glycosylation of cell surface proteins plays a crucial role in a multitude of biological processes, such as cell adhesion and recognition. To understand the process of protein glycosylation, the reaction mechanisms of the participating enzymes need to be known. However, the reaction mechanism of retaining glycosyltransferases has not yet been sufficiently explained. Here we investigated the catalytic mechanism of human isoform 2 of the retaining glycosyltransferase polypeptide UDP-GalNAc transferase by coupling two different QM/MM-based approaches, namely a potential energy surface scan in two distance difference dimensions and a minimum energy reaction path optimisation using the Nudged Elastic Band method. Potential energy scan studies often suffer from inadequate sampling of reactive processes due to a predefined scan coordinate system. At the same time, path optimisation methods enable the sampling of a virtually unlimited number of dimensions, but their results cannot be unambiguously interpreted without knowledge of the potential energy surface. By combining these methods, we have been able to eliminate the most significant sources of potential errors inherent to each of these approaches. The structural model is based on the crystal structure of human isoform 2. In the QM/MM method, the QM region consists of 275 atoms, the remaining 5776 atoms were in the MM region. We found that ppGalNAcT2 catalyzes a same-face nucleophilic substitution with internal return (SNi). The optimized transition state for the reaction is 13.8 kcal/mol higher in energy than the reactant while the energy of the product complex is 6.7 kcal/mol lower. During the process of nucleophilic attack, a proton is synchronously transferred to the leaving phosphate. The presence of a short-lived metastable oxocarbenium intermediate is likely, as indicated by the reaction energy profiles obtained using high-level density functionals.

Theoretical study of enzymatic catalysis explains why the trapped covalent intermediate in the E303C mutant of glycosyltransferase GTB was not detected in the wild-type enzyme.

Hybrid quantum mechanics/molecular mechanics calculations were used to study the catalytic mechanism of the retaining human α-(1,3)-galactosyltransferase (GTBWT) and its E303C mutant (GTBE303C). Both backside (via covalent glycosyl-enzyme intermediate, CGEI) and frontside SNi-like mechanisms (via oxocarbenium-ion intermediate, OCII) were investigated. The calculations suggest that both mechanisms are feasible in the enzymatic catalysis. The nucleophilic attack of the acceptor substrate to the anomeric carbon of OCII is the rate-determining step with an overall reaction barrier (ΔE(‡) = 19.5 kcal mol(-1)) in agreement with an experimental rate constant (kcat = 5.1 s(-1)). A calculated α-secondary kinetic isotope effect (α-KIE) of 1.27 (GTBWT) and 1.26 (GTBE303C) predicts dissociative character of the transition state in agreement with experimentally measured α-KIE of other retaining glycosyltransferases. Remarkably, stable CGEI in GTBE303C compared with its counterpart in GTBWT may explain why the CGEI has been detected by mass spectrometry only in GTBE303C ( Soya N, Fang Y, Palcic MM, Klassen JS. 2011. Trapping and characterization of covalent intermediates of mutant retaining glycosyltransferases. Glycobiology, 21: 547-552).

Ab initio modelling of the anomeric and exo anomeric effects in 2-methoxytetrahydropyran and 2-methoxythiane corrected for intramolecular BSSE.

Accurate ab initio calculations including basis set limit (BSL) extrapolations, removal of intramolecular basis set superposition error (BSSE), solvent effect corrections, and thermal effects have been carried out to compare the structure and the anomeric and exo-anomeric effect in 2-methoxytetrahydropyran and 2-methoxythiane. The effect of intramolecular BSSE on the energetics was outlined for the first time in these types of compounds. It was found that both title compounds show comparable behaviour with respect to BSSE. The energy gap between the axial and equatorial form of 2-methoxythiane is reduced by 0.23 kcal mol(-1) due to the BSSE correction at the MP2/aug-cc-pVTZ level of theory, and in 2-methoxytetrahydropyran it is reduced by 0.21 kcal mol(-1). The intramolecular BSSE influenced also the energy differences between the gauche and trans conformers in both compounds. Energy decomposition analysis (EDA) reveals that the dominant destabilising interaction is repulsion and its primary stabilizing counterpart is the polarization interaction.

A theoretical study on the catalytic mechanism of the retaining α-1,2-mannosyltransferase Kre2p/Mnt1p: the impact of different metal ions on catalysis.

Glycosyltransferases are sugar-processing enzymes that require a specific metal ion cofactor for catalysis. In the presence of other ions the catalysis is often impaired. Here, for the first time, the enzymatic catalysis in the presence of various metal ions was modeled for a glycosyltransferase using a large enzymatic model. The catalytic mechanism of α-1,2-mannosyltransferase Kre2p/Mnt1p in the presence of Mn(2+) and other ions (Mg(2+), Zn(2+) and Ca(2+)) was modeled at the two hybrid DFT-QM/MM (M06-2X/OPLS2005 and B3LYP/OPLS2005) levels. Kinetic and structural parameters of transition states and intermediates, as well as kinetic isotope effects, were predicted and compared with available experimental and theoretical data. The catalysis in the presence of the metal ions is predicted as a stepwise SNi-like nucleophilic substitution reaction (DNint*AN(‡)DhAxh) via oxocarbenium ion intermediates. In the rate-determining step the leaving phosphate group of the donor substrate plays a role of the base catalyst. The predicted increased enzymatic reactivity (kcat: Zn(2+) ≈ Mg(2+) < Mn(2+) < Ca(2+)) correlated with the metal ion ability to polarize the Kre2p environment (Mg(2+) > Zn(2+) > Mn(2+) > Ca(2+)). The formation of the retained anomeric configuration in the product is controlled by a strict geometry of the active site of Kre2p. The 6-OH group of the attacking acceptor substrate may assist in protection of the anomeric carbon against unwanted hydrolysis by a through-space interaction with the electron deficient C1[double bond, length as m-dash]O5(+) moiety of the oxocarbenium-ion-like transition state.

A QM/MM investigation of the catalytic mechanism of metal-ion-independent core 2 ß1,6-N-acetylglucosaminyltransferase.

ß1,6-GlcNAc-transferase (C2GnT) is an important controlling factor of biological functions for many glycoproteins and its activity has been found to be altered in breast, colon, and lung cancer cells, in leukemia cells, in the lymhomonocytes of multiple sclerosis patients, leukocytes from diabetes patients, and in conditions causing an immune deficiency. The result of the action of C2GnT is the core 2 structure that is essential for the further elongation of the carbohydrate chains of O-glycans. The catalytic mechanism of this metal-ion-independent glycosyltransferase is of paramount importance and is investigated here by using quantum mechanical (QM) (density functional theory (DFT))/molecular modeling (MM) methods with different levels of theory. The structural model of the reaction site used in this report is based on the crystal structures of C2GnT. The entire enzyme-substrate system was subdivided into two different subsystems: the QM subsystem containing 206 atoms and the MM region containing 5914 atoms. Three predefined reaction coordinates were employed to investigate the catalytic mechanism. The calculated potential energy surfaces discovered the existence of a concerted SN 2-like mechanism. In this mechanism, a nucleophilic attack by O6 facilitated by proton transfer to the catalytic base and the separation of the leaving group all occur almost simultaneously. The transition state for the proposed reaction mechanism at the M06-2X/6-31G** (with diffuse functions on the O1', O5', OGlu , and O6 atoms) level was located at C1-O6=1.74 Šand C1-O1=2.86 Å. The activation energy for this mechanism was estimated to be between 20 and 29 kcal mol⁻¹, depending on the method used. These calculations also identified a low-barrier hydrogen bond between the nucleophile O6H and the catalytic base Glu320, and a hydrogen bond between the N-acetamino group and the glycosidic oxygen of the donor in the TS. It is proposed that these interactions contribute to a stabilization of TS and participate in the catalytic mechanism.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review.

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Binding of ß-amyloid (1-42) peptide to negatively charged phospholipid membranes in the liquid-ordered state: modeling and experimental studies.

To explore the initial stages of amyloid ß peptide (Aß42) deposition on membranes, we have studied the interaction of Aß42 in the monomeric form with lipid monolayers and with bilayers in either the liquid-disordered or the liquid-ordered (L(o)) state, containing negatively charged phospholipids. Molecular dynamics (MD) simulations of the system have been performed, as well as experimental measurements. For bilayers in the L(o) state, in the absence of the negatively charged lipids, interaction is weak and it cannot be detected by isothermal calorimetry. However, in the presence of phosphatidic acid, or of cardiolipin, interaction is detected by different methods and in all cases interaction is strongest with lower (2.5-5 mol%) than higher (10-20 mol%) proportions of negatively charged phospholipids. Liquid-disordered bilayers consistently allowed a higher Aß42 binding than L(o) ones. Thioflavin T assays and infrared spectroscopy confirmed a higher proportion of ß-sheet formation under conditions when higher peptide binding was measured. The experimental results were supported by MD simulations. We used 100 ns MD to examine interactions between Aß42 and three different 512 lipid bilayers consisting of palmitoylsphingomyelin, dimyristoyl phosphatidic acid, and cholesterol in three different proportions. MD pictures are different for the low- and high-charge bilayers, in the former case the peptide is bound through many contact points to the bilayer, whereas for the bilayer containing 20 mol% anionic phospholipid only a small fragment of the peptide appears to be bound. The MD results indicate that the binding and fibril formation on the membrane surface depends on the composition of the bilayer, and is the result of a subtle balance of many inter- and intramolecular interactions between the Aß42 and membrane.

Substrate-assisted catalytic mechanism of O-GlcNAc transferase discovered by quantum mechanics/molecular mechanics investigation.

In higher eukaryotes, a variety of proteins are post-translationally modified by adding O-linked N-acetylglucosamine (GlcNAc) residue to serine or threonine residues. Misregulation of O-GlcNAcylation is linked to a wide variety of diseases, such as diabetes, cancer, and neurodegenerative diseases, including Alzheimer's disease. GlcNAc transfer is catalyzed by an inverting glycosyltransferase O-GlcNAc transferase (uridine diphospho-N-acetylglucosamine:polypeptide ß-N-acetylaminyltransferase, OGT) that belongs to the GT-B superfamily. The catalytic mechanism of this metal-independent glycosyltransferase is of primary importance and is investigated here using QM(DFT)/MM methods. The structural model of the reaction site used in this paper is based on the crystal structures of OGT. The entire enzyme-substrate system was partitioned into two different subsystems: the QM subsystem containing 198 atoms, and the MM region containing 11,326 atoms. The catalytic mechanism was monitored by means of three two-dimensional potential energy maps calculated as a function of three predefined reaction coordinates at different levels of theory. These potential energy surfaces revealed the existence of a concerted S(N)2-like mechanism, in which a nucleophilic attack by O(Ser), facilitated by proton transfer to the catalytic base, and the dissociation of the leaving group occur almost simultaneously. The transition state for the proposed reaction mechanism at the MPW1K level was located at C1-O(Ser) = 1.92 Å and C1-O1 = 3.11 Å. The activation energy for this passage was estimated to be ~20 kcal mol(-1). These calculations also identified, for the first time for glycosyltransferases, the substrate-assisted mechanism in which the N-acetamino group of the donor participates in the catalytic mechanism.

Conformational free energy modeling of druglike molecules by metadynamics in the WHIM space.

Protein-ligand affinities can be significantly influenced not only by the interaction itself but also by conformational equilibrium of both binding partners, free ligand and free protein. Identification of important conformational families of a ligand and prediction of their thermodynamics is important for efficient ligand design. Here we report conformational free energy modeling of nine small-molecule drugs in explicitly modeled water by metadynamics with a bias potential applied in the space of weighted holistic invariant molecular (WHIM) descriptors. Application of metadynamics enhances conformational sampling compared to unbiased molecular dynamics simulation and allows to predict relative free energies of key conformations. Selected free energy minima and one example of transition state were tested by a series of unbiased molecular dynamics simulation. Comparison of free energy surfaces of free and target-bound Imatinib provides an estimate of free energy penalty of conformational change induced by its binding to the target.

N-(4-Substituted-benzoyl)-N'-(ß-d-glucopyranosyl)ureas as inhibitors of glycogen phosphorylase: Synthesis and evaluation by kinetic, crystallographic, and molecular modelling methods.

N-(4-Substituted-benzoyl)-N'-(ß-d-glucopyranosyl) ureas (substituents: Me, Ph, Cl, OH, OMe, NO(2), NH(2), COOH, and COOMe) were synthesised by ZnCl(2) catalysed acylation of O-peracetylated ß-d-glucopyranosyl urea as well as in reactions of O-peracetylated or O-unprotected glucopyranosylamines and acyl-isocyanates. O-deprotections were carried out by base or acid catalysed transesterifications where necessary. Kinetic studies revealed that most of these compounds were low micromolar inhibitors of rabbit muscle glycogen phosphorylase b (RMGPb). The best inhibitor was the 4-methylbenzoyl compound (K(i)=2.3µM). Crystallographic analyses of complexes of several of the compounds with RMGPb showed that the analogues exploited, together with water molecules, the available space at the ß-pocket subsite and induced a more extended shift of the 280s loop compared to RMGPb in complex with the unsubstituted benzoyl urea. The results suggest the key role of the water molecules in ligand binding and structure-based ligand design. Molecular docking study of selected inhibitors was done to show the ability of the binding affinity prediction. The binding affinity of the highest scored docked poses was calculated and correlated with experimentally measured K(i) values. Results show that correlation is high with the R-squared (R(2)) coefficient over 0.9.

Binding of ß-D-glucosides and ß-D-mannosides by rice and barley ß-D-glycosidases with distinct substrate specificities.

Predominantly, rice Os3BGlu7 operates as a ß-d-glucosidase (EC 3.2.1.21), while barley HvBII acts as a ß-d-mannosidase (EC 3.2.1.25). Saturation transfer difference nuclear magnetic resonance (STD NMR) and transferred nuclear Overhauser effect (trNOE) spectroscopy in conjunction with quantum mechanics/molecular mechanics (QM/MM) modeling and docking at the 6-31+G* level were used to investigate binding of S- and O-linked gluco- and manno-configured aryl-ß-d-glycosides to Os3BGlu7 and HvBII. Kinetic analyses with 4-nitrophenyl ß-d-thioglucoside (4NP-S-Glc) and 4-nitrophenyl ß-d-thiomannoside (4NP-S-Man) indicated that the inhibitions were competitive with apparent K(i) constants of 664 and 710 µM for Os3BGlu7 and 95 and 266 µM for HvBII, respectively. The STD NMR and trNOESY experiments revealed that 4NP-S-Glc and 4NP-S-Man bound weakly in (4)C(1) conformations to Os3BGlu7; 4NP-S-Glc adopted (3)S(5) (B(3,O)) or (1)S(3) ((1,4)B) conformations, and 4NP-S-Man preferred (4)C(1) geometry, when bound to HvBII. The QM modeling and docking, based on GLIDE scores, predicted that 4NP-O-Glc, 4NP-O-Man, and 4NP-S-Man bound preferentially in (1)S(3) geometries to both enzymes, contrary to 4NP-S-Glc that could also adopt a (4)C(1) conformation, although in a "flipped-down" ring position. The experimental and computational data suggested that in glycoside recognition and substrate specificity of Os3BGlu7 and HvBII, a combination of the following determinants is likely to play key roles: (i) the inherent conformational and spatial flexibilities of gluco- and manno-configured substrates in the enzymes' active sites, (ii) the subtle differences in the spatial disposition of active site residues and their capacities to form interactions with specific groups of substrates, and (iii) the small variations in the charge distributions and shapes of the catalytic sites.

Catalytic Mechanism of Processive GlfT2: Transition Path Sampling Investigation of Substrate Translocation.

We applied the transition path sampling (TPS) method to study the translocation step of the catalytic mechanism of galactofuranosyl transferase 2 (GlfT2). Using TPS in the field of enzymatic reactions is still relatively rare, and we show its effectiveness on this enzymatic system. We decipher an unknown mechanism of the translocation step and, thus, provide a complete understanding of the catalytic mechanism of GlfT2 at the atomistic level. The GlfT2 enzyme is involved in the formation of the mycobacterial cell wall and transfers galactofuranose (Galf) from UDP-Galf onto a growing acceptor Galf chain. The biosynthesis of the galactan chain is accomplished in a processive manner, with the growing acceptor substrate remaining bound to GlfT2. The glycosidic bond formed by GlfT2 between the two Galf residues alternates between ß-(1-6) and ß-(1-5) linkages. The translocation of the growing galactan between individual additions of Galf residues is crucial for the function of GlfT2. Analysis of unbiased trajectory ensembles revealed that the translocation proceeds differently depending on the glycosidic linkage between the last two Galf residues. We also showed that the protonation state of the catalytic residue Asp372 significantly influences the translocation. Approximate transition state structures and potential energy reaction barriers of the translocation process were determined. The calculated potential reaction barriers in the range of 6-14 kcal/mol show that the translocation process is not the rate-limiting step in galactan biosynthesis.

Conformational free energy surface of alpha-N-acetylneuraminic acid: an interplay between hydrogen bonding and solvation.

The conformational free energy surface of alpha-N-acetylneuraminic acid (Neu5Ac, sialic acid) in the space of ring-puckering coordinates was calculated using the metadynamics method. Free energy surfaces in vacuum and with an explicit solvent were calculated in GLYCAM 06 force field. In vacuum three structures are almost equivalently populated, namely, the (2)C(5) chair and the B(3,6)/(2)S(6) and (O)S(3) boat/skew-boat conformations. The B(3,6)/(2)S(6) structure is stabilized by an ionic hydrogen bond between the amide N-H bond and the carboxylic group. However, this structure is unfavorable in a water environment in which the experimentally observed (2)C(5) chair conformation is predicted to be more stable than the other structures. These results indicate that environment significantly influences conformation of Neu5Ac and that Neu5Ac-processing enzymes might modify a conformation of their substrates solely by a changing polarity of the environment. The structure of Neu5Ac bound in influenza neuraminidase ((4)S(2)/B(2,5)) belongs to conformations preferred in a water environment. The free energy penalty of this conformational change was calculated (relative to (2)C(5)) as 10.2 +/- 2.0 and 17.3 +/- 2.0 kJ/mol for (4,O)B/(O)S(3) and (4)S(2), respectively. This result indicates that mimicking of the enzyme-bound conformation is likely to be a viable strategy for the design of neuraminidase inhibitors.