1
|
Klein F, Soñora M, Helene Santos L, Nazareno Frigini E, Ballesteros-Casallas A, Rodrigo Machado M, Pantano S. The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels. J Struct Biol 2023; 215:107985. [PMID: 37331570 DOI: 10.1016/j.jsb.2023.107985] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 05/18/2023] [Accepted: 06/13/2023] [Indexed: 06/20/2023]
Abstract
The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.
Collapse
Affiliation(s)
- Florencia Klein
- Laboratoire de Biochimie Théorique, UPR9080, CNRS, Paris, France
| | - Martín Soñora
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay
| | | | - Ezequiel Nazareno Frigini
- Instituto Multidisciplinario de Investigaciones Biológicas de San Luis (IMIBIO-SL), Universidad Nacional de San Luis - CONICET, San Luis, Argentina
| | - Andrés Ballesteros-Casallas
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay; Area Bioinformática, DETEMA, Facultad de Química, Universidad de la República, General Flores 2124, Montevideo, 11600, Uruguay
| | | | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay; Area Bioinformática, DETEMA, Facultad de Química, Universidad de la República, General Flores 2124, Montevideo, 11600, Uruguay.
| |
Collapse
|
2
|
Dhabal D, Sankaranarayanan SKRS, Molinero V. Stability and Metastability of Liquid Water in a Machine-Learned Coarse-Grained Model with Short-Range Interactions. J Phys Chem B 2022; 126:9881-9892. [PMID: 36383428 DOI: 10.1021/acs.jpcb.2c06246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Coarse-grained water models are ∼100 times more efficient than all-atom models, enabling simulations of supercooled water and crystallization. The machine-learned monatomic model ML-BOP reproduces the experimental equation of state (EOS) and ice-liquid thermodynamics at 0.1 MPa on par with the all-atom TIP4P/2005 and TIP4P/Ice models. These all-atom models were parametrized using high-pressure experimental data and are either accurate for water's EOS (TIP4P/2005) or ice-liquid equilibrium (TIP4P/Ice). ML-BOP was parametrized from temperature-dependent ice and liquid experimental densities and melting data at 0.1 MPa; its only pressure training is from compression of TIP4P/2005 ice at 0 K. Here we investigate whether ML-BOP replicates the experimental EOS and ice-water thermodynamics along all pressures of ice I. We find that ML-BOP reproduces the temperature, enthalpy, entropy, and volume of melting of hexagonal ice up to 400 MPa and the EOS of water along the melting line with an accuracy that rivals that of both TIP4P/2005 and TIP4P/Ice. We interpret that the accuracy of ML-BOP originates from its ability to capture the shift between compact and open local structures to changes in pressure and temperature. ML-BOP reproduces the sharpening of the tetrahedral peak of the pair distribution function of water upon supercooling, and its pressure dependence. We characterize the region of metastability of liquid ML-BOP with respect to crystallization and cavitation. The accessibility of ice crystallization to simulations of ML-BOP, together with its accurate representation of the thermodynamics of water, makes it promising for investigating the interplay between anomalies, glass transition, and crystallization under conditions challenging to access through experiments.
Collapse
Affiliation(s)
- Debdas Dhabal
- Department of Chemistry, The University of Utah, Salt Lake City, Utah84112-0850, United States
| | - Subramanian K R S Sankaranarayanan
- Department of Mechanical and Industrial Engineering, University of Illinois, Chicago, Illinois60607, United States.,Center for Nanoscale Materials, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, Salt Lake City, Utah84112-0850, United States
| |
Collapse
|
3
|
Song J, Wan M, Yang Y, Gao L, Fang W. Development of accurate coarse-grained force fields for weakly polar groups by an indirect parameterization strategy. Phys Chem Chem Phys 2021; 23:6763-6774. [PMID: 33720253 DOI: 10.1039/d1cp00032b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Coarse-grained (CG) molecular dynamics simulations are widely used to predict morphological structures and interpret mechanisms of mesoscopic behavior between the scope of traditional experiments and all-atom simulations. However, most current CG force fields (FFs) are not precise enough, especially for polar molecules or functional groups. A main obstacle in developing accurate CG FFs for polar molecules is the freezing problem met at room temperature. In this work, we introduce an indirect parametrization strategy for weakly polar groups by considering their short-chain homologs to avoid freezing. Here, a polar group containing three to four heavy atoms is mapped into one CG bead that is connected to one alkyl bead composed of three or four carbons. The CG beads interact via 4-parameter nonbonded Morse potentials and harmonic bonded potentials. An efficient meta-multilinear interpolation parameterization algorithm, as recently developed by us, is used to rigorously optimize the force parameters. Satisfactory accuracy is witnessed in terms of the density, heat of vaporization, surface tension, and solvation free energy of the homologs of twelve polar molecules, all deviating from the experiment by less than 5%. The transferability of the current FF is indicated by the predicted density, heat of vaporization, and end-to-end distance distributions of fatty acid methyl esters composed of multiple functional groups parameterized in this work.
Collapse
Affiliation(s)
- Junjie Song
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, 19 Xin-Jie-Kou-Wai Street, Beijing 100875, China.
| | | | | | | | | |
Collapse
|
4
|
Lynch C, Rao S, Sansom MSP. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem Rev 2020; 120:10298-10335. [PMID: 32841020 PMCID: PMC7517714 DOI: 10.1021/acs.chemrev.9b00830] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Indexed: 12/18/2022]
Abstract
This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
Collapse
Affiliation(s)
- Charlotte
I. Lynch
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Shanlin Rao
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| |
Collapse
|
5
|
Kacar G, de With G. Parametrizing hydrogen bond interactions in dissipative particle dynamics simulations: The case of water, methanol and their binary mixtures. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.112581] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
|
6
|
de Souza RM, Ratochinski RH, Karttunen M, Dias LG. Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles. J Chem Inf Model 2020; 60:522-536. [PMID: 31714768 DOI: 10.1021/acs.jcim.9b00790] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.
Collapse
Affiliation(s)
- R M de Souza
- Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.,Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil.,The Center for Advanced Materials and Biomaterials Research , The University of Western Ontario , London , Ontario , Canada N6K 3K7
| | - R H Ratochinski
- Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil
| | - Mikko Karttunen
- Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.,The Center for Advanced Materials and Biomaterials Research , The University of Western Ontario , London , Ontario , Canada N6K 3K7.,Department of Applied Mathematics , The University of Western Ontario , London , Ontario , Canada N6A 5B7
| | - L G Dias
- Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil
| |
Collapse
|
7
|
Miyazaki Y, Okazaki S, Shinoda W. pSPICA: A Coarse-Grained Force Field for Lipid Membranes Based on a Polar Water Model. J Chem Theory Comput 2019; 16:782-793. [DOI: 10.1021/acs.jctc.9b00946] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Yusuke Miyazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Susumu Okazaki
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Wataru Shinoda
- Department of Materials Chemistry, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| |
Collapse
|
8
|
Machado MR, Zeida A, Darré L, Pantano S. From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field. Interface Focus 2019; 9:20180085. [PMID: 31065347 PMCID: PMC6501346 DOI: 10.1098/rsfs.2018.0085] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/11/2019] [Indexed: 12/11/2022] Open
Abstract
Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.
Collapse
Affiliation(s)
- Matías R. Machado
- Institut Pasteur de Montevideo, Group of Biomolecular Simulations, Mataojo 2020, CP 11400 Montevideo, Uruguay
| | - Ari Zeida
- Departamento de Bioquímica and Center for Free Radical and Biomedical Research, Facultad de Medicina, Universidad de la República, Montevideo, Uruguay
| | - Leonardo Darré
- Institut Pasteur de Montevideo, Group of Biomolecular Simulations, Mataojo 2020, CP 11400 Montevideo, Uruguay
- Institut Pasteur de Montevideo, Functional Genomics Unit, Mataojo 2020, CP 11400 Montevideo, Uruguay
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Group of Biomolecular Simulations, Mataojo 2020, CP 11400 Montevideo, Uruguay
| |
Collapse
|
9
|
Nieto-Draghi C, Rousseau B. Thermodynamically Consistent Force Field for Coarse-Grained Modeling of Aqueous Electrolyte Solution. J Phys Chem B 2019; 123:2424-2431. [DOI: 10.1021/acs.jpcb.8b11190] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Carlos Nieto-Draghi
- IFP Energies nouvelles, 1-4 Avenue de Bois Préau, 92852 Rueil-Malmaison, France
| | - Bernard Rousseau
- Laboratoire de Chimie Physique, UMR 8000 CNRS, Université Paris-Sud, 91405 Orsay, France
| |
Collapse
|
10
|
Ghoufi A, Malfreyt P. Calculation of the surface tension of water: 40 years of molecular simulations. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1513648] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Aziz Ghoufi
- Institut de Physique de Rennes, Université Rennes 1, Rennes, France
| | - Patrice Malfreyt
- Institut de Chimie de Clermont-Ferrand (ICCF), Université Clermont Auvergne, CNRS, SIGMA Clermont, Clermont-Ferrand, France
| |
Collapse
|
11
|
Bejagam KK, Singh S, An Y, Berry C, Deshmukh SA. PSO-Assisted Development of New Transferable Coarse-Grained Water Models. J Phys Chem B 2018; 122:1958-1971. [DOI: 10.1021/acs.jpcb.7b10542] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Karteek K. Bejagam
- Department
of Chemical Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Samrendra Singh
- CNH Industrial, Burr Ridge, Chicago, Illinois 60527, United States
| | - Yaxin An
- Department
of Chemical Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Carter Berry
- Department
of Chemical Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Sanket A. Deshmukh
- Department
of Chemical Engineering, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
12
|
Kacar G. Characterizing the structure and properties of dry and wet polyethylene glycol using multi-scale simulations. Phys Chem Chem Phys 2018; 20:12303-12311. [DOI: 10.1039/c8cp01802b] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Multi-scale simulations to study the structure and material properties of PEG in dry and wet conditions.
Collapse
Affiliation(s)
- Gokhan Kacar
- Department of Genetics and Bioengineering
- Faculty of Engineering
- Trakya University
- Edirne
- Turkey
| |
Collapse
|
13
|
Onufriev AV, Izadi S. Water models for biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1347] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Alexey V. Onufriev
- Department of Physics; Virginia Tech; Blacksburg VA USA
- Department of Computer Science; Virginia Tech; Blacksburg VA USA
- Center for Soft Matter and Biological Physics; Virginia Tech; Blacksburg VA USA
| | - Saeed Izadi
- Early Stage Pharmaceutical Development; Genentech Inc.; South San Francisco, CA USA
| |
Collapse
|
14
|
Rodríguez-López T, Khalak Y, Karttunen M. Non-conformal coarse-grained potentials for water. J Chem Phys 2017; 147:134108. [DOI: 10.1063/1.4985914] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Tonalli Rodríguez-López
- Departamento de Física, Universidad Autónoma Metropolitana, Iztapalapa, Apdo 55 534, Mexico DF 09340, Mexico
- Department of Chemistry and Waterloo Institute of Nanotechnology, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L3G1, Canada
| | - Yuriy Khalak
- Department of Mathematics and Computer Science and Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, MetaForum, 5600 MB Eindhoven, The Netherlands
| | - Mikko Karttunen
- Department of Mathematics and Computer Science and Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, MetaForum, 5600 MB Eindhoven, The Netherlands
- Department of Chemistry, Western University, 1151 Richmond Street, London, Ontario N6A5B7, Canada
- Department of Applied Mathematics, Western University, 1151 Richmond Street, London, Ontario N6A5B7, Canada
| |
Collapse
|
15
|
Machado MR, González HC, Pantano S. MD Simulations of Viruslike Particles with Supra CG Solvation Affordable to Desktop Computers. J Chem Theory Comput 2017; 13:5106-5116. [DOI: 10.1021/acs.jctc.7b00659] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matı́as R. Machado
- Biomolecular Simulations
Group, Institut Pasteur de Montevideo, Mataojo 2020, Montevideo CP 11400, Uruguay
| | - Humberto C. González
- Biomolecular Simulations
Group, Institut Pasteur de Montevideo, Mataojo 2020, Montevideo CP 11400, Uruguay
| | - Sergio Pantano
- Biomolecular Simulations
Group, Institut Pasteur de Montevideo, Mataojo 2020, Montevideo CP 11400, Uruguay
| |
Collapse
|
16
|
Properties of the polarizable MARTINI water model: A comparative study for aqueous electrolyte solutions. J Mol Liq 2015. [DOI: 10.1016/j.molliq.2015.08.062] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
17
|
Pluhackova K, Böckmann RA. Biomembranes in atomistic and coarse-grained simulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015. [PMID: 26194872 DOI: 10.1088/0953-8984/27/32/323103] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.
Collapse
Affiliation(s)
- Kristyna Pluhackova
- Computational Biology, Department of Biology, Friedrich-Alexander Universität Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| | | |
Collapse
|
18
|
Darré L, Machado MR, Brandner AF, González HC, Ferreira S, Pantano S. SIRAH: A Structurally Unbiased Coarse-Grained Force Field for Proteins with Aqueous Solvation and Long-Range Electrostatics. J Chem Theory Comput 2015; 11:723-39. [DOI: 10.1021/ct5007746] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Leonardo Darré
- Institut Pasteur de Montevideo, Montevideo, Uruguay
- Department
of Chemistry, King’s College, London, United Kingdom
| | | | | | | | | | | |
Collapse
|
19
|
Peter EK, Pivkin IV. A polarizable coarse-grained water model for dissipative particle dynamics. J Chem Phys 2014; 141:164506. [DOI: 10.1063/1.4899317] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Emanuel K. Peter
- Institute of Computational Science, Faculty of Informatics, University of Lugano, Lugano, Switzerland
| | - Igor V. Pivkin
- Institute of Computational Science, Faculty of Informatics, University of Lugano, Lugano, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| |
Collapse
|
20
|
Shen L, Hu H. Resolution-Adapted All-Atomic and Coarse-Grained Model for Biomolecular Simulations. J Chem Theory Comput 2014; 10:2528-36. [DOI: 10.1021/ct401029k] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Lin Shen
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
| | - Hao Hu
- Department
of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China
- The University of Hong Kong-Shenzhen Institute of Research and Innovation, Kejizhong Second Road, Shenzhen, China
| |
Collapse
|
21
|
Ingólfsson HI, Lopez CA, Uusitalo JJ, de Jong DH, Gopal SM, Periole X, Marrink SJ. The power of coarse graining in biomolecular simulations. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2014; 4:225-248. [PMID: 25309628 PMCID: PMC4171755 DOI: 10.1002/wcms.1169] [Citation(s) in RCA: 325] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Computational modeling of biological systems is challenging because of the multitude of spatial and temporal scales involved. Replacing atomistic detail with lower resolution, coarse grained (CG), beads has opened the way to simulate large-scale biomolecular processes on time scales inaccessible to all-atom models. We provide an overview of some of the more popular CG models used in biomolecular applications to date, focusing on models that retain chemical specificity. A few state-of-the-art examples of protein folding, membrane protein gating and self-assembly, DNA hybridization, and modeling of carbohydrate fibers are used to illustrate the power and diversity of current CG modeling.
Collapse
Affiliation(s)
- Helgi I Ingólfsson
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Cesar A Lopez
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Jaakko J Uusitalo
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Djurre H de Jong
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Srinivasa M Gopal
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Xavier Periole
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of GroningenGroningen, The Netherlands
| |
Collapse
|
22
|
|
23
|
Abstract
By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
24
|
Fogarty JC, Chiu SW, Kirby P, Jakobsson E, Pandit SA. Automated optimization of water-water interaction parameters for a coarse-grained model. J Phys Chem B 2014; 118:1603-11. [PMID: 24460506 PMCID: PMC3983376 DOI: 10.1021/jp409545x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have developed an automated parameter optimization software framework (ParOpt) that implements the Nelder-Mead simplex algorithm and applied it to a coarse-grained polarizable water model. The model employs a tabulated, modified Morse potential with decoupled short- and long-range interactions incorporating four water molecules per interaction site. Polarizability is introduced by the addition of a harmonic angle term defined among three charged points within each bead. The target function for parameter optimization was based on the experimental density, surface tension, electric field permittivity, and diffusion coefficient. The model was validated by comparison of statistical quantities with experimental observation. We found very good performance of the optimization procedure and good agreement of the model with experiment.
Collapse
Affiliation(s)
- Joseph C Fogarty
- Department of Physics, University of South Florida , Tampa, Florida 33620, United States
| | | | | | | | | |
Collapse
|
25
|
Gonzalez HC, Darré L, Pantano S. Transferable Mixing of Atomistic and Coarse-Grained Water Models. J Phys Chem B 2013; 117:14438-48. [DOI: 10.1021/jp4079579] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Leonardo Darré
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo 11400, Uruguay
- Department
of Chemistry, King’s College London, London, United Kingdom
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo 11400, Uruguay
| |
Collapse
|
26
|
Geometric isotope effects on small chloride ion water clusters with path integral molecular dynamics simulations. Chem Phys 2013. [DOI: 10.1016/j.chemphys.2013.10.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
27
|
Howard JC, Tschumper GS. Wavefunction methods for the accurate characterization of water clusters. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1168] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry University of Mississippi, University Mississippi USA
| |
Collapse
|
28
|
Darré L, Tek A, Baaden M, Pantano S. Mixing Atomistic and Coarse Grain Solvation Models for MD Simulations: Let WT4 Handle the Bulk. J Chem Theory Comput 2012; 8:3880-94. [DOI: 10.1021/ct3001816] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Leonardo Darré
- Institut Pasteur de Montevideo,
Mataojo 2020, CP 11400, Uruguay
| | - Alex Tek
- Laboratoire de Biochimie Théorique,
CNRS, UPR9080, Univ Paris Diderot, Sorbonne Paris Cité. 13
rue Pierre et Marie Curie, 75005, Paris, France
- Université Pierre et Marie
Curie, UPMC Sorbonne Universités, 4 place Jussieu 75005 Paris,
France
| | - Marc Baaden
- Laboratoire de Biochimie Théorique,
CNRS, UPR9080, Univ Paris Diderot, Sorbonne Paris Cité. 13
rue Pierre et Marie Curie, 75005, Paris, France
| | - Sergio Pantano
- Institut Pasteur de Montevideo,
Mataojo 2020, CP 11400, Uruguay
| |
Collapse
|