1
|
Ye Y, Zheng L, Hong L, García Sakai V, de Souza NR, Teng D, Wu B, Xu Y, Cai J, Liu Z. Direct Observation of the Mutual Coupling Effect in the Protein-Water-Glycerol Mixture by Combining Neutron Scattering and Selective Deuteration. J Phys Chem B 2024; 128:405-414. [PMID: 38183282 DOI: 10.1021/acs.jpcb.3c05135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2024]
Abstract
Numerous studies have discussed the impact of cosolvents on the structure, dynamics, and stability of proteins in aqueous solutions. However, the dynamics of cosolvents in the protein-water-cosolvent ternary system is largely unexplored in experiments due to technical difficulty. Consequently, a comprehensive understanding of the interplay among proteins, water, and cosolvents is still lacking. Here, we employed selective deuteration and neutron scattering techniques to characterize the individual motions of each component in the protein/water/glycerol (GLY) mixture across various temperatures. The consistent dynamic onset temperatures and the correlation between the MSD of the protein and the viscosity of solvents revealed the mutual coupling effects among the three components. Furthermore, our experimental and simulation results showed that the hydrogen bond relaxation energy barrier in the ternary system is ∼43 kJ/mol, whereas in the protein-water binary system it is merely ∼35 kJ/mol. Therefore, we suggest that GLY can enhance hydrogen bond interactions in the ternary system through the mutual coupling effect, thereby serving as one of the protective mechanisms of protein preservation by GLY.
Collapse
Affiliation(s)
- Yongfeng Ye
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Lirong Zheng
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Liang Hong
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
- Shanghai National Centre for Applied Mathematics (SJTU Center), MOE-LSC, Shanghai Jiao Tong University, Shanghai 200240, China
- Zhangjiang Institute for Advanced Study, Shanghai Jiao Tong University, Shanghai 201203, China
- Shanghai Artificial Intelligence Laboratory, Shanghai 200232, China
| | - Victoria García Sakai
- ISIS Pulsed Neutron and Muon Source, Rutherford Appleton Laboratory, Science & Technology Facilities Council, Didcot OX11 0QX, U.K
| | - Nicolas R de Souza
- Australian Nuclear Science and Technology (ANSTO), Locked Bag 2001, Kirrawee DC, New South Wales 2232, Australia
| | - Dahong Teng
- Organ Transplantation Center, Fujian Medical University Union Hospital, Fuzhou 350000, China
| | - Bin Wu
- Organ Transplantation Center, Fujian Medical University Union Hospital, Fuzhou 350000, China
| | - Yichao Xu
- Organ Transplantation Center, Fujian Medical University Union Hospital, Fuzhou 350000, China
| | - Jinzhen Cai
- Organ Transplantation Center, Fujian Medical University Union Hospital, Fuzhou 350000, China
| | - Zhuo Liu
- Institute of Natural Sciences, Shanghai Jiao Tong University, Shanghai 200240, China
- Shanghai National Centre for Applied Mathematics (SJTU Center), MOE-LSC, Shanghai Jiao Tong University, Shanghai 200240, China
- Zhangjiang Institute for Advanced Study, Shanghai Jiao Tong University, Shanghai 201203, China
- Shanghai Artificial Intelligence Laboratory, Shanghai 200232, China
| |
Collapse
|
2
|
Filianina M, Bin M, Berkowicz S, Reiser M, Li H, Timmermann S, Blankenburg M, Amann-Winkel K, Gutt C, Perakis F. Nanocrystallites Modulate Intermolecular Interactions in Cryoprotected Protein Solutions. J Phys Chem B 2023. [PMID: 37399586 DOI: 10.1021/acs.jpcb.3c02413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2023]
Abstract
Studying protein interactions at low temperatures has important implications for optimizing cryostorage processes of biological tissue, food, and protein-based drugs. One of the major issues is related to the formation of ice nanocrystals, which can occur even in the presence of cryoprotectants and can lead to protein denaturation. The presence of ice nanocrystals in protein solutions poses several challenges since, contrary to microscopic ice crystals, they can be difficult to resolve and can complicate the interpretation of experimental data. Here, using a combination of small- and wide-angle X-ray scattering (SAXS and WAXS), we investigate the structural evolution of concentrated lysozyme solutions in a cryoprotected glycerol-water mixture from room temperature (T = 300 K) down to cryogenic temperatures (T = 195 K). Upon cooling, we observe a transition near the melting temperature of the solution (T ≈ 245 K), which manifests both in the temperature dependence of the scattering intensity peak position reflecting protein-protein length scales (SAXS) and the interatomic distances within the solvent (WAXS). Upon thermal cycling, a hysteresis is observed in the scattering intensity, which is attributed to the formation of nanocrystallites in the order of 10 nm. The experimental data are well described by the two-Yukawa model, which indicates temperature-dependent changes in the short-range attraction of the protein-protein interaction potential. Our results demonstrate that the nanocrystal growth yields effectively stronger protein-protein attraction and influences the protein pair distribution function beyond the first coordination shell.
Collapse
Affiliation(s)
- Mariia Filianina
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Maddalena Bin
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Sharon Berkowicz
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Mario Reiser
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| | - Hailong Li
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
- Max Plank Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Sonja Timmermann
- Department of Physics, Universität Siegen, Walter-Flex-Strasse 3, 57072 Siegen, Germany
| | - Malte Blankenburg
- Deutsches Elektronen-Synchrotron (DESY), Notkestrasse 85, 22607 Hamburg, Germany
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
- Max Plank Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
- Institute of Physics, Johannes Gutenberg University, 55128 Mainz, Germany
| | - Christian Gutt
- Department of Physics, Universität Siegen, Walter-Flex-Strasse 3, 57072 Siegen, Germany
| | - Fivos Perakis
- Department of Physics, AlbaNova University Center, Stockholm University, S-106 91 Stockholm, Sweden
| |
Collapse
|
3
|
Ji X, Wang J, Wang T, Huang X, Li X, Wang N, Huang Y, Li R, Zhao B, Zhang T, Hao H. Understanding the Role of Solvent on Regulating Crystal Habit. CrystEngComm 2022. [DOI: 10.1039/d1ce01486b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Crystal habit is one of the most important properties for active pharmaceutical ingredients (API), which can significantly affect their downstream processing. In this study, to better understand the role of...
Collapse
|
4
|
Katava M, Stirnemann G, Pachetti M, Capaccioli S, Paciaroni A, Sterpone F. Specific Interactions and Environment Flexibility Tune Protein Stability under Extreme Crowding. J Phys Chem B 2021; 125:6103-6111. [PMID: 34100611 DOI: 10.1021/acs.jpcb.1c01511] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Macromolecular crowding influences protein mobility and stability in vivo. A precise description of the crowding effect on protein thermal stability requires the estimate of the combined effects of excluded volume, specific protein-environment interactions, as well as the thermal response of the crowders. Here, we explore an ideal model system, the lysozyme protein in powder state, to dissect the factors controlling the melting of the protein under extreme crowding. By deploying state-of-the art molecular simulations, supported by calorimetric experiments, we assess the role of the environment flexibility and of intermolecular electrostatic interactions. In particular, we show that the temperature-dependent flexibility of the macromolecular crowders, along with specific interactions, significantly alleviates the stabilizing contributions of the static volume effect.
Collapse
Affiliation(s)
- Marina Katava
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Guillaume Stirnemann
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Maria Pachetti
- Elettra-Sincrotrone Trieste, S.S. 14 km 163.5, Area Science Park, 34149 Trieste, Italy.,Department of Physics, University of Trieste, Via Valerio 2, 34127 Trieste, Italy
| | - Simone Capaccioli
- Dipartimento di Fisica, Universitá di Pisa, largo Pontecorvo 3, 56127 Pisa, Italy.,CISUP, Centro per l'Integrazione della Strumentazione dell'Università di Pisa, Lungarno Pacinotti 43, I-56127 Pisa, Italy
| | - Alessandro Paciaroni
- Dipartimento di Fisica e Geologia, Universitá di Perugia, via A. Pascoli, 06123 Perugia, Italy
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, IBPC, CNRS UPR9080, Université Paris Diderot, Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| |
Collapse
|
5
|
Quasi-native transition and self-diffusion of proteins in water-glycerol mixture. Biophys Chem 2020; 257:106274. [DOI: 10.1016/j.bpc.2019.106274] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2019] [Revised: 10/10/2019] [Accepted: 10/21/2019] [Indexed: 01/03/2023]
|
6
|
Yani Y, Kanaujia P, Chow PS, Tan RBH. Effect of API-Polymer Miscibility and Interaction on the Stabilization of Amorphous Solid Dispersion: A Molecular Simulation Study. Ind Eng Chem Res 2017. [DOI: 10.1021/acs.iecr.7b03187] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Yin Yani
- Institute of Chemical & Engineering Sciences, A*STAR (Agency for Science, Technology and Research), 1 Pesek Road, Jurong Island 627833, Singapore
| | - Parijat Kanaujia
- Institute of Chemical & Engineering Sciences, A*STAR (Agency for Science, Technology and Research), 1 Pesek Road, Jurong Island 627833, Singapore
| | - Pui Shan Chow
- Institute of Chemical & Engineering Sciences, A*STAR (Agency for Science, Technology and Research), 1 Pesek Road, Jurong Island 627833, Singapore
| | - Reginald B. H. Tan
- Institute of Chemical & Engineering Sciences, A*STAR (Agency for Science, Technology and Research), 1 Pesek Road, Jurong Island 627833, Singapore
- Department of Chemical & Biomolecular Engineering, National University of Singapore, 4 Engineering Drive 4, 117576 Singapore
| |
Collapse
|
7
|
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
|
8
|
Demuth D, Haase N, Malzacher D, Vogel M. Effects of solvent concentration and composition on protein dynamics: 13C MAS NMR studies of elastin in glycerol-water mixtures. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2015; 1854:995-1000. [PMID: 25917596 DOI: 10.1016/j.bbapap.2015.04.012] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/04/2015] [Revised: 03/27/2015] [Accepted: 04/13/2015] [Indexed: 11/25/2022]
Abstract
We use (13)C CP MAS NMR to investigate the dependence of elastin dynamics on the concentration and composition of the solvent at various temperatures. For elastin in pure glycerol, line-shape analysis shows that larger-scale fluctuations of the protein backbone require a minimum glycerol concentration of ~0.6 g/g at ambient temperature, while smaller-scale fluctuations are activated at lower solvation levels of ~0.2 g/g. Immersing elastin in various glycerol-water mixtures, we observe at room temperature that the protein mobility is higher for lower glycerol fractions in the solvent and, thus, lower solvent viscosity. When decreasing the temperature, the elastin spectra approach the line shape for the rigid protein at 245 K for all studied samples, indicating that the protein ceases to be mobile on the experimental time scale of ~10(-5) s. Our findings yield evidence for a strong coupling between elastin fluctuations and solvent dynamics and, hence, such interaction is not restricted to the case of protein-water mixtures. Spectral resolution of different carbon species reveals that the protein-solvent couplings can, however, be different for side chain and backbone units. We discuss these results against the background of the slaving model for protein dynamics.
Collapse
Affiliation(s)
- Dominik Demuth
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstraße 6, 64289 Darmstadt, Germany
| | - Nils Haase
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstraße 6, 64289 Darmstadt, Germany
| | - Daniel Malzacher
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstraße 6, 64289 Darmstadt, Germany
| | - Michael Vogel
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstraße 6, 64289 Darmstadt, Germany.
| |
Collapse
|
9
|
Li D, Liu L, Yu H, Zhai Z, Zhang Y, Guo B, Yang C, Liu B. A molecular simulation study of the protection of insulin bioactive structure by trehalose. J Mol Model 2014; 20:2496. [DOI: 10.1007/s00894-014-2496-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2013] [Accepted: 10/09/2014] [Indexed: 10/24/2022]
|
10
|
Hill JJ, Shalaev EY, Zografi G. The importance of individual protein molecule dynamics in developing and assessing solid state protein preparations. J Pharm Sci 2014; 103:2605-2614. [PMID: 24867196 DOI: 10.1002/jps.24021] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2014] [Revised: 05/05/2014] [Accepted: 05/06/2014] [Indexed: 11/09/2022]
Abstract
Processing protein solutions into the solid state is a common approach for generating stable amorphous protein mixtures that are suitable for long-term storage. Great care is typically given to protecting the protein native structure during the various drying steps that render it into the amorphous solid state. However, many studies illustrate that chemical and physical degradations still occur in spite of this amorphous material having good glassy properties and it being stored at temperatures below its glass transition temperature (Tg). Because of these persistent issues and recent biophysical studies that have refined the debate ascribing meaning to the molecular dynamical transition temperature and Tg of protein molecules, we provide an updated discussion on the impact of assessing and managing localized, individual protein molecule nondiffusive motions in the context of proteins being prepared into bulk amorphous mixtures. Our aim is to bridge the pharmaceutical studies addressing bulk amorphous preparations and their glassy behavior, with the biophysical studies historically focused on the nondiffusive internal protein dynamics and a protein's activity, along with their combined efforts in assessing the impact of solvent hydrogen-bonding networks on local stability. We also provide recommendations for future research efforts in solid-state formulation approaches.
Collapse
Affiliation(s)
- John J Hill
- Department of Bioengineering, University of Washington, Seattle, WA 98195.
| | | | - George Zografi
- School of Pharmacy, University of Wisconsin-Madison, Madison, Wisconsin 53705-2222
| |
Collapse
|
11
|
Vural D, Hong L, Smith JC, Glyde HR. Long-time mean-square displacements in proteins. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:052706. [PMID: 24329295 DOI: 10.1103/physreve.88.052706] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/28/2013] [Indexed: 06/03/2023]
Abstract
We propose a method for obtaining the intrinsic, long-time mean square displacement (MSD) of atoms and molecules in proteins from finite-time molecular dynamics (MD) simulations. Typical data from simulations are limited to times of 1 to 10 ns, and over this time period the calculated MSD continues to increase without a clear limiting value. The proposed method consists of fitting a model to MD simulation-derived values of the incoherent intermediate neutron scattering function, I(inc)(Q,t), for finite times. The infinite-time MSD, <r(2)>, appears as a parameter in the model and is determined by fits of the model to the finite-time I(inc)(Q,t). Specifically, the <r(2)> is defined in the usual way in terms of the Debye-Waller factor as I(Q,t=∞)=exp(-Q(2)<r(2)>/3). The method is illustrated by obtaining the intrinsic MSD <r(2)> of hydrated lysozyme powder (h=0.4 g water/g protein) over a wide temperature range. The intrinsic <r(2)> obtained from data out to 1 and to 10 ns is found to be the same. The intrinsic <r(2)> is approximately twice the value of the MSD that is reached in simulations after times of 1 ns which correspond to those observed using neutron instruments that have an energy resolution width of 1 μeV.
Collapse
Affiliation(s)
- Derya Vural
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716-2570, USA
| | - Liang Hong
- Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA and Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Jeremy C Smith
- Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA and Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Henry R Glyde
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716-2570, USA
| |
Collapse
|
12
|
He X, Lopes PEM, MacKerell AD. Polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator. Biopolymers 2013; 99:724-38. [PMID: 23703219 PMCID: PMC3902549 DOI: 10.1002/bip.22286] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2013] [Accepted: 05/05/2013] [Indexed: 01/12/2023]
Abstract
A polarizable empirical force field for acyclic polyalcohols based on the classical Drude oscillator is presented. The model is optimized with an emphasis on the transferability of the developed parameters among molecules of different sizes in this series and on the condensed-phase properties validated against experimental data. The importance of the explicit treatment of electronic polarizability in empirical force fields is demonstrated in the cases of this series of molecules with vicinal hydroxyl groups that can form cooperative intra- and intermolecular hydrogen bonds. Compared to the CHARMM additive force field, improved treatment of the electrostatic interactions avoids overestimation of the gas-phase dipole moments resulting in significant improvement in the treatment of the conformational energies and leads to the correct balance of intra- and intermolecular hydrogen bonding of glycerol as evidenced by calculated heat of vaporization being in excellent agreement with experiment. Computed condensed phase data, including crystal lattice parameters and volumes and densities of aqueous solutions are in better agreement with experimental data as compared to the corresponding additive model. Such improvements are anticipated to significantly improve the treatment of polymers in general, including biological macromolecules.
Collapse
Affiliation(s)
- Xibing He
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
| | - Pedro E. M. Lopes
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, MD 21201
| |
Collapse
|
13
|
GhattyVenkataKrishna PK, Carri GA. The effect of complex solvents on the structure and dynamics of protein solutions: The case of Lysozyme in trehalose/water mixtures. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2013; 36:14. [PMID: 23404569 DOI: 10.1140/epje/i2013-13014-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2012] [Revised: 10/15/2012] [Accepted: 01/17/2013] [Indexed: 06/01/2023]
Abstract
We present a Molecular Dynamics simulation study of the effect of trehalose concentration on the structure and dynamics of individual proteins immersed in trehalose/water mixtures. Hen egg-white Lysozyme is used in this study and trehalose concentrations of 0%, 10%, 20%, 30% and 100% by weight are explored. Surprisingly, we have found that changes in trehalose concentration do not change the global structural characteristics of the protein as measured by standard quantities like the mean square deviation, radius of gyration, solvent accessible surface area, inertia tensor and asphericity. Only in the limit of pure trehalose these metrics change significantly. Specifically, we found that the protein is compressed by 2% when immersed in pure trehalose. At the amino acid level there is noticeable rearrangement of the surface residues due to the change in polarity of the surrounding environment with the addition of trehalose. From a dynamic perspective, our computation of the Incoherent Intermediate Scattering Function shows that the protein slows down with increasing trehalose concentration; however, this slowdown is not monotonic. Finally, we also report in-depth results for the hydration layer around the protein including its structure, hydrogen-bonding characteristics and dynamic behavior at different length scales.
Collapse
|
14
|
Herbers CR, Sauer D, Vogel M. 2H NMR studies of glycerol dynamics in protein matrices. J Chem Phys 2012; 136:124511. [PMID: 22462878 DOI: 10.1063/1.3697448] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We use (2)H NMR spectroscopy to investigate the rotational motion of glycerol molecules in matrices provided by the connective tissue proteins elastin and collagen. Analyzing spin-lattice relaxation, line-shape properties, and stimulated-echo decays, we determine the rates and geometries of the motion as a function of temperature and composition. It is found that embedding glycerol in an elastin matrix leads to a mild slowdown of glycerol reorientation at low temperatures and glycerol concentrations, while the effect vanishes at ambient temperatures or high solvent content. Furthermore, it is observed that the nonexponential character of the rotational correlation functions is much more prominent in the elastin matrix than in the bulk liquid. Results from spin-lattice relaxation and line shape measurements indicate that, in the mixed systems, the strong nonexponentiality is in large part due to the existence of distributions of correlation times, which are broader on the long-time flank and, hence, more symmetric than in the neat system. Stimulated-echo analysis of slow glycerol dynamics reveals that, when elastin is added, the mechanism for the reorientation crosses over from small-angle jump dynamics to large-angle jump dynamics and the geometry of the motion changes from isotropic to anisotropic. The results are discussed against the background of present and previous findings for glycerol and water dynamics in various protein matrices and compared with observations for other dynamically highly asymmetric mixtures so as to ascertain in which way the viscous freezing of a fast component in the matrix of a slow component differs from the glassy slowdown in neat supercooled liquids.
Collapse
Affiliation(s)
- C R Herbers
- Institut für Physikalische Chemie, Westfälische Wilhelms-Universität Münster, Corrensstraße 28/30, 48149 Münster, Germany
| | | | | |
Collapse
|
15
|
Mamontov E, Chu XQ. Water–protein dynamic coupling and new opportunities for probing it at low to physiological temperatures in aqueous solutions. Phys Chem Chem Phys 2012; 14:11573-88. [DOI: 10.1039/c2cp41443k] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
16
|
Khodadadi S, Roh JH, Kisliuk A, Mamontov E, Tyagi M, Woodson SA, Briber RM, Sokolov AP. Dynamics of biological macromolecules: not a simple slaving by hydration water. Biophys J 2010; 98:1321-6. [PMID: 20371332 DOI: 10.1016/j.bpj.2009.12.4284] [Citation(s) in RCA: 95] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2009] [Revised: 11/24/2009] [Accepted: 12/04/2009] [Indexed: 11/29/2022] Open
Abstract
We studied the dynamics of hydrated tRNA using neutron and dielectric spectroscopy techniques. A comparison of our results with earlier data reveals that the dynamics of hydrated tRNA is slower and varies more strongly with temperature than the dynamics of hydrated proteins. At the same time, tRNA appears to have faster dynamics than DNA. We demonstrate that a similar difference appears in the dynamics of hydration water for these biomolecules. The results and analysis contradict the traditional view of slaved dynamics, which assumes that the dynamics of biological macromolecules just follows the dynamics of hydration water. Our results demonstrate that the dynamics of biological macromolecules and their hydration water depends strongly on the chemical and three-dimensional structures of the biomolecules. We conclude that the whole concept of slaving dynamics should be reconsidered, and that the mutual influence of biomolecules and their hydration water must be taken into account.
Collapse
Affiliation(s)
- S Khodadadi
- Department of Polymer Science, University of Akron, Akron, Ohio, USA
| | | | | | | | | | | | | | | |
Collapse
|
17
|
Uvarov MN, Kulik LV, Dzuba SA. Spin relaxation of fullerene C70 photoexcited triplet in molecular glasses: Evidence for onset of fast orientational motions of molecules in the matrix near 100 K. J Chem Phys 2009; 131:144501. [DOI: 10.1063/1.3244983] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Busselez R, Lefort R, Guendouz M, Frick B, Merdrignac-Conanec O, Morineau D. Molecular dynamics of glycerol and glycerol-trehalose bioprotectant solutions nanoconfined in porous silicon. J Chem Phys 2009; 130:214502. [DOI: 10.1063/1.3147222] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
|
19
|
Busselez R, Lefort R, Ji Q, Affouard F, Morineau D. Molecular dynamics simulation of nanoconfined glycerol. Phys Chem Chem Phys 2009; 11:11127-33. [DOI: 10.1039/b911859d] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
20
|
Khodadadi S, Pawlus S, Roh JH, Garcia Sakai V, Mamontov E, Sokolov AP. The origin of the dynamic transition in proteins. J Chem Phys 2008; 128:195106. [DOI: 10.1063/1.2927871] [Citation(s) in RCA: 146] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Tarek M, Tobias DJ. The role of protein–solvent hydrogen bond dynamics in the structural relaxation of a protein in glycerol versus water. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2008; 37:701-9. [DOI: 10.1007/s00249-008-0324-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2007] [Revised: 03/28/2008] [Accepted: 04/01/2008] [Indexed: 11/28/2022]
|
22
|
Sokolov A, Roh J, Mamontov E, García Sakai V. Role of hydration water in dynamics of biological macromolecules. Chem Phys 2008. [DOI: 10.1016/j.chemphys.2007.07.013] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
23
|
Yongye AB, Foley BL, Woods RJ. On achieving experimental accuracy from molecular dynamics simulations of flexible molecules: aqueous glycerol. J Phys Chem A 2008; 112:2634-9. [PMID: 18311953 PMCID: PMC4201037 DOI: 10.1021/jp710544s] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The rotational isomeric states (RIS) of glycerol at infinite dilution have been characterized in the aqueous phase via a 1 micros conventional molecular dynamics (MD) simulation, a 40 ns enhanced sampling replica exchange molecular dynamics (REMD) simulation, and a reevaluation of the experimental NMR data. The MD and REMD simulations employed the GLYCAM06/AMBER force field with explicit treatment of solvation. The shorter time scale of the REMD sampling method gave rise to RIS and theoretical scalar 3J(HH) coupling constants that were comparable to those from the much longer traditional MD simulation. The 3J(HH) coupling constants computed from the MD methods were in excellent agreement with those observed experimentally. Despite the agreement between the computed and the experimental J-values, there were variations between the rotamer populations computed directly from the MD data and those derived from the experimental NMR data. The experimentally derived populations were determined utilizing limiting J-values from an analysis of NMR data from substituted ethane molecules and may not be completely appropriate for application in more complex molecules, such as glycerol. Here, new limiting J-values have been derived via a combined MD and quantum mechanical approach and were used to decompose the experimental 3J(HH) coupling constants into population distributions for the glycerol RIS.
Collapse
Affiliation(s)
- Austin B. Yongye
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, Georgia 30602
| | - B. Lachele Foley
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, Georgia 30602
| | - Robert J. Woods
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, Georgia 30602
| |
Collapse
|
24
|
Heyes DJ, Sakuma M, Scrutton NS. Laser excitation studies of the product release steps in the catalytic cycle of the light-driven enzyme, protochlorophyllide oxidoreductase. J Biol Chem 2007; 282:32015-20. [PMID: 17848549 DOI: 10.1074/jbc.m706098200] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The latter stages of the catalytic cycle of the light-driven enzyme, protochlorophyllide oxidoreductase, have been investigated using novel laser photoexcitation methods. The formation of the ternary product complex was initiated with a 6-ns laser pulse, which allowed the product release steps to be kinetically accessed for the first time. Subsequent absorbance changes associated with the release of the NADP+ and chlorophyllide products from the enzyme could be followed on a millisecond timescale. This has facilitated a detailed kinetic and thermodynamic characterization for the interconversion of all the various bound and unbound product species. Initially, NADP+ is released from the enzyme in a biphasic process with rate constants of 1210 and 237 s(-1). The rates of both phases show a significant dependence on the viscosity of the solvent and become considerably slower at higher glycerol concentrations. The fast phase of this process exhibits no dependence on NADP+ concentration, suggesting that conformational changes are required prior to NADP+ release. Following NADP+ release, the NADPH rebinds to the enzyme with a maximum rate constant of approximately 72 s(-1). At elevated temperatures (>298 K) chlorophyllide is released from the enzyme to yield the free product with a maximum rate constant of 20 s(-1). The temperature dependencies of the rates of each of these steps were measured, and enthalpies and entropies of activation were calculated using the Eyring equation. A comprehensive kinetic and thermodynamic scheme for these final stages of the reaction mechanism is presented.
Collapse
Affiliation(s)
- Derren J Heyes
- Manchester Interdisciplinary Biocentre, Faculty of Life Sciences, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom.
| | | | | |
Collapse
|
25
|
Abstract
The dynamics of a folded protein is studied in water and glycerol at a series of temperatures below and above their respective dynamical transition. The system is modeled in two distinct states whereby the protein is decoupled from the bulk solvent at low temperatures, and communicates with it through a vicinal layer at physiological temperatures. A linear viscoelastic model elucidates the less-than-expected increase in the relaxation times observed in the backbone dynamics of the protein. The model further explains the increase in the flexibility of the protein once the transition takes place and the differences in the flexibility under the different solvent environments. Coupling between the vicinal layer and the protein fluctuations is necessary to interpret these observations. The vicinal layer is postulated to form once a threshold for the volumetric fluctuations in the protein to accommodate solvents of different sizes is reached. Compensation of entropic-energetic contributions from the protein-coupled vicinal layer quantifies the scaling of the dynamical transition temperatures in various solvents. The protein adapts different conformational routes for organizing the required coupling to a specific solvent, which is achieved by adjusting the amount of conformational jumps in the surface-group dihedrals.
Collapse
Affiliation(s)
- Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey.
| | | | | |
Collapse
|
26
|
Dirama TE, Curtis JE, Carri GA, Sokolov AP. Coupling between lysozyme and trehalose dynamics: microscopic insights from molecular-dynamics simulations. J Chem Phys 2007; 124:034901. [PMID: 16438608 DOI: 10.1063/1.2159471] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have carried out molecular-dynamics simulations on fully flexible all-atom models of the protein lysozyme immersed in trehalose, an effective biopreservative, with the purpose of exploring the nature and extent of the dynamical coupling between them. Our study shows a strong coupling over a wide range of temperatures. We found that the onset of anharmonic behavior was dictated by changes in the dynamics and relaxation processes in the trehalose glass. The physical origin of protein-trehalose coupling was traced to the hydrogen bonds formed at the interface between the protein and the solvent. Moreover, protein-solvent hydrogen bonding was found to control the structural relaxation of the protein. The dynamics of the protein was found to be heterogeneous; the motions of surface and core atoms had different dependencies on temperature and, in addition, the surface atoms were more sensitive to the dynamics of the solvent than the core atoms. From the solvent perspective we found that the dynamics near the protein surface showed an unexpected enhanced mobility compared to the bulk. These results shed some light on the microscopic origins of the dynamical coupling in protein-solvent systems.
Collapse
Affiliation(s)
- Taner E Dirama
- Department of Polymer Science, The University of Akron, Akron, Ohio 44325-3909, USA
| | | | | | | |
Collapse
|
27
|
Hunt NT, Kattner L, Shanks RP, Wynne K. The dynamics of water-protein interaction studied by ultrafast optical Kerr-effect spectroscopy. J Am Chem Soc 2007; 129:3168-72. [PMID: 17315992 DOI: 10.1021/ja066289n] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Changes in the ultrafast dynamics and terahertz Raman spectrum accompanying a helix-to-coil transition of a homo-polypeptide have been observed for the first time. Formation of the alpha-helix is associated with a shift to lower frequency of a broad Raman band attributable to solvent-peptide intermolecular hydrogen bonding. This band facilitates direct spectroscopic observation of so-called hydration water near a peptide and yields the first quantitative estimate of the time scale of the ultrafast dynamics in the solvation shell, which range from 0.18 to 0.33 ps (185-100 cm(-1)) depending on the secondary structure of the peptide. Such fast motions of solvent molecules have been referred to as the "lubricant of life" and are thought to play key roles in determining structure and activity of proteins.
Collapse
Affiliation(s)
- Neil T Hunt
- Department of Physics, SUPA, University of Strathclyde, Glasgow G4 0NG, Scotland, United Kingdom.
| | | | | | | |
Collapse
|
28
|
Curtis JE, Dirama TE, Carri GA, Tobias DJ. Inertial Suppression of Protein Dynamics in a Binary Glycerol−Trehalose Glass. J Phys Chem B 2006; 110:22953-6. [PMID: 17107124 DOI: 10.1021/jp0615499] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The traditional approach used to predict the ability of a glassy matrix to maximally preserve the activity of a protein solute is the glass transition temperature (T(g)) of the glass. Recently it has been shown that the addition of a low T(g) diluent (glycerol) can rigidify the structure of a high T(g) glassy matrix in binary glycerol-trehalose glasses. The optimal density of glycerol in trehalose minimizes the average mean square displacements of non-exchangeable protons in the glass samples. The amount of glycerol added to a trehalose glass coincides with the maximal recovery of biological activity in a separate study using similar binary glass samples. In this study, we use molecular dynamics (MD) simulations to investigate the dynamics of a hydrated protein encased in glycerol, unary trehalose and binary glycerol-trehalose glasses. We have found that we are able to reproduce the rigidification of the glycerol-trehalose glassy matrix and that there is a direct correlation between bulk glass dynamics and the extent of atomic fluctuation of protein atoms. The detailed microscopic picture that emerges is that protein dynamics are suppressed mainly by inertia of the bulk glass and to a lesser extent specific interactions at the protein-solvent interface. Thus, the inertia of the glassy matrix may be an influential factor in the determination of pharmaceutically relevant formulations.
Collapse
|
29
|
Roh JH, Curtis JE, Azzam S, Novikov VN, Peral I, Chowdhuri Z, Gregory RB, Sokolov AP. Influence of hydration on the dynamics of lysozyme. Biophys J 2006; 91:2573-88. [PMID: 16844746 PMCID: PMC1562387 DOI: 10.1529/biophysj.106.082214] [Citation(s) in RCA: 181] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Quasielastic neutron and light-scattering techniques along with molecular dynamics simulations were employed to study the influence of hydration on the internal dynamics of lysozyme. We identified three major relaxation processes that contribute to the observed dynamics in the picosecond to nanosecond time range: 1), fluctuations of methyl groups; 2), fast picosecond relaxation; and 3), a slow relaxation process. A low-temperature onset of anharmonicity at T approximately 100 K is ascribed to methyl-group dynamics that is not sensitive to hydration level. The increase of hydration level seems to first increase the fast relaxation process and then activate the slow relaxation process at h approximately 0.2. The quasielastic scattering intensity associated with the slow process increases sharply with an increase of hydration to above h approximately 0.2. Activation of the slow process is responsible for the dynamical transition at T approximately 200 K. The dependence of the slow process on hydration correlates with the hydration dependence of the enzymatic activity of lysozyme, whereas the dependence of the fast process seems to correlate with the hydration dependence of hydrogen exchange of lysozyme.
Collapse
Affiliation(s)
- J H Roh
- Department of Polymer Science, The University of Akron, Akron, Ohio 44325-3909, USA
| | | | | | | | | | | | | | | |
Collapse
|