Presentation of membrane-anchored glycosphingolipids determined from molecular dynamics simulations and NMR paramagnetic relaxation rate enhancement

Molecular dynamics simulations and solid-state nuclear magnetic resonance spectroscopy measurements of C-H bond order parameters and effective correlation times in a POPC-GM3 bilayer

Glycolipids such as gangliosides affect the properties of lipid membranes and in extension the interactions between membranes and other biomolecules like proteins. To better understand how the properties of individual lipid molecules can contribute to shape the functional aspects of a membrane, the spatial restriction and dynamics of C-H bond segments can be measured using nuclear magnetic resonance (NMR) spectroscopy. We combine solid-state NMR spectroscopy with all-atom molecular dynamics (MD) simulations to investigate how ganglioside GM3 affects the bilayer structure and dynamics of C-H bond segments. These two methods yield reorientational correlation functions, molecular profiles of C-H bond order parameters |S_CH| and effective correlation times τ_e, which we compare for lipids in POPC bilayers with and without 30 mol% GM3. Our results revealed that all C-H segments of POPC reorient slower in the presence of GM3 and that the defining features of the GM3-POPC bilayer lie in the GM3 headgroup; it gives the bilayer an extended headgroup layer with high order (|S_CH| up to 0.3-0.4) and slow dynamics (τ_e up to 100 ns), a character that may be mechanistically important in ganglioside interactions with other biomolecules.

The importance of side branches of glycosylphosphatidylinositol anchors: a molecular dynamics perspective

Many proteins are anchored to the cell surface of eukaryotes using a unique family of glycolipids called glycosylphosphatidylinositol (GPI) anchors. These glycolipids also exist without a covalently bound protein, in particular on the cell surfaces of protozoan parasites where they are densely populated. GPIs and GPI-anchored proteins participate in multiple cellular processes such as signal transduction, cell adhesion, protein trafficking and pathogenesis of Malaria, Toxoplasmosis, Trypanosomiasis and prion diseases, among others. All GPIs share a common conserved glycan core modified in a cell-dependent manner with additional side glycans or phosphoethanolamine residues. Here, we use atomistic molecular dynamic simulations and perform a systematic study to evaluate the structural properties of GPIs with different side chains inserted in lipid bilayers. Our results show a flop-down orientation of GPIs with respect to the membrane surface and the presentation of the side chain residues to the solvent. This finding agrees well with experiments showing the role of the side residues as active epitopes for recognition of GPIs by macrophages and induction of GPI-glycan-specific immune responses. Protein-GPI interactions were investigated by attaching parasitic GPIs to Green Fluorescent Protein. GPIs are observed to recline on the membrane surface and pull down the attached protein close to the membrane facilitating mutual contacts between protein, GPI and the lipid bilayer. This model is efficient in evaluating the interaction of GPIs and GPI-anchored proteins with membranes and can be extended to study other parasitic GPIs and proteins and develop GPI-based immunoprophylaxis to treat infectious diseases.

Oligosaccharide Presentation Modulates the Molecular Recognition of Glycolipids by Galectins on Membrane Surfaces

Galectins are a family of glycan binding proteins that stand out for the wide range of biological phenomena in which they are involved. Most galectin functions are associated with their glycan binding capacities, which are generally well characterized at the oligosaccharide level, but not at the glycoprotein or glycolipid level. Glycolipids form the part of cell membranes where they can act as galectin cellular receptors. In this scenario, glycan presentation as well as the membrane chemical and structural features are expected to have a strong impact in these molecular association processes. Herein, liposomes were used as membrane mimicking scaffolds for the presentation of glycosphingolipids (GSLs) and to investigate their interaction with Galectin-3 and the N-domain of Galectin-8 (Gal8N). The binding towards GM3 and GM1 and their non-silaylated GSLs was compared to the binding to the free glycans, devoid of lipid. The analysis was carried out using a combination of NMR methods, membrane perturbation studies, and molecular modeling. Our results showed a different tendency of the two galectins in their binding capacities towards the glycans, depending on whether they were free oligosaccharides or as part of GSL inserted into a lipid bilayer, highlighting the significance of GSL glycan presentation on membranes in lectin binding.

How to assess the structural dynamics of transcription factors by integrating sparse NMR and EPR constraints with molecular dynamics simulations

