Computation of the Electronic and Spectroscopic Properties of Carbohydrates Using Novel Density Functional and Vibrational Self-Consistent Field Methods

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.

Applications of terahertz spectroscopy in biosystems

Terahertz (THz) spectroscopic investigations of condensed-phase biological samples are reviewed ranging from the simple crystalline forms of amino acids, carbohydrates and polypeptides to the more complex aqueous forms of small proteins, DNA and RNA. Vibrationally resolved studies of crystalline samples have revealed the exquisite sensitivity of THz modes to crystalline order, temperature, conformational form, peptide sequence and local solvate environment and have given unprecedented measures of the binding force constants and anharmonic character of the force fields, properties necessary to improve predictability but not readily obtainable using any other method. These studies have provided benchmark vibrational data on extended periodic structures for direct comparisons with classical (CHARMm) and quantum chemical (density functional theory) theories. For the larger amorphous and/or aqueous phase samples, the THz modes form a continuum-like absorption that arises because of the full accessibility to conformational space and/or the rapid time scale for inter-conversion in these environments. Despite severe absorption by liquid water, detailed investigations have uncovered the photo- and hydration-induced conformational flexibility of proteins, the solvent shell depth of the water/biomolecule boundary layers and the solvent reorientation dynamics occurring in these interfacial layers that occur on sub-picosecond time scales. As such, THz spectroscopy has enhanced and extended the accessibility to intermolecular forces, length- and timescales important in biological structure and activity.

An ab initio study of S(N)2 reactivity at C6 in hexopyranose derivatives. II. Role of populations, barriers, and reaction path curvature

This paper continues our investigation into a simple dipole-dipole interaction model proposed to explain the dramatically reduced S(N)2 reactivity at the primary C6 position of galacto-configured pyranose systems relative to their gluco isomers. We present ab initio calculations (B3LYP/6-31+G(d,p)) on six model structures that show that this effect is not a major influence. Reactant rotameric equilibria as well as free-energy reaction barriers and reaction path curvature were evaluated. Results point to a number of other factors that could account for the observed reactivity differences. Our results cast doubt on the general relevance of transition structure dipole-dipole repulsions to S(N)2 reactivity.

Evaluation of carbohydrate molecular mechanical force fields by quantum mechanical calculations

A quantitative evaluation of 20 second-generation carbohydrate force fields was carried out using ab initio and density functional methods. Geometry-optimized structures (B3LYP/6-31G(d)) and relative energies using augmented correlation consistent basis sets were calculated in gas phase for monosaccharide carbohydrate benchmark systems. Selected results are: (i). The interaction energy of the alpha-d-glucopyranose.H(2)O heterodimer is estimated to be 4.9 kcal/mol, using a composite method including terms at highly correlated (CCSD(T)) level. Most molecular mechanics force fields are in error in this respect; (ii). The (3)E envelope (south) pseudorotational conformer of methyl 5-deoxy-beta-d-xylofuranoside is 0.66 kcal/mol more stable than the (3)E envelope (north) conformer and the alpha-anomer of methyl d-glucopyranoside is 0.82 kcal/mol more stable than the beta-anomer; (iii). The relative energies of the (gg, gt and tg) rotamers of methyl alpha-d-glucopyranoside and methyl alpha-d-galactopyranoside are (0.13, 0.00, 0.15) and (0.64, 0.00, 0.77) kcal/mol, respectively. The results of the quantum mechanical calculations are compared with the results of calculations using the 20 second-generation carbohydrate force fields. No single force field is consistently better than the others for all the test cases. A statistical assessment of the performance of the force fields indicates that CHEAT(95), CFF, certain versions of Amber and of MM3 have the best overall performance, for these gas phase monosaccharide systems.

Intramolecular interaction energies in model alanine and glycine tetrapeptides. Evaluation of anisotropy, polarization, and correlation effects. A parallel ab initio HF/MP2, DFT, and polarizable molecular mechanics study

An extension of the SIBFA polarizable molecular mechanics procedure to flexible oligopeptides is reported. The procedure is evaluated by computing the relative conformational energies, deltaE(conf), of the alanine tetrapeptide in 10 representative conformations, which were originally derived by Beachy et al. (J Am Chem Soc 1997, 119, 5908) to benchmark molecular mechanics procedures with respect to ab initio computations. In the present study, a particular emphasis is on the separable nature of the components of the energy and the particular impact of the polarization energy component on deltaE(conf). We perform comparisons with respect to single-point HF, DFT, LMP2, and MP2 computations done at the SIBFA-derived energy minima. Such comparisons are made first for the 10 conformers derived from phi/psi torsional angle energy-minimization (the rigid rotor approach), and, in a second step, after allowing additional relaxation of the C(alpha) centered valence angles. In both series of energy-minimization, the SIBFA deltaE(conf) compared best with the LMP2 results using the 6-311G** basis set, the rms being 1.3 kcal/mol. In the absence of the polarization component, the rms is 3.5 kcal/mol. In both series of minimizations, the magnitudes of deltaE(conf), computed as differences with respect to the most stable conformer taken as energy zero, decrease along the series: HF > DFT > LMP2 > SIBFA > MP2, indicative of increasing stabilization of the most highly folded conformers.