Molecular dynamics simulation of carboxy-myoglobin embedded in a trehalose-water matrix

A Long Journey into the Investigation of the Structure–Dynamics–Function Paradigm in Proteins through the Activities of the Palermo Biophysics Group

An overview of the biophysics activity at the Department of Physics and Chemistry Emilio Segrè of the University of Palermo is given. For forty years, the focus of the research has been on the protein structure–dynamics–function paradigm, with the aim of understanding the molecular basis of the relevant mechanisms and the key role of solvent. At least three research lines are identified; the main results obtained in collaboration with other groups in Italy and abroad are presented. This review is dedicated to the memory of Professors Massimo Ugo Palma, Maria Beatrice Palma Vittorelli, and Lorenzo Cordone, which were the founders of the Palermo School of Biophysics. We all have been, directly or indirectly, their pupils; we miss their enthusiasm for scientific research, their deep physical insights, their suggestions, their strict but always constructive criticisms, and, most of all, their friendship. This paper is dedicated also to the memory of Prof. Hans Frauenfelder, whose pioneering works on nonexponential rebinding kinetics, protein substates, and energy landscape have inspired a large part of our work in the field of protein dynamics.

Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (ν_CO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A₀, A₁, and A₃) and (ii) the modeling and the interpretation of the ν_CO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A₀ and A₁/A₃) spectroscopic regions; (ii) the characteristic bimodal shape of the A₁/A₃ spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the ν_CO energy relaxation, characterizing its kinetic efficiency.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Properties of Aqueous Trehalose Mixtures: Glass Transition and Hydrogen Bonding

Trehalose is a naturally occurring disaccharide known to remarkably stabilize biomacromolecules in the biologically active state. The stabilizing effect is typically observed over a large concentration range and affects many macromolecules including proteins, lipids, and DNA. Of special interest is the transition from aqueous solution to the dense and highly concentrated glassy state of trehalose that has been implicated in bioadaptation of different organisms toward desiccation stress. Although several mechanisms have been suggested to link the structure of the low water content glass with its action as an exceptional stabilizer, studies are ongoing to resolve which are most pertinent. Specifically, the role that hydrogen bonding plays in the formation of the glass is not well resolved. Here we model aqueous trehalose mixtures over a wide concentration range, using molecular dynamics simulations with two available force fields. Both force fields indicate glass transition temperatures and osmotic pressures that are close to experimental values, particularly at high trehalose contents. We develop and employ a methodology that allows us to analyze the thermodynamics of hydrogen bonds in simulations at different water contents and temperatures. Remarkably, this analysis is able to link the liquid to glass transition with changes in hydrogen bond characteristics. Most notably, the onset of the glassy state can be quantitatively related to the transition from weakly to strongly correlated hydrogen bonds. Our findings should help resolve the properties of the glass and the mechanisms of its formation in the presence of added macromolecules.

More than a Confinement: “Soft” and “Hard” Enzyme Entrapment Modulates Biological Catalyst Function

Catalysis makes chemical and biochemical reactions kinetically accessible. From a technological point of view, organic, inorganic, and biochemical catalysis is relevant for several applications, from industrial synthesis to biomedical, material, and food sciences. A heterogeneous catalyst, i.e., a catalyst confined in a different phase with respect to the reagents’ phase, requires either its physical confinement in an immobilization matrix or its physical adsorption on a surface. In this review, we will focus on the immobilization of biological catalysts, i.e., enzymes, by comparing hard and soft immobilization matrices and their effect on the modulation of the catalysts’ function. Indeed, unlike smaller molecules, the catalytic activity of protein catalysts depends on their structure, conformation, local environment, and dynamics, properties that can be strongly affected by the immobilization matrices, which, therefore, not only provide physical confinement, but also modulate catalysis.

Exploring Dynamics and Structure of Biomolecules, Cryoprotectants, and Water Using Molecular Dynamics Simulations: Implications for Biostabilization and Biopreservation

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.

The effects of pressure on the energy landscape of proteins

