This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.106.087130v1 92/12/4157    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Sonne, J. Articles by Peters, G. H. Search for Related Content PubMed PubMed Citation Articles by Sonne, J. Articles by Peters, G. H.

Reparameterization of All-Atom Dipalmitoylphosphatidylcholine Lipid Parameters Enables Simulation of Fluid Bilayers at Zero Tension

Jacob Sonne ^* , Morten Ø. Jensen , Flemming Y. Hansen ^*, Lars Hemmingsen  and Günther H. Peters ^* 

^* Department of Chemistry, Technical University of Denmark, Lyngby, Denmark; MEMPHYS-Center for Biomembrane Physics, University of Southern Denmark, Odense, Denmark; Department of Life Sciences and Chemistry, Roskilde University, Roskilde, Denmark; and Department of Natural Sciences, The Faculty of Life Sciences, University of Copenhagen, Copenhagen, Denmark

Correspondence: Address reprint requests to Günther H. Peters, E-mail: ghp{at}kemi.dtu.dk.

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of 48 Ų. To obtain fluid (L) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 ± 0.1 Ų. Compared to the 48 Ų, the new value of 60.4 Ų is in fair agreement with the experimental value of 64 Ų. In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.