Modeling ganglioside headgroups by conformational analysis and molecular dynamics

GM1 structural requirements to mediate neuronal functions

Since the 1980s, it has been known that the administration of ganglioside GM1 to cultured cells induced or enhanced neuronal differentiation. GM1 mechanism of action relies on its direct interaction and subsequent activation of the membrane tyrosine kinase receptor, TrkA, which naturally serves as NGF receptor. This process is mediated by the sole oligosaccharide portion of GM1, the pentasaccharide β-Gal-(1-3)-β-GalNAc-(1-4)-[α-Neu5Ac-(2-3)]-β-Gal-(1-4)-β-Glc. Here we detailed the minimum structural requirements of the oligosaccharide portion of GM1 for mediating the TrkA dependent neuritogenic processing. By in vitro and in silico biochemical approaches, we demonstrated that the minimal portion of GM1 required for the TrkA activation is the inner core of the ganglioside's oligosaccharide β-Gal-(1-3)-β-GalNAc-(1-4)-[α-Neu5Ac-(2-3)]-β-Gal. The addition of a sialic acid residue at position 3 of the outer galactose of the GM1 oligosaccharide, which forms the oligosaccharide of GD1a, prevented the interaction with TrkA and the resulting neuritogenesis. On the contrary, the addition of a fucose residue at position 2 of the outer galactose, forming the Fucosyl-GM1 oligosaccharide, did not prevent the TrkA-mediated neuritogenesis.

Structural Analysis of Oligosaccharides and Glycoconjugates Using NMR

Carbohydrate chains play critical roles in cellular recognition and subsequent signal transduction in the nervous system. Furthermore, gangliosides are targets for various amyloidogenic proteins associated with neurodegenerative disorders. To better understand the molecular mechanisms underlying these biological phenomena, atomic views are essential to delineate dynamic biomolecular interactions. Nuclear magnetic resonance (NMR) spectroscopy provides powerful tools for studying structures, dynamics, and interactions of biomolecules at the atomic level. This chapter describes the basics of solution NMR techniques and their applications to the analysis of 3D structures and interactions of glycoconjugates in the nervous system.

Molecular Recognition Insights of Sialic Acid Glycans by Distinct Receptors Unveiled by NMR and Molecular Modeling

All cells are decorated with a highly dense and complex structure of glycan chains, which are mostly attached to proteins and lipids. In this context, sialic acids are a family of nine-carbon acidic monosaccharides typically found at the terminal position of glycan chains, modulating several physiological and pathological processes. Sialic acids have many structural and modulatory roles due to their negative charge and hydrophilicity. In addition, the recognition of sialic acid glycans by mammalian cell lectins, such as siglecs, has been described as an important immunological checkpoint. Furthermore, sialic acid glycans also play a pivotal role in host-pathogen interactions. Various pathogen receptors exposed on the surface of viruses and bacteria are responsible for the binding to sialic acid sugars located on the surface of host cells, becoming a critical point of contact in the infection process. Understanding the molecular mechanism of sialic acid glycans recognition by sialic acid-binding proteins, present on the surface of pathogens or human cells, is essential to realize the biological mechanism of these events and paves the way for the rational development of strategies to modulate sialic acid-protein interactions in diseases. In this perspective, nuclear magnetic resonance (NMR) spectroscopy, assisted with molecular modeling protocols, is a versatile and powerful technique to investigate the structural and dynamic aspects of glycoconjugates and their interactions in solution at the atomic level. NMR provides the corresponding ligand and protein epitopes, essential for designing and developing potential glycan-based therapies. In this review, we critically discuss the current state of knowledge about the structural features behind the molecular recognition of sialic acid glycans by different receptors, naturally present on human cells or pathogens, disclosed by NMR spectroscopy and molecular modeling protocols.

Complement Factor H and Simian Virus 40 bind the GM1 ganglioside in distinct conformations

