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

The water association band as a marker of hydrogen bonds in trehalose amorphous matrices

The water association band is a suitable marker of residual water behavior in bioprotective trehalose matrices.

Resolution and Characterization of Chemical Steps in Enzyme Catalytic Sequences by Using Low-Temperature and Time-Resolved, Full-Spectrum EPR Spectroscopy in Fluid Cryosolvent and Frozen Solution Systems

Approaches to the resolution and characterization of individual chemical steps in enzyme catalytic sequences, by using temperatures in the cryogenic range of 190-250 K, and kinetics measured by time-resolved, full-spectrum electron paramagnetic resonance spectroscopy in fluid cryosolvent and frozen solution systems, are described. The preparation and performance of the adenosylcobalamin-dependent ethanolamine ammonia-lyase enzyme from Salmonella typhimurium in the two systems exemplifies the biochemical and spectroscopic methods. General advantages of low-temperature studies are (1) slowing of reaction steps, so that measurements can be made by using straightforward T-step kinetic methods and commercial instrumentation, (2) resolution of individual reaction steps, so that first-order kinetic analysis can be applied, and (3) accumulation of intermediates that are not detectable at room temperatures. The broad temperature range from room temperature to 190 K encompasses three regimes: (1) temperature-independent mean free energy surface (corresponding to native behavior); (2) the narrow temperature region of a glass-like transition in the protein, over which the free energy surface changes, revealing dependence of the native reaction on collective protein/solvent motions; and (3) the temperature range below the glass transition region, for which persistent reaction corresponds to nonnative, alternative reaction pathways, in the vicinity of the native configurational envelope. Representative outcomes of low-temperature kinetics studies are portrayed on Eyring and free energy surface (landscape) plots, and guidelines for interpretations are presented.

The importance of individual protein molecule dynamics in developing and assessing solid state protein preparations

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.

Reverse micelles as a tool for probing solvent modulation of protein dynamics: Reverse micelle encapsulated hemoglobin

Hydration waters impact protein dynamics. Dissecting the interplay between hydration waters and dynamics requires a protein that manifests a broad range of dynamics. Proteins in reverse micelles (RMs) have promise as tools to achieve this objective because the water content can be manipulated. Hemoglobin is an appropriate tool with which to probe hydration effects. We describe both a protocol for hemoglobin encapsulation in reverse micelles and a facile method using PEG and cosolvents to manipulate water content. Hydration properties are probed using the water-sensitive fluorescence from Hb bound pyranine and covalently attached Badan. Protein dynamics are probed through ligand recombination traces derived from photodissociated carbonmonoxy hemoglobin on a log scale that exposes the potential role of both α and β solvent fluctuations in modulating protein dynamics. The results open the possibility of probing hydration level phenomena in this system using a combination of NMR and optical probes.

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.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Comparative analysis of inner cavities and ligand migration in non-symbiotic AHb1 and AHb2

This study reports a comparative analysis of the topological properties of inner cavities and the intrinsic dynamics of non-symbiotic hemoglobins AHb1 and AHb2 from Arabidopsis thaliana. The two proteins belong to the 3/3 globin fold and have a sequence identity of about 60%. However, it is widely assumed that they have distinct physiological roles. In order to investigate the structure-function relationships in these proteins, we have examined the bis-histidyl and ligand-bound hexacoordinated states by atomistic simulations using in silico structural models. The results allow us to identify two main pathways to the distal cavity in the bis-histidyl hexacoordinated proteins. Nevertheless, a larger accessibility to small gaseous molecules is found in AHb2. This effect can be attributed to three factors: the mutation Leu35(AHb1)→Phe32(AHb2), the enhanced flexibility of helix B, and the more favorable energetic profile for ligand migration to the distal cavity. The net effect of these factors would be to facilitate the access of ligands, thus compensating the preference for the fully hexacoordination of AHb2, in contrast to the equilibrium between hexa- and pentacoordinated species in AHb1. On the other hand, binding of the exogenous ligand introduces distinct structural changes in the two proteins. A well-defined tunnel is formed in AHb1, which might be relevant to accomplish the proposed NO detoxification reaction. In contrast, no similar tunnel is found in AHb2, which can be ascribed to the reduced flexibility of helix E imposed by the larger number of salt bridges compared to AHb1. This feature would thus support the storage and transport functions proposed for AHb2. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Network-based models as tools hinting at nonevident protein functionality

