The protein-solvent glass transition

Temperature-dependent iron motion in extremophile rubredoxins - no need for 'corresponding states'

Extremophile organisms are known that can metabolize at temperatures down to - 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins - rubredoxins - from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of 'corresponding states' posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although 'corresponding states' would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

Do Ionic Liquids Slow Down in Stages?

High impact recent articles have reported on the existence of a liquid-liquid (L-L) phase transition as a function of both pressure and temperature in ionic liquids (ILs) containing the popular trihexyltetradecylphosphonium cation (P666,14+), sometimes referred to as the "universal liquifier". The work presented here reports on the structural-dynamic pathway from liquid to glass of the most well-studied IL comprising the P666,14+ cation. We present experimental and computational evidence that, on cooling, the path from the room-temperature liquid to the glass state is one of separate structural-dynamic changes. The first stage involves the slowdown of the charge network, while the apolar subcomponent is fully mobile. A second, separate stage entails the slowdown of the apolar domain. Whereas it is possible that these processes may be related to the liquid-liquid and glass transitions, more research is needed to establish this conclusively.

Molecular Dynamics Simulation of the Influence of Temperature and Salt on the Dynamic Hydration Layer in a Model Polyzwitterionic Polymer PAEDAPS

We investigate the hydration of poly(3-[2-(acrylamido) ethyldimethylammonio] propanesulfonate) over a range of temperatures in pure water and with the inclusion of 0.1 mol/L NaCl using atomistic molecular dynamics simulation. Drawing on concepts drawn from the field of glass-forming liquids, we use the Debye-Waller parameter () for describing the water mobility gradient around the polybetaine backbone extending to an overall distance ≈18 Å. The water mobility in this layer is defined through the mean-square water molecule displacement at a time on the order of water's β-relaxation time. The brushlike topology of polybetaines leads to two regions in the dynamic hydration layer. The inner region of ≈10.5 Å is explored by pendant group conformational motions, and the outer region of ≈7.5 Å represents an extended layer of reduced water mobility relative to bulk water. The dynamic hydration layer extends far beyond the static hydration layer, adjacent to the polymer.

Glass-like behavior of intercalated organic solvents in graphite oxide detected by spin-probe EPR

Membranes based on graphite oxide (GO) are promising materials for the separation of polar liquids and gases. Understanding the properties of solvents immersed in GO is important for the development of various technological applications. Here, the molecular motions of the TEMPO nitroxide spin probe in acetonitrile intercalated into the GO inter-plane space were studied using electron paramagnetic resonance (EPR), including its pulsed version, and electron spin echo (ESE). For a sample containing 75% acetonitrile relative to equilibrium sorption at room temperature, ESE-detected stochastic librations were observed for TEMPO molecules above 135 K. Since these librations are an inherent property of molecular glasses, this fact indicates that intercalated acetonitrile forms a two-dimensional glass state. Above 225 K, an acceleration of stochastic librations was observed, indicating the manifestation of a glass-like dynamical cross-over. Continuous wave (CW) EPR spectra of TEMPO showed the absence of overall tumbling motions in the entire investigated temperature range of up to 340 K, indicating that the intercalated acetonitrile does not behave as a bulk liquid (the melting point of acetonitrile is 229 K). Dynamical librations of TEMPO molecules detected by CW EPR were found to accelerate above 240 K.

The Role of Polymer Chain Stiffness and Guest Nanoparticle Loading in Improving the Glass Transition Temperature of Polymer Nanocomposites

The impact of polymer chain stiffness characterized by the bending modulus (k_θ) on the glass transition temperature (T_g) of pure polymer systems, as well as polymer nanocomposites (PNCs), is investigated using molecular dynamics simulations. At small k_θ values, the pure polymer system and respective PNCs are in an amorphous state, whereas at large k_θ values, both systems are in a semicrystalline state with a glass transition at low temperature. For the pure polymer system, T_g initially increases with k_θ and does not change obviously at large k_θ. However, the T_g of PNCs shows interesting behaviors with the increasing volume fraction of nanoparticles (f_NP) at different k_θ values. T_g tends to increase with f_NP at small k_θ, whereas it becomes suppressed at large k_θ.

Reorganization energy of electron transfer

The theory of electron transfer reactions establishes the conceptual foundation for redox solution chemistry, electrochemistry, and bioenergetics. Electron and proton transfer across the cellular membrane provide all energy of life gained through natural photosynthesis and mitochondrial respiration. Rates of biological charge transfer set kinetic bottlenecks for biological energy storage. The main system-specific parameter determining the activation barrier for a single electron-transfer hop is the reorganization energy of the medium. Both harvesting of light energy in natural and artificial photosynthesis and efficient electron transport in biological energy chains require reduction of the reorganization energy to allow fast transitions. This review article discusses mechanisms by which small values of the reorganization energy are achieved in protein electron transfer and how similar mechanisms can operate in other media, such as nonpolar and ionic liquids. One of the major mechanisms of reorganization energy reduction is through non-Gibbsian (nonergodic) sampling of the medium configurations on the reaction time. A number of alternative mechanisms, such as electrowetting of active sites of proteins, give rise to non-parabolic free energy surfaces of electron transfer. These mechanisms, and nonequilibrium population of donor-acceptor vibrations, lead to a universal phenomenology of separation between the Stokes shift and variance reorganization energies of electron transfer.

Effect of Water Deuteration on Protein Electron Transfer