Protein dynamics is characterized by fluctuations among different conformational substates, i.e. the different minima of their energy landscape. At temperatures above ~200 K, these fluctuations lead to a steep increase in the thermal dependence of all dynamical properties, phenomenon known as Protein Dynamical Transition. In spite of the intense studies, little is known about the effects of pressure on these processes, investigated mostly near room temperature. We studied by neutron scattering the dynamics of myoglobin in a wide temperature and pressure range. Our results show that high pressure reduces protein motions, but does not affect the onset temperature for the Protein Dynamical Transition, indicating that the energy differences and barriers among conformational substates do not change with pressure. Instead, high pressure values strongly reduce the average structural differences between the accessible conformational substates, thus increasing the roughness of the free energy landscape of the system.

Investigating Structure and Dynamics of Proteins in Amorphous Phases Using Neutron Scattering

In order to increase shelf life and minimize aggregation during storage, many biotherapeutic drugs are formulated and stored as either frozen solutions or lyophilized powders. However, characterizing amorphous solids can be challenging with the commonly available set of biophysical measurements used for proteins in liquid solutions. Therefore, some questions remain regarding the structure of the active pharmaceutical ingredient during freezing and drying of the drug product and the molecular role of excipients. Neutron scattering is a powerful technique to study structure and dynamics of a variety of systems in both solid and liquid phases. Moreover, neutron scattering experiments can generally be correlated with theory and molecular simulations to analyze experimental data. In this article, we focus on the use of neutron techniques to address problems of biotechnological interest. We describe the use of small-angle neutron scattering to study the solution structure of biological molecules and the packing arrangement in amorphous phases, that is, frozen glasses and freeze-dried protein powders. In addition, we discuss the use of neutron spectroscopy to measure the dynamics of glassy systems at different time and length scales. Overall, we expect that the present article will guide and prompt the use of neutron scattering to provide unique insights on many of the outstanding questions in biotechnology.

Structural and dynamical characteristics of trehalose and sucrose matrices at different hydration levels as probed by FTIR and high-field EPR

Some organisms can survive complete dehydration and high temperatures by adopting an anhydrobiotic state in which the intracellular medium contains large amounts of disaccharides, particularly trehalose and sucrose. Trehalose is most effective also in protecting isolated in vitro biostructures. In an attempt to clarify the molecular mechanisms of disaccharide bioprotection, we compared the structure and dynamics of sucrose and trehalose matrices at different hydration levels by means of high-field W-band EPR and FTIR spectroscopy. The hydration state of the samples was characterized by FTIR spectroscopy and the structural organization was probed by EPR using a nitroxide radical dissolved in the respective matrices. Analysis of the EPR spectra showed that the structure and dynamics of the dehydrated matrices as well as their evolution upon re-hydration differ substantially between trehalose and sucrose. The dehydrated trehalose matrix is homogeneous in terms of distribution of the residual water and spin-probe molecules. In contrast, dehydrated sucrose forms a heterogeneous matrix. It is comprised of sucrose polycrystalline clusters and several bulk water domains. The amorphous form was found only in 30% (volume) of the sucrose matrix. Re-hydration leads to a structural homogenization of the sucrose matrix, whilst in the trehalose matrix several domains develop differing in the local water/radical content and radical mobility. The molecular model of the matrices provides an explanation for the different protein-matrix dynamical coupling observed in dried ternary sucrose and trehalose matrices, and accounts for the superior efficacy of trehalose as a bioprotectant. Furthermore, for bacterial photosynthetic reaction centers it is shown that at low water content the protein-matrix coupling is modulated by the sugar/protein molar ratio in sucrose matrices only. This effect is suggested to be related to the preference for sucrose, rather than trehalose, as a bioprotective disaccharide in some anhydrobiotic organisms.

Proteins in amorphous saccharide matrices: structural and dynamical insights on bioprotection