Network-based models of proteins are popular tools employed to determine dynamic features related to the folded structure. They encompass all topological and geometric computational approaches idealizing proteins as directly interacting nodes. Topology makes use of neighborhood information of residues, and geometry includes relative placement of neighbors. Coarse-grained approaches efficiently predict alternative conformations because of inherent collectivity in the protein structure. Such collectivity is moderated by topological characteristics that also tune neighborhood structure: That rich residues have richer neighbors secures robustness toward random loss of interactions/nodes due to environmental fluctuations/mutations. Geometry conveys the additional information of force balance to network models, establishing the local shape of the energy landscape. Here, residue and/or bond perturbations are critically evaluated to suggest new experiments, as network-based computational techniques prove useful in capturing domain movements and conformational shifts resulting from environmental alterations. Evolutionarily conserved residues are optimally connected, defining a subnetwork that may be utilized for further coarsening.

Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework

Energy landscape theories have provided a common ground for understanding the protein folding problem, which once seemed to be overwhelmingly complicated. At the same time, the native state was found to be an ensemble of interconverting states with frustration playing a more important role compared to the folding problem. The landscape of the folded protein – the native landscape – is glassier than the folding landscape; hence, a general description analogous to the folding theories is difficult to achieve. On the other hand, the native basin phase volume is much smaller, allowing a protein to fully sample its native energy landscape on the biological timescales. Current computational resources may also be used to perform this sampling for smaller proteins, to build a ‘topographical map’ of the native landscape that can be used for subsequent analysis. Several major approaches to representing this topographical map are highlighted in this review, including the construction of kinetic networks, hierarchical trees and free energy surfaces with subsequent structural and kinetic analyses. In this review, we extensively discuss the important question of choosing proper collective coordinates characterizing functional motions. In many cases, the substates on the native energy landscape, which represent different functional states, can be used to obtain variables that are well suited for building free energy surfaces and analyzing the protein's functional dynamics. Normal mode analysis can provide such variables in cases where functional motions are dictated by the molecule's architecture. Principal component analysis is a more expensive way of inferring the essential variables from the protein's motions, one that requires a long molecular dynamics simulation. Finally, the two popular models for the allosteric switching mechanism, ‘preexisting equilibrium’ and ‘induced fit’, are interpreted within the energy landscape paradigm as extreme points of a continuum of transition mechanisms. Some experimental evidence illustrating each of these two models, as well as intermediate mechanisms, is presented and discussed.

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).

Glass matrix-facilitated thermal reduction: a tool for probing reactions of met hemoglobin with nitrite and nitric oxide

Isolating elemental steps that comprise a protein reaction in solution is a difficult process. In this study, the use of sugar-derived glass matrices is evaluated as a biophysical tool to help dissect out elemental steps and isolate intermediates. Two features of the glass are utilized in this endeavor: (i) the capacity of trehalose glass matrices to support thermal reduction over macroscopic distances; and (ii) the ability of glass matrices to significantly damp large amplitude protein dynamics. The focus of the study is on the reaction of nitric oxide (NO) with a nitrite ion coordinated to the heme iron of hemoglobin (Hb). The thermal reduction property of the glass is used to generate NO from nitrite within the glass, and the damping of protein dynamics is used to control entry of NO into the distal heme pocket of Hb, where it can either interact with bound nitrite or bind to the heme iron. The results not only relate to earlier controversial studies addressing the reactions of Hb with NO and nitrite but also raise the prospect that these properties of sugar-derived glassy matrices can be exploited as a new biophysical tool to modulate and probe reactions of NO with hemeproteins as well as a wide range of other metalloproteins.

