Conformational flexibility of N-glycans in solution studied by REMD simulations

Impact of Polarization on the Ring Puckering Dynamics of Hexose Monosaccharides

Analysis of crystal structures of hexose monosaccharides α-d-mannose (α-MAN), β-d-mannose (β-MAN), α-d-glucose (α-GLC), β-d-glucose (β-GLC), α-d-galactose (α-GAL), β-d-galactose (β-GAL), α-d-altrose (α-ALT), β-d-altrose (β-ALT), α-d-idose (α-IDO), and β-d-idose (β-IDO) reveals that the monosaccharide ring adopts multiple ring conformations. These ring conformations can be broadly classified as chair, half-chair, envelope, boat, and skew-boat conformations. The ability of the monosaccharide ring to adopt multiple conformations has been closely tied with their bioactivity. However, it has been difficult to capture the dynamic information of these conformations from experimental studies. Even from simulations, capturing these different conformations is challenging because of the energy barriers involved in the transitions between the stable ⁴C₁ and ¹C₄ chair forms. In this study, we analyze the influence of the polarizable force field on the ring dynamics of five major types of unsubstituted aldohexoses─glucose, mannose, galactose, altrose, and idose─and their anomers. We simulate microsecond trajectories to capture the influence of the CHARMM36 additive and polarizable carbohydrate force fields on the ring dynamics. The microsecond trajectories allow us to comment on the issues associated with equilibrium molecular dynamics simulations. Further, we use the extended system adaptive biasing force (eABF) method to compare the conformational sampling efficiencies of the additive and polarizable force fields. Our studies reveal that inclusion of polarization enhances the sampling of ring conformations and lowers the energy barriers between the ⁴C₁ and ¹C₄ conformations. Overall, the CHARMM36 additive force field is observed to be rigid and favor the ⁴C₁ conformations. Although the inclusion of polarizability results in enhancing ring flexibility, we observe sampling that does not agree with experimental results, warranting a revision of the polarizable Drude parameters.

Exploration, Representation, and Rationalization of the Conformational Phase Space of N-Glycans

Despite their fundamental biological relevance, structure-property relationships in N-glycans are fundamentally lacking, and their highly multidimensional compositional and conformational phase spaces remain largely unexplored. The torsional flexibility of the glycosidic linkages and the ring dynamics result in wide, rugged free-energy landscapes that are difficult to sample in molecular dynamics simulations. We show that a novel enhanced-sampling scheme combining replica exchange with solute and collective-variable tempering, enabling transitions over all relevant energy barriers, delivers converged distributions of solvated N-glycan conformers. Several dimensionality-reduction algorithms are compared and employed to generate conformational free-energy maps in two dimensions. Together with an originally developed conformation-based nomenclature scheme that uniquely identifies glycan conformers, our modeling procedure is applied to reveal the effect of chemical substitutions on the conformational ensemble of selected high-mannose-type and complex glycans. Moreover, the structure-prediction capabilities of two commonly used glycan force fields are assessed via the theoretical prediction of experimentally available nuclear magnetic resonance J-coupling constants. The results especially confirm the key role of ω and ψ torsion angles in discriminating between different conformational states and suggest an intriguing correlation between the torsional and ring-puckering degrees of freedom that may be biologically relevant.

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.

Conformational preferences of triantennary and tetraantennary hybrid N-glycans in aqueous solution: Insights from 20 μs long atomistic molecular dynamic simulations

In the current study, we have investigated the conformational dynamics of a triantennary (N-glycan1) and tetraantennary (N-glycan2) hybrid N-glycans found on the surface of the HIV glycoprotein using 20 μs long all-atom molecular dynamics (MD) simulations. The main objective of the present study is to elucidate the influence of adding a complex branch on the overall glycan structural dynamics. Our investigation suggests that the average RMSD value increases when a complex branch is added to N-glycan1. However, the RMSD distribution is relatively wider in the case of N-glycan1 compared to N-glycan2, which indicates that multiple complex branches restrict the conformational variability of glycans. A similar observation is obtained from the principal component analysis of both glycans. All the puckering states (⁴C₁ to ¹C₄) of each monosaccharide except mannose are sampled in our simulations, although the ⁴C₁ chair form is energetically more favorable than ¹C₄. In N-glycan1, the 1-6 linkage in the mannose branch [Man(9)-α(1-6)-Man(5)] stays in the gauche-gauche cluster, whereas it moves towards trans-gauche in N-glycan2. For both glycans, mannose branches are more flexible than the complex branches, and adding a complex branch does not influence the dynamics of the mannose branches. We have noticed that the end-to-end distance of the complex branch shortens by ∼ 10 Å in the presence of another complex branch. This suggests that in the presence of an additional complex branch, the other complex branch adopts a close folded structure. All these conformational changes involve the selective formation of inter-residue and water-mediated hydrogen-bond networks.

