Search

A mouse monoclonal antibody (mAb FL100A) previously prepared against Flavobacterium psychrophilum (Fp) CSF259‐93 has now been examined for binding to lipopolysaccharides (LPS) of this strain and Fp 950106‐1/1. The corresponding O‐polysaccharides (O‐PS) of these strains are formed by identical trisaccharide repeats composed of l‐Rhamnose (l‐Rha), 2‐acetamido‐2‐deoxy‐l‐fucose (l‐FucNAc) and 2‐acetamido‐4‐R1‐2,4‐dideoxy‐d‐quinovose (d‐Qui2NAc4NR1) where R1 represents a dihydroxyhexanamido moiety. The O‐PS loci of these strains are also identical except for the gene (wzy1 or wzy2) that encodes the polysaccharide polymerase. Accordingly, adjacent O‐PS repeats are joined through d‐Qui2NAc4NR1 and l‐Rha by wzy2‐dependent α(1–2) linkages in Fp CSF259‐93 versus wzy1‐dependent β(1–3) linkages in Fp 950106‐1/1. mAb FL100A reacted strongly with Fp CSF259‐93 O‐PS and LPS but weakly or not at all with Fp 950106‐1/1 LPS and O‐PS. Importantly, it also labelled cell surface blebs on the former but not the latter strain. Additionally, mAb binding was approximately 5‐times stronger to homologous Fp CSF259‐93 LPS than to LPS from a strain with a different R‐group gene. A conformational epitope for mAb FL100A binding was suggested from molecular dynamic simulations of each O‐PS. Thus, Fp CSF259‐93 O‐PS formed a stable well‐defined compact helix in which the R1 groups were displayed in a regular pattern on the helix exterior while unreactive Fp 950106‐1/1 O‐PS adopted a flexible extended linear conformation. Taken together, the findings establish the specificity of mAb FL100A for Wzy2‐linked F. psychrophilum O‐PS and LPS. [ABSTRACT FROM AUTHOR]

Siglecs are important lectins found in different types of immune cells and function as regulatory molecules by recognizing self-associated glycans and converting extracellular interactions into signals for inhibiting immune cell functions. Although many Siglecs have been found to show broad specificities and recognize different types of sulfated oligosaccharides, Siglec-8 and Siglec-9 displayed a high degree of specificity for sialyl N-acetyllactosamine (sLacNAc) with sulfations at O6-positions of the galactose (6’-sulfation) and N-acetylglucosamine (6-sulfation), respectively. Siglec-3 was recently discovered to bind sLacNAc both sulfations. In addition to a conserved arginine residue for binding to sialic acid residue, the sequence variety in the CC’ loop may provide binding specificities to sulfated oligosaccharides in Siglecs. Thus, the present study employed molecular models to study the impact of different residues in the CC’ loops of Siglec-8/9/3 to the recognitions of 6-sulfations in Gal and/or GlcNAc of sLacNAc. The negatively charged residues in the CC’ loop of Siglec-9 formed unfavorable electrostatic repulsions with the 6-sulfate in Gal and resulted no recognitions, in contrast to the favorable interactions formed between the positively charged residues in the CC’ loop of Siglec-8 and the 6-sulfate in Gal resulting strong specificity. A two-state binding model was proposed for Siglec-3 recognizing 6-sulfations in Gal and GlcNAc of sLacNAc, as the neutral residues in the CC’ loop of Siglec-3 could not form strong favorable interactions to lock the 6-sulfate in Gal within a single binding pose or strong unfavorable interactions to repel the 6-sulfate in Gal. The oligosaccharide adopted two distinctive binding poses and oriented the sulfate groups to form interactions with residues in the CC’ loop and G-strand. The present study provided a structural mechanism for the sequence variety in the CC’ loop of Siglec-8/9/3 determining the recognitions to the sulfated oligosaccharides and offered insights into the binding specificities for Siglecs.