Bioprotection by sugars, and in particular trehalose peculiarity, is a relevant topic due to the implications in several fields. The underlying mechanisms are not yet clearly elucidated, and remain the focus of current investigations. Here we revisit data obtained at our lab on binary sugar/water and ternary protein/sugar/water systems, in wide ranges of water content and temperature, in the light of the current literature. The data here discussed come from complementary techniques (Infrared Spectroscopy, Molecular Dynamics simulations, Small Angle X-ray Scattering and Calorimetry), which provided a consistent description of the bioprotection by sugars from the atomistic to the macroscopic level. We present a picture, which suggests that protein bioprotection can be explained in terms of a strong coupling of the biomolecule surface to the matrix via extended hydrogen-bond networks, whose properties are defined by all components of the systems, and are strongly dependent on water content. Furthermore, the data show how carbohydrates having similar hydrogen-bonding capabilities exhibit different efficiency in preserving biostructures.

How strongly does trehalose interact with lysozyme in the solid state? Insights from molecular dynamics simulation and inelastic neutron scattering

Therapeutic proteins are usually conserved in glassy matrixes composed of stabilizing excipients and a small amount of water, which both control their long-term stability, and thus their potential use in medical treatments. To shed some light on the protein-matrix interactions in such systems, we performed molecular dynamics (MD) simulations on matrixes of (i) the model globular protein lysozyme (L), (ii) the well-known bioprotectant trehalose (T), and (iii) the 1:1 (in weight) lysozyme/trehalose mixture (LT), at hydration levels h of 0.0, 0.075, and 0.15 (in g of water/g of protein or sugar). We also supplemented these simulations with complementary inelastic neutron scattering (INS) experiments on the L, T, and LT lyophilized (freeze-dried) samples. The densities and free volume distributions indicate that trehalose improves the molecular packing of the LT glass with respect to the L one. Accordingly, the low-frequency vibrational densities of states (VDOS) and the mean square displacements (MSDs) of lysozyme reveal that it is less flexible-and thus less likely to unfold-in the presence of trehalose. Furthermore, at low contents (h = 0.075), water systematically stiffens the vibrational motions of lysozyme and trehalose, whereas it increases their MSDs on the nanosecond (ns) time scale. This stems from the hydrogen bonds (HBs) that lysozyme and trehalose form with water, which, interestingly, are stronger than the ones they form with each other but which, nonetheless, relax faster on the ns time scale, given the larger mobility of water. Moreover, lysozyme interacts preferentially with water in the hydrated LT mixtures, and trehalose appears to slow down significantly the relaxation of lysozyme-water HBs. Overall, our results suggest that the stabilizing efficiency of trehalose arises from its ability to (i) increase the number of HBs formed by proteins in the dry state and (ii) make the HBs formed by water with proteins stable on long (>ns) time scales.

SAXS study on myoglobin embedded in amorphous saccharide matrices

We report on Small Angle X-ray Scattering (SAXS) measurements performed on samples of carboxy-myoglobin and met-myoglobin embedded in low hydrated matrices of four different saccharides (trehalose, sucrose, maltose and lactose). Results confirm the already reported occurrence of inhomogeneities, which are not peculiar of trehalose samples, but appear also in maltose and lactose, and in some cases also sucrose, being dependent on the sample hydration and on the presence of sodium dithionite. This behaviour confirms our previous interpretation about the nature of the inhomogeneities, and prompt it as a possible general behaviour for highly concentrated sugar matrices.

Protein thermal denaturation and matrix glass transition in different protein-trehalose-water systems

Biopreservation by saccharides is a widely studied issue due to its scientific and technological importance; in particular, ternary amorphous protein-saccharide-water systems are extensively exploited to model the characteristics of the in vivo biopreservation process. We present here a differential scanning calorimetry (DSC) study on amorphous trehalose-water systems with embedded different proteins (myoglobin, lysozyme, BSA, hemoglobin), which differ for charge, surface, and volume properties. In our study, the protein/trehalose molar ratio is kept constant at 1/40, while the water/sugar molar ratio is varied between 2 and 300; results are compared with those obtained for binary trehalose-water systems. DSC upscans offer the possibility of investigating, in the same measurement, the thermodynamic properties of the matrix (glass transition, T(g)) and the functional properties of the encapsulated protein (thermal denaturation, T(den)). At high-to-intermediate hydration, the presence of the proteins increases the glass transition temperature of the encapsulating matrix. The effect mainly depends on size properties, and it can be ascribed to confinement exerted by the protein on the trehalose-water solvent. Conversely, at low hydration, lower T(g) values are measured in the presence of proteins: the lack of water promotes sugar-protein interactions, thus weakening the confinement effect and softening the matrix with respect to the binary system. A parallel T(den) increase is also observed; remarkably, this stabilization can reach ∼70 K at low hydration, a finding potentially of high biotechnological relevance. A linear relationship between T(g) and T(den) is also observed, in line with previous results; this finding suggests that collective water-trehalose interactions, responsible for the glass transition, also influence the protein denaturation.