Novel diagnostic and prognostic factors for the advanced melanoma based on the glycosylation-related changes studied by biophysical profiling methods

Melanoma is a life-threatening disease due to the early onset of metastasis and frequent resistance to the applied treatment. For now, no single histological, immunohistochemical or serological biomarker was able to provide a precise predictive value for the aggressive behavior in melanoma patients. Thus, the search for quantifying methods allowing a simultaneous diagnosis and prognosis of melanoma patients is highly desirable. By investigating specific molecular interactions with some biosensor-based techniques, one can determine novel prognostic factors for this tumor. In our previous study, we have shown the possibility of a qualitative in vitro distinguishing the commercially available melanoma cells at different progression stages based on the measurements of the lectin Concanavalin A interacting with surface glycans present on cells. Here, we present the results of the quantitative diagnostic and prognostic study of both commercial and patient-derived melanoma cells based on the evaluation of two novel factors: lectin affinity and glycan viscoelastic index obtained from the quartz crystal microbalance with dissipation monitoring (QCM-D) measurements. Two approaches to the QCM-D measurements were applied, the first uses the ability of melanoma cells to grow as a monolayer of cells on the sensor (cell-based sensors), and the second shortens the time of the analysis (suspension cell based-sensors). The results were confirmed by the complementary label-free (atomic force microscopy, AFM; and surface plasmon resonance, SPR) and labeling (lectin-ELISA; and microscale thermophoresis, MST) techniques. This new approach provides additional quantitative diagnosis and a personalized prognosis which can be done simultaneously to the traditional histopathological analysis.

2-Dimensional ultra-high performance liquid chromatography and DMT-MM derivatization paired with tandem mass spectrometry for comprehensive serum N-glycome characterization

Glycosylation is a prominent co- and post-translational modification which contributes to a variety of important biological functions. Protein glycosylation characteristics, particularly N-glycosylation, are influenced by changes in one's pathological state, such as through the presence of disease, and as such, there is great interest in N-glycans as potential disease biomarkers. Human serum is an attractive source for N-glycan based biomarker studies as circulatory proteins are representative of one's physiology, with many serum proteins containing N-glycosylation. The difficulty in comprehensively characterizing the serum N-glycome arises from the absence of a biosynthetic template resulting in great structural heterogeneity and complexity. To help overcome these challenges we developed a 2-dimensional liquid chromatography platform which utilizes offline weak anion exchange (WAX) chromatography in the first dimension and hydrophilic interaction liquid chromatography (HILIC) in the second dimension to separate N-glycans by charge, corresponding to degree of sialylation, and size, respectively. Performing these separations offline enables subsequent derivatization with 4-(4,6-Dimethoxy-1,3,5-triazin-2-yl)-4-methylmorpholinium chloride (DMT-MM) for sialic acid linkage determination and the identification of sialic acid linkage isomers. Subsequent tandem mass spectrometry analysis revealed the identification of 212 complete and partial N-glycan structures including low abundant N-glycans containing acetyl and sulphate modifications. The identifications obtained through this platform were then applied to N-glycans released from a set of stage 3 gastric cancer serum samples obtained from patients before (pre-op) and after (post-op) tumour resection to investigate how the serum N-glycome can facilitate differentiation between the two pathological states.

Biomolecular Recognition of the Glycan Neoantigen CA19-9 by Distinct Antibodies

Glycans decorate the cell surface, secreted glycoproteins and glycolipids, and altered glycans are often found in cancers. Despite their high diagnostic and therapeutic potential, however, glycans are polar and flexible molecules that are quite challenging for the development and design of high-affinity binding antibodies. To understand the mechanisms by which glycan neoantigens are specifically recognized by antibodies, we analyze the biomolecular recognition of the tumor-associated carbohydrate antigen CA19-9 by two distinct antibodies using X-ray crystallography. Despite the potential plasticity of glycans and the very different antigen-binding surfaces presented by the antibodies, both structures reveal an essentially identical extended CA19-9 conformer, suggesting that this conformer's stability selects the antibodies. Starting from the bound structure of one of the antibodies, we use the AbLIFT computational algorithm to design a variant with seven core mutations in the variable domain's light-heavy chain interface that exhibits tenfold improved affinity for CA19-9. The results reveal strategies used by antibodies to specifically recognize glycan antigens and show how automated antibody-optimization methods may be used to enhance the clinical potential of existing antibodies.

