Conformation of a trimannoside bound to mannose-binding protein by nuclear magnetic resonance and molecular dynamics simulations

Computational Chemistry Tools in Glycobiology: Modelling of Carbohydrate–Protein Interactions

Molecular modelling provides a major impact in the field of glycosciences, helping in the characterisation of the molecular basis of the recognition between lectins from pathogens and human glycoconjugates, and in the design of glycocompounds with anti-infectious properties. The conformational properties of oligosaccharides are complex, and therefore, the simulation of these properties is a challenging task. Indeed, the development of suitable force fields is required for the proper simulation of important problems in glycobiology, such as the interatomic interactions responsible for oligosaccharide and glycoprotein dynamics, including O-linkages in oligo- and polysaccharides, and N- and O-linkages in glycoproteins. The computational description of representative examples is discussed, herein, related to biologically active oligosaccharides and their interaction with lectins and other proteins, and the new routes open for the design of glycocompounds with promising biological activities.

Molecular mechanisms of inhibition of influenza by surfactant protein D revealed by large-scale molecular dynamics simulation

Surfactant protein D (SP-D), a mammalian C-type lectin, is the primary innate inhibitor of influenza A virus (IAV) in the lung. Interactions of SP-D with highly branched viral N-linked glycans on hemagglutinin (HA), an abundant IAV envelope protein and critical virulence factor, promote viral aggregation and neutralization through as yet unknown molecular mechanisms. Two truncated human SP-D forms, wild-type (WT) and double mutant D325A+R343V, representing neck and carbohydrate recognition domains are compared in this study. Whereas both WT and D325A+R343V bind to isolated glycosylated HA, WT does not inhibit IAV in neutralization assays; in contrast, D325A+R343V neutralization compares well with that of full-length native SP-D. To elucidate the mechanism for these biochemical observations, we have determined crystal structures of D325A+R343V in the presence and absence of a viral nonamannoside (Man9). On the basis of the D325A+R343V-Man9 structure and other crystallographic data, models of complexes between HA and WT or D325A+R343V were produced and subjected to molecular dynamics. Simulations reveal that whereas WT and D325A+R343V both block the sialic acid receptor site of HA, the D325A+R343V complex is more stable, with stronger binding caused by additional hydrogen bonds and hydrophobic interactions with HA residues. Furthermore, the blocking mechanism of HA differs for WT and D325A+R343V because of alternate glycan binding modes. The combined results suggest a mechanism through which the mode of SP-D-HA interaction could significantly influence viral aggregation and neutralization. These studies provide the first atomic-level molecular view of an innate host defense lectin inhibiting its viral glycoprotein target.

Shaping up for structural glycomics: a predictive protocol for oligosaccharide conformational analysis applied to N-linked glycans

•

Aqueous 10 μs simulations of N-linked mannosyl cores and sialyl Lewis (sLe) antennae are validated.

•

Sequence dependent glycosidic linkage and pyranose ring μs motions are implicated in bioactivity.

•

Stacked pyranoses in sLe^a and sLe^x are predicted to be atypically rigid on μs timescales.

•

In a 25 μs simulation of sLe^x, all known conformers were sampled within the initial 10 μs of dynamics.

•

Unbiased 10 μs simulations are proposed as a route to systematic and accurate glycomic 3D-analysis.

The human glycome comprises a vast untapped repository of 3D-structural information that holds the key to glycan recognition and a new era of rationally designed mimetic chemical probes, drugs, and biomaterials. Toward routine prediction of oligosaccharide conformational populations and exchange rates at thermodynamic equilibrium, we apply hardware-accelerated aqueous molecular dynamics to model μs motions in N-glycans that underpin inflammation and immunity. In 10 μs simulations, conformational equilibria of mannosyl cores, sialyl Lewis (sLe) antennae, and constituent sub-sequences agreed with prior refinements (X-ray and NMR). Glycosidic linkage and pyranose ring flexing were affected by branching, linkage position, and secondary structure, implicating sequence dependent motions in glycomic functional diversity. Linkage and ring conformational transitions that have eluded precise quantification by experiment and conventional (ns) simulations were predicted to occur on μs timescales. All rings populated non-chair shapes and the stacked galactose and fucose pyranoses of sLe^a and sLe^x were rigidified, suggesting an exploitable 3D-signature of cell adhesion protein binding. Analyses of sLe^x dynamics over 25 μs revealed that only 10 μs were sufficient to explore all aqueous conformers. This simulation protocol, which yields conformational ensembles that are independent of initial 3D-structure, is proposed as a route to understanding oligosaccharide recognition and structure–activity relationships, toward development of carbohydrate-based novel chemical entities.