Myoglobin embedded in saccharide amorphous matrices: water-dependent domains evidenced by small angle X-ray scattering

We report Small Angle X-ray Scattering (SAXS) measurements performed on samples of carboxy-myoglobin (MbCO) embedded in low-water trehalose glasses. Results showed that, in such samples, "low-protein" trehalose-water domains are present, surrounded by a protein-trehalose-water background; such finding is supported by Infrared Spectroscopy (FTIR) measurements. These domains, which do not appear in the absence of the protein and in analogous sucrose systems, preferentially incorporate the incoming water at the onset of rehydration, and disappear following large hydration. This observation suggests that, in organisms under anhydrobiosis, analogous domains could play a buffering role against the daily variations of the atmospheric moisture. The reported results are rationalized by assuming sizably different protein-matrix coupling in trehalose with respect to sucrose, analogous to the one suggested for the photosynthetic reaction centre from Rhodobacter sphaeroides (F. Francia et al., J. Am. Chem. Soc., 2008, 130, 10240-10246).

Thermal denaturation of myoglobin in water--disaccharide matrixes: relation with the glass transition of the system

Proteins embedded in glassy saccharide systems are protected against adverse environmental conditions [Crowe et al. Annu. Rev. Physiol. 1998, 60, 73-103]. To further characterize this process, we studied the relationship between the glass transition temperature of the protein-containing saccharide system (T(g)) and the temperature of thermal denaturation of the embedded protein (T(den)). To this end, we studied by differential scanning calorimetry the thermal denaturation of ferric myoglobin in water/disaccharide mixtures containing nonreducing (trehalose, sucrose) or reducing (maltose, lactose) disaccharides. All the samples studied are, at room temperature, liquid systems whose viscosity varies from very low to very large values, depending on the water content. At a high water/saccharide mole ratio, homogeneous glass formation does not occur; regions of glass form, whose T(g) does not vary by varying the saccharide content, and the disaccharide barely affects the myoglobin denaturation temperature. At a suitably low water/saccharide mole ratio, by lowering the temperature, the systems undergo transition to the glassy state whose T(g) is determined by the water content; the Gordon-Taylor relationship between T(g) and the water/disaccharide mole ratio is obeyed; and T(den) increases by decreasing the hydration regardless of the disaccharide, such effect being entropy-driven. The presence of the protein was found to lower the T(g). Furthermore, for nonreducing disaccharides, plots of T(den) vs T(g) give linear correlations, whereas for reducing disaccharides, data exhibit an erratic behavior below a critical water/disaccharide ratio. We ascribe this behavior to the likelihood that in the latter samples, proteins have undergone Maillard reaction before thermal denaturation.

Mixing properties of lyophilized protein systems: a spectroscopic and calorimetric study