Mammalian cell surfaces are decorated with a variety of glycan chains that orchestrate development and defense and are exploited by pathogens for cellular attachment and entry. While glycosidic linkages are, in principle, flexible, the conformational space that a given glycan can sample is subject to spatial and electrostatic restrictions imposed by its overall chemical structure. Here, we show how the glycan moiety of the GM1 ganglioside, a branched, monosialylated pentasaccharide that serves as a ligand for various proteins, undergoes differential conformational selection in its interactions with different lectins. Using STD NMR and X-ray crystallography, we found that the innate immune regulator complement Factor H (FH) binds a previously not reported GM1 conformation that is not compatible with the GM1-binding sites of other structurally characterized GM1-binding lectins such as the Simian Virus 40 (SV40) capsid. Molecular dynamics simulations of the free glycan in explicit solvent on the 10 μs timescale reveal that the FH-bound conformation nevertheless corresponds to a minimum in the Gibbs free energy plot. In contrast to the GM1 conformation recognized by SV40, the FH-bound GM1 conformation is associated with poor NOE restraints, explaining how it escaped(1)H-(1)H NOE-restrained modeling in the past and highlighting the necessity for ensemble representations of glycan structures.

Application of paramagnetic NMR-validated molecular dynamics simulation to the analysis of a conformational ensemble of a branched oligosaccharide

Oligosaccharides of biological importance often exhibit branched covalent structures and dynamic conformational multiplicities. Here we report the application of a method that we developed, which combined molecular dynamics (MD) simulations and lanthanide-assisted paramagnetic NMR spectroscopy, to evaluate the dynamic conformational ensemble of a branched oligosaccharide. A lanthanide-chelating tag was attached to the reducing end of the branched tetrasaccharide of GM2 ganglioside to observe pseudocontact shifts as the source of long distance information for validating the conformational ensemble derived from MD simulations. By inspecting the results, the conformational space of the GM2 tetrasaccharide was compared with that of its nonbranched derivative, the GM3 trisaccharide.

NMR characterization of the interactions between lyso-GM1 aqueous micelles and amyloid beta

Gangliosides are targets for a variety of pathologically relevant proteins, including amyloid beta (Abeta), an important component implicated in Alzheimer's disease (AD). To provide a structural basis for this pathogenic interaction associated with AD, we conducted NMR analyses of the Abeta interactions with gangliosides using lyso-GM1 micelles as a model system. Our NMR data revealed that the sugar-lipid interface is primarily perturbed upon binding of Abeta to the micelles, underscoring the importance of the inner part of the ganglioside cluster for accommodating Abeta in comparison with the outer carbohydrate branches that provide microbial toxin- and virus-binding sites.

Structural aspects of ganglioside-containing membranes

The demand for understanding the physical role of gangliosides in membranes is pressing, due to the high number of diverse and crucial biological functions in which they are involved, needing a unifying thread. To this purpose, model systems including gangliosides have been subject of extensive structural studies. Although showing different levels of complication, all models share the need for simplicity, in order to allow for physico-chemical clarity, so they keep far from the extreme complexity of the true biological systems. Nonetheless, as widely agreed, they provide a basic hint on the structural contribution specific molecules can pay to the complex aggregate. This topic we address in the present review. Gangliosides are likely to play their physical role through metamorphism, cooperativity and demixing, that is, they tend to segregate and identify regions where they can dictate and modulate the geometry and the topology of the structure, and its mechanical properties. Strong three-dimensional organisation and cooperativity are exploited to scale up the local arrangement hierarchically from the nano- to the mesoscale, influencing the overall morphology of the structure.

Effect of structural modifications of ganglioside GM2 on intra-molecular carbohydrate-to-carbohydrate interaction and enzymatic susceptibility

