1
|
Kubincová A, Riniker S, Hünenberger PH. Solvent-scaling as an alternative to coarse-graining in adaptive-resolution simulations: The adaptive solvent-scaling (AdSoS) scheme. J Chem Phys 2021; 155:094107. [PMID: 34496576 DOI: 10.1063/5.0057384] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new approach termed Adaptive Solvent-Scaling (AdSoS) is introduced for performing simulations of a solute embedded in a fine-grained (FG) solvent region itself surrounded by a coarse-grained (CG) solvent region, with a continuous FG ↔ CG switching of the solvent resolution across a buffer layer. Instead of relying on a distinct CG solvent model, the AdSoS scheme is based on CG models defined by a dimensional scaling of the FG solvent by a factor s, accompanied by an s-dependent modulation of the atomic masses and interaction parameters. The latter changes are designed to achieve an isomorphism between the dynamics of the FG and CG models, and to preserve the dispersive and dielectric solvation properties of the solvent with respect to a solute at FG resolution. This scaling approach offers a number of advantages compared to traditional coarse-graining: (i) the CG parameters are immediately related to those of the FG model (no need to parameterize a distinct CG model); (ii) nearly ideal mixing is expected for CG variants with similar s-values (ideal mixing holding in the limit of identical s-values); (iii) the solvent relaxation timescales should be preserved (no dynamical acceleration typical for coarse-graining); (iv) the graining level NG (number of FG molecules represented by one CG molecule) can be chosen arbitrarily (in particular, NG = s3 is not necessarily an integer); and (v) in an adaptive-resolution scheme, this level can be varied continuously as a function of the position (without requiring a bundling mechanism), and this variation occurs at a constant number of particles per molecule (no occurrence of fractional degrees of freedom in the buffer layer). By construction, the AdSoS scheme minimizes the thermodynamic mismatch between the different regions of the adaptive-resolution system, leading to a nearly homogeneous scaled solvent density s3ρ. Residual density artifacts in and at the surface of the boundary layer can easily be corrected by means of a grid-based biasing potential constructed in a preliminary pure-solvent simulation. This article introduces the AdSoS scheme and provides an initial application to pure atomic liquids (no solute) with Lennard-Jones plus Coulomb interactions in a slab geometry.
Collapse
Affiliation(s)
- Alžbeta Kubincová
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zurich, Vladimir Prelog-Weg 2, CH-8093 Zürich, Switzerland
| |
Collapse
|
2
|
Thaler S, Praprotnik M, Zavadlav J. Back-mapping augmented adaptive resolution simulation. J Chem Phys 2020; 153:164118. [PMID: 33138420 DOI: 10.1063/5.0025728] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Concurrent multiscale techniques such as Adaptive Resolution Scheme (AdResS) can offer ample computational advantages over conventional atomistic (AT) molecular dynamics simulations. However, they typically rely on aphysical hybrid regions to maintain numerical stability when high-resolution degrees of freedom (DOFs) are randomly re-inserted at the resolution interface. We propose an Energy Minimized AT (DOF) Insertion (EMATI) method that uses an informed rather than random AT DOF insertion to tackle the root cause of the issue, i.e., overlapping AT potentials. EMATI enables us to directly couple AT and coarse-grained resolutions without any modifications of the interaction potentials. We exemplify AdResS-EMATI in a system of liquid butane and show that it yields improved structural and thermodynamic properties at the interface compared to competing AdResS approaches. Furthermore, our approach extends the applicability of the AdResS without a hybrid region to systems for which force capping is inadequate.
Collapse
Affiliation(s)
- S Thaler
- Professorship of Multiscale Modeling of Fluid Materials, Department of Mechanical Engineering, Technical University of Munich, Munich, Germany
| | - M Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, SI-1001 Ljubljana, Slovenia
| | - J Zavadlav
- Professorship of Multiscale Modeling of Fluid Materials, Department of Mechanical Engineering, Technical University of Munich, Munich, Germany
| |
Collapse
|
3
|
Taylor PA, Jayaraman A. Molecular Modeling and Simulations of Peptide–Polymer Conjugates. Annu Rev Chem Biomol Eng 2020; 11:257-276. [DOI: 10.1146/annurev-chembioeng-092319-083243] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Peptide–polymer conjugates are a class of soft materials composed of covalently linked blocks of protein/polypeptides and synthetic/natural polymers. These materials are practically useful in biological applications, such as drug delivery, DNA/gene delivery, and antimicrobial coatings, as well as nonbiological applications, such as electronics, separations, optics, and sensing. Given their broad applicability, there is motivation to understand the molecular and macroscale structure, dynamics, and thermodynamic behavior exhibited by such materials. We focus on the past and ongoing molecular simulation studies aimed at obtaining such fundamental understanding and predicting molecular design rules for the target function. We describe briefly the experimental work in this field that validates or motivates these computational studies. We also describe the various models (e.g., atomistic, coarse-grained, or hybrid) and simulation methods (e.g., stochastic versus deterministic, enhanced sampling) that have been used and the types of questions that have been answered using these computational approaches.
Collapse
Affiliation(s)
- Phillip A. Taylor
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, USA
| | - Arthi Jayaraman
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, USA
- Department of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
4
|
Liu Y, De Vries AH, Barnoud J, Pezeshkian W, Melcr J, Marrink SJ. Dual Resolution Membrane Simulations Using Virtual Sites. J Phys Chem B 2020; 124:3944-3953. [PMID: 32314586 PMCID: PMC7232679 DOI: 10.1021/acs.jpcb.0c01842] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
![]()
All-atomistic
(AA) and coarse-grain (CG) simulations have been
successfully applied to investigate a broad range of biomolecular
processes. However, the accessible time and length scales of AA simulation
are limited and the specific molecular details of CG simulation are
simplified. Here, we propose a virtual site (VS) based hybrid scheme
that can concurrently couple AA and CG resolutions in a single membrane
simulation, mitigating the shortcomings of either representation.
With some adjustments to make the AA and CG force fields compatible,
we demonstrate that lipid bilayer properties are well kept in our
hybrid approach. Our VS hybrid method was also applied to simulate
a small lipid vesicle, with the inner leaflet and interior solvent
represented in AA, and the outer leaflet together with exterior solvent
at the CG level. Our multiscale method opens the way to investigate
biomembrane properties at increased computational efficiency, in particular
applications involving large solvent filled regions.
Collapse
Affiliation(s)
- Yang Liu
- Groningen Biomolecular Sciences and Biotechnology Institute and the Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Alex H De Vries
- Groningen Biomolecular Sciences and Biotechnology Institute and the Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Jonathan Barnoud
- Groningen Biomolecular Sciences and Biotechnology Institute and the Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Weria Pezeshkian
- Groningen Biomolecular Sciences and Biotechnology Institute and the Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Josef Melcr
- Groningen Biomolecular Sciences and Biotechnology Institute and the Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and the Zernike Institute for Advanced Material, University of Groningen, 9747 AG Groningen, The Netherlands
| |
Collapse
|
5
|
Renevey A, Riniker S. Benchmarking Hybrid Atomistic/Coarse-Grained Schemes for Proteins with an Atomistic Water Layer. J Phys Chem B 2019; 123:3033-3042. [DOI: 10.1021/acs.jpcb.8b12149] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Annick Renevey
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
6
|
Friedman R, Khalid S, Aponte-Santamaría C, Arutyunova E, Becker M, Boyd KJ, Christensen M, Coimbra JTS, Concilio S, Daday C, van Eerden FJ, Fernandes PA, Gräter F, Hakobyan D, Heuer A, Karathanou K, Keller F, Lemieux MJ, Marrink SJ, May ER, Mazumdar A, Naftalin R, Pickholz M, Piotto S, Pohl P, Quinn P, Ramos MJ, Schiøtt B, Sengupta D, Sessa L, Vanni S, Zeppelin T, Zoni V, Bondar AN, Domene C. Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology. J Membr Biol 2018; 251:609-631. [PMID: 30350011 PMCID: PMC6244758 DOI: 10.1007/s00232-018-0050-y] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2018] [Accepted: 10/03/2018] [Indexed: 12/23/2022]
Affiliation(s)
- Ran Friedman
- Department of Chemistry and Biomedical Sciences and Centre of Excellence "Biomaterials Chemistry", Linnæus University, Kalmar, Sweden.
| | - Syma Khalid
- University of Southampton, Southampton, SO17 1BJ, UK
| | - Camilo Aponte-Santamaría
- Max Planck Tandem Group in Computational Biophysics, University of Los Andes, Bogotá, Colombia.,Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Heidelberg, Germany
| | - Elena Arutyunova
- Department of Biochemistry, University of Alberta, Edmonton, Canada
| | | | - Kevin J Boyd
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, USA
| | - Mikkel Christensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark.,Interdisciplinary Nanoscience center (iNANO), Aarhus University, Aarhus, Denmark.,Sino-Danish Center for Education and Research, Beijing, China
| | - João T S Coimbra
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Porto, Portugal
| | - Simona Concilio
- Department of Industrial Engineering, University of Salerno, Fisciano, SA, Italy
| | - Csaba Daday
- Heidelberg Institute for Theoretical Studies, Heidelberg, Germany
| | | | - Pedro A Fernandes
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Porto, Portugal
| | - Frauke Gräter
- Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Heidelberg, Germany.,Heidelberg Institute for Theoretical Studies, Heidelberg, Germany
| | | | | | - Konstantina Karathanou
- Department of Physics, Theoretical Molecular Biophysics Group, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | | | - M Joanne Lemieux
- Department of Biochemistry, University of Alberta, Edmonton, Canada
| | | | - Eric R May
- Department of Molecular and Cell Biology, University of Connecticut, Storrs, CT, USA
| | - Antara Mazumdar
- GBB Institute, University of Groningen, Groningen, The Netherlands
| | - Richard Naftalin
- Physiology and Vascular Biology Departments, King's College London School of Medicine, London, UK
| | - Mónica Pickholz
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, CONICET-Universidad de Buenos Aires, IFIBA, Buenos Aires, Argentina
| | - Stefano Piotto
- Department of Pharmacy, University of Salerno, Fisciano, SA, Italy
| | - Peter Pohl
- Institute of Biophysics, Johannes Kepler University, Linz, Austria
| | - Peter Quinn
- Biochemistry Department, King's College London, London, UK
| | - Maria J Ramos
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, Porto, Portugal
| | - Birgit Schiøtt
- Department of Chemistry, Aarhus University, Aarhus, Denmark.,Interdisciplinary Nanoscience center (iNANO), Aarhus University, Aarhus, Denmark
| | - Durba Sengupta
- Physical Chemistry Division, National Chemical Laboratory, Pune, India
| | - Lucia Sessa
- Department of Pharmacy, University of Salerno, Fisciano, SA, Italy
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Fribourg, Switzerland
| | - Talia Zeppelin
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Valeria Zoni
- Department of Biology, University of Fribourg, Fribourg, Switzerland
| | - Ana-Nicoleta Bondar
- Department of Physics, Theoretical Molecular Biophysics Group, Freie Universität Berlin, Arnimallee 14, 14195, Berlin, Germany
| | - Carmen Domene
- Department of Chemistry, University of Bath, Claverton Down Bath, BA2 7AY, UK.,Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford, OX1 3TA, UK
| |
Collapse
|
7
|
Fiorentini R, Kremer K, Potestio R, Fogarty AC. Using force-based adaptive resolution simulations to calculate solvation free energies of amino acid sidechain analogues. J Chem Phys 2018; 146:244113. [PMID: 28668024 DOI: 10.1063/1.4989486] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
The calculation of free energy differences is a crucial step in the characterization and understanding of the physical properties of biological molecules. In the development of efficient methods to compute these quantities, a promising strategy is that of employing a dual-resolution representation of the solvent, specifically using an accurate model in the proximity of a molecule of interest and a simplified description elsewhere. One such concurrent multi-resolution simulation method is the Adaptive Resolution Scheme (AdResS), in which particles smoothly change their resolution on-the-fly as they move between different subregions. Before using this approach in the context of free energy calculations, however, it is necessary to make sure that the dual-resolution treatment of the solvent does not cause undesired effects on the computed quantities. Here, we show how AdResS can be used to calculate solvation free energies of small polar solutes using Thermodynamic Integration (TI). We discuss how the potential-energy-based TI approach combines with the force-based AdResS methodology, in which no global Hamiltonian is defined. The AdResS free energy values agree with those calculated from fully atomistic simulations to within a fraction of kBT. This is true even for small atomistic regions whose size is on the order of the correlation length, or when the properties of the coarse-grained region are extremely different from those of the atomistic region. These accurate free energy calculations are possible because AdResS allows the sampling of solvation shell configurations which are equivalent to those of fully atomistic simulations. The results of the present work thus demonstrate the viability of the use of adaptive resolution simulation methods to perform free energy calculations and pave the way for large-scale applications where a substantial computational gain can be attained.
Collapse
Affiliation(s)
- Raffaele Fiorentini
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Aoife C Fogarty
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
8
|
Renevey A, Riniker S. Improved accuracy of hybrid atomistic/coarse-grained simulations using reparametrised interactions. J Chem Phys 2018; 146:124131. [PMID: 28388132 DOI: 10.1063/1.4979128] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Reducing the number of degrees of freedom in molecular models-so-called coarse-graining-is a popular approach to increase the accessible time scales and system sizes in molecular dynamics simulations. It involves, however, per se a loss of information. In order to retain a high accuracy in the region of interest, hybrid methods that combine two levels of resolution in a single system are an attractive trade-off. Hybrid atomistic (AT)/coarse-grained (CG) simulations have previously been shown to preserve the secondary structure elements of AT proteins in CG water but to cause an artificial increase in intramolecular hydrogen bonds, resulting in a reduced flexibility of the proteins. Recently, it was found that the AT-CG interactions employed in these simulations were too favourable for apolar solutes and not favourable enough for polar solutes. Here, the AT-CG interactions are reparametrised to reproduce the solvation free energy of a series of AT alkanes and side-chain analogues in CG water, while retaining the good mixing behaviour of AT water with CG water. The new AT-CG parameters are tested in hybrid simulations of four proteins in CG water. Structural and dynamic properties are compared to those obtained in fully AT simulations and, if applicable, to experimental data. The results show that the artificial increase of intramolecular hydrogen bonds is drastically reduced, leading to a better reproduction of the structural properties and flexibility of the proteins in atomistic water, without the need for an atomistic solvent layer.
Collapse
Affiliation(s)
- Annick Renevey
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
9
|
Zavadlav J, Marrink SJ, Praprotnik M. Multiscale Simulation of Protein Hydration Using the SWINGER Dynamical Clustering Algorithm. J Chem Theory Comput 2018; 14:1754-1761. [PMID: 29439560 DOI: 10.1021/acs.jctc.7b01129] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
To perform computationally efficient concurrent multiscale simulations of biological macromolecules in solution, where the all-atom (AT) models are coupled to supramolecular coarse-grained (SCG) solvent models, previous studies resorted to modified AT water models, such as the bundled-simple point charge (SPC) models, that use semiharmonic springs to restrict the relative movement of water molecules within a cluster. Those models can have a significant impact on the simulated biomolecules and can lead, for example, to a partial unfolding of a protein. In this work, we employ the recently developed alternative approach with a dynamical clustering algorithm, SWINGER, which enables a direct coupling of original unmodified AT and SCG water models. We perform an adaptive resolution molecular dynamics simulation of a Trp-Cage miniprotein in multiscale water, where the standard SPC water model is interfaced with the widely used MARTINI SCG model, and demonstrate that, compared to the corresponding full-blown AT simulations, the structural and dynamic properties of the solvated protein and surrounding solvent are well reproduced by our approach.
Collapse
Affiliation(s)
- Julija Zavadlav
- Computational Science & Engineering Laboratory , ETH Zurich , Clausiusstrasse 33 , CH-8092 Zurich , Switzerland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials , University of Groningen , Nijenborgh 7 , 9747 AG Groningen , The Netherlands
| | - Matej Praprotnik
- Laboratory for Molecular Modeling , National Institute of Chemistry , Hajdrihova 19 , SI-1001 Ljubljana , Slovenia.,Department of Physics, Faculty of Mathematics and Physics , University of Ljubljana , Jadranska 19 , SI-1000 Ljubljana , Slovenia
| |
Collapse
|
10
|
Sinitskiy AV, Voth GA. Quantum mechanics/coarse-grained molecular mechanics (QM/CG-MM). J Chem Phys 2018; 148:014102. [DOI: 10.1063/1.5006810] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anton V. Sinitskiy
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
11
|
Zavadlav J, Bevc S, Praprotnik M. Adaptive resolution simulations of biomolecular systems. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2017; 46:821-835. [PMID: 28905203 DOI: 10.1007/s00249-017-1248-0] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2017] [Revised: 06/12/2017] [Accepted: 08/15/2017] [Indexed: 10/18/2022]
Abstract
In this review article, we discuss and analyze some recently developed hybrid atomistic-mesoscopic solvent models for multiscale biomolecular simulations. We focus on the biomolecular applications of the adaptive resolution scheme (AdResS), which allows solvent molecules to change their resolution back and forth between atomistic and coarse-grained representations according to their positions in the system. First, we discuss coupling of atomistic and coarse-grained models of salt solution using a 1-to-1 molecular mapping-i.e., one coarse-grained bead represents one water molecule-for development of a multiscale salt solution model. In order to make use of coarse-grained molecular models that are compatible with the MARTINI force field, one has to resort to a supramolecular mapping, in particular to a 4-to-1 mapping, where four water molecules are represented with one coarse-grained bead. To this end, bundled atomistic water models are employed, i.e., the relative movement of water molecules that are mapped to the same coarse-grained bead is restricted by employing harmonic springs. Supramolecular coupling has recently also been extended to polarizable coarse-grained water models with explicit charges. Since these coarse-grained models consist of several interaction sites, orientational degrees of freedom of the atomistic and coarse-grained representations are coupled via a harmonic energy penalty term. The latter aligns the dipole moments of both representations. The reviewed multiscale solvent models are ready to be used in biomolecular simulations, as illustrated in a few examples.
Collapse
Affiliation(s)
- Julija Zavadlav
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia.,Chair of Computational Science, ETH Zurich, Clausiusstrasse 33, 8092, Zurich, Switzerland
| | - Staš Bevc
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia
| | - Matej Praprotnik
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia. .,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia.
| |
Collapse
|
12
|
Delgado-Buscalioni R. Thermodynamics of adaptive molecular resolution. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2016; 374:rsta.2016.0152. [PMID: 27698044 DOI: 10.1098/rsta.2016.0152] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 07/05/2016] [Indexed: 05/27/2023]
Abstract
A relatively general thermodynamic formalism for adaptive molecular resolution (AMR) is presented. The description is based on the approximation of local thermodynamic equilibrium and considers the alchemic parameter λ as the conjugate variable of the potential energy difference between the atomistic and coarse-grained model Φ=U(1)-U(0) The thermodynamic formalism recovers the relations obtained from statistical mechanics of H-AdResS (Español et al, J. Chem. Phys. 142, 064115, 2015 (doi:10.1063/1.4907006)) and provides relations between the free energy compensation and thermodynamic potentials. Inspired by this thermodynamic analogy, several generalizations of AMR are proposed, such as the exploration of new Maxwell relations and how to treat λ and Φ as 'real' thermodynamic variablesThis article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.
Collapse
Affiliation(s)
- R Delgado-Buscalioni
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, Campus de Cantoblanco, 28049 Madrid, Spain Institute for Condensed Matter Physics, IFIMAC, Campus de Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
13
|
Li ZL, Ding HM, Ma YQ. Interaction of peptides with cell membranes: insights from molecular modeling. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:083001. [PMID: 26828575 DOI: 10.1088/0953-8984/28/8/083001] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The investigation of the interaction of peptides with cell membranes is the focus of active research. It can enhance the understanding of basic membrane functions such as membrane transport, fusion, and signaling processes, and it may shed light on potential applications of peptides in biomedicine. In this review, we will present current advances in computational studies on the interaction of different types of peptides with the cell membrane. Depending on the properties of the peptide, membrane, and external environment, the peptide-membrane interaction shows a variety of different forms. Here, on the basis of recent computational progress, we will discuss how different peptides could initiate membrane pores, translocate across the membrane, induce membrane endocytosis, produce membrane curvature, form fibrils on the membrane surface, as well as interact with functional membrane proteins. Finally, we will present a conclusion summarizing recent progress and providing some specific insights into future developments in this field.
Collapse
Affiliation(s)
- Zhen-lu Li
- National Laboratory of Solid State Microstructures and Department of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, People's Republic of China
| | | | | |
Collapse
|
14
|
Yan XC, Tirado-Rives J, Jorgensen WL. Hydration Properties and Solvent Effects for All-Atom Solutes in Polarizable Coarse-Grained Water. J Phys Chem B 2016; 120:8102-14. [PMID: 26901452 DOI: 10.1021/acs.jpcb.6b00399] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Due to the importance of water in chemical and biological systems, a coarse-grained representation of the solvent can greatly simplify the description of the system while retaining key thermodynamic properties of the medium. A multiscale solvation model that couples all-atom solutes and polarizable Martini coarse-grained water (AAX/CGS) is developed to reproduce free energies of hydration of organic solutes. Using Monte Carlo/free energy perturbation (MC/FEP) calculations, results from multiscale and all-atom simulations are compared. Improved accuracy is obtained with the AAX/CGS approach for hydrophobic and sulfur- or halogen-containing solutes, but larger deviations are found for polar solute molecules where hydrogen bonding is featured. Furthermore, solvent effects on conformational and tautomeric equilibria of AA solutes were investigated using AA, CG, and GB/SA solvent models. It is found that the CG solvent model can reproduce well the medium effects from experiment and AA simulations; however, the GB/SA solvent model fails in some cases. A 7-30-fold reduction in computational cost is found for the present AAX/CGS multiscale simulations compared to the AA alternative.
Collapse
Affiliation(s)
- Xin Cindy Yan
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - William L Jorgensen
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| |
Collapse
|