We review recent advances in modeling structural ensembles of transcription factors from nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopic data, integrated with molecular dynamics (MD) simulations. We focus on approaches that confirm computed conformational ensembles by sparse constraints obtained from magnetic resonance. This combination enables the deduction of functional and structural protein models even if nuclear Overhauser effects (NOEs) are too scarce for conventional structure determination. We highlight recent insights into the folding-upon-DNA binding transitions of intrinsically disordered transcription factors that could be assessed using such integrative approaches.

The role of surfaces on amyloid formation

Interfaces can strongly accelerate or inhibit protein aggregation, destabilizing proteins that are stable in solution or, conversely, stabilizing proteins that are aggregation-prone. Although this behaviour is well-known, our understanding of the molecular mechanisms underlying surface-induced protein aggregation is still largely incomplete. A major challenge is represented by the high number of physico-chemical parameters involved, which are highly specific to the considered combination of protein, surface properties, and solution conditions. The key aspect determining the role of interfaces is the relative propensity of the protein to aggregate at the surface with respect to bulk. In this review, we discuss the multiple molecular determinants that regulate this balance. We summarize current experimental techniques aimed at characterizing protein aggregation at interfaces, and highlight the need to complement experimental analysis with theoretical modelling. In particular, we illustrate how chemical kinetic analysis can be combined with experimental methods to provide insights into the molecular mechanisms underlying surface-induced protein aggregation, under both stagnant and agitation conditions. We summarize recent progress in the study of important amyloids systems, focusing on selected relevant interfaces.

Thermal Response in Cellulose Iβ Based on Molecular Dynamics

The structural details of cellulose I β were discussed according to molecular dynamics simulations with the GLYCAM-06 force field. The simulation outcomes were in agreement with previous experimental data, including structural parameters and hydrogen bond pattern at 298 K. We found a new conformation of cellulose Iβ existed at the intermediate temperature that is between the low and high temperatures. Partial chain rotations along the backbone direction were found and conformations of hydroxymethyl groups that alternated from tg to either gt or gg were observed when the temperature increased from 298 K to 400 K. In addition, the gg conformation is preferred than gt. For the structure adopted at high temperature of 500 K, major chains were twisted and two chains detached from each plain. In contrast to the observation under intermediate temperature, the population of hydroxymethyl groups in gt exceeded that in gg conformation at high temperature. In addition, three patterns of hydrogen bonding were identified at low, intermediate and high temperatures in the simulations. The provided structural information indicated the transitions occurred around 350 K and 450 K, considered as the transitional temperatures of cellulose Iβ in this work.

Force Field Comparison of GM1 in a DOPC Bilayer Validated with AFM and FRET Experiments

The great physiological relevance of glycolipids is being increasingly recognized, and glycolipid interactions have been shown to be central to cell-cell recognition, neuronal plasticity, protein-ligand recognition, and other important processes. However, detailed molecular-level understanding of these processes remains to be fully resolved. Molecular dynamics simulations could reveal the details of the glycolipid interactions, but the results may be influenced by the choice of the employed force field. Here, we have compared the behavior and properties of GM1, a common, biologically important glycolipid, using the CHARMM36, OPLS, GROMOS, and Amber99-GLYCAM06 (in bilayers comprising SLIPIDS and LIPID14 lipids) force fields in bilayers comprising 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids and compared the results to atomic force microscopy and fluorescence resonance energy transfer experiments. We found discrepancies within the GM1 behavior displayed between the investigated force fields. Based on a direct comparison with complementary experimental results derived from fluorescence and AFM measurements, we recommend using the Amber99-GLYCAM force field in bilayers comprising LIPID14 or SLIPIDS lipids followed by CHARMM36 and OPLS force fields in simulations. The GROMOS force field is not recommended for reproducing the properties of the GM1 head group.

β-Amyloid aggregation and heterogeneous nucleation

In this article, we consider the role of heterogeneous nucleation in β-amyloid aggregation. Heterogeneous nucleation is more common and occurs at lower levels of supersaturation than homogeneous nucleation. The nucleation period is also the stage at which most of the polymorphism of amyloids arises, this being one of the defining features of amyloids. We focus on several well-known heterogeneous nucleators of β-amyloid, including lipid surfaces, especially those enriched in gangliosides and cholesterol, and divalent metal ions. These two broad classes of nucleators affect β-amyloid particularly in light of the amphiphilicity of these peptides: the N-terminal region, which is largely polar and charged, contains the metal binding site, whereas the C-terminal region is aliphatic and is important in lipid binding. Notably, these two classes of nucleators can interact cooperatively, aggregation begetting greater aggregation.

CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and Lipoglycans

