Dynamics of hemoglobin in human erythrocytes and in solution: influence of viscosity studied by ultrafast vibrational echo experiments

Blood-Blister Aneurysms of the Internal Carotid Artery in Tibetan and Han Populations : A Retrospective Observational Study

OBJECTIVE

Blood-blister aneurysms (BBAs) of the internal carotid artery (ICA) are challenging lesions with high morbidity and mortality rates. Although research on BBAs is well documented in different populations, the study of BBAs in the Tibetan population is extremely rare. This study aimed to evaluate the characteristics of BBAs and analyze the treatment modalities and long-term outcomes in the Tibetan population in comparison with the Han population.

METHODS

The characteristics of patients with BBAs of the ICA from January 2009 to January 2021 at our institution were reviewed. The features of aneurysms, treatment modalities, complications, and follow-up outcomes were retrospectively analyzed.

RESULTS

A total of 130 patients (41 Tibetan and 89 Han patients) with BBAs of the ICA who underwent treatment were enrolled. Compared with the Han group, the Tibetan group significantly demonstrated a high ratio of BBAs among ICAs (8.6%, 41/477 vs. 1.6%, 89/5563; p<0.05), a high ratio of vasospasm (34.1%, 14/41 vs. 6.7%, 6/89; p=0.001), a high risk of ischemic events (43.9%, 18/41 vs. 22.5%, 20/89; p<0.05), and a low ratio of good outcomes (modified Rankin scale, 0-2) at the 1-year follow-up (51.2%, 21/41 vs. 74.2%, 66/89; p<0.05). The multivariate regression model showed that ischemic events significantly contributed to the prediction of outcomes at 1 year. Further analysis revealed that microsurgery and vasospasm were associated with ischemic events.

CONCLUSION

In comparison with Han patients, the Tibetan population had a high ratio of BBA occurrence, a high incidence of ischemic events, and a high ratio of poor outcomes. The endovascular approach showed more benefits in BBA patients.

Transparent window 2D IR spectroscopy of proteins

Proteins are complex, heterogeneous macromolecules that exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires experimental tools to characterize them with high spatial and temporal precision. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution. Two-dimensional (2D) IR methods that provide richer information are becoming more routine but remain challenging to apply to proteins. Spectral congestion typically prevents selective investigation of native vibrations; however, the problem can be overcome by site-specific introduction of amino acid side chains that have vibrational groups with frequencies in the "transparent window" of protein spectra. This Perspective provides an overview of the history and recent progress in the development of transparent window 2D IR of proteins.

The unique structural features of carbonmonoxy hemoglobin from the sub-Antarctic fish Eleginops maclovinus

Tetrameric hemoglobins (Hbs) are prototypical systems for the investigations of fundamental properties of proteins. Although the structure of these proteins has been known for nearly sixty years, there are many aspects related to their function/structure that are still obscure. Here, we report the crystal structure of a carbonmonoxy form of the Hb isolated from the sub-Antarctic notothenioid fish Eleginops maclovinus characterised by either rare or unique features. In particular, the distal site of the α chain results to be very unusual since the distal His is displaced from its canonical position. This displacement is coupled with a shortening of the highly conserved E helix and the formation of novel interactions at tertiary structure level. Interestingly, the quaternary structure is closer to the T-deoxy state of Hbs than to the R-state despite the full coordination of all chains. Notably, these peculiar structural features provide a rationale for some spectroscopic properties exhibited by the protein in solution. Finally, this unexpected structural plasticity of the heme distal side has been associated with specific sequence signatures of various Hbs.

Elevated hemoglobin is associated with poor prognosis in Tibetans with poor-grade aneurysmal subarachnoid hemorrhage after clipping: A Retrospective Case-Control Study

OBJECTIVES