Traditional theories of long-range protein electron transfer describe the reaction rate in terms of the tunneling distance and the reaction free energy. They do not recognize two physical effects: (i) local wetting of the active site by hydration water and (ii) protein identity affecting the rate through dynamics and flexibility. We find, by molecular dynamics simulations, a significant, ∼25 times, slowing down of the rate of protein electron transfer upon deuteration. H/D substitution changes the rate constant pre-exponential factor in the regime of electron transfer controlled by medium dynamics. Switching from light to heavy water increases the effective medium relaxation time. The effect is caused by both a global change in the flexibility of the protein backbone and locally stronger hydrogen bonds to charged residues.

Determining biomolecular structures near room temperature using X-ray crystallography: concepts, methods and future optimization

For roughly two decades, cryocrystallography has been the overwhelmingly dominant method for determining high-resolution biomolecular structures. Competition from single-particle cryo-electron microscopy and micro-electron diffraction, increased interest in functionally relevant information that may be missing or corrupted in structures determined at cryogenic temperature, and interest in time-resolved studies of the biomolecular response to chemical and optical stimuli have driven renewed interest in data collection at room temperature and, more generally, at temperatures from the protein-solvent glass transition near 200 K to ∼350 K. Fischer has recently reviewed practical methods for room-temperature data collection and analysis [Fischer (2021), Q. Rev. Biophys. 54, e1]. Here, the key advantages and physical principles of, and methods for, crystallographic data collection at noncryogenic temperatures and some factors relevant to interpreting the resulting data are discussed. For room-temperature data collection to realize its potential within the structural biology toolkit, streamlined and standardized methods for delivering crystals prepared in the home laboratory to the synchrotron and for automated handling and data collection, similar to those for cryocrystallography, should be implemented.

Vibrational Analysis of Hydration-Layer Water around Ubiquitin, Unpeeled Layer by Layer: Molecular-Dynamics Perceptions

Classical molecular-dynamics simulations have been performed to examine the interplay between ubiquitin and its hydration-water sub-layers, chiefly from a vibrational-mode and IR viewpoint-where we analyse individual sub-layers characteristics. The vibrational Density of States (VDOS) revealed that the first solvation sub-shell indicates a confined character therein. For layers of increasing distance from the surface, the adoption of greater bulk-like spectral behaviour was evident, suggesting that vibrational harmonisation to bulk occurs within 6-7 Å of the surface.

Exploring the Limits of Biological Complexity Amenable to Studies by Incoherent Neutron Spectroscopy

The wavelengths of neutrons available at neutron scattering facilities are comparable with intra- and inter-molecular distances, while their energies are comparable with molecular vibrational energies, making such neutrons highly suitable for studies of molecular-level dynamics. The unmistakable trend in neutron spectroscopy has been towards measurements of systems of greater complexity. Several decades of studies of dynamics using neutron scattering have witnessed a progression from measurements of solids to liquids to protein complexes and biomembranes, which may exhibit properties characteristic of both solids and liquids. Over the last two decades, the frontier of complexity amenable to neutron spectroscopy studies has reached the level of cells. Considering this a baseline for neutron spectroscopy of systems of the utmost biological complexity, we briefly review what has been learned to date from neutron scattering studies at the cellular level and then discuss in more detail the recent strides into neutron spectroscopy of tissues and whole multicellular organisms.

Resolution and characterization of confinement- and temperature-dependent dynamics in solvent phases that surround proteins in frozen aqueous solution by using spin-probe EPR spectroscopy

Spin probe electron paramagnetic resonance spectroscopy is applied to characterize the dynamics of concentric hydration and mesophase solvent domains that surround proteins within the ice boundary in frozen aqueous solutions. The solvent dynamics are tuned by variation of temperature (190-265K) and by the degree of ice boundary confinement, which is modulated by the volume of added cryosolvent (0-~50Å separation distance from protein surface). Goals are to: (1) characterize the protein-coupled solvent dynamics on correlation time scales of ~10^-10<τ<10^-7s, and spatial scales from protein surface to periphery of the surrounding solution, from the perspective of a free, small-molecule (~7Å diameter) probe, and (2) reveal properties of the solvent-protein coupling that can be correlated with protein functions, that are measureable under the same conditions. Rotational mobility of the nitroxide spin probe, TEMPOL, resolves and tracks two solvent components, the protein-associated domain (PAD; akin to hydration layer) and surrounding mesodomain, through their distinct temperature- and confinement-dependent values of τ and normalized weight. Detailed protocols are described for simulation of two-component nitroxide EPR spectra, which are categorized by line shape regime and guided by a library of template spectra and simulation parameters derived from two model soluble globular proteins. The order-disorder transition in the PAD, which is a universal feature of protein-coupled solvent dynamics, provides a well-defined, tunable property for elucidating mechanism in solvent-protein-function dynamical coupling. The low-temperature mesodomain system and EPR spin probe method are generally applicable to reveal solvent contributions to a broad range of macromolecule-mediated biological processes.

Effects of cryo-EM cooling on structural ensembles