Molecular modeling and NMR studies of benzyl substituted mannosyl trisaccharide binding to two mannose-specific lectins: Allium sativam agglutinin I and Concanavalin A

The interaction of trimannoside, α-benzyl 3, 6-di-O-(α-D-mannopyranosyl)-α-D-mannopyranoside, 1 with ASAI (Allium sativam agglutinin I, garlic lectin) was studied to reveal the conformational preferences of this ligand in bound-state and detailed binding mode at atomic level. The binding phenomenon was then compared with another well-known mannose-binding lectin, ConA (Concanavalin A). Structural studies of the ligand in free state were done using NMR spectroscopy and Molecular Dynamics simulations. It is found that the substituted-trimannoside can undergo conformational transitions in solution, with one major and one minor conformation per glycosidic linkage (α 1→3 and α 1→6). On the other hand in the bound-state only one of the two major conformations was significantly populated. The role of phenyl ring in the binding process was explored. An extended binding site was observed for the trimannoside in ASAI utilizing the aromatic substituent, which is not seen in ConA. Binding data from difference absorption spectroscopy supported this fact that the binding of benzyl-substituted ligand is tighter with ASAI than ConA. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 952-967, 2010.

Use of residual dipolar couplings in structural analysis of protein-ligand complexes by solution NMR spectroscopy

Investigation of structure-function relationships in protein complexes, specifically protein-ligand interactions, carry great significance in elucidating the structural and mechanistic bases of molecular recognition events and their role in regulating cell processes. Nuclear magnetic resonance (NMR) spectroscopy is one of the leading structural and analytical techniques in in-depth studies of protein-ligand interactions. Recent advances in NMR methodology such as transverse relaxation-optimized spectroscopy (TROSY) and residual dipolar couplings (RDCs) measured in liquid crystalline alignment medium, offer a viable alternative to traditional nuclear Overhauser enhancement (NOE)-based approaches for structure determination of large protein complexes. RDCs provide a way to constrain the relative orientation of two molecules in complex with each other by aligning their independently determined order tensors. The potential for utilization of RDCs can be extended to proteins with multiple ligands or even multimeric protein-ligand complexes, where symmetry properties of the protein can be taken advantage of. Availability of effective RDC data collection and analysis protocols can certainly aid this process by their incorporation into structure calculation protocols using intramolecular and intermolecular orientational restraints. This chapter discusses in detail some of these protocols including methods for sample preparation in liquid crystalline media, NMR experiments for RDC data collection, as well as software tools for RDC data analysis and protein-ligand complex structure determination.

Docking, synthesis, and NMR studies of mannosyl trisaccharide ligands for DC-SIGN lectin

DC-SIGN, a lectin, which presents at the surface of immature dendritic cells, constitutes nowadays a promising target for the design of new antiviral drugs. This lectin recognizes highly glycosylated proteins present at the surface of several pathogens such as HIV, Ebola virus, Candida albicans, Mycobacterium tuberculosis, etc. Understanding the binding mode of this lectin is a topic of tremendous interest and will permit a rational design of new and more selective ligands. Here, we present computational and experimental tools to study the interaction of di- and trisaccharides with DC-SIGN. Docking analysis of complexes involving mannosyl di- and trisaccharides and the carbohydrate recognition domain (CRD) of DC-SIGN have been performed. Trisaccharides Manalpha1,2[Manalpha1,6]Man 1 and Manalpha1,3[Manalpha1,6]Man 2 were synthesized from an orthogonally protected mannose as a common intermediate. Using these ligands and the soluble extracellular domain (ECD) of DC-SIGN, NMR experiments based on STD and transfer-NOE were performed providing additional information. Conformational analysis of the mannosyl ligands in the free and bound states was done. These studies have demonstrated that terminal mannoses at positions 2 or 3 in the trisaccharides are the most important moiety and present the strongest contact with the binding site of the lectin. Multiple binding modes could be proposed and therefore should be considered in the design of new ligands.

