1
|
Lu D, Luo D, Zhang Y, Wang B. A Robust Induced Fit Docking Approach with the Combination of the Hybrid All-Atom/United-Atom/Coarse-Grained Model and Simulated Annealing. J Chem Theory Comput 2024; 20:6414-6423. [PMID: 38966989 DOI: 10.1021/acs.jctc.4c00653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/06/2024]
Abstract
Molecular docking remains an indispensable tool in computational biology and structure-based drug discovery. However, the correct prediction of binding poses remains a major challenge for molecular docking, especially for target proteins where a substrate binding induces significant reorganization of the active site. Here, we introduce an Induced Fit Docking (IFD) approach named AA/UA/CG-SA-IFD, which combines a hybrid All-Atom/United-Atom/Coarse-Grained model with Simulated Annealing. In this approach, the core region is represented by the All-Atom(AA) model, while the protein environment beyond the core region and the solvent are treated with either the United-Atom (UA) or the Coarse-Grained (CG) model. By combining the Elastic Network Model (ENM) for the CG region, the hybrid model ensures a reasonable description of ligand binding and the environmental effects of the protein, facilitating highly efficient and reliable sampling of ligand binding through Simulated Annealing (SA) at a high temperature. Upon validation with two testing sets, the AA/UA/CG-SA-IFD approach demonstrates remarkable accuracy and efficiency in induced fit docking, even for challenging cases where the docked poses significantly deviate from crystal structures.
Collapse
Affiliation(s)
- Dexin Lu
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Ding Luo
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| | - Yuwei Zhang
- Jiangsu Key Laboratory of New Power Batteries, Jiangsu Collaborative Innovation Centre of Biomedical Functional Materials, School of Chemistry and Materials Science, Nanjing Normal University, Nanjing 210023, P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, P. R. China
| |
Collapse
|
2
|
Li S, Wu B, Luo YL, Han W. Simulations of Functional Motions of Super Large Biomolecules with a Mixed-Resolution Model. J Chem Theory Comput 2024; 20:2228-2245. [PMID: 38374639 PMCID: PMC10938502 DOI: 10.1021/acs.jctc.3c01046] [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: 09/22/2023] [Revised: 01/18/2024] [Accepted: 01/29/2024] [Indexed: 02/21/2024]
Abstract
Many large protein machines function through an interplay between large-scale movements and intricate conformational changes. Understanding functional motions of these proteins through simulations becomes challenging for both all-atom and coarse-grained (CG) modeling techniques because neither approach alone can readily capture the full details of these motions. In this study, we develop a multiscale model by employing the popular MARTINI CG model to represent a heterogeneous environment and structurally stable proteins and using the united-atom (UA) model PACE to describe proteins undergoing subtle conformational changes. PACE was previously developed to be compatible with the MARTINI solvent and membrane. Here, we couple the protein descriptions of the two models by directly mixing UA and CG interaction parameters to greatly simplify parameter determination. Through extensive validations with diverse protein systems in solution or membrane, we demonstrate that only additional parameter rescaling is needed to enable the resulting model to recover the stability of native structures of proteins under mixed representation. Moreover, we identify the optimal scaling factors that can be applied to various protein systems, rendering the model potentially transferable. To further demonstrate its applicability for realistic systems, we apply the model to a mechanosensitive ion channel Piezo1 that has peripheral arms for sensing membrane tension and a central pore for ion conductance. The model can reproduce the coupling between Piezo1's large-scale arm movement and subtle pore opening in response to membrane stress while consuming much less computational costs than all-atom models. Therefore, our model shows promise for studying functional motions of large protein machines.
Collapse
Affiliation(s)
- Shu Li
- Centre
for Artificial Intelligence Driven Drug Discovery, Faculty of Applied
Sciences, Macao Polytechnic University, Macao 999078, China
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Bohua Wu
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Yun Lyna Luo
- Department
of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| | - Wei Han
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Department
of Chemistry, Faculty of Science, Hong Kong
Baptist University, Hong Kong SAR 999077, China
- Shenzhen
Bay Laboratory, Institute of Chemical Biology, Shenzhen 518132, China
| |
Collapse
|
3
|
Ge Y, Wang X, Zhu Q, Yang Y, Dong H, Ma J. Machine Learning-Guided Adaptive Parametrization for Coupling Terms in a Mixed United-Atom/Coarse-Grained Model for Diphenylalanine Self-Assembly in Aqueous Ionic Liquids. J Chem Theory Comput 2023; 19:6718-6732. [PMID: 37725682 DOI: 10.1021/acs.jctc.3c00809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
Precise regulation of the peptide self-assembly into ordered nanostructures with intriguing properties has attracted intense attention. However, predicting peptide assembly at atomic resolution is a challenge due to both the structural flexibility of peptides and the associated huge computational costs. A machine learning-guided adaptive parametrization method was proposed for developing a mixed atomic and coarse-grained (CG) model through a multiobjective optimization strategy. Our model incorporates the united-atom (UA) model for diphenylalanine (P) and the polarizable electrostatic-variable coarse-grained (VaCG) model for aqueous ionic liquid [BMIM]+[BF4]- solution. In this mixed model, the coupling van der Waals (vdW) interaction is addressed by introducing virtual sites (VS) in the UA model to interact with solvent CG beads. The coupling parameters, including the electrostatic parameter and vdW parameters, are automatically optimized through ML-guided adaptive parametrization. The performance of this model was tested by some microstructural properties, e.g., the average number of P-P intermolecular hydrogen bonds (HBs) and radius distribution functions (RDFs) between P and different fragments of IL, in comparison with all-atom (AA) simulations. The computational cost is significantly reduced using such a parametrization scheme, which could search tens of thousands of force-field parameter sets, while needing only a small fraction of them to be assessed with molecular dynamics (MD) simulations. We used such a mixed resolution model to investigate the self-assembly in IL-water mixtures with variants of IL concentration (X). The long-range-ordered fibril structure is formed in a pure water system (X = 0). With an increase of IL concentrations, the formation of an ordered self-assembly nanostructure is prohibited, instead forming branched fibril at X = 2 mol % or amorphous aggregates when X > 10 mol %, resulting from the interplay between π-stacking and HB interactions between P and IL. The qualitative agreement between the simulated structures and the observed morphologies in experiments indicates the applicability of ML-guided parametrization strategy in the study of complex systems, such as polymers, lipid bilayers, and polysaccharides.
Collapse
Affiliation(s)
- Yang Ge
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Xueping Wang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Qiang Zhu
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yuqin Yang
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
| | - Hao Dong
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China
- State Key Laboratory of Analytical Chemistry for Life Science, Nanjing University, Nanjing 210023, China
- Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| | - Jing Ma
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
4
|
Bag S, Meinel MK, Müller-Plathe F. Toward a Mobility-Preserving Coarse-Grained Model: A Data-Driven Approach. J Chem Theory Comput 2022; 18:7108-7120. [PMID: 36449362 DOI: 10.1021/acs.jctc.2c00898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
Coarse-grained molecular dynamics (MD) simulation is a promising alternative to all-atom MD simulation for the fast calculation of system properties, which is imperative in designing materials with a specific target property. There have been several coarse-graining strategies developed over the past few years that provide accurate structural properties of the system. However, these coarse-grained models share a major drawback in that they introduce an artificial acceleration in molecular mobility. In this paper, we report a data-driven approach to generate coarse-grained models that preserve the all-atom molecular mobility. We designed a machine learning model in the form of an artificial neural network, which directly predicts the simulation-ready mobility-preserving coarse-grained potential as an output given the all-atom force field (FF) parameters as inputs. As a proof of principle, we took 2,3,4-trimethylpentane as a model system and described the development of machine learning models in detail. We quantify the artificial acceleration in molecular mobility by defining the acceleration factor as the ratio of the coarse-grained and the all-atom diffusion coefficient. The predicted coarse-grained potential generated by the best machine learning model can bring down the acceleration factor to a value of ∼2, which could be otherwise as large as 7 for a typical value of 3 × 10-9 m2 s-1 for the all-atom diffusion coefficient. We believe our method will be of interest in the community as a route to generating coarse-grained potentials with accurate dynamics.
Collapse
Affiliation(s)
- Saientan Bag
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287Darmstadt, Germany
| | - Melissa K Meinel
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287Darmstadt, Germany
| |
Collapse
|
5
|
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
|
6
|
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
|
7
|
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
|
8
|
Meinel MK, Müller-Plathe F. Loss of Molecular Roughness upon Coarse-Graining Predicts the Artificially Accelerated Mobility of Coarse-Grained Molecular Simulation Models. J Chem Theory Comput 2020; 16:1411-1419. [DOI: 10.1021/acs.jctc.9b00943] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Melissa K. Meinel
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Profile Area Thermofluids and Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Strasse 8, D-64287 Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Profile Area Thermofluids and Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Strasse 8, D-64287 Darmstadt, Germany
| |
Collapse
|
9
|
Quantum mechanics/molecular mechanics multiscale modeling of biomolecules. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2020. [DOI: 10.1016/bs.apoc.2020.08.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
10
|
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
|
11
|
Webb MA, Delannoy JY, de Pablo JJ. Graph-Based Approach to Systematic Molecular Coarse-Graining. J Chem Theory Comput 2018; 15:1199-1208. [PMID: 30557028 DOI: 10.1021/acs.jctc.8b00920] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A novel methodology is introduced here to generate coarse-grained (CG) representations of molecular models for simulations. The proposed strategy relies on basic graph-theoretic principles and is referred to as graph-based coarse-graining (GBCG). It treats a given system as a molecular graph and derives a corresponding CG representation by using edge contractions to combine nodes in the graph, which correspond to atoms in the molecule, into CG sites. A key element of this methodology is that the nodes are combined according to well-defined protocols that rank-order nodes based on the underlying chemical connectivity. By iteratively performing these operations, successively coarser representations of the original atomic system can be produced to yield a systematic set of CG mappings with hierarchical resolution in an automated fashion. These capabilities are demonstrated in the context of several test systems, including toluene, pentadecane, a polysaccharide dimer, and a rhodopsin protein. In these examples, GBCG yields multiple, intuitive structures that naturally preserve the chemical topology of the system. Importantly, these representations are rendered from algorithmic implementation rather than an arbitrary ansatz, which, until now, has been the conventional approach for defining CG mapping schemes. Overall, the results presented here indicate that GBCG is efficient, robust, and unambiguous in its application, making it a valuable tool for future CG modeling.
Collapse
Affiliation(s)
- Michael A Webb
- Institute for Molecular Engineering , University of Chicago , Chicago , Illinois 60637 , United States
| | - Jean-Yves Delannoy
- Simulation, Modeling and Artificial Intelligence team , Solvay , Bristol , Pennsylvania 19007 , United States
| | - Juan J de Pablo
- Institute for Molecular Engineering , University of Chicago , Chicago , Illinois 60637 , United States.,Institute for Molecular Engineering and Materials Science Division , Argonne National Laboratory , Lemont , Illinois 60439 , United States
| |
Collapse
|
12
|
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
|
13
|
Kar P, Feig M. Hybrid All-Atom/Coarse-Grained Simulations of Proteins by Direct Coupling of CHARMM and PRIMO Force Fields. J Chem Theory Comput 2017; 13:5753-5765. [PMID: 28992696 DOI: 10.1021/acs.jctc.7b00840] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Hybrid all-atom/coarse-grained (AA/CG) simulations of proteins offer a computationally efficient compromise where atomistic details are only applied to biologically relevant regions while benefiting from the speedup of treating the remaining parts of a given system at the CG level. The recently developed CG model, PRIMO, allows a direct coupling with an atomistic force field with no additional modifications or coupling terms and the ability to carry out dynamic simulations without any restraints on secondary or tertiary structures. A hybrid AA/CG scheme based on combining all-atom CHARMM and coarse-grained PRIMO representations was validated via molecular dynamics and replica exchange simulations of soluble and membrane proteins. The AA/CG scheme was also tested in the calculation of the free energy profile for the transition from the closed to the open state of adenylate kinase via umbrella sampling molecular dynamics method. The overall finding is that the AA/CG scheme generates dynamics and energetics that are qualitatively and quantitatively comparable to AA simulations while offering the computational advantages of coarse-graining. This model opens the door to challenging applications where high accuracy is required only in parts of large biomolecular complexes.
Collapse
Affiliation(s)
- Parimal Kar
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University , East Lansing, Michigan 48824, United States
| |
Collapse
|
14
|
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
|
15
|
Charchar P, Christofferson AJ, Todorova N, Yarovsky I. Understanding and Designing the Gold-Bio Interface: Insights from Simulations. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2016; 12:2395-418. [PMID: 27007031 DOI: 10.1002/smll.201503585] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2015] [Revised: 02/01/2016] [Indexed: 05/20/2023]
Abstract
Gold nanoparticles (AuNPs) are an integral part of many exciting and novel biomedical applications, sparking the urgent need for a thorough understanding of the physicochemical interactions occurring between these inorganic materials, their functional layers, and the biological species they interact with. Computational approaches are instrumental in providing the necessary molecular insight into the structural and dynamic behavior of the Au-bio interface with spatial and temporal resolutions not yet achievable in the laboratory, and are able to facilitate a rational approach to AuNP design for specific applications. A perspective of the current successes and challenges associated with the multiscale computational treatment of Au-bio interfacial systems, from electronic structure calculations to force field methods, is provided to illustrate the links between different approaches and their relationship to experiment and applications.
Collapse
Affiliation(s)
- Patrick Charchar
- School of Engineering, RMIT University, Melbourne, Victoria, 3001, Australia
| | | | - Nevena Todorova
- School of Engineering, RMIT University, Melbourne, Victoria, 3001, Australia
| | - Irene Yarovsky
- School of Engineering, RMIT University, Melbourne, Victoria, 3001, Australia
| |
Collapse
|
16
|
Stanzione F, Jayaraman A. Hybrid Atomistic and Coarse-Grained Molecular Dynamics Simulations of Polyethylene Glycol (PEG) in Explicit Water. J Phys Chem B 2016; 120:4160-73. [DOI: 10.1021/acs.jpcb.6b02327] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Francesca Stanzione
- Department
of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States
| | - Arthi Jayaraman
- Department
of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States
- Department
of Materials Science and Engineering, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
17
|
Szklarczyk OM, Bieler NS, Hünenberger PH, van Gunsteren WF. Flexible Boundaries for Multiresolution Solvation: An Algorithm for Spatial Multiscaling in Molecular Dynamics Simulations. J Chem Theory Comput 2015; 11:5447-63. [PMID: 26574333 DOI: 10.1021/acs.jctc.5b00406] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An algorithm is proposed for performing molecular dynamics (MD) simulations of a biomolecular solute represented at atomistic resolution surrounded by a surface layer of atomistic fine-grained (FG) solvent molecules within a bulk represented by coarse-grained (CG) solvent beads. The method, called flexible boundaries for multiresolution solvation (FBMS), is based on: (i) a three-region layering of the solvent around the solute, involving an FG layer surrounded by a mixed FG-CG buffer layer, itself surrounded by a bulk CG region; (ii) a definition of the layer boundary that relies on an effective distance to the solute surface and is thus adapted to the shape of the solute as well as adjusts to its conformational changes. The effective surface distance is defined by inverse-nth power averaging over the distances to all non-hydrogen solute atoms, and the layering is enforced by means of half-harmonic distance restraints, attractive for the FG molecules and repulsive for the CG beads. A restraint-free region at intermediate distances enables the formation of the buffer layer, where the FG and CG solvents can mix freely. The algorithm is tested and validated using the GROMOS force field and the associated FG (SPC) and CG (polarizable CGW) water models. The test systems include pure-water systems, where one FG molecule plays the role of a solute, and a deca-alanine peptide with two widely different solute shapes considered, α-helical and fully extended. In particular, as the peptide unfolds, the number of FG molecules required to fill its close-range solvation layer increases, with the additional molecules being provided by the buffer layer. Further validation involves simulations of four proteins in multiresolution FG/CG mixtures. The resulting structural, energetic, and solvation properties are found to be similar to those observed in corresponding pure FG simulations.
Collapse
Affiliation(s)
- Oliwia M Szklarczyk
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, Swiss Federal Institute of Technology ETH , 8093 Zürich, Switzerland
| | - Noah S Bieler
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, Swiss Federal Institute of Technology ETH , 8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, Swiss Federal Institute of Technology ETH , 8093 Zürich, Switzerland
| | - Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, Swiss Federal Institute of Technology ETH , 8093 Zürich, Switzerland
| |
Collapse
|
18
|
Zavadlav J, Podgornik R, Praprotnik M. Adaptive Resolution Simulation of a DNA Molecule in Salt Solution. J Chem Theory Comput 2015; 11:5035-44. [DOI: 10.1021/acs.jctc.5b00596] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Julija Zavadlav
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Hajdrihova
19, SI-1001 Ljubljana, Slovenia
| | - Rudolf Podgornik
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, SI-1000 Ljubljana, Slovenia
- Theoretical
Physics Department, J. Stefan Institute, Jamova c. 39, SI-1000 Ljubljana, Slovenia
| | - Matej Praprotnik
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Hajdrihova
19, SI-1001 Ljubljana, Slovenia
| |
Collapse
|
19
|
Kuhn AB, Gopal SM, Schäfer LV. On Using Atomistic Solvent Layers in Hybrid All-Atom/Coarse-Grained Molecular Dynamics Simulations. J Chem Theory Comput 2015; 11:4460-72. [DOI: 10.1021/acs.jctc.5b00499] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alexander B. Kuhn
- Lehrstuhl für
Theoretische
Chemie, Ruhr-University Bochum, D-44780 Bochum, Germany
| | - Srinivasa M. Gopal
- Lehrstuhl für
Theoretische
Chemie, Ruhr-University Bochum, D-44780 Bochum, Germany
| | - Lars V. Schäfer
- Lehrstuhl für
Theoretische
Chemie, Ruhr-University Bochum, D-44780 Bochum, Germany
| |
Collapse
|
20
|
Sokkar P, Boulanger E, Thiel W, Sanchez-Garcia E. Hybrid Quantum Mechanics/Molecular Mechanics/Coarse Grained Modeling: A Triple-Resolution Approach for Biomolecular Systems. J Chem Theory Comput 2015; 11:1809-18. [DOI: 10.1021/ct500956u] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Pandian Sokkar
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Eliot Boulanger
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| | - Elsa Sanchez-Garcia
- Max-Planck-Institut für
Kohlenforschung, Kaiser-Wilhelm-Platz
1, 45470 Mülheim
an der Ruhr, Germany
| |
Collapse
|
21
|
Spiriti J, Zuckerman DM. Tunable Coarse Graining for Monte Carlo Simulations of Proteins via Smoothed Energy Tables: Direct and Exchange Simulations. J Chem Theory Comput 2014; 10:5161-5177. [PMID: 25400525 PMCID: PMC4230378 DOI: 10.1021/ct500622z] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2014] [Indexed: 12/03/2022]
Abstract
Many commonly used coarse-grained models for proteins are based on simplified interaction sites and consequently may suffer from significant limitations, such as the inability to properly model protein secondary structure without the addition of restraints. Recent work on a benzene fluid (Lettieri S.; Zuckerman D. M.J. Comput. Chem.2012, 33, 268-275) suggested an alternative strategy of tabulating and smoothing fully atomistic orientation-dependent interactions among rigid molecules or fragments. Here we report our initial efforts to apply this approach to the polar and covalent interactions intrinsic to polypeptides. We divide proteins into nearly rigid fragments, construct distance and orientation-dependent tables of the atomistic interaction energies between those fragments, and apply potential energy smoothing techniques to those tables. The amount of smoothing can be adjusted to give coarse-grained models that range from the underlying atomistic force field all the way to a bead-like coarse-grained model. For a moderate amount of smoothing, the method is able to preserve about 70-90% of the α-helical structure while providing a factor of 3-10 improvement in sampling per unit computation time (depending on how sampling is measured). For a greater amount of smoothing, multiple folding-unfolding transitions of the peptide were observed, along with a factor of 10-100 improvement in sampling per unit computation time, although the time spent in the unfolded state was increased compared with less smoothed simulations. For a β hairpin, secondary structure is also preserved, albeit for a narrower range of the smoothing parameter and, consequently, for a more modest improvement in sampling. We have also applied the new method in a "resolution exchange" setting, in which each replica runs a Monte Carlo simulation with a different degree of smoothing. We obtain exchange rates that compare favorably to our previous efforts at resolution exchange (Lyman E.; Zuckerman D. M.J. Chem. Theory Comput.2006, 2, 656-666).
Collapse
Affiliation(s)
- Justin Spiriti
- Department of Computational
and Systems Biology, University of Pittsburgh, Pittsburgh, Pennsylvania 15213, United States
| | - Daniel M. Zuckerman
- Department of Computational
and Systems Biology, University of Pittsburgh, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
22
|
Kar P, Feig M. Recent advances in transferable coarse-grained modeling of proteins. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2014; 96:143-80. [PMID: 25443957 PMCID: PMC5366245 DOI: 10.1016/bs.apcsb.2014.06.005] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Computer simulations are indispensable tools for studying the structure and dynamics of biological macromolecules. Biochemical processes occur on different scales of length and time. Atomistic simulations cannot cover the relevant spatiotemporal scales at which the cellular processes occur. To address this challenge, coarse-grained (CG) modeling of the biological systems is employed. Over the last few years, many CG models for proteins continue to be developed. However, many of them are not transferable with respect to different systems and different environments. In this review, we discuss those CG protein models that are transferable and that retain chemical specificity. We restrict ourselves to CG models of soluble proteins only. We also briefly review recent progress made in the multiscale hybrid all-atom/CG simulations of proteins.
Collapse
Affiliation(s)
- Parimal Kar
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA
| | - Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan, USA; Department of Chemistry, Michigan State University, East Lansing, Michigan, USA.
| |
Collapse
|
23
|
Pezeshki S, Lin H. Recent developments in QM/MM methods towards open-boundary multi-scale simulations. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.911870] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
24
|
Zavadlav J, Melo MN, Cunha AV, de Vries AH, Marrink SJ, Praprotnik M. Adaptive Resolution Simulation of MARTINI Solvents. J Chem Theory Comput 2014; 10:2591-8. [DOI: 10.1021/ct5001523] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Julija Zavadlav
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Hajdrihova
19, SI-1001 Ljubljana, Slovenia
| | - Manuel N. Melo
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Ana V. Cunha
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Alex H. de Vries
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Matej Praprotnik
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Hajdrihova
19, SI-1001 Ljubljana, Slovenia
| |
Collapse
|
25
|
Zavadlav J, Melo MN, Marrink SJ, Praprotnik M. Adaptive resolution simulation of an atomistic protein in MARTINI water. J Chem Phys 2014; 140:054114. [DOI: 10.1063/1.4863329] [Citation(s) in RCA: 69] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
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
|