Performance of the general amber force field in modeling aqueous POPC membrane bilayers

Balázs Jójárt, Tamás A. Martinek

Research output: Article

97 Citations (Scopus)


The aim of this work was to answer the question of whether the general amber force field (GAFF) is good enough to simulate fully hydrated POPC membrane bilayers. The test system contained 128 POPC and 2985 TIP3P water molecules. The equilibration was carried out in a nonarbitrary manner to reach the stable liquid-crystalline phase. The simulations were performed by using particle mesh Ewald electrostatics implemented in molecular dynamics packages Amber8 (NPT ensembles) and NAMD2 (NPγT ensembles). The computational results were assessed against the following experimental membrane properties: (i) area per lipid, (ii) area compressibility modulus, (iii) order parameter, (iv) gauche conformations per acyl chain, (v) lateral diffusion coefficients, (vi) electron density profile, and (vii) bound water at the lipid/water interface. The analyses revealed that the tested force field combination approximates the experimental values at an unexpectedly good level when the NPγT ensemble is applied with a surface tension of 60 mN m-1 per bilayer. It is concluded that the GAFF/TIP3P combination can be utilized for aqueous membrane bilayer simulations, as it provides acceptable accuracy for biomolecular modeling.

Original languageEnglish
Pages (from-to)2051-2058
Number of pages8
JournalJournal of Computational Chemistry
Issue number12
Publication statusPublished - szept. 1 2007


ASJC Scopus subject areas

  • Chemistry(all)
  • Computational Mathematics

Cite this