Conformation of the branched O-specific polysaccharide of Shigella dysenteriae type 2: molecular mechanics calculations show a compact helical structure exposing an epitope which potentially mimics galabiose

Substituted Oligosaccharides as Protein Mimics: Deep Learning Free Energy Landscapes

Protein-protein complexes power the majority of cellular processes. Interfering with the formation of such complexes using well-designed mimics is a difficult, yet actively pursued, research endeavor. Due to the limited availability of results on the conformational preferences of oligosaccharides compared to polypeptides, the former have been much less explored than the latter as protein mimics, despite interesting ADMET characteristics. In this work, the conformational landscapes of a series of 956 substituted glucopyranose oligomers of lengths 3 to 12 designed as protein interface mimics are revealed using microsecond-time-scale, enhanced-sampling molecular dynamics simulations. Deep convolutional networks are trained on these large conformational ensembles, to predict the stability of longer oligosaccharide structures from those of their constituent trimer motifs. Deep generative adversarial networks are then designed to suggest plausible conformations for oligosaccharide mimics of arbitrary length and substituent sequences that can subsequently be used as input to docking simulations. Analyzing the performance of the neural networks also yields insights into the intricate collective effects that dominate oligosaccharide conformational dynamics.

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.

NMR and conformational studies of linear and cyclic oligo-(1→6)-β-D-glucosamines

The conformational behavior of a series of linear and cyclic oligo-(1→6)-β-D-glucosamines and their N-acetylated derivatives, which are related to fragments of natural poly-N-acetylglucosamine, was studied by theoretical molecular modeling and experimental determination of transglycosidic vicinal coupling constants (3)J(C,H) and (3)J(H,H). Molecular dynamics simulations were performed under several types of conditions varying in the consideration of ionization of amino groups, solvent effect, and temperature. Neural network clustering and asphericity calculations were performed on the basis of molecular dynamics data. It was shown that disaccharide fragments in the studied linear oligosaccharides were not rigid, and tended to have several conformers, thus determining the overall twisted shape with helical elements. In addition, it was found that the behavior of C5-C6 bond depended significantly upon the simulation conditions. The cyclic di-, tri-, and tetrasaccharides mostly had symmetrical ring-shaped conformations. The larger cycles tended to adopt more complicated shapes, and the conformational behavior of their disaccharide fragments was close to that in the linear oligosaccharides.

Shape: automatic conformation prediction of carbohydrates using a genetic algorithm

Background

Detailed experimental three dimensional structures of carbohydrates are often difficult to acquire. Molecular modelling and computational conformation prediction are therefore commonly used tools for three dimensional structure studies. Modelling procedures generally require significant training and computing resources, which is often impractical for most experimental chemists and biologists. Shape has been developed to improve the availability of modelling in this field.

Results

The Shape software package has been developed for simplicity of use and conformation prediction performance. A trivial user interface coupled to an efficient genetic algorithm conformation search makes it a powerful tool for automated modelling. Carbohydrates up to a few hundred atoms in size can be investigated on common computer hardware. It has been shown to perform well for the prediction of over four hundred bioactive oligosaccharides, as well as compare favourably with previously published studies on carbohydrate conformation prediction.

Conclusion

The Shape fully automated conformation prediction can be used by scientists who lack significant modelling training, and performs well on computing hardware such as laptops and desktops. It can also be deployed on computer clusters for increased capacity. The prediction accuracy under the default settings is good, as it agrees well with experimental data and previously published conformation prediction studies. This software is available both as open source and under commercial licenses.

α2,3/α2,6-Sialylation of N-glycans: non-synonymous signals with marked developmental regulation in bovine reproductive tracts