Glycolipids (such as glycoglycerolipids, glycosphingolipids, and glycosylphosphatidylinositol) and lipoglycans (such as lipopolysaccharides (LPS), lipooligosaccharides (LOS), mycobacterial lipoarabinomannan, and mycoplasma lipoglycans) are typically found on the surface of cell membranes and play crucial roles in various cellular functions. Characterizing their structure and dynamics at the molecular level is essential to understand their biological roles, but systematic generation of glycolipid and lipoglycan structures is challenging because of great variations in lipid structures and glycan sequences (i.e., carbohydrate types and their linkages). To facilitate the generation of all-atom glycolipid/LPS/LOS structures, we have developed Glycolipid Modeler and LPS Modeler in CHARMM-GUI ( http://www.charmm-gui.org ), a web-based interface that simplifies building of complex biological simulation systems. In addition, we have incorporated these modules into Membrane Builder so that users can readily build a complex symmetric or asymmetric biological membrane system with various glycolipids and LPS/LOS. These tools are expected to be useful in innovative and novel glycolipid/LPS/LOS modeling and simulation research by easing tedious and intricate steps in modeling complex biological systems and shall provide insight into structures, dynamics, and underlying mechanisms of complex glycolipid-/LPS-/LOS-containing biological membrane systems.

Identification of Rare Lewis Oligosaccharide Conformers in Aqueous Solution Using Enhanced Sampling Molecular Dynamics

Determining the conformations accessible to carbohydrate ligands in aqueous solution is important for understanding their biological action. In this work, we evaluate the conformational free-energy surfaces of Lewis oligosaccharides in explicit aqueous solvent using a multidimensional variant of the swarm-enhanced sampling molecular dynamics (msesMD) method; we compare with multi-microsecond unbiased MD simulations, umbrella sampling, and accelerated MD approaches. For the sialyl Lewis A tetrasaccharide, msesMD simulations in aqueous solution predict conformer landscapes in general agreement with the other biased methods and with triplicate unbiased 10 μs trajectories; these simulations find a predominance of closed conformer and a range of low-occupancy open forms. The msesMD simulations also suggest closed-to-open transitions in the tetrasaccharide are facilitated by changes in ring puckering of its GlcNAc residue away from the ⁴C₁ form, in line with previous work. For sialyl Lewis X tetrasaccharide, msesMD simulations predict a minor population of an open form in solution corresponding to a rare lectin-bound pose observed crystallographically. Overall, from comparison with biased MD calculations, we find that triplicate 10 μs unbiased MD simulations may not be enough to fully sample glycan conformations in aqueous solution. However, the computational efficiency and intuitive approach of the msesMD method suggest potential for its application in glycomics as a tool for analysis of oligosaccharide conformation.

Ganglioside-Mediated Assembly of Amyloid β-Protein: Roles in Alzheimer's Disease

Assembly and deposition of amyloid β-protein (Aβ) is an early and invariable pathological event of Alzheimer's disease (AD), a chronic neurodegenerative disease affecting the neurons in the brain of aging population. Thus, clarification of the molecular mechanism underlying Aβ assembly is crucial not only for understanding the pathogenesis of AD, but also for developing disease-modifying remedies. In 1995, ganglioside-bound Aβ (GAβ), with unique molecular characteristics, including its altered immunoreactivity and its conspicuous ability to accelerate Aβ assembly, was discovered in an autopsied brain showing early pathological changes of AD. Based on these findings, it was hypothesized that GAβ is an endogenous seed for amyloid fibril formation in the AD brain. A body of evidence that supports the GAβ hypothesis has been growing for over 20years as follows. First, the conformational changes of Aβ from a random coil to an α-helix, and then to a β-sheet in the presence of ganglioside were validated by several techniques. Second, the seed activity of GAβ to accelerate the assembly of soluble Aβ into amyloid fibrils was confirmed by various in vitro and in vivo experiments. Third, it was found that the Aβ binding to ganglioside to form GAβ occurs under limited conditions, which were provided by the lipid environment surrounding ganglioside. Fourth, the region-specific Aβ deposition in the brain appeared to be dependent on the presence of the lipid environment that was in favor of GAβ generation. In this chapter, further progress of the study of ganglioside-mediated Aβ assembly, especially from the aspects of physicochemistry, structural biology, and neuropathology, is reviewed.

GENESIS 1.1: A hybrid-parallel molecular dynamics simulator with enhanced sampling algorithms on multiple computational platforms

GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.

Influence of Ganglioside GM1 Concentration on Lipid Clustering and Membrane Properties and Curvature

Gangliosides are a class of glycosphingolipids (GSLs) with amphiphilic character that are found at the outer leaflet of the cell membranes, where their ability to organize into special domains makes them vital cell membrane components. However, a molecular understanding of GSL-rich membranes in terms of their clustered organization, stability, and dynamics is still elusive. To gain molecular insight into the organization and dynamics of GSL-rich membranes, we performed all-atom molecular-dynamics simulations of bicomponent ganglioside GM1 in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) phospholipid bilayers with varying concentrations of GM1 (10%, 20%, and 30%). Overall, the simulations show very good agreement with available experimental data, including x-ray electron density profiles along the membrane normal, NMR carbohydrate proton-proton distances, and x-ray crystal structures. This validates the quality of our model systems for investigating GM1 clustering through an ordered-lipid-cluster analysis. The increase in GM1 concentration induces tighter lipid packing, driven mainly by inter-GM1 carbohydrate-carbohydrate interactions, leading to a greater preference for the positive curvature of GM1-containing membranes and larger cluster sizes of ordered-lipid clusters (with a composite of GM1 and POPC). These clusters tend to segregate and form a large percolated cluster at a 30% GM1 concentration at 293 K. At a higher temperature of 330 K, however, the segregation is not maintained.

