A Consensus View of Protein Dynamics

Total Page:16

File Type:pdf, Size:1020Kb

A Consensus View of Protein Dynamics A consensus view of protein dynamics Manuel Rueda*†, Carles Ferrer-Costa*†, Tim Meyer*†‡, Alberto Pe´ rez*†, Jordi Camps†§, Adam Hospital*†§, Josep Lluis Gelpi´*†‡, and Modesto Orozco*†‡§¶ *Molecular Modelling and Bioinformatics Unit and §Structural Biology Node, Institut de Recerca Biome`dica, Parc Cientı´ficde Barcelona, Josep Samitier 1-5, 08028 Barcelona, Spain; †Computational Biology Program, Barcelona Supercomputing Center, Jordi Girona 31, Edifici Nexus II, 08028 Barcelona, Spain; and ‡Departament de Bioquı´micai Biologia Molecular, Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 645, 08028 Barcelona, Spain Edited by Harold A. Scheraga, Cornell University, Ithaca, NY, and approved November 10, 2006 (received for review July 6, 2006) The dynamics of proteins in aqueous solution has been investi- Table 1. Structures representative of protein metafolds gated through a massive approach based on ‘‘state of the art’’ Exp No. of Atoms in molecular dynamics simulations performed for all protein PDB ID code structure disulfide bridges simulation metafolds using the four most popular force fields (OPLS, CHARMM, AMBER, and GROMOS). A detailed analysis of the 1AGI X-ray 3 23,403͞22,734 massive database of trajectories (>1.5 terabytes of data obtained 1BFG X-ray 0 25,212͞24,509 ͞ using Ϸ50 years of CPU) allowed us to obtain a robust-consensus 1BJ7 X-ray 2 19,371 18,518 1BSN NMR 0 33,753͞32,948 picture of protein dynamics in aqueous solution. 1CHN X-ray 0 18,022͞17,282 1CQY X-ray 0 28,856͞28,310 ͉ ͉ ͉ force field molecular dynamics molecular modeling protein structure 1CSP X-ray 0 13,293͞12,949 1CZT X-ray 1 31,072͞30,167 lexibility is a key determinant of the biological functionality of 1EMR X-ray 0 32,808͞31,885 1FAS X-ray 4 15,709͞15,398 Fproteins. A significant percentage of proteins are unfolded in ͞ the absence of ligand, and many others change their conformation 1FVQ NMR 0 16,322 15,921 1GND X-ray 0 68,647͞66,213 as a result of the presence of other molecules or changes in the 1I6F NMR 4 15,837͞15,552 environment (1–4). Experimental representation of flexibility, 1IL6 NMR 4 27,464͞26,486 even aprioripossible, is in general difficult, which makes atomistic 1J5D NMR 0 23,184͞22,682 simulation the only viable alternative to study this important 1JLI NMR 1 36,307͞35,615 phenomenon for many proteins. Among the different theoretical 1K40 X-ray 0 37,686͞36,883 methods available for description of protein flexibility, molecular 1KTE X-ray 1 19,490͞18,869 dynamics is probably the most powerful (5–13). Since it was first 1KXA X-ray 0 24,969͞24,113 ͞ applied to proteins in the late 1970s, molecular dynamics (MD) has 1LIT X-ray 3 26,289 25,651 1LKI X-ray 3 33,388͞32,403 been largely used to study the dynamics of proteins (14, 15). 1NSO NMR 0 31,774͞31,134 Unfortunately, due to the cost of simulations and the diversity of 1OOI X-ray 3 21,754͞21,057 force fields, a consensus view of protein dynamics has not been yet 1OPC X-ray 0 22,833͞22,245 obtained using this technique. In this paper, we present a systematic 1PDO X-ray 0 30,947͞30,199 study of the most populated protein metafolds [see supporting 1PHT X-ray 0 19,995͞19,554 information (SI) Data Set] using state of the art atomistic MD 1SDF NMR 2 32,579͞32,168 simulation conditions and the four most widely used force fields 1SP2 NMR 0 17,093͞16,930 [OPLS (O), CHARMM (C), AMBER (A), and GROMOS (G)] 1SUR X-ray 0 35,945͞34,691 ͞ (16–25). The result of this massive supercomputer effort (Ͼ1.5 2HVM X-ray 3 32,956 31,550 terabytes of data and a computational equivalent to 50 CPU years) For each protein, we quote the number of disulfide bridges, the origin of is a consensus picture of protein dynamics under conditions close experimental structure, and the number of atoms in the simulation box (A, C, to the physiological ones. O͞G). PDB, Protein Data Bank. Results and Discussion Supporting Information. For further details, see SI Data Set and (see Fig. 1). Thus, differences in radii of gyration between MD SI Figs. 5–17 samplings and experimental structures are Ͻ1%, and average differences in solvent accessible surface are Ϸ5%. The average Force Field Convergence. Previous to any dynamic study, we need to backbone rmsds between simulated and experimental structures determine whether force fields are providing a similar picture of are Ϸ2.0 Å (A, 1.9; C, 2.0; O, 1.9; G, 2.5 Å), close to the thermal protein structure and whether such a picture is similar to that noise of MD simulations. Quite surprisingly, the presence or derived experimentally. Analysis of collected samplings indicates an absence of disulfide bridges does not modify the deviation of our average divergence (␣ in Eq. 1)ofϷ2 Å between force fields samplings from experimental structures, which is, however, depen- (slightly larger deviations are found in G simulations), which, dent on the origin of the experimental structure. Thus, proteins considering the thermal noise of the simulations (⍀ in Eq. 2), suggests a similarity of Ϸ70% between the four samplings and an Ϫ1 average ‘‘effective distance (⍀ ) between them of only 1.4 Å Author contributions: M.R., C.F.-C., and T.M. contributed equally to this work; M.R., C.F.-C., (slightly larger values are always obtained for G-simulations; see A.H., J.L.G., and M.O. designed research; M.R., C.F.-C., T.M., A.P., and J.C. performed Table 2). This finding indicates that all simulations are, in fact, research; M.R., C.F.-C., T.M., and A.P. analyzed data; and M.O. wrote the paper. sampling a similar region of the conformational space. This sug- The authors declare no conflict of interest. gestion is supported by the analysis of the radii of gyration and This article is a PNAS direct submission. solvent accessible surface (differences of Ϸ1% in radii of gyration Abbreviations: MD, molecular dynamics; O, OPLS; C, CHARMM; A, AMBER; G, GROMOS. and 2% in solvent-accessible surface between the four force fields; ¶To whom correspondence should be addressed. E-mail: [email protected]. see Fig. 1). This article contains supporting information online at www.pnas.org/cgi/content/full/ Not only are the structures sampled in the four simulations 0605534104/DC1. similar, but they are also close to the experimental conformations © 2007 by The National Academy of Sciences of the USA 796–801 ͉ PNAS ͉ January 16, 2007 ͉ vol. 104 ͉ no. 3 www.pnas.org͞cgi͞doi͞10.1073͞pnas.0605534104 Downloaded by guest on September 28, 2021 BIOPHYSICS Fig. 1. Average values for rmsd (Upper Left), rRMSD (based on TM-score; Upper Right) with respect to experimental structures, radii of gyration (Lower Left), and solvent accessible surface (Lower Right) for A, C, O, and G force fields. COMPUTER SCIENCES showing the largest displacements with experiment correspond intramolecular hydrogen bonds found in experimental structures is always to structures solved by NMR spectroscopy (see Fig. 1 and very well preserved (see Fig. 2), but obviously, hydrogen bond Table 1). The deviations in these cases are due to displacements in donor and acceptors are more promiscuous in MD than inferred regions of large flexibility, where NMR structures are in general from single experimental structures (see below). All MD simula- poorly defined. This is clear when calculations are repeated using tions preserve very well (90%; see Fig. 2) secondary structure the TM-score, a more robust measure that reduces the weight of (␣-helices and ␤-sheets). In summary: (i) different force fields flexible entities in the rms-fit. In this case, no system is found for provide similar picture of protein structure, and (ii) MD simulations which deviation from experimental data are Ͼ2.5 Å (see Fig. 1). sample regions close to experimental structures. Divergences be- Moderate deviations between MD trajectories and experimental tween simulations and experimental structure signal proteins with model might then reflect either lack of quality in some parts of the very flexible moieties whose exact conformation is often not well experimental model, or just the intrinsic differences between defined. The two prerequisites noted above are then fulfilled, and simulation and experimental conditions. Clearly, the commonly we can then safely analyze the dynamics properties of proteins. made assumption that any deviation from experimental structure of B-factors are the simplest method to analyze local deformability MD simulations is always due to simulation errors is not correct and and have the advantage to be accessible from x-ray data, which can ignore important physical characteristics of proteins. provides us with an additional test on the quality of MD- Native protein–protein contacts are well preserved during sim- trajectories. As noted in Fig. 3, there is a very good agreement ulations (in average, 90% correct predictions). The total number of between experimental and MD-derived B-factors (see Fig. 3). Detailed analysis (see SI Fig. 5) show that such an agreement is especially good for A, C, and O simulations, whereas B-factor Table 2. Similarity indexes distribution is slightly displaced to higher values in G trajectories. In AMBER CHARMM OPLS GROMOS all of the cases, ␤-sheet segments are the most rigid, whereas turns AMBER 1.0 0.7 0.7 0.6 are in general the most flexible ones, in good agreement to what is 1.0 1.0 0.8 0.8 deduced from experimental data. In fact, the only significant CHARMM 1.0 0.7 0.6 difference between MD and x-ray B-factors is that, although 1.0 0.8 0.9 experimental B-factors above 60 Å2 are rare, they are not so OPLS 1.0 0.6 uncommon in MD simulations.
Recommended publications
  • Fast Pressure-Jump All-Atom Simulations and Experiments Reveal Site-Specific Protein Dehydration-Folding Dynamics
    Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics Maxim B. Prigozhina,1,2, Yi Zhangb,1, Klaus Schultenb,c, Martin Gruebelea,b,c,3, and Taras V. Pogorelova,b,d,e,3 aDepartment of Chemistry, University of Illinois at Urbana–Champaign, Urbana, IL 61801; bBeckman Institute for Advanced Science and Technology, Center for Biophysics and Quantitative Biology, University of Illinois at Urbana–Champaign, Urbana, IL 61801; cDepartment of Physics, University of Illinois at Urbana–Champaign, Urbana, IL 61801; dNational Center for Supercomputing Applications, University of Illinois at Urbana–Champaign, Urbana, IL 61801; and eSchool of Chemical Sciences, University of Illinois at Urbana–Champaign, Urbana, IL 61801 Edited by José N. Onuchic, Rice University, Houston, TX, and approved February 8, 2019 (received for review August 30, 2018) As theory and experiment have shown, protein dehydration is a and simulated using MD (18). Desolvated cavities in the protein major contributor to protein folding. Dehydration upon folding core (15) are restored by downward pressure jumps, which fa- can be characterized directly by all-atom simulations of fast cilitate refolding without a change in solvent or thermal motion. pressure drops, which create desolvated pockets inside the In addition, pressure perturbation reduces the effects of protein nascent hydrophobic core. Here, we study pressure-drop refolding aggregation, common in temperature-jump (T-jump) experi- of three λ-repressor fragment (λ6–85) mutants computationally and ments performed near the denaturation temperature (16, 17). experimentally. The three mutants report on tertiary structure for- Finally, pressure-jump (P-jump) experiments can utilize intrinsic mation via different fluorescent helix–helix contact pairs.
    [Show full text]
  • Studying Backbone Torsional Dynamics of Intrinsically Disordered Proteins Using fluorescence Depolarization Kinetics
    J Biosci Vol. 43, No. 3, July 2018, pp. 455–462 Ó Indian Academy of Sciences DOI: 10.1007/s12038-018-9766-1 Studying backbone torsional dynamics of intrinsically disordered proteins using fluorescence depolarization kinetics 1,3 1,2,3 DEBAPRIYA DAS and SAMRAT MUKHOPADHYAY * 1Centre for Protein Science, Design and Engineering, Indian Institute of Science Education and Research (IISER) Mohali, Punjab, India 2Department of Biological Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Punjab, India 3Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Punjab, India *Corresponding author (Email, [email protected]) Published online: 7 June 2018 Intrinsically disordered proteins (IDPs) do not autonomously adopt a stable unique 3D structure and exist as an ensemble of rapidly interconverting structures. They are characterized by significant conformational plasticity and are associated with several biological functions and dysfunctions. The rapid conformational fluctuation is governed by the backbone segmental dynamics arising due to the dihedral angle fluctuation on the Ramachandran /–w conformational space. We discovered that the intrinsic backbone torsional mobility can be monitored by a sensitive fluorescence readout, namely fluorescence depolarization kinetics, of tryptophan in an archetypal IDP such as a-synuclein. This methodology allows us to map the site-specific torsional mobility in the dihedral space within picosecond-nanosecond time range at a low protein concen- tration under the native condition. The characteristic timescale of * 1.4 ns, independent of residue position, represents collective torsional dynamics of dihedral angles (/ and w) of several residues from tryptophan and is independent of overall global tumbling of the protein.
    [Show full text]
  • Linking Predictions of Protein Structure and Disorder Through Molecular Simulation
    The Order-Disorder Continuum: Linking Predictions of Protein Structure and Disorder through Molecular Simulation The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation Hsu, Claire C., Markus J. Buehler, and Anna Tarakanova. "The Order-Disorder Continuum: Linking Predictions of Protein Structure and Disorder through Molecular Simulation." Scientific Reports, 10 (February 2020): 2068. © 2020, The Author(s). As Published http://dx.doi.org/10.1038/s41598-020-58868-w Publisher Springer Science and Business Media LLC Citable link https://hdl.handle.net/1721.1/125166 Terms of Use Creative Commons Attribution 4.0 International license Detailed Terms https://creativecommons.org/licenses/by/4.0/ www.nature.com/scientificreports OPEN The Order-Disorder Continuum: Linking Predictions of Protein Structure and Disorder through Molecular Simulation Claire C. Hsu1, Markus J. Buehler2 & Anna Tarakanova3,4* Intrinsically disordered proteins (IDPs) and intrinsically disordered regions within proteins (IDRs) serve an increasingly expansive list of biological functions, including regulation of transcription and translation, protein phosphorylation, cellular signal transduction, as well as mechanical roles. The strong link between protein function and disorder motivates a deeper fundamental characterization of IDPs and IDRs for discovering new functions and relevant mechanisms. We review recent advances in experimental techniques that have improved identifcation of disordered regions
    [Show full text]
  • The Conformation of Serum Albumin in Solution: a Combined Phosphorescence Depolarization-Hydrodynamic Modeling Study
    2422 Biophysical Journal Volume 80 May 2001 2422–2430 The Conformation of Serum Albumin in Solution: A Combined Phosphorescence Depolarization-Hydrodynamic Modeling Study M. Luisa Ferrer,* Ricardo Duchowicz,* Beatriz Carrasco,† Jose´ Garcı´a de la Torre,† and A. Ulises Acun˜a* *Instituto Quı´mica-Fı´sica “Rocasolano”, Consejo Superior de Investigaciones Cientificas, 28006 Madrid, and †Departamento de Quı´mica-Fı´sica, Universidad de Murcia, 30071 Murcia, Spain ABSTRACT There is a striking disparity between the heart-shaped structure of human serum albumin (HSA) observed in single crystals and the elongated ellipsoid model used for decades to interpret the protein solution hydrodynamics at neutral pH. These two contrasting views could be reconciled if the protein were flexible enough to change its conformation in solution from that found in the crystal. To investigate this possibility we recorded the rotational motions in real time of an erythrosin- bovine serum albumin complex (Er-BSA) over an extended time range, using phosphorescence depolarization techniques. These measurements are consistent with the absence of independent motions of large protein segments in solution, in the time range from nanoseconds to fractions of milliseconds, and give a single rotational correlation time ␾(BSA, 1 cP, 20°C) ϭ 40 Ϯ 2 ns. In addition, we report a detailed analysis of the protein hydrodynamics based on two bead-modeling methods. In the first, BSA was modeled as a triangular prismatic shell with optimized dimensions of 84 ϫ 84 ϫ 84 ϫ 31.5 Å, whereas in the second, the atomic-level structure of HSA obtained from crystallographic data was used to build a much more refined rough-shell model.
    [Show full text]
  • Numerical Evaluation of Protein Global Vibrations at Terahertz Frequencies by Means of Elastic Lattice Models †
    Proceedings Numerical Evaluation of Protein Global Vibrations at Terahertz Frequencies by Means of Elastic † Lattice Models Domenico Scaramozzino *, Giuseppe Lacidogna, Gianfranco Piana and Alberto Carpinteri Department of Structural, Geotechnical and Building Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129 Torino, Italy; [email protected] (G.L.); [email protected] (G.P.); [email protected] (A.C.) * Correspondence: [email protected] † Presented at 1st International Electronic Conference on Applied Sciences, 10–30 November 2020; Available online: https://asec2020.sciforum.net/. Published: 9 November 2020 Abstract: Proteins represent one of the most important building blocks for most biological processes. Their biological mechanisms have been found to correlate significantly with their dynamics, which is commonly investigated through molecular dynamics (MD) simulations. However, important insights on protein dynamics and biological mechanisms have also been obtained via much simpler and computationally efficient calculations based on elastic lattice models (ELMs). The application of structural mechanics approaches, such as modal analysis, to the protein ELMs has allowed to find impressive results in terms of protein dynamics and vibrations. The low-frequency vibrations extracted from the protein ELM are usually found to occur within the terahertz (THz) frequency range and correlate fairly accurately with the observed functional motions. In this contribution, the global vibrations of lysozyme will be investigated by means of a finite element (FE) truss model, and we will show that there exists complete consistency between the proposed FE approach and one of the more well-known ELMs for protein dynamics, the anisotropic network model (ANM). The proposed truss model can consequently be seen as a simple method, easily accessible to the structural mechanics community members, to analyze protein vibrations and their connections with the biological activity.
    [Show full text]
  • Dynamics of Proteins in Solution Biophysics
    Quarterly Reviews of Dynamics of proteins in solution Biophysics Marco Grimaldo1,2, Felix Roosen-Runge1,3, Fajun Zhang2, Frank Schreiber2 cambridge.org/qrb and Tilo Seydel1 1Institut Max von Laue - Paul Langevin, 71 avenue des Martyrs, 38042 Grenoble, France; 2Institut für Angewandte Physik, Universität Tübingen, Auf der Morgenstelle 10, 72076 Tübingen, Germany and 3Division for Physical Review Chemistry, Lund University, Naturvetarvägen 14, 22100 Lund, Sweden Cite this article: Grimaldo M, Roosen-Runge F, Abstract Zhang F, Schreiber F, Seydel T (2019). Dynamics of proteins in solution. Quarterly The dynamics of proteins in solution includes a variety of processes, such as backbone and Reviews of Biophysics 52, e7, 1–63. https:// side-chain fluctuations, interdomain motions, as well as global rotational and translational doi.org/10.1017/S0033583519000027 (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essen- tial transport processes, a detailed mechanistic understanding and monitoring of protein Received: 15 February 2019 Revised: 27 February 2019 dynamics in solution is highly desirable. The hierarchical character of protein dynamics Accepted: 4 March 2019 requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and Key words: focus in particular on results from neutron spectroscopy. We outline the underlying principles Quasielastic neutron spectroscopy; protein diffusion; protein
    [Show full text]
  • Comparisons of Protein Dynamics from Experimental Structure
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Digital Repository @ Iowa State University Biochemistry, Biophysics and Molecular Biology Biochemistry, Biophysics and Molecular Biology Publications 1-27-2018 Comparisons of Protein Dynamics from Experimental Structure Ensembles, Molecular Dynamics Ensembles, and Coarse-Grained Elastic Network Models Kannan Sankar Iowa State University Sambit K. Mishra Iowa State University Robert L. Jernigan Iowa State University, [email protected] Follow this and additional works at: https://lib.dr.iastate.edu/bbmb_ag_pubs Part of the Biochemistry Commons, Bioinformatics Commons, Biophysics Commons, Molecular Biology Commons, and the Structural Biology Commons The ompc lete bibliographic information for this item can be found at https://lib.dr.iastate.edu/ bbmb_ag_pubs/197. For information on how to cite this item, please visit http://lib.dr.iastate.edu/ howtocite.html. This Article is brought to you for free and open access by the Biochemistry, Biophysics and Molecular Biology at Iowa State University Digital Repository. It has been accepted for inclusion in Biochemistry, Biophysics and Molecular Biology Publications by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Comparisons of Protein Dynamics from Experimental Structure Ensembles, Molecular Dynamics Ensembles, and Coarse-Grained Elastic Network Models Abstract Predicting protein motions is important for bridging the gap between protein structure and function. With growing numbers of structures of the same, or closely related proteins becoming available, it is now possible to understand more about the intrinsic dynamics of a protein with principal component analysis (PCA) of the motions apparent within ensembles of experimental structures.
    [Show full text]
  • Nucleolar Proteome Dynamics Analysis
    letters to nature GLT1 activity and immunoblotting against Borrelia spp. in the mouse brain and other sites. Antimicrob. Agents Chemother. 38, 2632–2636 Levels of GLT1 protein were quantified by immunoblots3. Functional glutamate transport (1996). was measured by accumulation of 3H-glutamate in spinal cord slice or crude cortical 15. Chen, W. et al. Expression of a variant form of the glutamate transporter GLT1 in neuronal cultures 25 and in neurons and astrocytes in the rat brain. J. Neurosci. 22, 2142–2152 (2002). synaptosomal membranes . Measurement of total glutamate uptake actually reflects the þ combined physiological activity of all transporter subtypes. GLT1 protein is uniquely 16. Schlag, B. D. et al. Regulation of the glial Na -dependent glutamate transporters by cyclic AMP sensitive to transport inhibition by dihydrokainate (DHK). To estimate the contribution analogs and neurons. Mol. Pharmacol. 53, 355–369 (1998). of GLT1 to transport, aliquots of tissue homogenates were also incubated with 300 mM 17. Guo, H. et al. Increased expression of the glial glutamate transporter EAAT2 modulates excitotoxicity DHK. Non-specific uptake was determined in the presence of 300 mM threo-b- and delays the onset but not the outcome of ALS in mice. Hum. Mol. Genet. 12, 2519–2532 (2003). 18. Romera, C. et al. In vitro ischemic tolerance involves upregulation of glutamate transport partly hydroxyaspartate (THA), at 0 8C and in sodium-free homogenates. mediated by the TACE/ADAM17-tumor necrosis factor-alpha pathway. J. Neurosci. 24, 1350–1357 Generation of GLT1 BAC eGFP transgenic mouse (2004). 19. Spalloni, A. et al. Cu/Zn-superoxide dismutase (GLY93 ! ALA) mutation alters AMPA receptor 26 The BAC transgenic mice were generated as described previously with a shuttle vector subunit expression and function and potentiates kainate-mediated toxicity in motor neurons in provided by N.
    [Show full text]
  • Fluorescence Methods to Study Protein Dynamics in Living Cells
    Traffic 2000 1: 755–762 Munksgaard International Publishers Review Watching Proteins in the Wild: Fluorescence Methods to Study Protein Dynamics in Living Cells Chester Chamberlain and Klaus M. Hahn* advantage of the exploding interest in fluorescent analog cytochemistry. This in turn led to tremendous improvements Department of Cell Biology, BCC 162, The Scripps in the capabilities of the equipment, bringing the technique Research Institute, 10550 N Torrey Pines Rd, La Jolla, CA very much into the mainstream, where it has become an 92037, USA important means to study the mechanisms of many funda- * Corresponding author: K. Hahn, [email protected] mental cellular processes. The original fluorescent analogs were purposefully designed The advent of GFP imaging has led to a revolution in the with fluorophores that would not respond to their environ- study of live cell protein dynamics. Ease of access to ment, enabling precise quantitation of subcellular concentra- fluorescently tagged proteins has led to their wide- tions. The field is now reaching beyond such tagging to spread application and demonstrated the power of analogs whose fluorescence changes to report protein activ- studying protein dynamics in living cells. This has ity. This was initially made possible through covalent labeling spurred development of next generation approaches en- of proteins with dyes that shift their fluorescence spectra in abling not only the visualization of protein movements, response to changes in their protein environment (7,8), but but correlation of a protein’s dynamics with its changing has since been extended through the application of fluores- structural state or ligand binding. Such methods make use of fluorescence resonance energy transfer and dyes cence resonance energy transfer (FRET) using GFP or dyes that report changes in their environment, and take ad- (9–11).
    [Show full text]
  • Protein Dynamics and Entropy: Implications for Protein-Ligand Binding
    University of Pennsylvania ScholarlyCommons Publicly Accessible Penn Dissertations 2015 Protein Dynamics and Entropy: Implications for Protein-Ligand Binding Kyle William Harpole University of Pennsylvania, [email protected] Follow this and additional works at: https://repository.upenn.edu/edissertations Part of the Biochemistry Commons, and the Biophysics Commons Recommended Citation Harpole, Kyle William, "Protein Dynamics and Entropy: Implications for Protein-Ligand Binding" (2015). Publicly Accessible Penn Dissertations. 1756. https://repository.upenn.edu/edissertations/1756 This paper is posted at ScholarlyCommons. https://repository.upenn.edu/edissertations/1756 For more information, please contact [email protected]. Protein Dynamics and Entropy: Implications for Protein-Ligand Binding Abstract The nature of macromolecular interactions has been an area of deep interest for understanding many facets of biology. While a great deal of insight has been gained from structural knowledge, the contribution of protein dynamics to macromolecular interactions is not fully appreciated. This plays out from a thermodynamic perspective as the conformational entropy. The role of conformational entropy in macromolecular interactions has been difficulto t address experimentally. Recently, an empirical calibration has been developed to quantify the conformational entropy of proteins using solution NMR relaxation methods. This method has been demonstrated in two distinct protein systems. The goal of this work is to expand this calibration to assess whether conformational entropy can be effectively quantified from NMR-derived protein dynamics. First, we demonstrate that NMR dynamics do not correlate well between the solid and solution states, suggesting that the relationship between the conformational entropy of proteins is limited to solution state-derived NMR dynamics.
    [Show full text]
  • Protein Dynamics 11/16/16 and 11/18/16
    Protein Dynamics 11/16/16 and 11/18/16 Valerie Daggett [email protected] © Daggett Proteins Many shapes, many sizes, >100,000 in PDB è 807 unique autonomous protein folds Protein Structure • First protein crystal structure solved in the 1950s • Myoglobin • Nobel prize to Perutz and Kendrew Myoglobin Myoglobin Tightly packed! No pathway for O2 → heme! Protein Dynamics • Proteins are not static • Motion is an incontrovertible consequence of existing @ room temperature (or any T > 0 K) • Kinetic energy per atom is ~ 1 kcal/mole @ 298K (25°C) ⇒ several Å/ps • Motion recognized to be important since first crystal structure solved ...everything that living things do can be understood in terms of the jigglings and wigglings of atoms. Richard Feynman Function from Static Structure Dynamics necessary for function Dynamics Necessary for Function Outline of Lecture • Prevalent Protein Motions – Experiment Problems w/extracting detailed structural – Theory information regarding dynamics from expt – Biological Relevance » Is motion just a result of thermal energy and weak interactions, or is such motion functionally important and designed into a protein? Prevalent Protein Motions Type of Motion Extent of Motion (Å) Time Scale Bond Vibration 0.01-0.1 0.01-0.1 ps Side Chain 5-10 E 10-100 ps E Rotation 5 B 108-1012 ps B 1 (0.1-103 ms) Breathing 0.5-2 10 ps-1 ns 2 Relative Motion 1-15 10 ps-1 s 3 - hinge bending - allosteric tran. Local 5-10 < ms, µs, ns Denaturation deadtime 1012 ps/s 109 ns/s 1 Side Chain Rotations • Ring flips • 1H NMR – Rotational freedom
    [Show full text]
  • Proteomics of Phosphorylation and Protein Dynamics During Fertilization
    Proteomics of phosphorylation and protein dynamics PNAS PLUS during fertilization and meiotic exit in the Xenopus egg Marc Preslera, Elizabeth Van Italliea, Allon M. Kleina, Ryan Kunzb, Margaret L. Coughlina, Leonid Peshkina, Steven P. Gygib, Martin Wühra,b,c,d,1, and Marc W. Kirschnera,1 aDepartment of Systems Biology, Harvard Medical School, Boston, MA 02115; bDepartment of Cell Biology, Harvard Medical School, Boston, MA 02115; cDepartment of Molecular Biology, Princeton University, Princeton, NJ 08544; and dThe Lewis–Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08544 Contributed by Marc W. Kirschner, October 27, 2017 (sent for review June 6, 2017; reviewed by William M. Bement, James E. Ferrell, and Matthias Mann) Fertilization releases the meiotic arrest and initiates the events that studies were either qualitative or did not measure phosphorylation. prepare the egg for the ensuing developmental program. Protein With an en-richment step (21–23), it is feasible to measure rel- degradation and phosphorylation are known to regulate protein ative changes in a large number of phosphoproteins. However, activity during this process. However, the full extent of protein loss such studies are of limited utility without measuring site stoi- and phosphoregulation is still unknown. We examined absolute chiometry. Relying on fold change alone for phosphorylation protein and phosphosite dynamics of the fertilization response by studies will not distinguish between 10-fold relative changes from mass spectrometry-based proteomics in electroactivated eggs. To do 1% occupancy to 10% versus 10 to 100%. Several approaches are this, we developed an approach for calculating the stoichiometry of available (24–26), but none of these is able to determine occupan- phosphosites from multiplexed proteomics that is compatible cies of peptides measured with multiplexed MS or that have mul- with dynamic, stable, and multisite phosphorylation.
    [Show full text]