The effect of inter-molecular carbohydrate-to-carbohydrate interaction on basic cell biological processes has been well documented and appreciated. In contrast, very little is known about the intra-molecular carbohydrate-to-carbohydrate interaction. The presence of an interaction between the GalNAc and the Neu5Ac in GM2 detected by NMR spectroscopy represents a well-defined intra-molecular carbohydrate-to-carbohydrate interaction. This intriguing interaction is responsible for the GM2-epitope, GalNAcbeta1-->4(Neu5Acalpha2-->3)Gal-, to exhibit a rigid and compact conformation. We hypothesized that this compact conformation may be the cause for both the GalNAc and the Neu5Ac in GM2 to be refractory to enzymatic hydrolysis and the GM2 activator protein is able to interact with the compact trisaccharide GM2-epitope, rendering the GalNAc and the Neu5Ac accessible to beta-hexosaminidase A and sialidase. We have used a series of structurally modified GM2 to study the effect of modifications of sugar chains on the conformation and enzymatic susceptibility of this ganglioside. Our hypothesis was borne out by the fact that when the GalNAcbeta1-->4Gal linkage in GM2 was converted to the GalNAcbeta1-->6Gal, both the GalNAc and the Neu5Ac became susceptible to beta-hexosaminidase A and sialidase, respectively, without GM2 activator protein. We hope our work will engender interest in identifying other intra-molecular carbohydrate-to-carbohydrate interactions in glycoconjugates.

N-glycolyl GM1 ganglioside as a receptor for simian virus 40

Carbohydrate microarrays have emerged as powerful tools in analyses of microbe-host interactions. Using a microarray with 190 sequence-defined oligosaccharides in the form of natural glycolipids and neoglycolipids representative of diverse mammalian glycans, we examined interactions of simian virus 40 (SV40) with potential carbohydrate receptors. While the results confirmed the high specificity of SV40 for the ganglioside GM1, they also revealed that N-glycolyl GM1 ganglioside [GM1(Gc)], which is characteristic of simian species and many other nonhuman mammals, is a better ligand than the N-acetyl analog [GM1(Ac)] found in mammals, including humans. After supplementing glycolipid-deficient GM95 cells with GM1(Ac) and GM1(Gc) gangliosides and the corresponding neoglycolipids with phosphatidylethanolamine lipid groups, it was found that GM1(Gc) analogs conferred better virus binding and infectivity. Moreover, we visualized the interaction of NeuGc with VP1 protein of SV40 by molecular modeling and identified a conformation for GM1(Gc) ganglioside in complex with the virus VP1 pentamer that is compatible with its presentation as a membrane receptor. Our results open the way not only to detailed studies of SV40 infection in relation to receptor expression in host cells but also to the monitoring of changes that may occur with time in receptor usage by the virus.

Synthesis and conformational analysis of an α-d-mannopyranosyl-(1→2)-α-d-mannopyranosyl-(1→6)-α-d-mannopyranose mimic

A mimic of a (1-->2),(1-->6)-mannotrioside was synthesized by replacing the central mannose unit with an enantiomerically pure, conformationally stable trans-diaxial cyclohexanediol. The three-dimensional structure of the molecule was investigated by NMR spectroscopy supported by molecular modelling and was compared to the known features of the natural mannotrioside.

Characterization of the conformational and orientational dynamics of ganglioside GM1 in a dipalmitoylphosphatidylcholine bilayer by molecular dynamics simulations

The structure and dynamics of a single GM1 (Gal5-beta1,3-GalNAc4-beta1,4-(NeuAc3-alpha2,3)-Gal2-beta1,4-Glc1-beta1,1-Cer) embedded in a DPPC bilayer have been studied by MD simulations. Eleven simulations, each of 10 ns productive run, were performed with different initial conformations of GM1. Simulations of GM1-Os in water and of a DPPC bilayer were also performed to delineate the effects of the bilayer and GM1 on the conformational and orientational dynamics of each other. The conformation of the GM1 headgroup observed in the simulations is in agreement with those reported in literature; but the headgroup is restricted when embedded in the bilayer. NeuAc3 is the outermost saccharide towards the water phase. Glc1 and Gal2 prefer a parallel, and NeuAc3, GalNac4 and Gal5 prefer a perpendicular, orientation with respect to the bilayer normal. The overall characteristics of the bilayer are not affected by the presence of GM1; however, GM1 does influence the DPPC molecules in its immediate vicinity. The implications of these observations on the specific recognition and binding of GM1 embedded in a lipid bilayer by exogenous proteins as well as proteins embedded in lipids have been discussed.