Structure determination by cryo electron microscopy (cryo-EM) provides information on structural heterogeneity and ensembles at atomic resolution. To obtain cryo-EM images of macromolecules, the samples are first rapidly cooled down to cryogenic temperatures. To what extent the structural ensemble is perturbed during cooling is currently unknown. Here, to quantify the effects of cooling, we combined continuum model calculations of the temperature drop, molecular dynamics simulations of a ribosome complex before and during cooling with kinetic models. Our results suggest that three effects markedly contribute to the narrowing of the structural ensembles: thermal contraction, reduced thermal motion within local potential wells, and the equilibration into lower free-energy conformations by overcoming separating free-energy barriers. During cooling, barrier heights below 10 kJ/mol were found to be overcome, which is expected to reduce B-factors in ensembles imaged by cryo-EM. Our approach now enables the quantification of the heterogeneity of room-temperature ensembles from cryo-EM structures.

The rapid temperature drop during plunge-freezing affects the structural ensembles obtained by cryo-EM. To quantify the extent of perturbation, Bock and Grubmüller combined continuum calculations, MD simulations, and kinetic models.

Probing Small-Angle Molecular Motions with EPR Spectroscopy: Dynamical Transition and Molecular Packing in Disordered Solids

Disordered molecular solids present a rather broad class of substances of different origin—amorphous polymers, materials for photonics and optoelectronics, amorphous pharmaceutics, simple molecular glass formers, and others. Frozen biological media in many respects also may be referred to this class. Theoretical description of dynamics and structure of disordered solids still does not exist, and only some phenomenological models can be developed to explain results of particular experiments. Among different experimental approaches, electron paramagnetic resonance (EPR) applied to spin probes and labels also can deliver useful information. EPR allows probing small-angle orientational molecular motions (molecular librations), which intrinsically are inherent to all molecular solids. EPR is employed in its conventional continuous wave (CW) and pulsed—electron spin echo (ESE)—versions. CW EPR spectra are sensitive to dynamical librations of molecules while ESE probes stochastic molecular librations. In this review, different manifestations of small-angle motions in EPR of spin probes and labels are discussed. It is shown that CW-EPR-detected dynamical librations provide information on dynamical transition in these media, similar to that explored with neutron scattering, and ESE-detected stochastic librations allow elucidating some features of nanoscale molecular packing. The possible EPR applications are analyzed for gel-phase lipid bilayers, for biological membranes interacting with proteins, peptides and cryoprotectants, for supercooled ionic liquids (ILs) and supercooled deep eutectic solvents (DESs), for globular proteins and intrinsically disordered proteins (IDPs), and for some other molecular solids.

Millisecond mix-and-quench crystallography (MMQX) enables time-resolved studies of PEPCK with remote data collection

Time-resolved crystallography of biomolecules in action has advanced rapidly as methods for serial crystallography have improved, but the large number of crystals and the complex experimental infrastructure that are required remain serious obstacles to its widespread application. Here, millisecond mix-and-quench crystallography (MMQX) has been developed, which yields millisecond time-resolved data using far fewer crystals and routine remote synchrotron data collection. To demonstrate the capabilities of MMQX, the conversion of oxaloacetic acid to phosphoenolpyruvate by phosphoenolpyruvate carboxy-kinase (PEPCK) is observed with a time resolution of 40 ms. By lowering the entry barrier to time-resolved crystallography, MMQX should enable a broad expansion in structural studies of protein dynamics.

Photo-Switching of Protein Dynamical Collectivity

We examine changes in the picosecond structural dynamics with irreversible photobleaching of red fluorescent proteins (RFP) mCherry, mOrange2 and TagRFP-T. Measurements of the protein dynamical transition using terahertz time-domain spectroscopy show in all cases an increase in the turn-on temperature in the bleached state. The result is surprising given that there is little change in the protein surface, and thus, the solvent dynamics held responsible for the transition should not change. A spectral analysis of the measurements guided by quasiharmonic calculations of the protein absorbance reveals that indeed the solvent dynamical turn-on temperature is independent of the thermal stability/photostate however the protein dynamical turn-on temperature shifts to higher temperatures. This is the first demonstration of switching the protein dynamical turn-on temperature with protein functional state. The observed shift in protein dynamical turn-on temperature relative to the solvent indicates an increase in the required mobile waters necessary for the protein picosecond motions, that is, these motions are more collective. Melting-point measurements reveal that the photobleached state is more thermally stable, and structural analysis of related RFP’s shows that there is an increase in internal water channels as well as a more uniform atomic root mean squared displacement. These observations are consistent with previous suggestions that water channels form with extended light excitation providing O2 access to the chromophore and subsequent fluorescence loss. We report that these same channels increase internal coupling enhancing thermal stability and collectivity of the picosecond protein motions. The terahertz spectroscopic characterization of the protein and solvent dynamical onsets can be applied generally to measure changes in collectivity of protein motions.

The Dynamics of Hydrated Proteins Are the Same as Those of Highly Asymmetric Mixtures of Two Glass-Formers

