•Influenza specificity is influenced by differences in glycan conformational entropy.•Biantennary α2-6 sialosides can boost avidity by forming bidentate HA complexes.•All reported HA-oligosaccharide complexes can be grouped into one of four topologies.•Interpretations of specificity must consider the impact of changes in glycan substructure.Understanding the molecular origin of influenza receptor specificity is complicated by the paucity of quantitative affinity measurements, and the qualitative and variable nature of glycan array data. Further obstacles arise from the varied impact of viral glycosylation and the relatively narrow spectrum of biologically relevant receptors present on glycan arrays. A survey of receptor conformational properties is presented, leading to the conclusion that conformational entropy plays a key role in defining specificity, as does the newly reported ability of biantennary receptors that terminate in Siaα2-6Gal sequences to form bidentate interactions to two binding sites in a hemagglutinin trimer. Bidentate binding provides a functional explanation for the observation that Siaα2-6 receptors adopt an open-umbrella topology when bound to hemagglutinins from human-infective viruses, and calls for a reassessment of virus avidity and tissue tropism.

Molecular docking programs are primarily designed to align rigid, drug-like fragments into the binding sites of macromolecules and frequently display poor performance when applied to flexible carbohydrate molecules. A critical source of flexibility within an oligosaccharide is the glycosidic linkages. Recently, Carbohydrate Intrinsic (CHI) energy functions were reported that attempt to quantify the glycosidic torsion angle preferences. In the present work, the CHI-energy functions have been incorporated into the AutoDock Vina (ADV) scoring function, subsequently termed Vina-Carb (VC). Two user-adjustable parameters have been introduced, namely, a CHI- energy weight term (chi_coeff) that affects the magnitude of the CHI-energy penalty and a CHI-cutoff term (chi_cutoff) that negates CHI-energy penalties below a specified value. A data set consisting of 101 protein–carbohydrate complexes and 29 apoprotein structures was used in the development and testing of VC, including antibodies, lectins, and carbohydrate binding modules. Accounting for the intramolecular energies of the glycosidic linkages in the oligosaccharides during docking led VC to produce acceptable structures within the top five ranked poses in 74% of the systems tested, compared to a success rate of 55% for ADV. An enzyme system was employed in order to illustrate the potential application of VC to proteins that may distort glycosidic linkages of carbohydrate ligands upon binding. VC represents a significant step toward accurately predicting the structures of protein–carbohydrate complexes. Furthermore, the described approach is conceptually applicable to any class of ligands that populate well-defined conformational states.

We report a new classification method for pyranose ring conformations called Best-fit, Four-Membered Plane (BFMP), which describes pyranose ring conformations based on reference planes defined by four atoms. The method is able to characterize all asymmetrical and symmetrical shapes of a pyran ring, is readily automated, easy to interpret, and maps trivially to IUPAC definitions. It also provides a qualitative measurement of the distortion of the ring. Example applications include the analysis of data from crystal structures and molecular dynamics simulations.

Previous studies of calculated diffraction patterns for cellulose crystallites suggest that distortions that arise once models have been subjected to molecular dynamics (MD) simulation are the result of both microfibril twisting and changes in unit cell dimensions induced by the empirical force field; to date, it has not been possible to separate the individual contributions of these effects. To provide a better understanding of how twisting manifests in diffraction data, the present study demonstrates a method for generating twisted and linear cellulose structures that can be compared without the bias of dimensional changes, allowing assessment of the impact of twisting alone. Analysis of unit cell dimensions, microfibril volume, hydrogen bond patterns, glycosidic torsion angles, and hydroxymethyl group orientations confirmed that the twisted and linear structures collected with this method were internally consistent, and theoretical powder diffraction patterns for the two were shown to be effectively indistinguishable. These results indicate that differences between calculated patterns for the crystal coordinates and twisted structures from MD simulation can result entirely from changes in unit cell dimensions, and not from microfibril twisting. Although powder diffraction patterns for models in the 81-chain size regime were shown to be unaffected by twisting, suggesting that a modest degree of twist is not inconsistent with available crystallographic data, it may be that other diffraction techniques are capable of detecting this structural difference. Until such time as definitive experimental evidence comes to light, the results of this study suggest that both twisted and linear microfibrils may represent an appropriate model for cellulose Iβ.