Glycan Cluster Shielding and Antibody Epitopes on Lassa Virus Envelop Protein

An understanding of how an antiviral monoclonal antibody recognizes its target is vital for the development of neutralizing antibodies and vaccines. The extensive glycosylation of viral proteins almost certainly affects the antibody response, but the investigation of such effects is hampered by the huge range of structures and interactions of surface glycans through their inherent complexity and flexibility. Here, we built an atomistic model of a fully glycosylated envelope protein complex of the Lassa virus and performed molecular dynamics simulations to characterize the impact of surface glycans on the antibody response. The simulations attested to the variety of conformations and interactions of surface glycans. The results show that glycosylation nonuniformly shields the surface of the complex and only marginally affects protein dynamics. The glycans gather in distinct clusters through interaction with protein residues, and only a few regions are left accessible by an antibody. We successfully recovered known protein epitopes by integrating the simulation results with existing sequence- and structure-based epitope prediction methods. The results emphasize the rich structural environment of glycans and demonstrate that shielding is not merely envelopment by a uniform blanket of sugars. This work provides a molecular basis for integrating otherwise elusive structural properties of glycans into vaccine and neutralizing antibody developments.

Identification of Isomeric N-Glycans by Conformer Distribution Fingerprinting using Ion Mobility Mass Spectrometry

Glycans possess unparalleled structural complexity arising from chemically similar monosaccharide building blocks, configurations of anomeric linkages and different branching patterns, potentially giving rise to many isomers. This level of complexity is one of the main reasons that identification of exact glycan structures in biological samples still lags behind that of other biomolecules. Here, we introduce a methodology to identify isomeric N-glycans by determining gas phase conformer distributions (CDs) by measuring arrival time distributions (ATDs) using drift-tube ion mobility spectrometry-mass spectrometry. Key to the approach is the use of a range of well-defined synthetic glycans that made it possible to investigate conformer distributions in the gas phase of isomeric glycans in a systematic manner. In addition, we have computed CD fingerprints by molecular dynamics (MD) simulation, which compared well with experimentally determined CDs. It supports that ATDs resemble conformational populations in the gas phase and offer the prospect that such an approach can contribute to generating a library of CCS distributions (CCSDs) for structure identification.

Multiscale Simulations Examining Glycan Shield Effects on Drug Binding to Influenza Neuraminidase

Influenza neuraminidase is an important drug target. Glycans are present on neuraminidase and are generally considered to inhibit antibody binding via their glycan shield. In this work, we studied the effect of glycans on the binding kinetics of antiviral drugs to the influenza neuraminidase. We created all-atom in silico systems of influenza neuraminidase with experimentally derived glycoprofiles consisting of four systems with different glycan conformations and one system without glycans. Using Brownian dynamics simulations, we observe a two- to eightfold decrease in the rate of ligand binding to the primary binding site of neuraminidase due to the presence of glycans. These glycans are capable of covering much of the surface area of neuraminidase, and the ligand binding inhibition is derived from glycans sterically occluding the primary binding site on a neighboring monomer. Our work also indicates that drugs preferentially bind to the primary binding site (i.e., the active site) over the secondary binding site, and we propose a binding mechanism illustrating this. These results help illuminate the complex interplay between glycans and ligand binding on the influenza membrane protein neuraminidase.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Investigating Conformational Dynamics of Lewis Y Oligosaccharides and Elucidating Blood Group Dependency of Cholera Using Molecular Dynamics