Customarily, the studies of dynamics of hydrated proteins are focused on the fast hydration water ν-relaxation, the slow structural α-relaxation responsible for a single glass transition, and the protein dynamic transition (PDT). Guided by the analogy with the dynamics of highly asymmetric mixtures of molecular glass-formers, we explore the possibility that the dynamics of hydrated proteins are richer than presently known. By providing neutron scattering, dielectric relaxation, calorimetry, and deuteron NMR data in two hydrated globular proteins, myoglobin and BSA, and the fibrous elastin, we show the presence of two structural α-relaxations, α1 and α2, and the hydration water ν-relaxation, all coupled together with interconnecting properties. There are two glass transition temperatures T _g ^α1and T _g ^α2 corresponding to vitrification of the α1 and α2 processes. Relaxation time τ_α2(T) of the α2-relaxation changes its Arrhenius temperature dependence to super-Arrhenius on crossing T _g ^α1 from below. The ν-relaxation responds to the two vitrifications by changing the T-dependence of its relaxation time τ_ν(T) on crossing consecutively T _g ^α2 and T _g ^α1. It generates the PDT at T _d where τ_ν(T _d) matches about five times the experimental instrument timescale τ_exp, provided that T _d > T _g ^α1. This condition is satisfied by the hydrated globular proteins considered in this paper, and the ν-relaxation is in the liquid state with τ_ν(T) having the super-Arrhenius temperature dependence. However, if T _d < T _g ^α1, the ν-relaxation fails to generate the PDT because it is in the glassy state and τ_ν(T) has Arrhenius T-dependence, as in the case of hydrated elastin. Overall, the dynamics of hydrated proteins are the same as the dynamics of highly asymmetric mixtures of glass-formers. The results from this study have expanded the knowledge of the dynamic processes and their properties in hydrated proteins, and impact on research in this area is expected.

Temperature Dependence of Structural Dynamics at the Catalytic Cofactor of [FeFe]-hydrogenase

[FeFe]-hydrogenases are nature's blueprint for efficient hydrogen turnover. Understanding their enzymatic mechanism may improve technological H₂ fuel generation. The active-site cofactor (H-cluster) consists of a [4Fe-4S] cluster ([4Fe]_H), cysteine-linked to a diiron site ([2Fe]_H) carrying an azadithiolate (adt) group, terminal cyanide and carbon monoxide ligands, and a bridging carbon monoxide (μCO) in the oxidized protein (Hox). Recently, the debate on the structure of reduced H-cluster states was intensified by the assignment of new species under cryogenic conditions. We investigated temperature effects (4-280 K) in infrared (IR) and X-ray absorption spectroscopy (XAS) data of [FeFe]-hydrogenases using fit analyses and quantum-chemical calculations. IR data from our laboratory and literature sources were evaluated. At ambient temperatures, reduced H-cluster states with a bridging hydride (μH^-, in Hred and Hsred) or with an additional proton at [4Fe]_H (Hred') or at the distal iron of [2Fe]_H (Hhyd) prevail. At cryogenic temperatures, these species are largely replaced by states that hold a μCO, lack [4Fe]_H protonation, and bind an additional proton at the adt nitrogen (HredH⁺ and HsredH⁺). XAS revealed the atomic coordinate dispersion (i.e., the Debye-Waller parameter, 2σ²) of the iron-ligand bonds and Fe-Fe distances in the oxidized and reduced H-cluster. 2σ² showed a temperature dependence typical for the so-called protein-glass transition, with small changes below ∼200 K and a pronounced increase above this "breakpoint". This behavior is attributed to the freezing-out of larger-scale anharmonic motions of amino acid side chains and water species. We propose that protonation at [4Fe]_H as well as ligand rearrangement and μH^- binding at [2Fe]_H are impaired because of restricted molecular mobility at cryogenic temperatures so that protonation can be biased toward adt. We conclude that a H-cluster with a μCO, selective [4Fe]_H or [2Fe]_H protonation, and catalytic proton transfer via adt facilitates efficient H₂ conversion in [FeFe]-hydrogenase.

Collective excitations in α-helical protein structures interacting with the water environment

Low-frequency vibrational excitations of protein macromolecules in the terahertz frequency region are suggested to contribute to many biological processes such as enzymatic catalysis, intra-protein energy/charge transport, recognition, and allostery. To explain high effectiveness of these processes, two possible mechanisms of the long-lived excitation were proposed by H. Fröhlich and A.S. Davydov, which relate to either vibrational modes or solitary waves, respectively. In this paper, we developed a quantum dynamic model of vibrational excitation in α-helical proteins interacting with the aqueous environment. In the model, we distinguished three coupled subsystems, i.e., (i) a chain of hydrogen-bonded peptide groups (PGs), interacting with (ii) the subsystem of the side-chain residuals which in turn interact with (iii) the environment, surrounding water responsible for dissipation and fluctuation in the system. It was shown that the equation of motion for phonon variables of the PG chain can be transformed to nonlinear Schrodinger equation which admits bifurcation into the solution corresponding to the weak-damped vibrational modes (Fröhlich-type regime) and Davydov solitons. A bifurcation parameter is derived through the strength of phonon-phonon interaction between the side-chains and hydration-shell water molecules. As shown, the energy of these excited states is pumped through the interaction of the side-chains with fluctuating water environment of the proteins. The suggested mechanism of the collective vibrational mode excitation is discussed in connection with the recent experiments on the long-lived collective protein excitations in the terahertz frequency region and vibrational energy transport pathways in proteins.

Control of electron transfer by protein dynamics in photosynthetic reaction centers

Trehalose and glycerol are low molecular mass sugars/polyols that have found widespread use in the protection of native protein states, in both short- and long-term storage of biological materials, and as a means of understanding protein dynamics. These myriad uses are often attributed to their ability to form an amorphous glassy matrix. In glycerol, the glass is formed only at cryogenic temperatures, while in trehalose, the glass is formed at room temperature, but only upon dehydration of the sample. While much work has been carried out to elucidate a mechanistic view of how each of these matrices interact with proteins to provide stability, rarely have the effects of these two independent systems been directly compared to each other. This review aims to compile decades of research on how different glassy matrices affect two types of photosynthetic proteins: (i) the Type II bacterial reaction center from Rhodobacter sphaeroides and (ii) the Type I Photosystem I reaction center from cyanobacteria. By comparing aggregate data on electron transfer, protein structure, and protein dynamics, it appears that the effects of these two distinct matrices are remarkably similar. Both seem to cause a "tightening" of the solvation shell when in a glassy state, resulting in severely restricted conformational mobility of the protein and associated water molecules. Thus, trehalose appears to be able to mimic, at room temperature, nearly all of the effects on protein dynamics observed in low temperature glycerol glasses.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

