Two-bead polarizable water models combined with a two-bead multipole force field (TMFF) for coarse-grained simulation of proteins

2017 ◽  
Vol 19 (10) ◽  
pp. 7410-7419 ◽  
Author(s):  
Min Li ◽  
John Z. H. Zhang

(a) Four water molecules contained in the polarizable CG water models in (b) two-bead polarizable water model 1 (TPW1) and (c) two-bead polarizable water model 2 (TPW2).

Author(s):  
Min Li ◽  
John Z.H. Zhang

The development of coarse-grained (CG) water model is increasingly important in CG studies of biological processes. In this work, we developed a generic CG force field of liquid water on...


2017 ◽  
Vol 146 (5) ◽  
pp. 054501 ◽  
Author(s):  
Julian Michalowsky ◽  
Lars V. Schäfer ◽  
Christian Holm ◽  
Jens Smiatek

2018 ◽  
Author(s):  
Fatima Sajadi ◽  
Christopher Rowley

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 TIP3P 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 DPPC and POPC 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. The water permeability predicted by these models is significantly different; the TIP3P-model diffusion in solution and at the lipid--water interface is anomalously fast due to the spuriously low viscosity of TIP3P-model water, but the PMF of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials.


2011 ◽  
Vol 110-116 ◽  
pp. 3883-3887
Author(s):  
N. Maftouni ◽  
M. Amininassab ◽  
M. N. Mello ◽  
S. Marink

It is very essential to know mechanical properties in different regions of nanobio membrane as one of the most important parts of living systems. Here the coarse-grained (CG) simulations method have been used to study the pressure profile in a system including nanobio membrane and water. CG simulations have become an important tool to study many biomolecular processes, exploring scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization has been ignored in the initial CG water molecules model. Given the important role of water as a solvent in biological systems, its treatment is very important to the properties derived from simulation studies. Till now two models have been parameterized to simulate water: i) standard MARTINI water and ii) polarizable coarse-grained water model. Both of mentioned water models are proper to be used in combination with the CG MARTINI force field. In this work both of these models have been used for simulation. One micro second CG molecular dynamics simulation has been done for two separate systems. Each system includes water and hydrated 1-palmitoyl-2-oleoyl-1-sn-3-phosphatidylcholine (POPC) lipid nanobio membrane. The difference between two systems is in simulated water models that one system has standard MARTINI water and the other one has polarizable water. In each case pressure profile calculation has been done via Virial pressure theorem. Results indicate that using polarizable water model leads to higher picks in pressure profile in water region near surface of nanobio membrane. This can be related to density of polarizable water and also may play role as a small barrier.


2010 ◽  
Vol 6 (6) ◽  
pp. e1000810 ◽  
Author(s):  
Semen O. Yesylevskyy ◽  
Lars V. Schäfer ◽  
Durba Sengupta ◽  
Siewert J. Marrink

2019 ◽  
Vol 16 (1) ◽  
pp. 782-793 ◽  
Author(s):  
Yusuke Miyazaki ◽  
Susumu Okazaki ◽  
Wataru Shinoda

Author(s):  
Fatima Sajadi ◽  
Christopher Rowley

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 TIP3P 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 DPPC and POPC 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. The water permeability predicted by these models is significantly different; the TIP3P-model diffusion in solution and at the lipid--water interface is anomalously fast due to the spuriously low viscosity of TIP3P-model water, but the PMF of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials.


PeerJ ◽  
2018 ◽  
Vol 6 ◽  
pp. e5472 ◽  
Author(s):  
Fatima Sajadi ◽  
Christopher N. Rowley

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.


Sign in / Sign up

Export Citation Format

Share Document