The ability of ligands to displace conserved water molecules in protein binding sites is of significant interest in drug design and is particularly pertinent in the case of glycomimetic drugs. This concept was explored in previous work [Clarke et al. J. Am. Chem. Soc. 2001, 123, 12238−12247 and Kadirvelraj et al. J. Am. Chem. Soc. 2008, 130, 16933–16942] for a highly conserved water molecule located in the binding site of the prototypic carbohydrate-binding protein Concanavalin A (Con A). A synthetic ligand was designed with the aim of displacing such water. While the synthetic ligand bound to Con A in an analogous manner to that of the natural ligand, crystallographic analysis demonstrated that it did not displace the conserved water. In order to quantify the affinity of this particular water for the Con A surface, we report here the calculated standard binding free energy for this water in both ligand-bound and free Con A, employing three popular water models: TIP3P, TIP4P, and TIP5P. Although each model was developed to perform well in simulations of bulk-phase water, the computed binding energies for the isolated water molecule displayed a high sensitivity to the model. Both molecular dynamics simulation and free energy results indicate that the choice of water model may greatly influence the characterization of surface water molecules as conserved (TIP5P) or not (TIP3P) in protein binding sites, an observation of considerable significance to rational drug design. Structural and theoretical aspects at the basis of the different behaviors are identified and discussed.

Challenges for structural characterization of membrane-bound glycosphingolipids include their high internal dynamic motions and their physical proximity to membrane surfaces. Here we demonstrate that NMR paramagnetic relaxation enhancement can be used, alongside independent molecular dynamics simulations and an outer-sphere relaxation model, to quantitatively characterize the presentation (insertion depth and orientation relative to a membrane surface) of ganglioside GM1 in biologically relevant membrane environments. Longitudinal and transverse paramagnetic relaxation enhancement effects were measured for GM1, anchored to phospholipid bicelles, using both water-soluble and membrane-anchored paramagnetic probes, respectively. A method was developed to rapidly calculate paramagnetic relaxation enhancement effects from thousands of structures taken from a simulation of GM1 in a phospholipid bilayer. The combined computational and experimental approach yielded experimentally verified atomic-resolution 3D models of a highly plastic membrane-bound biomolecule.

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.

Biomolecular surface mapping methods offer an important alternative method for characterizing protein-protein and protein-ligand interactions in cases in which it is not possible to determine high-resolution three-dimensional (3D) structures of complexes. Hydroxyl radical footprinting offers a significant advance in footprint resolution compared with traditional chemical derivatization. Here we present results of footprinting performed with hydroxyl radicals generated on the nanosecond time scale by laser-induced photodissociation of hydrogen peroxide. We applied this emerging method to a carbohydrate-binding protein, galectin-1. Since galectin-1 occurs as a homodimer, footprinting was employed to characterize the interface of the monomeric subunits. Efficient analysis of the mass spectrometry data for the oxidized protein was achieved with the recently developed ByOnic (Palo Alto, CA) software that was altered to handle the large number of modifications arising from side-chain oxidation. Quantification of the level of oxidation has been achieved by employing spectral intensities for all of the observed oxidation states on a per-residue basis. The level of accuracy achievable from spectral intensities was determined by examination of mixtures of synthetic peptides related to those present after oxidation and tryptic digestion of galectin-1. A direct relationship between side-chain solvent accessibility and level of oxidation emerged, which enabled the prediction of the level of oxidation given the 3D structure of the protein. The precision of this relationship was enhanced through the use of average solvent accessibilities computed from 10 ns molecular dynamics simulations of the protein.

Bacterial surface capsular polysaccharides (CPS) that are similar in carbohydrate sequence may differ markedly in immunogenicity and antigenicity. The structural origin of these phenomena is poorly understood. Such a case is presented by the Gram-positive bacteria Streptococcus agalactiae (Group B Streptococcus; GBS) type III (GBSIII) and Streptococcus pneumoniae (Pn) type 14 (Pn14), which share closely related CPS sequences. Nevertheless, antibodies (Abs) against GBSIII rarely cross-react with the CPS from Pn14. To establish the origin for the variation in CPS antigenicity, models for the immune complexes of CPS fragments from GBSIII and Pn14, with the variable fragment (Fv) of a GBS-specific mAb (mAb 1B1), are presented. The complexes are generated through a combination of comparative Ab modeling and automated ligand docking, followed by explicitly solvated 10-ns molecular dynamics simulations. The relationship between carbohydrate sequence and antigenicity is further quantified through the computation of interaction energies using the Molecular Mechanics–Generalized Born Surface Area (MM-GBSA) method, augmented by conformational entropy estimates. Despite the electrostatic differences between Pn14 and GBSIII CPS, analysis indicates that entropic penalties are primarily responsible for the loss of affinity of the highly flexible Pn14 CPS for mAb 1B1. The similarity of the solution conformation of the relatively rigid GBSIII CPS with that in the immune complex characterizes the previously undescribed 3D structure of the conformational epitope. The analysis provides a comprehensive interpretation for a large body of biochemical and immunological data related to Ab recognition of bacterial polysaccharides and should be applicable to other Ab–carbohydrate interactions.