High hemoglobin (HGB) concentration is frequently seen in Tibetans in clinical practice; however, the impact on postsurgical prognosis in patients with poor-grade aneurysmal subarachnoid hemorrhage (aSAH) is not precisely known. Thus, we sought to understand the association between high HGB level and postoperative outcome in Tibetans with poor-grade aSAH.

PATIENTS AND METHODS

Results of clipping in consecutive Tibetan patients with poor-grade aSAH were analyzed retrospectively for the period January 2012 to January 2017. Based on the upper limit (160g/L) of normal hemoglobin levels, patients were divided to a high (HHC) and a normal (NHC) HGB-level cohort according to the first routine blood result on admission. Propensity score matching was used for baseline matching in the 2 cohorts. Postoperative complications in the 2 groups were compared. Prognosis after ictus, including 6-month neurological functional status and mortality at 30 days and 6 months were also assessed.

RESULTS

Risk of ischemia, pulmonary embolism and lower-limb deep venous thrombosis (DVT) was higher in HHC than NHC (62.88% vs. 21.64%; P<0.001; 10.30% vs. 1.31%, P<0.005; 24.74% vs. 7.21%, P<0.001, respectively). Hospital stay also differed significantly (15.82±3.86 vs. 10.37±4.80 days; P<0.001). Out of the 194 patients, 150 survived at 6 months. At 6-month neurological functional follow-up, 8 NHC patients had favorable modified Rankin scale (mRS) scores ≤2 at discharge, versus only 1 HHC patient, showing better outcome in NHC than HHC (8.25% vs. 1.03%; P=0.035). In-hospital mortality was significantly greater in HHC than NHC (17.52% vs. 7.22%; P=0.029). 30-day post-ictus mortality was 30.93% in HHC versus 14.43% in NHC (P<0.006). There was also a significant difference in mortality at 6 months post-ictus (47.42% vs. 18.56%; P<0.001).

CONCLUSION

High HGB level was associated with increased risk of postsurgical cerebral ischemia, pulmonary embolism and lower-limb DVT and poor prognosis in poor-grade aSAH patients. Preoperative hemodilution therapy might be beneficial in reducing operative complications, reducing hospital stay and improving short-term prognosis for neurological functional recovery in aSAH patients with high HGB concentration, but further detailed research is needed.

Multispectral multidimensional spectrometer spanning the ultraviolet to the mid-infrared

Multidimensional spectroscopy is the optical analog to nuclear magnetic resonance, probing dynamical processes with ultrafast time resolution. At optical frequencies, the technical challenges of multidimensional spectroscopy have hindered its progress until recently, where advances in laser sources and pulse-shaping have removed many obstacles to its implementation. Multidimensional spectroscopy in the visible and infrared (IR) regimes has already enabled respective advances in our understanding of photosynthesis and the structural rearrangements of liquid water. A frontier of ultrafast spectroscopy is to extend and combine multidimensional techniques and frequency ranges, which have been largely restricted to operating in the distinct visible or IR regimes. By employing two independent amplifiers seeded by a single oscillator, it is straightforward to span a wide range of time scales (femtoseconds to seconds), all of which are often relevant to the most important energy conversion and catalysis problems in chemistry, physics, and materials science. Complex condensed phase systems have optical transitions spanning the ultraviolet (UV) to the IR and exhibit dynamics relevant to function on time scales of femtoseconds to seconds and beyond. We describe the development of the Multispectral Multidimensional Nonlinear Spectrometer (MMDS) to enable studies of dynamical processes in atomic, molecular, and material systems spanning femtoseconds to seconds, from the UV to the IR regimes. The MMDS employs pulse-shaping methods to provide an easy-to-use instrument with an unprecedented spectral range that enables unique combination spectroscopies. We demonstrate the multispectral capabilities of the MMDS on several model systems.

Viscosity Control of Protein Solution by Small Solutes: A Review