The glycan part endows cellular glycoconjugates with significant potential for biological recognition. N-Glycan branches often end with alpha2,3/alpha2,6-sialylation, posing the question whether and how placement of the sialic acid at 3 - or 6 -acceptor positions of galactose has cell biological relevance. As attractive model to study developmental regulation we monitored the expression of alpha2,3/alpha2,6-sialylated determinants in fetal and adult bovine testes and ovaries by lectin histochemistry. Distinct expression patterns were detected in both organ types. Oocyte staining, as a prominent example, was restricted to the presence of alpha2,6-sialylated glycans. Treatment with sialidase abolished binding and thus excluded sulfate esters as lectin targets. We added computer simulations to rationalize the observed evidence for non-random expression of the two closely related sialylgalactose isomers. Extensive molecular mechanics and molecular dynamics calculations reveal that the seemingly minor shift of the glycosidic bond from the alpha2,3 position to the alpha2,6 configuration causes significant shape and flexibility changes. They give each disaccharide its own characteristic meaning as signal in the sugar code.

The use of a genetic algorithm search for molecular mechanics (MM3)-based conformational analysis of oligosaccharides

We have implemented a system called glygal that can perform conformational searches on oligosaccharides using several different genetic algorithm (GA) search methods. The searches are performed in the torsion angle conformational space, considering both the primary glycosidic linkages as well as the pendant groups (C-5-C-6 and hydroxyl groups) where energy calculations are performed using the MM3(96) force field. The system includes a graphical user interface for setting calculation parameters and incorporates a 3D molecular viewer. The system was tested using dozens of structures and we present two case studies for two previously investigated O-specific oligosaccharides of the Shigella dysenteriae type 2 and 4. The results obtained using glygal show a significant reduction in the number of structures that need to be sampled in order to find the best conformation, as compared to filtered systematic search.

The conformations of the O-specific polysaccharides of Shigella dysenteriae type 4 and Escherichia coli O159 studied with molecular mechanics (MM3) filtered systematic search

The branched O-antigens of Escherichia coli O159 and Shigella dysenteriae type 4 are structurally related and are known to show cross-reactivity with antibodies. In the present study, conformational analyses were performed on these two O-antigens using molecular mechanics MM3(96) with filtered systematic search. The results show very strong steric restrictions for the trisaccharide at the branch point of the E. coli O159 antigen, especially for the beta-d-GlcNAc-(1-->3)-beta-d-GlcNAc linkage of the main chain. For the type 4 O-antigen the calculations show essentially a single conformation with respect to the alpha-d-GlcNAc-(1-->3)-alpha-d-GlcNAc linkage of the main chain and three different favoured conformations for the fucose branch. Consecutive repeating units of the S. dysenteriae type 4 and E. coli O159 O-antigens form linear extended chains with significant flexibility between the branches. Comparative calculations carried out with the SWEET server indicate that our method of filtered systematic search is a superior method in the case of branched, constrained oligosaccharides. Based on the results of the MM3 calculations, we propose that the common epitope explaining the cross-reactivity comprises the fucose branch, the downstream GlcNAc and part of the uronic acid.

Dynamics of the Escherichia coli O91 O-Antigen Polysaccharide in Solution as Studied by Carbon-13 NMR Relaxation

The dynamics of the O-antigen part of the lipopolysaccharide from the enterohemorrhagic Escherichia coli O91 has been determined in solution using (13)C NMR relaxation measurements at two magnetic field strengths, 9.4 and 14.1 T, thereby facilitating the testing of several dynamical models. The biological repeating unit, consisting of five sugar residues and substituents, could be determined by spectral analysis of different (1)H,(13)C correlations and corroborated by the relaxation data. The site specifically (13)C-labeled material was shown to have approximately 10 repeating units with a narrow distribution. A model-free analysis of the relaxation data revealed a complex dynamical behavior where the sugar residues could be described by a global correlation time (tau(m) = 5.4 ns), generalized order parameters (S(2) approximately 0.63), and different correlation times for internal motions related to their position in the repeating unit along the polymer (tau(e) approximately 360-520 ps). One of the sugar residues showed, in addition, a chemical exchange contribution. Furthermore, a substituent on another sugar residue was described by two order parameters (S(f)(2) = 0.51 and S(s)(2) = 0.21). The solution dynamics of the polysaccharide are thus described by highly intricate motions, both in amplitude and time scales. These results are of significance in the general description of polysaccharides surrounding bacterial cell surfaces and in the presentation of antigenic epitopes to the immune system of an invaded host.