The purpose of this study was to investigate the solid-state properties of lyophilized formulations of protein (ribonuclease A) containing sucrose or trehalose across a wide range of compositions, both in the presence or absence of hydroxyethylstarch (HES). Infrared spectroscopy reveals that the protein forms hydrogen bonds to sugars (sucrose or trehalose) as water is removed from the sample. The strength and/or number of hydrogen bonds in dried samples increase as the weight fraction of sugar increases. Significant deviations of glass transition temperatures (T(g)'s) from those predicted by free volume theory are seen in both protein-sugar systems. The behavior can be explained by formation of protein-sugar hydrogen bonds at the expense of self-interactions between the sugars. Attractive interactions between lyophilized ribonuclease A and HES were detected spectroscopically and from thermodynamic analysis of T(g) values, contrary to the view that HES is sterically hindered from interacting with the protein surface. Sucrose-HES interactions were much less favorable than trehalose-HES interactions, suggesting that phase separation in sugar/HES/protein mixtures would be more likely in the presence of sucrose than trehalose. Finally, the thermodynamics of mixing were investigated using differential scanning calorimetry (DSC) providing some of the first data for such solid protein sugar formulations with and without HES. In nearly all samples, positive excess enthalpy, excess entropy and excess free energy were observed, with the excess free energy being greater for samples containing sucrose rather than trehalose. Analysis of Flory-Huggins chi parameters suggests that phase separation between protein and excipients may be thermodynamically favored in these dried solid preparations.

Biomolecular motion characterization by a self-distribution-function procedure in elastic incoherent neutron scattering

In the present paper we first focus on the role of the instrumental resolution in elastic incoherent neutron scattering (EINS) where the connection between the self-distribution function (SDF) and the measured EINS intensity profile is highlighted. Second we show how the SDF procedure, previously introduced, allows both the total and the partial mean-square displacement evaluations through the total and the partial SDFs. Finally, we compare the SDF and the Gaussian procedures, by applying the two approaches to EINS data collected, by the IN13 backscattering spectrometer (Institute Laue-Langevin, Grenoble), on aqueous mixtures of two homologous disaccharides, i.e., sucrose and trehalose, and on myoglobin.

Effect of trehalose on a phospholipid membrane under mechanical stress

Explicit solvent molecular dynamics simulations were used to investigate at atomic resolution the effect of trehalose on a hydrated phospholipid bilayer under mechanical stress (stretching forces imposed in the form of negative lateral pressure). Simulations were performed in the absence or presence of trehalose at 325 K, and with different values for negative lateral pressure. In the concentration regime (2 molal) and range of lateral pressures (1 to -250 bar) investigated, trehalose was found to interact directly with the membrane, partially replacing water molecules in the formation of hydrogen bonds with the lipid headgroups. Similar to previous findings in the context of thermal stress, the number, degree of bridging, and reaching depth of these hydrogen bonds increased with the magnitude of perturbation. However, at the concentration considered, trehalose was not sufficient to preserve the integrity of the membrane structure and to prevent its extreme elongation (and possible disruption) under the effect of stretching forces.

Trehalose and trehalose-based polymers for environmentally benign, biocompatible and bioactive materials

Trehalose is a non-reducing disaccharide that is found in many organisms but not in mammals. This sugar plays important roles in cryptobiosis of selaginella mosses, tardigrades (water bears), and other animals which revive with water from a state of suspended animation induced by desiccation. The interesting properties of trehalose are due to its unique symmetrical low-energy structure, wherein two glucose units are bonded face-to-face by 1→1-glucoside links. The Hayashibara Co. Ltd., is credited for developing an inexpensive, environmentally benign and industrial-scale process for the enzymatic conversion of α-1,4-linked polyhexoses to α,α-d-trehalose, which made it easy to explore novel food, industrial, and medicinal uses for trehalose and its derivatives. Trehalose-chemistry is a relatively new and emerging field, and polymers of trehalose derivatives appear environmentally benign, biocompatible, and biodegradable. The discriminating properties of trehalose are attributed to its structure, symmetry, solubility, kinetic and thermodynamic stability and versatility. While syntheses of trehalose-based polymer networks can be straightforward, syntheses and characterization of well defined linear polymers with tailored properties using trehalose-based monomers is challenging, and typically involves protection and deprotection of hydroxyl groups to attain desired structural, morphological, biological, and physical and chemical properties in the resulting products. In this review, we will overview known literature on trehalose’s fascinating involvement in cryptobiology; highlight its applications in many fields; and then discuss methods we used to prepare new trehalose-based monomers and polymers and explain their properties.

Dynamics of C-phycocyanin in various deuterated trehalose/water environments measured by quasielastic and elastic neutron scattering

The molecular understanding of protein stabilization by the disaccharide trehalose in extreme temperature or hydration conditions is still debated. In the present study, we investigated the role of trehalose on the dynamics of the protein C-phycocyanin (C-PC) by neutron scattering. To single out the motions of C-PC hydrogen (H) atoms in various trehalose/water environments, measurements were performed in deuterated trehalose and heavy water (D2O). We report that trehalose decreases the internal C-PC dynamics, as shown by a reduced diffusion coefficient of protein H atoms. By fitting the Elastic Incoherent Structure Factor--which gives access to the "geometry" of the internal proton motions--with the model of diffusion inside a sphere, we found that the presence of trehalose induces a significantly higher proportion of immobile C-PC hydrogens. We investigated, by elastic neutron scattering, the mean square displacements (MSDs) of deuterated trehalose/D2O-embedded C-PC as a function of temperature in the range of 40-318 K. Between 40 and approximately 225 K, harmonic MSDs of C-PC are slightly smaller in samples containing trehalose. Above a transition temperature of approximately 225 K, we observed anharmonic motions in all trehalose/water-coated C-PC samples. In the hydrated samples, MSDs are not significantly changed by addition of 15% trehalose but are slightly reduced by 30% trehalose. In opposition, no dynamical transition was detected in dry trehalose-embedded C-PC, whose hydrogen motions remain harmonic up to 318 K. These results suggest that a role of trehalose would be to stabilize proteins by inhibiting some fluctuations at the origin of protein unfolding and denaturation.

Ligand recombination and a hierarchy of solvent slaved dynamics: the origin of kinetic phases in hemeproteins

Ligand recombination studies play a central role both for characterizing different hemeproteins and their conformational states but also for probing fundamental biophysical processes. Consequently, there is great importance to providing a foundation from which one can understand the physical processes that give rise to and modulate the large range of kinetic patterns associated with ligand recombination in myoglobins and hemoglobins. In this work, an overview of cryogenic and solution phase recombination phenomena for COMb is first reviewed and then a new paradigm is presented for analyzing the temperature and viscosity dependent features of kinetic traces in terms of multiple phases that reflect which tier(s) of solvent slaved protein dynamics is (are) operative on the photoproduct population during the time course of the measurement. This approach allows for facile inclusion of both ligand diffusion among accessible cavities and conformational relaxation effects. The concepts are illustrated using kinetic traces and MEM populations derived from the CO recombination process for wild type and mutant myoglobins either in sol-gel matrices bathed in glycerol or in trehalose-derived glassy matrices.

GFP-mut2 proteins in trehalose-water matrixes: spatially heterogeneous protein-water-sugar structures

We report investigations on the properties of nanoenvironments around single-GFP-mut2 proteins in trehalose-water matrixes. Single-GFPmut2 molecules embedded in thin trehalose-water films were characterized in terms of their fluorescence brightness, bleaching dynamics, excited state lifetime, and fluorescence polarization. For each property, sets of approximately 100-150 single molecules have been investigated as a function of trehalose content and hydration. Three distinct and interconverting families of proteins have been found which differ widely in terms of bleaching dynamics, brightness, and fluorescence polarization, whose relative populations sizably depend on sample hydration. The reported results evidence the simultaneous presence of different protein-trehalose-water nanostructures whose rigidity increases by lowering the sample hydration. Such spatial inhomogeneity is in line with the well-known heterogeneous dynamics in supercooled fluids and in nonsolid carbohydrate glasses and gives a pictorial representation of the sharp, sudden reorganization of the above structures after uptake <==>release of water molecules.

A comparative study of carboxy myoglobin in saccharide-water systems by molecular dynamics simulation

Results from room-temperature molecular dynamics simulation on a system containing carboxy-myoglobin, water, and maltose molecules are reported. Protein atomic fluctuations, protein-solvent and solvent-solvent hydrogen bonding have been analyzed and compared to the ones in trehalose-water and sucrose-water systems (Proteins 2005, 59, 291-302). Results help in rationalizing, at a molecular level, the effects of homologues disaccharides on protein structure/dynamics experimentally observed. Furthermore, the effectiveness of disaccharides in bioprotection in terms of peculiar protein-matrix coupling is also discussed.

Cytochrome C in a dry trehalose matrix: structural and dynamical effects probed by x-ray absorption spectroscopy

We report on the structure and dynamics of the Fe ligand cluster of reduced horse heart cytochrome c in solution, in a dried polyvinyl alcohol (PVA) film, and in two trehalose matrices characterized by different contents of residual water. The effect of the solvent/matrix environment was studied at room temperature using Fe K-edge x-ray absorption fine structure (XAFS) spectroscopy. XAFS data were analyzed by combining ab initio simulations and multi-parameter fitting in an attempt to disentangle structural from disorder parameters. Essentially the same structural and disorder parameters account adequately for the XAFS spectra measured in solution, both in the absence and in the presence of glycerol, and in the PVA film, showing that this polymer interacts weakly with the embedded protein. Instead, incorporation in trehalose leads to severe structural changes, more prominent in the more dried matrix, consisting of 1), an increase up to 0.2 A of the distance between Fe and the imidazole N atom of the coordinating histidine residue and 2), an elongation up to 0.16 A of the distance between Fe and the fourth-shell C atoms of the heme pyrrolic units. These structural distortions are accompanied by a substantial decrease of the relative mean-square displacements of the first ligands. In the extensively dried trehalose matrix, extremely low values of the Debye Waller factors are obtained for the pyrrolic and for the imidazole N atoms. This finding is interpreted as reflecting a drastic hindering in the relative motions of the Fe ligand cluster atoms and an impressive decrease in the static disorder of the local Fe structure. It appears, therefore, that the dried trehalose matrix dramatically perturbs the energy landscape of cytochrome c, giving rise, at the level of local structure, to well-resolved structural distortions and restricting the ensemble of accessible conformational substates.