The relationship between the three-dimensional structures of oligosaccharides and polysaccharides and their biological properties has been the focus of many recent studies. The overall conformation of an oligosaccharide depends primarily on the orientation of the torsion angles (φ, ψ, and ω) between glycosyl residues. Numerous experimental studies have shown that in glucopyranosides the ω-torsion angle (O6-C6-C5-O5) displays a preference for gauche orientations, in disagreement with predictions based on gas-phase quantum mechanics calculations. In contrast, the ω-angle in galactopyranosides displays a high proportion of the anti-orientation. For oligosaccharides containing glycosidic linkages at the 6-position (1→6 linked), variations in rotamer population have a direct effect on the oligosaccharides' structure and function, and yet the physical origin of these conformational preferences remains unclear. Although it is generally recognized that the gauche effect in carbohydrates is a solvent-dependent phenomenon, the mechanism through which solvent induces the gauche preference is not understood. In the present work, quantum mechanics and solvated molecular dynamics calculations were performed on two representative carbohydrates, methyl α-d-glucopyranoside and methyl α-d-galactopyranoside. We show that correct reproduction of the experimental rotamer distributions about the ω-angles is obtained only after explicit water is included in the molecular dynamics simulations. The primary role of the water appears to be to disrupt the hydrogen bonding within the carbohydrate, thereby allowing the rotamer populations to be determined by internal electronic and steric repulsions between the oxygen atoms. The results reported here provide a quantitative explanation of the conformational behavior of (1→6)-linked carbohydrates.

The human immune response to an infection by Gram-negative bacteria involves detection of lipopolysaccharides (LPS), also known as endotoxins, which comprise the bacterial outer cell wall. Distinct from mammalian glycolipid structures, LPS have a conserved chemical pattern that is recognized by the pattern recognition receptor complex formed by myeloid differentiation protein 2 (MD-2) and toll-like receptor 4 (TLR4). A remarkable immune-mediated structure-toxicity relationship has been defined that relates to the number of acyl chains in the endotoxin. While there is a clear correlation between endotoxin acylation and elicited agonist or antagonist responses, the 3D structural basis of this relationship remains unclear. In order to explore, at atomic-resolution, the effects of a range of chemically distinct endotoxins on the structure and dynamics of their MD-2·endotoxin complexes, we examined a series of variably acylated lipid A molecules from Escherichia coli and Neisseria meningitidis in complex with human MD-2. Through the application of molecular dynamics simulations, in concert with experimental data, we have identified specific structural and dynamic features of the MD-2-endotoxin complexes that may control dimerization of TLR4 molecules. As dimerization is central to the release of downstream chemical mediators, the results provide a structural foundation for the ability of endotoxins to act as either agonists or antagonists of the TLR4 pathway.Graphical abstractDownload full-size imageHighlights► MD-2 complexes with variably acylated endotoxins share similar 2° and 3° structure. ► These complexes have markedly different quaternary structure. ► Placement of endotoxin in binding cavity is regulated by volume of the lipid domain. ► Both dynamic and structural properties of the complex are functionally important. ► TLR4 activation correlates with exposed apolar surface in MD-2·endotoxin complexes.

Biomolecular surface mapping methods offer an important alternative method for characterizing protein–protein and protein–ligand interactions in cases in which it is not possible to determine high-resolution three-dimensional (3D) structures of complexes. Hydroxyl radical footprinting offers a significant advance in footprint resolution compared with traditional chemical derivatization. Here we present results of footprinting performed with hydroxyl radicals generated on the nanosecond time scale by laser-induced photodissociation of hydrogen peroxide. We applied this emerging method to a carbohydrate-binding protein, galectin-1. Since galectin-1 occurs as a homodimer, footprinting was employed to characterize the interface of the monomeric subunits. Efficient analysis of the mass spectrometry data for the oxidized protein was achieved with the recently developed ByOnic (Palo Alto, CA) software that was altered to handle the large number of modifications arising from side-chain oxidation. Quantification of the level of oxidation has been achieved by employing spectral intensities for all of the observed oxidation states on a per-residue basis. The level of accuracy achievable from spectral intensities was determined by examination of mixtures of synthetic peptides related to those present after oxidation and tryptic digestion of galectin-1. A direct relationship between side-chain solvent accessibility and level of oxidation emerged, which enabled the prediction of the level of oxidation given the 3D structure of the protein. The precision of this relationship was enhanced through the use of average solvent accessibilities computed from 10 ns molecular dynamics simulations of the protein.A combination of oxidative surface footprinting and molecular dynamics simulation allows identification of the interfacial residues in dimeric galectin-1.Download high-res image (185KB)Download full-size image