Cholera is caused by Vibrio cholerae and is an example of a blood-group-dependent disease. Recent studies suggest that the receptor-binding B subunit of the cholera toxin (CT) binds histo-blood group antigens at a secondary binding site. Herein, we studied the conformational dynamics of Lewis Y (Le^Y) oligosaccharides, H-tetrasaccharides and A-pentasaccharides, in aqueous solution by conducting accelerated molecular dynamics (aMD) simulations. The flexible nature of both oligosaccharides was displayed in aMD simulations. Furthermore, aMD simulations revealed that for both oligosaccharides in the free form, ⁴C₁ and ¹C₄ puckers were sampled for all but GalNAc monosaccharides, while either the ⁴C₁ (GlcNAc, Gal, GalNAc) or ¹C₄ (Fuc2, Fuc3) pucker was sampled in the CT-bound forms. In aMD, the complete transition from the ⁴C₁ to ¹C₄ pucker was sampled for GlcNAc and Gal in both oligosaccharides. Further, we have observed a transition from the open to closed conformer in the case of A-pentasaccharide, while H-tetrasaccharide remains in the open conformation throughout the simulation. Both oligosaccharides adopted an open conformation in the CT binding site. Moreover, we have investigated the molecular basis of recognition of Le^Y oligosaccharides by the B subunit of the cholera toxin of classical and El Tor biotypes using the molecular mechanics generalized Born surface area (MM/GBSA) scheme. The O blood group determinant, H-tetrasaccharide, exhibits a stronger affinity to both biotypes compared to the A blood group determinant, A-pentasaccharide, which agrees with the experimental data. The difference in binding free energy between O and A blood group determinants mainly arises due to the increased entropic cost and desolvation energy in the case of A-pentasaccharide compared to that of H-tetrasaccharide. Our study also reveals that the terminal Fuc3 contributes most to the binding free energy compared to other carbohydrate residues as it forms multiple hydrogen bonds with CT. Overall, our study might help in designing glycomimetic drugs targeting the cholera toxin.

Structural and Functional Insights into GluK3-kainate Receptor Desensitization and Recovery

GluK3-kainate receptors are atypical members of the iGluR family that reside at both the pre- and postsynapse and play a vital role in the regulation of synaptic transmission. For a better understanding of structural changes that underlie receptor functions, GluK3 receptors were trapped in desensitized and resting/closed states and structures analyzed using single particle cryo-electron microscopy. While the desensitized GluK3 has domain organization as seen earlier for another kainate receptor-GluK2, antagonist bound GluK3 trapped a resting state with only two LBD domains in dimeric arrangement necessary for receptor activation. Using structures as a guide, we show that the N-linked glycans at the interface of GluK3 ATD and LBD likely mediate inter-domain interactions and attune receptor-gating properties. The mutational analysis also identified putative N-glycan interacting residues. Our results provide a molecular framework for understanding gating properties unique to GluK3 and exploring the role of N-linked glycosylation in their modulation.

Replica-Exchange Methods for Biomolecular Simulations

In this study, a replica-exchange method was developed to overcome conformational sampling difficulties in computer simulations of spin glass or other systems with rugged free-energy landscapes. This method was then applied to the protein-folding problem in combination with molecular dynamics (MD) simulation. Owing to its simplicity and sampling efficiency, the replica-exchange method has been applied to many other biological problems and has been continuously improved. The method has often been combined with other sampling techniques, such as umbrella sampling, free-energy perturbation, metadynamics, and Gaussian accelerated MD (GaMD). In this chapter, we first summarize the original replica-exchange molecular dynamics (REMD) method and discuss how new algorithms related to the original method are implemented to add new features. Heterogeneous and flexible structures of an N-glycan in a solution are simulated as an example of applications by REMD, replica exchange with solute tempering, and GaMD. The sampling efficiency of these methods on the N-glycan system and the convergence of the free-energy changes are compared. REMD simulation protocols and trajectory analysis using the GENESIS software are provided to facilitate the practical use of advanced simulation methods.

Molecular Dynamics of Gangliosides

Computational methodologies have immense potential to delineate the dynamic conformations of glycoconjugates including gangliosides, thereby characterizing the conformational adaptability of their glycans upon interacting with various target proteins. Replica-exchange molecular dynamics simulations have been employed to effectively explore the vast conformational spaces of large, branched carbohydrate moieties. When experimentally validated using NMR, molecular simulations can provide dynamical views of molecular recognition events involving the ganglioside glycans.

Glycan Reader is improved to recognize most sugar types and chemical modifications in the Protein Data Bank

Motivation

Glycans play a central role in many essential biological processes. Glycan Reader was originally developed to simplify the reading of Protein Data Bank (PDB) files containing glycans through the automatic detection and annotation of sugars and glycosidic linkages between sugar units and to proteins, all based on atomic coordinates and connectivity information. Carbohydrates can have various chemical modifications at different positions, making their chemical space much diverse. Unfortunately, current PDB files do not provide exact annotations for most carbohydrate derivatives and more than 50% of PDB glycan chains have at least one carbohydrate derivative that could not be correctly recognized by the original Glycan Reader.

Results