Inertial Suppression of Protein Dynamics in a Binary Glycerol−Trehalose Glass

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.

Infrared absorption study of the heme pocket dynamics of carbonmonoxyheme proteins

The temperature dependencies of the infrared absorption CO bands of carboxy complexes of horseradish peroxidase (HRP(CO)) in glycerol/water mixture at pH 6.0 and 9.3 are interpreted using the theory of optical absorption bandshape. The bands' anharmonic behavior is explained assuming that there is a higher-energy set of conformational substates (CSS(h)), which are populated upon heating and correspond to the protein substates with disordered water molecules in the heme pocket. Analysis of the second moments of the CO bands of the carboxy complexes of myoglobin (Mb(CO)) and hemoglobin (Hb(CO)), and of HRP(CO) with benzohydroxamic acid (HRP(CO)+BHA), shows that the low energy CSS(h) exists also in the open conformation of Mb(CO), where the heme pocket is spacious enough to accommodate a water molecule. In the HRP(CO)+BHA and closed conformations of Mb(CO) and Hb(CO), the heme pocket is packed with BHA and different amino acids, the CSS(h) has much higher energy and is hardly populated even at the highest temperatures. Therefore only motions of these amino acids contribute to the band broadening. These motions are linked to the protein surface and frozen in the glassy matrix, whereas in the liquid solvent they are harmonic. Thus the second moment of the CO band is temperature-independent in glass and is proportional to the temperature in liquid. The temperature dependence of the second moment of the CO peak of HRP(CO) in the trehalose glass exhibits linear coupling to an oscillator. This oscillator can be a moving water molecule locked in the heme pocket in the whole interval of temperatures or a trehalose molecule located in the heme pocket.