A Simple Model System for the Study of Carbohydrate−Aromatic Interactions

A molecular scaffold was identified which enables the establishment of intramolecular interactions between a monosaccharide and a nearby phenyl ring. A group of molecules containing four different monosaccharides (glucose, galactose, N-acetyl-glucosamine, and N-acetyl-galactosamine) was synthesized and used to investigate the extent and nature of this carbohydrate-arene interaction, as well as the effect on the overall 3D structure of the molecules involved. The sugar-aromatic distance was evaluated by rigorous NMR studies supported by molecular modeling and found to be constant throughout the series, independent of the nature of the sugar and of the conformational behavior of the fragment connecting the two elements. Ab initio calculations at the B3LYP/DZV(2d,p) level of theory enable the analysis of the electronic nature of the interaction. The study shows that, given the opportunity, persistent intramolecular aromatic-sugar interactions can be established and can significantly influence overall molecular shape and energetics. These results have important implications in the design of structural mimics of oligosaccharides.

Conformations of higher gangliosides and their binding with cholera toxin - investigation by molecular modeling, molecular mechanics, and molecular dynamics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface higher gangliosides (GT1A and GT1B) and their interaction with Cholera Toxin. The water mediated hydrogen bonding network exists between sugar residues in gangliosides. An integrated molecular modeling, molecular mechanics, and molecular dynamics calculation of cholera toxin complexed with GT1A and GT1B reveal that, the active site of cholera toxin can accommodate these higher gangliosides. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different higher ganglioside towards cholera toxin. This study identifies that the binding site of cholera toxin is shallow and can accommodate a maximum of two NeuNAc residues. The NeuNAc binding site of cholera toxin may be crucial for the design of inhibitors that can prevent the infection of cholera.

Conformational analyses of mycothiol, a critical intracellular glycothiol in Mycobacteria

Intracellular thiols are essential biomolecules, which play several critical roles in living organisms including controlling intracellular redox potential and acting as cofactors for several vital detoxification enzymes including S-transferases and formaldehyde dehydrogenases. The tripeptide gamma-L-glutamyl-L-cysteinylglycine, more commonly known as glutathione, is well known as the major intracellular thiol in eukaryotes and in some bacteria. However, glutathione is absent in the Actinomycetales bacteria such as Mycobacteria and Streptomyces and is believed to be replaced by 1-D-myo-inosityl-2-(N-acetyl-L-cysteinyl)amido-2-deoxy-alpha-D-glucopyranoside, mycothiol, in these organisms. Although much is known about the chemistry and biochemistry of glutathione, currently much less is known concerning mycothiol and its properties. The structure of mycothiol is composed of a glycoside linkage between myo-inositol and D-glucosamine with an N-acetyl-L-cysteine linked to the 2'-amino group of the d-glucosamine moiety. Mycothiol is currently of intense interest due to its essential role in the cellular physiology of Mycobacteria, such as Mycobacterium tuberculosis, and its possible role in antimycobacterial drug resistance. A detailed investigation of its chemistry is therefore essential in ameliorating our knowledge of this key glycothiol, and in shedding additional light on its biochemical role in these pathogenic organisms. This report presents a detailed conformational analysis of mycothiol utilizing a variety of force fields and stochastic search protocols. Cluster analyses of energetically low lying conformations have indicated the presence of several key conformations that are populated in the gas phase and with implicit water solvation. These conformations are compared to recent NMR studies on a derivative of mycothiol. This information should be an important contribution to our basic understanding of the chemistry of this glycothiol and critical in the design of novel inhibitors of pathogen enzymes that require it.