Molecular modeling of biomolecules by paramagnetic NMR and computational hybrid methods

The 3D atomic structures of biomolecules and their complexes are key to our understanding of biomolecular function, recognition, and mechanism. However, it is often difficult to obtain structures, particularly for systems that are complex, dynamic, disordered, or exist in environments like cell membranes. In such cases sparse data from a variety of paramagnetic NMR experiments offers one possible source of structural information. These restraints can be incorporated in computer modeling algorithms that can accurately translate the sparse experimental data into full 3D atomic structures. In this review, we discuss various types of paramagnetic NMR/computational hybrid modeling techniques that can be applied to successful modeling of not only the atomic structure of proteins but also their interacting partners. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Histo-Blood Group Antigen Presentation Is Critical for Binding of Norovirus VLP to Glycosphingolipids in Model Membranes

Virus entry depends on biomolecular recognition at the surface of cell membranes. In the case of glycolipid receptors, these events are expected to be influenced by how the glycan epitope close to the membrane is presented to the virus. This presentation of membrane-associated glycans is more restricted than that of glycans in solution, particularly because of orientational constraints imposed on the glycolipid through its lateral interactions with other membrane lipids and proteins. We have developed and employed a total internal reflection fluorescence microscopy-based binding assay and a scheme for molecular dynamics (MD) membrane simulations to investigate the consequences of various glycan presentation effects. The system studied was histo-blood group antigen (HBGA) epitopes of membrane-bound glycosphingolipids (GSLs) derived from small intestinal epithelium of humans (type 1 chain) and dogs (type 2 chain) interacting with GII.4 norovirus-like particles. Our experimental results showed strong binding to all lipid-linked type 1 chain HBGAs but no or only weak binding to the corresponding type 2 chain HBGAs. This is in contrast to results derived from STD experiments with free HBGAs in solution where binding was observed for Lewis x. The MD data suggest that the strong binding to type 1 chain glycolipids was due to the well-exposed (1,2)-linked α-l-Fucp and (1,4)-linked α-l-Fucp residues, while the weaker binding or lack of binding to type 2 chain HBGAs was due to the very restricted accessibility of the (1,3)-linked α-l-Fucp residue when the glycolipid is embedded in a phospholipid membrane. Our results not only contribute to a general understanding of protein-carbohydrate interactions on model membrane surfaces, particularly in the context of virus binding, but also suggest a possible role of human intestinal GSLs as potential receptors for norovirus uptake.

"Rules of Engagement" of Protein-Glycoconjugate Interactions: A Molecular View Achievable by using NMR Spectroscopy and Molecular Modeling

Understanding the dynamics of protein-ligand interactions, which lie at the heart of host-pathogen recognition, represents a crucial step to clarify the molecular determinants implicated in binding events, as well as to optimize the design of new molecules with therapeutic aims. Over the last decade, advances in complementary biophysical and spectroscopic methods permitted us to deeply dissect the fine structural details of biologically relevant molecular recognition processes with high resolution. This Review focuses on the development and use of modern nuclear magnetic resonance (NMR) techniques to dissect binding events. These spectroscopic methods, complementing X-ray crystallography and molecular modeling methodologies, will be taken into account as indispensable tools to provide a complete picture of protein-glycoconjugate binding mechanisms related to biomedicine applications against infectious diseases.