Viscosity of protein solution is one of the most troublesome issues for the high-concentration formulation of protein drugs. In this review, we summarize the practical methods that suppress the viscosity of protein solution using small molecular additives. The small amount of salts decreases the viscosity that results from electrostatic repulsion and attraction. The chaotrope suppresses the hydrophobic attraction and cluster formation, which can lower the solution viscosity. Arginine hydrochloride (ArgHCl) also suppresses the solution viscosity due to the hydrophobic and aromatic interactions between protein molecules. The small molecular additives are the simplest resolution of the high viscosity of protein solution as well as understanding of the primary cause in complex phenomena of protein interactions.

Photophysical behavior of a potential drug candidate, trans-[2-(4-methoxystyryl)]quinoline-1-oxide tuned by environment effects

Styryl and/or quinoline structural fragments, present in a large number of bioactive substances, inspired the design of various new drug candidates. In this paper, we describe the photophysical behavior of trans-[2-(4-methoxystyryl)]quinoline-1-oxide (trans-MSQNO) on the basis of X-ray analysis data, theoretical calculations as well as steady state and time-resolved spectroscopy experiments in various media. The molecule crystallizes in orthorhombic unit cell containing eight molecules of N-oxide, space group Pbca. The NO bond is substantially shorter in comparison with the NO bond in the ZnTPP unit [1.3052(11) Å vs. 1.335(2) Å]. Variation of emission colors from the violet (~450nm) through blue (480nm), green (525nm) and yellow (575nm) is observed in different environments. Comparable values of lifetimes estimated both at ambient temperature and at 77K suggest that excited state dynamics in this case is viscosity independent. DFT and TD DFT B3LYP/6-31G(d, p) calculations performed for four different trans-MSQNO rotamers in the gas phase, as well as nonpolar and polar media (PCM model) suggest that an equilibrium between them can be significantly altered even by a relatively weak interactions with the environment. It is suggested that varying intensity ratios of experimental absorption bands in different media may be due to the dominant share of one or more rotamers of the excited trans-MSQNO molecule. Gas phase calculations show also that the vertical ππ*, S₀→S₁, transition resulting from the HOMO→LUMO electronic configuration exhibits only a partial CT nature. On the other hand, in polar media, a substantial increment of excited state dipole moment of all rotamers compared to the ground state, its increase with increasing solvent polarity and a significant red shift in the absorption and emission spectra, point to the enhanced CT nature of the S₁ excited state. Hence, the trans-MSQNO molecule may be considered a subsequent styrylquinoline drug candidate where the CT drug-receptor interactions are of a high importance.

Histidine Orientation Modulates the Structure and Dynamics of a de Novo Metalloenzyme Active Site

The ultrafast dynamics of a de novo metalloenzyme active site is monitored using two-dimensional infrared spectroscopy. The homotrimer of parallel, coiled coil α-helices contains a His3-Cu(I) metal site where CO is bound and serves as a vibrational probe of the hydrophobic interior of the self-assembled complex. The ultrafast spectral dynamics of Cu-CO reveals unprecedented ultrafast (2 ps) nonequilibrium structural rearrangements launched by vibrational excitation of CO. This initial rapid phase is followed by much slower ∼40 ps vibrational relaxation typical of metal-CO vibrations in natural proteins. To identify the hidden coupled coordinate, small molecule analogues and the full peptide were studied by QM and QM/MM calculations, respectively. The calculations show that variation of the histidines' dihedral angles in coordinating Cu controls the coupling between the CO stretch and the Cu-C-O bending coordinates. Analysis of different optimized structures with significantly different electrostatic field magnitudes at the CO ligand site indicates that the origin of the stretch-bend coupling is not directly due to through-space electrostatics. Instead, the large, ∼3.6 D dipole moments of the histidine side chains effectively transduce the electrostatic environment to the local metal coordination orientation. The sensitivity of the first coordination sphere to the protein electrostatics and its role in altering the potential energy surface of the bound ligands suggests that long-range electrostatics can be leveraged to fine-tune function through enzyme design.

