101
|
Ariyaratne A, Wu C, Tseng CY, Zocchi G. Dissipative dynamics of enzymes. PHYSICAL REVIEW LETTERS 2014; 113:198101. [PMID: 25415926 DOI: 10.1103/physrevlett.113.198101] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/08/2014] [Indexed: 06/04/2023]
Abstract
We explore enzyme conformational dynamics at sub-Å resolution, specifically, temperature effects. The ensemble-averaged mechanical response of the folded enzyme is viscoelastic in the whole temperature range between the warm and cold denaturation transitions. The dissipation parameter γ of the viscoelastic description decreases by a factor of 2 as the temperature is raised from 10 to 45 °C; the elastic parameter K shows a similar decrease. Thus, when probed dynamically, the enzyme softens for increasing temperature. Equilibrium mechanical experiments with the DNA spring (and a different enzyme) also show, qualitatively, a small softening for increasing temperature.
Collapse
Affiliation(s)
- Amila Ariyaratne
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, USA
| | - Chenhao Wu
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, USA
| | - Chiao-Yu Tseng
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, USA
| | - Giovanni Zocchi
- Department of Physics and Astronomy, University of California, Los Angeles, California 90095, USA
| |
Collapse
|
102
|
Moritsugu K, Terada T, Kidera A. Multiscale enhanced sampling driven by multiple coarse-grained models. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.10.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
103
|
Fuglebakk E, Tiwari SP, Reuter N. Comparing the intrinsic dynamics of multiple protein structures using elastic network models. Biochim Biophys Acta Gen Subj 2014; 1850:911-922. [PMID: 25267310 DOI: 10.1016/j.bbagen.2014.09.021] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Revised: 09/15/2014] [Accepted: 09/16/2014] [Indexed: 12/15/2022]
Abstract
BACKGROUND Elastic network models (ENMs) are based on the simple idea that a protein can be described as a set of particles connected by springs, which can then be used to describe its intrinsic flexibility using, for example, normal mode analysis. Since the introduction of the first ENM by Monique Tirion in 1996, several variants using coarser protein models have been proposed and their reliability for the description of protein intrinsic dynamics has been widely demonstrated. Lately an increasing number of studies have focused on the meaning of slow dynamics for protein function and its potential conservation through evolution. This leads naturally to comparisons of the intrinsic dynamics of multiple protein structures with varying levels of similarity. SCOPE OF REVIEW We describe computational strategies for calculating and comparing intrinsic dynamics of multiple proteins using elastic network models, as well as a selection of examples from the recent literature. MAJOR CONCLUSIONS The increasing interest for comparing dynamics across protein structures with various levels of similarity, has led to the establishment and validation of reliable computational strategies using ENMs. Comparing dynamics has been shown to be a viable way for gaining greater understanding for the mechanisms employed by proteins for their function. Choices of ENM parameters, structure alignment or similarity measures will likely influence the interpretation of the comparative analysis of protein motion. GENERAL SIGNIFICANCE Understanding the relation between protein function and dynamics is relevant to the fundamental understanding of protein structure-dynamics-function relationship. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Edvin Fuglebakk
- Department of Molecular Biology, University of Bergen, Pb. 7803, N-5020 Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Pb. 7803, N-5020 Bergen, Norway.
| | - Sandhya P Tiwari
- Department of Molecular Biology, University of Bergen, Pb. 7803, N-5020 Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Pb. 7803, N-5020 Bergen, Norway.
| | - Nathalie Reuter
- Department of Molecular Biology, University of Bergen, Pb. 7803, N-5020 Bergen, Norway; Computational Biology Unit, Department of Informatics, University of Bergen, Pb. 7803, N-5020 Bergen, Norway.
| |
Collapse
|
104
|
Tavanti F, Tozzini V. A multi-scale-multi-stable model for the rhodopsin photocycle. Molecules 2014; 19:14961-78. [PMID: 25237751 PMCID: PMC6271392 DOI: 10.3390/molecules190914961] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2014] [Revised: 08/28/2014] [Accepted: 09/08/2014] [Indexed: 11/19/2022] Open
Abstract
We report a multi-scale simulation study of the photocycle of the rhodopsins. The quasi-atomistic representation ("united atoms" UA) of retinal is combined with a minimalist coarse grained (CG, one-bead-per amino acid) representation of the protein, in a hybrid UA/CG approach, which is the homolog of QM/MM, but at lower resolution. An accurate multi-stable parameterization of the model allows simulating each state and transition among them, and the combination of different scale representation allows addressing the entire photocycle. We test the model on bacterial rhodopsin, for which more experimental data are available, and then also report results for mammalian rhodopsins. In particular, the analysis of simulations reveals the spontaneous appearance of meta-stable states in quantitative agreement with experimental data.
Collapse
Affiliation(s)
- Francesco Tavanti
- NEST-Istituto Nanoscienze, CNR, Piazza San Silvestro 12, 56127 Pisa, Italy
| | - Valentina Tozzini
- NEST-Istituto Nanoscienze, CNR, Piazza San Silvestro 12, 56127 Pisa, Italy.
| |
Collapse
|
105
|
Noel JK, Whitford PC. How Simulations Reveal Dynamics, Disorder, and the Energy Landscapes of Biomolecular Function. Isr J Chem 2014. [DOI: 10.1002/ijch.201400018] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
106
|
Energy landscape views for interplays among folding, binding, and allostery of calmodulin domains. Proc Natl Acad Sci U S A 2014; 111:10550-5. [PMID: 25002491 DOI: 10.1073/pnas.1402768111] [Citation(s) in RCA: 110] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Ligand binding modulates the energy landscape of proteins, thus altering their folding and allosteric conformational dynamics. To investigate such interplay, calmodulin has been a model protein. Despite much attention, fully resolved mechanisms of calmodulin folding/binding have not been elucidated. Here, by constructing a computational model that can integrate folding, binding, and allosteric motions, we studied in-depth folding of isolated calmodulin domains coupled with binding of two calcium ions and associated allosteric conformational changes. First, mechanically pulled simulations revealed coexistence of three distinct conformational states: the unfolded, the closed, and the open states, which is in accord with and augments structural understanding of recent single-molecule experiments. Second, near the denaturation temperature, we found the same three conformational states as well as three distinct binding states: zero, one, and two calcium ion bound states, leading to as many as nine states. Third, in terms of the nine-state representation, we found multiroute folding/binding pathways and shifts in their probabilities with the calcium concentration. At a lower calcium concentration, "combined spontaneous folding and induced fit" occurs, whereas at a higher concentration, "binding-induced folding" dominates. Even without calcium binding, we observed that the folding pathway of calmodulin domains can be modulated by the presence of metastable states. Finally, full-length calmodulin also exhibited an intriguing coupling between two domains when applying tension.
Collapse
|
107
|
Leioatts N, Suresh P, Romo TD, Grossfield A. Structure-based simulations reveal concerted dynamics of GPCR activation. Proteins 2014; 82:2538-51. [PMID: 24889093 DOI: 10.1002/prot.24617] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 05/06/2014] [Accepted: 05/20/2014] [Indexed: 11/08/2022]
Abstract
G protein-coupled receptors (GPCRs) are a vital class of proteins that transduce biological signals across the cell membrane. However, their allosteric activation mechanism is not fully understood; crystal structures of active and inactive receptors have been reported, but the functional pathway between these two states remains elusive. Here, we use structure-based (Gō-like) models to simulate activation of two GPCRs, rhodopsin and the β₂ adrenergic receptor (β₂AR). We used data-derived reaction coordinates that capture the activation mechanism for both proteins, showing that activation proceeds through quantitatively different paths in the two systems. Both reaction coordinates are determined from the dominant concerted motions in the simulations so the technique is broadly applicable. There were two surprising results. First, the main structural changes in the simulations were distributed throughout the transmembrane bundle, and not localized to the obvious areas of interest, such as the intracellular portion of Helix 6. Second, the activation (and deactivation) paths were distinctly nonmonotonic, populating states that were not simply interpolations between the inactive and active structures. These transitions also suggest a functional explanation for β₂AR's basal activity: it can proceed through a more broadly defined path during the observed transitions.
Collapse
Affiliation(s)
- Nicholas Leioatts
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, New York, 14642
| | | | | | | |
Collapse
|
108
|
Zhou Y, Pu J. Reaction Path Force Matching: A New Strategy of Fitting Specific Reaction Parameters for Semiempirical Methods in Combined QM/MM Simulations. J Chem Theory Comput 2014; 10:3038-54. [PMID: 26588275 DOI: 10.1021/ct4009624] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a general strategy of reparametrizing semiempirical (SE) methods against ab initio (AI) methods for combined quantum mechanical and molecular mechanical (QM/MM) simulations of specific chemical reactions in condensed phases. The resulting approach, designated Reaction Path Force Matching (RP-FM), features cycles of sampling configurations along a reaction path on an efficient SE/MM potential energy surface (PES) and adjusting specific reaction parameters (SRPs) in the SE method such that the atomic forces computed at the target AI/MM level are reproduced. Iterative applications of the RP-FM cycle make possible achieving the accuracy of AI/MM simulations without explicitly sampling the computationally expensive AI/MM PES. The bypassed sampling, nevertheless, is implicitly accomplished through the aid of the efficient SE-SRP/MM PES, on which the target-level reaction path is expected to be obtained upon convergence. We demonstrate the effectiveness of the RP-FM procedure for a symmetric proton transfer reaction in the gas phase and in solution. The remarkable agreements between the RP-FM optimized SE-SRP methods and the target AI method on various properties, including energy profiles, potential of mean force free energy profiles, atomic forces, charge populations, and solvation effects, suggest that RP-FM can be used as an efficient and reliable strategy for simulating condensed-phase chemical reactions.
Collapse
Affiliation(s)
- Yan Zhou
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis , 402 N. Blackford St., LD326, Indianapolis, Indiana 46202, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis , 402 N. Blackford St., LD326, Indianapolis, Indiana 46202, United States
| |
Collapse
|
109
|
Seyler SL, Beckstein O. Sampling large conformational transitions: adenylate kinase as a testing ground. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.919497] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
110
|
Tempkin JOB, Qi B, Saunders MG, Roux B, Dinner AR, Weare J. Using multiscale preconditioning to accelerate the convergence of iterative molecular calculations. J Chem Phys 2014; 140:184114. [PMID: 24832260 PMCID: PMC11450774 DOI: 10.1063/1.4872021] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2013] [Accepted: 04/09/2014] [Indexed: 11/14/2022] Open
Abstract
Iterative procedures for optimizing properties of molecular models often converge slowly owing to the computational cost of accurately representing features of interest. Here, we introduce a preconditioning scheme that allows one to use a less expensive model to guide exploration of the energy landscape of a more expensive model and thus speed the discovery of locally stable states of the latter. We illustrate our approach in the contexts of energy minimization and the string method for finding transition pathways. The relation of the method to other multilevel simulation techniques and possible extensions are discussed.
Collapse
Affiliation(s)
- Jeremy O B Tempkin
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Bo Qi
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Marissa G Saunders
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Benoit Roux
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Aaron R Dinner
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, USA
| | - Jonathan Weare
- James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
111
|
Knott M, Best RB. Discriminating binding mechanisms of an intrinsically disordered protein via a multi-state coarse-grained model. J Chem Phys 2014; 140:175102. [PMID: 24811666 PMCID: PMC4032430 DOI: 10.1063/1.4873710] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2014] [Accepted: 04/17/2014] [Indexed: 12/31/2022] Open
Abstract
Many proteins undergo a conformational transition upon binding to their cognate binding partner, with intrinsically disordered proteins (IDPs) providing an extreme example in which a folding transition occurs. However, it is often not clear whether this occurs via an "induced fit" or "conformational selection" mechanism, or via some intermediate scenario. In the first case, transient encounters with the binding partner favour transitions to the bound structure before the two proteins dissociate, while in the second the bound structure must be selected from a subset of unbound structures which are in the correct state for binding, because transient encounters of the incorrect conformation with the binding partner are most likely to result in dissociation. A particularly interesting situation involves those intrinsically disordered proteins which can bind to different binding partners in different conformations. We have devised a multi-state coarse-grained simulation model which is able to capture the binding of IDPs in alternate conformations, and by applying it to the binding of nuclear coactivator binding domain (NCBD) to either ACTR or IRF-3 we are able to determine the binding mechanism. By all measures, the binding of NCBD to either binding partner appears to occur via an induced fit mechanism. Nonetheless, we also show how a scenario closer to conformational selection could arise by choosing an alternative non-binding structure for NCBD.
Collapse
Affiliation(s)
- Michael Knott
- Department of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Robert B Best
- Department of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
112
|
Gur M, Madura JD, Bahar I. Global transitions of proteins explored by a multiscale hybrid methodology: application to adenylate kinase. Biophys J 2014; 105:1643-52. [PMID: 24094405 DOI: 10.1016/j.bpj.2013.07.058] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Revised: 06/25/2013] [Accepted: 07/11/2013] [Indexed: 12/19/2022] Open
Abstract
Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.
Collapse
Affiliation(s)
- Mert Gur
- Department of Computational and Systems Biology, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania
| | | | | |
Collapse
|
113
|
Li CH, Yang YX, Su JG, Liu B, Tan JJ, Zhang XY, Wang CX. Allosteric transitions of the maltose transporter studied by an elastic network model. Biopolymers 2014; 101:758-68. [DOI: 10.1002/bip.22455] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2013] [Revised: 12/02/2013] [Accepted: 12/09/2013] [Indexed: 11/11/2022]
Affiliation(s)
- Chun Hua Li
- College of Life Science and Bioengineering; Beijing University of Technology; Beijing 100124 China
| | - Yong Xiao Yang
- College of Life Science and Bioengineering; Beijing University of Technology; Beijing 100124 China
| | - Ji Guo Su
- College of Science; Yanshan University; Qinhuangdao 066004 China
| | - Bin Liu
- College of Life Science and Bioengineering; Beijing University of Technology; Beijing 100124 China
| | - Jian Jun Tan
- College of Life Science and Bioengineering; Beijing University of Technology; Beijing 100124 China
| | - Xiao Yi Zhang
- College of Life Science and Bioengineering; Beijing University of Technology; Beijing 100124 China
| | - Cun Xin Wang
- College of Life Science and Bioengineering; Beijing University of Technology; Beijing 100124 China
| |
Collapse
|
114
|
Bastolla U. Computing protein dynamics from protein structure with elastic network models. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1186] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Ugo Bastolla
- Centro de Biologa Molecular Severo Ochoa (CSIC‐UAM)Universidad Autónoma de MadridMadridSpain
| |
Collapse
|
115
|
C(α) torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation. J Mol Model 2014; 20:2196. [PMID: 24728650 DOI: 10.1007/s00894-014-2196-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2014] [Accepted: 03/02/2014] [Indexed: 02/02/2023]
Abstract
Given the increasing complexity of simulated molecular systems, and the fact that simulation times have now reached milliseconds to seconds, immense amounts of data (in the gigabyte to terabyte range) are produced in current molecular dynamics simulations. Manual analysis of these data is a very time-consuming task, and important events that lead from one intermediate structure to another can become occluded in the noise resulting from random thermal fluctuations. To overcome these problems and facilitate a semi-automated data analysis, we introduce in this work a measure based on C(α) torsion angles: torsion angles formed by four consecutive C(α) atoms. This measure describes changes in the backbones of large systems on a residual length scale (i.e., a small number of residues at a time). Cluster analysis of individual C(α) torsion angles and its fuzzification led to continuous time patches representing (meta)stable conformations and to the identification of events acting as transitions between these conformations. The importance of a change in torsion angle to structural integrity is assessed by comparing this change to the average fluctuations in the same torsion angle over the complete simulation. Using this novel measure in combination with other measures such as the root mean square deviation (RMSD) and time series of distance measures, we performed an in-depth analysis of a simulation of the open form of DNA polymerase I. The times at which major conformational changes occur and the most important parts of the molecule and their interrelations were pinpointed in this analysis. The simultaneous determination of the time points and localizations of major events is a significant advantage of the new bottom-up approach presented here, as compared to many other (top-down) approaches in which only the similarity of the complete structure is analyzed.
Collapse
|
116
|
Abstract
Proteins are fascinating supramolecular structures, which are able to recognize ligands transforming binding information into chemical signals. They can transfer information across the cell, can catalyse complex chemical reactions, and are able to transform energy into work with much more efficiency than any human engine. The unique abilities of proteins are tightly coupled with their dynamic properties, which are coded in a complex way in the sequence and carefully refined by evolution. Despite its importance, our experimental knowledge of protein dynamics is still rather limited, and mostly derived from theoretical calculations. I will review here, in a systematic way, the current state-of-the-art theoretical approaches to the study of protein dynamics, emphasizing the most recent advances, examples of use and the expected lines of development in the near future.
Collapse
Affiliation(s)
- Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), Baldiri i Reixac 8, Barcelona 08028, Spain.
| |
Collapse
|
117
|
Zheng W, Tekpinar M. Analysis of protein conformational transitions using elastic network model. Methods Mol Biol 2014; 1084:159-72. [PMID: 24061921 DOI: 10.1007/978-1-62703-658-0_9] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/20/2023]
Abstract
In this chapter, we demonstrate the usage of a coarse-grained elastic network model to analyze protein conformational transitions in the NS3 helicase (NS3hel) of Hepatitis C virus (HCV). This analysis allows us to identify and visualize collective domain motions involved in the conformational transitions and predict the order of structural events during the transitions. It is highly efficient and applicable to many multi-domain protein structures which undergo large conformational changes to fulfill their functions. This method is made available through a Web server ( http://enm.lobos.nih.gov ).
Collapse
Affiliation(s)
- Wenjun Zheng
- Department of Physics, University at Buffalo, Buffalo, NY, USA
| | | |
Collapse
|
118
|
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
|
119
|
Abstract
By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| |
Collapse
|
120
|
Identifying essential pairwise interactions in elastic network model using the alpha shape theory. J Comput Chem 2014; 35:1111-21. [DOI: 10.1002/jcc.23587] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2013] [Revised: 01/24/2014] [Accepted: 02/26/2014] [Indexed: 11/07/2022]
|
121
|
Zheng W. All-atom and coarse-grained simulations of the forced unfolding pathways of the SNARE complex. Proteins 2014; 82:1376-86. [DOI: 10.1002/prot.24505] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2013] [Revised: 12/11/2013] [Accepted: 01/06/2014] [Indexed: 01/03/2023]
Affiliation(s)
- Wenjun Zheng
- Department of Physics; University at Buffalo, State University of New York; New York
| |
Collapse
|
122
|
Hills RD. Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics. Methods Mol Biol 2014; 1084:123-140. [PMID: 24061919 DOI: 10.1007/978-1-62703-658-0_7] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.
Collapse
Affiliation(s)
- Ronald D Hills
- Department of Pharmaceutical Sciences, University of New England, Portland, ME, USA
| |
Collapse
|
123
|
Abstract
Morphing was initially developed as a cinematic effect, where one image is seamlessly transformed into another image. The technique was widely adopted by biologists to visualize the transition between protein conformational states, generating an interpolated pathway from an initial to a final protein structure. Geometric morphing seeks to create visually suggestive movies that illustrate structural changes between conformations but do not necessarily represent a biologically relevant pathway, while minimum energy path (MEP) interpolations aim at describing the true transition state between the crystal structure minima in the energy landscape.
Collapse
Affiliation(s)
- Dahlia R Weiss
- Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, CA, USA
| | | |
Collapse
|
124
|
|
125
|
Huang MJ, Kapral R, Mikhailov AS, Chen HY. Coarse-grain simulations of active molecular machines in lipid bilayers. J Chem Phys 2013; 138:195101. [PMID: 23697442 DOI: 10.1063/1.4803507] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
A coarse-grain method for simulations of the dynamics of active protein inclusions in lipid bilayers is described. It combines the previously proposed hybrid simulations of bilayers [M.-J. Huang, R. Kapral, A. S. Mikhailov, and H.-Y. Chen, J. Chem. Phys. 137, 055101 (2012)], based on molecular dynamics for the lipids and multi-particle collision dynamics for the solvent, with an elastic-network description of active proteins. The method is implemented for a model molecular machine which performs active conformational motions induced by ligand binding and its release after reaction. The situation characteristic for peripheral membrane proteins is considered. Statistical investigations of the effects of single active or passive inclusions on the shape of the membrane are carried out. The results show that the peripheral machine produces asymmetric perturbations of the thickness of two leaflets of the membrane. It also produces a local saddle in the midplane height of the bilayer. Analysis of the power spectrum of the fluctuations of the membrane midplane shows that the conformational motion of the machine perturbs these membrane fluctuations. The hydrodynamic lipid flows induced by cyclic conformational changes in the machine are analyzed. It is shown that such flows are long-ranged and should provide an additional important mechanism for interactions between active inclusions in biological membranes.
Collapse
Affiliation(s)
- Mu-Jie Huang
- Department of Physics, National Central University, Jhongli 32001, Taiwan
| | | | | | | |
Collapse
|
126
|
Al-Bluwi I, Vaisset M, Siméon T, Cortés J. Modeling protein conformational transitions by a combination of coarse-grained normal mode analysis and robotics-inspired methods. BMC STRUCTURAL BIOLOGY 2013; 13 Suppl 1:S2. [PMID: 24564964 PMCID: PMC3953241 DOI: 10.1186/1472-6807-13-s1-s2] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
BACKGROUND Obtaining atomic-scale information about large-amplitude conformational transitions in proteins is a challenging problem for both experimental and computational methods. Such information is, however, important for understanding the mechanisms of interaction of many proteins. METHODS This paper presents a computationally efficient approach, combining methods originating from robotics and computational biophysics, to model protein conformational transitions. The ability of normal mode analysis to predict directions of collective, large-amplitude motions is applied to bias the conformational exploration performed by a motion planning algorithm. To reduce the dimension of the problem, normal modes are computed for a coarse-grained elastic network model built on short fragments of three residues. Nevertheless, the validity of intermediate conformations is checked using the all-atom model, which is accurately reconstructed from the coarse-grained one using closed-form inverse kinematics. RESULTS Tests on a set of ten proteins demonstrate the ability of the method to model conformational transitions of proteins within a few hours of computing time on a single processor. These results also show that the computing time scales linearly with the protein size, independently of the protein topology. Further experiments on adenylate kinase show that main features of the transition between the open and closed conformations of this protein are well captured in the computed path. CONCLUSIONS The proposed method enables the simulation of large-amplitude conformational transitions in proteins using very few computational resources. The resulting paths are a first approximation that can directly provide important information on the molecular mechanisms involved in the conformational transition. This approximation can be subsequently refined and analyzed using state-of-the-art energy models and molecular modeling methods.
Collapse
Affiliation(s)
- Ibrahim Al-Bluwi
- CNRS, LAAS, 7 avenue du colonel Roche, F-31400 Toulouse, France
- Univ de Toulouse, LAAS, F-31400 Toulouse, France
| | - Marc Vaisset
- CNRS, LAAS, 7 avenue du colonel Roche, F-31400 Toulouse, France
- Univ de Toulouse, LAAS, F-31400 Toulouse, France
| | - Thierry Siméon
- CNRS, LAAS, 7 avenue du colonel Roche, F-31400 Toulouse, France
- Univ de Toulouse, LAAS, F-31400 Toulouse, France
| | - Juan Cortés
- CNRS, LAAS, 7 avenue du colonel Roche, F-31400 Toulouse, France
- Univ de Toulouse, LAAS, F-31400 Toulouse, France
| |
Collapse
|
127
|
Hodak H. The Nobel Prize in chemistry 2013 for the development of multiscale models of complex chemical systems: a tribute to Martin Karplus, Michael Levitt and Arieh Warshel. J Mol Biol 2013; 426:1-3. [PMID: 24184197 DOI: 10.1016/j.jmb.2013.10.037] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Hélène Hodak
- Journal of Molecular Biology-Elsevier, 600 Technology Square, Cambridge, MA 02139, USA.
| |
Collapse
|
128
|
Tripathi S, Portman JJ. Allostery and Folding of the N-terminal Receiver Domain of Protein NtrC. J Phys Chem B 2013; 117:13182-93. [DOI: 10.1021/jp403181p] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Swarnendu Tripathi
- Department
of Physics, University of Houston, Houston, Texas 77204, United States
| | - John J. Portman
- Department
of Physics, Kent State University, Kent, Ohio 44242, United States
| |
Collapse
|
129
|
|
130
|
NMR spectroscopy on domain dynamics in biomacromolecules. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2013; 112:58-117. [DOI: 10.1016/j.pbiomolbio.2013.05.001] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2013] [Revised: 05/06/2013] [Accepted: 05/07/2013] [Indexed: 12/22/2022]
|
131
|
Xia F, Tong D, Lu L. Robust Heterogeneous Anisotropic Elastic Network Model Precisely Reproduces the Experimental B-factors of Biomolecules. J Chem Theory Comput 2013; 9:3704-14. [PMID: 26584122 DOI: 10.1021/ct4002575] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A computational method called the progressive fluctuation matching (PFM) is developed for constructing robust heterogeneous anisotropic network models (HANMs) for biomolecular systems. An HANM derived through the PFM approach consists of harmonic springs with realistic positive force constants, and yields the calculated B-factors that are basically identical to the experimental ones. For the four tested protein systems including crambin, trypsin inhibitor, HIV-1 protease, and lysozyme, the root-mean-square deviations between the experimental and the computed B-factors are only 0.060, 0.095, 0.247, and 0.049 Å(2), respectively, and the correlation coefficients are 0.99 for all. By comparing the HANM/ANM normal modes to their counterparts derived from both an atomistic force field and an NMR structure ensemble, it is found that HANM may provide more accurate results on protein dynamics.
Collapse
Affiliation(s)
- Fei Xia
- School of Biological Sciences, Nanyang Technological University , 60 Nanyang Drive, Singapore, 637551
| | - Dudu Tong
- School of Biological Sciences, Nanyang Technological University , 60 Nanyang Drive, Singapore, 637551
| | - Lanyuan Lu
- School of Biological Sciences, Nanyang Technological University , 60 Nanyang Drive, Singapore, 637551
| |
Collapse
|
132
|
Ping J, Hao P, Li YX, Wang JF. Molecular dynamics studies on the conformational transitions of adenylate kinase: a computational evidence for the conformational selection mechanism. BIOMED RESEARCH INTERNATIONAL 2013; 2013:628536. [PMID: 23936827 PMCID: PMC3712241 DOI: 10.1155/2013/628536] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/05/2013] [Accepted: 06/13/2013] [Indexed: 12/22/2022]
Abstract
Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme that catalyzes reversible transfer of phosphoryl group from ATP to AMP with a large-scale domain motion. The detailed mechanism for this conformational transition remains unknown. In the current study, we performed long time-scale molecular dynamics simulations on both open and closed states of ADK. Based on the structural analyses of the simulation trajectories, we detected over 20 times conformational transitions between the open and closed states of ADK and identified two novel conformations as intermediate states in the catalytic processes. With these findings, we proposed a possible mechanism for the large-scale domain motion of Escherichia coli ADK and its catalytic process: (1) the substrate free ADK adopted an open conformation; (2) ATP bound with LID domain closure; (3) AMP bound with NMP domain closure; (4) phosphoryl transfer occurred with ATP, and AMP converted into two ADPs, and no conformational transition was detected in the enzyme; (5) LID domain opened with one ADP released; (6) another ADP released with NMP domain open. As both open and closed states sampled a wide range of conformation transitions, our simulation strongly supported the conformational selection mechanism for Escherichia coli ADK.
Collapse
Affiliation(s)
- Jie Ping
- Pathogen Diagnostic Center, Institut Pasteur of Shanghai Chinese Academy of Sciences, Shanghai 200025, China
| | - Pei Hao
- Pathogen Diagnostic Center, Institut Pasteur of Shanghai Chinese Academy of Sciences, Shanghai 200025, China
- Shanghai Center for Bioinformation Technology, 100 Qinzhou Road, Shanghai 200235, China
| | - Yi-Xue Li
- Shanghai Center for Bioinformation Technology, 100 Qinzhou Road, Shanghai 200235, China
- Bioinformatics Center, Key Laboratory of Systems Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China
| | - Jing-Fang Wang
- Shanghai Center for Bioinformation Technology, 100 Qinzhou Road, Shanghai 200235, China
- Key Laboratory of Systems Biomedicine (Ministry of Education), Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
133
|
Terada TP, Kimura T, Sasai M. Entropic mechanism of allosteric communication in conformational transitions of dihydrofolate reductase. J Phys Chem B 2013; 117:12864-77. [PMID: 23705773 DOI: 10.1021/jp402071m] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The mechanism of allosteric conformational transitions of Escherichia coli dihydrofolate reductase (DHFR) is investigated theoretically by applying a newly developed coarse-grained model. Functional forms of interaction potentials in the model depend on the local structural environments around those interactions to represent the many-residue effects due to atomic packing in each local region, and hence, this model is called "the chameleon model". The chameleon model consistently describes the free-energy landscape of two conformational transitions in the catalytic cycle of DHFR, which we call conformational transition 1 (CT1) and conformational transition 2 (CT2); CT1 is accompanied by the hydride transfer reaction, and CT2 is accompanied by the product ligand release. The transition state of CT1 is entropically stabilized by the disordering of loops at the peripheral regions of the protein, which enhances the positively correlated fluctuations at the center part of the protein, showing that the allosteric communication between distant regions through the central region is intrinsically associated with the entropic stabilization of the transition state. The transition state of CT2 is entropically stabilized through the mechanism that enhances the breathing motion of two domains, showing that the difference in the distribution of interactions brings about the difference in the transition mechanism between CT1 and CT2. The chameleon model opens a way to consistently describe the dynamical energy landscape of enzymatic reactions.
Collapse
Affiliation(s)
- Tomoki P Terada
- Department of Computational Science and Engineering and ‡Department of Applied Physics, Nagoya University , Nagoya 464-8603, Japan
| | | | | |
Collapse
|
134
|
Sfriso P, Hospital A, Emperador A, Orozco M. Exploration of conformational transition pathways from coarse-grained simulations. ACTA ACUST UNITED AC 2013; 29:1980-6. [PMID: 23740746 DOI: 10.1093/bioinformatics/btt324] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
MOTIVATION A new algorithm to trace conformational transitions in proteins is presented. The method uses discrete molecular dynamics as engine to sample protein conformational space. A multiple minima Go-like potential energy function is used in combination with several enhancing sampling strategies, such as metadynamics, Maxwell Demon molecular dynamics and essential dynamics. The method, which shows an unprecedented computational efficiency, is able to trace a wide range of known experimental transitions. Contrary to simpler methods our strategy does not introduce distortions in the chemical structure of the protein and is able to reproduce well complex non-linear conformational transitions. The method, called GOdMD, can easily introduce additional restraints to the transition (presence of ligand, known intermediate, known maintained contacts, …) and is freely distributed to the community through the Spanish National Bioinformatics Institute (http://mmb.irbbarcelona.org/GOdMD). AVAILABILITY Freely available on the web at http://mmb.irbbarcelona.org/GOdMD.
Collapse
Affiliation(s)
- Pedro Sfriso
- Institute for Research in Biomedicine (IRB Barcelona), Joint IRB-BSC Program in Computational Biology, Baldiri Reixac 10, Barcelona, Spain
| | | | | | | |
Collapse
|
135
|
Affiliation(s)
- Marissa G. Saunders
- Department of Chemistry, Institute for Biophysical Dynamics, James Franck Institute, and Computation Institute, University of Chicago, Chicago, Illinois 60637;
| | - Gregory A. Voth
- Department of Chemistry, Institute for Biophysical Dynamics, James Franck Institute, and Computation Institute, University of Chicago, Chicago, Illinois 60637;
| |
Collapse
|
136
|
Sahún-Roncero M, Rubio-Ruiz B, Saladino G, Conejo-García A, Espinosa A, Velázquez-Campoy A, Gervasio FL, Entrena A, Hurtado-Guerrero R. The mechanism of allosteric coupling in choline kinase α1 revealed by the action of a rationally designed inhibitor. Angew Chem Int Ed Engl 2013; 52:4582-6. [PMID: 23441033 DOI: 10.1002/anie.201209660] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2012] [Revised: 01/16/2013] [Indexed: 01/21/2023]
Abstract
Applying a CHOK hold: Combined experimental and computational studies of the binding mode of a rationally designed inhibitor of the dimeric choline kinase α1 (CHOKα1) explain the molecular mechanism of negative cooperativity (see scheme) and how the monomers are connected. The results give insight into how the symmetry of the dimer can be partially conserved despite a lack of conservation in the static crystal structures.
Collapse
Affiliation(s)
- María Sahún-Roncero
- Institute of Biocomputation and Physics of Complex Systems (BIFI), University of Zaragoza, BIFI-IQFR (CSIC) Joint Unit, Mariano Esquillor s/n, Campus Rio Ebro, Edificio I+D; Fundacion ARAID, Edificio Pignatelli 36, Spain
| | | | | | | | | | | | | | | | | |
Collapse
|
137
|
Sahún-Roncero M, Rubio-Ruiz B, Saladino G, Conejo-García A, Espinosa A, Velázquez-Campoy A, Gervasio FL, Entrena A, Hurtado-Guerrero R. The Mechanism of Allosteric Coupling in Choline Kinase α1 Revealed by the Action of a Rationally Designed Inhibitor. Angew Chem Int Ed Engl 2013. [DOI: 10.1002/ange.201209660] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
138
|
Mohazab AR, Plotkin SS. Polymer uncrossing and knotting in protein folding, and their role in minimal folding pathways. PLoS One 2013; 8:e53642. [PMID: 23365638 PMCID: PMC3554774 DOI: 10.1371/journal.pone.0053642] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2012] [Accepted: 11/30/2012] [Indexed: 11/19/2022] Open
Abstract
We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold [Formula: see text], [Formula: see text], [Formula: see text], and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation "alignment". The consensus minimal pathway is constructed and shown schematically for representative cases of an [Formula: see text], [Formula: see text], and knotted protein. An overlap parameter is defined between pathways; we find that [Formula: see text] proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and [Formula: see text] proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.
Collapse
Affiliation(s)
- Ali R. Mohazab
- Department of Physics and Astronomy, University of British Columbia, Vancouver, B.C, Canada
| | - Steven S. Plotkin
- Department of Physics and Astronomy, University of British Columbia, Vancouver, B.C, Canada
| |
Collapse
|
139
|
Elber R, Kirmizialtin S. Molecular machines. Curr Opin Struct Biol 2013; 23:206-11. [PMID: 23305848 DOI: 10.1016/j.sbi.2012.12.002] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2012] [Accepted: 12/03/2012] [Indexed: 12/12/2022]
Abstract
Molecular machines (MM) are essential components of living cells. They conduct mechanical work, transport materials into and out of cells, assist in processing enzymatic reactions, and more. Their operations are frequently combined with significant conformational transitions. Computational studies of these conformational transitions and their coupling to molecular functions are discussed. It is argued that coarse descriptions of these molecules which are based on mass density and shape provide useful information on directions of action. It is further argued that MM are likely to have well focused and narrow reaction pathways. The proposal for such pathways is supported by evolutionary analyses of homologous machines. Finally, these observations are used to build atomically detailed models of these systems that are making the link from structure to functions (kinetics and thermodynamics). For that purpose enhanced sampling techniques are required.
Collapse
Affiliation(s)
- Ron Elber
- Department of Chemistry and Biochemistry, University of Texas at Austin, 105 East 24th St., Stop A5300 Austin, TX 78712-0165, USA.
| | | |
Collapse
|
140
|
Daily MD, Yu H, Phillips GN, Cui Q. Allosteric activation transitions in enzymes and biomolecular motors: insights from atomistic and coarse-grained simulations. Top Curr Chem (Cham) 2013; 337:139-64. [PMID: 23468286 PMCID: PMC3976962 DOI: 10.1007/128_2012_409] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The chemical step in enzymes is usually preceded by a kinetically distinct activation step that involves large-scale conformational transitions. In "simple" enzymes this step corresponds to the closure of the active site; in more complex enzymes, such as biomolecular motors, the activation step is more complex and may involve interactions with other biomolecules. These activation transitions are essential to the function of enzymes and perturbations in the scale and/or rate of these transitions are implicated in various serious human diseases; incorporating key flexibilities into engineered enzymes is also considered a major remaining challenge in rational enzyme design. Therefore it is important to understand the underlying mechanism of these transitions. This is a significant challenge to both experimental and computational studies because of the allosteric and multi-scale nature of such transitions. Using our recent studies of two enzyme systems, myosin and adenylate kinase (AK), we discuss how atomistic and coarse-grained simulations can be used to provide insights into the mechanism of activation transitions in realistic systems. Collectively, the results suggest that although many allosteric transitions can be viewed as domain displacements mediated by flexible hinges, there are additional complexities and various deviations. For example, although our studies do not find any evidence for "cracking" in AK, our results do underline the contribution of intra-domain properties (e.g., dihedral flexibility) to the rate of the transition. The study of mechanochemical coupling in myosin highlights that local changes important to chemistry require stabilization from more extensive structural changes; in this sense, more global structural transitions are needed to activate the chemistry in the active site. These discussions further emphasize the importance of better understanding factors that control the degree of co-operativity for allosteric transitions, again hinting at the intimate connection between protein stability and functional flexibility. Finally, a number of topics of considerable future interest are briefly discussed.
Collapse
Affiliation(s)
- Michael D Daily
- Pacific Northwest National Laboratory, Richland, Washington, 99352, USA
| | | | | | | |
Collapse
|
141
|
Schrank TP, Wrabl JO, Hilser VJ. Conformational heterogeneity within the LID domain mediates substrate binding to Escherichia coli adenylate kinase: function follows fluctuations. Top Curr Chem (Cham) 2013; 337:95-121. [PMID: 23543318 DOI: 10.1007/128_2012_410] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Proteins exist as dynamic ensembles of molecules, implying that protein amino acid sequences evolved to code for both the ground state structure as well as the entire energy landscape of excited states. Accumulating theoretical and experimental evidence suggests that enzymes use such conformational fluctuations to facilitate allosteric processes important for substrate binding and possibly catalysis. This phenomenon can be clearly demonstrated in Escherichia coli adenylate kinase, where experimentally observed local unfolding of the LID subdomain, as opposed to a more commonly postulated rigid-body opening motion, is related to substrate binding. Because "entropy promoting" glycine mutations designed to increase specifically the local unfolding of the LID domain also affect substrate binding, changes in the excited energy landscape effectively tune the function of this enzyme without changing the ground state structure or the catalytic site. Thus, additional thermodynamic information, above and beyond the single folded structure of an enzyme-substrate complex, is likely required for a full and quantitative understanding of how enzymes work.
Collapse
Affiliation(s)
- Travis P Schrank
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch at Galveston, 301 University Boulevard, Galveston, TX, 77555-1068, USA,
| | | | | |
Collapse
|
142
|
Wang Y, Gan L, Wang E, Wang J. Exploring the Dynamic Functional Landscape of Adenylate Kinase Modulated by Substrates. J Chem Theory Comput 2012; 9:84-95. [PMID: 26589012 DOI: 10.1021/ct300720s] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Adenylate kinase (ADK) has been explored widely, through both experimental and theoretical studies. However, still less is known about how the functional dynamics of ADK is modulated explicitly by its natural substrates. Here, we report a quantitative study of the dynamic energy landscape for ADK responding to the substrate binding by integrating both experimental investigations and theoretical modeling. We make theoretical predictions which are in remarkable agreement with the single molecule experiments on the substrate-bound complex. With our combined models of ADK in its apo form, in the presence of AMP or ATP, and in complex with both substrates, we specifically address the following key questions: (1) Are there intermediate state(s) during their catalytic cycle and if so how many? (2) How many pathways are there along the open-to-closed transitions and what are their corresponding weights? (3) How do substrates influence the pathway weights and the stability of the intermediates? (4) Which lid's motion is rate-limiting along the turnover cycle, the NMP or the LID domain? Our models predict two major parallel stepwise pathways and two on-pathway intermediates which are denoted as IN (NMP domain open while LID domain closed) and IL (LID domain open and NMP domain closed), respectively. Further investigation of temperature effects suggests that the IN pathway is dominant at room temperature, but the IL pathway is dominant at the optimal temperature. This leads us to propose that the IL pathway is more dominant by entropy and IN pathway by enthalpy. Remarkably, our results show that even with maximum concentrations of natural substrates, ADK still fluctuates between multiple functional states, reflecting an intrinsic capability of large-scale conformational fluctuations which may be essential to its biological function. The results based on the dual-ligands model provide the theoretical validation of random bisubstrate biproducts (Bi-Bi) mechanism for the enzymatic reaction of ADK. Additionally, the pathway flux analysis strongly suggests that the motion of the NMP domain is the rate-determining step for the conformational cycle (opening and closing).
Collapse
Affiliation(s)
- Yong Wang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, P.R. China
| | - Linfeng Gan
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, P.R. China
| | - Erkang Wang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, P.R. China
| | - Jin Wang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, P.R. China.,College of Physics, Jilin University, Changchun, Jilin, P.R. China.,Department of Chemistry and Physics, State University of New York at Stony Brook, Stony Brook, New York 11794-3400, United States
| |
Collapse
|
143
|
Martin DR, Matyushov DV. Solvated dissipative electro-elastic network model of hydrated proteins. J Chem Phys 2012; 137:165101. [DOI: 10.1063/1.4759105] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
|
144
|
Angarica VE, Sancho J. Protein dynamics governed by interfaces of high polarity and low packing density. PLoS One 2012; 7:e48212. [PMID: 23110216 PMCID: PMC3482218 DOI: 10.1371/journal.pone.0048212] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2012] [Accepted: 09/24/2012] [Indexed: 11/18/2022] Open
Abstract
The folding pathway, three-dimensional structure and intrinsic dynamics of proteins are governed by their amino acid sequences. Internal protein surfaces with physicochemical properties appropriate to modulate conformational fluctuations could play important roles in folding and dynamics. We show here that proteins contain buried interfaces of high polarity and low packing density, coined as LIPs: Light Interfaces of high Polarity, whose physicochemical properties make them unstable. The structures of well-characterized equilibrium and kinetic folding intermediates indicate that the LIPs of the corresponding native proteins fold late and are involved in local unfolding events. Importantly, LIPs can be identified using very fast and uncomplicated computational analysis of protein three-dimensional structures, which provides an easy way to delineate the protein segments involved in dynamics. Since LIPs can be retained while the sequences of the interacting segments diverge significantly, proteins could in principle evolve new functional features reusing pre-existing encoded dynamics. Large-scale identification of LIPS may contribute to understanding evolutionary constraints of proteins and the way protein intrinsic dynamics are encoded.
Collapse
Affiliation(s)
- Vladimir Espinosa Angarica
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Universidad de Zaragoza, Zaragoza, Spain
- Biocomputation and Complex Systems Physics Institute (BIFI), Joint Unit BIFI-IQFR, CSIC, Universidad de Zaragoza, Zaragoza, Spain
| | - Javier Sancho
- Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Universidad de Zaragoza, Zaragoza, Spain
- Biocomputation and Complex Systems Physics Institute (BIFI), Joint Unit BIFI-IQFR, CSIC, Universidad de Zaragoza, Zaragoza, Spain
- * E-mail:
| |
Collapse
|
145
|
Tekpinar M, Zheng W. Coarse-grained and all-atom modeling of structural states and transitions in hemoglobin. Proteins 2012; 81:240-52. [PMID: 22987685 DOI: 10.1002/prot.24180] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2012] [Revised: 08/27/2012] [Accepted: 09/10/2012] [Indexed: 11/08/2022]
Abstract
Hemoglobin (Hb), an oxygen-binding protein composed of four subunits (α1, α2, β1, and β2), is a well-known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse-grained (CG) modeling, all-atom simulation, and structural data from X-ray crystallography and wide-angle X-ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot-spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all-atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all-atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling.
Collapse
Affiliation(s)
- Mustafa Tekpinar
- Physics Department, University at Buffalo, Buffalo, New York 14260, USA
| | | |
Collapse
|
146
|
Düttmann M, Mittnenzweig M, Togashi Y, Yanagida T, Mikhailov AS. Complex intramolecular mechanics of G-actin--an elastic network study. PLoS One 2012; 7:e45859. [PMID: 23077498 PMCID: PMC3471905 DOI: 10.1371/journal.pone.0045859] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2012] [Accepted: 08/17/2012] [Indexed: 11/30/2022] Open
Abstract
Systematic numerical investigations of conformational motions in single actin molecules were performed by employing a simple elastic-network (EN) model of this protein. Similar to previous investigations for myosin, we found that G-actin essentially behaves as a strain sensor, responding by well-defined domain motions to mechanical perturbations. Several sensitive residues within the nucleotide-binding pocket (NBP) could be identified, such that the perturbation of any of them can induce characteristic flattening of actin molecules and closing of the cleft between their two mobile domains. Extending the EN model by introduction of a set of breakable links which become effective only when two domains approach one another, it was observed that G-actin can possess a metastable state corresponding to a closed conformation and that a transition to this state can be induced by appropriate perturbations in the NBP region. The ligands were roughly modeled as a single particle (ADP) or a dimer (ATP), which were placed inside the NBP and connected by elastic links to the neighbors. Our approximate analysis suggests that, when ATP is present, it stabilizes the closed conformation of actin. This may play an important role in the explanation why, in the presence of ATP, the polymerization process is highly accelerated.
Collapse
Affiliation(s)
- Markus Düttmann
- Department of Physical Chemistry, Fritz Haber Institute of the Max Planck Society, Berlin, Germany.
| | | | | | | | | |
Collapse
|
147
|
Abstract
A multiscale coarse-graining method called the normal-mode analysis based fluctuation matching (NMA-FM) is developed for constructing coarse-grained models of biomolecular systems. In the framework of normal-mode analysis, an arbitrary fine-grained model can be systematically converted to a more coarse-grained model, while the crucial low-frequency motions of the fine-grained system are able to be reproduced in the coarse-grained modeling. The method relies on the technique of fluctuation matching that has been devised earlier for parametrizing heterogeneous elastic network models based on data from atomistic molecular dynamics simulations. The new approach is quite efficient since it avoids expensive atomistic molecular dynamics simulations and can start from already coarse-grained elastic network models. In the practical aspect, the method is suitable for conformational analyses of large biomacromolecules and calculations of mechanical properties of biomaterials, which is demonstrated by the studied systems including an amyloid dimer, lysozyme and adenylate kinase proteins, and the S2 subdomain of myosin.
Collapse
Affiliation(s)
- Fei Xia
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore
| | - Lanyuan Lu
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, Singapore
| |
Collapse
|
148
|
Dans PD, Pérez A, Faustino I, Lavery R, Orozco M. Exploring polymorphisms in B-DNA helical conformations. Nucleic Acids Res 2012; 40:10668-78. [PMID: 23012264 PMCID: PMC3510489 DOI: 10.1093/nar/gks884] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The traditional mesoscopic paradigm represents DNA as a series of base-pair steps whose energy response to equilibrium perturbations is elastic, with harmonic oscillations (defining local stiffness) around a single equilibrium conformation. In addition, base sequence effects are often analysed as a succession of independent XpY base-pair steps (i.e. a nearest-neighbour (NN) model with only 10 unique cases). Unfortunately, recent massive simulations carried out by the ABC consortium suggest that the real picture of DNA flexibility may be much more complex. The paradigm of DNA flexibility therefore needs to be revisited. In this article, we explore in detail one of the most obvious violations of the elastic NN model of flexibility: the bimodal distributions of some helical parameters. We perform here an in-depth statistical analysis of a very large set of MD trajectories and also of experimental structures, which lead to very solid evidence of bimodality. We then suggest ways to improve mesoscopic models to account for this deviation from the elastic regime.
Collapse
Affiliation(s)
- Pablo D Dans
- Joint IRB-BSC Program on Computational Biology, Institute for Research in Biomedicine, Parc Cientific de Barcelona, Josep Samitier 1-5, Barcelona 08028, Spain
| | | | | | | | | |
Collapse
|
149
|
Schofield J, Inder P, Kapral R. Modeling of solvent flow effects in enzyme catalysis under physiological conditions. J Chem Phys 2012; 136:205101. [PMID: 22667589 DOI: 10.1063/1.4719539] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
A stochastic model for the dynamics of enzymatic catalysis in explicit, effective solvents under physiological conditions is presented. Analytically-computed first passage time densities of a diffusing particle in a spherical shell with absorbing boundaries are combined with densities obtained from explicit simulation to obtain the overall probability density for the total reaction cycle time of the enzymatic system. The method is used to investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The direct simulation of the enzyme-substrate binding and reaction is carried out using an elastic network model for the protein, and the solvent motions are described by multiparticle collision dynamics which incorporates hydrodynamic flow effects. Systems where solvent-enzyme coupling occurs through explicit intermolecular interactions, as well as systems where this coupling is taken into account by including the protein and substrate in the multiparticle collision step, are investigated and compared with simulations where hydrodynamic coupling is absent. It is demonstrated that the flow of solvent particles around the enzyme facilitates the large-scale hinge motion of the enzyme with bound substrates, and has a significant impact on the shape of the probability densities and average time scales of substrate binding for substrates near the enzyme, the closure of the enzyme after binding, and the overall time of completion of the cycle.
Collapse
Affiliation(s)
- Jeremy Schofield
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada.
| | | | | |
Collapse
|
150
|
Gur M, Erman B. Quasi-harmonic fluctuations of two bound peptides. Proteins 2012; 80:2769-79. [DOI: 10.1002/prot.24160] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2012] [Revised: 07/27/2012] [Accepted: 08/06/2012] [Indexed: 11/10/2022]
|