1
|
Pirnia A, Maqdisi R, Mittal S, Sener M, Singharoy A. Perspective on Integrative Simulations of Bioenergetic Domains. J Phys Chem B 2024; 128:3302-3319. [PMID: 38562105 DOI: 10.1021/acs.jpcb.3c07335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.
Collapse
Affiliation(s)
- Adam Pirnia
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Ranel Maqdisi
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Sumit Mittal
- VIT Bhopal University, Sehore 466114, Madhya Pradesh, India
| | - Melih Sener
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| |
Collapse
|
2
|
Zanardi I, Venturi S, Panesi M. Adaptive physics-informed neural operator for coarse-grained non-equilibrium flows. Sci Rep 2023; 13:15497. [PMID: 37726349 PMCID: PMC10509218 DOI: 10.1038/s41598-023-41039-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Accepted: 08/21/2023] [Indexed: 09/21/2023] Open
Abstract
This work proposes a new machine learning (ML)-based paradigm aiming to enhance the computational efficiency of non-equilibrium reacting flow simulations while ensuring compliance with the underlying physics. The framework combines dimensionality reduction and neural operators through a hierarchical and adaptive deep learning strategy to learn the solution of multi-scale coarse-grained governing equations for chemical kinetics. The proposed surrogate's architecture is structured as a tree, with leaf nodes representing separate neural operator blocks where physics is embedded in the form of multiple soft and hard constraints. The hierarchical attribute has two advantages: (i) It allows the simplification of the training phase via transfer learning, starting from the slowest temporal scales; (ii) It accelerates the prediction step by enabling adaptivity as the surrogate's evaluation is limited to the necessary leaf nodes based on the local degree of non-equilibrium of the gas. The model is applied to the study of chemical kinetics relevant for application to hypersonic flight, and it is tested here on pure oxygen gas mixtures. In 0-[Formula: see text] scenarios, the proposed ML framework can adaptively predict the dynamics of almost thirty species with a maximum relative error of 4.5% for a wide range of initial conditions. Furthermore, when employed in 1-[Formula: see text] shock simulations, the approach shows accuracy ranging from 1% to 4.5% and a speedup of one order of magnitude compared to conventional implicit schemes employed in an operator-splitting integration framework. Given the results presented in the paper, this work lays the foundation for constructing an efficient ML-based surrogate coupled with reactive Navier-Stokes solvers for accurately characterizing non-equilibrium phenomena in multi-dimensional computational fluid dynamics simulations.
Collapse
Affiliation(s)
- Ivan Zanardi
- Center for Hypersonics and Entry Systems Studies, Department of Aerospace Engineering, University of Illinois Urbana-Champaign, Urbana, 61801, IL, USA
| | - Simone Venturi
- Center for Hypersonics and Entry Systems Studies, Department of Aerospace Engineering, University of Illinois Urbana-Champaign, Urbana, 61801, IL, USA
| | - Marco Panesi
- Center for Hypersonics and Entry Systems Studies, Department of Aerospace Engineering, University of Illinois Urbana-Champaign, Urbana, 61801, IL, USA.
| |
Collapse
|
3
|
Polymer brushes for friction control: Contributions of molecular simulations. Biointerphases 2023; 18:010801. [PMID: 36653299 DOI: 10.1116/6.0002310] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
When polymer chains are grafted to solid surfaces at sufficiently high density, they form brushes that can modify the surface properties. In particular, polymer brushes are increasingly being used to reduce friction in water-lubricated systems close to the very low levels found in natural systems, such as synovial joints. New types of polymer brush are continually being developed to improve with lower friction and adhesion, as well as higher load-bearing capacities. To complement experimental studies, molecular simulations are increasingly being used to help to understand how polymer brushes reduce friction. In this paper, we review how molecular simulations of polymer brush friction have progressed from very simple coarse-grained models toward more detailed models that can capture the effects of brush topology and chemistry as well as electrostatic interactions for polyelectrolyte brushes. We pay particular attention to studies that have attempted to match experimental friction data of polymer brush bilayers to results obtained using molecular simulations. We also critically look at the remaining challenges and key limitations to overcome and propose future modifications that could potentially improve agreement with experimental studies, thus enabling molecular simulations to be used predictively to modify the brush structure for optimal friction reduction.
Collapse
|
4
|
Zhang Y, Oestreicher J, Binns WJ, Briggs S, Kim CS, Béland LK. A Coarse-Grained Interaction Model for Sodium Dominant Montmorillonite. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:13226-13237. [PMID: 36256513 DOI: 10.1021/acs.langmuir.2c02233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Montmorillonite is the main crystalline mineral present in bentonite. It is an absorbent, swelling material; the physical chemistry underlying its ability to absorb water and swell occurs at the nanoscale, governed by electrical double-layer interactions. In turn, absorption and swelling lead to important changes in the macroscopic transport properties of the clay. Mesoscale models can help us establish a link between these nanoscale processes and macroscale properties, notably by providing a detailed description of its pore network. Models on the scale of hundreds to thousands of nanometers are required, which cannot realistically be handled using traditional all-atom molecular dynamics simulations. This work presents a coarse-grained (CG) mesoscale model of sodium montmorillonite. In our model, montmorillonite platelets are represented by two types of particles: central nonhydrogen-bonded particles and edge hydrogen-bonding particles. The particle interactions are described by two-body potentials, which were optimized based on all-atom molecular dynamics simulations. Specifically, several potential mean force calculations involving dry and hydrated montmorillonite were performed, using the ClayFF potential to calculate interatomic forces. The CG model was validated by testing the scalability of the model, testing its ability to reproduce potentials of mean force reported elsewhere in the literature, and by comparing the calculated elastic properties of a system containing 1000 Na montmorillonite platelets to experimentally measured elastic properties of bentonite. The simulated elastic properties obtained using our mesoscale model agree with these experimental values.
Collapse
Affiliation(s)
- Yaoting Zhang
- Department of Mechanical & Materials Engineering, Queen's University, Nicol Hall, 60 Union Street, Kingston, K7L 3N6Ontario, Canada
- Nuclear Waste Management Organization, 22 St. Clair Avenue East Fourth Floor, Toronto, M4T 2S3Ontario, Canada
| | - Jeremy Oestreicher
- Department of Mechanical & Materials Engineering, Queen's University, Nicol Hall, 60 Union Street, Kingston, K7L 3N6Ontario, Canada
| | - W Jeffrey Binns
- Nuclear Waste Management Organization, 22 St. Clair Avenue East Fourth Floor, Toronto, M4T 2S3Ontario, Canada
| | - Scott Briggs
- Nuclear Waste Management Organization, 22 St. Clair Avenue East Fourth Floor, Toronto, M4T 2S3Ontario, Canada
| | - Chang Seok Kim
- Nuclear Waste Management Organization, 22 St. Clair Avenue East Fourth Floor, Toronto, M4T 2S3Ontario, Canada
| | - Laurent K Béland
- Department of Mechanical & Materials Engineering, Queen's University, Nicol Hall, 60 Union Street, Kingston, K7L 3N6Ontario, Canada
| |
Collapse
|
5
|
Styles KM, Brown AT, Sagona AP. A Review of Using Mathematical Modeling to Improve Our Understanding of Bacteriophage, Bacteria, and Eukaryotic Interactions. Front Microbiol 2021; 12:724767. [PMID: 34621252 PMCID: PMC8490754 DOI: 10.3389/fmicb.2021.724767] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 08/27/2021] [Indexed: 12/27/2022] Open
Abstract
Phage therapy, the therapeutic usage of viruses to treat bacterial infections, has many theoretical benefits in the ‘post antibiotic era.’ Nevertheless, there are currently no approved mainstream phage therapies. One reason for this is a lack of understanding of the complex interactions between bacteriophage, bacteria and eukaryotic hosts. These three-component interactions are complex, with non-linear or synergistic relationships, anatomical barriers and genetic or phenotypic heterogeneity all leading to disparity between performance and efficacy in in vivo versus in vitro environments. Realistic computer or mathematical models of these complex environments are a potential route to improve the predictive power of in vitro studies for the in vivo environment, and to streamline lab work. Here, we introduce and review the current status of mathematical modeling and highlight that data on genetic heterogeneity and mutational stochasticity, time delays and population densities could be critical in the development of realistic phage therapy models in the future. With this in mind, we aim to inform and encourage the collaboration and sharing of knowledge and expertise between microbiologists and theoretical modelers, synergising skills and smoothing the road to regulatory approval and widespread use of phage therapy.
Collapse
Affiliation(s)
- Kathryn M Styles
- School of Life Sciences, University of Warwick, Coventry, United Kingdom
| | - Aidan T Brown
- School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom
| | - Antonia P Sagona
- School of Life Sciences, University of Warwick, Coventry, United Kingdom
| |
Collapse
|
6
|
Shen H, Wu Z, Lu C. Extension of the CAVS model to the simulation of helical peptides in a membrane environment. Phys Chem Chem Phys 2021; 23:12850-12863. [PMID: 34060579 DOI: 10.1039/d1cp01238j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Considering the effect of peptide insertion on the dipole potential of the lipid membrane, we extend the CAVS coarse-grained (CG) model to the simulation of helical peptides in a membrane environment. In this approach, the CG scheme for a peptide backbone is similar to the treatment in the united-atom model, while we treated the side chain of an amino acid by grouping 1-3 heavy atoms into a CG unit. The CAVS CG force field for peptides is optimized by reproducing the experimental results for the backbone (φ, ψ) distribution and predicting the PMF profiles of transferring organic molecules in a lipid bilayer membrane obtained from all-atom simulations. The CAVS simulation of a helical peptide in a phosphatidylcholine (PC) lipid bilayer revealed that the insertion of a peptide increases the dipole potential of the PC lipid bilayer, in which the peptide and its neutralized ions make a significant contribution. Finally, we carried out the CAVS simulation for five different helical peptides in the PC lipid bilayer to explore the behavior of peptide tilt, showing excellent agreement with the all-atom simulations. Our work suggests that the peptide tilt should relieve the deformation stress from the lipid bilayer, and the peptide aggregation could reduce the peptide tilt by resisting the deformation stress from the surrounding lipids.
Collapse
Affiliation(s)
- Hujun Shen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Education University, Guiyang, Guizhou 550018, P. R. China.
| | - Zhenhua Wu
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Chan Lu
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Education University, Guiyang, Guizhou 550018, P. R. China.
| |
Collapse
|
7
|
Maximiano P, Durães L, Simões P. Overview of Multiscale Molecular Modeling and Simulation of Silica Aerogels. Ind Eng Chem Res 2019. [DOI: 10.1021/acs.iecr.9b03781] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Pedro Maximiano
- CIEPQPF, Department of Chemical Engineering, University of Coimbra, Rua Sı́lvio de Lima, 3030-790 Coimbra, Portugal
| | - Luísa Durães
- CIEPQPF, Department of Chemical Engineering, University of Coimbra, Rua Sı́lvio de Lima, 3030-790 Coimbra, Portugal
| | - Pedro Simões
- CIEPQPF, Department of Chemical Engineering, University of Coimbra, Rua Sı́lvio de Lima, 3030-790 Coimbra, Portugal
| |
Collapse
|
8
|
Al-Shar'i NA, Al-Balas QA. Molecular Dynamics Simulations of Adenosine Receptors: Advances, Applications and Trends. Curr Pharm Des 2019; 25:783-816. [DOI: 10.2174/1381612825666190304123414] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2019] [Accepted: 02/26/2019] [Indexed: 01/09/2023]
Abstract
:
Adenosine receptors (ARs) are transmembrane proteins that belong to the G protein-coupled receptors
(GPCRs) superfamily and mediate the biological functions of adenosine. To date, four AR subtypes are known,
namely A1, A2A, A2B and A3 that exhibit different signaling pathways, tissue localization, and mechanisms of
activation. Moreover, the widespread ARs and their implication in numerous physiological and pathophysiological
conditions had made them pivotal therapeutic targets for developing clinically effective agents.
:
The crystallographic success in identifying the 3D crystal structures of A2A and A1 ARs has dramatically enriched
our understanding of their structural and functional properties such as ligand binding and signal transduction.
This, in turn, has provided a structural basis for a larger contribution of computational methods, particularly molecular
dynamics (MD) simulations, toward further investigation of their molecular properties and designing
bioactive ligands with therapeutic potential. MD simulation has been proved to be an invaluable tool in investigating
ARs and providing answers to some critical questions. For example, MD has been applied in studying ARs
in terms of ligand-receptor interactions, molecular recognition, allosteric modulations, dimerization, and mechanisms
of activation, collectively aiding in the design of subtype selective ligands.
:
In this review, we focused on the advances and different applications of MD simulations utilized to study the
structural and functional aspects of ARs that can foster the structure-based design of drug candidates. In addition,
relevant literature was briefly discussed which establishes a starting point for future advances in the field of drug
discovery to this pivotal group of drug targets.
Collapse
Affiliation(s)
- Nizar A. Al-Shar'i
- Department of Medicinal Chemistry and Pharmacognosy, Faculty of Pharmacy, Jordan University of Science and Technology, P.O. Box 3030, Irbid 22110, Jordan
| | - Qosay A. Al-Balas
- Department of Medicinal Chemistry and Pharmacognosy, Faculty of Pharmacy, Jordan University of Science and Technology, P.O. Box 3030, Irbid 22110, Jordan
| |
Collapse
|
9
|
Sahoo A, Matysiak S. Computational insights into lipid assisted peptide misfolding and aggregation in neurodegeneration. Phys Chem Chem Phys 2019; 21:22679-22694. [DOI: 10.1039/c9cp02765c] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
An overview of recent advances in computational investigation of peptide–lipid interactions in neurodegeneration – Alzheimer's, Parkinson's and Huntington's disease.
Collapse
Affiliation(s)
- Abhilash Sahoo
- Biophysics Program
- Institute of Physical Science and Technology
- University of Maryland
- College Park
- USA
| | - Silvina Matysiak
- Biophysics Program
- Institute of Physical Science and Technology
- University of Maryland
- College Park
- USA
| |
Collapse
|
10
|
Davtyan A, Voth GA, Andersen HC. Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics. J Chem Phys 2016; 145:224107. [DOI: 10.1063/1.4971430] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Aram Davtyan
- Department of Chemistry, The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Hans C. Andersen
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
11
|
Vögele M, Holm C, Smiatek J. Coarse-grained simulations of polyelectrolyte complexes: MARTINI models for poly(styrene sulfonate) and poly(diallyldimethylammonium). J Chem Phys 2016; 143:243151. [PMID: 26723636 DOI: 10.1063/1.4937805] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present simulations of aqueous polyelectrolyte complexes with new MARTINI models for the charged polymers poly(styrene sulfonate) and poly(diallyldimethylammonium). Our coarse-grained polyelectrolyte models allow us to study large length and long time scales with regard to chemical details and thermodynamic properties. The results are compared to the outcomes of previous atomistic molecular dynamics simulations and verify that electrostatic properties are reproduced by our MARTINI coarse-grained approach with reasonable accuracy. Structural similarity between the atomistic and the coarse-grained results is indicated by a comparison between the pair radial distribution functions and the cumulative number of surrounding particles. Our coarse-grained models are able to quantitatively reproduce previous findings like the correct charge compensation mechanism and a reduced dielectric constant of water. These results can be interpreted as the underlying reason for the stability of polyelectrolyte multilayers and complexes and validate the robustness of the proposed models.
Collapse
Affiliation(s)
- Martin Vögele
- Institute for Computational Physics, University of Stuttgart, Stuttgart, Germany
| | - Christian Holm
- Institute for Computational Physics, University of Stuttgart, Stuttgart, Germany
| | - Jens Smiatek
- Institute for Computational Physics, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
12
|
Li G, Shen H, Zhang D, Li Y, Wang H. Coarse-Grained Modeling of Nucleic Acids Using Anisotropic Gay-Berne and Electric Multipole Potentials. J Chem Theory Comput 2016; 12:676-93. [PMID: 26717419 DOI: 10.1021/acs.jctc.5b00903] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
In this work, we attempt to apply a coarse-grained (CG) model, which is based on anisotropic Gay-Berne and electric multipole (EMP) potentials, to the modeling of nucleic acids. First, a comparison has been made between the CG and atomistic models (AMBER point-charge model) in the modeling of DNA and RNA hairpin structures. The CG results have demonstrated a good quality in maintaining the nucleic acid hairpin structures, in reproducing the dynamics of backbone atoms of nucleic acids, and in describing the hydrogen-bonding interactions between nucleic acid base pairs. Second, the CG and atomistic AMBER models yield comparable results in modeling double-stranded DNA and RNA molecules. It is encouraging that our CG model is capable of reproducing many elastic features of nucleic acid base pairs in terms of the distributions of the interbase pair step parameters (such as shift, slide, tilt, and twist) and the intrabase pair parameters (such as buckle, propeller, shear, and stretch). Finally, The GBEMP model has shown a promising ability to predict the melting temperatures of DNA duplexes with different lengths.
Collapse
Affiliation(s)
- Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, Liaoning Province, People's Republic of China
| | - Hujun Shen
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, Liaoning Province, People's Republic of China
| | - Dinglin Zhang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, Liaoning Province, People's Republic of China
| | - Yan Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, Liaoning Province, People's Republic of China
| | - Honglei Wang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences , Dalian 116023, Liaoning Province, People's Republic of China
| |
Collapse
|
13
|
Davtyan A, Dama JF, Voth GA, Andersen HC. Dynamic force matching: A method for constructing dynamical coarse-grained models with realistic time dependence. J Chem Phys 2016; 142:154104. [PMID: 25903863 DOI: 10.1063/1.4917454] [Citation(s) in RCA: 69] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
Coarse-grained (CG) models of molecular systems, with fewer mechanical degrees of freedom than an all-atom model, are used extensively in chemical physics. It is generally accepted that a coarse-grained model that accurately describes equilibrium structural properties (as a result of having a well constructed CG potential energy function) does not necessarily exhibit appropriate dynamical behavior when simulated using conservative Hamiltonian dynamics for the CG degrees of freedom on the CG potential energy surface. Attempts to develop accurate CG dynamic models usually focus on replacing Hamiltonian motion by stochastic but Markovian dynamics on that surface, such as Langevin or Brownian dynamics. However, depending on the nature of the system and the extent of the coarse-graining, a Markovian dynamics for the CG degrees of freedom may not be appropriate. In this paper, we consider the problem of constructing dynamic CG models within the context of the Multi-Scale Coarse-graining (MS-CG) method of Voth and coworkers. We propose a method of converting a MS-CG model into a dynamic CG model by adding degrees of freedom to it in the form of a small number of fictitious particles that interact with the CG degrees of freedom in simple ways and that are subject to Langevin forces. The dynamic models are members of a class of nonlinear systems interacting with special heat baths that were studied by Zwanzig [J. Stat. Phys. 9, 215 (1973)]. The properties of the fictitious particles can be inferred from analysis of the dynamics of all-atom simulations of the system of interest. This is analogous to the fact that the MS-CG method generates the CG potential from analysis of equilibrium structures observed in all-atom simulation data. The dynamic models generate a non-Markovian dynamics for the CG degrees of freedom, but they can be easily simulated using standard molecular dynamics programs. We present tests of this method on a series of simple examples that demonstrate that the method provides realistic dynamical CG models that have non-Markovian or close to Markovian behavior that is consistent with the actual dynamical behavior of the all-atom system used to construct the CG model. Both the construction and the simulation of such a dynamic CG model have computational requirements that are similar to those of the corresponding MS-CG model and are good candidates for CG modeling of very large systems.
Collapse
Affiliation(s)
- Aram Davtyan
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - James F Dama
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Hans C Andersen
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
14
|
Wang Y, Wang B, Liao H, Song X, Wu H, Wang H, Shen H, Ma X, Tan M. Liposomal nanohybrid cerasomes for mitochondria-targeted drug delivery. J Mater Chem B 2015; 3:7291-7299. [PMID: 32262837 DOI: 10.1039/c5tb01197c] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Mitochondrial dysfunctions cause numerous human disorders and the development of mitochondria-targeted nanocarriers for drug delivery has aroused great attention. Herein, we report the synthesis of a liposomal nanohybrid cerasome modified with triphosphonium (TPP) for drug delivery to the mitochondrial matrix. The cerasomes were observed to possess an average size of about 38 nm in diameter, and the theoretical simulation of GBEMP mapping demonstrated that the amphiphilic organotrialkoxysilanes were stable as a bilayer equilibrium conformation after self-assembly. The cerasomes showed good stability, excellent biocompatibility and sustainable drug release behavior. Moreover, the TPP-targeted cerasomes resulted in greater drug accumulation in mitochondria, thus leading to a greater antitumor effect as compared to non-targeted cerasomes by using doxorubicin as a modal drug. The specific accumulation of TPP-targeted cerasomes within mitochondria was also confirmed by using JC-1 as the fluorescent probe to analyze the mitochondrial transmembrane potential change.
Collapse
Affiliation(s)
- Yanfang Wang
- Division of Biotechnology, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | | | | | | | | | | | | | | | | |
Collapse
|
15
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
16
|
Shen H, Li Y, Xu P, Li X, Chu H, Zhang D, Li G. An anisotropic coarse-grained model based on Gay-Berne and electric multipole potentials and its application to simulate a DMPC bilayer in an implicit solvent model. J Comput Chem 2015; 36:1103-13. [DOI: 10.1002/jcc.23895] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2014] [Accepted: 02/08/2015] [Indexed: 01/12/2023]
Affiliation(s)
- Hujun Shen
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| | - Yan Li
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| | - Peijun Xu
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| | - Xiaofang Li
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| | - Huiying Chu
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| | - Dinglin Zhang
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design; State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences; Dalian 116023, Liaoning Province Peoples Republic of China
| |
Collapse
|
17
|
Davtyan A, Dama JF, Sinitskiy AV, Voth GA. The Theory of Ultra-Coarse-Graining. 2. Numerical Implementation. J Chem Theory Comput 2014; 10:5265-75. [PMID: 26583210 DOI: 10.1021/ct500834t] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The increasing interest in the modeling of complex macromolecular systems in recent years has spurred the development of numerous coarse-graining (CG) techniques. However, many of the CG models are constructed assuming that all details beneath the resolution of CG degrees of freedom are fast and average out, which sets limits on the resolution of feasible coarse-grainings and on the range of applications of the CG models. Ultra-coarse-graining (UCG) makes it possible to construct models at any desired resolution while accounting for discrete conformational or chemical changes within the CG sites that can modulate the interactions between them. Here, we discuss the UCG methodology and its numerical implementation. We pay particular attention to the numerical mechanism for including state transitions between different conformations within CG sites because this has not been discussed previously. Using a simple example of 1,2-dichloroethane, we demonstrate the ability of the UCG model to reproduce the multiconfigurational behavior of this molecular liquid, even when each molecule is modeled with only one CG site. The methodology can also be applied to other molecular liquids and macromolecular systems with time scale separation between conformational transitions and other intramolecular motions and rotations.
Collapse
Affiliation(s)
- Aram Davtyan
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago , Chicago, Illinois 60637, United States
| | - James F Dama
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago , Chicago, Illinois 60637, United States
| | - Anton V Sinitskiy
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago , Chicago, Illinois 60637, United States
| | - Gregory A Voth
- Department of Chemistry, The James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago , Chicago, Illinois 60637, United States
| |
Collapse
|
18
|
Braun D, Boresch S, Steinhauser O. Transport and dielectric properties of water and the influence of coarse-graining: Comparing BMW, SPC/E, and TIP3P models. J Chem Phys 2014; 140:064107. [DOI: 10.1063/1.4864117] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
|
19
|
Shen H, Li Y, Ren P, Zhang D, Li G. An Anisotropic Coarse-Grained Model for Proteins Based On Gay-Berne and Electric Multipole Potentials. J Chem Theory Comput 2014; 10:731-750. [PMID: 24659927 PMCID: PMC3958967 DOI: 10.1021/ct400974z] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
![]()
Gay–Berne
anisotropic potential has been widely used to
evaluate the nonbonded interactions between coarse-grained particles
being described as elliptical rigid bodies. In this paper, we are
presenting a coarse-grained model for twenty kinds of amino acids
and proteins, based on the anisotropic Gay–Berne and point
electric multipole (EMP) potentials. We demonstrate that the anisotropic
coarse-grained model, namely GBEMP model, is able to reproduce many
key features observed from experimental protein structures (Dunbrack
Library), as well as from atomistic force field simulations (using
AMOEBA, AMBER, and CHARMM force fields), while saving the computational
cost by a factor of about 10–200 depending on specific cases
and atomistic models. More importantly, unlike other coarse-grained
approaches, our framework is based on the fundamental intermolecular
forces with explicit treatment of electrostatic and repulsion-dispersion
forces. As a result, the coarse-grained protein model presented an
accurate description of nonbonded interactions (particularly electrostatic
component) between hetero/homodimers (such as peptide–peptide,
peptide–water). In addition, the encouraging performance of
the model was reflected by the excellent correlation between GBEMP
and AMOEBA models in the calculations of the dipole moment of peptides.
In brief, the GBEMP model given here is general and transferable,
suitable for simulating complex biomolecular systems.
Collapse
Affiliation(s)
- Hujun Shen
- Laboratory of Molecular Modeling and Design, State key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Rd. Dalian 116023, PR China
| | - Yan Li
- Laboratory of Molecular Modeling and Design, State key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Rd. Dalian 116023, PR China
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA
| | - Dinglin Zhang
- Laboratory of Molecular Modeling and Design, State key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Rd. Dalian 116023, PR China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, 457 Zhongshan Rd. Dalian 116023, PR China
| |
Collapse
|
20
|
Page AJ, Isomoto T, Knaup JM, Irle S, Morokuma K. Effects of Molecular Dynamics Thermostats on Descriptions of Chemical Nonequilibrium. J Chem Theory Comput 2012; 8:4019-28. [PMID: 26605569 DOI: 10.1021/ct3004639] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The performance of popular molecular dynamics (MD) thermostat algorithms in constant temperature simulations of equilibrium systems is well-known. This is not the case, however, in the context of nonequilibrium chemical systems, such as chemical reactions or nanoscale self-assembly processes. In this work, we investigate the effect of popular thermostat algorithms on the "natural" (i.e., Hamiltonian) dynamics of a nonequilibrium, chemically reacting system. By comparing constant-temperature quantum mechanical MD (QM/MD) simulations of carbon vapor condensation using velocity scaling, Berendsen, Andersen, Langevin, and Nosé-Hoover chain thermostat algorithms with natural NVE simulations, we show that efficient temperature control and reliable reaction dynamics are mutually exclusive in such a system. This problem may be circumvented, however, by placing the reactive system in an inert He atmosphere, which is itself described using NVT MD. We demonstrate that both realistic temperature control and dynamics consistent with natural NVE dynamics can then be obtained simultaneously. In essence, the thermal energy created by the natural dynamics of the NVE subsystem is drained by the thermostat acting on the NVT atmosphere, without adversely affecting the dynamics of the reactive system itself.
Collapse
Affiliation(s)
- Alister J Page
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Tetsushi Isomoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| | - Jan M Knaup
- Department of Physics, Harvard University, Cambridge, Masachusetts 02138, United States.,Bremen Center for Computational Materials Science - BCCMS, University of Bremen , 28359 Bremen, Germany
| | - Stephan Irle
- Department of Chemistry, Graduate School of Science, Nagoya University, Nagoya 464-4602, Japan
| | - Keiji Morokuma
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan.,Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
21
|
Shen H, Xia Z, Li G, Ren P. A Review of Physics-Based Coarse-Grained Potentials for the Simulations of Protein Structure and Dynamics. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY VOLUME 8 2012. [DOI: 10.1016/b978-0-444-59440-2.00005-3] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|