Osmolyte-induced perturbations of hydrogen bonding between hydration layer waters: correlation with protein conformational changes

Gadolinium vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe osmolyte-induced changes in the hydrogen bond strength between first and second shell waters on the surface of free Gd(3+) and Gd(3+) coordinated to EDTA and to structured calcium binding peptides in solution. In parallel, Raman is used to probe the corresponding impact of the same set of osmolytes on hydrogen bonding among waters in the bulk phase. Increasing concentration of added urea is observed to progressively weaken the hydrogen bonding within the hydration layer but has minimal observed impact on bulk water. In contrast, polyols are observed to enhance hydrogen bonding in both the hydration layer and the bulk with the amplitude being polyol dependent with trehalose being more effective than sucrose, glucose, or glycerol. The observed patterns indicate that the size and properties of the osmolyte as well as the local architecture of the specific surface site of hydration impact preferential exclusion effects and local hydrogen bond strength. Correlation of the vibronic spectra with CD measurements on the peptides as a function of added osmolytes shows an increase in secondary structure with added polyols and that the progressive weakening of the hydrogen bonding upon addition of urea first increases water occupancy within the peptide and only subsequently does the peptide unfold. The results support models in which the initial steps in the unfolding process involve osmolyte-induced enhancement of water occupancy within the interior of the protein.

Nanosecond motions in proteins impose bounds on the timescale distributions of local dynamics

We elucidate the physics of protein dynamical transition via 10-100-ns molecular dynamics simulations at temperatures spanning 160-300 K. By tracking the energy fluctuations, we show that the protein dynamical transition is marked by a crossover from nonstationary to stationary processes that underlie the dynamics of protein motions. A two-timescale function captures the nonexponential character of backbone structural relaxations. One timescale is attributed to the collective segmental motions and the other to local relaxations. The former is well defined by a single-exponential, nanosecond decay, operative at all temperatures. The latter is described by a set of processes that display a distribution of timescales. Although their average remains on the picosecond timescale, the distribution is markedly contracted at the onset of the transition. It is shown that the collective motions impose bounds on timescales spanned by local dynamical processes. The nonstationary character below the transition implicates the presence of a collection of substates whose interactions are restricted. At these temperatures, a wide distribution of local-motion timescales, extending beyond that of nanoseconds, is observed. At physiological temperatures, local motions are confined to timescales faster than nanoseconds. This relatively narrow window makes possible the appearance of multiple channels for the backbone dynamics to operate.

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.

Reactivity of glass-embedded met hemoglobin derivatives toward external NO: implications for nitrite-mediated production of bioactive NO

Many protein reactions are exceedingly difficult to dissect under standard conditions due to low concentrations of reactants and intermediates. A case in point are several proposed reactions of hemoglobin with both nitrite and nitric oxide. In the present work, glassy matrices are used to dynamically control the rate at which externally introduced gaseous NO accesses and reacts with several different met Hb derivatives including the nitrite, nitrate, and aquomet forms. This novel yet general approach reveals a clear difference between nitrite and other ligands including nitrate, water, and an internal imidazole. For nitrate, water, and the internal distal imidazole, the observed spectral changes indicate that NO entering the distal heme pocket is effective in displacing these ligands from the ferric heme iron. In contrast, when the ligand is nitrite, the resulting initial spectra indicate the formation of an intermediate that has distinctly ferrous-like properties. The spectrum and the response of DAF fluorescence to the presence of the intermediate are consistent with a recently proposed nitrite anhydrase reaction. This proposed intermediate is especially significant in that it represents a pathway for a nitrite-dependent catalytic process whereby Hb generates relatively long-lived bioactive forms of NO such as S-nitrosoglutathione. The failure to form this intermediate either at low pH or when the glass is extensively dried is consistent with the requirement for a specific conformation of reactants and residue side chains within the distal heme pocket.

