Simulations of Lipid Bilayers Using the CHARMM36 Force Field with the TIP3P-FB and TIP4P-FB Water Models

Overview

Journal PeerJ

Specialties Biology
Environmental Health
General Medicine

Date 2018 Aug 22

PMID 30128211

Citations 14

Authors

Fatima Sajadi

Christopher N Rowley

Affiliations

Soon will be listed here.

Abstract

The CHARMM36 force field for lipids is widely used in simulations of lipid bilayers. The CHARMM family of force fields were developed for use with the mTIP3P water model. This water model has an anomalously high dielectric constant and low viscosity, which limits its accuracy in the calculation of quantities like permeability coefficients. The TIP3P-FB and TIP4P-FB water models are more accurate in terms of the dielectric constant and transport properties, which could allow more accurate simulations of systems containing water and lipids. To test whether the CHARMM36 lipid force field is compatible with the TIP3P-FB and TIP4P-FB water models, we have performed simulations of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. The calculated headgroup area, compressibility, order parameters, and X-ray form factors are in good agreement with the experimental values, indicating that these improved water models can be used with the CHARMM36 lipid force field without modification when calculating membrane physical properties. The water permeability predicted by these models is significantly different; the mTIP3P-model diffusion in solution and at the lipid-water interface is anomalously fast due to the spuriously low viscosity of mTIP3P-model water, but the potential of mean force of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials. As a result, the rates of water permeation calculated the FB water models are slower than the experimental value by a factor of 15-17, while simulations with the mTIP3P model only underestimate the water permeability by a factor of 3.

Citing Articles

Exploring Marine Natural Compounds: Innovative Therapeutic Candidates Against Chagas Disease Through Virtual Screening and Molecular Dynamics.

Maya-Ramirez C, Saih A, Mendez Tenorio A, Wong Baeza C, Nogueda Torres B, Santiago Hernandez J Life (Basel). 2025; 15(2).

PMID: 40003601 PMC: 11856606. DOI: 10.3390/life15020192.

Insights into Dermal Permeation of Skin Oil Oxidation Products from Enhanced Sampling Molecular Dynamics Simulation.

Thomas R, Prabhakar P, Tobias D, von Domaros M J Phys Chem B. 2025; 129(6):1784-1794.

PMID: 39901666 PMC: 11831647. DOI: 10.1021/acs.jpcb.4c08090.

Impacts of Nutlin-3a and exercise on murine double minute 2-enriched glioma treatment.

Chen Y, Fan Z, Luo Z, Kang X, Wan R, Li F Neural Regen Res. 2024; 20(4):1135-1152.

PMID: 38989952 PMC: 11438351. DOI: 10.4103/NRR.NRR-D-23-00875.

Cryptotanshinone-Induced Permeabilization of Model Phospholipid Membranes: A Biophysical Study.

Ortiz J, Aranda F, Teruel J, Ortiz A Membranes (Basel). 2024; 14(6).

PMID: 38921485 PMC: 11205401. DOI: 10.3390/membranes14060118.

Identification of a potent and selective LAPTc inhibitor by RapidFire-Mass Spectrometry, with antichagasic activity.

Izquierdo M, Lin D, ONeill S, Webster L, Paterson C, Thomas J PLoS Negl Trop Dis. 2024; 18(2):e0011956.

PMID: 38359089 PMC: 10901353. DOI: 10.1371/journal.pntd.0011956.

References

Jambeck J, Lyubartsev A . Derivation and systematic validation of a refined all-atom force field for phosphatidylcholine lipids. J Phys Chem B. 2012; 116(10):3164-79. PMC: 3320744. DOI: 10.1021/jp212503e. View

Braun D, Boresch S, Steinhauser O . Transport and dielectric properties of water and the influence of coarse-graining: comparing BMW, SPC/E, and TIP3P models. J Chem Phys. 2014; 140(6):064107. DOI: 10.1063/1.4864117. View

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