Sulfated glycans have been found to be associated with various diseases and therefore have significant potential in molecular pathology as biomarkers. Although lectins are useful reagents for detecting glycans, there is a paucity of sulfate-recognizing lectins, and those that exist, such as from Maackia amurensis, display mixed specificities. Recombinant lectin engineering offers an emerging tool for creating novel glycan recognition by altering and/or enhancing endogenous specificities. The present study demonstrated the use of computational approaches in the engineering of a mutated form of E-selectin that displayed highly specific recognition of 60 -sulfo-sialyl Lewis X (60 -sulfosLex ), with negligible binding to its endogenous nonsulfated ligand, sLex . This new specificity mimics that of the unrelated protein Siglec-8, for which 60 -sulfo-sLex is its preferred ligand. Molecular dynamics simulations and energy calculations predicted that two point mutations (E92A/E107A) would be required to stabilize binding to the sulfated oligosaccharide with E-selectin. In addition to eliminating putative repulsions between the negatively charged side chains and the sulfate moiety, the mutations also abolished favorable interactions with the endogenous ligand. Glycan microarray screening of the recombinantly expressed proteins confirmed the predicted specificity change but also identified the introduction of unexpected affinity for the unfucosylated form of 60 -sulfo-sLex (60 -sulfo-sLacNAc). Three key requirements were demonstrated in this case for engineering specificity for sulfated oligosaccharide: 1) removal of unfavorable interactions with the 60 -sulfate, 2) introduction of favorable interactions for the sulfate, and 3) removal of favorable interactions with the endogenous ligand. [ABSTRACT FROM AUTHOR]

Background: Most licensed antiviral drugs are nucleoside analogs. A recent research focuses on blocking a virus from entering the cells in the viral cell adsorption/entry stage. In this entry mechanism, the glycans present on the viral surface play a fundamental role. Homochiral L-peptides acting this fusion mechanism have shown some inhibition of viral infection. Peptides with regularly alternating enantiomeric sequence (L,D-peptides) can assume structures that are not accessible to the corresponding homochiral molecules. Furthermore, L,D-peptides are less sensitive to enzymatic digestion. Aim: In silico design of an L, D-peptide with a high affinity for the viral surface glycans and consequently able to interfere with its fusion mechanism. Methods: We used a 3 α , 6 α -Mannopentaose (3-6MP) molecule to simulate a viral surface glycan. We performed molecular dynamics (MD) simulations of 3-6MP and D,L-peptide in water using the force field AMBER12-GLYCAM06i. We evaluated the binding constant from trajectories. The D,L-peptide molecule was modified over the sequence, the length, the terminals and finally glycosylated to attain a very high binding constant value for 3-6MP. In addition, we studied the specific interaction between T lymphocyte CD4 glycoprotein and HIV envelope glycoprotein gp120 through MD simulations between a D,L-peptide, bounded to a typical CD4 glycan, and a highly conserved HIV gp120 glycan. Results: Concerning the interaction with 3-6MP molecule, the very effective molecule obtained was H-d-Trp-l-Pro-d-Asn-l-Pro-d-Trp-l-Pro-d-Asn-l-Pro-OH where the Asn residues in position 3 and 7 are glycosylated with alpha-D-Mannopyranosyl-(1 → 3)-[alpha-D-mannopyranosyl-(1 → 6)]-alpha-D-mannopyranosyl-(1 → 4)-N-acetyl-beta-D-glucopyranosyl-(1 → 4)-N-acetyl-beta-D-glucopyranosyl-1-OH. As far as the interaction with HIV envelope is considered instead, the very effective molecule obtained, able to antagonize the CD4 glycoprotein, was H-d-Trp-l-Pro-d-Asn-l-Pro-d-Trp-l-Pro-d-Asn-l-Pro-OH where the Asn residue in position 3 is glycosylated with alpha-D-galactopyranosyl-(1 → 4)- N-acetyl-beta-D-glucopyranosyl-(1 → 2)-alpha-D-Mannopyranosyl-(1 → 3)- [alpha-D-galactopyranosyl-(1 → 4)- N-acetyl-beta-D-glucopyranosyl-(1 → 2)- alpha-D-Mannopyranosyl-(1 → 6)]- beta-D-mannopyranosyl-(1 → 4)-N-acetyl-beta-D-glucopyranosyl-(1 → 4)-N-acetyl-beta-D-glucopyranosyl-1-OH Conclusion: The above optimized glycosylated D,L-peptide molecules could be very promising representatives of a new powerful class of antiviral agents. In silico design of glyco-D,L-peptides with high affinity for the viral surface glycans and consequently able to interfere with its fusion mechanism. Molecular Dynamics (MD) simulations are performed using the GROMACS suite and the force field AMBER12-GLYCAM06i. A very high binding is realized between: the glyco-D,L-peptide TrpProAsnPenta and the common viral surface glycan 3-6MP [1]; the glyco-D,L-peptide TrpProAsn9glyco and the conserved gp120 HIV glycan N262glyco [2]. [ABSTRACT FROM AUTHOR]