Role of solvent on protein-matrix coupling in MbCO embedded in water-saccharide systems: a Fourier transform infrared spectroscopy study

Embedding protein in sugar systems of low water content enables one to investigate the protein dynamic-structure function in matrixes whose rigidity is modulated by varying the content of residual water. Accordingly, studying the dynamics and structure thermal evolution of a protein in sugar systems of different hydration constitutes a tool for disentangling solvent rigidity from temperature effects. Furthermore, studies performed using different sugars may give information on how the detailed composition of the surrounding solvent affects the internal protein dynamics and structural evolution. In this work, we compare Fourier transform infrared spectroscopy measurements (300-20 K) on MbCO embedded in trehalose, sucrose, maltose, raffinose, and glucose matrixes of different water content. At all the water contents investigated, the protein-solvent coupling was tighter in trehalose than in the other sugars, thus suggesting a molecular basis for the trehalose peculiarity. These results are in line with the observation that protein-matrix phase separation takes place in lysozyme-lactose, whereas it is absent in lysozyme-trehalose systems; indeed, these behaviors may respectively be due to the lack or presence of suitable water-mediated hydrogen-bond networks, which match the protein surface to the surroundings. The above processes might be at the basis of pattern recognition in crowded living systems; indeed, hydration shells structural and dynamic matching is first needed for successful come together of interacting biomolecules.

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.