Charge density-dependent modifications of hydration shell waters by Hofmeister ions

Gadolinium (Gd(3+)) vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe, as a function of added Hofmeister series salts, changes in the OH stretching frequency derived from first-shell waters of aqueous Gd(3+) and of Gd(3+) coordinated to three different types of molecules: (i) a chelate (EDTA), (ii) structured peptides (mSE3/SE2) of the lanthanide-binding tags (LBTs) family with a single high-affinity binding site, and (iii) a calcium-binding protein (calmodulin) with four binding sites. The vibronic sideband (VSB) corresponding to the OH stretching mode of waters coordinated to Gd(3+), whose frequency is inversely correlated with the strength of the hydrogen bonding to neighboring waters, exhibits an increase in frequency when Gd(3+) becomes coordinated to either EDTA, calmodulin, or mSE3 peptide. In all of these cases, the addition of cation chloride or acetate salts to the solution increases the frequency of the vibronic band originating from the OH stretching mode of the coordinated waters in a cation- and concentration-dependent fashion. The cation dependence of the frequency increase scales with charge density of the cations, giving rise to an ordering consistent with the Hofmeister ordering. On the other hand, water Raman spectroscopy shows no significant change upon addition of these salts. Additionally, it is shown that the cation effect is modulated by the specific anion used. The results indicate a mechanism of action for Hofmeister series ions in which hydrogen bonding among hydration shell waters is modulated by several factors. High charge density cations sequester waters in a configuration that precludes strong hydrogen bonding to neighboring waters. Under such conditions, anion effects emerge as anions compete for hydrogen-bonding sites with the remaining free waters on the surface of the hydration shell. The magnitude of the anion effect is both cation and Gd(3+)-binding site specific.

Ultrafast catalytic processes and conformational changes in the light-driven enzyme protochlorophyllide oxidoreductase (POR)

The enzyme POR (protochlorophyllide oxidoreductase), from the family of alcohol dehydrogenases, reduces protochlorophyllide into chlorophyllide on the absorption of light. The reduction involves the transfer of two protons and two electrons and is an important regulatory step in the biosynthesis of chlorophyll. In recent years, due to the availability of large quantities of the pure enzyme, much of the catalytic reaction has been unravelled by using a variety of spectroscopic methods, including ultrafast initial events in catalysis. In addition, it has been demonstrated that a light-activated conformational change of the protein is necessary to activate catalysis. This makes POR a very important model system to study the relationship between structural changes of enzymes and functionality.

Ligand binding to truncated hemoglobin N from Mycobacterium tuberculosis is strongly modulated by the interplay between the distal heme pocket residues and internal water

