Author(s): Klauda JB, Brooks BR, MacKerell AD Jr, Venable RM, Pastor RW
Abstract Share this page
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.
This article was published in J Phys Chem B
and referenced in Journal of Computer Science & Systems Biology