Comprehensive characterization of oligosaccharide conformational ensembles with conformer classification by free-energy landscape via reproductive kernel Hilbert space

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.

The European Polysaccharide Network of Excellence (EPNOE) research roadmap 2040: Advanced strategies for exploiting the vast potential of polysaccharides as renewable bioresources

Polysaccharides are among the most abundant bioresources on earth and consequently need to play a pivotal role when addressing existential scientific challenges like climate change and the shift from fossil-based to sustainable biobased materials. The Research Roadmap 2040 of the European Polysaccharide Network of Excellence (EPNOE) provides an expert's view on how future research and development strategies need to evolve to fully exploit the vast potential of polysaccharides as renewable bioresources. It is addressed to academic researchers, companies, as well as policymakers and covers five strategic areas that are of great importance in the context of polysaccharide related research: (I) Materials & Engineering, (II) Food & Nutrition, (III) Biomedical Applications, (IV) Chemistry, Biology & Physics, and (V) Skills & Education. Each section summarizes the state of research, identifies challenges that are currently faced, project achievements and developments that are expected in the upcoming 20 years, and finally provides outlines on how future research activities need to evolve.

One-Bond 13C-1H and 13C-13C Spin-Coupling Constants as Constraints in MA'AT Analysis of Saccharide Conformation

MA'AT analysis uses ensembles of redundant experimental NMR spin-coupling constants, parametrized J-coupling equations obtained from density functional theory (DFT) calculations, and circular statistics to produce probability distributions of molecular torsion angles in solution and information on librational motions about these angles (Meredith et al., J. Chem. Info. Model. 2022, 62, 3135-3141). Current DFT methods give nearly quantitative two- and three-bond J_HH, J_CH, and ¹J_CC values for use in MA'AT analysis of saccharides. In contrast, the accuracy of DFT-calculated one-bond ¹J_CH and ¹J_CC values is more difficult to determine, preventing their use in MA'AT modeling. This report describes experimental and computational studies that address this problem using two approaches (Strategies 1 and 2). Differences [¹J_CH^calc - ¹J_CH^exp] (Strategy 1) ranged from -1.2 to 2.5 Hz, giving an average difference of 0.8 ± 1.7 Hz. Percent differences ranged from -0.8% to 1.6%, giving an average % difference of 0.5 ± 1.1%. In comparison, [¹J_CH^MA'AT - ¹J_CH^exp] (Strategy 2) ranged from -1.8 to 0.2 Hz, giving an average difference of -1.2 ± 0.7 Hz. Percent differences ranged from -1.2% to 0.1%, giving an average % difference of -0.8 ± 0.5%. Strategy 1 gave an average difference of 2.1 Hz between calculated and experimental ¹J_CC values, with an average % difference of 5.1 ± 0.2%. Calculated ¹J_CC values were consistently larger than experimental values. Strategy 2 also gave calculated ¹J_CC values that were larger than the experimental values, with an average difference of 2.3 ± 0.6 Hz, and an average % difference of 5.6 ± 1.6%. The findings of both strategies are similar and indicate that ¹J_CH values in saccharides can be calculated nearly quantitatively, but ¹J_CC values appear to be consistently overestimated by ∼5% using current DFT methods.

Between Two Chairs: Combination of Theory and Experiment for the Determination of the Conformational Dynamics of Xylosides

The conformational properties of monosaccharides constitute fundamental features of oligosaccharides. While the energy landscape of monosaccharides can be altered by a specific biochemical environment or by chemical modifications, the analysis of resulting dynamic conformational equilibria is not feasible by experimental means alone. In this work, a series of β-d-xylopyranosides is used to outline how a combination of experimental NMR parameters and computed molecular properties can be used to determine conformers and quantify the composition of conformational equilibria. We demonstrate that identifying the most stable conformers using energy calculations is challenging and computing of NMR shieldings is typically not sensitive enough. On the other hand, computed spin-spin coupling constants for the xyloside ring can be used to unambiguously assign experimental NMR data of dynamic conformational equilibria and quantify the ratio of different conformers in the mixture. As a proof of principle, this procedure allowed to analyze a hitherto unknown dynamic equilibrium of a diamino-xyloside as a precursor of a molecular switch.

Nonconventional NMR Spin-Coupling Constants in Oligosaccharide Conformational Modeling: Structural Dependencies Determined from Density Functional Theory Calculations

Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA'AT modeling of the O-glycosidic linkages of oligosaccharides. Four (¹ J _C1',H1', ¹ J _C1',C2', ² J _C1',H2', and ² J _C2',H1') and eight (¹ J _C4,H4, ¹ J _C3,C4, ¹ J _C4,C5, ² J _C3,H4, ² J _C4,H3, ² J _C5,H4, ² J _C4,H5, and ² J _C3,C5) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C-O torsion angles, ϕ and ψ, respectively. Long-range ⁴ J _H1',H4 was also examined as a potential conformational constraint of either ϕ or ψ. Secondary effects of exocyclic (hydroxyl) C-O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA'AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C-O bonds. The findings also provide new insights into the general effects of exocyclic C-O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.

Experimental and computational characterization of dynamic biomolecular interaction systems involving glycolipid glycans

On cell surfaces, carbohydrate chains that modify proteins and lipids mediate various biological functions, which are exerted not only through carbohydrate-protein interactions but also through carbohydrate-carbohydrate interactions. These glycans exhibit considerable degrees of conformational variability and often form clusters providing multiple binding sites. The integration of nuclear magnetic resonance spectroscopy and molecular dynamics simulation has made it possible to delineate the dynamical structures of carbohydrate chains. This approach has facilitated the remodeling of oligosaccharide conformational space in the prebound state to improve protein-binding affinity and has been applied to visualize dynamic carbohydrate-carbohydrate interactions that control glycoprotein-glycoprotein complex formation. Functional glycoclusters have been characterized by experimental and computational approaches applied to various model membranes and artificial self-assembling systems. This line of investigation has provided dynamic views of molecular assembling on glycoclusters, giving mechanistic insights into physiological and pathological molecular events on cell surfaces as well as clues for the design and creation of molecular systems exerting improved glycofunctions. Further development and accumulation of such studies will allow detailed understanding and artificial control of the "glycosynapse" foreseen by Dr. Sen-itiroh Hakomori.

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.