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.

Asparagine-linked glycosylation, the co-translational covalent attachment of carbohydrates to asparagine side chains, has a major effect on the folding, stability, and function of many proteins. The carbohydrate composition in mature glycoproteins is heterogeneous due to modification of the initial oligosaccharide by glycosidases and glycosyltransferases during the glycoprotein passage through the endoplasmic reticulum and Golgi apparatus. Despite the diversity of carbohydrate structures, the core β-D-(GlcNAc) 2 remains conserved in all N-linked glycoproteins. Previously, results from our laboratory showed that the molecular composition of the core disaccharide has a critical and unique conformational effect on the peptide backbone. Herein, we employ a synergistic experimental and computational approach to study the effect of the stereochemistry of the carbohydrate-peptide linkage on glycopeptide structure. A glycopeptide derived from a hemagglutinin protein fragment was synthesized, with the carbohydrate attached to the peptide with an α-linked stereochemistry. Computational and biophysical analyses reveal that the conformations of the peptide and α-and β-linked glycopeptides are uniquely influenced by the attached saccharide. The value of computational approaches for probing the influence of attached saccharides on polypeptide conformation is highlighted.

The hydration behavior of two model disaccharides, methyl-alpha-D-maltoside (1) and methyl-alpha-D-isomaltoside (2), has been investigated by a comparative 10 ns molecular dynamics study. The detailed hydration of the two disaccharides was described using three force fields especially developed for modeling of carbohydrates in explicit solvent. To validate the theoretical results the two compounds were synthesized and subjected to 500 MHz NMR spectroscopy, including pulsed field gradient diffusion measurements (1: 4.0. 10(-6) cm(2). s(-1); 2: 4.2. 10(-6) cm(2). s(-1)). In short, the older CHARMM-based force field exhibited a more structured carbohydrate-water interaction leading to better agreement with the diffusional properties of the two compounds, whereas especially the alpha-(1-->6) linkage and the primary hydroxyl groups were inaccurately modeled. In contrast, the new generation of the CHARMM-based force field (CSFF) and the most recent version of the AMBER-based force field (GLYCAM-2000a) exhibited less structured carbohydrate-water interactions with the result that the diffusional properties of the two disaccharides were underestimated, whereas the simulations of the alpha-(1-->6) linkage and the primary hydroxyl groups were significantly improved and in excellent agreement with homo- and heteronuclear coupling constants. The difference between the two classes of force field (more structured and less structured carbohydrate-water interaction) was underlined by calculation of the isotropic hydration as calculated by radial pair distributions. At one extreme, the radial O em leader O pair distribution function yielded a peak density of 2.3 times the bulk density in the first hydration shell when using the older CHARMM force field, whereas the maximum density observed in the GLYCAM force field was calculated to be 1.0, at the other extreme.

We have performed a series of first-principles electronic structure calculations to examine the reaction pathways and the corresponding free energy barriers for the ester hydrolysis of protonated cocaine in its chair and boat conformations. The calculated free energy barriers for the benzoyl ester hydrolysis of protonated chair cocaine are close to the corresponding barriers calculated for the benzoyl ester hydrolysis of neutral cocaine. However, the free energy barrier calculated for the methyl ester hydrolysis of protonated cocaine in its chair conformation is significantly lower than for the methyl ester hydrolysis of neutral cocaine and for the dominant pathway of the benzoyl ester hydrolysis of protonated cocaine. The significant decrease of the free energy barrier, ~ 4 kcal/mol, is attributed to the intramolecular acid catalysis of the methyl ester hydrolysis of protonated cocaine, because the transition state structure is stabilized by the strong hydrogen bond between the carbonyl oxygen of the methyl ester moiety and the protonated tropane N. The relative magnitudes of the free energy barriers calculated for different pathways of the ester hydrolysis of protonated chair cocaine are consistent with the experimental kinetic data for cocaine hydrolysis under physiologic conditions. Similar intramolecular acid catalysis also occurs for the benzoyl ester hydrolysis of (protonated) boat cocaine in the physiologic condition, although the contribution of the intramolecular hydrogen bonding to transition state stabilization is negligible. Nonetheless, the predictability of the intramolecular hydrogen bonding could be useful in generating antibody-based catalysts that recruit cocaine to the boat conformation and an analog that elicited antibodies to approximate the protonated tropane N and the benzoyl O more closely than the natural boat conformer might increase the contribution from hydrogen bonding. Such a stable analog of the transition state for intramolecular catalysis of cocaine benzoyl-ester hydrolysis was synthesized and used to successfully elicit a number of anti-cocaine catalytic antibodies.