Myelin Oligodendrocyte Glycoprotein (MOG) is a transmembrane protein in the myelin sheath. It acts as an auto-antigen under certain unknown conditions causing demyelination, thus resulting in Myelin Oligodendrocyte Glycoprotein Antibody-associated Disease (MOGAD). The significance of glycosylation in the conformational dynamics of the extracellular region (EC1) of the MOG were evident from the previous computational studies. Here, in this study, we performed the molecular dynamics simulation of the entire human MOG in the myelin sheath for 100 ns using the NAMD program. The results indicated that the EC1 and cytoplasmic region (CP) dominate the conformational rigidity of the protein, and enhance its interaction with lipids. This in turn helps in maintaining the myelin integrity in the presence of glycan. The transmembrane regions have reduced interaction with lipids in the glycosylated system. Moreover, the C-terminal extracellular region 2 (EC2) behaves exactly opposite to that of EC1 in the glycan presence. This may be attributed to the glycosylation site in the EC1 region. Hence, not only the region EC1 (having 3 crucial epitopes) but even the CP region were important for understanding the proper function of MOG in the glycan presence.Communicated by Ramaswamy H. Sarma.

While the construction of a dependable force field for performing classical molecular dynamics (MD) simulation is crucial for elucidating the structure and function of biomolecular systems, the attempts to do this for glycans are relatively sparse compared to those for proteins and nucleic acids. Currently, the use of GLYCAM06 force field is the most popular, but there have been a number of concerns about its accuracy in the systematic description of structural changes. In the present work, we focus on the improvement of the GLYCAM06 force field for β-d-glucose, a simple and the most abundant monosaccharide molecule, with the aid of machine learning techniques implemented with the TensorFlow library. Following the pre-sampling over a wide range of configuration space generated by MD simulation, the atomic charge and dihedral angle parameters in the GLYCAM06 force field were re-optimized to accurately reproduce the relative energies of β-d-glucose obtained by the density functional theory (DFT) calculations according to the structural changes. The validation for the newly proposed force-field parameters was then carried out by verifying that the relative energy errors compared to the DFT value were significantly reduced and that some inconsistencies with experimental (e.g., NMR) results observed in the GLYCAM06 force field were resolved relevantly.

Glycosaminoglycans (GAGs) are an important class of carbohydrates that serve critical roles in blood clotting, tissue repair, cell migration and adhesion, and lubrication. The variable sulfation pattern and iduronate ring conformations in GAGs influence their polymeric structure and nature of interaction. This study characterizes several heparin-like GAG disaccharides and tetrasaccharides using NMR and molecular dynamics simulations to assist in the development of parameters for GAGs within the GLYCAM06 force field. The force field additions include parameters and charges for a transferable sulfate group for O- and N-sulfation, neutral (COOH) forms of iduronic and glucuronic acid, and Δ4,5-unsaturated uronate (ΔUA) residues. ΔUA residues frequently arise from the enzymatic digestion of heparin and heparin sulfate. Simulations of disaccharides containing ΔUA reveal that the presence of sulfation on this residue alters the relative populations of 1H2 and 2H1 ring conformations. Simulations of heparin tetrasaccharides containing N-sulfation in place of N-acetylation on glucosamine residues influence the ring conformations of adjacent iduronate residues. [ABSTRACT FROM AUTHOR]