The survival of Mycobacterium tuberculosis requires detoxification of host *NO. Oxygenated Mycobacterium tuberculosis truncated hemoglobin N catalyzes the rapid oxidation of nitric oxide to innocuous nitrate with a second-order rate constant (k'(NOD) approximately 745 x 10(6) m(-1) x s(-1)), which is approximately 15-fold faster than the reaction of horse heart myoglobin. We ask what aspects of structure and/or dynamics give rise to this enhanced reactivity. A first step is to expose what controls ligand/substrate binding to the heme. We present evidence that the main barrier to ligand binding to deoxy-truncated hemoglobin N (deoxy-trHbN) is the displacement of a distal cavity water molecule, which is mainly stabilized by residue Tyr(B10) but not coordinated to the heme iron. As observed in the Tyr(B10)/Gln(E11) apolar mutants, once this kinetic barrier is lowered, CO and O(2) binding is very rapid with rates approaching 1-2 x 10(9) m(-1) x s(-1). These large values almost certainly represent the upper limit for ligand binding to a heme protein and also indicate that the iron atom in trHbN is highly reactive. Kinetic measurements on the photoproduct of the *NO derivative of met-trHbN, where both the *NO and water can be directly followed, revealed that water rebinding is quite fast (approximately 1.49 x 10(8) s(-1)) and is responsible for the low geminate yield in trHbN. Molecular dynamics simulations, performed with trHbN and its distal mutants, indicated that in the absence of a distal water molecule, ligand access to the heme iron is not hindered. They also showed that a water molecule is stabilized next to the heme iron through hydrogen-bonding with Tyr(B10) and Gln(E11).

Influence of distal residue B10 on CO dynamics in myoglobin and neuroglobin

For many years, myoglobin has served as a paradigm for structure-function studies in proteins. Ligand binding and migration within myoglobin has been studied in great detail by crystallography and spectroscopy, showing that gaseous ligands such as O(2), CO, and NO not only bind to the heme iron but may also reside transiently in three internal ligand docking sites, the primary docking site B and secondary sites C and D. These sites affect ligand association and dissociation in specific ways. Neuroglobin is another vertebrate heme protein that also binds small ligands. Ligand migration pathways in neuroglobin have not yet been elucidated. Here, we have used Fourier transform infrared temperature derivative spectroscopy at cryogenic temperatures to compare the influence of the side chain volume of amino acid residue B10 on ligand migration to and rebinding from docking sites in myoglobin and neuroglobin.

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of the dimeric hemoglobin of Scapharca inaequivalvis (HbI) after photodissociation. By combining these studies with X-ray crystallography, three transient ligand docking sites were identified: a primary docking site B in close vicinity to the heme iron, and two secondary docking sites C and D corresponding to the Xe4 and Xe2 cavities of myoglobin. To assess the relevance of these findings for physiological binding, we also performed flash photolysis experiments on HbICO at room temperature and equilibrium binding studies with dioxygen. Our results show that the Xe4 and Xe2 cavities serve as transient docking sites for unbound ligands in the protein, but not as way stations on the entry/exit pathway. For HbI, the so-called histidine gate mechanism proposed for other globins appears as a plausible entry/exit route as well.

Ligand migration in nonsymbiotic hemoglobin AHb1 from Arabidopsis thaliana

AHb1 is a hexacoordinated type 1 nonsymbiotic hemoglobin recently discovered in Arabidopsis thaliana. To gain insight into the ligand migration inside the protein, we studied the CO rebinding kinetics of AHb1 encapsulated in silica gels, in the presence of glycerol. The CO rebinding kinetics after nanosecond laser flash photolysis exhibits complex ligand migration patterns, consistent with the existence of discrete docking sites in which ligands can temporarily be stored before rebinding to the heme at different times. This finding may be of relevance to the physiological NO dioxygenase activity of this protein, which requires sequential binding of two substrates, NO and O2, to the heme.

Conformational dependence of hemoglobin reactivity under high viscosity conditions: the role of solvent slaved dynamics

The concept of protein dynamic states is introduced. This concept is based on (i) protein dynamics being organized hierarchically with respect to solvent slaving and (ii) which tier of dynamics is operative over the time window of a given measurement. The protein dynamic state concept is used to analyze the kinetic phases derived from the recombination of carbon monoxide to sol-gel-encapsulated human adult hemoglobin (HbA) and select recombinant mutants. The temperature-dependent measurements are made under very high viscosity conditions obtained by bathing the samples in an excess of glycerol. The results are consistent with a given tier of solvent slaved dynamics becoming operative at a time delay (with respect to the onset of the measurement) that is primarily solvent- and temperature-dependent. However, the functional consequences of the dynamics are protein- and conformation-specific. The kinetic traces from both equilibrium populations and trapped allosteric intermediates show a consistent progression that exposes the role of both conformation and hydration in the control of reactivity. Iron-zinc symmetric hybrid forms of HbA are used to show the dramatic difference between the kinetic patterns for T state alpha and beta subunits. The overall results support a model for allostery in HbA in which the ligand-binding-induced transition from the deoxy T state to the high -affinity R state proceeds through a progression of T state intermediates.

How a vicinal layer of solvent modulates the dynamics of proteins

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.