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.

The rotational preferences of the hydroxymethyl group in pyranosides is known to depend on the local environment, whether in solid, solution, or gas phase. By combining molecular dynamics (MD) simulations with NMR spectroscopy the rotational preferences for the ω angle in methyl 2,3-di--methyl-α-D-glucopyranoside () and methyl 2,3-di--methyl-α-D-galactopyranoside () in a variety of solvents, with polarities ranging from 80 to 2.3 D have been determined. The effects of solvent polarity on intramolecular hydrogen bonding have been identified and quantified. In water, the internal hydrogen bonding networks are disrupted by competition with hydrogen bonds to the solvent. When the internal hydrogen bonds are differentially disrupted, the rotamer populations associated with the ω angle may be altered. In the case of in water, the preferential disruption of the interaction between HO6 and O4 destabilizes the rotamer, leading to the observed preference for gauche rotamers. Without the hydrogen bond enhancement offered by a low polarity environment, both and display rotamer populations that are consistent with expectations based on the minimization of repulsive intramolecular oxygen-oxygen interactions. In a low polarity environment, HO6 prefers to interact with O4, however, in water these interactions are markedly weakened, indicating that HO6 acts as a hydrogen bond donor to water.

No abstract

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. S treptococcus agalactiae [Group B Streptococcus (GBS)] andStreptococcus pneumoniae (Pn) are responsible for the majority of life-threatening cases of septicemia, meningitis, and pneumonia in neonates (1, 2). Gram-positive bacteria, such as GBS and Pn, are classified into serotypes according to the unique carbohydrate sequence of the bacterial surface capsular polysaccharide (CPS) and protein antigens. Serotypes vary in antigenicity, immunogenicity, virulence, and geographical distribution (3). Quantification of the structural and dynamic properties responsible for the affinity and specificity of antigenic oligosaccharide-antibody (Ab) interactions is a crucial step in furthering the understanding of the immune response to bacterial and fungal pathogens. In GBS, the CPS is a high-molecularweight polymer composed of varying sequences of ␤-D-, and sometimes L-rhamnopyranose. The glyceryl side chain of the Neu5Ac residues may also be Oacetylated (4). In all GBS strains identified to date, the Neu5Ac residues occur in the terminal position on the side-chain branches of the polymeric repeat unit of the CPS. They play an important role in defining the antigenicity and immunogenicity of the CPS (5). Variations within the CPS sequence result in the nine kno...

An NMR study of the diaxial/diequatorial chair equilibrium in a range of silylated derivatives of trans-1,4- and trans-1,2-dihydroxycyclohexane is reported and discussed with a view to explaining unusually large populations of chair conformations with axial substituents, noted previously for some monosilyloxycyclohexanes and in some silylated sugars. X-ray diffraction studies of three bis-triphenylsilyloxycyclohexanes are reported and show both axial and equatorial silyloxy groups with the exocyclic bonds eclipsed. Eclipsing is also suggested by molecular mechanics (MM3) calculations on such derivatives. Both axial and equatorial tertiary silyl groups have 1,3-repulsive interactions with whatever substituents or hydrogen atoms are at the two adjacent equatorial positions, and these are relieved by rotation toward the eclipsed conformation of the exocyclic C-O bond. The three substituents on silicon interact attractively with the nine atoms at the 3, 4, and 5-positions of the cyclohexane ring and calculations suggest that these stabilizing interactions are significantly greater in the axial than in the equatorial conformation. An equatorial C-OSiR(3) bond with one or two equatorial neighbors has a restricted potential energy well that becomes much broader when the bond is axial without any equatorial neighbors in the alternative chair. Adjacent silyl groups in the 1,2-disubstituted series interact in a stabilizing way overall in all conformations, this being particularly marked in the diaxial conformation of the more complex ethers. These factors lead to unusually large axial populations.