A two-state model is commonly used for interpreting ring conformations of furanoses based on NMR scalar J-coupling constants, with the ring populating relatively narrow distributions in the North and the South of the pseudorotation itinerary. The validity of this simple approach has been questioned, and is examined here in detail employing molecular dynamics (MD) simulations with a new GLYCAM force field parameter set for furanoses. Theoretical J-coupling constants derived from unrestrained MD simulations with the new furanose-specific parameters agreed with the experimental coupling constants to within 1 Hz on average. The results confirm that a two state model is a reasonable description for the ring conformation in the majority of methyl furanosides. However, in the case of methyl α-D-arabinofuranoside the ring populates a continuum of states from North to South via the eastern side of the pseudorotational itinerary. Two key properties are responsible for these differences. Firstly, East and West regions in β- and α-anomers, respectively, are destabilized by the absence of the anomeric effect. And, secondly, East or West conformations can be further destabilized by repulsive interactions among vicinal hydroxyl groups and ring oxygen atoms when the vicinal hydroxyl groups are in syn-configurations (such as in ribose and lyxose) more so than when in anti (arabinose, xylose). [ABSTRACT FROM AUTHOR]

The ABO blood group system is the most important blood type system in human transfusion medicine. Here, we explore the specificity of antibody recognition towards ABO blood group antigens using computational modeling and biolayer interferometry. Automated docking and molecular dynamics (MD) simulations were used to explore the origin of the specificity of an anti-blood group A antibody variable fragment (Fv AC1001). The analysis predicts a number of Fv-antigen interactions that contribute to affinity, including a hydrogen bond between a HisL49 and the carbonyl moiety of the GalNAc in antigen A. This interaction was consistent with the dependence of affinity on pH, as measured experimentally; at lower pH there is an increase in binding affinity. Binding energy calculations provide unique insight into the origin of interaction energies at a per-residue level in both the scFv and the trisaccharide antigen. The calculations indicate that while the antibody can accommodate both blood group A and B antigens in its combining site, the A antigen is preferred by approximately 4 kcal/mol, consistent with the lack of binding observed for the B antigen.

The ABO blood group system is the most important blood type system in human transfusion medicine. Here, we explore the specificity of antibody recognition toward ABO blood group antigens using computational modeling and bio layer interferometry. Automated docking and molecular dynamics simulations were used to explore the origin of the specificity of an anti-blood group A antibody variable fragment (Fiv AC1001). The analysis predicts a number of Fv-antigen interactions that contribute to affinity, including a hydrogen bond between a HisL49 and the carbonyl moiety of the GalNAc in antigen A. This interaction was consistent with the dependence of affinity on pH, as measured experimentally; at lower pH there is an increase in binding affinity. Binding energy calculations provide unique insight into the origin of interaction energies at a per-residue level in both the scFv and the trisaccharide antigen. The calculations indicate that while the antibody can accommodate both blood group A and B antigens in its combining site, the A antigen is preferred by 4 kcal/mol, consistent with the lack of binding observed for the B antigen. [ABSTRACT FROM AUTHOR]

Protein glycosylation with O-linked N-acetylglucosamine (O-GlcNAc) is a post-translational modification of serine/threonine residues in nucleocytoplasmic proteins. O-GlcNAc has been shown to play a role in many different cellular processes and O-GlcNAcylation is often found at sites that are also known to be phosphorylated. Unlike phosphorylation, O-GlcNAc levels are regulated by only two enzymes, O-GlcNAc transferase (OGT) and O-GlcNAc hydrolase (O-GlcNAcase or OGA). So far, no obvious consensus sequence has been found for sites of O-GlcNAcylation. Additionally, O-GlcNAcase recognizes and cleaves all O-GlcNAcylated proteins, independent of their sequence. In this work, we generate and analyze five models of O-GlcNAcylated peptides in complex with a bacterial OGA. Each of the five glycopeptides bind to OGA in a similar fashion, with OGA–peptide interactions primarily, but not exclusively, involving the peptide backbone atoms, thus explaining the lack of sensitivity to peptide sequence. Nonetheless, differences in peptide sequences, particularly at the −1 to −4 positions, lead to variations in predicted affinity, consistent with observed experimental variations in enzyme kinetics. The potential exists, therefore, to employ the present analysis to guide the development glycopeptide-specific inhibitors, or conversely, the conversion of OGA into a reagent that could target specific O-GlcNAcylated peptide sequences. [ABSTRACT FROM PUBLISHER]