Glycan Reader has been improved and now identifies most sugar types and chemical modifications (including various glycolipids) in the PDB, and both PDB and PDBx/mmCIF formats are supported. CHARMM-GUI Glycan Reader is updated to generate the simulation system and input of various glycoconjugates with most sugar types and chemical modifications. It also offers a new functionality to edit the glycan structures through addition/deletion/modification of glycosylation types, sugar types, chemical modifications, glycosidic linkages, and anomeric states. The simulation system and input files can be used for CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Glycan Fragment Database in GlycanStructure.Org is also updated to provide an intuitive glycan sequence search tool for complex glycan structures with various chemical modifications in the PDB.

Availability and implementation

http://www.charmm-gui.org/input/glycan and http://www.glycanstructure.org.

Contact

wonpil@lehigh.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Conformational effects of N-glycan core fucosylation of immunoglobulin G Fc region on its interaction with Fcγ receptor IIIa

Antibody-dependent cellular cytotoxicity (ADCC) is promoted through interaction between the Fc region of immunoglobulin G1 (IgG1) and Fcγ receptor IIIa (FcγRIIIa), depending on N-glycosylation of these glycoproteins. In particular, core fucosylation of IgG1-Fc N-glycans negatively affects this interaction and thereby compromises ADCC activity. To address the mechanisms of this effect, we performed replica-exchange molecular dynamics simulations based on crystallographic analysis of a soluble form of FcγRIIIa (sFcγRIIIa) in complex with IgG1-Fc. Our simulation highlights increased conformational fluctuation of the N-glycan at Asn162 of sFcγRIIIa upon fucosylation of IgG1-Fc, consistent with crystallographic data giving no interpretable electron density for this N-glycan, except for the innermost part. The fucose residue disrupts optimum intermolecular carbohydrate-carbohydrate interactions, rendering this sFcγRIIIa glycan distal from the Fc glycan. Moreover, our simulation demonstrates that core fucosylation of IgG1-Fc affects conformational dynamics and rearrangements of surrounding amino acid residues, typified by Tyr296 of IgG1-Fc, which was more extensively involved in the interaction with sFcγRIIIa without Fc core fucosylation. Our findings offer a structural foundation for designing and developing therapeutic antibodies with improved ADCC activity.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Conformational Heterogeneity of the HIV Envelope Glycan Shield

To better understand the conformational properties of the glycan shield covering the surface of the HIV gp120/gp41 envelope (Env) trimer, and how the glycan shield impacts the accessibility of the underlying protein surface, we performed enhanced sampling molecular dynamics (MD) simulations of a model glycosylated HIV Env protein and related systems. Our simulation studies revealed a conformationally heterogeneous glycan shield with a network of glycan-glycan interactions more extensive than those observed to date. We found that partial preorganization of the glycans potentially favors binding by established broadly neutralizing antibodies; omission of several specific glycans could increase the accessibility of other glycans or regions of the protein surface to antibody or CD4 receptor binding; the number of glycans that can potentially interact with known antibodies is larger than that observed in experimental studies; and specific glycan conformations can maximize or minimize interactions with individual antibodies. More broadly, the enhanced sampling MD simulations described here provide a valuable tool to guide the engineering of specific Env glycoforms for HIV vaccine design.

Conformational Analysis of a High-Mannose-Type Oligosaccharide Displaying Glucosyl Determinant Recognised by Molecular Chaperones Using NMR-Validated Molecular Dynamics Simulation

Exploration of the conformational spaces of flexible oligosaccharides is essential to gain deeper insights into their functional mechanisms. Here we characterised dynamic conformation of a high-mannose-type dodecasaccharide with a terminal glucose residue, a critical determinant recognised by molecular chaperones. The dodecasaccharide was prepared by our developed chemoenzymatic technique, which uses ¹³ C labelling and lanthanide tagging to detect conformation-dependent paramagnetic effects by NMR spectroscopy. The NMR-validated molecular dynamics simulation produced the dynamic conformational ensemble of the dodecasaccharide. This determined its spatial distribution as well as the glycosidic linkage conformation of the terminal glucose determinant. Moreover, comparison of our results with previously reported crystallographic data indicates that the chaperone binding to its target oligosaccharides involves an induced-fit mechanism.

Delineating the conformational flexibility of trisaccharides from NMR spectroscopy experiments and computer simulations

