Carbohydrate NMR chemical shift prediction by GeqShift employing E(3) equivariant graph neural networks.

Carbohydrates, vital components of biological systems, are well-known for their structural diversity. Nuclear Magnetic Resonance (NMR) spectroscopy plays a crucial role in understanding their intricate molecular arrangements and is essential in assessing and verifying the molecular structure of organic molecules. An important part of this process is to predict the NMR chemical shift from the molecular structure. This work introduces a novel approach that leverages E(3) equivariant graph neural networks to predict carbohydrate NMR spectral data. Notably, our model achieves a substantial reduction in mean absolute error, up to threefold, compared to traditional models that rely solely on two-dimensional molecular structure. Even with limited data, the model excels, highlighting its robustness and generalization capabilities. The model is dubbed GeqShift (geometric equivariant shift) and uses equivariant graph self-attention layers to learn about NMR chemical shifts, in particular since stereochemical arrangements in carbohydrate molecules are characteristics of their structures.

Enzymatic synthesis of N-formylated sialosides via a five-enzyme cascade.

Here we report an enzymatic approach to synthesize N-formylneuraminic acid (Neu5Fo) containing sialosides, through a five-enzyme cascade. This method stands as an alternative to traditional chemical syntheses, aiming for precision and efficiency in generating sialosides with a tailored N-formyl group generated directly from formic acid. The newly synthesized Neu5Fo was characterized using various NMR techniques revealing a conformational equilibrium at the amide bond of the formyl group in slow exchange on the NMR time scale with a trans : cis ratio of ∼2 : 1. This work not only suggests potential for exploring the biological roles of sialosides but also points to the possibility of developing novel therapeutic agents.

Nonreducing Sugar Scaffold Enables the Development of Immunomodulatory TLR4-specific LPS Mimetics with Picomolar Potency.

Innate immune defense mechanisms against infection and cancer encompass the modulation of pattern recognition receptor (PRR)-mediated inflammation, including upregulation of various transcription factors and the activation of pro-inflammatory pathways important for immune surveillance. Dysfunction of PRRs-mediated signaling has been implicated in cancer and autoimmune diseases, while the overactivation of PRRs-driven responses during infection can lead to devastating consequences such as acute lung injury or sepsis. We used crystal structure-based design to develop immunomodulatory lipopolysaccharide (LPS) mimetics targeting one of the ubiquitous PRRs, Toll-like Receptor 4 (TLR4). Taking advantage of an exo-anomeric conformation and specific molecular shape of synthetic nonreducing ß,ß-diglucosamine, which was investigated by NMR, we developed two sets of lipid A mimicking glycolipids capable of either potently activating innate immune responses or inhibiting pro-inflammatory signaling. Stereoselective 1,1'-glycosylation towards fully orthogonally protected nonreducing GlcNß(1↔1')ßGlcN followed by stepwise assembly of differently functionalised phosphorylated glycolipids provided biologically active molecules that were evaluated for their ability to trigger or to inhibit cellular innate immune responses. Two LPS mimetics, identified as potent TLR4-specific inducers of the intracellular signaling pathways, serve as vaccine adjuvant- and immunotherapy candidates, while anionic glycolipids with TLR4-inhibitory potential hold therapeutic promise for the management of acute or chronic inflammation.

Structure and mechanism of biosynthesis of Streptococcus mutans cell wall polysaccharide.