ABSTRACT Molecular dynamics (MD) simulations of cellulose microfibrils are pertinent to the paper, textile, and biofuels industries for their unique capacity to characterize dynamic behavior and atomic-level interactions with solvent molecules and cellulase enzymes. While high-resolution crystallographic data have established a solid basis for computational analysis of cellulose, previous work has demonstrated a tendency for modeled microfibrils to diverge from the linear experimental structure and adopt a twisted conformation. Here, we investigate the dependence of this twisting behavior on computational approximations and establish the theoretical basis for its occurrence. We examine the role of solvent, the effect of nonbonded force field parameters [partial charges and van der Waals (vdW) contributions], and the use of explicitly modeled oxygen lone pairs in both the solute and solvent. Findings suggest that microfibril twisting is favored by vdW interactions, and counteracted by both intrachain hydrogen bonds and solvent effects at the microfibril surface. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 746-756, 2013. [ABSTRACT FROM AUTHOR]

Glycan chains that terminate in sialic acid (Neu5Ac) are frequently the receptors targeted by pathogens for initial adhesion. Carbohydrate-binding proteins (lectins) with specificity for Neu5Ac are particularly useful in the detection and isolation of sialylated glycoconjugates, such as those associated with pathogen adhesion as well as those characteristic of several diseases including cancer. Structural studies of lectins are essential in order to understand the origin of their specificity, which is particularly important when employing such reagents as diagnostic tools. Here, we report a crystallographic and molecular dynamics (MD) analysis of a lectin from Polyporus squamosus (PSL) that is specific for glycans terminating with the sequence Neu5Acα2-6Galβ. Because of its importance as a histological reagent, the PSL structure was solved (to 1.7 Å) in complex with a trisaccharide, whose sequence (Neu5Acα2-6Galβ1-4GlcNAc) is exploited by influenza A hemagglutinin for viral adhesion to human tissue. The structural data illuminate the origin of the high specificity of PSL for the Neu5Acα2-6Gal sequence. Theoretical binding free energies derived from the MD data confirm the key interactions identified crystallographically and provide additional insight into the relative contributions from each amino acid, as well as estimates of the importance of entropic and enthalpic contributions to binding. [ABSTRACT FROM PUBLISHER]

Abstract: The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°. [Copyright &y& Elsevier]

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5–6 of PECAM-1 and a region of high affinity in domains 2–3, consistent with experimental data and ligand–protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules. [ABSTRACT FROM PUBLISHER]