The Roles of Histidines and Charged Residues as Potential Triggers of a Conformational Change in the Fusion Loop of Ebola Virus Glycoprotein

Ebola virus (EBOV) enters cells from late endosomes/lysosomes under mildly acidic conditions. Entry by fusion with the endosomal membrane requires the fusion loop (FL, residues 507–560) of the EBOV surface glycoprotein to undergo a pH-dependent conformational change. To find the pH trigger for this reaction we mutated multiple conserved histidines and charged and uncharged hydrophilic residues in the FL and measured their activity by liposome fusion and cell entry of virus-like particles. The FL location in the membrane was assessed by NMR using soluble and lipid-bound paramagnetic relaxation agents. While we could not identify a single residue to be alone responsible for pH triggering, we propose that a distributed pH effect over multiple residues induces the conformational change that enhances membrane insertion and triggers the fusion activity of the EBOV FL.

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.

Paramagnetic NMR probes for characterization of the dynamic conformations and interactions of oligosaccharides

Paramagnetism-assisted nuclear magnetic resonance (NMR) techniques have recently been applied to a wide variety of biomolecular systems, using sophisticated immobilization methods to attach paramagnetic probes, such as spin labels and lanthanide-chelating groups, at specific sites of the target biomolecules. This is also true in the field of carbohydrate NMR spectroscopy. NMR analysis of oligosaccharides is often precluded by peak overlap resulting from the lack of variability of local chemical structures, by the insufficiency of conformational restraints from nuclear Overhauser effect (NOE) data due to low proton density, and moreover, by the inherently flexible nature of carbohydrate chains. Paramagnetic probes attached to the reducing ends of oligosaccharides cause paramagnetic relaxation enhancements (PREs) and/or pseudocontact shifts (PCSs) resolve the peak overlap problem. These spectral perturbations can be sources of long-range atomic distance information, which complements the local conformational information derived from J couplings and NOEs. Furthermore, paramagnetic NMR approaches, in conjunction with computational methods, have opened up possibilities for the description of dynamic conformational ensembles of oligosaccharides in solution. Several applications of paramagnetic NMR techniques are presented to demonstrate their utility for characterizing the conformational dynamics of oligosaccharides and for probing the carbohydrate-recognition modes of proteins. These techniques can be applied to the characterization of transient, non-stoichiometric interactions and will contribute to the visualization of dynamic biomolecular processes involving sugar chains.

Coupling between overall rotational diffusion and domain motions in proteins and its effect on dielectric spectra

In this work, we formulate a closed-form solution of the model of a semirigid molecule for the case of fluctuating and reorienting molecular electric dipole moment. We illustrate with numeric calculations the impact of protein domain motions on dielectric spectra using the example of the 128 kDa protein dimer of Enzyme I. We demonstrate that the most drastic effect occurs for situations when the characteristic time of protein domain dynamics is comparable to the time of overall molecular rotational diffusion. We suggest that protein domain motions could be a possible explanation for the high-frequency contribution that accompanies the major relaxation dispersion peak in the dielectric spectra of protein aqueous solutions. We propose that the presented computational methodology could be used for the simultaneous analysis of dielectric spectroscopy and nuclear magnetic resonance data. Proteins 2015; 83:1571-1581. © 2015 Wiley Periodicals, Inc.

The Greater Affinity of JC Polyomavirus Capsid for α2,6-Linked Lactoseries Tetrasaccharide c than for Other Sialylated Glycans Is a Major Determinant of Infectivity

The human JC polyomavirus (JCPyV) establishes an asymptomatic, persistent infection in the kidneys of the majority of the population and is the causative agent of the fatal demyelinating disease progressive multifocal leukoencephalopathy (PML) in immunosuppressed individuals. The Mad-1 strain of JCPyV, a brain isolate, was shown earlier to require α2,6-linked sialic acid on the lactoseries tetrasaccharide c (LSTc) glycan for attachment to host cells. In contrast, a JCPyV kidney isolate type 3 strain, WT3, has been reported to interact with sialic acid-containing gangliosides, but the role of these glycans in JCPyV infection has remained unclear. To help rationalize these findings and probe the effects of strain-specific differences on receptor binding, we performed a comprehensive analysis of the glycan receptor specificities of these two representative JCPyV strains using high-resolution X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy, and correlated these data with the results of infectivity assays. We show here that capsid proteins of Mad-1 and WT3 JCPyV can both engage LSTc as well as multiple sialylated gangliosides. However, the binding affinities exhibit subtle differences, with the highest affinity observed for LSTc. Engagement of LSTc is a prerequisite for functional receptor engagement, while the more weakly binding gangliosides are not required for productive infection. Our findings highlight the complexity of virus-carbohydrate interactions and demonstrate that subtle differences in binding affinities, rather than the binding event alone, help determine tissue tropism and viral pathogenesis.

