1
|
van Gunsteren WF, Oostenbrink C. Methods for Classical-Mechanical Molecular Simulation in Chemistry: Achievements, Limitations, Perspectives. J Chem Inf Model 2024; 64:6281-6304. [PMID: 39136351 DOI: 10.1021/acs.jcim.4c00823] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/27/2024]
Abstract
More than a half century ago it became feasible to simulate, using classical-mechanical equations of motion, the dynamics of molecular systems on a computer. Since then classical-physical molecular simulation has become an integral part of chemical research. It is widely applied in a variety of branches of chemistry and has significantly contributed to the development of chemical knowledge. It offers understanding and interpretation of experimental results, semiquantitative predictions for measurable and nonmeasurable properties of substances, and allows the calculation of properties of molecular systems under conditions that are experimentally inaccessible. Yet, molecular simulation is built on a number of assumptions, approximations, and simplifications which limit its range of applicability and its accuracy. These concern the potential-energy function used, adequate sampling of the vast statistical-mechanical configurational space of a molecular system and the methods used to compute particular properties of chemical systems from statistical-mechanical ensembles. During the past half century various methodological ideas to improve the efficiency and accuracy of classical-physical molecular simulation have been proposed, investigated, evaluated, implemented in general simulation software or were abandoned. The latter because of fundamental flaws or, while being physically sound, computational inefficiency. Some of these methodological ideas are briefly reviewed and the most effective methods are highlighted. Limitations of classical-physical simulation are discussed and perspectives are sketched.
Collapse
Affiliation(s)
- Wilfred F van Gunsteren
- Institute for Molecular Physical Science, Swiss Federal Institute of Technology, ETH, CH-8093 Zurich, Switzerland
| | - Chris Oostenbrink
- Institute of Molecular Modelling and Simulation, BOKU University, 1190 Vienna, Austria
- Christian Doppler Laboratory for Molecular Informatics in the Biosciences, BOKU University, Muthgasse 18, 1190 Vienna, Austria
| |
Collapse
|
2
|
Capone M, Romanelli M, Castaldo D, Parolin G, Bello A, Gil G, Vanzan M. A Vision for the Future of Multiscale Modeling. ACS PHYSICAL CHEMISTRY AU 2024; 4:202-225. [PMID: 38800726 PMCID: PMC11117712 DOI: 10.1021/acsphyschemau.3c00080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/30/2023] [Revised: 01/31/2024] [Accepted: 02/01/2024] [Indexed: 05/29/2024]
Abstract
The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.
Collapse
Affiliation(s)
- Matteo Capone
- Department
of Physical and Chemical Sciences, University
of L’Aquila, L’Aquila 67010, Italy
| | - Marco Romanelli
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Davide Castaldo
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Giovanni Parolin
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
| | - Alessandro Bello
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
- Department
of Physics, Informatics and Mathematics, University of Modena and Reggio Emilia, Modena 41125, Italy
| | - Gabriel Gil
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
- Instituto
de Cibernética, Matemática y Física (ICIMAF), La Habana 10400, Cuba
| | - Mirko Vanzan
- Department
of Chemical Sciences, University of Padova, Padova 35131, Italy
- Department
of Physics, University of Milano, Milano 20133, Italy
| |
Collapse
|
3
|
Yue S, Zhao J, Sun Y, Niu H, Li H, Jing Y, Aluru NR. Multi-scale simulation of proton diffusion in dislocation cores in BaZrO 3. Phys Chem Chem Phys 2022; 24:21440-21451. [PMID: 36047850 DOI: 10.1039/d2cp02989h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Dislocations are important for their effects on the chemical, electrical, magnetic, and transport properties of oxide materials, especially for electrochemical devices such as solid fuel cells and resistive memories, but these effects are still under-studied at the atomic level. We have developed a quantum mechanical/molecular mechanical (QM/MM)-based multiscale simulation program to reveal the diffusion properties of protons on 〈100〉 edge dislocations in BaZrO3 perovskite oxide. We find that the large free space and the presence of hydrogen bonds in the dislocation core structure lead to significant trapping of protons. The diffusion properties of protons in dislocation cores were investigated, and no evidence of pipeline diffusion was found from the calculated migration energy barriers, which not only did not accelerate ion diffusion but rather decreases the conductivity of ions. The proton diffusion properties of Y-doped BaZrO3 (BZY), with a dislocation core structure (BZY-D) and with a grain boundary structure (BZY-GB) were also compared. In all three structures, local lattice deformation occupies an essential part in the proton transfer and rotation processes. The change in bond order is calculated and it is found that the interaction with oxygen and Zr ions during proton transfer and rotation controls the energy barrier for local lattice deformation of the O-B-O motion, which affects the proton diffusion in the structure. Our study provides insight into proton diffusion in dislocations in terms of mechanical behavior, elucidates the origin of the energy barrier associated with proton diffusion in dislocations, and provides guidance for the preparation and application of proton conductors.
Collapse
Affiliation(s)
- Shaofeng Yue
- Department of Astronautical Science and Mechanics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, P. R. China.
| | - Junqing Zhao
- Department of Astronautical Science and Mechanics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, P. R. China.
| | - Yi Sun
- Department of Astronautical Science and Mechanics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, P. R. China.
| | - Hongwei Niu
- Department of Astronautical Science and Mechanics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, P. R. China.
| | - Huyang Li
- Department of Astronautical Science and Mechanics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, P. R. China.
| | - Yuhang Jing
- Department of Astronautical Science and Mechanics, Harbin Institute of Technology, Harbin, Heilongjiang 150001, P. R. China.
| | - N R Aluru
- Walker Department of Mechanical Engineering, Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA
| |
Collapse
|
4
|
Kekenes-Huskey PM, Burgess DE, Sun B, Bartos DC, Rozmus ER, Anderson CL, January CT, Eckhardt LL, Delisle BP. Mutation-Specific Differences in Kv7.1 ( KCNQ1) and Kv11.1 ( KCNH2) Channel Dysfunction and Long QT Syndrome Phenotypes. Int J Mol Sci 2022; 23:7389. [PMID: 35806392 PMCID: PMC9266926 DOI: 10.3390/ijms23137389] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2022] [Revised: 06/22/2022] [Accepted: 06/24/2022] [Indexed: 11/16/2022] Open
Abstract
The electrocardiogram (ECG) empowered clinician scientists to measure the electrical activity of the heart noninvasively to identify arrhythmias and heart disease. Shortly after the standardization of the 12-lead ECG for the diagnosis of heart disease, several families with autosomal recessive (Jervell and Lange-Nielsen Syndrome) and dominant (Romano-Ward Syndrome) forms of long QT syndrome (LQTS) were identified. An abnormally long heart rate-corrected QT-interval was established as a biomarker for the risk of sudden cardiac death. Since then, the International LQTS Registry was established; a phenotypic scoring system to identify LQTS patients was developed; the major genes that associate with typical forms of LQTS were identified; and guidelines for the successful management of patients advanced. In this review, we discuss the molecular and cellular mechanisms for LQTS associated with missense variants in KCNQ1 (LQT1) and KCNH2 (LQT2). We move beyond the "benign" to a "pathogenic" binary classification scheme for different KCNQ1 and KCNH2 missense variants and discuss gene- and mutation-specific differences in K+ channel dysfunction, which can predispose people to distinct clinical phenotypes (e.g., concealed, pleiotropic, severe, etc.). We conclude by discussing the emerging computational structural modeling strategies that will distinguish between dysfunctional subtypes of KCNQ1 and KCNH2 variants, with the goal of realizing a layered precision medicine approach focused on individuals.
Collapse
Affiliation(s)
- Peter M. Kekenes-Huskey
- Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, IL 60153, USA
| | - Don E. Burgess
- Department of Physiology, College of Medicine, University of Kentucky, Lexington, KY 40536, USA; (D.E.B.); (E.R.R.)
| | - Bin Sun
- Department of Pharmacology, Harbin Medical University, Harbin 150081, China;
| | | | - Ezekiel R. Rozmus
- Department of Physiology, College of Medicine, University of Kentucky, Lexington, KY 40536, USA; (D.E.B.); (E.R.R.)
| | - Corey L. Anderson
- Cellular and Molecular Arrythmias Program, Division of Cardiovascular Medicine, Department of Medicine, University of Wisconsin-Madison, Madison, WI 53705, USA; (C.L.A.); (C.T.J.); (L.L.E.)
| | - Craig T. January
- Cellular and Molecular Arrythmias Program, Division of Cardiovascular Medicine, Department of Medicine, University of Wisconsin-Madison, Madison, WI 53705, USA; (C.L.A.); (C.T.J.); (L.L.E.)
| | - Lee L. Eckhardt
- Cellular and Molecular Arrythmias Program, Division of Cardiovascular Medicine, Department of Medicine, University of Wisconsin-Madison, Madison, WI 53705, USA; (C.L.A.); (C.T.J.); (L.L.E.)
| | - Brian P. Delisle
- Department of Physiology, College of Medicine, University of Kentucky, Lexington, KY 40536, USA; (D.E.B.); (E.R.R.)
| |
Collapse
|
5
|
Wang H, Zheng X. Theoretical Study of Macrocyclic Host Molecules: From Supramolecular Recognition to Self-Assembly. Phys Chem Chem Phys 2022; 24:19011-19028. [DOI: 10.1039/d2cp02152h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Supramolecular chemistry focuses on molecular recognition and self-assembly of various building blocks through weak non-covalent interactions, including anion-π, hydrogen bond (HB), hydrophobic interactions, van der Waals (vdW) interactions, etc, which...
Collapse
|
6
|
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
|
7
|
Chan C, Du S, Dong Y, Cheng X. Computational and Experimental Approaches to Investigate Lipid Nanoparticles as Drug and Gene Delivery Systems. Curr Top Med Chem 2021; 21:92-114. [PMID: 33243123 PMCID: PMC8191596 DOI: 10.2174/1568026620666201126162945] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Revised: 10/16/2020] [Accepted: 10/22/2020] [Indexed: 02/06/2023]
Abstract
Lipid nanoparticles (LNPs) have been widely applied in drug and gene delivery. More than twenty years ago, DoxilTM was the first LNPs-based drug approved by the US Food and Drug Administration (FDA). Since then, with decades of research and development, more and more LNP-based therapeutics have been used to treat diverse diseases, which often offer the benefits of reduced toxicity and/or enhanced efficacy compared to the active ingredients alone. Here, we provide a review of recent advances in the development of efficient and robust LNPs for drug/gene delivery. We emphasize the importance of rationally combining experimental and computational approaches, especially those providing multiscale structural and functional information of LNPs, to the design of novel and powerful LNP-based delivery systems.
Collapse
Affiliation(s)
- Chun Chan
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy, The Ohio State University, Columbus, OH 43210, USA
| | - Shi Du
- Division of Pharmaceutics and Pharmacology, College of Pharmacy, The Ohio State University, Columbus, OH 43210, USA
| | - Yizhou Dong
- Division of Pharmaceutics and Pharmacology, College of Pharmacy, The Ohio State University, Columbus, OH 43210, USA
- Department of Biomedical Engineering; The Center for Clinical and Translational Science; The Comprehensive Cancer Center; Dorothy M. Davis Heart & Lung Research Institute; Department of Radiation Oncology, The Ohio State University, Columbus, OH 43210, USA
| | - Xiaolin Cheng
- Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy, The Ohio State University, Columbus, OH 43210, USA
- Biophysics Graduate Program, Translational Data Analytics Institute, The Ohio State University, Columbus, OH 43210, USA
| |
Collapse
|
8
|
van Gunsteren WF, Daura X, Fuchs PFJ, Hansen N, Horta BAC, Hünenberger PH, Mark AE, Pechlaner M, Riniker S, Oostenbrink C. On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues. Chemphyschem 2020; 22:264-282. [PMID: 33377305 DOI: 10.1002/cphc.202000968] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Indexed: 12/14/2022]
Abstract
Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.
Collapse
Affiliation(s)
- Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Xavier Daura
- Institute of Biotechnology and Biomedicine, Universitat Autonoma de Barcelona (UAB), 08193, Barcelona, Spain.,Catalan Institution for Research and Advanced Studies (ICREA), 08010, Barcelona, Spain
| | - Patrick F J Fuchs
- Sorbonne Université, Ecole Normale Supérieure, PSL Research University, CNRS, Laboratoire des Biomolécules (LBM), F-75005, Paris, France.,Université de Paris, UFR Sciences du Vivant, F-75013, Paris, France
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569, Stuttgart, Germany
| | - Bruno A C Horta
- Instituto de Química, Universidade Federal de Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| | - Philippe H Hünenberger
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Alan E Mark
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, QLD, 4072, Australia
| | - Maria Pechlaner
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Chris Oostenbrink
- Institute of Molecular Modelling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| |
Collapse
|
9
|
Bunker A, Róg T. Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery. Front Mol Biosci 2020; 7:604770. [PMID: 33330633 PMCID: PMC7732618 DOI: 10.3389/fmolb.2020.604770] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 11/02/2020] [Indexed: 12/12/2022] Open
Abstract
In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.
Collapse
Affiliation(s)
- Alex Bunker
- Division of Pharmaceutical Biosciences, Drug Research Program, Faculty of Pharmacy, University of Helsinki, Helsinki, Finland
| | - Tomasz Róg
- Department of Physics, University of Helsinki, Helsinki, Finland
| |
Collapse
|
10
|
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]
|
11
|
Bastida A, Zúñiga J, Requena A, Cerezo J. Energetic Self-Folding Mechanism in α-Helices. J Phys Chem B 2019; 123:8186-8194. [PMID: 31498638 DOI: 10.1021/acs.jpcb.9b05860] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
A novel energetic route driving the folding of a polyalanine peptide from an extended conformation to its α-helix native conformation is described, supported by a new method to compute mean potential energy surfaces accurately in terms of the dihedral angles of the peptide chain from extensive molecular dynamics simulations. The energetic self-folding (ESF) route arises specifically from the balance between the intrinsic propensity of alanine residues toward the αR conformation and two, opposite, effects: the destabilizing interaction with neighbor residues and the stabilizing formation of native hydrogen bonds, with the latter being dominant for large peptide lengths. The ESF mechanism provides simple but robust support to the nucleation-elongation or zipper models and offers a quantitative energetic funnel picture of the folding process. The mechanism is validated by the reasonable agreement between the computed folding energies and the experimental values.
Collapse
Affiliation(s)
- Adolfo Bastida
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| | - José Zúñiga
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| | - Alberto Requena
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| | - Javier Cerezo
- Departamento de Química Física , Universidad de Murcia , 30100 Murcia , Spain
| |
Collapse
|
12
|
The Roots of Bio‐Molecular Simulation: The Eight‐Week CECAM Workshop ‘Models for Protein Dynamics’ of 1976. Helv Chim Acta 2019. [DOI: 10.1002/hlca.201800239] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
13
|
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
|
14
|
Hahn DF, Milić JV, Hünenberger PH. Vase
‐
Kite
Equilibrium of Resorcin[4]arene Cavitands Investigated Using Molecular Dynamics Simulations with Ball‐and‐Stick Local Elevation Umbrella Sampling. Helv Chim Acta 2019. [DOI: 10.1002/hlca.201900060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- David F. Hahn
- Laboratory of Physical Chemistry, Department of Chemistry and Applied BiosciencesETH Zürich Vladimir-Prelog-Weg 2 CH-8093 Zürich Switzerland
| | - Jovana V. Milić
- Laboratory of Photonics and InterfacesÉcole Polytechnique Fédérale de Lausanne, EPFL SB ISIC LPI, Station 6 CH-1015 Lausanne Switzerland
| | - Philippe H. Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied BiosciencesETH Zürich Vladimir-Prelog-Weg 2 CH-8093 Zürich Switzerland
| |
Collapse
|
15
|
Wei W, Ma J, Xie D, Zhou Y. Linking inhibitor motions to proteolytic stability of sunflower trypsin inhibitor-1. RSC Adv 2019; 9:13776-13786. [PMID: 35519558 PMCID: PMC9063939 DOI: 10.1039/c9ra02114k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2019] [Accepted: 04/25/2019] [Indexed: 12/20/2022] Open
Abstract
The remarkable capability of an enzyme isn't only determined by its active site but also controlled by the environment. To unravel the environment role in catalysis, the dynamic motions as well as the static mechanism need to be studied. In this work, QM/MM MD simulations were employed to study the proteolysis process of SFTI-1 and BiKF, which revealed that a combination of static non-bonded interactions and dynamic motions along the reaction coordinate can account for the different hydrolysis rates between them. A comparison among SFTI-1 and three analogs with similar non-bonded interactions further revealed a positive correlation between the mobility of inhibitors and the hydrolysis rates. Apart from the cyclic backbone and disulfide bond, intramolecular hydrogen bonds also increase the rigidity of the backbone of inhibitors, and therefore hinder inhibitor motions to resist proteolysis. These new detailed mechanistic insights suggest the need to consider inhibitor motions in the rational design of peptide inhibitors. Besides the non-bonded interactions, inhibitor motions especially rotation of the scissile bond also influence proteolytic stability.![]()
Collapse
Affiliation(s)
- Wanqing Wei
- Institute of Theoretical and Computational Chemistry, Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Jing Ma
- Institute of Theoretical and Computational Chemistry, Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| | - Yanzi Zhou
- Institute of Theoretical and Computational Chemistry, Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University Nanjing 210023 China
| |
Collapse
|
16
|
Renevey A, Riniker S. Benchmarking Hybrid Atomistic/Coarse-Grained Schemes for Proteins with an Atomistic Water Layer. J Phys Chem B 2019; 123:3033-3042. [DOI: 10.1021/acs.jpcb.8b12149] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Annick Renevey
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
17
|
Wei W, Chen Y, Xie D, Zhou Y. Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation. Phys Chem Chem Phys 2019; 21:5049-5058. [PMID: 30762035 DOI: 10.1039/c8cp07784c] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.
Collapse
Affiliation(s)
- Wanqing Wei
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China.
| | | | | | | |
Collapse
|
18
|
Duster AW, Garza CM, Aydintug BO, Negussie MB, Lin H. Adaptive Partitioning QM/MM for Molecular Dynamics Simulations: 6. Proton Transport through a Biological Channel. J Chem Theory Comput 2019; 15:892-905. [DOI: 10.1021/acs.jctc.8b01128] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Adam W. Duster
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Christina M. Garza
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Baris O. Aydintug
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Mikias B. Negussie
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| | - Hai Lin
- Chemistry Department, CB 194, University of Colorado, Denver, Colorado 80217, United States
| |
Collapse
|
19
|
Freindorf M, Tao Y, Sethio D, Cremer D, Kraka E. New mechanistic insights into the Claisen rearrangement of chorismate – a Unified Reaction Valley Approach study. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1530464] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Daniel Sethio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Dieter Cremer
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| |
Collapse
|
20
|
Duster AW, Wang CH, Lin H. Adaptive QM/MM for Molecular Dynamics Simulations: 5. On the Energy-Conserved Permuted Adaptive-Partitioning Schemes. Molecules 2018; 23:E2170. [PMID: 30154373 PMCID: PMC6225285 DOI: 10.3390/molecules23092170] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Revised: 08/23/2018] [Accepted: 08/24/2018] [Indexed: 11/16/2022] Open
Abstract
In combined quantum-mechanical/molecular-mechanical (QM/MM) dynamics simulations, the adaptive-partitioning (AP) schemes reclassify atoms on-the-fly as QM or MM in a smooth manner. This yields a mobile QM subsystem with contents that are continuously updated as needed. Here, we tailor the Hamiltonian adaptive many-body correction (HAMBC) proposed by Boreboom et al. [J. Chem. Theory Comput.2016, 12, 3441] to the permuted AP (PAP) scheme. The treatments lead to the HAMBC-PAP method (HPAP), which both conserves energy and produces accurate solvation structures in the test of "water-in-water" model system.
Collapse
Affiliation(s)
- Adam W Duster
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| | - Chun-Hung Wang
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| | - Hai Lin
- Department of Chemistry, University of Colorado Denver, Denver, CO 80217, USA.
| |
Collapse
|
21
|
Sinitskiy AV, Voth GA. Quantum mechanics/coarse-grained molecular mechanics (QM/CG-MM). J Chem Phys 2018; 148:014102. [DOI: 10.1063/1.5006810] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anton V. Sinitskiy
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A. Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
22
|
Ainsley J, Lodola A, Mulholland AJ, Christov CZ, Karabencheva-Christova TG. Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2018; 113:1-32. [PMID: 30149903 DOI: 10.1016/bs.apcsb.2018.07.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.
Collapse
Affiliation(s)
- Jon Ainsley
- Department of Applied Sciences, Faculty of Health and Life Sciences, Northumbria University, Newcastle upon Tyne, United Kingdom
| | | | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | - Christo Z Christov
- Department of Chemistry, Michigan Technological University, Houghton, MI, United States.
| | | |
Collapse
|
23
|
Ding HM, Ma YQ. Computational approaches to cell-nanomaterial interactions: keeping balance between therapeutic efficiency and cytotoxicity. NANOSCALE HORIZONS 2018; 3:6-27. [PMID: 32254106 DOI: 10.1039/c7nh00138j] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Owing to their unique properties, nanomaterials have been widely used in biomedicine since they have obvious inherent advantages over traditional ones. However, nanomaterials may also cause dysfunction in proteins, genes and cells, resulting in cytotoxic and genotoxic responses. Recently, more and more attention has been paid to these potential toxicities of nanomaterials, especially to the risks of nanomaterials to human health and safety. Therefore, when using nanomaterials for biomedical applications, it is of great importance to keep the balance between therapeutic efficiency and cytotoxicity (i.e., increase the therapeutic efficiency as well as decrease the potential toxicity). This requires a deeper understanding of the interactions between various types of nanomaterials and biological systems at the nano/bio interface. In this review, from the point of view of theoretical researchers, we will present the current status regarding the physical mechanism of cytotoxicity caused by nanomaterials, mainly based on recent simulation results. In addition, the strategies for minimizing the nanotoxicity naturally and artificially will also be discussed in detail. Furthermore, we should notice that toxicity is not always bad for clinical use since causing the death of specific cells is the main way of treating disease. Enhancing the targeting ability of nanomaterials to diseased cells and minimizing their side effects on normal cells will always be hugely challenging issues in nanomedicine. By combining the latest computational studies with some experimental verifications, we will provide special insights into recent advances regarding these problems, especially for the design of novel environment-responsive nanomaterials.
Collapse
Affiliation(s)
- Hong-Ming Ding
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, Soochow University, Suzhou 215006, China
| | | |
Collapse
|
24
|
van Gunsteren WF, Daura X, Hansen N, Mark AE, Oostenbrink C, Riniker S, Smith LJ. Validation of Molecular Simulation: An Overview of Issues. Angew Chem Int Ed Engl 2017; 57:884-902. [PMID: 28682472 DOI: 10.1002/anie.201702945] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2017] [Indexed: 12/14/2022]
Abstract
Computer simulation of molecular systems enables structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic, or supra-molecular level. To interpret results of such simulations appropriately, the quality of the calculated properties must be evaluated. This depends on the way the simulations are performed and on the way they are validated by comparison to values Qexp of experimentally observable quantities Q. One must consider 1) the accuracy of Qexp , 2) the accuracy of the function Q(rN ) used to calculate a Q-value based on a molecular configuration rN of N particles, 3) the sensitivity of the function Q(rN ) to the configuration rN , 4) the relative time scales of the simulation and experiment, 5) the degree to which the calculated and experimental properties are equivalent, and 6) the degree to which the system simulated matches the experimental conditions. Experimental data is limited in scope and generally corresponds to averages over both time and space. A critical analysis of the various factors influencing the apparent degree of (dis)agreement between simulations and experiment is presented and illustrated using examples from the literature. What can be done to enhance the validation of molecular simulation is also discussed.
Collapse
Affiliation(s)
- Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Xavier Daura
- Institute of Biotechnology and Biomedicine, Universitat Autonoma de Barcelona, UAB, 08193, Barcelona, Spain.,Catalan Institution for Research and Advanced Studies, ICREA, 08010, Barcelona, Spain
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569, Stuttgart, Germany
| | - Alan E Mark
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, QLD 4072, Australia
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Sereina Riniker
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH, 8093, Zurich, Switzerland
| | - Lorna J Smith
- Department of Chemistry, Inorganic Chemistry, Laboratory, University of Oxford, South Parks Road, Oxford, OX1 3QR, UK
| |
Collapse
|
25
|
van Gunsteren WF, Daura X, Hansen N, Mark AE, Oostenbrink C, Riniker S, Smith LJ. Validierung von molekularen Simulationen: eine Übersicht verschiedener Aspekte. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201702945] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Wilfred F. van Gunsteren
- Laboratorium für Physikalische Chemie; Eidgenössische Technische Hochschule Zürich; 8093 Zürich Schweiz
| | - Xavier Daura
- Institute of Biotechnology and Biomedicine; Universitat Autònoma de Barcelona; Spanien
- Catalan Institution for Research and Advanced Studies, ICREA; Barcelona Spanien
| | - Niels Hansen
- Institut für Technische Thermodynamik und Thermische Verfahrenstechnik; Universität Stuttgart; Deutschland
| | - Alan E. Mark
- School of Chemistry and Molecular Biosciences; University of Queensland; St. Lucia Australien
| | - Chris Oostenbrink
- Institut für Molekulare Modellierung und Simulation; Universität für Bodenkultur Wien; Österreich
| | - Sereina Riniker
- Laboratorium für Physikalische Chemie; Eidgenössische Technische Hochschule Zürich; 8093 Zürich Schweiz
| | - Lorna J. Smith
- Inorganic Chemistry Laboratory; Department of Chemistry; University of Oxford; Großbritannien
| |
Collapse
|
26
|
Grimme S, Bannwarth C, Caldeweyher E, Pisarek J, Hansen A. A general intermolecular force field based on tight-binding quantum chemical calculations. J Chem Phys 2017; 147:161708. [DOI: 10.1063/1.4991798] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Jana Pisarek
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie der Universität Bonn, Beringstr. 4, D-53115 Bonn,
Germany
| |
Collapse
|
27
|
Feliu N, Hassan M, Garcia Rico E, Cui D, Parak W, Alvarez-Puebla R. SERS Quantification and Characterization of Proteins and Other Biomolecules. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:9711-9730. [PMID: 28826207 DOI: 10.1021/acs.langmuir.7b01567] [Citation(s) in RCA: 89] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Changes in protein expression levels and protein structure may indicate genomic mutations and may be related to some diseases. Therefore, the precise quantification and characterization of proteins can be used for disease diagnosis. Compared with several other alternative methods, surface-enhanced Raman scattering (SERS) spectroscopy is regarded as an excellent choice for the quantification and structural characterization of proteins. Herein, we review the main advance of using plasmonic nanostructures as SERS sensing platform for this purpose. Three design approaches, including direct SERS, indirect SERS, and SERS-encoded nanoparticles, are discussed in the direction of developing new precise approaches of quantification and characterization of proteins. While this Review is focused on proteins, in order to highlight concepts of SERS-based sensors also detection of other biomolecules will be discussed.
Collapse
Affiliation(s)
- Neus Feliu
- Fachbereich Physik, Philipps Universität Marburg , 35037 Marburg, Germany
- Experimental Cancer Medicine, Department of Laboratory Medicine, Karolinska Institutet , Stockholm, 141 86 Sweden
| | - Moustapha Hassan
- Experimental Cancer Medicine, Department of Laboratory Medicine, Karolinska Institutet , Stockholm, 141 86 Sweden
| | - Eduardo Garcia Rico
- Fundacion de Investigacion HM Hospitales , San Bernardo 101, 28015 Madrid, Spain
- Centro Integral Oncologico Clara Campal (CIOCC) , Oña 10, 28050 Madrid, Spain
- Servicio de Oncologia Clinica, Hospital Universitario HM Torrelodones , Castillo de Olivares s/n, 28250 Torrelodones, Spain
- School of Medicine, San Pablo CEU , Calle Julián Romea, 18, 28003 Madrid, Spain
| | - Daxiang Cui
- Institute of Nano Biomedicine and Engineering, Key Laboratory for Thin Film and Microfabrication Technology of the Ministry of Education, Department of Instrument Science and Engineering, School of Electronic Information and Electrical Engineering, National Center for Translational Medicine, Shanghai Jiao Tong University , 200240 Shanghai, China
| | - Wolfgang Parak
- Fachbereich Physik, Philipps Universität Marburg , 35037 Marburg, Germany
- Institute of Nano Biomedicine and Engineering, Key Laboratory for Thin Film and Microfabrication Technology of the Ministry of Education, Department of Instrument Science and Engineering, School of Electronic Information and Electrical Engineering, National Center for Translational Medicine, Shanghai Jiao Tong University , 200240 Shanghai, China
- Fachbereich Physik und Chemie, Universität Hamburg , 20146 Harmburg, Germany
| | - Ramon Alvarez-Puebla
- Departamento de Química Física e Inorgánica, Universitat Rovira i Virgili , Carrer de Marcellí Domingo s/n, 43007 Tarragona, Spain
- ICREA , Passeig Lluís Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
28
|
Frush EH, Sekharan S, Keinan S. In Silico Prediction of Ligand Binding Energies in Multiple Therapeutic Targets and Diverse Ligand Sets—A Case Study on BACE1, TYK2, HSP90, and PERK Proteins. J Phys Chem B 2017; 121:8142-8148. [DOI: 10.1021/acs.jpcb.7b07224] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Elizabeth Hatcher Frush
- Cloud Pharmaceuticals, Inc., 6 Davis Drive,
Research Triangle Park, North Carolina 27709, United States
| | - Sivakumar Sekharan
- Cloud Pharmaceuticals, Inc., 6 Davis Drive,
Research Triangle Park, North Carolina 27709, United States
| | - Shahar Keinan
- Cloud Pharmaceuticals, Inc., 6 Davis Drive,
Research Triangle Park, North Carolina 27709, United States
| |
Collapse
|
29
|
Culka M, Gisdon FJ, Ullmann GM. Computational Biochemistry-Enzyme Mechanisms Explored. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2017; 109:77-112. [PMID: 28683923 DOI: 10.1016/bs.apcsb.2017.04.004] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Understanding enzyme mechanisms is a major task to achieve in order to comprehend how living cells work. Recent advances in biomolecular research provide huge amount of data on enzyme kinetics and structure. The analysis of diverse experimental results and their combination into an overall picture is, however, often challenging. Microscopic details of the enzymatic processes are often anticipated based on several hints from macroscopic experimental data. Computational biochemistry aims at creation of a computational model of an enzyme in order to explain microscopic details of the catalytic process and reproduce or predict macroscopic experimental findings. Results of such computations are in part complementary to experimental data and provide an explanation of a biochemical process at the microscopic level. In order to evaluate the mechanism of an enzyme, a structural model is constructed which can be analyzed by several theoretical approaches. Several simulation methods can and should be combined to get a reliable picture of the process of interest. Furthermore, abstract models of biological systems can be constructed combining computational and experimental data. In this review, we discuss structural computational models of enzymatic systems. We first discuss various models to simulate enzyme catalysis. Furthermore, we review various approaches how to characterize the enzyme mechanism both qualitatively and quantitatively using different modeling approaches.
Collapse
Affiliation(s)
- Martin Culka
- Computational Biochemistry, University of Bayreuth, Bayreuth, Germany
| | - Florian J Gisdon
- Computational Biochemistry, University of Bayreuth, Bayreuth, Germany
| | | |
Collapse
|
30
|
Duster AW, Wang C, Garza CM, Miller DE, Lin H. Adaptive quantum/molecular mechanics: what have we learned, where are we, and where do we go from here? WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1310] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Adam W. Duster
- Chemistry Department University of Colorado Denver Denver CO USA
| | - Chun‐Hung Wang
- Chemistry Department University of Colorado Denver Denver CO USA
| | | | | | - Hai Lin
- Chemistry Department University of Colorado Denver Denver CO USA
| |
Collapse
|
31
|
Taylor WR. SimGen: A General Simulation Method for Large Systems. J Mol Biol 2017; 429:408-415. [PMID: 27771481 PMCID: PMC5282398 DOI: 10.1016/j.jmb.2016.10.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2016] [Revised: 09/28/2016] [Accepted: 10/06/2016] [Indexed: 11/23/2022]
Abstract
SimGen is a stand-alone computer program that reads a script of commands to represent complex macromolecules, including proteins and nucleic acids, in a structural hierarchy that can then be viewed using an integral graphical viewer or animated through a high-level application programming interface in C++. Structural levels in the hierarchy range from α-carbon or phosphate backbones through secondary structure to domains, molecules, and multimers with each level represented in an identical data structure that can be manipulated using the application programming interface. Unlike most coarse-grained simulation approaches, the higher-level objects represented in SimGen can be soft, allowing the lower-level objects that they contain to interact directly. The default motion simulated by SimGen is a Brownian-like diffusion that can be set to occur across all levels of representation in the hierarchy. Links can also be defined between objects, which, when combined with large high-level random movements, result in an effective search strategy for constraint satisfaction, including structure prediction from predicted pairwise distances. The implementation of SimGen makes use of the hierarchic data structure to avoid unnecessary calculation, especially for collision detection, allowing it to be simultaneously run and viewed on a laptop computer while simulating large systems of over 20,000 objects. It has been used previously to model complex molecular interactions including the motion of a myosin-V dimer "walking" on an actin fibre, RNA stem-loop packing, and the simulation of cell motion and aggregation. Several extensions to this original functionality are described.
Collapse
|
32
|
Zhang X, Zhao Y, Duan X, Zhang HN, Cao Z, Mo Y. Mechanisms for the deamination reaction of 8-oxoguanine catalyzed by 8-oxoguanine deaminase: A combined QM/MM molecular dynamics study. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1142/s0219633616500668] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The deamination reaction of 8-oxoguanine (8-oxoG) catalyzed by 8-oxoguanine deaminase (8-oxoGD) plays a critically important role in the DNA repair activity for oxidative damage. In order to elucidate the complete enzymatic catalysis mechanism at the stages of 8-oxoguanine binding, departure of 2-hydroxy-1H-purine-6,8(7H,9H)-dione from the active site, and formation of 8-oxoxanthine, extensive combined QM(PM3)/MM molecular dynamics simulations have been performed. Computations show that the rate-limiting step corresponds to the nucleophilic attack from zinc-coordinate hydroxide group to free 8-oxoguanine. Through conformational analyses, we demonstrate that Trp115, Trp123 and Leu119 connect to O8@8-oxoguanine with hydrogen bonds, and we suggest that mutations of tryptophan (115 and 123) to histidine or phenylalanine and mutation of leucine (119) to alanine could potentially lead to a mutant with enhanced activity. On this ground, a proton transfer mechanism for the formation of 8-oxoxanthine was further discussed. Both Glu218 and water molecule could be used as proton shuttles, and water molecule plays a major role in proton transfer in substrate. On the other hand, comparative simulations on the deamination of guanine and isocytosine reveal that, for the helping of hydrogen bonds between O8@8-oxoguanine and enzyme, O8@8-oxoguanine is the fastest to be deaminated among the three substrates which are also supported by the experimental kinetic constants.
Collapse
Affiliation(s)
- Xin Zhang
- State Key Laboratory of Chemical Resource Engineering, Institute of Materia Medica, College of Science, Beijing University of Chemical Technology, Beijing 100029, P. R. China
- Department of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, USA
| | - Yuan Zhao
- Department of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, USA
- 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
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng, 475004, P. R. China
| | - Xinli Duan
- State Key Laboratory of Chemical Resource Engineering, Institute of Materia Medica, College of Science, Beijing University of Chemical Technology, Beijing 100029, P. R. China
| | - Hui N. Zhang
- State Key Laboratory of Chemical Resource Engineering, Institute of Materia Medica, College of Science, Beijing University of Chemical Technology, Beijing 100029, P. R. China
| | - Zexing Cao
- 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
| | - Yirong Mo
- Department of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, USA
| |
Collapse
|
33
|
Fogarty AC, Potestio R, Kremer K. A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site. Proteins 2016; 84:1902-1913. [PMID: 27699855 DOI: 10.1002/prot.25173] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2016] [Revised: 09/07/2016] [Accepted: 09/16/2016] [Indexed: 11/11/2022]
Abstract
In multi-resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902-1913. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Aoife C Fogarty
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| |
Collapse
|
34
|
Dong G, Lu J, Lai W. Insights into the Mechanism of Aromatic Ring Cleavage of Noncatecholic Compound 2-Aminophenol by Aminophenol Dioxygenase: A Quantum Mechanics/Molecular Mechanics Study. ACS Catal 2016. [DOI: 10.1021/acscatal.6b00372] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Geng Dong
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Jiarui Lu
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Wenzhen Lai
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| |
Collapse
|
35
|
Zhou Y, Xie D, Zhang Y. Amide Rotation Hindrance Predicts Proteolytic Resistance of Cystine-Knot Peptides. J Phys Chem Lett 2016; 7:1138-42. [PMID: 26958702 PMCID: PMC4824663 DOI: 10.1021/acs.jpclett.6b00373] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Cystine-knot peptides have remarkable stability against protease degradation and are attractive scaffolds for peptide-based therapeutic and diagnostic agents. In this work, by studying the hydrolysis reaction of a cystine-knot inhibitor MCTI-A and its variants with ab initio QM/MM molecular dynamics simulations, we have elucidated an amide rotation hindrance mechanism for proteolysis resistance: The proteolysis of MCTI-A is retarded due to the higher free energy cost during the rotation of NH group around scissile peptide bond at the tetrahedral intermediate of acylation, and covalent constraint provided by disulfide bonds is the key factor to hinder this rotation. A nearly linear correlation has been revealed between free energy barriers of the peptide hydrolysis reaction and the amide rotation free energy changes at the protease-peptide Michaelis complex state. This suggests that amide rotation hindrance could be one useful feature to estimate peptide proteolysis stability.
Collapse
Affiliation(s)
- Yanzi Zhou
- Laboratory of Mesoscopic Chemistry, Collaborative Innovation Center of Chemistry for Life Sciences, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
- Corresponding Author: 1) , 2)
| | - Daiqian Xie
- Laboratory of Mesoscopic Chemistry, Collaborative Innovation Center of Chemistry for Life Sciences, Institute of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003 USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Corresponding Author: 1) , 2)
| |
Collapse
|
36
|
Shen L, Yang W. Quantum Mechanics/Molecular Mechanics Method Combined with Hybrid All-Atom and Coarse-Grained Model: Theory and Application on Redox Potential Calculations. J Chem Theory Comput 2016; 12:2017-27. [PMID: 26930454 DOI: 10.1021/acs.jctc.5b01107] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
We developed a new multiresolution method that spans three levels of resolution with quantum mechanical, atomistic molecular mechanical, and coarse-grained models. The resolution-adapted all-atom and coarse-grained water model, in which an all-atom structural description of the entire system is maintained during the simulations, is combined with the ab initio quantum mechanics and molecular mechanics method. We apply this model to calculate the redox potentials of the aqueous ruthenium and iron complexes by using the fractional number of electrons approach and thermodynamic integration simulations. The redox potentials are recovered in excellent accordance with the experimental data. The speed-up of the hybrid all-atom and coarse-grained water model renders it computationally more attractive. The accuracy depends on the hybrid all-atom and coarse-grained water model used in the combined quantum mechanical and molecular mechanical method. We have used another multiresolution model, in which an atomic-level layer of water molecules around redox center is solvated in supramolecular coarse-grained waters for the redox potential calculations. Compared with the experimental data, this alternative multilayer model leads to less accurate results when used with the coarse-grained polarizable MARTINI water or big multipole water model for the coarse-grained layer.
Collapse
Affiliation(s)
- Lin Shen
- Department of Chemistry, Duke University , Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University , Durham, North Carolina 27708, United States
| |
Collapse
|
37
|
Zhang X, Zhao Y, Yan H, Cao Z, Mo Y. Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD). J Comput Chem 2016; 37:1163-74. [DOI: 10.1002/jcc.24306] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Revised: 01/04/2016] [Accepted: 01/05/2016] [Indexed: 12/16/2022]
Affiliation(s)
- Xin Zhang
- State Key Laboratory of Chemical Resource Engineering, Institute of Materia Medica, College of Science, Beijing University of Chemical Technology; Beijing 100029 China
- Department of Chemistry; Western Michigan University; Kalamazoo Michigan 49008
| | - Yuan Zhao
- 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 China
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University; Kaifeng 475004 China
| | - Honggao Yan
- Department of Biochemistry; The Center for Biological Modeling, Michigan State University; East Lansing Michigan 48824
- Department of Chemistry; The Center for Biological Modeling, Michigan State University; East Lansing Michigan 48824
| | - Zexing Cao
- 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 China
| | - Yirong Mo
- Department of Chemistry; Western Michigan University; Kalamazoo Michigan 49008
| |
Collapse
|
38
|
Lonsdale R, Reetz MT. Reduction of α,β-Unsaturated Ketones by Old Yellow Enzymes: Mechanistic Insights from Quantum Mechanics/Molecular Mechanics Calculations. J Am Chem Soc 2015; 137:14733-42. [PMID: 26521678 DOI: 10.1021/jacs.5b08687] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Enoate reductases catalyze the reduction of activated C═C bonds with high enantioselectivity. The oxidative half-reaction, which involves the addition of a hydride and a proton to opposite faces of the C═C bond, has been studied for the first time by hybrid quantum mechanics/molecular mechanics (QM/MM). The reduction of 2-cyclohexen-1-one by YqjM from Bacillus subtilis was selected as the model system. Two-dimensional QM/MM (B3LYP-D/OPLS2005) reaction pathways suggest that the hydride and proton are added as distinct steps, with the former step preceding the latter. Furthermore, we present interesting insights into the reactivity of this enzyme, including the weak binding of the substrate in the active site, the role of the two active site histidine residues for polarization of the substrate C═O bond, structural details of the transition states to hydride and proton transfer, and the role of Tyr196 as proton donor. The information presented here will be useful for the future design of enantioselective YqjM mutants for other substrates.
Collapse
Affiliation(s)
- Richard Lonsdale
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.,Fachbereich Chemie der Philipps-Universität , Hans-Meerwein-Strasse, 35032 Marburg, Germany
| | - Manfred T Reetz
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.,Fachbereich Chemie der Philipps-Universität , Hans-Meerwein-Strasse, 35032 Marburg, Germany
| |
Collapse
|
39
|
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
|
40
|
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
|
41
|
Dumont E, Monari A. Understanding DNA under oxidative stress and sensitization: the role of molecular modeling. Front Chem 2015; 3:43. [PMID: 26236706 PMCID: PMC4500984 DOI: 10.3389/fchem.2015.00043] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2015] [Accepted: 06/29/2015] [Indexed: 12/12/2022] Open
Abstract
DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.
Collapse
Affiliation(s)
- Elise Dumont
- Laboratoire de Chimie, UMR 5182 Centre National de la Recherche Scientifique, École Normale Supérieure de Lyon Lyon, France
| | - Antonio Monari
- Université de Lorraine - Nancy, Theory-Modeling-Simulation, Structure et Réactivité des Systèmes Moléculaires Complexes (SRSMC) Vandoeuvre-les-Nancy, France ; Centre National de la Recherche Scientifique, Theory-Modeling-Simulation, Structure et Réactivité des Systèmes Moléculaires Complexes (SRSMC) Vandoeuvre-les-Nancy, France
| |
Collapse
|
42
|
Lipid insertion domain unfolding regulates protein orientational transition behavior in a lipid bilayer. Biophys Chem 2015; 206:22-39. [PMID: 26164502 DOI: 10.1016/j.bpc.2015.06.011] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2015] [Revised: 06/24/2015] [Accepted: 06/25/2015] [Indexed: 01/02/2023]
Abstract
We have used coarse-grained (CG) and united atom (UA) molecular dynamics simulations to explore the mechanisms of protein orientational transition of a model peptide (Aβ42) in a phosphatidylcholine/cholesterol (PC/CHO) lipid bilayer. We started with an inserted state of Aβ42 containing a folded (I) or unfolded (II) K28-A42 lipid insertion domain (LID), which was stabilized by the K28-snorkeling and A42-anchoring to the PC polar groups in the lipid bilayer. After a UA-to-CG transformation and a 1000ns-CG simulation for enhancing the sampling of protein orientations, we discovered two transitions: I-to-"deep inserted" state with disrupted K28-snorkeling and II-to-"deep surface" state with disrupted A42-anchoring. The new states remained stable after a CG-to-UA transformation and a 200ns-UA simulation relaxation. Significant changes in the cholesterol-binding domain of Aβ42 and protein-induced membrane disruptions were evident after the transitions. We propose that the conformation of the LID regulates protein orientational transitions in the lipid membrane.
Collapse
|
43
|
Preise der Deutschen Bunsen-Gesellschaft: S. Riniker, C. Bräuchle, H.-J. Freund und S. W. Hell / John Stauffer Lectureship: F. Diederich. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/ange.201504944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
44
|
Deutsche Bunsen-Gesellschaft Awards: S. Riniker, C. Bräuchle, H.-J. Freund, and S. W. Hell / John Stauffer Lectureship: F. Diederich. Angew Chem Int Ed Engl 2015. [DOI: 10.1002/anie.201504944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
45
|
Szklarczyk OM, Arvaniti E, van Gunsteren WF. Polarizable coarse-grained models for molecular dynamics simulation of liquid cyclohexane. J Comput Chem 2015; 36:1311-21. [DOI: 10.1002/jcc.23929] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2014] [Revised: 03/09/2015] [Accepted: 03/12/2015] [Indexed: 11/08/2022]
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
| | - Eirini Arvaniti
- 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
|
46
|
Fogarty AC, Potestio R, Kremer K. Adaptive resolution simulation of a biomolecule and its hydration shell: Structural and dynamical properties. J Chem Phys 2015; 142:195101. [DOI: 10.1063/1.4921347] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Aoife C. Fogarty
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
47
|
Brunk E, Rothlisberger U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem Rev 2015; 115:6217-63. [PMID: 25880693 DOI: 10.1021/cr500628b] [Citation(s) in RCA: 301] [Impact Index Per Article: 33.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Elizabeth Brunk
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,‡Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California 94618, United States
| | - Ursula Rothlisberger
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,§National Competence Center of Research (NCCR) MARVEL-Materials' Revolution: Computational Design and Discovery of Novel Materials, 1015 Lausanne, Switzerland
| |
Collapse
|
48
|
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
|
49
|
Ding HM, Ma YQ. Theoretical and computational investigations of nanoparticle-biomembrane interactions in cellular delivery. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2015; 11:1055-71. [PMID: 25387905 DOI: 10.1002/smll.201401943] [Citation(s) in RCA: 189] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2014] [Revised: 09/05/2014] [Indexed: 05/18/2023]
Abstract
With the rapid development of nanotechnology, nanoparticles have been widely used in many applications such as phototherapy, cell imaging, and drug/gene delivery. A better understanding of how nanoparticles interact with bio-system (especially cells) is of great importance for their potential biomedical applications. In this review, the current status and perspective of theoretical and computational investigations is presented on the nanoparticle-biomembrane interactions in cellular delivery. In particular, the determining parameters (including the properties of nanoparticles, cell membranes and environments) that govern the cellular uptake of nanoparticles (direct penetration and endocytosis) are discussed. Further, some special attention is paid to their interactions beyond the translocation of nanoparticles across membranes (e.g., nanoparticles escaping from endosome and entering into nucleus). Finally, a summary is given, and the challenging problems of this field in the future are identified.
Collapse
Affiliation(s)
- Hong-ming Ding
- National Laboratory of Solid State Microstructures and Department of Physics, Nanjing University, Nanjing, 210093, China
| | | |
Collapse
|
50
|
Lei J, Zhou Y, Xie D, Zhang Y. Mechanistic insights into a classic wonder drug--aspirin. J Am Chem Soc 2015; 137:70-3. [PMID: 25514511 PMCID: PMC4309031 DOI: 10.1021/ja5112964] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2014] [Indexed: 01/24/2023]
Abstract
Aspirin, one of the oldest and most common anti-inflammatory agents, has recently been shown to reduce cancer risks. The principal pharmacological effects of aspirin are known to arise from its covalent modification of cyclooxygenase-2 (COX-2) through acetylation of Ser530, but the detailed mechanism of its biochemical action and specificity remains to be elucidated. In this work, we have filled this gap by employing a state-of-the-art computational approach, Born-Oppenheimer molecular dynamics simulations with ab initio quantum mechanical/molecular mechanical potential and umbrella sampling. Our studies have characterized a substrate-assisted inhibition mechanism for aspirin acetylating COX: it proceeds in two successive stages with a metastable tetrahedral intermediate, in which the carboxyl group of aspirin serves as the general base. The computational results confirmed that aspirin would be 10-100 times more potent against COX-1 than against COX-2, and revealed that this inhibition specificity between the two COX isoforms can be attributed mainly to the difference in kinetics rate of the covalent inhibition reaction, not the aspirin-binding step. The structural origin of this differential inhibition of the COX enzymes by aspirin has also been elucidated.
Collapse
Affiliation(s)
- Jinping Lei
- Institute
of Theoretical and Computational Chemistry, Laboratory of Mesoscopic
Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Yanzi Zhou
- Institute
of Theoretical and Computational Chemistry, Laboratory of Mesoscopic
Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Daiqian Xie
- Institute
of Theoretical and Computational Chemistry, Laboratory of Mesoscopic
Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
- Synergetic
Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026,China
| | - Yingkai Zhang
- Department
of Chemistry, New York University, New York, New York 10003 United States
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|