Recognition of tandem sialic acid residues by CD38: a theoretical study

The electronic structures of gangliosides are described using semiempirical and ab inito molecular orbital theories as well as the density functional theory to clarify the causative factors of the differences in inhibitory effects and to elucidate the recognition mechanisms of the enzyme. Our results suggest that CD38 is likely to recognize the two phosphate groups in NAD and the two carboxyl groups in tandem sialic acid residues of gangliosides. The recognition mechanisms of the substrate are proposed based on the good correlation found between the orbital energy of the highest occupied molecular orbital of the gangliosides and the degree of the inhibitory effect.

A new type of anti-ganglioside antibodies present in neurological patients

High titers of anti-GA1 antibodies have been associated with neurological syndromes. In most cases, these antibodies cross-react with the structurally related glycolipids GM1 and GD1b, although specific anti-GA1 antibodies have also been reported. The role of specific anti-GA1 antibodies is uncertain since the presence of GA1 in the human nervous system has not been clarified. A rabbit was immunized with GD1a and its sera were screened for antibody reactivity by standard immunoassay methods (HPTLC-immunostaining and ELISA). Anti-GD1a antibodies were not detected but, unexpectedly, anti-GA1 IgG-antibodies were found. Antibody binding to GA1 was inhibited by soluble GA1 but also by GD1a. These results indicate that the rabbit produced antibodies that recognize epitopes present on the glycolipids, that are absent or not exposed on solid phase adsorbed GD1a. We investigated the presence of these unusual anti-ganglioside antibodies in normal and neurological patient sera. Approximately, 10% of normal human sera contained low titer of specific anti-GA1 IgG-antibodies but none of them recognized soluble GD1a. High titers of IgG-antibodies reacting only with GA1 were detected in 12 patient sera out of 325 analyzed. Of these, 6 sera showed binding that was inhibited by soluble GD1a and four of them also by GM1. This new type of anti-ganglioside antibodies should be considered important elements for understanding of the pathogenesis of these diseases as well as their diagnosis.

Disialogangliosides and Their Interaction with Cholera Toxin—Investigation by Molecular Modeling, Molecular Mechanics and Molecular Dynamics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface disialogangliosides (GD1A, GD1B and GD3) in aqueous environment. The molecular mechanics calculation reveals that water mediated hydrogen bonding network plays a significant role in the structural stabilization of GD1A, GD1B and GD3. These water mediated hydrogen bonds not only exist between neighboring residues but also exist between residues that are separated by 2 to 3 residues in between. The conformational energy difference between different conformational states of gangliosides correlates very well with the number of water mediated and direct hydrogen bonds. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin is worked out. The NeuNAc has a limited allowed eulerian space at the binding site of Cholera Toxin (2.4%). The molecular modeling, molecular mechanics and molecular dynamics of disialoganglioside-cholera toxin complex reveal that cholera toxin can accommodate the disialoganglioside GD1A in three different modes. A single mode of binding is permissible for GD1B and GD3. Direct and water mediated hydrogen bonding interactions stabilizes these binding modes and play an essential role in defining the order of specificity for different disialogangliosides towards cholera toxin. This study not only provides models for the disialoganglioside-cholera toxin complexes but also identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.

Intramolecular Carbohydrate-Aromatic Interactions and Intermolecular van der Waals Interactions Enhance the Molecular Recognition Ability of GM1 Glycomimetics for Cholera Toxin