Structural glycobiology: a game of snakes and ladders

Oligo- and polysaccharides are infamous for being extremely flexible molecules, populating a series of well-defined rotational isomeric states under physiological conditions. Characterization of this heterogeneous conformational ensemble has been a major obstacle impeding high-resolution structure determination of carbohydrates and acting as a bottleneck in the effort to understand the relationship between the carbohydrate structure and function. This challenge has compelled the field to develop and apply theoretical and experimental methods that can explore conformational ensembles by both capturing and deconvoluting the structural and dynamic properties of carbohydrates. This review focuses on computational approaches that have been successfully used in combination with experiment to detail the three-dimensional structure of carbohydrates in a solution and in a complex with proteins. In addition, emerging experimental techniques for three-dimensional structural characterization of carbohydrate-protein complexes and future challenges in the field of structural glycobiology are discussed. The review is divided into five sections: (1) The complexity and plasticity of carbohydrates, (2) Predicting carbohydrate-protein interactions, (3) Calculating relative and absolute binding free energies for carbohydrate-protein complexes, (4) Emerging and evolving techniques for experimental characterization of carbohydrate-protein structures, and (5) Current challenges in structural glycoscience.

GLYCAM06: a generalizable biomolecular force field. Carbohydrates

A new derivation of the GLYCAM06 force field, which removes its previous specificity for carbohydrates, and its dependency on the AMBER force field and parameters, is presented. All pertinent force field terms have been explicitly specified and so no default or generic parameters are employed. The new GLYCAM is no longer limited to any particular class of biomolecules, but is extendible to all molecular classes in the spirit of a small-molecule force field. The torsion terms in the present work were all derived from quantum mechanical data from a collection of minimal molecular fragments and related small molecules. For carbohydrates, there is now a single parameter set applicable to both alpha- and beta-anomers and to all monosaccharide ring sizes and conformations. We demonstrate that deriving dihedral parameters by fitting to QM data for internal rotational energy curves for representative small molecules generally leads to correct rotamer populations in molecular dynamics simulations, and that this approach removes the need for phase corrections in the dihedral terms. However, we note that there are cases where this approach is inadequate. Reported here are the basic components of the new force field as well as an illustration of its extension to carbohydrates. In addition to reproducing the gas-phase properties of an array of small test molecules, condensed-phase simulations employing GLYCAM06 are shown to reproduce rotamer populations for key small molecules and representative biopolymer building blocks in explicit water, as well as crystalline lattice properties, such as unit cell dimensions, and vibrational frequencies.

Density functional study of the conformational space of 4C1 D-glucuronic acid

The conformational space of (4)C(1) alpha- and beta-d-glucuronic acid was scanned by HF/3-21G(p) calculations followed by optimization of the 15 most stable structures for each, using the B3LYP density functional theory method in conjunction with a diffuse polarized valence triple-zeta basis set. We found a general preference of the alpha anomers in the isolated molecules in agreement with the large endo-anomeric hyperconjugation effects in these structures. From the other intramolecular interactions (exo-anomeric hyperconjugation, hydrogen-bonding, dipole-dipole, and steric interactions), the effect of the hydrogen bonding is the most pronounced and plays a major role in determining the stability order within the alpha and beta series. The most stable conformer of both alpha and beta (4)C(1) d-glucuronic acid is the structure with the maximum number (5) of intramolecular hydrogen bonds. Introduction of solvent (water) effects by the SCI-PCM model resulted in two characteristic changes of the energetic properties: the gas-phase stability order changed considerably, and the energy range of the 15 most stable conformers decreased from 30 to 15 kJ/mol. The geometrical parameters reflect well the superimposed effects of hyperconjugation and hydrogen-bonding interactions. Most characteristics are the variations of the C-O bond distances (within a range of 0.04 A) upon the combined intramolecular effects.