Ultrafast Dynamics of Carboxy-Hemoglobin: Two-Dimensional Infrared Spectroscopy Experiments and Simulations

This Letter presents a comparison between experimental and simulated 2D mid-infrared spectra of carboxy-hemoglobin in the spectral region of the carbon monoxide stretching mode. The simulations rely on a fluctuating potential energy surface that includes both the effect of heme and the protein surroundings computed from molecular dynamics simulations. A very good agreement between theory and experiment is obtained with no adjustable parameters. The simulations show that the effect of the distal histidine through the hydrogen bond is strong and is directly responsible for the slow decay of the frequency-frequency correlation function on a 10 ps time scale. This study confirms that fluctuations in carboxy-hemoglobin are more inhomogeneous than those in the more frequently studied carboxy-myoglobin. The comparison between simulations and experiments brings valuable information on the complex relation between protein structure and spectral diffusion.

Crowding induced collective hydration of biological macromolecules over extended distances

Ultrafast two-dimensional infrared (2D-IR) spectroscopy reveals picosecond protein and hydration dynamics of crowded hen egg white lysozyme (HEWL) labeled with a metal-carbonyl vibrational probe covalently attached to a solvent accessible His residue. HEWL is systematically crowded alternatively with polyethylene glycol (PEG) or excess lysozyme in order to distinguish the chemically inert polymer from the complex electrostatic profile of the protein crowder. The results are threefold: (1) A sharp dynamical jamming-like transition is observed in the picosecond protein and hydration dynamics that is attributed to an independent-to-collective hydration transition induced by macromolecular crowding that slows the hydration dynamics up to an order of magnitude relative to bulk water. (2) The interprotein distance at which the transition occurs suggests collective hydration of proteins over distances of 30-40 Å. (3) Comparing the crowding effects of PEG400 to our previously reported experiments using glycerol exposes fundamental differences between small and macromolecular crowding agents.

Two-dimensional infrared spectroscopy of azido-nicotinamide adenine dinucleotide in water

Mid-IR active analogs of enzyme cofactors have the potential to be important spectroscopic reporters of enzyme active site dynamics. Azido-nicotinamide adenine dinucleotide (NAD(+)), which has been recently synthesized in our laboratory, is a mid-IR active analog of NAD(+), a ubiquitous redox cofactor in biology. In this study, we measure the frequency-frequency time correlation function for the antisymmetric stretching vibration of the azido group of azido-NAD(+) in water. Our results are consistent with previous studies of pseudohalides in water. We conclude that azido-NAD(+) is sensitive to local environmental fluctuations, which, in water, are dominated by hydrogen-bond dynamics of the water molecules around the probe. Our results demonstrate the potential of azido-NAD(+) as a vibrational probe and illustrate the potential of substituted NAD(+)-analogs as reporters of local structural dynamics that could be used for studies of protein dynamics in NAD-dependent enzymes.

Characterizing the dynamics of functionally relevant complexes of formate dehydrogenase

The potential for femtosecond to picosecond time-scale motions to influence the rate of the intrinsic chemical step in enzyme-catalyzed reactions is a source of significant controversy. Among the central challenges in resolving this controversy is the difficulty of experimentally characterizing thermally activated motions at this time scale in functionally relevant enzyme complexes. We report a series of measurements to address this problem using two-dimensional infrared spectroscopy to characterize the time scales of active-site motions in complexes of formate dehydrogenase with the transition-state-analog inhibitor azide (N(3)(-)). We observe that the frequency-frequency time correlation functions (FFCF) for the ternary complexes with NAD(+) and NADH decay completely with slow time constants of 3.2 ps and 4.6 ps, respectively. This result suggests that in the vicinity of the transition state, the active-site enzyme structure samples a narrow and relatively rigid conformational distribution indicating that the transition-state structure is well organized for the reaction. In contrast, for the binary complex, we observe a significant static contribution to the FFCF similar to what is seen in other enzymes, indicating the presence of the slow motions that occur on time scales longer than our measurement window.