The design and synthesis of two GM1 glycomimetics, 6 and 7, and analysis of their conformation in the free state and when complexed to cholera toxin is described. These compounds, which include an (R)-cyclohexyllactic acid and an (R)-phenyllactic acid fragment, respectively, display significant affinity for cholera toxin. A detailed NMR spectroscopy study of the toxin/glycomimetic complexes, assisted by molecular modeling techniques, has allowed their interactions with the toxin to be explained at the atomic level. It is shown that intramolecular van der Waals and CH-pi carbohydrate-aromatic interactions define the conformational properties of 7, which adopts a three-dimensional structure significantly preorganized for proper interaction with the toxin. The exploitation of this kind of sugar-aromatic interaction, which is very well described in the context of carbohydrate/protein complexes, may open new avenues for the rational design of sugar mimics.

Conformational analysis and dynamics of mannobiosides and mannotriosides using Monte Carlo/stochastic dynamics simulations

The conformation and dynamics of alpha-(1-->2)-mannobioside, alpha-(1-->6)-mannobioside, and of the trisaccharide alpha-Man-(1-->2)-alpha-Man-(1-->6)-alpha- Man-OMe were studied using Monte Carlo/stochastic dynamics (MC/SD) simulations, the AMBER* force field, and the GB/SA implicit water solvation model. The results are in agreement with available experimental data.

Structural analysis of lipid complexes of GM2-activator protein

The GM2-activator protein (GM2-AP) is a small lysosomal lipid transfer protein essential for the hydrolytic conversion of ganglioside GM2 to GM3 by beta-hexosaminidase A. The crystal structure of human apo-GM2-AP is known to consist of a novel beta-cup fold with a spacious hydrophobic interior. Here, we present two new structures of GM2-AP with bound lipids, showing two different lipid-binding modes within the apolar pocket. The 1.9A structure with GM2 bound shows the position of the ceramide tail and significant conformational differences among the three molecular copies in the asymmetric unit. The tetrasaccharide head group is not visible and is presumed to be disordered. However, its general position could be established through modeling. The structure of a low-pH crystal, determined at 2.5A resolution, has a significantly enlarged hydrophobic channel that merges with the apolar pocket. Electron density inside the pocket and channel suggests the presence of a trapped phospholipid molecule. Structure alignments among the four crystallographically unique monomers provide information on the potential role for lipid binding of flexible chain segments at the rim of the cavity opening. Two discrete orientations of the S130-T133 loop define an open and a closed configuration of the hydrophobic channel that merges with the apolar pocket. We propose: (i) that the low-pH structure represents an active membrane-binding conformation; (ii) that the mobile S130-T133 loop serves as a gate for passage of ligand into the apolar pocket; and (iii) that this loop and the adjacent apolar V59-W63 loop form a surface patch with two exposed tryptophan residues that could interface with lipid bilayers.

Mimicking gangliosides by design: mimics of GM1 headgroup

Detailed knowledge of the three-dimensional structure of ganglioside headgroups has allowed the successful design of structural and functional mimics of ganglioside GMI oligosaccharide. Our recent work in this area is reviewed in this paper.

The origin of anti-GM1 antibodies in neuropathies: the "binding site drift" hypothesis

Elevated titers of serum antibodies against GM1-ganglioside are associated with a variety of autoimmune neuropathies. The origin of these autoantibodies is still unknown, although there is evidence that they are produced by CD5+ B-lymphocytes and that antigen mimicry is involved. Anti-GM, IgM-antibodies in the normal human immunological repertoire are low affinity antibodies that cross-react with other glycoconjugates carrying Gal beta1-3GalNAc and probably do not have GM1-mediated biological activity. Other anti-GM1 IgM-antibodies with higher affinity and/or different fine specificity are present in patients with motor syndromes. Based on our studies of structural requirement for binding, we hypothesize that disease-associated anti-GM1 antibodies originate at random by mutations affecting the binding site of naturally-occurring ones. The hypothesis is conceptually similar to the established phenomenon of "genetic drift" in species evolutionary biology and is therefore termed "binding site drift".