The conformation of saccharides in solution is challenging to characterize in the context of a single well-defined three-dimensional structure. Instead, they are better represented by an ensemble of conformations associated with their structural diversity and flexibility. In this study, we delineate the conformational heterogeneity of five trisaccharides via a combination of experimental and computational techniques. Experimental NMR measurements target conformationally sensitive parameters, including J couplings and effective distances around the glycosidic linkages, while the computational simulations apply the well-calibrated additive CHARMM carbohydrate force field in combination with efficient enhanced sampling molecular dynamics simulation methods. Analysis of conformational heterogeneity is performed based on sampling of discreet states as defined by dihedral angles, on root-mean-square differences of Cartesian coordinates and on the extent of volume sampled. Conformational clustering, based on the glycosidic linkage dihedral angles, shows that accounting for the full range of sampled conformations is required to reproduce the experimental data, emphasizing the utility of the molecular simulations in obtaining an atomic detailed description of the conformational properties of the saccharides. Results show the presence of differential conformational preferences as a function of primary sequence and glycosidic linkage types. Significant differences in conformational ensembles associated with the anomeric configuration of a single glycosidic linkage reinforce the impact of such changes on the conformational properties of carbohydrates. The present structural insights of the studied trisaccharides represent a foundation for understanding the range of conformations adopted in larger oligosaccharides and how these molecules encode their conformational heterogeneity into the monosaccharide sequence.

Conformational sampling of oligosaccharides using Hamiltonian replica exchange with two-dimensional dihedral biasing potentials and the weighted histogram analysis method (WHAM)

Oligosaccharides and polysaccharides exert numerous functional roles in biology through their structural diversity and conformational properties. To investigate their conformational properties using computational methods, Hamiltonian replica exchange (H-REX) combined with two-dimensional grid-based correction maps as biasing potentials (bpCMAP) significantly improves the sampling efficiency about glycosidic linkages. In the current study, we extend the application of H-REX with bpCMAP to complex saccharides and establish systematic procedures for bpCMAP construction, determination of replica distribution, and data analysis. Our main findings are that (1) the bpCMAP for each type of glycosidic linkage can be constructed from the corresponding disaccharide using gas-phase umbrella sampling simulations, (2) the replica distribution can be conveniently determined following the exact definition of the average acceptance ratio based on the assigned distribution of biasing potentials, and (3) the extracted free energy surface (or potential of mean force (PMF)) can be improved using the weighted histogram analysis method (WHAM) allowing for the inclusion of data from the excited state replicas in the calculated probability distribution. The method is applied to a branched N-glycan found on the HIV gp120 protein, and a linear N-glycan. Considering the general importance of N-glycans and the wide appreciation of the sampling problem, the present method represents an efficient procedure for the conformational sampling of complex oligo- and polysaccharides under explicit solvent conditions. More generally, the use of WHAM is anticipated to be of general utility for the calculation of PMFs from H-REX simulations in a wide range of macromolecular systems.

Exploration of conformational spaces of high-mannose-type oligosaccharides by an NMR-validated simulation

Exploration of the conformational spaces of flexible biomacromolecules is essential for quantitatively understanding the energetics of their molecular recognition processes. We employed stable isotope- and lanthanide-assisted NMR approaches in conjunction with replica-exchange molecular dynamics (REMD) simulations to obtain atomic descriptions of the conformational dynamics of high-mannose-type oligosaccharides, which harbor intracellular glycoprotein-fate determinants in their triantennary structures. The experimentally validated REMD simulation provided quantitative views of the dynamic conformational ensembles of the complicated, branched oligosaccharides, and indicated significant expansion of the conformational space upon removal of a terminal mannose residue during the functional glycan-processing pathway.

Enhanced conformational sampling of carbohydrates by Hamiltonian replica-exchange simulation

Knowledge of the structure and conformational flexibility of carbohydrates in an aqueous solvent is important to improving our understanding of how carbohydrates function in biological systems. In this study, we extend a variant of the Hamiltonian replica-exchange molecular dynamics (MD) simulation to improve the conformational sampling of saccharides in an explicit solvent. During the simulations, a biasing potential along the glycosidic-dihedral linkage between the saccharide monomer units in an oligomer is applied at various levels along the replica runs to enable effective transitions between various conformations. One reference replica runs under the control of the original force field. The method was tested on disaccharide structures and further validated on biologically relevant blood group B, Lewis X and Lewis A trisaccharides. The biasing potential-based replica-exchange molecular dynamics (BP-REMD) method provided a significantly improved sampling of relevant conformational states compared with standard continuous MD simulations, with modest computational costs. Thus, the proposed BP-REMD approach adds a new dimension to existing carbohydrate conformational sampling approaches by enhancing conformational sampling in the presence of solvent molecules explicitly at relatively low computational cost.