MM3 potential energy surfaces of trisaccharide models of lambda-, mu-, and nu-carrageenans

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.

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.

Force-induced structural changes in non-sulfated carrageenan based oligosaccharides - a theoretical study

In this work we use the Enforced Geometry Optimization (EGO) approach to simulate force-induced structural changes in the monomer, and di- up to pentameric oligomers of neutral and non-sulfated carrageenan based oligosaccharides. Our results indicate that the monosaccharide unit sequence in the oligomeric structure determines: (i) the type(s) of the enforced conformational transition(s), and (ii) the mechanical resistance to external forces. It is a direct consequence of the different glycosidic bonding types in the examined carregeenan based oligosaccharides.

Density functional conformational study of 2-O-sulfated 3,6 anhydro-α-D-galactose and of neo-κ- and ι-carrabiose molecules in gas phase and water

We examined the conformational preferences of the 2-O-sulfated-3,6-α-D-anhydrogalactose (compound I) and two 1,3 linked disaccharides constituting-κ or ι-carrageenans using density functional and ab initio methods in gas phase and aqueous solution. Systematic modifications of two torsion angles leading to 324 and 144 starting geometries for the compound I and each disaccharide were used to generate adiabatic maps using B3LYP/6-31G(d). The lower energy conformers were then fully optimized using B3LYP, B3PW91 and MP2 with several basis sets. Overall, we discuss the impact of full relaxation on the energy and structure of the dominant conformations, present the performance comparison with previous molecular mechanics calculations if available, and determine whether our results are impacted, when polarization and diffuse functions are added to the 6-31G(d) basis set, or when the MP2 level of theory is used.

NMR study on hydroxy protons of κ- and κ/μ-hybrid carrageenan oligosaccharides: experimental evidence of hydrogen bonding and chemical exchange interactions in κ/μ oligosaccharides

The hydroxy protons of κ- and κ/μ-hybrid carrageenan oligosaccharides have been studied by NMR spectroscopy in 85% H(2)O/15% acetone-d(6). Hydration and hydrogen bonding interactions in di- (κ), tetra- (κκ), hexa (κκκ), and octa- (κκκκ) κ-oligosaccharides and hexa- (κμκ), octa- (κμμκ), and deca- (κμμμκ) κ/μ-oligosaccharides have been investigated by measuring the chemical shifts, temperature coefficients, and chemical exchange of the hydroxy protons. These NMR parameters indicate that no strong and persistent intramolecular hydrogen bonds involving hydroxy protons stabilize the structure of κ-carrageenan oligosaccharides in aqueous solution. In the κ/μ-oligosaccharides, the presence of chemical exchange between OH3 of α-d-Gal-6-sulfate (D6S) and OH2 of β-d-Gal-4-sulfate (G4S) across the β-d-Gal-4-S-(1→4)-α-d-Gal-6-S linkage reveals the existence of a weak hydrogen bond interaction between the two hydroxyl groups. The smaller temperature coefficients of OH2_D6S and OH3_D6S indicate reduced hydration, interpreted as spatial proximity to the 4-sulfate group and O5 ring oxygen of the neighboring G4S residues, respectively. These first experimental results on the conformation of κ/μ-carrageenan oligosaccharides shine light on the potential role of "kinks" in the properties of the three-dimensional carrageenan gel network.

Relaxed energetic maps of kappa-carrabiose: a DFT study

The B3LYP density function was used with the 6-31G(d) basis set to perform relaxed energetic contour maps of the charged form of kappa-carrabiose in the gas phase and for the neutral form first in the gas phase and then by simulating the presence of water as solvent using the Onsager model. Only one starting conformation has been considered to perform all the calculations. Rigid energetic maps have been then constructed either by addition of diffuse or polarization functions to the basis set obtaining in that way 6-31+G(d)//6-31G(d), 6-31+G(d,p)//6-31G(d), and 6-311++G(d,p)//6-31G(d) energetic maps that have been carefully examined. The obtained structures corresponding to the lower energy conformers have been then fully optimized using different basis sets with the B3LYP method, a reversion in term of energy has been observed for the two first minima in the case of the charged disaccharide in the gas phase, this was attributed to the large grid of 30 degrees that could lead to the exclusion of an intermediate value corresponding to the real minimum of energy. We thus suggest that after establishing potential energy maps it is essential to proceed to full optimizations of the lower energy conformers. Calculations using the more accurate correlated method MP2 with the 6-31G(d) basis set have also been performed for conformers of the two disaccharides in the gas phase.