Polarizable Water Model for the Coarse-Grained MARTINI Force Field
Top Cited Papers
Open Access
- 10 June 2010
- journal article
- research article
- Published by Public Library of Science (PLoS) in PLoS Computational Biology
- Vol. 6 (6), e1000810
- https://doi.org/10.1371/journal.pcbi.1000810
Abstract
Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial 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, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models. Many biomolecular processes involve charged species moving between regions of high polarity, such as the water phase, and regions of lower polarity, such as the lipid membrane. Due to the change in electrostatic screening between these two environments, the strength of the interactions between the moving charge and the surrounding molecules also changes. This has important consequences for the way biological activity is controlled. To help understand the forces driving the movement of biomolecules, we developed a computational model which is capable of describing these processes at near-atomic detail. To do so efficiently, we use a coarse-grained description of the molecules, in which some of the atomistic detail is averaged out. To capture the inhomogeneous nature of the dielectric response, we re-introduce some detail in the water model; the new model effectively mimics the orientational polarizability of real water molecules, and screens electrostatic interactions realistically. This enables the study of a number of important biological processes that were hitherto considered challenging for coarse-grained models, such as the permeation of ions across a lipid membrane and the rupture of membranes due to an electrostatic field, at relatively low computational cost.Keywords
This publication has 58 references indexed in Scilit:
- Current Status of the AMOEBA Polarizable Force FieldThe Journal of Physical Chemistry B, 2010
- Multiscale simulations of protein landscapes: Using coarse‐grained models as reference potentials to full explicit modelsProteins: Structure, Function, and Bioinformatics, 2009
- Free Energy for the Permeation of Na+ and Cl− Ions and Their Ion-Pair through a Zwitterionic Dimyristoyl Phosphatidylcholine Lipid Bilayer by Umbrella Integration with Harmonic Fourier BeadsJournal of the American Chemical Society, 2009
- Coarse-Grained Simulations of the Membrane-Active Antimicrobial Peptide Maculatin 1.1Biophysical Journal, 2008
- Coarse-Grained Molecular Dynamics Studies of the Concentration and Size Dependence of Fifth- and Seventh-Generation PAMAM Dendrimers on Pore Formation in DMPC BilayerThe Journal of Physical Chemistry B, 2008
- GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular SimulationJournal of Chemical Theory and Computation, 2008
- GROMACS: Fast, flexible, and freeJournal of Computational Chemistry, 2005
- A Multiscale Coarse-Graining Method for Biomolecular SystemsThe Journal of Physical Chemistry B, 2005
- VMD: Visual molecular dynamicsJournal of Molecular Graphics, 1996
- Molecular dynamics with coupling to an external bathThe Journal of Chemical Physics, 1984