Streptococcus mutans, the causative agent of human dental caries, expresses a cell wall attached Serotype c- specific Carbohydrate (SCC) that is critical for cell viability. SCC consists of a repeating →3)α-Rha(1→2)α-Rha(1→ polyrhamnose backbone, with glucose (Glc) side-chains and glycerol phosphate (GroP) decorations. This study reveals that SCC has one major and two minor Glc modifications. The major Glc modification, α-Glc, attached to position 2 of 3-rhamnose, is installed by SccN and SccM glycosyltransferases and is the site of the GroP addition. The minor Glc modifications are ß-Glc linked to position 4 of 3-rhamnose installed by SccP and SccQ glycosyltransferases, and α-Glc attached to position 4 of 2-rhamnose installed by SccN working in tandem with an unknown enzyme. Both the major and the minor ß-Glc modifications control bacterial morphology, but only the GroP and major Glc modifications are critical for biofilm formation.

On the influence of solvent on the stereoselectivity of glycosylation reactions.

Methodology development in carbohydrate chemistry entails the stereoselective formation of C-O bonds as a key step in the synthesis of oligo- and polysaccharides. The anomeric selectivity of a glycosylation reaction is affected by a multitude of parameters, such as the nature of the donor and acceptor, activator/promotor system, temperature and solvent. The influence of different solvents on the stereoselective outcome of glycosylation reactions employing thioglucopyranosides as glycosyl donors with a non-participating protecting group at position 2 has been studied. A large change in selectivity as a function of solvent was observed and a correlation between selectivity and the Kamlet-Taft solvent parameter π* was found. Furthermore, molecular modeling using density functional theory methodology was conducted to decipher the role of the solvent and possible reaction pathways were investigated.

Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy.

Glycans in the form of oligosaccharides, polysaccharides, and glycoconjugates are ubiquitous in nature, and their structures range from linear assemblies to highly branched and decorated constructs. Solution state NMR spectroscopy facilitates elucidation of preferred conformations and shapes of the saccharides, motions, and dynamic aspects related to processes over time as well as the study of transient interactions with proteins. Identification of intermolecular networks at the atomic level of detail in recognition events by carbohydrate-binding proteins known as lectins, unraveling interactions with antibodies, and revealing substrate scope and action of glycosyl transferases employed for synthesis of oligo- and polysaccharides may efficiently be analyzed by NMR spectroscopy. By utilizing NMR active nuclei present in glycans and derivatives thereof, including isotopically enriched compounds, highly detailed information can be obtained by the experiments. Subsequent analysis may be aided by quantum chemical calculations of NMR parameters, machine learning-based methodologies and artificial intelligence. Interpretation of the results from NMR experiments can be complemented by extensive molecular dynamics simulations to obtain three-dimensional dynamic models, thereby clarifying molecular recognition processes involving the glycans.

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures.

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-ß-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

NMR chemical shift prediction and structural elucidation of linker-containing oligo- and polysaccharides using the computer program CASPER.

Carbohydrate structures containing alkyl groups as aglycones are useful for investigating enzyme activity and glycan-protein interactions. Moreover, linker-containing oligosaccharides with a spacer group are commonly used to print glycan microarrays or to prepare protein-conjugates as vaccine candidates. The structural accuracy of these synthesized glycans are essential for interpretation of results from biological experiments in which the compounds have been used and NMR spectroscopy can unravel and confirm their structures. An approach for efficient 1H and 13C NMR chemical shift assignments employed a parallel NOAH-10 measurement followed by NMR spin-simulation to refine the 1H NMR chemical shifts, as exemplified for a disaccharide with an azidoethyl group as an aglycone, the NMR chemical shifts of which have been used to enhance the quality of CASPER (http://www.casper.organ.su.se/casper/). The CASPER program has been further developed to aid characterization of linker-containing oligo- and polysaccharides, either by chemical shift prediction for comparison to experimental NMR data or as structural investigation of synthesized glycans based on acquired unassigned NMR data. The ability of CASPER to elucidate structures of linker-containing oligosaccharides is demonstrated and comparisons to assigned or unassigned NMR data show the utility of CASPER in supporting a proposed oligosaccharide structure. Prediction of NMR chemical shifts of an oligosaccharide, corresponding to the repeating unit of an O-antigen polysaccharide, having a linker as an aglycone and a non-natural substituent derivative thereof are presented to exemplify the diversity of structures handled. Furthermore, NMR chemical shift predictions of synthesized polysaccharides, corresponding to bacterial polysaccharides, containing a linker are described showing that in addition to oligosaccharide structures also polysaccharide structures having an aglycone spacer group can be analyzed by CASPER.

Structural and functional insights into the Pseudomonas aeruginosa glycosyltransferase WaaG and the implications for lipopolysaccharide biosynthesis.

The glycosyltransferase WaaG in Pseudomonas aeruginosa (PaWaaG) is involved in the synthesis of the core region of lipopolysaccharides. It is a promising target for developing adjuvants that could help in the uptake of antibiotics. Herein, we have determined structures of PaWaaG in complex with the nucleotide-sugars UDP-glucose, UDP-galactose, and UDP-GalNAc. Structural comparison with the homolog from Escherichia coli (EcWaaG) revealed five key differences in the sugar-binding pocket. Solution-state NMR analysis showed that WT PaWaaG specifically hydrolyzes UDP-GalNAc and unlike EcWaaG, does not hydrolyze UDP-glucose. Furthermore, we found that a PaWaaG mutant (Y97F/T208R/N282A/T283A/T285I) designed to resemble the EcWaaG sugar binding site, only hydrolyzed UDP-glucose, underscoring the importance of the identified amino acids in substrate specificity. However, neither WT PaWaaG nor the PaWaaG mutant capable of hydrolyzing UDP-glucose was able to complement an E. coli ΔwaaG strain, indicating that more remains to be uncovered about the function of PaWaaG in vivo. This structural and biochemical information will guide future structure-based drug design efforts targeting PaWaaG.

Conformational flexibility of the disaccharide ß-L-Fucp-(1→4)-α-D-Glcp-OMe as deduced from NMR spectroscopy experiments and computer simulations.

Carbohydrates in biological systems are referred to as glycans and modification of their structures is a hallmark indicator of disease. Analysis of the three-dimensional structure forms the basis for further insight into how they function and comparison of crystal structure with solution-state conformation(s) is particularly relevant, which has been performed for the disaccharide ß-L-Fucp-(1→4)-α-D-Glcp-OMe. In water solution the conformational space at the glycosidic linkage between the two sugar residues is identified from molecular dynamics (MD) simulations as having a low-energy exo-syn conformation, deviating somewhat from the solid-state conformation, and two anti-conformational states, i.e., anti-ϕ and anti-ψ, indicating flexibility at the glycosidic linkage. NMR data were obtained from 1D 1H,1H-NOESY and STEP-NOESY experiments, measurement of transglycosidic 3JCH coupling constants and NMR spin-simulation. The free energy profile of the ω torsion angle computed from MD simulation was in excellent agreement with the rotamer distribution from NMR experiment being for gt:gg:tg 38 : 53 : 9, respectively, with a proposed inter-residue O5'⋯HO6 hydrogen bond being predominant in the gg rotamer. Quantum mechanics methodology was used to calculate transglycosidic NMR 3JCH coupling constants, averaged over a conformational ensemble of structures representing various rotamers of exocyclic groups, in good to excellent agreement with Karplus-type relationships previously developed. Furthermore, 1H and 13C NMR chemical shifts were calculated using the same methodology and were found to be in excellent agreement with experimental data.

Alternating L4 loop architecture of the bacterial polysaccharide co-polymerase WzzE.

Lipopolysaccharides such as the enterobacterial common antigen are important components of the enterobacterial cell envelope that act as a protective barrier against the environment and are often polymerized by the inner membrane bound Wzy-dependent pathway. By employing cryo-electron microscopy we show that WzzE, the co-polymerase component of this pathway that is responsible for the length modulation of the enterobacterial common antigen, is octameric with alternating up-down conformations of its L4 loops. The alternating up-down nature of these essential loops, located at the top of the periplasmic bell, are modulated by clashing helical faces between adjacent protomers that flank the L4 loops around the octameric periplasmic bell. This alternating arrangement and a highly negatively charged binding face create a dynamic environment in which the polysaccharide chain is extended, and suggest a ratchet-type mechanism for polysaccharide elongation.

Resolving the complexity in human milk oligosaccharides using pure shift NMR methods and CASPER.

Human milk oligosaccharides belong to an important class of bioactive molecules with diverse effects on the development of infants. NMR is capable of providing vital structural information about oligosaccharides which can aid in determining structure-function relationships. However, this information is often concealed by signal overlap in 1H spectra, due to the narrow chemical shift range and signal multiplicity. Signal overlap in oligosaccharide spectra can be greatly reduced, and resolution improved, by utilising pure shift methods. Here the benefits of combining pure shift methods with the CASPER computational approach to resonance assignment in oligosaccharides are demonstrated.

Ultra-selective, ultra-clean 1D rotating-frame Overhauser effect spectroscopy.

An ultra-selective 1D NMR experiment - GEMSTONE-ROESY - enables clear, unambiguous assignment of ROE signals in the not uncommon situation that traditional selective methods fail. Its usefulness is demonstrated in the analysis of the natural products cyclosporin and lacto-N-difucohexaose I, providing detailed insight into the structures and conformations of these molecules.

Glycan Stability and Flexibility: Thermodynamic and Kinetic Characterization of Nonconventional Hydrogen Bonding in Lewis Antigens.

We provide evidence for CH-based nonconventional hydrogen bonds (H-bonds) for 10 Lewis antigens and two of their rhamnose analogues. We also characterize the thermodynamics and kinetics of the H-bonds in these molecules and present a plausible explanation for the presence of nonconventional H-bonds in Lewis antigens. Using an alternative method to simultaneously fit a series of temperature-dependent fast exchange nuclear magnetic resonance (NMR) spectra, we determined that the H-bonded conformation is favored by ∼1 kcal/mol over the non-H-bonded conformation. Additionally, a comparison of temperature-dependent 13C linewidths in various Lewis antigens and the two rhamnose analogues reveals H-bonds between the carbonyl oxygen of the N-acetyl group of N-acetylglucosamine and the OH2 group of galactose/fucose. The data presented herein provide insight into the contribution of nonconventional H-bonding to molecular structure and could therefore be used for the rational design of therapeutics.

Modeling and Simulation of Bacterial Outer Membranes with Lipopolysaccharides and Capsular Polysaccharides.

Capsule is one of the common virulence factors in Gram-negative bacteria protecting pathogens from host defenses and consists of long-chain capsular polysaccharides (CPS) anchored in the outer membrane (OM). Elucidating structural properties of CPS is important to understand its biological functions as well as the OM properties. However, the outer leaflet of the OM in current simulation studies is represented exclusively by LPS due to the complexity and diversity of CPS. In this work, representative Escherichia coli CPS, KLPS (a lipid A-linked form) and KPG (a phosphatidylglycerol-linked form), are modeled and incorporated into various symmetric bilayers with co-existing LPS in different ratios. All-atom molecular dynamics simulations of these systems have been conducted to characterize various bilayer properties. Incorporation of KLPS makes the acyl chains of LPS more rigid and ordered, while incorporation of KPG makes them less ordered and flexible. These results are consistent with the calculated area per lipid (APL) of LPS, in which the APL of LPS becomes smaller when KLPS is incorporated, whereas it gets larger when KPG is included. Torsional analysis reveals that the influence of the CPS presence on the conformational distributions of the glycosidic linkages of LPS is small, and minor differences are also detected for the inner and outer regions of the CPS. Combined with previously modeled enterobacterial common antigens (ECAs) in the form of mixed bilayers, this work provides more realistic OM models as well as the basis for characterization of interactions between the OM and OM proteins.

Glycosidic α-linked mannopyranose disaccharides: an NMR spectroscopy and molecular dynamics simulation study employing additive and Drude polarizable force fields.

D-Mannose is a structural component in N-linked glycoproteins from viruses and mammals as well as in polysaccharides from fungi and bacteria. Structural components often consist of D-Manp residues joined via α-(1→2)-, α-(1→3)-, α-(1→4)- or α-(1→6)-linkages. As models for these oligo- and polysaccharides, a series of mannose-containing disaccharides have been investigated with respect to conformation and dynamics. Translational diffusion NMR experiments were performed to deduce rotational correlation times for the molecules, 1D 1H,1H-NOESY and 1D 1H,1H-T-ROESY NMR experiments were carried out to obtain inter-residue proton-proton distances and one-dimensional long-range and 2D J-HMBC experiments were acquired to gain information about conformationally dependent heteronuclear coupling constants across glycosidic linkages. To attain further spectroscopic data, the doubly 13C-isotope labeled α-D-[1,2-13C2]Manp-(1→4)-α-D-Manp-OMe was synthesized thereby facilitating conformational analysis based on 13C,13C coupling constants as interpreted by Karplus-type relationships. Molecular dynamics simulations were carried out for the disaccharides with explicit water as solvent using the additive CHARMM36 and Drude polarizable force fields for carbohydrates, where the latter showed broader population distributions. Both simulations sampled conformational space in such a way that inter-glycosidic proton-proton distances were very well described whereas in some cases deviations were observed between calculated conformationally dependent NMR scalar coupling constants and those determined from experiment, with closely similar root-mean-square differences for the two force fields. However, analyses of dipole moments and radial distribution functions with water of the hydroxyl groups indicate differences in the underlying physical forces dictating the wider conformational sampling with the Drude polarizable versus additive C36 force field and indicate the improved utility of the Drude polarizable model in investigating complex carbohydrates.

Primary Structure of Glycans by NMR Spectroscopy.

Glycans, carbohydrate molecules in the realm of biology, are present as biomedically important glycoconjugates and a characteristic aspect is that their structures in many instances are branched. In determining the primary structure of a glycan, the sugar components including the absolute configuration and ring form, anomeric configuration, linkage(s), sequence, and substituents should be elucidated. Solution state NMR spectroscopy offers a unique opportunity to resolve all these aspects at atomic resolution. During the last two decades, advancement of both NMR experiments and spectrometer hardware have made it possible to unravel carbohydrate structure more efficiently. These developments applicable to glycans include, inter alia, NMR experiments that reduce spectral overlap, use selective excitations, record tilted projections of multidimensional spectra, acquire spectra by multiple receivers, utilize polarization by fast-pulsing techniques, concatenate pulse-sequence modules to acquire several spectra in a single measurement, acquire pure shift correlated spectra devoid of scalar couplings, employ stable isotope labeling to efficiently obtain homo- and/or heteronuclear correlations, as well as those that rely on dipolar cross-correlated interactions for sequential information. Refined computer programs for NMR spin simulation and chemical shift prediction aid the structural elucidation of glycans, which are notorious for their limited spectral dispersion. Hardware developments include cryogenically cold probes and dynamic nuclear polarization techniques, both resulting in enhanced sensitivity as well as ultrahigh field NMR spectrometers with a 1H NMR resonance frequency higher than 1 GHz, thus improving resolution of resonances. Taken together, the developments have made and will in the future make it possible to elucidate carbohydrate structure in great detail, thereby forming the basis for understanding of how glycans interact with other molecules.

Elucidation of the O-antigen structure of Escherichia coli O93 and characterization of its biosynthetic genes.

The structure of the O-antigen from the international reference strain Escherichia coli O93:-:H16 has been determined. A nonrandom modal chain-length distribution was observed for the lipopolysaccharide, a pattern which is typical when long O-specific polysaccharides are expressed. By a combination of (i) bioinformatics information on the gene cluster related to O-antigen synthesis including putative function on glycosyl transferases, (ii) the magnitude of NMR coupling constants of anomeric protons, and (iii) unassigned 2D 1H, 13C-HSQC, and 1H,1H-TOCSY NMR spectra it was possible to efficiently elucidate the structure of the carbohydrate polymer in an automated fashion using the computer program CASPER. The polysaccharide also carries O-acetyl groups and their locations were determined by 2D NMR experiments showing that ~½ of the population was 2,6-di-O-acetylated, ~» was 2-O-acetylated, whereas ~» did not carry O-acetyl group(s) in the 3-O-substituted mannosyl residue of the repeating unit. The structure of the tetrasaccharide repeating unit of the O-antigen is given by: →2)-ß-d-Manp-(1→3)-ß-d-Manp2Ac6Ac-(1→4)-ß-d-GlcpA-(1→3)-α-d-GlcpNAc-(1→, which should also be the biological repeating unit and it shares structural elements with capsular polysaccharides from E. coli K84 and K50. The structure of the acidic O-specific polysaccharide from Cellulophaga baltica strain NN015840T differs to that of the O-antigen from E. coli O93 by lacking the O-acetyl group at O6 of the O-acetylated mannosyl residue.

Structural elucidation of the O-antigen polysaccharide from Escherichia coli O125ac and biosynthetic aspects thereof.

Enteropathogenic Escherichia coli O125, the cause of infectious diarrheal disease, is comprised of two serogroups, viz., O125ab and O125ac, which display the aggregative adherence pattern with epithelial cells. Herein, the structure of the O-antigen polysaccharide from E. coli O125ac:H6 has been elucidated. Sugar analysis revealed the presence of fucose, mannose, galactose and N-acetyl-galactosamine as major components. Unassigned 1H and 13C NMR data from one- and two-dimensional NMR experiments of the O125ac O-antigen in conjunction with sugar components were used as input to the CASPER program, which can determine polysaccharide structure in a fully automated way, and resulted in the following branched pentasaccharide structure of the repeating unit: →4)[ß-d-Galp-(1 â†’ 3)]-ß-d-GalpNAc-(1 â†’ 2)-α-d-Manp-(1 â†’ 3)-α-l-Fucp-(1 â†’ 3)-α-d-GalpNAc-(1→, where the side chain is denoted by square brackets. The proposed O-antigen structure was confirmed by 1H and 13C NMR chemical shift assignments and determination of interresidue connectivities. Based on this structure, that of the O125ab O-antigen, which consists of hexasaccharide repeating units with an additional glucosyl group, was possible to establish in a semi-automated fashion by CASPER. The putative existence of gnu and gne in the gene clusters of the O125 serogroups is manifested by N-acetyl-d-galactosamine residues as the initial sugar residue of the biological repeating unit as well as within the repeating unit. The close similarity between O-antigen structures is consistent with the presence of two subgroups in the E. coli O125 serogroup.

Antibiotics and Carbohydrate-Containing Drugs Targeting Bacterial Cell Envelopes: An Overview.

Certain bacteria constitute a threat to humans due to their ability to escape host defenses as they easily develop drug resistance. Bacteria are classified into gram-positive and gram-negative according to the composition of the cell membrane structure. Gram-negative bacteria have an additional outer membrane (OM) that is not present in their gram-positive counterpart; the latter instead hold a thicker peptidoglycan (PG) layer. This review covers the main structural and functional properties of cell wall polysaccharides (CWPs) and PG. Drugs targeting CWPs are discussed, both noncarbohydrate-related (ß-lactams, fosfomycin, and lipopeptides) and carbohydrate-related (glycopeptides and lipoglycopeptides). Bacterial resistance to these drugs continues to evolve, which calls for novel antibacterial approaches to be developed. The use of carbohydrate-based vaccines as a valid strategy to prevent bacterial infections is also addressed.