Conformation of 3'CMP bound to RNase A using TrNOESY

The conditions for accurately determining distance constraints from TrNOESY data on a small ligand (3'CMP) bound to a small protein (RNase A, <14 kDa) are described. For small proteins, normal TrNOESY conditions of 10:1 ligand:protein or greater can lead to inaccurate structures for the ligand-bound conformation due to the contribution of the free ligand to the TrNOESY signals. By using two ligand:protein ratios (2:1 and 5:1), which give the same distance constraints, a conformation of 3'CMP bound to RNase A was determined (glycosidic torsion angle, chi=-166 degrees ; pseudorotational phase angle, 0 degrees < or = P < or =36 degrees ). Ligand-protein NOESY cross peaks were also observed and used to dock 3'CMP into the binding pocket of the apo-protein (7rsa). After energy minimization, the conformation of the 3'CMP:RNase A complex was similar to the X-ray structure (1 rpf) except that a C3'-endo conformation for the ribose ring (rather than C2'-exo conformation) was found in the TrNOESY structure.

High-resolution four-dimensional carbon-correlated 1H-1H ROESY experiments employing isotags and the filter diagonalization method for effective assignment of glycosidic linkages in oligosaccharides

Four-dimensional nuclear magnetic resonance spectroscopy of oligosaccharides that correlates 1H-1H ROESY cross peaks to two additional 13C frequency dimensions is reported. The 13C frequencies were introduced by derivatization of all free hydroxyl groups with doubly 13C-labeled acetyl isotags. Pulse sequences were optimized for processing with the filter diagonalization method. The extensive overlap typically observed in 2D ROESY 1H-1H planes was alleviated by resolution of ROESY cross peaks in the two added dimensions associated with the carbon frequencies of the isotags. This enabled the interresidue 1H-1H ROESY cross peaks to be unambiguously assigned hence spatially proximate sugar spin systems across glycosidic bonds could be effectively ascertained. An experiment that selectively amplifies interresidue ROESY 1H-1H cross peaks is also reported. It moves the magnetization of an intraresidue proton normally correlated to a sugar H-1 signal orthogonally along the z axis prior to a Tr-ROESY mixing sequence. This virtually eliminates the incoherent intraresidue ROESY transfer, suppresses coherent TOCSY transfer, and markedly enhances the intensity of interresidue ROESY cross peaks.

New carbohydrate specificity and HIV-1 fusion blocking activity of the cyanobacterial protein MVL: NMR, ITC and sedimentation equilibrium studies

Carbohydrate-binding proteins that bind their carbohydrate ligands with high affinity are rare and therefore of interest because they expand our understanding of carbohydrate specificity and the structural requirements that lead to high-affinity interactions. Here, we use NMR and isothermal titration calorimetry techniques to determine carbohydrate specificity and affinities for a novel cyanobacterial protein, MVL, and show that MVL binds oligomannosides such as Man(6)GlcNAc(2) with sub-micromolar affinities. The amino acid sequence of MVL contains two homologous repeats, each comprising 54 amino acid residues. Using multi-dimensional NMR techniques, we show that MVL contains two novel carbohydrate recognition domains composed of four non-contiguous regions comprising approximately 15 amino acid residues each, and that these residues make numerous intermolecular contacts with their carbohydrate ligands. NMR screening of a comprehensive panel of di-, tri-, and high-mannose oligosaccharides establish that high-affinity binding requires at least the presence of a discrete conformation presented by Manbeta(1-->4)GlcNAc in the context of larger oligomannosides. As shown by sedimentation equilibrium and gel-filtration experiments, MVL is a monodisperse dimer in solution, and NMR data establish that the three-dimensional structure must be symmetric. MVL inhibits HIV-1 Envelope-mediated cell fusion with an IC(50) value of approximately 30 nM.

Molecular dynamics simulations of galectin-1-oligosaccharide complexes reveal the molecular basis for ligand diversity

