Energy Landscape in Protein Folding and Unfolding

Total Page:16

File Type:pdf, Size:1020Kb

Energy Landscape in Protein Folding and Unfolding Energy landscape in protein folding and unfolding Francesco Mallamacea,b,c,1, Carmelo Corsaroa,d, Domenico Mallamacee, Sebastiano Vasid, Cirino Vasia, Piero Baglionif, Sergey V. Buldyrevg, Sow-Hsin Chenb, and H. Eugene Stanleyc,1 aCNR-Istituto per i Processi Chimico Fisici Messina, I-98166 Messina, Italy; bDepartment of Nuclear Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139; cCenter for Polymer Studies and Department of Physics, Boston University, Boston, MA 02215; dDipartimento di Fisica e di Scienze della Terra, Università di Messina, I-98166 Messina, Italy; eConsorzio per lo Sviluppo dei Sistemi a Grande Interfase, Unità di Catania, I-95125 Catania, Italy; fDipartimento di Chimica, Università di Firenze and Consorzio per lo Sviluppo dei Sistemi a Grande Interfase, I-50019 Florence, Italy; and gDepartment of Physics, Yeshiva University, New York, NY 10033 Contributed by H. Eugene Stanley, December 22, 2015 (sent for review March 17, 2015; reviewed by Anders Nilsson and Michele Parrinello) We use 1H NMR to probe the energy landscape in the protein fold- simulation (9, 10) and different experimental techniques (11–16). ing and unfolding process. Using the scheme ⇄ reversible unfolded Whereas MD simulations directly model peptide conformational (intermediate) → irreversible unfolded (denatured) state, we study transitions in terms of the energy landscape, experiments supply the thermal denaturation of hydrated lysozyme that occurs when useful but limited information, revealing some details in the the temperature is increased. Using thermal cycles in the range structure and the collective dynamics of protein, both dry and in 295 < T < 365 K and following different trajectories along the protein solution. Examples of this include the spectroscopic techniques energy surface, we observe that the hydrophilic (the amide NH) and (e.g., NMR, neutron, X-ray, Raman, and FTIR) that supply data hydrophobic (methyl CH3 and methine CH) peptide groups evolve on structure and dynamical modes of the protein, and the and exhibit different behaviors. We also discuss the role of water calorimetric measurements that follow the reversible folding– and hydrogen bonding in the protein configurational stability. unfolding as far as irreversible denaturation (9, 10). The protein dynamics have been studied from the glass state in the deep protein folding | proton NMR | energy landscape | hydration water supercooled regime (T < 200 K) to the completely denatured state (T ’ 350 K) by proving the essential role of the hydration n intriguing problem of statistical physics concerns the water. As for bulk water, hydrogen bond (HB) interactions Aevolutionary pathways that molecular systems follow as they strongly determine the properties of these systems, for which SCIENCES form mesoscale structures and exhibit new functional behaviors water is not simply a solvent but is also an integral and active APPLIED PHYSICAL (1). An example of this problem is the self-organization of bio- component, i.e., it is itself an important “biomolecule” that plays systems that evolve from basic molecules. This challenging sub- both a dynamic and structural role (17). Hence, HB interactions ject is studied by using a variety of theoretical methods (2–4). are the key to understanding water’s properties and how water The free-energy landscape model is nowadays the most used to functions in biological environments (18). describe such phenomena and especially the aging of the protein Among the different experimental methods, calorimetry, by folding mechanism (1, 5, 6), i.e., the way in which proteins fold to monitoring the process reaction rates, focuses directly the energetic their native state and then unfold (protein denaturation) (6, 7). properties of the hydrated proteins folding–unfolding mechanism BIOPHYSICS AND The model is based on the idea that in complex materials and and measures the enthalpy and entropy behaviors. Unlike common COMPUTATIONAL BIOLOGY systems there are many thermodynamical configurations in which chemical reaction rates, which increase as the temperature is in- the free-energy surface exhibits a number of local minima sepa- creased, the rate constant of the protein folding reaction initially rated by barriers, i.e., as the system explores its phase space the increases on increasing T, by following an Arrhenius law goes to a trajectory of its evolution is an alternating sequence of local maximum, and then decreases as the temperature continues to energy minima and saddle points (transition states), which are increase (19). This latter situation is also well-described in MD associated with the positions of all of the system particles. A simulations (20). The activation enthalpy is thus T-dependent and trajectory thus specifies the path of the system as it evolves by the corresponding protein specific heat, CPðTÞ, exhibits large moving across its energy landscape. A peptide is a linear chain of amino acids, and globular pro- Significance teins are polypeptide chains that fold into their native confor- mation. During the folding process a polypeptide undergoes many conformational changes and there is a significant decrease Protein folding represents an open question in science, and the in the system configurational entropy as the native state is free-energy landscape framework is one way to describe it. In particular, the role played by water in the processes is of spe- approached. To understand folding we focus on how proteins – search conformational space. The process is accompanied by cial interest. To clarify these issues we study, during folding many microscopic reactions, the nature of which is determined unfolding, the temperature evolution of the magnetization for by the specifics of the energy surface. Thus, the characteristics of hydrophilic and hydrophobic groups of hydrated lysozyme the energy surface of a polypeptide chain are the key to a quan- using NMR spectroscopy. Our findings confirm the validity of titative understanding of folding. Although the degrees of freedom the theoretical scenario of a process dominated by different energetic routes, also explaining the water role in the protein of a polypeptide chain allow an enormously large number of configuration stability. We also highlight that the protein native possible configurations, “constraints” on the energy decrease these state limit is represented by the water singular temperature that configurations visited in the folding reaction to a limited number characterizes its compressibility and expansivity and is the origin (8). Understanding the free-energy surface (“landscape”)enables of the thermodynamical anomalies of its liquid state. us to understand the folding process. A balance between the potential energy and the configurational entropy leads to a free- Author contributions: F.M., C.C., P.B., S.V.B., and H.E.S. designed research; F.M., C.C., D.M., energy barrier that generates the two-state folding behavior usu- S.V., and C.V. performed research; F.M., D.M., S.V., P.B., S.V.B., and S.-H.C. analyzed data; ally observed in small proteins. The potential energy decreases as and F.M., C.C., and H.E.S. wrote the paper. the native state is approached and favors folding, but decreasing Reviewers: A.N., SLAC National Accelerator Laboratory; and M.P., ETH Zurich. the entropy of the configuration is unfavorable to folding. The authors declare no conflict of interest. The thermodynamics and kinetics of folding and unfolding 1To whom correspondence may be addressed. Email: [email protected] or have been intensively studied by molecular dynamics (MD) [email protected]. www.pnas.org/cgi/doi/10.1073/pnas.1524864113 PNAS | March 22, 2016 | vol. 113 | no. 12 | 3159–3163 Downloaded by guest on September 25, 2021 changes: first it increases with T, reaches a maximum, and then precisely, to quantify the behavior of the protein hydrophilic (the decreases. Hence a CPðTÞ plot against T exhibits an endothermic amide NH) and hydrophobic (methyl and methine) groups, during peak whose area is related to the enthalpy of transformation of the folding–unfolding process, we use as control variable the the protein. Such a behavior is also strongly influenced by the temperature, which slows the process kinetics to many hours. protein hydrophobic side chains and hydration water (19, 21, 22). The goal is to determine the topography of the protein energy NMR spectroscopy instead can be used to follow folding (or landscape by following different trajectories along the energy unfolding) processes by probing interactions as they form at surface. This allows us also to consider explicitly the role of water the level of individual residues, and to compare the findings and hydrogen bonding in protein configurational stability. obtained with simulations (13, 23–25). NMR studies involve in- vestigation of the equilibrium conversion between native and Results denatured states resulting from a protein being subjected to heat, Thermal denaturation of lysozyme occurs according to the extremes of pH, or chemical denaturants. The kinetics of the scheme, native state (N) ⇄ reversible unfolded (intermediate) process may limit the observations. If the folding occurs on a state (RU) → irreversible unfolded (denatured) state (IU); N ⇄ time scale for which spectra can be measured and recorded se- RU → IU (21, 22). This is consistent with the general view that quentially, the spectral changes accompanying the reaction can the first step in the denaturation of small one-domain globular be
Recommended publications
  • DNA Energy Landscapes Via Calorimetric Detection of Microstate Ensembles of Metastable Macrostates and Triplet Repeat Diseases
    DNA energy landscapes via calorimetric detection of microstate ensembles of metastable macrostates and triplet repeat diseases Jens Vo¨ lkera, Horst H. Klumpb, and Kenneth J. Breslauera,c,1 aDepartment of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, 610 Taylor Rd, Piscataway, NJ 08854; bDepartment of Molecular and Cell Biology, University of Cape Town, Private Bag, Rondebosch 7800, South Africa; and cCancer Institute of New Jersey, New Brunswick, NJ 08901 Communicated by I. M. Gelfand, Rutgers, The State University of New Jersey, Piscataway, NJ, October 15, 2008 (received for review September 8, 2008) Biopolymers exhibit rough energy landscapes, thereby allowing biological role(s), if any, of kinetically stable (metastable) mi- biological processes to access a broad range of kinetic and ther- crostates that make up the time-averaged, native state ensembles modynamic states. In contrast to proteins, the energy landscapes of macroscopic nucleic acid states remains to be determined. of nucleic acids have been the subject of relatively few experimen- As part of an effort to address this deficiency, we report here tal investigations. In this study, we use calorimetric and spectro- experimental evidence for the presence of discrete microstates scopic observables to detect, resolve, and selectively enrich ener- in metastable triplet repeat bulge looped ⍀-DNAs of potential getically discrete ensembles of microstates within metastable DNA biological significance. The specific triplet repeat bulge looped structures. Our results are consistent with metastable, ‘‘native’’ ⍀-DNA species investigated mimic slipped DNA structures DNA states being composed of an ensemble of discrete and corresponding to intermediates in the processes that lead to kinetically stable microstates of differential stabilities, rather than DNA expansion in triplet repeat diseases (36).
    [Show full text]
  • Ruggedness in the Free Energy Landscape Dictates Misfolding of the Prion Protein
    Article Ruggedness in the Free Energy Landscape Dictates Misfolding of the Prion Protein Roumita Moulick 1, Rama Reddy Goluguri 1 and Jayant B. Udgaonkar 1,2 1 - National Centre for Biological Sciences, Tata Institute of Fundamental Research, Bengaluru 560065, India 2 - Indian Institute of Science Education and Research, Pashan, Pune 411008, India Correspondence to Jayant B. Udgaonkar: Indian Institute of Science Education and Research, Pashan, Pune 411008, India. [email protected] https://doi.org/10.1016/j.jmb.2018.12.009 Edited by Sheena Radford Abstract Experimental determination of the key features of the free energy landscapes of proteins, which dictate their adeptness to fold correctly, or propensity to misfold and aggregate and which are modulated upon a change from physiological to aggregation-prone conditions, is a difficult challenge. In this study, sub-millisecond kinetic measurements of the folding and unfolding of the mouse prion protein reveal how the free energy landscape becomes more complex upon a shift from physiological (pH 7) to aggregation-prone (pH 4) conditions. Folding and unfolding utilize the same single pathway at pH 7, but at pH 4, folding occurs on a pathway distinct from the unfolding pathway. Moreover, the kinetics of both folding and unfolding at pH 4 depend not only on the final conditions but also on the conditions under which the processes are initiated. Unfolding can be made to switch to occur on the folding pathway by varying the initial conditions. Folding and unfolding pathways appear to occupy different regions of the free energy landscape, which are separated by large free energy barriers that change with a change in the initial conditions.
    [Show full text]
  • Introducing the Levinthal's Protein Folding Paradox and Its Solution
    Article pubs.acs.org/jchemeduc Introducing the Levinthal’s Protein Folding Paradox and Its Solution Leandro Martínez* Department of Physical Chemistry, Institute of Chemistry, University of Campinas, Campinas, Saõ Paulo 13083-970, Brazil *S Supporting Information ABSTRACT: The protein folding (Levinthal’s) paradox states that it would not be possible in a physically meaningful time to a protein to reach the native (functional) conformation by a random search of the enormously large number of possible structures. This paradox has been solved: it was shown that small biases toward the native conformation result in realistic folding times of realistic-length sequences. This solution of the paradox is, however, not amenable to most chemistry or biology students due to the demanding mathematics. Here, a simplification of the study of the paradox and its solution is provided so that it is accessible to chemists and biologists at an undergraduate or graduate level. Despite its simplicity, the model captures some fundamental aspects of the protein folding mechanism and allows students to grasp the actual significance of the popular folding funnel representation of the protein energy landscape. The analysis of the folding model provides a rich basis for a discussion of the relationships between kinetics and thermodynamics on a fundamental level and on student perception of the time-scales of molecular phenomena. KEYWORDS: Upper-Division Undergraduate, Graduate Education, Biochemistry, Biophysical Chemistry, Molecular Biology, Proteins/Peptides,
    [Show full text]
  • Chapter 6 Protein Structure and Folding
    Chapter 6 Protein Structure and Folding 1. Secondary Structure 2. Tertiary Structure 3. Quaternary Structure and Symmetry 4. Protein Stability 5. Protein Folding Myoglobin Introduction 1. Proteins were long thought to be colloids of random structure 2. 1934, crystal of pepsin in X-ray beam produces discrete diffraction pattern -> atoms are ordered 3. 1958 first X-ray structure solved, sperm whale myoglobin, no structural regularity observed 4. Today, approx 50’000 structures solved => remarkable degree of structural regularity observed Hierarchy of Structural Layers 1. Primary structure: amino acid sequence 2. Secondary structure: local arrangement of peptide backbone 3. Tertiary structure: three dimensional arrangement of all atoms, peptide backbone and amino acid side chains 4. Quaternary structure: spatial arrangement of subunits 1) Secondary Structure A) The planar peptide group limits polypeptide conformations The peptide group ha a rigid, planar structure as a consequence of resonance interactions that give the peptide bond ~40% double bond character The trans peptide group The peptide group assumes the trans conformation 8 kJ/mol mire stable than cis Except Pro, followed by cis in 10% Torsion angles between peptide groups describe polypeptide chain conformations The backbone is a chain of planar peptide groups The conformation of the backbone can be described by the torsion angles (dihedral angles, rotation angles) around the Cα-N (Φ) and the Cα-C bond (Ψ) Defined as 180° when extended (as shown) + = clockwise, seen from Cα Not
    [Show full text]
  • Protein Misfolding in Human Diseases
    Linköping Studies in Science and Technology Dissertation No. 1239 Protein Misfolding in Human Diseases Karin Almstedt Biochemistry Department of Physics, Chemistry and Biology Linköping University, SE-581 83 Linköping, Sweden Linköping 2009 The cover shows a Himalayan mountain range, symbolizing a protein folding landscape. During the course of the research underlying this thesis, Karin Almstedt was enrolled in Forum Scientium, a multidisciplinary doctoral program at Linköping University, Sweden. Copyright © 2009 Karin Almstedt ISBN: 978-91-7393-698-9 ISSN: 0345-7524 Printed in Sweden by LiU-Tryck Linköping 2009 Never Stop Exploring 83 ABSTRACT The studies in this thesis are focused on misfolded proteins involved in human disease. There are several well known diseases that are due to aberrant protein folding. These types of diseases can be divided into three main categories: 1. Loss-of-function diseases 2. Gain-of-toxic-function diseases 3. Infectious misfolding diseases Most loss-of-function diseases are caused by aberrant folding of important proteins. These proteins often misfold due to inherited mutations. The rare disease carbonic anhydrase II deficiency syndrome (CADS) can manifest in carriers of point mutations in the human carbonic anhydrase II (HCA II) gene. One mutation associated with CADS entails the His107Tyr (H107Y) substitution. We have demonstrated that the H107Y mutation is a remarkably destabilizing mutation influencing the folding behavior of HCA II. A mutational survey of position H107 and a neighboring conserved position E117 has been performed entailing the mutants H107A, H107F, H107N, E117A and the double mutants H107A/E117A and H107N/E117A. We have also shown that the binding of specific ligands can stabilize the disease causing mutant, and shift the folding equilibrium towards the native state, providing a starting point for small molecule drugs for CADS.
    [Show full text]
  • Complex Energy Landscapes” (IPAM Program, Fall 2017)
    Whitepaper: “Complex Energy Landscapes” (IPAM program, fall 2017) Authors list (alphabetically) Nestor F. Aguirre, Chris Anderson, Lorenzo Boninsegna, Gábor Csányi, Mauricio del Razo Sarmina, Marco Di Gennaro, Florent Hédin, Graeme Henkelman, Richard G. Hennig, Jan Janßen, Tony Lelièvre, Hao Li, Mitchell Luskin, Noa Marom, Jörg Neugebauer, Feliks Nüske, Joshua Paul, Danny Perez, Giovanni Pinamonti, Petr Plechac, Biswas Rijal, Gideon Simpson, Justin C. Smith, Thomas Swinburne, Anne Marie Z. Tan, Mira Todorova, Dallas R. Trinkle, Stephen Xie, Ping Yang Abstract Recent advances in computational resources and the development of high-throughput frameworks enable the efficient sampling of complicated multivariate functions. This includes energy and electronic property landscapes of inorganic, organic, biomolecular, and hybrid ​ ​ materials and functional nanostructures. Combined with new developments in data science, this ​ leads to a rapidly growing need for numerical methods and a fundamental mathematical ​ ​ understanding of efficient sampling approaches, optimization techniques, hierarchical surrogate models, coarse-graining techniques, and methods for uncertainty quantification. The complexity of these energy and property landscapes originates from their simultaneous dependence on discrete degrees of freedom – i.e., number of atoms and species types – and continuous ones – i.e., the position of atoms. Moreover, dynamical behavior governed by complex landscapes involves a rich hierarchy of timescales and is characterized by rare events
    [Show full text]
  • Protein Functional Landscapes, Dynamics, Allostery 297
    Quarterly Reviews of Biophysics 43, 3 (2010), pp. 295–332. f Cambridge University Press 2010 295 doi:10.1017/S0033583510000119 Printed in the United States of America Protein functional landscapes,dynamics, allostery: a tortuouspathtowards a universal theoretical framework Pavel I. Zhuravlev and Garegin A. Papoian* Department of Chemistry and Biochemistry, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, USA Abstract. Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein’s functional dynamics.
    [Show full text]
  • Evolution, Energy Landscapes and the Paradoxes of Protein Folding
    Biochimie 119 (2015) 218e230 Contents lists available at ScienceDirect Biochimie journal homepage: www.elsevier.com/locate/biochi Review Evolution, energy landscapes and the paradoxes of protein folding Peter G. Wolynes Department of Chemistry and Center for Theoretical Biological Physics, Rice University, Houston, TX 77005, USA article info abstract Article history: Protein folding has been viewed as a difficult problem of molecular self-organization. The search problem Received 15 September 2014 involved in folding however has been simplified through the evolution of folding energy landscapes that are Accepted 11 December 2014 funneled. The funnel hypothesis can be quantified using energy landscape theory based on the minimal Available online 18 December 2014 frustration principle. Strong quantitative predictions that follow from energy landscape theory have been widely confirmed both through laboratory folding experiments and from detailed simulations. Energy Keywords: landscape ideas also have allowed successful protein structure prediction algorithms to be developed. Folding landscape The selection constraint of having funneled folding landscapes has left its imprint on the sequences of Natural selection Structure prediction existing protein structural families. Quantitative analysis of co-evolution patterns allows us to infer the statistical characteristics of the folding landscape. These turn out to be consistent with what has been obtained from laboratory physicochemical folding experiments signaling a beautiful confluence of ge- nomics and chemical physics. © 2014 Elsevier B.V. and Societ e française de biochimie et biologie Moleculaire (SFBBM). All rights reserved. Contents 1. Introduction . ....................... 218 2. General evidence that folding landscapes are funneled . ....................... 221 3. How rugged is the folding landscape e physicochemical approach . ....................... 222 4.
    [Show full text]
  • Behind the Folding Funnel Diagram
    commentary Behind the folding funnel diagram Martin Karplus This Commentary clarifies the meaning of the funnel diagram, which has been widely cited in papers on protein folding. To aid in the analysis of the funnel diagram, this Commentary reviews historical approaches to understanding the mechanism of protein folding. The primary role of free energy in protein folding is discussed, and it is pointed out that the decrease in the configurational entropy as the native state is approached hinders folding, rather than guiding it. Diagrams are introduced that provide a less ambiguous representation of the factors governing the protein folding reaction than the funnel diagram. nderstanding the mechanism fewer) generally fold in times on the order to that shown in Figure 1. It is a schematic by which proteins fold to their of milliseconds to seconds (except in cases representation of how the effective Unative state remains a problem of with special factors that slow the folding, potential energy, implicitly averaged fundamental interest in biology, in spite of such as proline isomerization), there was over solvent interactions (in the vertical the fact that it has been studied for many indeed a paradox. The ultimate statement direction), and the configurational entropy years1. Moreover, now that misfolding has of the paradox was given in the language (in the horizontal direction) of a protein been shown to be the source of a range of of computational complexity5. Many decrease as the native state is approached; diseases, a knowledge of the factors that phenomenological models were proposed picturesque three-dimensional funnels are determine whether a polypeptide chain to show how the conformational space that also being used11.
    [Show full text]
  • Chemical Physics of Protein Folding
    Chemical physics of protein folding Peter G. Wolynesa,1, William A. Eatonb, and Alan R. Fershtc aDepartment of Chemistry and Center for Theoretical Biological Physics, Rice University, Houston, TX 77005; bNational Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892; and cMedical Research Council Laboratory of Molecular Biology, Cambridge CB2 0QH, United Kingdom n 2008, to the consternation of some, logically distinct traps. In other words, the Many of the most intriguing questions of one of the editors of this special issue energy landscape of evolved proteins ap- mechanism can only be addressed feebly at on the “Chemical Physics of Protein pears to be funneled (6, 7). the level of ensemble-averaged experi- I ” “ Folding, was quoted as saying, What The overall funneled nature of the folding ments; thus, folding science has called forth was called the protein folding problem 20 landscape provides a first guess of how some extremely challenging experiments in years ago is solved” (1). One purpose of this folding begins and continues: Proteins fold which molecules are studied individually special issue is to drive home this point. The by assembling primarily native substructures, oneatatime(17–19). On the computa- other, more important purpose is to illus- whereas they only transiently sample mis- tional side, although effective theories of trate how workers on the protein folding folds. This insight explains the success of folding can often use simplified models that problem, by moving beyond their early ob- protein engineering in providing detailed make interesting and surprising predictions session with seeming paradoxes (2), are structural information on the transition state (20, 21), many questions of detail can developing a quantitative understanding ensemble, the so-called “φ-value analysis” now be answered satisfactorily with fully of how the simpler biological structures (8, 9).
    [Show full text]
  • Modeling the Free Energy Landscape of Biomolecules Via Dihedral Angle Principal Component Analysis of Molecular Dynamics Simulations
    Modeling the Free Energy Landscape of Biomolecules via Dihedral Angle Principal Component Analysis of Molecular Dynamics Simulations Dissertation zur Erlangung des Doktorgrades der Naturwissenschaften vorgelegt beim Fachbereich Biochemie, Chemie und Pharmazie der Goethe-Universit¨at in Frankfurt am Main von Alexandros Altis aus Frankfurt am Main Frankfurt am Main 2008 (D 30) vom Fachbereich Biochemie, Chemie und Pharmazie der Goethe-Universit¨at Frankfurt am Main als Dissertation angenommen. Dekan: Prof. Dr. Dieter Steinhilber 1. Gutachter: Prof. Dr. Gerhard Stock 2. Gutachter: JProf. Dr. Karin Hauser Datum der Disputation: ............................................. Contents 1 Introduction 1 2 Dihedral Angle Principal Component Analysis 7 2.1 Introduction to molecular dynamics simulation . ........ 9 2.2 Definition and derivation of principal components . ......... 11 2.3 Circularstatistics ................................ 12 2.4 Dihedral angle principal component analysis (dPCA) . ....... 18 2.5 Asimpleexample-trialanine . 19 2.6 Interpretation of eigenvectors . ..... 22 2.7 ComplexdPCA ................................. 25 2.8 Energylandscapeofdecaalanine. 27 2.9 CartesianPCA ................................. 31 2.10 DirectangularPCA............................... 35 2.11 Correlationanalysis. 39 2.12 Nonlinear principal component analysis . ...... 42 2.13Conclusions ................................... 43 3 Free Energy Landscape 47 3.1 Introduction................................... 47 3.2 Clustering ...................................
    [Show full text]
  • Energy Landscapes for Proteins: from Single Funnels to Multifunctional Systems
    Energy Landscapes for Proteins: From Single Funnels to Multifunctional Systems K. R¨oder,1,† J. A. Joseph,1,† Brooke E. Husic1,2 and D. J. Wales1,∗ 1: Department of Chemistry, University of Cambridge, Lensfield Road, CB2 1EW, Cambridge, UK 2: Current affiliation: Chemistry Department, Stanford University, 333 Campus Drive, Stanford, CA 94305, USA ∗ Corresponding author: [email protected] † These authors contributed equally. Article type: Progress Report Keywords: Energy landscapes, protein folding, multifunnel paradigm In this report we advance the hypothesis that multifunctional systems may be associated with multifunnel potential and free energy landscapes, with particular focus on biomolecules. We compare systems that exhibit single, double, and multiple competing structures, and contrast multifunnel landscapes associated with misfolded amyloidogenic oligomers, which presumably do not arise as an evolutionary target. In this context, intrinsically disordered proteins could be considered intrinsically multifunctional molecules, associated with multifunnel landscapes. Potential energy landscape theory enables biomolecules to be treated in a common framework together with self-organising and multifunctional systems based on inorganic materials, atomic and molecular clusters, crystal polymorphs, and soft matter. 1 Introduction The field of protein science underwent dramatic transformations with the groundbreaking ex- periments of Sanger[1] and Kendrew,[2] characterising the sequence of insulin and the three- dimensional structure of myoglobin,
    [Show full text]