Molecular dynamics simulations of alpha2 --> 8-linked disialoside: conformational analysis and implications for binding to proteins

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, K_LPS (a lipid A-linked form) and K_PG (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 K_LPS makes the acyl chains of LPS more rigid and ordered, while incorporation of K_PG 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 K_LPS is incorporated, whereas it gets larger when K_PG 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.

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.

Complete structural elucidation of an oxidized polysialic acid drug intermediate by nuclear magnetic resonance spectroscopy

Polysialic acid (PSA) is a high molecular weight glycan composed of repeat units of α(2→8) linked 5-N-acetyl-neuraminic acid. Mild periodate oxidation of PSA selectively targets the end sialic acid ring containing three adjacent alcohols generating a putative aldehyde, which can be used for terminal attachment of PSA to therapeutic proteins. The work presented here permitted complete NMR peak assignments of not only the repeat units, but also the two terminal units at each end of oxidized PSA, an intermediate, which can be used to improve drug performance. The assignments were made using a variety of NMR techniques on oligomers of sialic acid as well as oxidized PSA with molecular masses of 4 and 20 kDa. This enabled structure elucidation that showed the actual moiety formed was not the expected aldehyde or its hydrate, but is a hemiacetal between the oxidation site on the terminal sialic acid ring and the penultimate ring. The existence of a hemiacetal structure has major implications on stability, reactivity, and conjugation chemistry of oxidized PSA. The assignment process also revealed deuterium exchange of the axial hydrogen at the 3- (methylene) position of the ring, which was in agreement with the literature.

Theoretical investigation on the binding specificity of sialyldisaccharides with hemagglutinins of influenza A virus by molecular dynamics simulations

Recognition of cell-surface sialyldisaccharides by influenza A hemagglutinin (HA) triggers the infection process of influenza. The changes in glycosidic torsional linkage and the receptor conformations may alter the binding specificity of HAs to the sialylglycans. In this study, 10-ns molecular dynamics simulations were carried out to examine the structural and dynamic behavior of the HAs bound with sialyldisaccharides Neu5Acα(2-3)Gal (N23G) and Neu5Acα(2-6)Gal (N26G). The analysis of the glycosidic torsional angles and the pair interaction energy between the receptor and the interacting residues of the binding site reveal that N23G has two binding modes for H1 and H5 and a single binding mode for H3 and H9. For N26G, H1 and H3 has two binding modes, and H5 and H9 has a single binding mode. The direct and water-mediated hydrogen bonding interactions between the receptors and HAs play dominant roles in the structural stabilization of the complexes. It is concluded from pair interaction energy and Molecular Mechanic-Poisson-Boltzmann Surface Area calculations that N26G is a better receptor for H1 when compared with N23G. N23G is a better receptor for H5 when compared with N26G. However, H3 and H9 can recognize N23G and N26G in equal binding specificity due to the marginal energy difference (≈2.5 kcal/mol). The order of binding specificity of N23G is H3 > H5 > H9 > H1 and N26G is H1 > H3 > H5 > H9, respectively. The proposed conformational models will be helpful in designing inhibitors for influenza virus.

The conformational properties of methyl alpha-(2,8)-di/trisialosides and their N-acyl analogues: implications for anti-Neisseria meningitidis B vaccine design

The conformational properties of di- and trisaccharide fragments of the polysialic acid O-antigen capsular polysaccharide (CPS) of Neisseria meningitidis B (NmB) have been investigated by a combination of solution phase NMR spectroscopy and explicit-solvent molecular dynamics (MD) simulations. Simulations employing 100 ns of conventional MD, as well as 160 ns of replica exchange MD (REMD), with the GLYCAM06 force field were shown to be in agreement with experimental NMR scalar J-coupling and NOE values. The presence of conformational families has been determined by monitoring interglycosidic torsion angles, by comparing structural superimpositions, as well as via a Bayesian statistical analysis of the torsional data. Attempts to augment the immunogenicity of NmB CPS often involve chemical modifications of the N-acetyl moiety. Here the effects of these chemical group modifications on the conformational properties of the trisialoside have been probed via REMD simulations of the N-glycolyl, N-propionyl, N-propyl and N-butanoyl analogues. Although there were conformational families unique to each non-native analogue, the chemical modifications resulted in largely equivalent overall conformational phase-spaces compared to the native trisialoside. On the basis of the conformational distributions, these shared conformational properties suggest that a recurrent global conformational epitope may be present in both the native and chemically modified CPS fragments. Explanations are therefore provided for monoclonal antibody cross-reactivity, in terms of recognition of a shared global CPS conformation, as well as for lack of cross-reactivity, in terms of fine structural differences associated with the N-acyl groups, which may be dominant in highly matured antibody responses.

Siglec-7 undergoes a major conformational change when complexed with the alpha(2,8)-disialylganglioside GT1b

The siglecs are a group of mammalian sialic acid binding receptors expressed predominantly in the immune system. The CD33-related siglecs show complex recognition patterns for sialylated glycans. Siglec-7 shows a preference for alpha(2,8)-disialylated ligands and provides a structural template for studying the key interactions that drive this selectivity. We have co-crystallized Siglec-7 with a synthetic oligosaccharide corresponding to the alpha(2,8)-disialylated ganglioside GT1b. The crystal structure of the complex offers a first glimpse into how this important family of lectins binds the structurally diverse gangliosides. The structure reveals that the C-C' loop, a region implicated in previous studies as driving siglec specificity, undergoes a dramatic conformational shift, allowing it to interact with the underlying neutral glycan core of the ganglioside. The structural data in combination with mutagenesis studies show that binding of the ganglioside is driven by extensive hydrophobic contacts together with key polar interactions and that the binding site structure is complementary to preferred solution conformations of GT1b.

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.

α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.

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.

Conformations of terminal sialyloligosaccharide fragments--a molecular dynamics study

Molecular dynamics simulations have been performed to understand the conformational features of the terminal sialyloligosaccharide fragments NeuNAc alpha(2-3)Gal, NeuNAc alpha(2-6)Gal, NeuNAc alpha(2-8)NeuNAc and NeuNAc alpha(2-9)NeuNAc. The conformational regions A(i), B(i) and C(i) were identified in the Ramachandran plot. Analysis of the 1000 ps trajectories collected through simulation (2000 ps in the case of NeuNAc alpha (2-9)NeuNAc) revealed that these molecules have conformational propensity in region B(i). The occurrence of these molecules in the common conformational space leads to a structural similarity between them. This structural similarity may be an essential requirement for the neuraminidase activity towards sialyloligosaccharides. The local change in the conformation of the active site residues of neuraminidases may contribute for the specificity differences between different linkages of sialyloligosaccharides. A highly conserved water-mediated hydrogen bond observed in these structures between the sugar residues, acts as an additional stabilizing force.

Anomeric configuration, glycosidic linkage, and the solution conformational entropy of O-linked disaccharides

Oligosaccharides perform a large number of biological roles, as dictated by their chemical structure and spatial arrangement. While conformational entropies are usually determined in vacuo by computer modeling, molecular recognition processes normally take place in solution. Here I show results of experiments using size-exclusion chromatography (SEC), an entropically driven solution technique. These clearly differentiate the individual contributions of the alpha and beta anomeric configurations and of the (1 --> 4) and (1 --> 6) glycosidic linkages to the solution conformational entropy of O-linked disaccharides. I also distinguish between the members of the epimeric disaccharide pair isomaltose-melibiose and trace the difference to that between their constituent monosaccharides, alpha-glucose and alpha-galactose.