CHARMM36 Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Phosphatidylethanolamine, Phosphatidylglycerol, and Ether Lipids

Overview

Journal J Chem Theory Comput

Publisher American Chemical Society

Specialties Biochemistry
Chemistry

Date 2021 Feb 23

PMID 33620194

Citations 27

Authors

Yalun Yu

Andreas Kramer

Richard M Venable

Bernard R Brooks

Jeffery B Klauda

Richard W Pastor

Affiliations

Soon will be listed here.

Abstract

Long-range Lennard-Jones (LJ) interactions have been incorporated into the CHARMM36 (C36) lipid force field (FF) using the LJ particle-mesh Ewald (LJ-PME) method in order to remove the inconsistency of bilayer and monolayer properties arising from the exclusion of long-range dispersion [Yu, Y.; Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. 2021, 10.1021/acs.jctc.0c01326. (preceding article in this issue)]. The new FF is denoted C36/LJ-PME. While the first optimization was based on three phosphatidylcholines (PCs), this work extends the validation and parametrization to more lipids including PC, phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and ether lipids. The agreement with experimental structure data is excellent for PC, PE, and ether lipids. C36/LJ-PME also compares favorably with scattering data of PG bilayers but less so with NMR deuterium order parameters of 1,2-dimyristoyl--glycero-3-phospho-(1'-rac-glycerol) (DMPG) at 303.15 K, indicating a need for future optimization regarding PG-specific parameters. Frequency dependence of NMR spin-lattice relaxation times is well-described by C36/LJ-PME, and the overall agreement with experiment is comparable to C36. Lipid diffusion is slower than C36 due to the added long-range dispersion causing a higher viscosity, although it is still too fast compared to experiment after correction for periodic boundary conditions. When using a 10 Å real-space cutoff, the simulation speed of C36/LJ-PME is roughly equal to C36. While more lipids will be incorporated into the FF in the future, C36/LJ-PME can be readily used for common lipids and extends the capability of the CHARMM FF by supporting monolayers and eliminating the cutoff dependence.

Citing Articles

β-Barrel proteins dictate the effect of core oligosaccharide composition on outer membrane mechanics.

Fitzmaurice D, Amador A, Starr T, Hocky G, Rojas E Biophys J. 2025; 124(5):765-777.

PMID: 39863924 PMC: 11897544. DOI: 10.1016/j.bpj.2025.01.017.

Computationally Efficient Polarizable MD Simulations: A Simple Water Model for the Classical Drude Oscillator Polarizable Force Field.

Teng X, Yu W, MacKerell Jr A J Phys Chem Lett. 2025; 16(4):1016-1023.

PMID: 39841123 PMC: 11895599. DOI: 10.1021/acs.jpclett.4c03451.

Cell membranes sustain phospholipid imbalance via cholesterol asymmetry.

Doktorova M, Symons J, Zhang X, Wang H, Schlegel J, Lorent J bioRxiv. 2024; .

PMID: 39713443 PMC: 11661119. DOI: 10.1101/2023.07.30.551157.

Atomistic Simulations and Analysis of Peripheral Membrane Proteins with Model Lipid Bilayers.

Campbell O, Allsopp R, Klauda J, Monje V Methods Mol Biol. 2024; 2888:281-303.

PMID: 39699738 DOI: 10.1007/978-1-0716-4318-1_19.

Rotary mechanism of the prokaryotic V motor driven by proton motive force.

Kishikawa J, Nishida Y, Nakano A, Kato T, Mitsuoka K, Okazaki K Nat Commun. 2024; 15(1):9883.

PMID: 39567487 PMC: 11579504. DOI: 10.1038/s41467-024-53504-x.

References

Yu Y, Klauda J . Update of the CHARMM36 United Atom Chain Model for Hydrocarbons and Phospholipids. J Phys Chem B. 2020; 124(31):6797-6812. DOI: 10.1021/acs.jpcb.0c04795. View

Venable R, Ingolfsson H, Lerner M, Perrin Jr B, Camley B, Marrink S . Lipid and Peptide Diffusion in Bilayers: The Saffman-Delbrück Model and Periodic Boundary Conditions. J Phys Chem B. 2016; 121(15):3443-3457. PMC: 6326097. DOI: 10.1021/acs.jpcb.6b09111. View

Pan J, Tristram-Nagle S, Nagle J . Effect of cholesterol on structural and mechanical properties of membranes depends on lipid chain saturation. Phys Rev E Stat Nonlin Soft Matter Phys. 2009; 80(2 Pt 1):021931. PMC: 2756665. DOI: 10.1103/PhysRevE.80.021931. View

Klauda J, Brooks B, MacKerell Jr A, Venable R, Pastor R . An ab initio study on the torsional surface of alkanes and its effect on molecular simulations of alkanes and a DPPC bilayer. J Phys Chem B. 2006; 109(11):5300-11. DOI: 10.1021/jp0468096. View

Shaikh S, Brzustowicz M, Gustafson N, Stillwell W, Wassall S . Monounsaturated PE does not phase-separate from the lipid raft molecules sphingomyelin and cholesterol: role for polyunsaturation?. Biochemistry. 2002; 41(34):10593-602. DOI: 10.1021/bi025712b. View

Nagle J, Tristram-Nagle S . Structure of lipid bilayers. Biochim Biophys Acta. 2000; 1469(3):159-95. PMC: 2747654. DOI: 10.1016/s0304-4157(00)00016-2. View

Kucerka N, Nieh M, Katsaras J . Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature. Biochim Biophys Acta. 2011; 1808(11):2761-71. DOI: 10.1016/j.bbamem.2011.07.022. View

Venable R, Luo Y, Gawrisch K, Roux B, Pastor R . Simulations of anionic lipid membranes: development of interaction-specific ion parameters and validation using NMR data. J Phys Chem B. 2013; 117(35):10183-92. PMC: 3813009. DOI: 10.1021/jp401512z. View

Kucerka N, Liu Y, Chu N, Petrache H, Tristram-Nagle S, Nagle J . Structure of fully hydrated fluid phase DMPC and DLPC lipid bilayers using X-ray scattering from oriented multilamellar arrays and from unilamellar vesicles. Biophys J. 2005; 88(4):2626-37. PMC: 1305359. DOI: 10.1529/biophysj.104.056606. View

10.

Kucerka N, Katsaras J, Nagle J . Comparing membrane simulations to scattering experiments: introducing the SIMtoEXP software. J Membr Biol. 2010; 235(1):43-50. PMC: 2876336. DOI: 10.1007/s00232-010-9254-5. View

11.

Venable R, Kramer A, Pastor R . Molecular Dynamics Simulations of Membrane Permeability. Chem Rev. 2019; 119(9):5954-5997. PMC: 6506413. DOI: 10.1021/acs.chemrev.8b00486. View

12.

Zhuang X, Makover J, Im W, Klauda J . A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties. Biochim Biophys Acta. 2014; 1838(10):2520-9. DOI: 10.1016/j.bbamem.2014.06.010. View

13.

Camley B, Lerner M, Pastor R, Brown F . Strong influence of periodic boundary conditions on lateral diffusion in lipid bilayer membranes. J Chem Phys. 2016; 143(24):243113. PMC: 4608966. DOI: 10.1063/1.4932980. View

14.

Wennberg C, Murtola T, Hess B, Lindahl E . Lennard-Jones Lattice Summation in Bilayer Simulations Has Critical Effects on Surface Tension and Lipid Properties. J Chem Theory Comput. 2015; 9(8):3527-37. DOI: 10.1021/ct400140n. View

15.

Botan A, Favela-Rosales F, Fuchs P, Javanainen M, Kanduc M, Kulig W . Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions. J Phys Chem B. 2015; 119(49):15075-88. PMC: 4677354. DOI: 10.1021/acs.jpcb.5b04878. View

16.

Henin J, Shinoda W, Klein M . United-atom acyl chains for CHARMM phospholipids. J Phys Chem B. 2008; 112(23):7008-15. PMC: 8451178. DOI: 10.1021/jp800687p. View

17.

Pan J, Tristram-Nagle S, Kucerka N, Nagle J . Temperature dependence of structure, bending rigidity, and bilayer interactions of dioleoylphosphatidylcholine bilayers. Biophys J. 2007; 94(1):117-24. PMC: 2134881. DOI: 10.1529/biophysj.107.115691. View

18.

Sundaralihgam M . DISCUSSION PAPER: MOLECULAR STRUCTURES AND CONFORMATIONS OF THE PHOSPHOLIPIDS AND SPHINGOMYELINS. Ann N Y Acad Sci. 2017; 195(1):324-355. DOI: 10.1111/j.1749-6632.1972.tb54814.x. View

19.

Pan J, Marquardt D, Heberle F, Kucerka N, Katsaras J . Revisiting the bilayer structures of fluid phase phosphatidylglycerol lipids: Accounting for exchangeable hydrogens. Biochim Biophys Acta. 2014; 1838(11):2966-9. DOI: 10.1016/j.bbamem.2014.08.009. View

20.

Somerharju P, Virtanen J, Eklund K, Vainio P, Kinnunen P . 1-Palmitoyl-2-pyrenedecanoyl glycerophospholipids as membrane probes: evidence for regular distribution in liquid-crystalline phosphatidylcholine bilayers. Biochemistry. 1985; 24(11):2773-81. DOI: 10.1021/bi00332a027. View