The Influence of Aqueous versus Glassy Solvents on Protein Dynamics: Vibrational Echo Experiments and Molecular Dynamics Simulations

Spectrally resolved infrared stimulated vibrational echo measurements are used to measure the vibrational dephasing of the CO stretching mode of carbonmonoxy-hemoglobin (HbCO), a myoglobin mutant (H64V), and a bacterial cytochrome c(552) mutant (Ht-M61A) in aqueous solution and trehalose glasses. The vibrational dephasing of the heme-bound CO is significantly slower for all three proteins embedded in trehalose glasses compared to that of aqueous protein solutions. All three proteins exhibit persistent but notably slower spectral diffusion when the protein surface is fixed by the glassy solvent. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to reveal that the structural dynamics, as sensed by the CO, of the three proteins in trehalose and aqueous solution are dominated by fast (tens of femtoseconds), motionally narrowed fluctuations. MD simulations of H64V in dynamic and "static" water are presented as models of the aqueous and glassy environments. FFCFs are calculated from the H64V simulations and qualitatively reproduce the important features of the experimentally extracted FFCFs. The suppression of long time scale (picoseconds to tens of picoseconds) frequency fluctuations (spectral diffusion) in the glassy solvent is the result of a damping of atomic displacements throughout the protein structure and is not limited to structural dynamics that occur only at the protein surface. The analysis provides evidence that some dynamics are coupled to the hydration shell of water, supporting the idea that the bioprotection offered by trehalose is due to its ability to immobilize the protein surface through a thin, constrained layer of water.

Internal dynamics and protein–matrix coupling in trehalose-coated proteins

We review recent studies on the role played by non-liquid, water-containing matrices on the dynamics and structure of embedded proteins. Two proteins were studied, in water-trehalose matrices: a water-soluble protein (carboxy derivative of horse heart myoglobin) and a membrane protein (reaction centre from Rhodobacter sphaeroides). Several experimental techniques were used: Mossbauer spectroscopy, elastic neutron scattering, FTIR spectroscopy, CO recombination after flash photolysis in carboxy-myoglobin, kinetic optical absorption spectroscopy following pulsed and continuous photoexcitation in Q(B) containing or Q(B) deprived reaction centre from R. sphaeroides. Experimental results, together with the outcome of molecular dynamics simulations, concurred to give a picture of how water-containing matrices control the internal dynamics of the embedded proteins. This occurs, in particular, via the formation of hydrogen bond networks that anchor the protein surface to the surrounding matrix, whose stiffness increases by lowering the sample water content. In the conclusion section, we also briefly speculate on how the protein-matrix interactions observed in our samples may shed light on the protein-solvent coupling also in liquid aqueous solutions.

Stability of proteins: Temperature, pressure and the role of the solvent

We focus on the various aspects of the physics related to the stability of proteins. We review the pure thermodynamic aspects of the response of a protein to pressure and temperature variations and discuss the respective stability phase diagram. We relate the experimentally observed shape of this diagram to the low degree of correlation between the fluctuations of enthalpy and volume changes associated with the folding-denaturing transition and draw attention to the fact that one order parameter is not enough to characterize the transition. We discuss in detail microscopic aspects of the various contributions to the free energy gap of proteins and put emphasis on how a cosolvent may either enlarge or diminish this gap. We review briefly the various experimental approaches to measure changes in protein stability induced by cosolvents, denaturants, but also by pressure and temperature. Finally, we discuss in detail our own molecular dynamics simulations on cytochrome c and show what happens under high pressure, how glycerol influences structure and volume fluctuations, and how all this compares with experiments.