IMPORTANCE

Viral infection is initiated by attachment to receptors on host cells, and this event plays an important role in viral disease. We investigated the receptor-binding properties of human JC polyomavirus (JCPyV), a virus that resides in the kidneys of the majority of the population and can cause the fatal demyelinating disease progressive multifocal leukoencephalopathy (PML) in the brains of immunosuppressed individuals. JCPyV has been reported to interact with multiple carbohydrate receptors, and we sought to clarify how the interactions between JCPyV and cellular carbohydrate receptors influenced infection. Here we demonstrate that JCPyV can engage numerous sialylated carbohydrate receptors. However, the virus displays preferential binding to LSTc, and only LSTc mediates a productive infection. Our findings demonstrate that subtle differences in binding affinity, rather than receptor engagement alone, are a key determinant of viral infection.

Molecular dynamics simulations of membrane- and protein-bound glycolipids using GLYCAM

This chapter outlines protocols for the preparation, execution, and analysis of molecular dynamics (MD) simulations of glycolipids in biologically relevant environments, i.e., imbedded in lipid bilayers or bound to proteins, with the goal of generating biologically relevant structural and dynamic information. Also included is a description of ensemble average (EA) charge set development consistent with the GLYCAM06 force field and its implementation using the AMBER molecular dynamics software suite.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Exploration of conformational spaces of high-mannose-type oligosaccharides by an NMR-validated simulation

Exploration of the conformational spaces of flexible biomacromolecules is essential for quantitatively understanding the energetics of their molecular recognition processes. We employed stable isotope- and lanthanide-assisted NMR approaches in conjunction with replica-exchange molecular dynamics (REMD) simulations to obtain atomic descriptions of the conformational dynamics of high-mannose-type oligosaccharides, which harbor intracellular glycoprotein-fate determinants in their triantennary structures. The experimentally validated REMD simulation provided quantitative views of the dynamic conformational ensembles of the complicated, branched oligosaccharides, and indicated significant expansion of the conformational space upon removal of a terminal mannose residue during the functional glycan-processing pathway.

Structures of B-lymphotropic polyomavirus VP1 in complex with oligosaccharide ligands

B-Lymphotropic Polyomavirus (LPyV) serves as a paradigm of virus receptor binding and tropism, and is the closest relative of the recently discovered Human Polyomavirus 9 (HPyV9). LPyV infection depends on sialic acid on host cells, but the molecular interactions underlying LPyV-receptor binding were unknown. We find by glycan array screening that LPyV specifically recognizes a linear carbohydrate motif that contains α2,3-linked sialic acid. High-resolution crystal structures of the LPyV capsid protein VP1 alone and in complex with the trisaccharide ligands 3′-sialyllactose and 3′-sialyl-N-acetyl-lactosamine (3SL and 3SLN, respectively) show essentially identical interactions. Most contacts are contributed by the sialic acid moiety, which is almost entirely buried in a narrow, preformed cleft at the outer surface of the capsid. The recessed nature of the binding site on VP1 and the nature of the observed glycan interactions differ from those of related polyomaviruses and most other sialic acid-binding viruses, which bind sialic acid in shallow, more exposed grooves. Despite their different modes for recognition, the sialic acid binding sites of LPyV and SV40 are half-conserved, hinting at an evolutionary strategy for diversification of binding sites. Our analysis provides a structural basis for the observed specificity of LPyV for linear glycan motifs terminating in α2,3-linked sialic acid, and links the different tropisms of known LPyV strains to the receptor binding site. It also serves as a useful template for understanding the ligand-binding properties and serological crossreactivity of HPyV9.

Viruses must engage specific receptors on host cells in order to initiate infection. The type of receptor and its concentration on cells determine viral spread and tropism, but for many viruses, the receptor and the mode of recognition by the virus are not known. We have characterized the structural requirements for receptor binding of B-lymphotropic polyomavirus (LPyV). This virus was originally isolated from African Green Monkey lymph node cultures and attracted interest because of its narrow tropism for a human tumor cell line. LPyV is also the closest relative of the recently discovered Human Polyomavirus 9 (HPyV9). We screened the LPyV coat protein VP1 on an carbohydrate microarray and found that it specifically recognizes a linear sugar motif that terminates in α2,3-linked sialic acid. We then determined the structures LPyV VP1 bound to these carbohydrates. The protein has a preformed, deeply recessed binding site for sialic acid. The binding site differs in both architecture and mode of recognition from the binding sites of other viruses. LPyV only binds linear carbohydrates that are able to penetrate into the binding slot.