Advances in studying protein disorder with solid-state NMR

Solution NMR is a key tool to study intrinsically disordered proteins (IDPs), whose importance for biological function is widely accepted. However, disordered proteins are not limited to solution and are also found in non-soluble systems such as fibrils and membrane proteins. In this Trends article, I will discuss how solid-state NMR can be used to study disorder in non-soluble proteins. Techniques based on dipolar couplings can study static protein disorder which either occurs naturally as e.g. in spider silk or can be induced by freeze trapping IDPs or unfolded proteins. In this case, structural ensembles are directly reflected by a static distribution of dihedral angels that can be determined by the distribution of chemical shifts or other methods. Techniques based on J-couplings can detect dynamic protein disorder under MAS. In this case, only average chemical shifts are measured but disorder can be characterized with a variety of data including secondary chemical shifts, relaxation rates, paramagnetic relaxation enhancements, or residual dipolar couplings. I describe both technical aspects and examples of solid-state NMR on protein disorder and end the article with a discussion of challenges and opportunities of this emerging field.

Direct observation of ligand migration within human hemoglobin at work

Human hemoglobin is the textbook example of the stereochemistry of an allosteric protein and of the exquisite control that a protein can exert over ligand binding. However, the fundamental basis by which the protein facilitates the ligand movement remains unknown. In this study, we used cryogenic X-ray crystallography and a high-repetition pulsed laser irradiation technique to elucidate the atomic details of ligand migration processes in hemoglobin after photolysis of the bound CO. Our data clarify the distinct CO migration pathways in the individual subunits of hemoglobin and unravel the functional roles of the internal cavities and neighboring amino acid residues in ligand exit and entry. Our results also demonstrate the high gas permeability and porosity of hemoglobin, facilitating O₂ delivery.

Hemoglobin is one of the best-characterized proteins with respect to structure and function, but the internal ligand diffusion pathways remain obscure and controversial. Here we captured the CO migration processes in the tense (T), relaxed (R), and second relaxed (R2) quaternary structures of human hemoglobin by crystallography using a high-repetition pulsed laser technique at cryogenic temperatures. We found that in each quaternary structure, the photodissociated CO molecules migrate along distinct pathways in the α and β subunits by hopping between the internal cavities with correlated side chain motions of large nonpolar residues, such as α14Trp(A12), α105Leu(G12), β15Trp(A12), and β71Phe(E15). We also observe electron density evidence for the distal histidine [α58/β63His(E7)] swing-out motion regardless of the quaternary structure, although less evident in α subunits than in β subunits, suggesting that some CO molecules have escaped directly through the E7 gate. Remarkably, in T-state Fe(II)-Ni(II) hybrid hemoglobins in which either the α or β subunits contain Ni(II) heme that cannot bind CO, the photodissociated CO molecules not only dock at the cavities in the original Fe(II) subunit, but also escape from the protein matrix and enter the cavities in the adjacent Ni(II) subunit even at 95 K, demonstrating the high gas permeability and porosity of the hemoglobin molecule. Our results provide a comprehensive picture of ligand movements in hemoglobin and highlight the relevance of cavities, nonpolar residues, and distal histidines in facilitating the ligand migration.

Solvent flows, conformation changes and lattice reordering in a cold protein crystal

When protein crystals are abruptly cooled, the unit-cell, protein and solvent-cavity volumes all contract, but the volume of bulk-like internal solvent may expand. Outflow of this solvent from the unit cell and its accumulation in defective interior crystal regions has been suggested as one cause of the large increase in crystal mosaicity on cooling. It is shown that when apoferritin crystals are abruptly cooled to temperatures between 220 and 260 K, the unit cell contracts, solvent is pushed out and the mosaicity grows. On temperature-dependent timescales of 10 to 200 s, the unit-cell and solvent-cavity volume then expand, solvent flows back in, and the mosaicity and B factor both drop. Expansion and reordering at fixed low temperature are associated with small-amplitude but large-scale changes in the conformation and packing of apoferritin. These results demonstrate that increases in mosaicity on cooling arise due to solvent flows out of or into the unit cell and to incomplete, arrested relaxation of protein conformation. They indicate a critical role for time in variable-temperature crystallographic studies, and the feasibility of probing interactions and cooperative conformational changes that underlie cold denaturation in the presence of liquid solvent at temperatures down to ∼200 K.

Lipid chain mobility and packing in DOPC bilayers at cryogenic temperatures