Galectin-1 is a member of a protein family historically characterized by its ability to bind carbohydrates containing a terminal galactosyl residue. Galectin-1 is found in a variety of mammalian tissues as a homodimer of 14.5-kDa subunits. A number of developmental and regulatory processes have been attributed to the ability of galectin-1 to bind a variety of oligosaccharides containing the Gal-beta-(1,4)-GlcNAc (LacNAc(II)) sequence. To probe the origin of this permissive binding, solvated molecular dynamics (MD) simulations of several representative galectin-1-ligand complexes have been performed. Simulations of structurally defined complexes have validated the computational approach and expanded upon data obtained from X-ray crystallography and surface plasmon resonance measurements. The MD results indicate that a set of anchoring interactions between the galectin-1 carbohydrate recognition domain (CRD) and the LacNAc core are maintained for a diverse set of ligands and that substituents at the nonreducing terminus of the oligosaccharide extend into the remainder of a characteristic surface groove. The anionic nature of ligands exhibiting relatively high affinities for galectin-1 implicates electrostatic interactions in ligand selectivity, which is confirmed by a generalized Born analysis of the complexes. The results suggest that the search for a single endogenous ligand or function for this lectin may be inappropriate and instead support a more general role for galectin-1, in which the lectin is able to crosslink heterogeneous oligosaccharides displayed on a variety of cell surfaces. Such binding promiscuity provides an explanation for the variety of adhesion phenomena mediated by galectin-1.

Structural characterization of a mannose-binding protein-trimannoside complex using residual dipolar couplings

The ligand-binding properties of a 53 kDa homomultimeric trimer from mannose-binding protein (MBP) have been investigated using residual dipolar couplings (RDCs) that are easily measured from NMR spectra of the ligand and isotopically labeled protein. Using a limited set of 1H-15N backbone amide NMR assignments for MBP and orientational information derived from the RDC measurements in aligned media, an order tensor for MBP has been determined that is consistent with symmetry-based predictions of an axially symmetric system. 13C-1H couplings for a bound trisaccharide ligand, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside (trimannoside) have been determined at natural abundance and used as orientational constraints. The bound ligand geometry and orientational constraints allowed docking of the trimannoside ligand in the binding site of MBP to produce a structural model for MBP-oligosaccharide interactions.

Suitability of binary mixtures of water with aprotic solvents to turn hydroxyl protons of carbohydrate ligands into conformational sensors in NOE and transferred NOE experiments

The structural analysis of protein-carbohydrate interactions is essential for the long-range aim to sort out entropic/ enthalpic factors in the binding process. Of conspicuous clinical interest, this work can also offer the perspective to devise new classes of therapeuticals which interfere with disease-related glycan recognition. We have shown that it is possible to use exchangeable hydroxyl protons of carbohydrate ligands as conformational sensors for defining their bound-state topology by measurements in dimethyl sulfoxide(d6) (Siebert et al. (2000) ChemBioChem, 1, 181-195). However, the proteins are required to maintain binding capacity in the aprotic solvent. To define conditions to limit its harmful effect on sensitive protein structures while still being able to pick up solvent-exchangeable hydroxyl signals we systematically tested binary solvent mixtures of dimethyl sulfoxide and acetone with water. These solvent mixtures did not preclude to monitor hydroxyl protons of carbohydrate ligands even at temperatures well above 0 degrees C. Notably, hydrogen bonding of the two tested disaccharides (Galbeta1-4Glcalpha/beta and Galalpha1-3Galalpha/beta or Galalpha1-3Galbeta1-OCH(3)), which are common lectin ligands, resembled the situation under physiological conditions. Also, a refined topological description for hydroxyl positioning could be achieved for Galalpha1-3Gal. At least equally important, this approach worked for elucidation of the mistletoe-lectin-bound topology of lactose in its syn-conformation with indication for formation of a characteristic interresidual hydrogen bond. These measurements were performed in a binary dimethyl sulfoxide(d6):water mixture (6:4 ratio, v/v) at -12 degrees C and encourage to pursue this line of investigation by monitoring in the course of stepwise temperature increases. Our experiments reveal that binary mixtures have favorable properties for the conformational analysis of the free- and bound-state topologies of bioactive ligands.