Application of Metabolic 13C Labeling in Conjunction with High-Field Nuclear Magnetic Resonance Spectroscopy for Comparative Conformational Analysis of High Mannose-Type Oligosaccharides

High mannose-type oligosaccharides are enzymatically trimmed in the endoplasmic reticulum, resulting in various processing intermediates with exposed glycotopes that are recognized by a series of lectins involved in glycoprotein fate determination in cells. Although recent crystallographic data have provided the structural basis for the carbohydrate recognition of intracellular lectins, atomic information of dynamic oligosaccharide conformations is essential for a quantitative understanding of the energetics of carbohydrate–lectin interactions. Carbohydrate NMR spectroscopy is useful for characterizing such conformational dynamics, but often hampered by poor spectral resolution and lack of recombinant techniques required to produce homogeneous glycoforms. To overcome these difficulties, we have recently developed a methodology for the preparation of a homogeneous high mannose-type oligosaccharide with ¹³C labeling using a genetically engineered yeast strain. We herein successfully extended this method to result in the overexpression of ¹³C-labeled Man₉GlcNAc₂ (M9) with a newly engineered yeast strain with the deletion of four genes involved in N-glycan processing. This enabled high-field NMR analyses of ¹³C-labeled M9 in comparison with its processing product lacking the terminal mannose residue ManD2. Long-range NOE data indicated that the outer branches interact with the core in both glycoforms, and such foldback conformations are enhanced upon the removal of ManD2. The observed conformational variabilities might be significantly associated with lectins and glycan-trimming enzymes.

The use of glycoinformatics in glycochemistry

Glycoinformatics is a small but growing branch of bioinformatics and chemoinformatics. Various resources are now available that can be of use to glycobiologists, but also to chemists who work on the synthesis or analysis of carbohydrates. This article gives an overview of existing glyco-specific databases and tools, with a focus on their application to glycochemistry: Databases can provide information on candidate glycan structures for synthesis, or on glyco-enzymes that can be used to synthesize carbohydrates. Statistical analyses of glycan databases help to plan glycan synthesis experiments. 3D-Structural data of protein-carbohydrate complexes are used in targeted drug design, and tools to support glycan structure analysis aid with quality control. Specific problems of glycoinformatics compared to bioinformatics for genomics or proteomics, especially concerning integration and long-term maintenance of the existing glycan databases, are also discussed.

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.

Paramagnetic Cu(II) for probing membrane protein structure and function: inhibition mechanism of the influenza M2 proton channel

Paramagnetic Cu(II) ions enhance nuclear spin relaxation in a distance-dependent fashion and can be used as a structural probe of proteins. Cu(II) can also serve as a functionally important ligand in proteins. Here we investigate the structural basis of Cu(II) inhibition of the influenza M2 proton channel through Cu(II)-induced paramagnetic relaxation enhancement (PRE). (13)C T(1) relaxation rates of the central residues of the transmembrane (TM) domain of M2 are significantly enhanced by Cu(II), and pronounced spectral broadening is observed for the proton-selective residue, His37. These data yielded quantitative distances of (13)C spins to the Cu(II) center and identified the Cu(II) binding site to be Nε2 of His37. This binding site is surrounded by four imidazole rings from the top and four indole rings of Trp41 from the bottom, thus explaining the high affinity of Cu(II) binding. Bound at this location, Cu(II) can inhibit proton currents by perturbing histidine-water proton exchange, preventing histidine conformational dynamics, and interfering with His-Trp cation-π interaction. The Cu(II) binding site is distinct from the binding site of the hydrophobic drug amantadine, which is about 10 Å N-terminal to His37. Consistently, Cu(II) and amantadine induce distinct conformational changes at several key residues, suggesting the possibility of designing new drugs that target the His37 site to inhibit amantadine-resistant mutant M2 proteins. In addition to the high-affinity His37 binding site, we also examined the weaker and nonspecific binding of Cu(II) to membrane-surface lipid phosphates and the extent of the resulting PRE to surface-proximal protein residues. This study demonstrates the feasibility of NMR studies of paramagnetic-ion-complexed membrane proteins, where the ion serves as both a functional ligand and a distance probe.