Hermans J, Pathiaseril A, Anderson A . Excess free energy of liquids from molecular dynamics simulations. Application to water models. J Am Chem Soc. 2011; 110(18):5982-6. DOI: 10.1021/ja00226a009. View

Gaalswyk K, Awoonor-Williams E, Rowley C . Generalized Langevin Methods for Calculating Transmembrane Diffusivity. J Chem Theory Comput. 2016; 12(11):5609-5619. DOI: 10.1021/acs.jctc.6b00747. View

Wu E, Cheng X, Jo S, Rui H, Song K, Davila-Contreras E . CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. J Comput Chem. 2014; 35(27):1997-2004. PMC: 4165794. DOI: 10.1002/jcc.23702. View

Jo S, Kim T, Iyer V, Im W . CHARMM-GUI: a web-based graphical user interface for CHARMM. J Comput Chem. 2008; 29(11):1859-65. DOI: 10.1002/jcc.20945. View

Mendes Ferreira T, Coreta-Gomes F, Ollila O, Moreno M, Vaz W, Topgaard D . Cholesterol and POPC segmental order parameters in lipid membranes: solid state 1H-13C NMR and MD simulation studies. Phys Chem Chem Phys. 2012; 15(6):1976-89. DOI: 10.1039/c2cp42738a. View

Mathai J, Tristram-Nagle S, Nagle J, Zeidel M . Structural determinants of water permeability through the lipid membrane. J Gen Physiol. 2008; 131(1):69-76. PMC: 2174160. DOI: 10.1085/jgp.200709848. View

10.

Guixa-Gonzalez R, Rodriguez-Espigares I, Ramirez-Anguita J, Carrio-Gaspar P, Martinez-Seara H, Giorgino T . MEMBPLUGIN: studying membrane complexity in VMD. Bioinformatics. 2014; 30(10):1478-80. DOI: 10.1093/bioinformatics/btu037. View

11.

Venable R, Brown F, Pastor R . Mechanical properties of lipid bilayers from molecular dynamics simulation. Chem Phys Lipids. 2015; 192:60-74. PMC: 4684433. DOI: 10.1016/j.chemphyslip.2015.07.014. View

12.

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

13.

Kucerka N, Tristram-Nagle S, Nagle J . Closer look at structure of fully hydrated fluid phase DPPC bilayers. Biophys J. 2006; 90(11):L83-5. PMC: 1459514. DOI: 10.1529/biophysj.106.086017. View

14.

Lee J, Cheng X, Swails J, Yeom M, Eastman P, Lemkul J . CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J Chem Theory Comput. 2015; 12(1):405-13. PMC: 4712441. DOI: 10.1021/acs.jctc.5b00935. View

15.

Douliez J, Leonard A, Dufourc E . Restatement of order parameters in biomembranes: calculation of C-C bond order parameters from C-D quadrupolar splittings. Biophys J. 1995; 68(5):1727-39. PMC: 1282076. DOI: 10.1016/S0006-3495(95)80350-4. View

16.

Kucerka N, Tristram-Nagle S, Nagle J . Structure of fully hydrated fluid phase lipid bilayers with monounsaturated chains. J Membr Biol. 2006; 208(3):193-202. DOI: 10.1007/s00232-005-7006-8. View

17.

Benz R, Castro-Roman F, Tobias D, White S . Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach. Biophys J. 2004; 88(2):805-17. PMC: 1305157. DOI: 10.1529/biophysj.104.046821. View

18.

Shirts M, Pande V . Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. J Chem Phys. 2005; 122(13):134508. DOI: 10.1063/1.1877132. View

19.

Lee C, Comer J, Herndon C, Leung N, Pavlova A, Swift R . Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds. J Chem Inf Model. 2016; 56(4):721-33. PMC: 5280572. DOI: 10.1021/acs.jcim.6b00022. View

20.

Ballal D, Venkataraman P, Fouad W, Cox K, Chapman W . Isolating the non-polar contributions to the intermolecular potential for water-alkane interactions. J Chem Phys. 2014; 141(6):064905. DOI: 10.1063/1.4892341. View