Lipid Configurations from Molecular Dynamics Simulations

Area of Science:

• Computational Biophysics
• Molecular Dynamics Simulations
• Lipid Bilayer Modeling

Background:

• Accurate force fields are essential for simulating biomembrane conformational properties.
• Understanding lipid behavior in bilayers is key to deciphering membrane function.
• Existing force fields' ability to replicate phospholipid crystal structures is debated.

Purpose of the Study:

• To evaluate the accuracy of two prominent lipid force fields (Berger and CHARMM36) in reproducing lipid conformational properties.
• To compare simulation results with established structural principles of phospholipids.
• To identify force field-dependent artifacts in biomembrane simulations.

Main Methods:

• Molecular dynamics simulations of palmitoyl-oleoyl phosphatidylcholine in hydrated fluid bilayers.
• Analysis of dihedral angle distributions for lipid chains and headgroups.
• Comparison of simulation outputs with experimental data and crystal structures.

Main Results:

• CHARMM36 accurately captures lipid chain inequivalence in fluid bilayers, unlike the Berger force field.
• Berger force field exhibits artifactual Na+ ion binding due to exposed carbonyls.
• CHARMM36 better reproduces NMR data and headgroup conformations across various lipid types compared to Berger.

Conclusions:

• CHARMM36 is superior for simulating lipid conformational properties in fluid bilayers.
• The Berger force field introduces artifacts, limiting its reliability for biomembrane simulations.
• Further experimental validation of choline headgroup conformations is warranted.