1
|
Maddock RMA, Marsh CO, Johns ST, Rooms LD, Duke PW, van der Kamp MW, Stach JEM, Race PR. Molecular basis of hyper-thermostability in the thermophilic archaeal aldolase MfnB. Extremophiles 2024; 28:42. [PMID: 39215799 PMCID: PMC11365854 DOI: 10.1007/s00792-024-01359-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2024] [Accepted: 08/20/2024] [Indexed: 09/04/2024]
Abstract
Methanogenic archaea are chemolithotrophic prokaryotes that can reduce carbon dioxide with hydrogen gas to form methane. These microorganisms make a significant contribution to the global carbon cycle, with methanogenic archaea from anoxic environments estimated to contribute > 500 million tons of global methane annually. Archaeal methanogenesis is dependent on the methanofurans; aminomethylfuran containing coenzymes that act as the primary C1 acceptor molecule during carbon dioxide fixation. Although the biosynthetic pathway to the methanofurans has been elucidated, structural adaptations which confer thermotolerance to Mfn enzymes from extremophilic archaea are yet to be investigated. Here we focus on the methanofuran biosynthetic enzyme MfnB, which catalyses the condensation of two molecules of glyceralde-3-phosphate to form 4‑(hydroxymethyl)-2-furancarboxaldehyde-phosphate. In this study, MfnB enzymes from the hyperthermophile Methanocaldococcus jannaschii and the mesophile Methanococcus maripaludis have been recombinantly overexpressed and purified to homogeneity. Thermal unfolding studies, together with steady-state kinetic assays, demonstrate thermoadaptation in the M. jannaschii enzyme. Molecular dynamics simulations have been used to provide a structural explanation for the observed properties. These reveal a greater number of side chain interactions in the M. jannaschii enzyme, which may confer protection from heating effects by enforcing spatial residue constraints.
Collapse
Affiliation(s)
- Rosie M A Maddock
- School of Biochemistry, Biomedical Sciences Building, University of Bristol, University Walk, Bristol, BS8 1TD, UK
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, Tyndall Avenue, Bristol, BS8 1TQ, UK
| | - Carl O Marsh
- School of Biochemistry, Biomedical Sciences Building, University of Bristol, University Walk, Bristol, BS8 1TD, UK
| | - Samuel T Johns
- School of Biochemistry, Biomedical Sciences Building, University of Bristol, University Walk, Bristol, BS8 1TD, UK
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, Tyndall Avenue, Bristol, BS8 1TQ, UK
| | - Lynden D Rooms
- School of Biochemistry, Biomedical Sciences Building, University of Bristol, University Walk, Bristol, BS8 1TD, UK
| | - Phillip W Duke
- Defence Science and Technology Laboratory, Porton Down, Salisbury, SP4 0JQ, UK
| | - Marc W van der Kamp
- School of Biochemistry, Biomedical Sciences Building, University of Bristol, University Walk, Bristol, BS8 1TD, UK
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, University of Bristol, Tyndall Avenue, Bristol, BS8 1TQ, UK
| | - James E M Stach
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK
| | - Paul R Race
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK.
| |
Collapse
|
2
|
Fujii Y, Ioka H, Minamoto C, Kurisaki I, Tanaka S, Ohta K, Tominaga K. Vibrational frequency fluctuations of poly(N,N-diethylacrylamide) in the vicinity of coil-to-globule transition studied by two-dimensional infrared spectroscopy and molecular dynamics simulations. J Chem Phys 2024; 161:064903. [PMID: 39120037 DOI: 10.1063/5.0218180] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2024] [Accepted: 07/24/2024] [Indexed: 08/10/2024] Open
Abstract
Poly(N,N-diethylacrylamide) (PdEA), one of the thermoresponsive polymers, in aqueous solutions has attracted much attention because of its characteristic properties, such as coil-to-globule (CG) transition. We performed two-dimensional infrared spectroscopy and molecular dynamics (MD) simulations to understand the hydration dynamics in the vicinity of the CG transition at the molecular level via vibrational frequency fluctuations of the carbonyl stretching modes in the side chains of PdEA. Furthermore, N,N-diethylpropionamide, a repeating monomer unit of PdEA, is also investigated for comparison. From decays of the frequency-frequency time correlation functions (FFTCFs) of the carbonyl stretching modes, we consider that inhomogeneity of the hydration environments originates from various backbone configurations of PdEA. The degree of the inhomogeneity depends on temperature. Hydration water molecules near the carbonyl groups are influenced by the confinements of the polymers. The restricted reorientation of the embedded water, the local torsions of the backbone, and the rearrangement of the whole structure contribute to the slow spectral diffusion. By performing MD simulations, we calculated the FFTCFs and dynamical quantities, such as fluctuations of the dihedral angles of the backbone and the orientation of the hydration water molecules. The simulated FFTCFs match well with the experimental results, indicating that the retarded water reorientations via the excluded volume effect play an important role in the vibrational frequency fluctuations of the carbonyl stretching mode. It is also found the embedded water molecules are influenced by the local torsions of the backbone structure within the time scales of the spectral diffusion.
Collapse
Affiliation(s)
- Yuki Fujii
- Department of Chemistry, Graduate School of Science, Kobe University, Rokkodai-cho 1-1, Nada, Kobe 657-8501, Japan
| | - Hikaru Ioka
- Department of Chemistry, Graduate School of Science, Kobe University, Rokkodai-cho 1-1, Nada, Kobe 657-8501, Japan
| | - Chihiro Minamoto
- Department of Applied Chemistry and Biotechnology, Niihama National College of Technology, Yakumo-cho 7-1, Niihama, Ehime 792-8580, Japan
| | - Ikuo Kurisaki
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
| | - Shigenori Tanaka
- Department of Computational Science, Graduate School of System Informatics, Kobe University, Rokkodai-cho 1-1, Nada, Kobe 657-8501, Japan
| | - Kaoru Ohta
- Molecular Photoscience Research Center, Kobe University, Rokkodai-cho 1-1, Nada, Kobe 657-8501, Japan
| | - Keisuke Tominaga
- Department of Chemistry, Graduate School of Science, Kobe University, Rokkodai-cho 1-1, Nada, Kobe 657-8501, Japan
- Molecular Photoscience Research Center, Kobe University, Rokkodai-cho 1-1, Nada, Kobe 657-8501, Japan
| |
Collapse
|
3
|
Kruchinin SE, Kislinskaya EE, Chuev GN, Fedotova MV. Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor. Int J Mol Sci 2022; 23:ijms232314785. [PMID: 36499117 PMCID: PMC9737982 DOI: 10.3390/ijms232314785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 11/20/2022] [Accepted: 11/23/2022] [Indexed: 11/29/2022] Open
Abstract
Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.
Collapse
Affiliation(s)
- Sergey E. Kruchinin
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Ekaterina E. Kislinskaya
- Department of Fundamental and Applied Chemistry, Institute of Mathematics, Information Technology and Science, Ivanovo State University, Ermak St., 39, 153025 Ivanovo, Russia
| | - Gennady N. Chuev
- Institute of Theoretical and Experimental Biophysics, The Russian Academy of Sciences, Institutskaya St., Pushchino, 142290 Moscow, Russia
- Correspondence: (G.N.C.); (M.V.F.)
| | - Marina V. Fedotova
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
- Correspondence: (G.N.C.); (M.V.F.)
| |
Collapse
|
4
|
Sun Q, Fu Y, Wang W. Temperature effects on hydrophobic interactions: Implications for protein unfolding. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111550] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
5
|
Jin T, Long F, Zhang Q, Zhuang W. Site-Specific Water Dynamics in the First Hydration Layer of an Anti-Freeze Glyco-Protein: A Simulation Study. Phys Chem Chem Phys 2022; 24:21165-21177. [DOI: 10.1039/d2cp00883a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Antifreeze glycoproteins (AFGPs) inhibit ice recrystallization by a mechanism remaining largely elusive. Dynamics of AFGPs’ hydration water and its involvement in the antifreeze activity, for instance, have not been identified...
Collapse
|
6
|
Mesophilic Pyrophosphatase Function at High Temperature: A Molecular Dynamics Simulation Study. Biophys J 2020; 119:142-150. [PMID: 32533942 DOI: 10.1016/j.bpj.2020.05.021] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2020] [Revised: 05/06/2020] [Accepted: 05/15/2020] [Indexed: 02/06/2023] Open
Abstract
The mesophilic inorganic pyrophosphatase from Escherichia coli (EcPPase) retains function at 353 K, the physiological temperature of hyperthermophilic Thermococcus thioreducens, whereas the homolog protein (TtPPase) from this hyperthermophilic organism cannot function at room temperature. To explain this asymmetric behavior, we examined structural and dynamical properties of the two proteins using molecular dynamics simulations. The global flexibility of TtPPase is significantly higher than its mesophilic homolog at all tested temperature/pressure conditions. However, at 353 K, EcPPase reduces its solvent-exposed surface area and increases subunit compaction while maintaining flexibility in its catalytic pocket. In contrast, TtPPase lacks this adaptability and has increased rigidity and reduced protein/water interactions in its catalytic pocket at room temperature, providing a plausible explanation for its inactivity near room temperature.
Collapse
|
7
|
Timr S, Madern D, Sterpone F. Protein thermal stability. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:239-272. [PMID: 32145947 DOI: 10.1016/bs.pmbts.2019.12.007] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.
Collapse
Affiliation(s)
- Stepan Timr
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | | | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
8
|
Trofimov YA, Krylov NA, Efremov RG. Confined Dynamics of Water in Transmembrane Pore of TRPV1 Ion Channel. Int J Mol Sci 2019; 20:ijms20174285. [PMID: 31480555 PMCID: PMC6747475 DOI: 10.3390/ijms20174285] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2019] [Revised: 08/28/2019] [Accepted: 08/29/2019] [Indexed: 11/23/2022] Open
Abstract
Solvation effects play a key role in chemical and biological processes. The microscopic properties of water near molecular surfaces are radically different from those in the bulk. Furthermore, the behavior of water in confined volumes of a nanometer scale, including transmembrane pores of ion channels, is especially nontrivial. Knowledge at the molecular level of structural and dynamic parameters of water in such systems is necessary to understand the mechanisms of ion channels functioning. In this work, the results of molecular dynamics (MD) simulations of water in the pore and selectivity filter domains of TRPV1 (Transient Receptor Potential Vanilloid type 1) membrane channel are considered. These domains represent nanoscale volumes with strongly amphiphilic walls, where physical behavior of water radically differs from that of free hydration (e.g., at protein interfaces) or in the bulk. Inside the pore and filter domains, water reveals a very heterogeneous spatial distribution and unusual dynamics: It forms compact areas localized near polar groups of particular residues. Residence time of water molecules in such areas is at least 1.5 to 3 times larger than that observed for similar groups at the protein surface. Presumably, these water “blobs” play an important role in the functional activity of TRPV1. In particular, they take part in hydration of the hydrophobic TRPV1 pore by localizing up to six waters near the so-called “lower gate” of the channel and reducing by this way the free energy barrier for ion and water transport. Although the channel is formed by four identical protein subunits, which are symmetrically packed in the initial experimental 3D structure, in the course of MD simulations, hydration of the same amino acid residues of individual subunits may differ significantly. This greatly affects the microscopic picture of the distribution of water in the channel and, potentially, the mechanism of its functioning. Therefore, reconstruction of the full picture of TRPV1 channel solvation requires thorough atomistic simulations and analysis. It is important that the naturally occurring porous volumes, like ion-conducting protein domains, reveal much more sophisticated and fine-tuned regulation of solvation than, e.g., artificially designed carbon nanotubes.
Collapse
Affiliation(s)
- Yury A Trofimov
- M.M. Shemyakin & Yu.A. Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Miklukho-Maklaya Street, 16/10, 117997 Moscow, Russia
- National Research University Higher School of Economics, Myasnitskaya ul. 20, 101000 Moscow, Russia
- National Research Nuclear University Moscow Engineering Physics Institute, Kashirskoe Shosse, 31, 115409 Moscow, Russia
| | - Nikolay A Krylov
- M.M. Shemyakin & Yu.A. Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Miklukho-Maklaya Street, 16/10, 117997 Moscow, Russia
- National Research University Higher School of Economics, Myasnitskaya ul. 20, 101000 Moscow, Russia
| | - Roman G Efremov
- M.M. Shemyakin & Yu.A. Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Miklukho-Maklaya Street, 16/10, 117997 Moscow, Russia.
- National Research University Higher School of Economics, Myasnitskaya ul. 20, 101000 Moscow, Russia.
- Moscow Institute of Physics and Technology (State University), Dolgoprudny, 141701 Moscow, Russia.
| |
Collapse
|
9
|
Arya S, Singh AK, Bhasne K, Dogra P, Datta A, Das P, Mukhopadhyay S. Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein. Biophys J 2019; 114:2540-2551. [PMID: 29874605 DOI: 10.1016/j.bpj.2018.04.028] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2017] [Revised: 04/16/2018] [Accepted: 04/17/2018] [Indexed: 10/14/2022] Open
Abstract
Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson's disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.
Collapse
Affiliation(s)
- Shruti Arya
- Centre for Protein Science, Design and Engineering, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India; Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India
| | - Avinash K Singh
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, India
| | - Karishma Bhasne
- Centre for Protein Science, Design and Engineering, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India; Department of Biological Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India
| | - Priyanka Dogra
- Centre for Protein Science, Design and Engineering, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India; Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India
| | - Anindya Datta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai, India.
| | - Payel Das
- Data Science Department, IBM Thomas J. Watson Research Center, Yorktown Heights, New York.
| | - Samrat Mukhopadhyay
- Centre for Protein Science, Design and Engineering, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India; Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India; Department of Biological Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, Punjab, India.
| |
Collapse
|
10
|
Abstract
Based on molecular dynamics simulations of four globular proteins in dilute aqueous solution, with three different water models, we examine several, essentially geometrical, aspects of the protein-water interface that remain controversial or incompletely understood. First, we compare different hydration shell definitions, based on spatial or topological proximity criteria. We find that the best method for constructing monolayer shells with nearly complete coverage is to use a 5 Å water-carbon cutoff and a 4 Å water-water cutoff. Using this method, we determine a mean interfacial water area of 11.1 Å2 which appears to be a universal property of the protein-water interface. We then analyze the local coordination and packing density of water molecules in the hydration shells and in subsets of the first shell. The mean polar water coordination number in the first shell remains within 1% of the bulk-water value, and it is 5% lower in the nonpolar part of the first shell. The local packing density is obtained from additively weighted Voronoi tessellation, arguably the most physically realistic method for allocating space between protein and water. We find that water in all parts of the first hydration shell, including the nonpolar part, is more densely packed than in the bulk, with a shell-averaged density excess of 6% for all four proteins. We suggest reasons why this value differs from previous experimental and computational results, emphasizing the importance of a realistic placement of the protein-water dividing surface and the distinction between spatial correlation and packing density. The protein-induced perturbation of water coordination and packing density is found to be short-ranged, with an exponential decay "length" of 0.6 shells. We also compute the protein partial volume, analyze its decomposition, and argue against the relevance of electrostriction.
Collapse
Affiliation(s)
- Filip Persson
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | - Pär Söderhjelm
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | - Bertil Halle
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| |
Collapse
|
11
|
Relative Contributions of Core Protein and Solvation Shell in the Terahertz Dielectric Properties of Protein Solutions. J Phys Chem B 2017; 121:9508-9512. [DOI: 10.1021/acs.jpcb.7b06442] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
12
|
Abstract
The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.
Collapse
Affiliation(s)
- Damien Laage
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
| | - Thomas Elsaesser
- Max-Born-Institut
für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - James T. Hynes
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
- Department
of Chemistry and Biochemistry, University
of Colorado, Boulder, Colorado 80309, United
States
| |
Collapse
|
13
|
Atamas N, Bardik V, Bannikova A, Grishina O, Lugovskoi E, Lavoryk S, Makogonenko Y, Korolovych V, Nerukh D, Paschenko V. The effect of water dynamics on conformation changes of albumin in pre-denaturation state: photon correlation spectroscopy and simulation. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.01.017] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
14
|
Fisette O, Päslack C, Barnes R, Isas JM, Langen R, Heyden M, Han S, Schäfer LV. Hydration Dynamics of a Peripheral Membrane Protein. J Am Chem Soc 2016; 138:11526-35. [PMID: 27548572 DOI: 10.1021/jacs.6b07005] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Water dynamics in the hydration shell of the peripheral membrane protein annexin B12 were studied using MD simulations and Overhauser DNP-enhanced NMR. We show that retardation of water motions near phospholipid bilayers is extended by the presence of a membrane-bound protein, up to around 10 Å above that protein. Near the membrane surface, electrostatic interactions with the lipid head groups strongly slow down water dynamics, whereas protein-induced water retardation is weaker and dominates only at distances beyond 10 Å from the membrane surface. The results can be understood from a simple model based on additive contributions from the membrane and the protein to the activation free energy barriers of water diffusion next to the biomolecular surfaces. Furthermore, analysis of the intermolecular vibrations of the water network reveals that retarded water motions near the membrane shift the vibrational modes to higher frequencies, which we used to identify an entropy gradient from the membrane surface toward the bulk water. Our results have implications for processes that take place at lipid membrane surfaces, including molecular recognition, binding, and protein-protein interactions.
Collapse
Affiliation(s)
- Olivier Fisette
- Center for Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University , 44780 Bochum, Germany
| | - Christopher Päslack
- Center for Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University , 44780 Bochum, Germany.,Max-Planck Institut für Kohlenforschung , 45470 Mülheim an der Ruhr, Germany
| | - Ryan Barnes
- Department of Chemistry and Biochemistry and Department of Chemical Engineering, University of California, Santa Barbara , Santa Barbara, California 93106, United States
| | - J Mario Isas
- Department of Biochemistry and Molecular Biology, Zilkha Neurogenetic Institute, Keck School of Medicine, University of Southern California , Los Angeles, California 90089, United States
| | - Ralf Langen
- Department of Biochemistry and Molecular Biology, Zilkha Neurogenetic Institute, Keck School of Medicine, University of Southern California , Los Angeles, California 90089, United States
| | - Matthias Heyden
- Max-Planck Institut für Kohlenforschung , 45470 Mülheim an der Ruhr, Germany
| | - Songi Han
- Department of Chemistry and Biochemistry and Department of Chemical Engineering, University of California, Santa Barbara , Santa Barbara, California 93106, United States
| | - Lars V Schäfer
- Center for Theoretical Chemistry, Faculty of Chemistry and Biochemistry, Ruhr-University , 44780 Bochum, Germany
| |
Collapse
|
15
|
Brotzakis ZF, Groot CCM, Brandeburgo WH, Bakker HJ, Bolhuis PG. Dynamics of Hydration Water around Native and Misfolded α-Lactalbumin. J Phys Chem B 2016; 120:4756-66. [DOI: 10.1021/acs.jpcb.6b02592] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Z. F. Brotzakis
- Van’t
Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science
Park 904, 1098 XH Amsterdam, The Netherlands
| | - C. C. M. Groot
- FOM Institute AMOLF, Science
Park 104, 1098 XG Amsterdam, The Netherlands
| | - W. H. Brandeburgo
- Van’t
Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science
Park 904, 1098 XH Amsterdam, The Netherlands
| | - H. J. Bakker
- FOM Institute AMOLF, Science
Park 104, 1098 XG Amsterdam, The Netherlands
| | - P. G. Bolhuis
- Van’t
Hoff Institute for Molecular Sciences, Universiteit van Amsterdam, Science
Park 904, 1098 XH Amsterdam, The Netherlands
| |
Collapse
|
16
|
Chakraborty D, Taly A, Sterpone F. Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins. J Phys Chem B 2015; 119:12760-70. [PMID: 26335353 DOI: 10.1021/acs.jpcb.5b05791] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.
Collapse
Affiliation(s)
- Debashree Chakraborty
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité , 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Antoine Taly
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité , 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité , 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
17
|
Khodadadi S, Sokolov AP. Protein dynamics: from rattling in a cage to structural relaxation. SOFT MATTER 2015; 11:4984-4998. [PMID: 26027652 DOI: 10.1039/c5sm00636h] [Citation(s) in RCA: 91] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We present an overview of protein dynamics based mostly on results of neutron scattering, dielectric relaxation spectroscopy and molecular dynamics simulations. We identify several major classes of protein motions on the time scale from faster than picoseconds to several microseconds, and discuss the coupling of these processes to solvent dynamics. Our analysis suggests that the microsecond backbone relaxation process might be the main structural relaxation of the protein that defines its glass transition temperature, while faster processes present some localized secondary relaxations. Based on the overview, we formulate a general picture of protein dynamics and discuss the challenges in this field.
Collapse
Affiliation(s)
- S Khodadadi
- Faculty of Applied Sciences, Delft University of Technology, Delft, The Netherlands
| | | |
Collapse
|
18
|
Sterpone F, Derreumaux P, Melchionna S. Protein Simulations in Fluids: Coupling the OPEP Coarse-Grained Force Field with Hydrodynamics. J Chem Theory Comput 2015; 11:1843-53. [PMID: 26574390 PMCID: PMC5242371 DOI: 10.1021/ct501015h] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
A novel simulation framework that integrates the OPEP coarse-grained (CG) model for proteins with the Lattice Boltzmann (LB) methodology to account for the fluid solvent at mesoscale level is presented. OPEP is a very efficient, water-free and electrostatic-free force field that reproduces at quasi-atomistic detail processes like peptide folding, structural rearrangements, and aggregation dynamics. The LB method is based on the kinetic description of the solvent in order to solve the fluid mechanics under a wide range of conditions, with the further advantage of being highly scalable on parallel architectures. The capabilities of the approach are presented, and it is shown that the strategy is effective in exploring the role of hydrodynamics on protein relaxation and peptide aggregation. The end result is a strategy for modeling systems of thousands of proteins, such as in the case of dense protein suspensions. The future perspectives of the multiscale approach are also discussed.
Collapse
Affiliation(s)
- Fabio Sterpone
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | | |
Collapse
|
19
|
Kalimeri M, Derreumaux P, Sterpone F. Are coarse-grained models apt to detect protein thermal stability? The case of OPEP force field. JOURNAL OF NON-CRYSTALLINE SOLIDS 2015; 407:494-501. [PMID: 28100926 PMCID: PMC5238951 DOI: 10.1016/j.jnoncrysol.2014.07.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
We present the first investigation of the kinetic and thermodynamic stability of two homologous thermophilic and mesophilic proteins based on the coarse-grained model OPEP. The object of our investigation is a pair of G-domains of relatively large size, 200 amino acids each, with an experimental stability gap of about 40 K. The OPEP force field is able to maintain stable the fold of these relatively large proteins within the hundrend-nanosecond time scale without including external constraints. This makes possible to characterize the conformational landscape of the folded protein as well as to explore the unfolding. In agreement with all-atom simulations used as a reference, we show that the conformational landscape of the thermophilic protein is characterized by a larger number of substates with slower dynamics on the network of states and more resilient to temperature increase. Moreover, we verify the stability gap between the two proteins using replica-exchange simulations and estimate a difference between the melting temperatures of about 23 K, in fair agreement with experiment. The detailed investigation of the unfolding thermodynamics, allows to gain insight into the mechanism underlying the enhanced stability of the thermophile relating it to a smaller heat capacity of unfolding.
Collapse
Affiliation(s)
- Maria Kalimeri
- Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, France
- Institut Universitaire de France, 103 Boulevard Saint-Michel, 75005, Paris, France
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, France
- Corresponding author.
| |
Collapse
|
20
|
Rahaman O, Kalimeri M, Melchionna S, Hénin J, Sterpone F. Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains. J Phys Chem B 2014; 119:8939-49. [PMID: 25317828 DOI: 10.1021/jp507571u] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.
Collapse
Affiliation(s)
- Obaidur Rahaman
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Maria Kalimeri
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Simone Melchionna
- ‡CNR-IPCF, Consiglio Nazionale delle Ricerche, Physics Dept., Univ. La Sapienza, P.le A. Moro 2, 00185, Rome, Italy
| | - Jérôme Hénin
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| | - Fabio Sterpone
- †Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Univ. Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005, Paris, France
| |
Collapse
|
21
|
Vajda T, Perczel A. Role of water in protein folding, oligomerization, amyloidosis and miniprotein. J Pept Sci 2014; 20:747-59. [PMID: 25098401 DOI: 10.1002/psc.2671] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2014] [Revised: 06/03/2014] [Accepted: 06/06/2014] [Indexed: 01/02/2023]
Abstract
The essential involvement of water in most fundamental extra-cellular and intracellular processes of proteins is critically reviewed and evaluated in this article. The role of water in protein behavior displays structural ambivalence; it can protect the disordered peptide-chain by hydration or helps the globular chain-folding, but promotes also the protein aggregation, as well (see: diseases). A variety of amyloid diseases begins as benign protein monomers but develops then into toxic amyloid aggregates of fibrils. Our incomplete knowledge of this process emphasizes the essential need to reveal the principles of governing this oligomerization. To understand the biophysical basis of the simpler in vitro amyloid formation may help to decipher also the in vivo way. Nevertheless, to ignore the central role of the water's effect among these events means to receive an uncompleted picture of the true phenomenon. Therefore this review represents a stopgap role, because the most published studies--with a few exceptions--have been neglected the crucial importance of water in the protein research. The following questions are discussed from the water's viewpoint: (i) interactions between water and proteins, (ii) protein hydration/dehydration, (iii) folding of proteins and miniproteins, (iv) peptide/protein oligomerization, and (v) amyloidosis.
Collapse
Affiliation(s)
- Tamás Vajda
- MTA-ELTE Protein Modelling Research Group, Eötvös Loránd University and Laboratory of Structural Chemistry & Biology, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest, 1117, Hungary
| | | |
Collapse
|
22
|
Fogarty A, Laage D. Water dynamics in protein hydration shells: the molecular origins of the dynamical perturbation. J Phys Chem B 2014; 118:7715-29. [PMID: 24479585 PMCID: PMC4103960 DOI: 10.1021/jp409805p] [Citation(s) in RCA: 169] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2013] [Revised: 12/27/2013] [Indexed: 02/08/2023]
Abstract
Protein hydration shell dynamics play an important role in biochemical processes including protein folding, enzyme function, and molecular recognition. We present here a comparison of the reorientation dynamics of individual water molecules within the hydration shell of a series of globular proteins: acetylcholinesterase, subtilisin Carlsberg, lysozyme, and ubiquitin. Molecular dynamics simulations and analytical models are used to access site-resolved information on hydration shell dynamics and to elucidate the molecular origins of the dynamical perturbation of hydration shell water relative to bulk water. We show that all four proteins have very similar hydration shell dynamics, despite their wide range of sizes and functions, and differing secondary structures. We demonstrate that this arises from the similar local surface topology and surface chemical composition of the four proteins, and that such local factors alone are sufficient to rationalize the hydration shell dynamics. We propose that these conclusions can be generalized to a wide range of globular proteins. We also show that protein conformational fluctuations induce a dynamical heterogeneity within the hydration layer. We finally address the effect of confinement on hydration shell dynamics via a site-resolved analysis and connect our results to experiments via the calculation of two-dimensional infrared spectra.
Collapse
Affiliation(s)
- Aoife
C. Fogarty
- Department
of Chemistry, UMR ENS-CNRS-UPMC 8640, École
Normale Supérieure, 24 rue Lhomond, 75005 Paris, France
| | - Damien Laage
- Department
of Chemistry, UMR ENS-CNRS-UPMC 8640, École
Normale Supérieure, 24 rue Lhomond, 75005 Paris, France
| |
Collapse
|
23
|
Arthur EJ, King JT, Kubarych KJ, Brooks CL. Heterogeneous preferential solvation of water and trifluoroethanol in homologous lysozymes. J Phys Chem B 2014; 118:8118-27. [PMID: 24823618 PMCID: PMC4216199 DOI: 10.1021/jp501132z] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
![]()
Cytoplasmic
osmolytes can significantly alter the thermodynamic
and kinetic properties of proteins relative to those under dilute
solution conditions. Spectroscopic experiments of lysozymes in cosolvents
indicate that such changes may arise from the heterogeneous, site-specific
hydrophobic interactions between protein surface residues and individual
solvent molecules. In pursuit of an accurate and predictive model
for explaining biomolecular interactions, we study the averaged structural
characteristics of mixed solvents with homologous lysozyme solutes
using all-atom molecular dynamics. By observing the time-averaged
densities of different aqueous solutions of trifluoroethanol, we deduce
trends in the heterogeneous solvent interactions over each protein’s
surface, and investigate how the homology of protein structure does
not necessarily translate to similarities in solvent structure and
composition—even when observing identical side chains.
Collapse
Affiliation(s)
- Evan J Arthur
- Department of Chemistry and ‡Biophysics Program, University of Michigan , 930 N. University Avenue , Ann Arbor, Michigan 48109-1055, USA
| | | | | | | |
Collapse
|
24
|
Kalimeri M, Rahaman O, Melchionna S, Sterpone F. How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain. J Phys Chem B 2013; 117:13775-85. [PMID: 24087838 DOI: 10.1021/jp407078z] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Proteins from thermophilic organisms are stable and functional well above ambient temperature. Understanding the molecular mechanism underlying such a resistance is of crucial interest for many technological applications. For some time, thermal stability has been assumed to correlate with high mechanical rigidity of the protein matrix. In this work we address this common belief by carefully studying a pair of homologous G-domain proteins, with their melting temperatures differing by 40 K. To probe the thermal-stability content of the two proteins we use extensive simulations covering the microsecond time range and employ several different indicators to assess the salient features of the conformational landscape and the role of internal fluctuations at ambient condition. At the atomistic level, while the magnitude of fluctuations is comparable, the distribution of flexible and rigid stretches of amino-acids is more regular in the thermophilic protein causing a cage-like correlation of amplitudes along the sequence. This caging effect is suggested to favor stability at high T by confining the mechanical excitations. Moreover, it is found that the thermophilic protein, when folded, visits a higher number of conformational substates than the mesophilic homologue. The entropy associated with the occupation of the different substates and the thermal resilience of the protein intrinsic compressibility provide a qualitative insight on the thermal stability of the thermophilic protein as compared to its mesophilic homologue. Our findings potentially open the route to new strategies in the design of thermostable proteins.
Collapse
Affiliation(s)
- Maria Kalimeri
- Laboratoire de Biochimie Théorique, IBPC, CNRS, UPR9080, Université Paris Diderot , Sorbonne Paris Cité, France
| | | | | | | |
Collapse
|
25
|
Fogarty AC, Duboué-Dijon E, Sterpone F, Hynes JT, Laage D. Biomolecular hydration dynamics: a jump model perspective. Chem Soc Rev 2013; 42:5672-83. [DOI: 10.1039/c3cs60091b] [Citation(s) in RCA: 90] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|