1
|
Chen H, Ogden D, Pant S, Cai W, Tajkhorshid E, Moradi M, Roux B, Chipot C. A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and Pitfalls. J Chem Theory Comput 2022; 18:1406-1422. [PMID: 35138832 PMCID: PMC8904302 DOI: 10.1021/acs.jctc.1c01049] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The string method with swarms of trajectories (SMwST) is an algorithm that identifies a physically meaningful transition pathway─a one-dimensional curve, embedded within a high-dimensional space of selected collective variables. The SMwST algorithm leans on a series of short, unbiased molecular dynamics simulations spawned at different locations of the discretized path, from whence an average dynamic drift is determined to evolve the string toward an optimal pathway. However conceptually simple in both its theoretical formulation and practical implementation, the SMwST algorithm is computationally intensive and requires a careful choice of parameters for optimal cost-effectiveness in applications to challenging problems in chemistry and biology. In this contribution, the SMwST algorithm is presented in a self-contained manner, discussing with a critical eye its theoretical underpinnings, applicability, inherent limitations, and use in the context of path-following free-energy calculations and their possible extension to kinetics modeling. Through multiple simulations of a prototypical polypeptide, combining the search of the transition pathway and the computation of the potential of mean force along it, several practical aspects of the methodology are examined with the objective of optimizing the computational effort, yet without sacrificing accuracy. In light of the results reported here, we propose some general guidelines aimed at improving the efficiency and reliability of the computed pathways and free-energy profiles underlying the conformational transitions at hand.
Collapse
Affiliation(s)
- Haochuan Chen
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Nankai University, Tianjin 300071, China
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
| | - Dylan Ogden
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Shashank Pant
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Wensheng Cai
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Nankai University, Tianjin 300071, China
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Biochemistry and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Mahmoud Moradi
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
| | - Christophe Chipot
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
2
|
Abstract
![]()
The kinetics of
a dynamical system comprising two metastable states
is formulated in terms of a finite-time propagator in phase space
(position and velocity) adapted to the underdamped Langevin equation.
Dimensionality reduction to a subspace of collective variables yields
familiar expressions for the propagator, committor, and steady-state
flux. A quadratic expression for the steady-state flux between the
two metastable states can serve as a robust variational principle
to determine an optimal approximate committor expressed in terms of
a set of collective variables. The theoretical formulation is exploited
to clarify the foundation of the string method with swarms-of-trajectories,
which relies on the mean drift of short trajectories to determine
the optimal transition pathway. It is argued that the conditions for
Markovity within a subspace of collective variables may not be satisfied
with an arbitrary short time-step and that proper kinetic behaviors
appear only when considering the effective propagator for longer lag
times. The effective propagator with finite lag time is amenable to
an eigenvalue-eigenvector spectral analysis, as elaborated previously
in the context of position-based Markov models. The time-correlation
functions calculated by swarms-of-trajectories along the string pathway
constitutes a natural extension of these developments. The present
formulation provides a powerful theoretical framework to characterize
the optimal pathway between two metastable states of a system.
Collapse
Affiliation(s)
- Benoît Roux
- Department of Biochemistry and Molecular Biology, The University of Chicago, Chicago, Illinois 60637, United States.,Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Maschio MC, Fregoni J, Molteni C, Corni S. Proline isomerization effects in the amyloidogenic protein β2-microglobulin. Phys Chem Chem Phys 2021; 23:356-367. [DOI: 10.1039/d0cp04780e] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The protein β2-microglobulin can aggregate in insoluble amyloid fibrils. By relying on extensive sampling simulations, we study the Pro32 isomerization as a possible triggering factor leading to structural modifications in β2-m.
Collapse
Affiliation(s)
| | - Jacopo Fregoni
- CNR-Nano S3
- Modena
- Italy
- Department of Chemical Sciences
- University of Padova
| | - Carla Molteni
- Department of Physics
- King's College London
- Strand
- London WC2R 2LS
- UK
| | - Stefano Corni
- CNR-Nano S3
- Modena
- Italy
- Department of Chemical Sciences
- University of Padova
| |
Collapse
|
4
|
Lev B, Allen TW. Simulating ion channel activation mechanisms using swarms of trajectories. J Comput Chem 2020; 41:387-401. [PMID: 31743478 DOI: 10.1002/jcc.26102] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Revised: 10/16/2019] [Accepted: 10/17/2019] [Indexed: 12/14/2022]
Abstract
Atomic-level studies of protein activity represent a significant challenge as a result of the complexity of conformational changes occurring on wide-ranging timescales, often greatly exceeding that of even the longest simulations. A prime example is the elucidation of protein allosteric mechanisms, where localized perturbations transmit throughout a large macromolecule to generate a response signal. For example, the conversion of chemical to electrical signals during synaptic neurotransmission in the brain is achieved by specialized membrane proteins called pentameric ligand-gated ion channels. Here, the binding of a neurotransmitter results in a global conformational change to open an ion-conducting pore across the nerve cell membrane. X-ray crystallography has produced static structures of the open and closed states of the proton-gated GLIC pentameric ligand-gated ion channel protein, allowing for atomistic simulations that can uncover changes related to activation. We discuss a range of enhanced sampling approaches that could be used to explore activation mechanisms. In particular, we describe recent application of an atomistic string method, based on Roux's "swarms of trajectories" approach, to elucidate the sequence and interdependence of conformational changes during activation. We illustrate how this can be combined with transition analysis and Brownian dynamics to extract thermodynamic and kinetic information, leading to understanding of what controls ion channel function. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Bogdan Lev
- School of Science, RMIT University, Melbourne, Victoria, 3000, Australia
| | - Toby W Allen
- School of Science, RMIT University, Melbourne, Victoria, 3000, Australia
| |
Collapse
|
5
|
Flood E, Boiteux C, Lev B, Vorobyov I, Allen TW. Atomistic Simulations of Membrane Ion Channel Conduction, Gating, and Modulation. Chem Rev 2019; 119:7737-7832. [DOI: 10.1021/acs.chemrev.8b00630] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Emelie Flood
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Céline Boiteux
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Bogdan Lev
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| | - Igor Vorobyov
- Department of Physiology & Membrane Biology/Department of Pharmacology, University of California, Davis, 95616, United States
| | - Toby W. Allen
- School of Science, RMIT University, Melbourne, Victoria 3000, Australia
| |
Collapse
|
6
|
Fujisaki H, Moritsugu K, Matsunaga Y. Exploring Configuration Space and Path Space of Biomolecules Using Enhanced Sampling Techniques-Searching for Mechanism and Kinetics of Biomolecular Functions. Int J Mol Sci 2018; 19:E3177. [PMID: 30326661 PMCID: PMC6213965 DOI: 10.3390/ijms19103177] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 10/10/2018] [Accepted: 10/11/2018] [Indexed: 01/07/2023] Open
Abstract
To understand functions of biomolecules such as proteins, not only structures but their conformational change and kinetics need to be characterized, but its atomistic details are hard to obtain both experimentally and computationally. Here, we review our recent computational studies using novel enhanced sampling techniques for conformational sampling of biomolecules and calculations of their kinetics. For efficiently characterizing the free energy landscape of a biomolecule, we introduce the multiscale enhanced sampling method, which uses a combined system of atomistic and coarse-grained models. Based on the idea of Hamiltonian replica exchange, we can recover the statistical properties of the atomistic model without any biases. We next introduce the string method as a path search method to calculate the minimum free energy pathways along a multidimensional curve in high dimensional space. Finally we introduce novel methods to calculate kinetics of biomolecules based on the ideas of path sampling: one is the Onsager⁻Machlup action method, and the other is the weighted ensemble method. Some applications of the above methods to biomolecular systems are also discussed and illustrated.
Collapse
Grants
- JPMJPR1679 Japan Science and Technology Agency
- 16K00059 Ministry of Education, Culture, Sports, Science and Technology
- 17KT0101 Ministry of Education, Culture, Sports, Science and Technology
- 25840060 Ministry of Education, Culture, Sports, Science and Technology
- 15K18520 Ministry of Education, Culture, Sports, Science and Technology
- JP18am0101109 Japan Agency for Medical Research and Development
- 17gm0810012h0001 Japan Agency for Medical Research and Development
Collapse
Affiliation(s)
- Hiroshi Fujisaki
- Department of Physics, Nippon Medical School, 1-7-1 Kyonan-cho, Musashino, Tokyo 180-0023, Japan.
- AMED-CREST, Japan Agency for Medical Research and Development, 1-1-5 Sendagi, Bunkyo-ku, Tokyo 113-8603, Japan.
| | - Kei Moritsugu
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.
| | - Yasuhiro Matsunaga
- RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.
- JST PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan.
| |
Collapse
|
7
|
Matsunaga Y, Yamane T, Terada T, Moritsugu K, Fujisaki H, Murakami S, Ikeguchi M, Kidera A. Energetics and conformational pathways of functional rotation in the multidrug transporter AcrB. eLife 2018; 7:31715. [PMID: 29506651 PMCID: PMC5839741 DOI: 10.7554/elife.31715] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2017] [Accepted: 02/12/2018] [Indexed: 12/31/2022] Open
Abstract
The multidrug transporter AcrB transports a broad range of drugs out of the cell by means of the proton-motive force. The asymmetric crystal structure of trimeric AcrB suggests a functionally rotating mechanism for drug transport. Despite various supportive forms of evidence from biochemical and simulation studies for this mechanism, the link between the functional rotation and proton translocation across the membrane remains elusive. Here, calculating the minimum free energy pathway of the functional rotation for the complete AcrB trimer, we describe the structural and energetic basis behind the coupling between the functional rotation and the proton translocation at atomic resolution. Free energy calculations show that protonation of Asp408 in the transmembrane portion of the drug-bound protomer drives the functional rotation. The conformational pathway identifies vertical shear motions among several transmembrane helices, which regulate alternate access of water in the transmembrane as well as peristaltic motions that pump drugs in the periplasm.
Collapse
Affiliation(s)
- Yasuhiro Matsunaga
- RIKEN Advanced Institute for Computational Science, Kobe, Japan.,JST PRESTO, Kawaguchi, Japan
| | - Tsutomu Yamane
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Japan
| | - Tohru Terada
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Kei Moritsugu
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Japan
| | | | - Satoshi Murakami
- Graduate School of Bioscience & Biotechnology, Tokyo Institute of Technology, Yokohama, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Japan
| | | |
Collapse
|
8
|
Maximova T, Moffatt R, Ma B, Nussinov R, Shehu A. Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics. PLoS Comput Biol 2016; 12:e1004619. [PMID: 27124275 PMCID: PMC4849799 DOI: 10.1371/journal.pcbi.1004619] [Citation(s) in RCA: 132] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.
Collapse
Affiliation(s)
- Tatiana Maximova
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Ryan Moffatt
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
| | - Buyong Ma
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
| | - Ruth Nussinov
- Basic Science Program, Leidos Biomedical Research, Inc. Cancer and Inflammation Program, National Cancer Institute, Frederick, Maryland, United States of America
- Sackler Institute of Molecular Medicine, Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Tel Aviv University, Tel Aviv, Israel
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, Virginia, United States of America
- Department of Biongineering, George Mason University, Fairfax, Virginia, United States of America
- School of Systems Biology, George Mason University, Manassas, Virginia, United States of America
| |
Collapse
|
9
|
Fujisaki H, Moritsugu K, Matsunaga Y, Morishita T, Maragliano L. Extended Phase-Space Methods for Enhanced Sampling in Molecular Simulations: A Review. Front Bioeng Biotechnol 2015; 3:125. [PMID: 26389113 PMCID: PMC4558547 DOI: 10.3389/fbioe.2015.00125] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2015] [Accepted: 08/10/2015] [Indexed: 01/12/2023] Open
Abstract
Molecular Dynamics simulations are a powerful approach to study biomolecular conformational changes or protein-ligand, protein-protein, and protein-DNA/RNA interactions. Straightforward applications, however, are often hampered by incomplete sampling, since in a typical simulated trajectory the system will spend most of its time trapped by high energy barriers in restricted regions of the configuration space. Over the years, several techniques have been designed to overcome this problem and enhance space sampling. Here, we review a class of methods that rely on the idea of extending the set of dynamical variables of the system by adding extra ones associated to functions describing the process under study. In particular, we illustrate the Temperature Accelerated Molecular Dynamics (TAMD), Logarithmic Mean Force Dynamics (LogMFD), and Multiscale Enhanced Sampling (MSES) algorithms. We also discuss combinations with techniques for searching reaction paths. We show the advantages presented by this approach and how it allows to quickly sample important regions of the free-energy landscape via automatic exploration.
Collapse
Affiliation(s)
| | - Kei Moritsugu
- Graduate School of Medical Life Science, Yokohama City University, Yokohama, Japan
| | | | - Tetsuya Morishita
- National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan
| | - Luca Maragliano
- Center for Synaptic Neuroscience, Istituto Italiano di Tecnologia, Genova, Italy
| |
Collapse
|
10
|
Brancolini G, Corazza A, Vuano M, Fogolari F, Mimmi MC, Bellotti V, Stoppini M, Corni S, Esposito G. Probing the influence of citrate-capped gold nanoparticles on an amyloidogenic protein. ACS NANO 2015; 9:2600-13. [PMID: 25695203 DOI: 10.1021/nn506161j] [Citation(s) in RCA: 67] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Nanoparticles (NPs) are known to exhibit distinct physical and chemical properties compared with the same materials in bulk form. NPs have been repeatedly reported to interact with proteins, and this interaction can be exploited to affect processes undergone by proteins, such as fibrillogenesis. Fibrillation is common to many proteins, and in living organisms, it causes tissue-specific or systemic amyloid diseases. The nature of NPs and their surface chemistry is crucial in assessing their affinity for proteins and their effects on them. Here we present the first detailed structural characterization and molecular mechanics model of the interaction between a fibrillogenic protein, β2-microglobulin, and a NP, 5 nm hydrophilic citrate-capped gold nanoparticles. NMR measurements and simulations at multiple levels (enhanced sampling molecular dynamics, Brownian dynamics, and Poisson-Boltzmann electrostatics) explain the origin of the observed protein perturbations mostly localized at the amino-terminal region. Experiments show that the protein-NP interaction is weak in the physiological-like, conditions and do not induce protein fibrillation. Simulations reproduce these findings and reveal instead the role of the citrate in destabilizing the lower pH protonated form of β2-microglobulin. The results offer possible strategies for controlling the desired effect of NPs on the conformational changes of the proteins, which have significant roles in the fibrillation process.
Collapse
Affiliation(s)
- Giorgia Brancolini
- †Center S3, CNR Institute Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | - Alessandra Corazza
- ‡Dipartimento di Scienze Mediche e Biologiche (DSMB), University of Udine, Piazzale Kolbe 4, 33100 Udine, Italy
- §Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, 00136 Roma, Italy
| | - Marco Vuano
- ‡Dipartimento di Scienze Mediche e Biologiche (DSMB), University of Udine, Piazzale Kolbe 4, 33100 Udine, Italy
| | - Federico Fogolari
- ‡Dipartimento di Scienze Mediche e Biologiche (DSMB), University of Udine, Piazzale Kolbe 4, 33100 Udine, Italy
- §Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, 00136 Roma, Italy
| | - Maria Chiara Mimmi
- ‡Dipartimento di Scienze Mediche e Biologiche (DSMB), University of Udine, Piazzale Kolbe 4, 33100 Udine, Italy
| | - Vittorio Bellotti
- §Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, 00136 Roma, Italy
- ⊥Dipartimento di Medicina Molecolare, Universita' di Pavia, Via Taramelli 3, 27100 Pavia, Italy
- ∥Division of Medicine, University College of London, London NW3 2PF, U.K
| | - Monica Stoppini
- §Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, 00136 Roma, Italy
- ⊥Dipartimento di Medicina Molecolare, Universita' di Pavia, Via Taramelli 3, 27100 Pavia, Italy
| | - Stefano Corni
- †Center S3, CNR Institute Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | - Gennaro Esposito
- ‡Dipartimento di Scienze Mediche e Biologiche (DSMB), University of Udine, Piazzale Kolbe 4, 33100 Udine, Italy
- §Istituto Nazionale Biostrutture e Biosistemi, Viale medaglie d'Oro 305, 00136 Roma, Italy
- ¶Science and Math Division, New York University at Abu Dhabi, Abu Dhabi, UAE
| |
Collapse
|
11
|
Estácio SG, Krobath H, Vila-Viçosa D, Machuqueiro M, Shakhnovich EI, Faísca PFN. A simulated intermediate state for folding and aggregation provides insights into ΔN6 β2-microglobulin amyloidogenic behavior. PLoS Comput Biol 2014; 10:e1003606. [PMID: 24809460 PMCID: PMC4014404 DOI: 10.1371/journal.pcbi.1003606] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2013] [Accepted: 03/16/2014] [Indexed: 01/26/2023] Open
Abstract
A major component of ex vivo amyloid plaques of patients with dialysis-related amyloidosis (DRA) is a cleaved variant of β2-microglobulin (ΔN6) lacking the first six N-terminal residues. Here we perform a computational study on ΔN6, which provides clues to understand the amyloidogenicity of the full-length β2-microglobulin. Contrary to the wild-type form, ΔN6 is able to efficiently nucleate fibrillogenesis in vitro at physiological pH. This behavior is enhanced by a mild acidification of the medium such as that occurring in the synovial fluid of DRA patients. Results reported in this work, based on molecular simulations, indicate that deletion of the N-terminal hexapeptide triggers the formation of an intermediate state for folding and aggregation with an unstructured strand A and a native-like core. Strand A plays a pivotal role in aggregation by acting as a sticky hook in dimer assembly. This study further predicts that the detachment of strand A from the core is maximized at pH 6.2 resulting into higher aggregation efficiency. The structural mapping of the dimerization interface suggests that Tyr10, His13, Phe30 and His84 are hot-spot residues in ΔN6 amyloidogenesis. Dialysis-related amyloidosis (DRA) is a conformational disease that affects individuals undergoing long-term haemodialysis. In DRA the progressive accumulation of protein human β2-microglobulin (Hβ2m) in the osteoarticular system, followed by its assembly into amyloid fibrils, eventually leads to tissue erosion and destruction. Disclosing the aggregation mechanism of Hβ2m under physiologically relevant conditions represents a major challenge due to the inability of the protein to efficiently nucleate fibrillogenesis in vitro at physiological pH. On the other hand, ΔN6, a truncated variant of Hβ2m, which is also a major component of ex vivo amyloid deposits extracted from DRA patients, is able to efficiently form amyloid fibrils de novo in physiological conditions. This amyloidogenic behavior is dramatically enhanced in a slightly more acidic pH (6.2) compatible with the mild acidification that occurs in the synovial fluid of DRA patients. In this work, an innovative three-stage methodological approach, relying on an array of molecular simulations, spanning different levels of resolution is used to investigate the initial stage of the de novo aggregation mechanism of ΔN6 in a physiologically relevant pH range. We identify an intermediate state for folding and aggregation, whose potential to dimerize is enhanced at pH 6.2. Our results provide rationalizations for previous experimental observations and new insights into the molecular basis of DRA.
Collapse
Affiliation(s)
- Sílvia G. Estácio
- Centro de Física da Matéria Condensada & Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| | - Heinrich Krobath
- Centro de Física da Matéria Condensada & Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| | - Diogo Vila-Viçosa
- Centro de Química e Bioquímica & Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| | - Miguel Machuqueiro
- Centro de Química e Bioquímica & Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
| | - Eugene I. Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts, United States of America
- * E-mail: (EIS); (PFNF)
| | - Patrícia F. N. Faísca
- Centro de Física da Matéria Condensada & Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, Lisboa, Portugal
- * E-mail: (EIS); (PFNF)
| |
Collapse
|
12
|
Das A, Gur M, Cheng MH, Jo S, Bahar I, Roux B. Exploring the conformational transitions of biomolecular systems using a simple two-state anisotropic network model. PLoS Comput Biol 2014; 10:e1003521. [PMID: 24699246 PMCID: PMC3974643 DOI: 10.1371/journal.pcbi.1003521] [Citation(s) in RCA: 96] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2013] [Accepted: 02/01/2014] [Indexed: 11/19/2022] Open
Abstract
Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results. Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.
Collapse
Affiliation(s)
- Avisek Das
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America
| | - Mert Gur
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Mary Hongying Cheng
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Sunhwan Jo
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America
| | - Ivet Bahar
- Department of Computational & Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois, United States of America
- * E-mail:
| |
Collapse
|
13
|
Vashisth H, Skiniotis G, Brooks CL. Collective variable approaches for single molecule flexible fitting and enhanced sampling. Chem Rev 2014; 114:3353-65. [PMID: 24446720 PMCID: PMC3983124 DOI: 10.1021/cr4005988] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2013] [Indexed: 12/12/2022]
Affiliation(s)
- Harish Vashisth
- Department
of Chemical Engineering, University of New
Hampshire, Durham, New Hampshire 03824, United States
| | - Georgios Skiniotis
- Life Sciences Institute, Department
of Biological Chemistry, and
Biophysics Program, and Department of Chemistry and Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Charles Lee Brooks
- Life Sciences Institute, Department
of Biological Chemistry, and
Biophysics Program, and Department of Chemistry and Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
14
|
Maragliano L, Roux B, Vanden-Eijnden E. Comparison between Mean Forces and Swarms-of-Trajectories String Methods. J Chem Theory Comput 2014; 10:524-33. [PMID: 26580029 PMCID: PMC6980172 DOI: 10.1021/ct400606c] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The original formulation of the string method in collective variable space is compared with a recent variant called string method with swarms-of-trajectories. The assumptions made in the original method are revisited and the significance of the minimum free energy path (MFEP) is discussed in the context of reactive events. These assumptions are compared to those made in the string method with swarms-of-trajectories, and shown to be equivalent in a certain regime: in particular an expression for the path identified by the swarms-of-trajectories method is given and shown to be closely related to the MFEP. Finally, the algorithmic aspects of both methods are compared.
Collapse
Affiliation(s)
- Luca Maragliano
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
- Biosciences Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Eric Vanden-Eijnden
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, United States
| |
Collapse
|
15
|
Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration. ENTROPY 2013. [DOI: 10.3390/e16010163] [Citation(s) in RCA: 282] [Impact Index Per Article: 25.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
16
|
Zinovjev K, Ruiz-Pernía JJ, Tuñón I. Toward an Automatic Determination of Enzymatic Reaction Mechanisms and Their Activation Free Energies. J Chem Theory Comput 2013; 9:3740-9. [PMID: 26584125 DOI: 10.1021/ct400153r] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present a combination of the string method and a path collective variable for the exploration of the free energy surface associated to a chemical reaction in condensed environments. The on-the-fly string method is employed to find the minimum free energy paths on a multidimensional free energy surface defined in terms of interatomic distances, which is a convenient selection to study bond forming/breaking processes. Once the paths have been determined, a reaction coordinate is defined as a measure of the advance of the system along these paths. This reaction coordinate can be then used to trace the reaction Potential of Mean Force from which the activation free energy can be obtained. This combination of methodologies has been here applied to the study, by means of Quantum Mechanics/Molecular Mechanics simulations, of the reaction catalyzed by guanidinoacetate methyltransferase. This enzyme catalyzes the methylation of guanidinoacetate by S-adenosyl-l-methionine, a reaction that involves a methyl transfer and a proton transfer and for which different reaction mechanisms have been proposed.
Collapse
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física, Universitat de València , 46100 Burjassot, Spain.,Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - J Javier Ruiz-Pernía
- Departament de Química Física i Analítica, Universitat Jaume I , 12071 Castellón, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València , 46100 Burjassot, Spain
| |
Collapse
|