Low-temperature molecular mobility and packing in biological tissues are important for their survival upon cryopreservation. Electron paramagnetic resonance (EPR) in its pulsed version of electron spin echo (ESE) allows studying stochastic librations of spin-labeled molecules, the type of motion which dominates at low temperatures. These librations are characterized by the parameter <α²>τ_c where <α²> is the mean squared angular amplitude and τ_c is the correlation time for the motion. This parameter is known to be larger for higher temperature and for looser intermolecular structure. In this work, ESE data for the bilayers comprised of doubly-unsaturated DOPC (dioleoyl-glycero-phosphocholine) lipids and mono-unsaturated POPC (palmitoyl-oleoyl-glycero-phosphocholine) lipids with spin-labeled stearic acids added were obtained in the temperature range between 80 and 210 K; the results were compared also with the previously obtained data for fully-saturated DPPC (dipalmitoyl-glycero-phosphocholine) lipid bilayers [J. Phys. Chem. B2014, 118, 12,478-12,485; Appl. Magn. Reson. 2018, 49, 1369-1383]. It turned out that for DOPC bilayers the <α²>τ_c values are of intermediate magnitude between those for POPC and DPPC bilayers, which implies an intermediate density of lipid packing. A possible explanation of this result could be rearrangement at cryogenic temperatures of the DOPC lipid tails, with their terminal segments folding cooperatively. This interpretation is also in agreement with the known thermodynamic properties of gel-fluid transition for DOPC bilayer.

Solid-State Biology and Seed Longevity: A Mechanical Analysis of Glasses in Pea and Soybean Embryonic Axes

The cytoplasm of anhydrobiotes (organisms that persist in the absence of water) solidifies during drying. Despite this stabilization, anhydrobiotes vary in how long they persist while dry. In this paper, we call upon concepts currently used to explain stability of amorphous solids (also known as glasses) in synthetic polymers, foods, and pharmaceuticals to the understand variation in longevity of biological systems. We use embryonic axes of pea (Pisum sativum) and soybean (Glycine max) seeds as test systems that have relatively long and short survival times, respectively, but similar geometries and water sorption behaviors. We used dynamic mechanical analysis to gain insights on structural and mobility properties that relate to stability of other organic systems with controlled composition. Changes of elastic and loss moduli, and the dissipation function, tan δ, in response to moisture and temperature were compared in pea and soybean axes containing less than 0.2 g H2O g-1 dry mass. The work shows high complexity of structure-regulated molecular mobility within dried seed matrices. As was previously observed for pea cotyledons, there were multiple relaxations of structural constraints to molecular movement, which demonstrate substantial localized, "fast" motions within solidified cytoplasm. There was detected variation in the coordination among long-range slow, diffusive and short-range fast, vibrational motions in glasses of pea compared to soybean, which suggest higher constraints to motion in pea and greater "fragility" in soybean. We are suggesting that differences in fragility contribute to variation of seed longevity. Indeed, fragility and coordination of short and long range motions are linked to stability and physical aging of synthetic polymers.

Ice formation and solvent nanoconfinement in protein crystals

Ice formation within protein crystals is a major obstacle to the cryocrystallographic study of protein structure, and has limited studies of how the structural ensemble of a protein evolves with temperature in the biophysically interesting range from ∼260 K to the protein-solvent glass transition near 200 K. Using protein crystals with solvent cavities as large as ∼70 Å, time-resolved X-ray diffraction was used to study the response of protein and internal solvent during rapid cooling. Solvent nanoconfinement suppresses freezing temperatures and ice-nucleation rates so that ice-free, low-mosaicity diffraction data can be reliably collected down to 200 K without the use of cryoprotectants. Hexagonal ice (Ih) forms in external solvent, but internal crystal solvent forms stacking-disordered ice (Isd) with a near-random stacking of cubic and hexagonal planes. Analysis of powder diffraction from internal ice and single-crystal diffraction from the host protein structure shows that the maximum crystallizable solvent fraction decreases with decreasing crystal solvent-cavity size, and that an ∼6 Å thick layer of solvent adjacent to the protein surface cannot crystallize. These results establish protein crystals as excellent model systems for the study of nanoconfined solvent. By combining fast cooling, intense X-ray beams and fast X-ray detectors, complete structural data sets for high-value targets, including membrane proteins and large complexes, may be collected at ∼220-240 K that have much lower mosaicities and comparable B factors, and that may allow more confident identification of ligand binding than in current cryocrystallographic practice.

Motions of water and solutes-Slaving versus plasticization phenomena

It is well-accepted that hydration water is crucial for the structure, dynamics, and function of proteins. However, the exact role of water for the motions and functions of proteins is still debated. Experiments have shown that protein and water dynamics are strongly coupled but with water motions occurring on a considerably faster time scale (the so-called slaving behavior). On the other hand, water also reduces the conformational entropy of proteins and thereby acts as a plasticizer of them. In this work, we analyze the dynamics (using broadband dielectric spectroscopy) of some specific non-biological water solutions in a broad concentration range to elucidate the role of water in the dynamics of the solutes. Our results demonstrate that at low water concentrations (less than 5 wt. %), the plasticization phenomenon prevails for all the materials analyzed. However, at higher water concentrations, two different scenarios can be observed: the slaving phenomenon or plasticization, depending on the solute analyzed. These results generalize the slaving phenomenon to some, but not all, non-biological solutions and allow us to analyze the key factors for observing the slaving behavior in protein solutions as well as to reshaping the slaving concept.

Dynamical transition in molecular glasses and proteins observed by spin relaxation of nitroxide spin probes and labels

