1
|
Coste A, Slejko E, Zavadlav J, Praprotnik M. Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media. J Chem Theory Comput 2024; 20:411-420. [PMID: 38118122 PMCID: PMC10782447 DOI: 10.1021/acs.jctc.3c00984] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 12/04/2023] [Accepted: 12/04/2023] [Indexed: 12/22/2023]
Abstract
Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.
Collapse
Affiliation(s)
- Amaury Coste
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Ljubljana SI-1001, Slovenia
| | - Ema Slejko
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Ljubljana SI-1001, Slovenia
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana SI-1000, Slovenia
| | - Julija Zavadlav
- Professorship
of Multiscale Modeling of Fluid Materials, TUM School of Engineering
and Design, Technical University of Munich, Garching Near Munich DE-85748, Germany
| | - Matej Praprotnik
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Ljubljana SI-1001, Slovenia
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana SI-1000, Slovenia
| |
Collapse
|
2
|
Smith ER, Theodorakis PE. Multiscale simulation of fluids: coupling molecular and continuum. Phys Chem Chem Phys 2024; 26:724-744. [PMID: 38113114 DOI: 10.1039/d3cp03579d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2023]
Abstract
Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.
Collapse
Affiliation(s)
- Edward R Smith
- Department of Mechanical and Aerospace Engineering, Brunel University London, Uxbridge, Middlesex UB8 3PH, UK.
| | | |
Collapse
|
3
|
Erban R, Togashi Y. Asymmetric Periodic Boundary Conditions for All-Atom Molecular Dynamics and Coarse-Grained Simulations of Nucleic Acids. J Phys Chem B 2023; 127:8257-8267. [PMID: 37713594 PMCID: PMC10544013 DOI: 10.1021/acs.jpcb.3c03887] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2023] [Revised: 07/31/2023] [Indexed: 09/17/2023]
Abstract
Periodic boundary conditions are commonly applied in molecular dynamics simulations in the microcanonical (NVE), canonical (NVT), and isothermal-isobaric (NpT) ensembles. In their simplest application, a biological system of interest is placed in the middle of a solvation box, which is chosen 'sufficiently large' to minimize any numerical artifacts associated with the periodic boundary conditions. This practical approach brings limitations to the size of biological systems that can be simulated. Here, we study simulations of effectively infinitely long nucleic acids, which are solvated in the directions perpendicular to the polymer chain, while periodic boundary conditions are also applied along the polymer chain. We study the effects of these asymmetric periodic boundary conditions (APBC) on the simulated results, including the mechanical properties of biopolymers and the properties of the surrounding solvent. To get some further insights into the advantages of using the APBC, a coarse-grained worm-like chain model is first studied, illustrating how the persistence length can be extracted from the local properties of the polymer chain, which are less affected by the APBC than some global averages. This is followed by all-atom molecular dynamics simulations of DNA in ionic solutions, where we use the APBC to investigate sequence-dependent properties of DNA molecules and properties of the surrounding solvent.
Collapse
Affiliation(s)
- Radek Erban
- Mathematical
Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock
Road, Oxford OX2 6GG, U.K.
| | - Yuichi Togashi
- Department
of Bioinformatics, College of Life Sciences, Ritsumeikan University, 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
- RIKEN
Center for Biosystems Dynamics Research, 2-2-3 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
4
|
Fiorentini R, Tarenzi T, Potestio R. Fast, Accurate, and System-Specific Variable-Resolution Modeling of Proteins. J Chem Inf Model 2023; 63:1260-1275. [PMID: 36735551 PMCID: PMC9976289 DOI: 10.1021/acs.jcim.2c01311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Indexed: 02/04/2023]
Abstract
In recent years, a few multiple-resolution modeling strategies have been proposed, in which functionally relevant parts of a biomolecule are described with atomistic resolution, with the remainder of the system being concurrently treated using a coarse-grained model. In most cases, the parametrization of the latter requires lengthy reference all-atom simulations and/or the usage of off-shelf coarse-grained force fields, whose interactions have to be refined to fit the specific system under examination. Here, we overcome these limitations through a novel multiresolution modeling scheme for proteins, dubbed coarse-grained anisotropic network model for variable resolution simulations, or CANVAS. This scheme enables a user-defined modulation of the resolution level throughout the system structure; a fast parametrization of the potential without the necessity of reference simulations; and the straightforward usage of the model on the most commonly used molecular dynamics platforms. The method is presented and validated with two case studies, the enzyme adenylate kinase and the therapeutic antibody pembrolizumab, by comparing the results obtained with the CANVAS model against fully atomistic simulations. The modeling software, implemented in Python, is made freely available for the community on a collaborative github repository.
Collapse
Affiliation(s)
- Raffaele Fiorentini
- Department
of Physics, University of Trento, via Sommarive 14, I-38123 Trento, Italy
- INFN-TIFPA, Trento
Institute for Fundamental Physics and Applications, via Sommarive 14, I-38123 Trento, Italy
| | - Thomas Tarenzi
- Department
of Physics, University of Trento, via Sommarive 14, I-38123 Trento, Italy
- INFN-TIFPA, Trento
Institute for Fundamental Physics and Applications, via Sommarive 14, I-38123 Trento, Italy
| | - Raffaello Potestio
- Department
of Physics, University of Trento, via Sommarive 14, I-38123 Trento, Italy
- INFN-TIFPA, Trento
Institute for Fundamental Physics and Applications, via Sommarive 14, I-38123 Trento, Italy
| |
Collapse
|
5
|
Alipour M, Motavaf M, Abdolmaleki P, Zali A, Ashrafi F, Safari S, Hajipour-Verdom B. Structural Analysis and Conformational Dynamics of Short Helical Hyperphosphorylated Segments of Tau Protein (Sequence 254–290) in Alzheimer’s Disease: A Molecular Dynamics Simulation Study. Front Mol Biosci 2022; 9:884705. [PMID: 36003083 PMCID: PMC9393928 DOI: 10.3389/fmolb.2022.884705] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Accepted: 06/09/2022] [Indexed: 11/21/2022] Open
Abstract
Alzheimer’s disease (AD) is a progressive neurodegenerative disorder whose early diagnosis leads to a chance for successful treatment and decreases the side effects. Hyperphosphorylation of tau proteins is a pathological hallmark of AD that causes it to lose its attachment ability to the microtubules. Alteration of tau structure due to its hyperphosphorylation is an exciting challenge regarding AD treatments. Here, we aimed to examine the structural alterations of short helical segments of tau protein with one to three phosphorylated sites by molecular dynamics simulation. Results indicated that the interaction of two similar segments with three phosphorylated sites (P-Ser262, 285, and 289) formed a compact and more stable structure than the one phosphorylated site complex (P-Ser262). Moreover, due to the high dynamics of the P-Ser262 complex, several structures were made with different conformational dynamics, but there was only one stable cluster of the P-Ser262, 285, and 289 complex during simulation. It seems that the P-Ser262, 285, and 289 complex plays an important role in the formation of paired helical filaments (PHFs) by forming a stable dimer. Generally, it is important to identify how structural features of segments in tau protein change when the phosphorylated sites increase from one to three sites and their effects on the formation of PHFs for drug design and diagnostic biomarkers.
Collapse
Affiliation(s)
- Mozhgan Alipour
- Functional Neurosurgery Research Center, Shohada Tajrish Comprehensive Neurosurgical Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran
- Department of Biophysics, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran
| | - Mahsa Motavaf
- Functional Neurosurgery Research Center, Shohada Tajrish Comprehensive Neurosurgical Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Parviz Abdolmaleki
- Department of Biophysics, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran
| | - Alireza Zali
- Functional Neurosurgery Research Center, Shohada Tajrish Comprehensive Neurosurgical Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Farzad Ashrafi
- Functional Neurosurgery Research Center, Shohada Tajrish Comprehensive Neurosurgical Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Saeid Safari
- Functional Neurosurgery Research Center, Shohada Tajrish Comprehensive Neurosurgical Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran
- *Correspondence: Saeid Safari, ; Behnam Hajipour-Verdom,
| | - Behnam Hajipour-Verdom
- Functional Neurosurgery Research Center, Shohada Tajrish Comprehensive Neurosurgical Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran
- Department of Biophysics, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran
- *Correspondence: Saeid Safari, ; Behnam Hajipour-Verdom,
| |
Collapse
|
6
|
Transformation characteristics of A-DNA in salt solution revealed through molecular dynamics simulations. Biophys Chem 2022; 288:106845. [DOI: 10.1016/j.bpc.2022.106845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 05/25/2022] [Accepted: 06/06/2022] [Indexed: 11/17/2022]
|
7
|
Papež P, Praprotnik M. Dissipative Particle Dynamics Simulation of Ultrasound Propagation through Liquid Water. J Chem Theory Comput 2022; 18:1227-1240. [PMID: 35001631 PMCID: PMC8830050 DOI: 10.1021/acs.jctc.1c01020] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Ultrasound is widely
used as a noninvasive method in therapeutic
and diagnostic applications. These can be further optimized by computational
approaches, as they allow for controlled testing and rational optimization
of the ultrasound parameters, such as frequency and amplitude. Usually,
continuum numerical methods are used to simulate ultrasound propagating
through different tissue types. In contrast, ultrasound simulations
using particle description are less common, as the implementation
is challenging. In this work, a dissipative particle dynamics model
is used to perform ultrasound simulations in liquid water. The effects
of frequency and thermostat parameters are studied and discussed.
We show that frequency and thermostat parameters affect not only the
attenuation but also the computed speed of sound. The present study
paves the way for development and optimization of a virtual ultrasound
machine for large-scale biomolecular simulations.
Collapse
Affiliation(s)
- Petra Papež
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, Ljubljana, SI-1001, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana, SI-1000, Slovenia
| | - Matej Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, Ljubljana, SI-1001, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana, SI-1000, Slovenia
| |
Collapse
|
8
|
Giulini M, Rigoli M, Mattiotti G, Menichetti R, Tarenzi T, Fiorentini R, Potestio R. From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules. Front Mol Biosci 2021; 8:676976. [PMID: 34164432 PMCID: PMC8215203 DOI: 10.3389/fmolb.2021.676976] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 05/06/2021] [Indexed: 12/18/2022] Open
Abstract
The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.
Collapse
Affiliation(s)
- Marco Giulini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Marta Rigoli
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Giovanni Mattiotti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Roberto Menichetti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Thomas Tarenzi
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaele Fiorentini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaello Potestio
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| |
Collapse
|
9
|
Cortes-Huerto R, Praprotnik M, Kremer K, Delle Site L. From adaptive resolution to molecular dynamics of open systems. THE EUROPEAN PHYSICAL JOURNAL. B 2021; 94:189. [PMID: 34720711 PMCID: PMC8547219 DOI: 10.1140/epjb/s10051-021-00193-w] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 08/31/2021] [Indexed: 05/14/2023]
Abstract
ABSTRACT We provide an overview of the Adaptive Resolution Simulation method (AdResS) based on discussing its basic principles and presenting its current numerical and theoretical developments. Examples of applications to systems of interest to soft matter, chemical physics, and condensed matter illustrate the method's advantages and limitations in its practical use and thus settle the challenge for further future numerical and theoretical developments.
Collapse
Affiliation(s)
| | - Matej Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, Ljubljana, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana, Slovenia
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Luigi Delle Site
- Department of Mathematics and Computer Science, Institute for Mathematics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
10
|
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
|
11
|
Soni MN, Kumar SP, S. R. KJ, Rawal RM, Pandya HA. An integrated computational approach to identify GC minor groove binders using various molecular docking scoring functions, dynamics simulations and binding free energy calculations. J Biomol Struct Dyn 2020; 38:3838-3855. [DOI: 10.1080/07391102.2019.1664331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Mehul N. Soni
- Department of Botany, Bioinformatics, and Climate Change Impacts Management, University School of Sciences, Gujarat University, Ahmedabad, India
| | - Sivakumar Prasanth Kumar
- Department of Life Sciences, University School of Sciences, Gujarat University, Ahmedabad, India
| | - Kaid Johar S. R.
- Department of Zoology, Bio-Medical Technology and Human Genetics, University School of Sciences, Gujarat University, Ahmedabad, India
| | - Rakesh M. Rawal
- Department of Life Sciences, University School of Sciences, Gujarat University, Ahmedabad, India
| | - Himanshu A. Pandya
- Department of Botany, Bioinformatics, and Climate Change Impacts Management, University School of Sciences, Gujarat University, Ahmedabad, India
- Department of Life Sciences, University School of Sciences, Gujarat University, Ahmedabad, India
| |
Collapse
|
12
|
Heidari M, Kremer K, Golestanian R, Potestio R, Cortes-Huerto R. Open-boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations. J Chem Phys 2020; 152:194104. [DOI: 10.1063/1.5143268] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Maziar Heidari
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Ramin Golestanian
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), 37077 Göttingen, Germany
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3PU, United Kingdom
| | - Raffaello Potestio
- Physics Department, University of Trento, Via Sommarive, 14, I-38123 Trento, Italy
- INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, I-38123 Trento, Italy
| | | |
Collapse
|
13
|
Delle Site L, Praprotnik M, Bell JB, Klein R. Particle–Continuum Coupling and its Scaling Regimes: Theory and Applications. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.201900232] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Luigi Delle Site
- Freie Universität Berlin Institute of Mathematics Arnimallee 6, 14195 Berlin Germany
| | - Matej Praprotnik
- Laboratory for Molecular Modeling National Institute of Chemistry SI‐1001 Ljubljana, Slovenia & Department of Physics Faculty of Mathematics and Physics University of Ljubljana SI‐1000 Ljubljana Slovenia
| | - John B. Bell
- Lawrence Berkeley National Lab 1 Cyclotron Rd. Berkeley CA 94720 USA
| | - Rupert Klein
- Freie Universität Berlin Institute of Mathematics Arnimallee 6, 14195 Berlin Germany
| |
Collapse
|
14
|
Zavadlav J, Marrink SJ, Praprotnik M. SWINGER: a clustering algorithm for concurrent coupling of atomistic and supramolecular liquids. Interface Focus 2019; 9:20180075. [PMID: 31065343 PMCID: PMC6501350 DOI: 10.1098/rsfs.2018.0075] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/01/2019] [Indexed: 12/11/2022] Open
Abstract
In this contribution, we review recent developments and applications of a dynamic clustering algorithm SWINGER tailored for the multiscale molecular simulations of biomolecular systems. The algorithm on-the-fly redistributes solvent molecules among supramolecular clusters. In particular, we focus on its applications in combination with the adaptive resolution scheme, which concurrently couples atomistic and coarse-grained molecular representations. We showcase the versatility of our multiscale approach on a few applications to biomolecular systems coupling atomistic and supramolecular water models such as the well-established MARTINI and dissipative particle dynamics models and provide an outlook for future work.
Collapse
Affiliation(s)
- Julija Zavadlav
- Computational Science and Engineering Laboratory, ETH-Zurich, Clausiusstrasse 33, 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, 1001 Ljubljana, Slovenia
| |
Collapse
|
15
|
Gunaratne RS, Wilson DB, Flegg MB, Erban R. Multi-resolution dimer models in heat baths with short-range and long-range interactions. Interface Focus 2019; 9:20180070. [PMID: 31065341 PMCID: PMC6501348 DOI: 10.1098/rsfs.2018.0070] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/11/2019] [Indexed: 11/16/2022] Open
Abstract
This work investigates multi-resolution methodologies for simulating dimer models. The solvent particles which make up the heat bath interact with the monomers of the dimer either through direct collisions (short-range) or through harmonic springs (long-range). Two types of multi-resolution methodologies are considered in detail: (a) describing parts of the solvent far away from the dimer by a coarser approach; (b) describing each monomer of the dimer by using a model with different level of resolution. These methodologies are then used to investigate the effect of a shared heat bath versus two uncoupled heat baths, one for each monomer. Furthermore, the validity of the multi-resolution methods is discussed by comparison to dynamics of macroscopic Langevin equations.
Collapse
Affiliation(s)
- Ravinda S. Gunaratne
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Daniel B. Wilson
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Mark B. Flegg
- School of Mathematical Sciences, Monash University, 9 Rainforest walk, Clayton campus, Victoria 3168, Australia
| | - Radek Erban
- Mathematical Institute, University of Oxford, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| |
Collapse
|
16
|
Popadić A, Svenšek D, Podgornik R, Praprotnik M. Density–Nematic Coupling in Isotropic Linear Polymers: Acoustic and Osmotic Birefringence. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201900019] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Aleksandar Popadić
- Laboratory for Molecular ModelingNational Institute of ChemistrySI‐1001 Ljubljana Slovenia
| | - Daniel Svenšek
- Department of Physics, Faculty of Mathematics and PhysicsUniversity of LjubljanaSI‐1000 Ljubljana Slovenia
| | - Rudolf Podgornik
- Department of Physics, Faculty of Mathematics and PhysicsUniversity of LjubljanaSI‐1000 Ljubljana Slovenia
- Department of Theoretical PhysicsJ. Stefan InstituteSI‐1000 Ljubljana Slovenia
- School of Physical Sciences and Kavli Institute for Theoretical SciencesUniversity of Chinese Academy of SciencesBeijing 100049 China
- CAS Key Laboratory of Soft Matter Physics, Institute of PhysicsChinese Academy of SciencesBeijing 100190 China
| | - Matej Praprotnik
- Laboratory for Molecular ModelingNational Institute of ChemistrySI‐1001 Ljubljana Slovenia
| |
Collapse
|
17
|
Tarenzi T, Calandrini V, Potestio R, Carloni P. Open-Boundary Molecular Mechanics/Coarse-Grained Framework for Simulations of Low-Resolution G-Protein-Coupled Receptor-Ligand Complexes. J Chem Theory Comput 2019; 15:2101-2109. [PMID: 30763087 PMCID: PMC6433333 DOI: 10.1021/acs.jctc.9b00040] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2019] [Indexed: 12/18/2022]
Abstract
G-protein-coupled receptors (GPCRs) constitute as much as 30% of the overall proteins targeted by FDA-approved drugs. However, paucity of structural experimental information and low sequence identity between members of the family impair the reliability of traditional docking approaches and atomistic molecular dynamics simulations for in silico pharmacological applications. We present here a dual-resolution approach tailored for such low-resolution models. It couples a hybrid molecular mechanics/coarse-grained (MM/CG) scheme, previously developed by us for GPCR-ligand complexes, with a Hamiltonian-based adaptive resolution scheme (H-AdResS) for the solvent. This dual-resolution approach removes potentially inaccurate atomistic details from the model while building a rigorous statistical ensemble-the grand canonical one-in the high-resolution region. We validate the method on a well-studied GPCR-ligand complex, for which the 3D structure is known, against atomistic simulations. This implementation paves the way for future accurate in silico studies of low-resolution ligand/GPCRs models.
Collapse
Affiliation(s)
- Thomas Tarenzi
- Computation-based Science and Technology Research Center CaSToRC , The Cyprus Institute , 20 Konstaninou Kavafi Street , 2121 Aglantzia, Nicosia , Cyprus
- Departments of Physics , Faculty of Mathematics, Computer Science and Natural Sciences, Aachen University , Otto-Blumenthal Straße , 52062 Aachen , Germany
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9 , Forschungszentrum Jülich , 52428 Jülich , Germany
| | - Vania Calandrini
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9 , Forschungszentrum Jülich , 52428 Jülich , Germany
| | - Raffaello Potestio
- Department of Physics , University of Trento , via Sommarive 14 Povo , Trento 38123 , Italy
- INFN-TIFPA, Trento Institute for Fundamental Physics and Applications , I-38123 Trento , Italy
| | - Paolo Carloni
- Departments of Physics , Faculty of Mathematics, Computer Science and Natural Sciences, Aachen University , Otto-Blumenthal Straße , 52062 Aachen , Germany
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9 , Forschungszentrum Jülich , 52428 Jülich , Germany
- JARA-HPC, Jülich Supercomputing Center , Forschungszentrum Jülich , 52428 Jülich , Germany
| |
Collapse
|
18
|
Ciccotti G, Delle Site L. The physics of open systems for the simulation of complex molecular environments in soft matter. SOFT MATTER 2019; 15:2114-2124. [PMID: 30761396 DOI: 10.1039/c8sm02523a] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Molecular dynamics (MD) has become one of the most powerful tools of investigation in soft matter. Despite such success, simulations of large molecular environments are mostly run using the approximation of closed systems without the possibility of exchange of matter. Due to the molecular complexity of soft matter systems, an optimal simulation strategy would require the application of concurrent multiscale resolution approaches such that each part of a large system can be considered as an open subsystem at a high resolution embedded in a large coarser reservoir of energy and particles. This paper discusses the current capability and the future perspectives of multiscale adaptive resolution MD methods to satisfy the conceptual principles of open systems and to perform simulations of complex molecular environments in soft matter.
Collapse
Affiliation(s)
- Giovanni Ciccotti
- Instituto per le Applicazioni del Calcolo, Consiglio Nazionale delle Ricerche, and Universita' La Sapienza, Rome, Italy.
| | | |
Collapse
|
19
|
Friedman R. Simulations of Biomolecules in Electrolyte Solutions. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201800163] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ran Friedman
- Department of Chemistry and Biomedical SciencesLinnæus UniversityKalmar SE‐391 82 Sweden
| |
Collapse
|
20
|
Krekeler C, Agarwal A, Junghans C, Praprotnik M, Delle Site L. Adaptive resolution molecular dynamics technique: Down to the essential. J Chem Phys 2018; 149:024104. [PMID: 30007373 DOI: 10.1063/1.5031206] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
We investigate the role of the thermodynamic (TD) force as an essential and sufficient technical ingredient for an efficient and accurate adaptive resolution algorithm. Such a force applied in the coupling region of an adaptive resolution molecular dynamics setup assures thermodynamic equilibrium between atomistically resolved and coarse-grained regions, allowing the proper exchange of molecules. We numerically prove that indeed for systems as relevant as liquid water and 1,3-dimethylimidazolium chloride ionic liquid, the combined action of the TD force and thermostat allows for computationally efficient and numerically accurate simulations, beyond the current capabilities of adaptive resolution setups, which employ switching functions in the coupling region.
Collapse
Affiliation(s)
| | - Animesh Agarwal
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Christoph Junghans
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Matej Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1001 Ljubljana, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, SI-1000 Ljubljana, Slovenia
| | - Luigi Delle Site
- Institute for Mathematics, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|