jp401718k_si_001.pdf (858.02 kB)
Comparing Simulations of Lipid Bilayers to Scattering Data: The GROMOS 43A1-S3 Force Field
journal contributionposted on 2013-05-02, 00:00 authored by Anthony R. Braun, Jonathan N. Sachs, John F. Nagle
Simulations of DOPC at T = 303 K were performed using the united atom force field 43A1-S3 at six fixed projected areas, AP = 62, 64, 66, 68, 70, and 72 Å2, as well as a tensionless simulation that produced an average ANPT = 65.8 Å2. After a small undulation correction for the system size consisting of 288 lipids, results were compared to experimental data. The best, and excellent, fit to neutron scattering data occurs at an interpolated AN = 66.6 Å2 and the best, but not as good, fit to the more extensive X-ray scattering data occurs at AX = 68.7 Å2. The distance ΔDB–H between the Gibbs dividing surface for water and the peak in the electron density profile agrees with scattering experiments. The calculated area compressibility KA = 277 ± 10 mN/m is in excellent agreement with the micromechanical experiment. The volume per lipid VL is smaller than volume experiments which suggests a workaround that raises all the areas by about 1.5%. Although AX ≠ AN ≠ ANPT, this force field obtains acceptable agreement with experiment for AL = 67.5 Å2 (68.5 Å2 in the workaround), which we suggest is a better DOPC result from 43A1-S3 simulations than its value from the tensionless NPT simulation. However, nonsimulation modeling obtains better simultaneous fits to both kinds of scattering data, which suggests that the force fields can still be improved.
area compressibility KAvolume experimentsDOPC result288 lipidssystem sizeatom force field 43force fieldsScattering DataworkaroundANPT68.7 Å2.force field obtains303 Kmicromechanical experimentnonsimulation modeling obtains65.8 Å2.GROMOS 43NPT simulationDBLipid Bilayers66.6 Å2lipid VLundulation correctionALAPAXdataelectron density profile