In glassy substances and biological media, dynamical transitions are observed in neutron scattering that manifests itself as deviations of the translational mean-squared displacement, 〈x²〉, of hydrogen atoms from harmonic dynamics. In biological media, the deviation occurs at two temperature intervals, at ∼100-150 K and at ∼170-230 K, and it is attributed to the motion of methyl groups in the former case and to the transition from harmonic to anharmonic or diffusive motions in the latter case. In this work, electron spin echo (ESE) spectroscopy-a pulsed version of electron paramagnetic resonance-is applied to study the spin relaxation of nitroxide spin probes and labels introduced in molecular glass former o-terphenyl and in protein lysozyme. The anisotropic contribution to the rate of the two-pulse ESE decay, ΔW, is induced by spin relaxation appearing because of restricted orientational stochastic molecular motion; it is proportional to 〈α²〉τ_c, where 〈α²〉 is the mean-squared angle of reorientation of the nitroxide molecule around the equilibrium position and τ_c is the correlation time of reorientation. The ESE time window allows us to study motions with τ_c < 10^-7 s. For glassy o-terphenyl, the 〈α²〉τ_c temperature dependence shows a transition near 240 K, which is in agreement with the literature data on 〈x²〉. For spin probes of essentially different size, the obtained data were found to be close, which evidences that motion is cooperative, involving a nanocluster of several neighboring molecules. For the dry lysozyme, the 〈α²〉τ_c values below 260 K were found to linearly depend on the temperature in the same way as it was observed in neutron scattering for 〈x²〉. As spin relaxation is influenced only by stochastic motion, the harmonic motions seen in ESE must be overdamped. In the hydrated lysozyme, ESE data show transitions near 130 K for all nitroxides, near 160 K for the probe located in the hydration layer, and near 180 K for the label in the protein interior. For this system, the two latter transitions are not observed in neutron scattering. The ESE-detected transitions are suggested to be related with water dynamics in the nearest hydration shell: with water glass transition near 130 K and with the onset of overall water molecular reorientations near 180 K; the disagreement with neutron scattering is ascribed to the larger time window for ESE-detected motions.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

Low-Temperature Dynamics of Chain-Labeled Lipids in Ester- and Ether-Linked Phosphatidylcholine Membranes

Continuous wave electron paramagnetic resonance spectroscopy and two-pulse echo detected spectra of chain-labeled lipids are used to study the dynamics of frozen lipid membranes over the temperature range 77-260 K. Bilayers of ester-linked dihexadecanoylphosphatidylcholine (DPPC) with noninterdigitated chains and ether-linked dihexadecyl phosphatidylcholine (DHPC) with interdigitated chains are considered. Rapid stochastic librations of small angular amplitude are found in both lipid matrices. In noninterdigitated DPPC bilayers, the mean-square angular amplitude, [Formula: see text], of the motion increases with temperature and it is larger close to the chain termini than close to the polar/apolar interface. In contrast, in interdigitated DHPC lamellae, [Formula: see text] is small and temperature and label-position independent at low temperature and increases steeply at high temperature. The rotational correlation time, τ_c, of librations lies in the subnanosecond range for DPPC and in the nanosecond range for DHPC. In all membrane samples, the temperature dependence of [Formula: see text] resembles that of the mean-square atomic displacement revealed by neutron scattering and a dynamical transition is detected in the range 210-240 K. The results highlight the librational oscillations and the glass-like behavior in bilayer and interdigitated lipid membranes.

Solvent-Slaved Dynamic Processes Observed by Tryptophan Phosphorescence of Human Serum Albumin

Despite extensive experimental and computational efforts to understand the nature of the hierarchy of protein fluctuations and the modulating role of the protein hydration shell, a detailed microscopic description of the dynamics of the protein-solvent system has yet to be achieved. By using single tryptophan protein phosphorescence, we follow site-specific internal protein dynamics over a broad temperature range and demonstrate three independent dynamic processes. Process I is seen at temperatures below the bulk solvent T_g, has low activation energy, and is likely due to fast vibrations that may be enabled by water mobility on the protein surface. Process II is observed above 170 K, with activation energy typical of β relaxations in a glass; it has the same temperature dependence as fluctuations of hydration shell waters. Process III is observed at T > 200 K; it has super-Arrhenius temperature dependence and closely follows the primary relaxation of the bulk. The fluorescence of pyranine bound to the protein reports on the mobility of water in the hydration shell; it reveals a shift in emission spectra with increasing temperature, indicative of a changing H-bond network at the surface of the protein. These results support a model of solvent-slaved protein dynamics.

Gd3+-Gd3+ distances exceeding 3 nm determined by very high frequency continuous wave electron paramagnetic resonance

Electron paramagnetic resonance spectroscopy in combination with site-directed spin labeling is a very powerful tool for elucidating the structure and organization of biomolecules. Gd3+ complexes have recently emerged as a new class of spin labels for distance determination by pulsed EPR spectroscopy at Q- and W-band. We present CW EPR measurements at 240 GHz (8.6 Tesla) on a series of Gd-rulers of the type Gd-PyMTA-spacer-Gd-PyMTA, with Gd-Gd distances ranging from 1.2 nm to 4.3 nm. CW EPR measurements of these Gd-rulers show that significant dipolar broadening of the central |-1/2〉 → |1/2〉 transition occurs at 30 K for Gd-Gd distances up to ∼3.4 nm with Gd-PyMTA as the spin label. This represents a significant extension for distances accessible by CW EPR, as nitroxide-based spin labels at X-band frequencies can typically only access distances up to ∼2 nm. We show that this broadening persists at biologically relevant temperatures above 200 K, and that this method is further extendable up to room temperature by immobilizing the sample in glassy trehalose. We show that the peak-to-peak broadening of the central transition follows the expected 1/r3 dependence for the electron-electron dipolar interaction, from cryogenic temperatures up to room temperature. A simple procedure for simulating the dependence of the lineshape on interspin distance is presented, in which the broadening of the central transition is modeled as an S = 1/2 spin whose CW EPR lineshape is broadened through electron-electron dipolar interactions with a neighboring S = 7/2 spin.