Abstract: Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl α-d-glucopyranoside and methyl α-d-galactopyranoside, employing the glycam06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the ω (O-5–C-5–C-6–O-6) and χ (C-3–C-4–O-4–HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4–HO-4 and O-6–HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics. [Copyright &y& Elsevier]

GLYCAM06 is a generalisable biomolecular force field that is extendible to diverse molecular classes in the spirit of a small-molecule force field. Here we report parameters for lipids, lipid bilayers and glycolipids for use with GLYCAM06. Only three lipid-specific atom types have been introduced, in keeping with the general philosophy of transferable parameter development. Bond stretching, angle bending, and torsional force constants were derived by fitting to quantum mechanical data for a collection of minimal molecular fragments and related small molecules. Partial atomic charges were computed by fitting to ensemble-averaged quantum-computed molecular electrostatic potentials. In addition to reproducing quantum mechanical internal rotational energies and experimental valence geometries for an array of small molecules, condensed-phase simulations employing the new parameters are shown to reproduce the bulk physical properties of a DMPC lipid bilayer. The new parameters allow for molecular dynamics simulations of complex systems containing lipids, lipid bilayers, glycolipids, and carbohydrates, using an internally consistent force field. By combining the AMBER parameters for proteins with the GLYCAM06 parameters, it is also possible to simulate protein-lipid complexes and proteins in biologically relevant membrane-like environments. [ABSTRACT FROM AUTHOR]

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 α- and β-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. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2008 [ABSTRACT FROM AUTHOR]

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 GBSlll 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. [ABSTRACT FROM AUTHOR]

Abstract: The conformational properties of the capsular polysaccharide (CPS) from group B Streptococcus serotype III (GBS III) are derived from 50ns explicitly solvated molecular dynamics simulations of a 25-residue fragment of the CPS. The results from the simulations are shown to be consistent with experimental NMR homo- and heteronuclear J-coupling and NOE data for both the sialylated native CPS and for the chemically desialylated polysaccharide. A helical structure is predicted with a diameter of 29.3Å and a pitch 89.5Å, in which the sialylated side chains are arrayed on the exterior surface of the helix. The results provide an explanation for the observation that CPS antigenicity varies with carbohydrate chain length up to approximately 4 pentasaccharide repeat units. The conformation of the immunodominant region is established and shown to be independent of the presence of sialic acid. The data provide an explanation for the observation that the specificity of the determinant, associated with the major population of antibodies raised upon immunization of rabbits with GBS III, is dependent on the presence of sialic acid. In the sialylated native CPS, the antibody response is largely directed against the immunodominant core of the helix. From simulations of the desialylated CPS, a model emerges which suggests that the minor population of antibodies, whose determinant is not sialic acid dependent, recognizes the same immunodominant region, but that in the disordered CPS this region is not presented in a regular repeating motif. [Copyright &y& Elsevier]

A program for converting the different existing AMBER and GLYCAM force fields for use with commercial molecular modeling packages is presented, using the Molecular Simulations Inc. (MSI or Accelrys) software package as a case model. Called AmberFFC, the program creates AMBER and GLYCAM force field files suitable for use with the Accelrys molecular mechanics modules by transforming the amino acid, nucleotide, and monosaccharide topology databases and force field parameter files to the Accelrys file format. It is intended for any modeler who is interested in using the current AMBER and GLYCAM force fields with the Accelrys FDiscover and CDiscover programs. AmberFFC has been written entirely with the Perl programming language, making it highly flexible and portable. In order to compare the implementation of the force fields converted by AmberFFC in the Accelrys package with their corresponding execution in the AMBER software, and also to verify the efficiency of the AmberFFC program, results from single point energy calculations for 13 model molecules were obtained with the two programs. It is demonstrated that results obtained with the CDiscover and FDiscover modules compare well to those found using Sander_classic, thus showing that AmberFFC is a highly efficient program. Some energy differences between the AMBER and Accelrys software have been observed, and their origin has been characterized and discussed. [ABSTRACT FROM AUTHOR]

While the construction of a dependable force field for performing classical molecular dynamics (MD) simulation is crucial for elucidating the structure and function of biomolecular systems, the attempts to do this for glycans are relatively sparse compared to those for proteins and nucleic acids. Currently, the use of GLYCAM06 force field is the most popular, but there have been a number of concerns about its accuracy in the systematic description of structural changes. In the present work, we focus on the improvement of the GLYCAM06 force field for β -d-glucose, a simple and the most abundant monosaccharide molecule, with the aid of machine learning techniques implemented with the TensorFlow library. Following the pre-sampling over a wide range of configuration space generated by MD simulation, the atomic charge and dihedral angle parameters in the GLYCAM06 force field were re-optimized to accurately reproduce the relative energies of β -d-glucose obtained by the density functional theory (DFT) calculations according to the structural changes. The validation for the newly proposed force-field parameters was then carried out by verifying that the relative energy errors compared to the DFT value were significantly reduced and that some inconsistencies with experimental (e.g., NMR) results observed in the GLYCAM06 force field were resolved relevantly. [ABSTRACT FROM AUTHOR]

The conformational properties of oligosaccharides are important in determining their biological properties, such as recognition by proteins. The structural and dynamic properties of many oligosaccharides are poorly understood both because of a lack of experimental data (usually obtained from solution NMR parameters) and because of gross approximations frequently invoked in theoretical models. To characterise the oligomannose oligosaccharide Man9GlcNAc2 we have acquired a more extensive NMR data set and performed the first unrestrained molecular dynamics (MD) simulation in water of this large oligosaccharide (employing the GLYCAM—93 parameter set with the AMBER force field). Good agreement is seen between the computed dynamics data and the results of both an isolated spin pair (ISPA) analysis of short mixing time NOE data and NOE build-up curves for mixing times from 100 to 2000 ms. The number of experimental conformational constraints obtained in this study are in principle sufficient to fully define a rigid structure. The fact that this could not be done indicates a high degree of internal flexibility and/or the presence of multiple conformations about the glycosidic linkages. Independently, the same conclusions are reached from an analysis of the MD results. In addition, the theoretical results allow the overall topology of the molecule and its intra-molecular and solvent-mediated hydrogen bonding pattern to be defined. Extensive re-organisation of solvent and inter-residue hydrogen bonds is shown to be required for significant conformational changes to occur, resulting in relatively long life-times for distinct glycosidic linkage conformations, despite the high local flexibility of the glycosidic linkages. This factor is also seen in the overall topology of the molecule, where the considerable internal flexibility is not translated into gross changes in structure. The control exerted by the solvent over both the flexibility and... [ABSTRACT FROM AUTHOR]

Computational methods have had a long history of application to carbohydrate systems and their development in this regard is discussed. The conformational analysis of carbohydrates differs in several ways from that of other biomolecules. Many glycans appear to exhibit numerous conformations coexisting in solution at room temperature and a conformational analysis of a carbohydrate must address both spatial and temporal properties. When solution nuclear magnetic resonance data are used for comparison, the simulation must give rise to ensemble-averaged properties. In contrast, when comparing to experimental data obtained from crystal structures a simulation of a crystal lattice, rather than of an isolated molecule, is appropriate. Molecular dynamics simulations are well suited for such condensed phase modeling. Interactions between carbohydrates and other biological macromolecules are also amenable to computational approaches. Having obtained a three-dimensional structure of the receptor protein, it is possible to model with accuracy the conformation of the carbohydrate in the complex. An example of the application of free energy perturbation simulations to the prediction of carbohydrate-protein binding energies is presented. [ABSTRACT FROM AUTHOR]

Myelin Oligodendrocyte glycoprotein (MOG) is found to play an important role in providing structural integrity to myelin sheath at the same time it acts as an auto-antigen which might lead to Multiple Sclerosis (MS). What causes this specific property of being an auto-antigen is still not known. Here we present molecular dynamics simulation studies of unfolding and folding of the protein MOG in both the absence and presence of N-glycan in order to understand the role of glycosylation in the stability and flexibility of the protein. The main results from these studies show that the glycosylation increases the stability of the protein MOG and inhibits the complete unfolding of MOG in the SMD. From the folding studies using TMD, it was observed that the glycan helps the protein to attain the near-native folded conformation. However, it was also observed from the direct TMD studies that the pathway of protein folding was enhanced by the trace-back of intermediate states in the presence of glycan. Image 1 • Effect of glycan in Myelin Oligodendrocyte Glycoprotein unfolding and folding. • SMD and TMD were used to study the protein unfolding and folding. • The major results of our study were • Glycosylation maintains the stability of protein MOG. • Glycan enhances the protein intra-molecular interaction and helps in proper folding. [ABSTRACT FROM AUTHOR]

Myelin Oligodendrocyte Glycoprotein (MOG) is found on the external surface of the myelin sheath and plays an important role in neurodegenerative diseases. It was observed that the protein MOG acts as an autoantigen and results in demyelination. The cause for the sudden change of protein to be autoantigen is still unclear. Here we present the molecular dynamics simulation studies of MOG in both unbound and bound states with an antibody. Both these systems were studied in the absence and presence of N-glycan in order to understand the effect of glycosylation in the MOG conformational changes. The results indicate that the glycosylation decreases the flexibility of protein in both free and bound states. Glycan influence the interaction of the complex with the water molecules whereas free protein MOG interaction with water molecules was not affected by the glycosylation. Glycan changes the 3 10 helices adjacent to the antibody interacting epitope MOG 35-55 to turns.Communicated by Ramaswamy H. Sarma.

Computational methodologies have immense potential to delineate the dynamic conformations of glycoconjugates including gangliosides, thereby characterizing the conformational adaptability of their glycans upon interacting with various target proteins. Replica-exchange molecular dynamics simulations have been employed to effectively explore the vast conformational spaces of large, branched carbohydrate moieties. When experimentally validated using NMR, molecular simulations can provide dynamical views of molecular recognition events involving the ganglioside glycans.

Glycan structures on glycoproteins and glycolipids play critical roles in biological recognition, targeting, and modulation of functions in animal systems. Many classes of glycan structures are capped with terminal sialic acid residues, which contribute to biological functions by either forming or masking glycan recognition sites on the cell surface or secreted glycoconjugates. Sialylated glycans are synthesized in mammals by a single conserved family of sialyltransferases that have diverse linkage and acceptor specificities. We examined the enzymatic basis for glycan sialylation in animal systems by determining the crystal structures of rat ST6GAL1, an enzyme that creates terminal α2,6-sialic acid linkages on complex-type N-glycans, at 2.4 Å resolution. Crystals were obtained from enzyme preparations generated in mammalian cells. The resulting structure revealed an overall protein fold broadly resembling the previously determined structure of pig ST3GAL1, including a CMP-sialic acid-binding site assembled from conserved sialylmotif sequence elements. Significant differences in structure and disulfide bonding patterns were found outside the sialylmotif sequences, including differences in residues predicted to interact with the glycan acceptor. Computational substrate docking and molecular dynamics simulations were performed to predict and evaluate the CMP-sialic acid donor and glycan acceptor interactions, and the results were compared with kinetic analysis of active site mutants. Comparisons of the structure with pig ST3GAL1 and a bacterial sialyltransferase revealed a similar positioning of donor, acceptor, and catalytic residues that provide a common structural framework for catalysis by the mammalian and bacterial sialyltransferases.

In the calculation of partial atomic charges, for use in molecular mechanics or dynamics simulations, it is common practice to select only a single conformation for the molecule of interest. For molecules that contain rotatable bonds, it is preferable to compute the charges from several relevant conformations. We present here results from a charge derivation protocol that determines the partial charges by averaging charges computed for conformations selected from explicitly solvated MD simulations, performed under periodic boundary conditions. This approach leads to partial charges that are weighted by a realistic population of conformations and that are suitable for condensed phase simulations. This protocol can, in principle, be applied to any class of molecule and to nonaqueous solvation. Carbohydrates contain numerous hydroxyl groups that exist in an ensemble of orientations in solution, and in this report we apply ensemble averaging to a series of methyl glycosides. We report the extent to which ensemble averaging leads to charge convergence among the various monosaccharides and among the constituent atoms within a given monosaccharide. Due to the large number of conformations (200) in our ensembles, we are able to compute statistically relevant standard deviations for the partial charges. An analysis of the standard deviations allows us to assess the extent to which equivalent atom types may, nevertheless, require unique partial charges. The configurations of the hydroxyl groups exert considerable influence on internal energies, and the limits of ensemble averaged charges are discussed in terms of these properties.

Charges derived from fitting a classical Coulomb model to quantum mechanical molecular electrostatic potentials (so called ESP-charges) are frequently used in simulations of macromolecules. Simulational methods that use ESP-charges generally reproduce the geometries of hydrogen bonded complexes, despite the fact that these charges are known to overestimate the strengths of these interactions. Through the use of a restraint function during the fitting of the partial charges to the electrostatic potentials the magnitudes of the charges may be attenuated (so called RESP-charges). For the AMBER force field RESP-charges have been proposed for proteins and nucleic acids. Here we examine a novel approach for determining the RESP-charges for carbohydrates based on molecular dynamics (MD) simulations of crystal structures. During a simulation, the crystallographic unit cell geometry is sensitive to both inter-molecular non-bonded forces and internal torsional rotations. However, for polar molecules, and specifically carbohydrates, the crystal geometries are particularly sensitive to the set of partial atomic charges employed in the simulation. Thus, given a force field in which the van der Waals and torsion terms are well parameterized, it is possible to assess the suitability of a set of partial charges by monitoring the properties of the crystal during an MD simulation. We have examined several charge sets for use with the GLYCAM parameters for carbohydrate and glycoprotein simulations and found that a restraint weight of 0.01 gives the best agreement with the neutron diffraction structure of α-d-glucopyranose. Unrestrained ESP-charges performed poorly as did the charges obtained from Mulliken and distributed multipole analyses of the quantum mechanical HF/6-31G* wavefunctions.