Rheological and syringeability properties of highly concentrated human polyclonal immunoglobulin solutions

This study of highly concentrated polyvalent immunoglobulin solutions, IgG, aimed at analyzing the relationships between protein concentration and aggregation on the one hand and viscosity on the other hand. Viscosity variations as a function of IgG concentration showed two well-defined behaviours: a Newtonian behaviour for low-concentrated solutions and a shear-thinning behaviour for highly concentrated ones. The viscosity data fitted very well with the Mooney model, suggesting the absence of intermolecular interactions in the IgG solutions that behaved like a non-interacting suspension of hard particles. The polyclonal nature of IgG seems to prevent intermolecular interaction. The shape factor, determined from Mooney fitting, revealed a non-spherical shape of the polyclonal IgG molecules. The rheological properties were also correlated with the injection force (F) through hypodermic needles by syringeability tests. Here, F was mainly affected by three parameters: the solution viscosity, the injection flow rate, and the needle characteristics. In fact, syringeability tests showed that F increased with IgG concentration and flow rate and decreased with the internal diameter of the needle. A zone for optimal injection conditions was then identified taking into account the different affecting parameters and mainly a maximum force for manual injection, which was fixed at 30N.

Early events associated with the excited state proton transfer in 2-(2'-pyridyl)benzimidazole

2-(2(')-Pyridyl)benzimidazole (2PBI) undergoes excited state proton transfer (ESPT) in acidic solutions, leading to a tautomer emission at 460 nm. This photoprocess has been studied using ultrafast fluorescence spectroscopic techniques in acidic neat aqueous solutions, in viscous mixtures of glycerol with water, as well as in sucrose solutions. The tautomer is found to be stabilized in the more viscous medium, leading to a greater relative quantum yield as well as lifetime. The long rise time in tautomer emission is not affected by viscosity though. Rather, it appears to have the same value as the long component of the decay of the cationic excited state (C(*)). In addition to the subnanosecond lifetime reported earlier, C(*) is found to exhibit a decay time of 2 ps. This is assigned to its protonation to form the nonfluorescent dication in its excited state (D(*)) considering the ground and excited state pK(a) values reported earlier. An additional rising component of 100 ps is observed in the region of C(*) emission. This is likely to arise from a structural change or charge redistribution in C(*) immediately after its creation and before the phototautomerization.

Exploring the molecular origins of protein dynamics in the active site of human carbonic anhydrase II

We present three-pulse vibrational echo measurements of azide ion bound to the active site Zn of human carbonic anhydrase II (HCA II) and of two separate active-site mutants Thr199 --> Ala (T199A) and Leu198 --> Phe (L198F). Because structural motions of the protein active site influence the frequency of bound ligands, the differences in the time scales of the frequency-frequency correlation functions (FFCFs) obtained from global fits to each set of data allow us to make inferences about the time scales of the active site dynamics of HCA II. Surprisingly, the deletion of a potential electrostatic interaction in results in very little change in the FFCF, but the insertion of the bulky phenylalanine ring in causes much faster dynamics. We conclude that the fast, sub-picosecond time scale in the correlation function is attributable to hydrogen bond dynamics, and the slow, apparently static contribution is due to the conformational flexibility of Zn-bound azide in the active site.

Viscosity-dependent protein dynamics

Spectrally resolved stimulated vibrational echo spectroscopy is used to investigate the dependence of fast protein dynamics on bulk solution viscosity at room temperature in four heme proteins: hemoglobin, myoglobin, a myoglobin mutant with the distal histidine replaced by a valine (H64V), and a cytochrome c552 mutant with the distal methionine replaced by an alanine (M61A). Fructose is added to increase the viscosity of the aqueous protein solutions over many orders of magnitude. The fast dynamics of the four globular proteins were found to be sensitive to solution viscosity and asymptotically approached the dynamical behavior that was previously observed in room temperature sugar glasses. The viscosity-dependent protein dynamics are analyzed in the context of a viscoelastic relaxation model that treats the protein as a deformable breathing sphere. The viscoelastic model is in qualitative agreement with the experimental data but does not capture sufficient system detail to offer a quantitative description of the underlying fluctuation amplitudes and relaxation rates. A calibration method based on the near-infrared spectrum of water overtones was constructed to accurately determine the viscosity of small volumes of protein solutions.

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Interpreting DNA Vibrational Circular Dichroism Spectra Using a Coupling Model from Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy was recently used to measure the vibrational couplings between carbonyl bonds located on DNA nucleobases (Krummel, A. T.; Mukherjee, P.; Zanni, M. T. J. Phys. Chem. B 2003, 107, 9165 and Krummel, A. T.; Zanni, M. T. J. Phys. Chem. B 2006, 110, 13991). Here, we extend the coupling model derived from these 2D IR experiments to simulate the vibrational absorption and vibrational circular dichroism (VCD) spectra of three double-stranded DNA oligomers: poly(dG)-poly(dC), poly(dG-dC), and dGGCC. Using this model, we determine that the VCD spectrum of A-form poly(dG)-poly(dC) is dominated by interactions between stacked bases, whereas the coupling between base pairs and stacked bases carries equal importance in the VCD spectrum of B-form poly(dG-dC). We also simulate the absorption and VCD spectra of dGGCC, which is a combination of A- and B-form configurations. These simulations give insight into the structural interpretation of VCD and absorption spectroscopies that have long been used to monitor DNA secondary structure and kinetics.

Dynamics of proteins encapsulated in silica sol-gel glasses studied with IR vibrational echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo spectroscopy is used to measure the fast dynamics of heme-bound CO in carbonmonoxy-myoglobin (MbCO) and -hemoglobin (HbCO) embedded in silica sol-gel glasses. On the time scale of approximately 100 fs to several picoseconds, the vibrational dephasing of the heme-bound CO is measurably slower for both MbCO and HbCO relative to that of aqueous protein solutions. The fast structural dynamics of MbCO, as sensed by the heme-bound CO, are influenced more by the sol-gel environment than those of HbCO. Longer time scale structural dynamics (tens of picoseconds), as measured by the extent of spectral diffusion, are the same for both proteins encapsulated in sol-gel glasses compared to that in aqueous solutions. A comparison of the sol-gel experimental results to viscosity-dependent vibrational echo data taken on various mixtures of water and fructose shows that the sol-gel-encapsulated MbCO exhibits dynamics that are the equivalent of the protein in a solution that is nearly 20 times more viscous than bulk water. In contrast, the HbCO dephasing in the sol-gel reflects only a 2-fold increase in viscosity. Attempts to alter the encapsulating pore size by varying the molar ratio of silane precursor to water (R value) used to prepare the sol-gel glasses were found to have no effect on the fast or steady-state spectroscopic results. The vibrational echo data are discussed in the context of solvent confinement and protein-pore wall interactions to provide insights into the influence of a confined environment on the fast structural dynamics experienced by a biomolecule.

Effects of Solvent Viscosity on Ligand Interconversion Dynamics in the Primary Docking Site of Heme Proteins

Interconversion dynamics of the ligand in the primary docking site of myoglobin (Mb) and hemoglobin (Hb) in trehalose and glycerol/D2O mixtures at 283 K was investigated by probing time-resolved vibrational spectra of CO photolyzed from these proteins. The interconversion dynamics in viscous media are similar to those in aqueous solution, indicating that it is minimally coupled to the solvent-coupled large-scale protein motion. Interconversion rates in the heme pocket of Hb in water solution are slower than those of Mb in trehalose glass, suggesting that the interconversion barrier in Hb is intrinsically higher than that in Mb and is not modified by the solvent viscosity.

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.