Role of Hydration Layer in Dynamical Transition in Proteins: Insights from Translational Self-Diffusivity

Elucidation of the role of hydration water underpinning dynamical crossover in proteins has proven challenging. Indeed, many contradictory findings in the literature seek to establish either causal or correlative links between water and protein behavior. Here, via molecular dynamics, we compute the temperature dependence of mean-square displacement and translational self-diffusivities for both hen egg white lysozyme and its hydration layer from 190 to 300 K. We find that the protein's mobility increases sharply at ∼230 K, indicating dynamical onset; concerted motion with hydration-water molecules is evident up to ∼285 K, confirming dynamical correlation between them. Exploring underlying mechanisms of such concerted motion, we scrutinize the water-protein hydrogen-bonding network as a function of temperature, noting sharp deviation from linearity of the hydrogen bond number's profile with temperature originating near the protein dynamical transition. Our studies reveal a common temperature profile/dependence of self-diffusivity values of the protein, hydration water, and the bulk solvent, originating from a common dependence on the bulk solvent viscosity, η_S. The key mechanistic role adopted by the protein-water hydrogen bond network in relation to the onset of proteins' dynamical transition is also discussed.

Evidence of Coupling between the Motions of Water and Peptides

Studies of protein dynamics at low temperatures are generally performed on hydrated powders and not in biologically realistic solutions of water because of water crystallization. However, here we avoid the problem of crystallization by reducing the size of the biomolecules. We have studied oligomers of the amino acid l-lysine, fully dissolved in water, and our dielectric relaxation data show that the glass transition-related dynamics of the oligomers is determined by the water dynamics, in a way similar to that previously observed for solvated proteins. This implies that the crucial role of water for protein dynamics can be extended to other types of macromolecular systems, where water is also able to determine their conformational fluctuations. Using the energy landscape picture of macromolecules, the thermodynamic criterion for such solvent-slaved macromolecular motions may be that the macromolecules need the entropy contribution from the solvent to overcome the enthalpy barriers between different conformational substates.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Dynamics of hydrated proteins and bio-protectants: Caged dynamics, β-relaxation, and α-relaxation

BACKGROUND

The properties of the three dynamic processes, α-relaxation, ν-relaxation, and caged dynamics in aqueous mixtures and hydrated proteins are analogous to corresponding processes found in van der Waals and polymeric glass-formers apart from minor differences.

METHODS

Collection of various experimental data enables us to characterize the structural α-relaxation of the protein coupled to hydration water (HW), the secondary or ν-relaxation of HW, and the caged HW process.

RESULTS

From the T-dependence of the ν-relaxation time of hydrated myoglobin, lysozyme, and bovine serum albumin, we obtain Ton at which it enters the experimental time windows of Mössbauer and neutron scattering spectroscopies, coinciding with protein dynamical transition (PDT) temperature Td. However, for all systems considered, the α-relaxation time at Ton or Td is many orders of magnitude longer. The other step change of the mean-square-displacement (MSD) at Tg_alpha originates from the coupling of the nearly constant loss (NCL) of caged HW to density. The coupling of the NCL to density is further demonstrated by another step change at the secondary glass temperature Tg_beta in two bio-protectants, trehalose and sucrose.

CONCLUSIONS

The structural α-relaxation plays no role in PDT. Since PDT is simply due to the ν-relaxation of HW, the term PDT is a misnomer. NCL of caged dynamics is coupled to density and show transitions at lower temperature, Tg_beta and Tg_alpha.

GENERAL SIGNIFICANCE

The so-called protein dynamical transition (PDT) of hydrated proteins is not caused by the structural α-relaxation of the protein but by the secondary ν-relaxation of hydration water. "This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".

Atomistic details of protein dynamics and the role of hydration water

BACKGROUND

The importance of protein dynamics for their biological activity is now well recognized. Different experimental and computational techniques have been employed to study protein dynamics, hierarchy of different processes and the coupling between protein and hydration water dynamics. Yet, understanding the atomistic details of protein dynamics and the role of hydration water remains rather limited.

SCOOP OF REVIEW

Based on overview of neutron scattering, molecular dynamic simulations, NMR and dielectric spectroscopy results we present a general picture of protein dynamics covering time scales from faster than ps to microseconds and the influence of hydration water on different relaxation processes.

MAJOR CONCLUSIONS

Internal protein dynamics spread over a wide time range from faster than picosecond to longer than microseconds. We suggest that the structural relaxation in hydrated proteins appears on the microsecond time scale, while faster processes present mostly motion of side groups and some domains. Hydration water plays a crucial role in protein dynamics on all time scales. It controls the coupled protein-hydration water relaxation on 10-100ps time scale. This process defines the friction for slower protein dynamics. Analysis suggests that changes in amount of hydration water affect not only general friction, but also influence significantly the protein's energy landscape.

GENERAL SIGNIFICANCE

The proposed atomistic picture of protein dynamics provides deeper understanding of various relaxation processes and their hierarchy, similarity and differences between various biological macromolecules, including proteins, DNA and RNA. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo".