J. Phys. Chem. B, 109 (11), 5300 -5311, 2005. 10.1021/jp0468096 S1089-5647(04)06809-9
Web Release Date: February 24, 2005

Copyright © 2005 American Chemical Society

An ab Initio Study on the Torsional Surface of Alkanes and Its Effect on Molecular Simulations of Alkanes and a DPPC Bilayer

Jeffery B. Klauda and Bernard R. Brooks

Laboratory of Biophysical Chemistry, National Institutes of Health, Building 50, 50 South Drive, Bethesda, Maryland 20892

Alexander D. MacKerell, Jr.

Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201

Richard M. Venable and Richard W. Pastor*

Laboratory of Biophysics, Center for Biologics Evaluation and Research, FDA, 1401 Rockville Pike, Rockville, Maryland 20852-1448

Received: July 19, 2004

In Final Form: November 22, 2004

Abstract:

Energies of 119 conformations of normal alkanes from butane to heptane were calculated at approximately the CCSD(T)/cc-pVQZ level. Energies of gauche (g) conformers relative to trans (t) decrease as chain length increases. In what is termed the "positive pentane effect", adjacent gauche conformers of the same sign are stabilized compared to nonadjacent conformers; e.g., for hexane the energies of tgt, tgg, and gtg are 0.600, 0.930, and 1.18 kcal/mol, respectively. Torsional terms in the CHARMM27 (C27) force field were fit to the calculated QM energies to yield a revised potential, C27r. Molecular dynamics simulations of normal alkanes (heptane, decane, tridecane, and pentadecane) with C27r yield higher populations of gauche states, increased transition rates, and improved agreement with experiment as compared to C27. In addition, C27r simulations of a hydrated DPPC lipid bilayer yield improved agreement with the experimental NMR deuterium order parameters for the aliphatic chain ends.


Download the full text: PDF | HTML