Formation of GM1 ganglioside clusters on the lipid membrane containing sphingomyeline and cholesterol

GM1 gangliosides form a microdomain with sphingomyeline (SM) and cholesterol (Chol) and are deeply involved in the aggregation of amyloid beta (Aβ) peptides on neural membranes. We performed molecular dynamics simulations on two kinds of lipid bilayers containing GM1 ganglioside: GM1/SM/Chol and GM1/POPC. Both 10 and 100 ns simulations and another set of 10 ns simulations with different initial lipid arrangement essentially showed the same computational results. GM1 molecules in the GM1/SM/Chol membrane were condensed, whereas those in GM1/POPC membrane scattered. That is, the formation of GM1 cluster was observed only on the GM1/SM/Chol mixed membrane. There appeared numerous hydrogen bonds among glycan portions of the GM1 clusters due to the condensation. A comparison in distribution of lipid molecules between the two kinds of membranes suggested that cholesterol had important roles to prevent the membrane from interdigitation and to stabilize other lipids for interacting with each other. This property of cholesterol promotes the formation of GM1 clusters.

Receptor binding by influenza virus: using computational techniques to extend structural data

Influenza attaches to host cells via hemagglutinin binding of cell-surface glycans. These relatively low-affinity interactions involving flexible ligands are critical in determining tissue and host specificity, but their dynamic nature complicates structural characterization of hemagglutinin-receptor complexes. Molecular simulation can assist in analyzing glycan and protein flexibility in crystallized complexes, assessing how binding might change under mutation or altered glycosylation patterns, and evaluating how soluble ligands may relate to physiological presentation on the plasma membrane. Molecular dynamics simulation also has the potential to help integrate structural and dynamic data sources. Here we review recent progress from analysis of molecular dynamics simulation and outline challenges for the future.

Fast NMR data acquisition from bicelles containing a membrane-associated peptide at natural-abundance

In spite of recent technological advances in NMR spectroscopy, its low sensitivity continues to be a major limitation particularly for the structural studies of membrane proteins. The need for a large quantity of a membrane protein and acquisition of NMR data for a long duration are not desirable. Therefore, there is considerable interest in the development of methods to speed up the NMR data acquisition from model membrane samples. In this study, we demonstrate the feasibility of acquiring two-dimensional spectra of an antimicrobial peptide (MSI-78; also known as pexiganan) embedded in isotropic bicelles using natural-abundance (15)N nuclei. A copper-chelated lipid embedded in bicelles is used to speed-up the spin-lattice relaxation of protons without affecting the spectral resolution and thus enabling fast data acquisition. Our results suggest that even a 2D SOFAST-HMQC spectrum can be obtained four times faster using a very small amount (∼3 mM) of a copper-chelated lipid. These results demonstrate that this approach will be useful in the structural studies of membrane-associated peptides and proteins without the need for isotopic enrichment for solution NMR studies.

Optimal mutation sites for PRE data collection and membrane protein structure prediction

Nuclear magnetic resonance paramagnetic relaxation enhancement (PRE) measures long-range distances to isotopically labeled residues, providing useful constraints for protein structure prediction. The method usually requires labor-intensive conjugation of nitroxide labels to multiple locations on the protein, one at a time. Here a computational procedure, based on protein sequence and simple secondary structure models, is presented to facilitate optimal placement of a minimum number of labels needed to determine the correct topology of a helical transmembrane protein. Tests on DsbB (four helices) using just one label lead to correct topology predictions in four of five cases, with the predicted structures <6 Å to the native structure. Benchmark results using simulated PRE data show that we can generally predict the correct topology for five and six to seven helices using two and three labels, respectively, with an average success rate of 76% and structures of similar precision. The results show promise in facilitating experimentally constrained structure prediction of membrane proteins.

Structure and Dynamics of Glycosphingolipids in Lipid Bilayers: Insights from Molecular Dynamics Simulations

Glycolipids are important constituents of biological membranes, and understanding their structure and dynamics in lipid bilayers provides insights into their physiological and pathological roles. Experimental techniques have provided details into their behavior at model and biological membranes; however, computer simulations are needed to gain atomic level insights. This paper summarizes the insights obtained from MD simulations into the conformational and orientational dynamics of glycosphingolipids and their exposure, hydration, and hydrogen-bonding interactions in membrane environment. The organization of glycosphingolipids in raft-like membranes and their modulation of lipid membrane structure are also reviewed.