Solvent dependence of dynamic transitions in protein solutions

Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling

We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.

Persistence of Methionine Side Chain Mobility at Low Temperatures in a Nine-Residue Low Complexity Peptide, as Probed by 2 H Solid-State NMR

Methionine side chains are flexible entities which play important roles in defining hydrophobic interfaces. We utilize deuterium static solid-state NMR to assess rotameric inter-conversions and other dynamic modes of the methionine in the context of a nine-residue random-coil peptide (RC9) with the low-complexity sequence GGKGMGFGL. The measurements in the temperature range of 313 to 161 K demonstrate that the rotameric interconversions in the hydrated solid powder state persist to temperatures below 200 K. Removal of solvation significantly reduces the rate of the rotameric motions. We employed 2 H NMR line shape analysis, longitudinal and rotation frame relaxation, and chemical exchange saturation transfer methods and found that the combination of multiple techniques creates a significantly more refined model in comparison with a single technique. Further, we compare the most essential features of the dynamics in RC9 to two different methionine-containing systems, characterized previously. Namely, the M35 of hydrated amyloid-β1-40 in the three-fold symmetric polymorph as well as Fluorenylmethyloxycarbonyl (FMOC)-methionine amino acid with the bulky hydrophobic group. The comparison suggests that the driving force for the enhanced methionine side chain mobility in RC9 is the thermodynamic factor stemming from distributions of rotameric populations, rather than the increase in the rate constant.

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.

Protein conformational entropy is not slaved to water

Conformational entropy can be an important element of the thermodynamics of protein functions such as the binding of ligands. The observed role for conformational entropy in modulating molecular recognition by proteins is in opposition to an often-invoked theory for the interaction of protein molecules with solvent water. The "solvent slaving" model predicts that protein motion is strongly coupled to various aspects of water such as bulk solvent viscosity and local hydration shell dynamics. Changes in conformational entropy are manifested in alterations of fast internal side chain motion that is detectable by NMR relaxation. We show here that the fast-internal side chain dynamics of several proteins are unaffected by changes to the hydration layer and bulk water. These observations indicate that the participation of conformational entropy in protein function is not dictated by the interaction of protein molecules and solvent water under the range of conditions normally encountered.

Normal Mode Analysis as a Routine Part of a Structural Investigation

Normal mode analysis (NMA) is a technique that can be used to describe the flexible states accessible to a protein about an equilibrium position. These states have been shown repeatedly to have functional significance. NMA is probably the least computationally expensive method for studying the dynamics of macromolecules, and advances in computer technology and algorithms for calculating normal modes over the last 20 years have made it nearly trivial for all but the largest systems. Despite this, it is still uncommon for NMA to be used as a component of the analysis of a structural study. In this review, we will describe NMA, outline its advantages and limitations, explain what can and cannot be learned from it, and address some criticisms and concerns that have been voiced about it. We will then review the most commonly used techniques for reducing the computational cost of this method and identify the web services making use of these methods. We will illustrate several of their possible uses with recent examples from the literature. We conclude by recommending that NMA become one of the standard tools employed in any structural study.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Ultraslow Water-Mediated Transmembrane Interactions Regulate the Activation of A2A Adenosine Receptor

Water molecules inside a G-protein coupled receptor (GPCR) have recently been spotlighted in a series of crystal structures. To decipher the dynamics and functional roles of internal water molecules in GPCR activity, we studied the A2A adenosine receptor using microsecond molecular-dynamics simulations. Our study finds that the amount of water flux across the transmembrane (TM) domain varies depending on the receptor state, and that the water molecules of the TM channel in the active state flow three times more slowly than those in the inactive state. Depending on the location in solvent-protein interface as well as the receptor state, the average residence time of water in each residue varies from ∼O(10(2)) ps to ∼O(10(2)) ns. Especially, water molecules, exhibiting ultraslow relaxation (∼O(10(2)) ns) in the active state, are found around the microswitch residues that are considered activity hotspots for GPCR function. A continuous allosteric network spanning the TM domain, arising from water-mediated contacts, is unique in the active state, underscoring the importance of slow water molecules in the activation of GPCRs.

Elastic Scattering Spectroscopy (ESS): an Instrument-Concept for Dynamics of Complex (Bio-) Systems From Elastic Neutron Scattering

A new type of neutron-scattering spectroscopy is presented that is designed specifically to measure dynamics in bio-systems that are difficult to obtain in any other way. The temporal information is largely model-free and is analogous to relaxation processes measured with dielectric spectroscopy, but provides additional spacial and geometric aspects of the underlying dynamics. Numerical simulations of the basic instrument design show the neutron beam can be highly focussed, giving efficiency gains that enable the use of small samples. Although we concentrate on continuous neutron sources, the extension to pulsed neutron sources is proposed, both requiring minimal data-treatment and being broadly analogous with dielectric spectroscopy, they will open the study of dynamics to new areas of biophysics.

Practical Considerations for Determination of Glass Transition Temperature of a Maximally Freeze Concentrated Solution

Glass transition temperature is a unique thermal characteristic of amorphous systems and is associated with changes in physical properties such as heat capacity, viscosity, electrical resistance, and molecular mobility. Glass transition temperature for amorphous solids is referred as (T g), whereas for maximally freeze concentrated solution, the notation is (T g'). This article is focused on the factors affecting determination of T g' for application to lyophilization process design and frozen storage stability. Also, this review provides a perspective on use of various types of solutes in protein formulation and their effect on T g'. Although various analytical techniques are used for determination of T g' based on the changes in physical properties associated with glass transition, the differential scanning calorimetry (DSC) is the most commonly used technique. In this article, an overview of DSC technique is provided along with brief discussion on the alternate analytical techniques for T g' determination. Additionally, challenges associated with T g' determination, using DSC for protein formulations, are discussed. The purpose of this review is to provide a practical industry perspective on determination of T g' for protein formulations as it relates to design and development of lyophilization process and/or for frozen storage; however, a comprehensive review of glass transition temperature (T g, T g'), in general, is outside the scope of this work.

Microscopic relaxations in a protein sustained down to 160K in a non-glass forming organic solvent

BACKGROUND

We have studied microscopic dynamics of a protein in carbon disulfide, a non-glass forming solvent, down to its freezing temperature of ca. 160K.

METHODS

We have utilized quasielastic neutron scattering.

RESULTS

A comparison of lysozyme hydrated with water and dissolved in carbon disulfide reveals a stark difference in the temperature dependence of the protein's microscopic relaxation dynamics induced by the solvent. In the case of hydration water, the common protein glass-forming solvent, the protein relaxation slows down in response to a large increase in the water viscosity on cooling down, exhibiting a well-known protein dynamical transition. The dynamical transition disappears in non-glass forming carbon disulfide, whose viscosity remains a weak function of temperature all the way down to freezing at just below 160K. The microscopic relaxation dynamics of lysozyme dissolved in carbon disulfide is sustained down to the freezing temperature of its solvent at a rate similar to that measured at ambient temperature.

CONCLUSIONS

Our results demonstrate that protein dynamical transition is not merely solvent-assisted, but rather solvent-induced, or, more precisely, is a reflection of the temperature dependence of the solvent's glass-forming dynamics.

GENERAL SIGNIFICANCE

We hypothesize that, if the long debated idea regarding the direct link between the microscopic relaxations and the biological activity in proteins is correct, then not only the microscopic relaxations, but also the activity, could be sustained in proteins all the way down to the freezing temperature of a non-glass forming solvent with a weak temperature dependence of its viscosity. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Slowdown of Interhelical Motions Induces a Glass Transition in RNA

RNA function depends crucially on the details of its dynamics. The simplest RNA dynamical unit is a two-way interhelical junction. Here, for such a unit--the transactivation response RNA element--we present evidence from molecular dynamics simulations, supported by nuclear magnetic resonance relaxation experiments, for a dynamical transition near 230 K. This glass transition arises from the freezing out of collective interhelical motional modes. The motions, resolved with site-specificity, are dynamically heterogeneous and exhibit non-Arrhenius relaxation. The microscopic origin of the glass transition is a low-dimensional, slow manifold consisting largely of the Euler angles describing interhelical reorientation. Principal component analysis over a range of temperatures covering the glass transition shows that the abrupt slowdown of motion finds its explanation in a localization transition that traps probability density into several disconnected conformational pools over the low-dimensional energy landscape. Upon temperature increase, the probability density pools then flood a larger basin, akin to a lakes-to-sea transition. Simulations on transactivation response RNA are also used to backcalculate inelastic neutron scattering data that match previous inelastic neutron scattering measurements on larger and more complex RNA structures and which, upon normalization, give temperature-dependent fluctuation profiles that overlap onto a glass transition curve that is quasi-universal over a range of systems and techniques.

Thermodynamic properties of bulk and confined water

The thermodynamic response functions of water display anomalous behaviors. We study these anomalous behaviors in bulk and confined water. We use nuclear magnetic resonance (NMR) to examine the configurational specific heat and the transport parameters in both the thermal stable and the metastable supercooled phases. The data we obtain suggest that there is a behavior common to both phases: that the dynamics of water exhibit two singular temperatures belonging to the supercooled and the stable phase, respectively. One is the dynamic fragile-to-strong crossover temperature (T(L) ≃ 225 K). The second, T* ∼ 315 ± 5 K, is a special locus of the isothermal compressibility K(T)(T, P) and the thermal expansion coefficient α(P)(T, P) in the P-T plane. In the case of water confined inside a protein, we observe that these two temperatures mark, respectively, the onset of protein flexibility from its low temperature glass state (T(L)) and the onset of the unfolding process (T*).

Enrichment of druggable conformations from apo protein structures using cosolvent-accelerated molecular dynamics

Here we describe the development of an improved workflow for utilizing experimental and simulated protein conformations in the structure-based design of inhibitors for anti-apoptotic Bcl-2 family proteins. Traditional structure-based approaches on similar targets are often constrained by the sparsity of available structures and difficulties in finding lead compounds that dock against flat, flexible protein-protein interaction surfaces. By employing computational docking of known small molecule inhibitors, we have demonstrated that structural ensembles derived from either accelerated MD (aMD) or MD in the presence of an organic cosolvent generally give better scores than those assessed from analogous conventional MD. Furthermore, conformations obtained from combined cosolvent aMD simulations started with the apo-Bcl-xL structure yielded better average and minimum docking scores for known binders than an ensemble of 72 experimental apo- and ligand-bound Bcl-xL structures. A detailed analysis of the simulated conformations indicates that the aMD effectively enhanced conformational sampling of the flexible helices flanking the main Bcl-xL binding groove, permitting the cosolvent acting as small ligands to penetrate more deeply into the binding pocket and shape ligand-bound conformations not evident in conventional simulations. We believe this approach could be useful for identifying inhibitors against other protein-protein interaction systems involving highly flexible binding sites, particularly for targets with less accumulated structural data.

Temperature-dependent dynamics of dry and hydrated β-casein studied by quasielastic neutron scattering

β-Casein is a component of casein micelle with amphillic nature and is recognized as a "natively disordered" protein that lacks secondary structures. In this study, the temperature and hydration effects on the dynamics of β-casein are explored by quasielastic neutron scattering (QENS). An upturn in the mean square displacement (MSD) of hydrated β-casein indicates an increase of protein flexibility at a temperature of ~225 K. Another increase in MSD at ~100 K, observed in both dry and hydrated β-casein, is ascribed to the methyl group rotations, which are not sensitive to hydration. QENS analysis in the energy domain reveals that the fraction of hydrogen atoms participating in motion in a sphere of diffusion is highly hydration dependent and increases with temperature. In the time domain analysis, a logarithmic-like decay is observed in the range of picosecond to nanosecond (β-relaxation time) in the dynamics of hydrated β-casein. This dynamical behavior has been observed in hydrated globular and oligomeric proteins. Our temperature-dependent QENS experiments provide evidence that lack of a secondary structure in β-casein results in higher flexibility in its dynamics and easier reversible thermal unfolding compared to other rigid biomolecules.

Scaling analysis of bio-molecular dynamics derived from elastic incoherent neutron scattering experiments

Numerous neutron scattering studies of bio-molecular dynamics employ a qualitative analysis of elastic scattering data and atomic mean square displacements. We provide a new quantitative approach showing that the intensity at zero energy exchange can be a rich source of information of bio-structural fluctuations on a pico- to nano-second time scale. Elastic intensity scans performed either as a function of the temperature (back-scattering) and∕or by varying the instrumental resolution (time of flight spectroscopy) yield the activation parameters of molecular motions and the approximate structural correlation function in the time domain. The two methods are unified by a scaling function, which depends on the ratio of correlation time and instrumental resolution time. The elastic scattering concept is illustrated with a dynamic characterization of alanine-dipeptide, protein hydration water, and water-coupled protein motions of lysozyme, per-deuterated c-phycocyanin (CPC) and hydrated myoglobin. The complete elastic scattering function versus temperature, momentum exchange, and instrumental resolution is analyzed instead of focusing on a single cross-over temperature of mean square displacements at the apparent onset temperature of an-harmonic motions. Our method predicts the protein dynamical transition (PDT) at Td from the collective (α) structural relaxation rates of the solvation shell as input. By contrast, the secondary (β) relaxation enhances the amplitude of fast local motions in the vicinity of the glass temperature Tg. The PDT is specified by step function in the elastic intensity leading from elastic to viscoelastic dynamic behavior at a transition temperature Td.

The effect of complex solvents on the structure and dynamics of protein solutions: The case of Lysozyme in trehalose/water mixtures

We present a Molecular Dynamics simulation study of the effect of trehalose concentration on the structure and dynamics of individual proteins immersed in trehalose/water mixtures. Hen egg-white Lysozyme is used in this study and trehalose concentrations of 0%, 10%, 20%, 30% and 100% by weight are explored. Surprisingly, we have found that changes in trehalose concentration do not change the global structural characteristics of the protein as measured by standard quantities like the mean square deviation, radius of gyration, solvent accessible surface area, inertia tensor and asphericity. Only in the limit of pure trehalose these metrics change significantly. Specifically, we found that the protein is compressed by 2% when immersed in pure trehalose. At the amino acid level there is noticeable rearrangement of the surface residues due to the change in polarity of the surrounding environment with the addition of trehalose. From a dynamic perspective, our computation of the Incoherent Intermediate Scattering Function shows that the protein slows down with increasing trehalose concentration; however, this slowdown is not monotonic. Finally, we also report in-depth results for the hydration layer around the protein including its structure, hydrogen-bonding characteristics and dynamic behavior at different length scales.

Protein surface and core dynamics show concerted hydration-dependent activation

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Hydrophobic Binding Hot Spots of Bcl-xL Protein-Protein Interfaces by Cosolvent Molecular Dynamics Simulation

Identifying binding hot spots in protein-protein interfaces is important for understanding the binding specificity and for the design of nonpeptide, small molecule inhibitors. Molecular dynamics simulation in the isopropanol/water cosolvent environment and in water was employed to investigate Bcl-xL protein, which has a highly flexible, large, and primarily hydrophobic binding site. Simulations of either the apo- or holocrystal structures of the Bcl-xL in pure water fail to generate conformations found in the cocrystal structures of Bcl-xL in complex with its binding partners due to hydrophobic collapse. In contrast, simulations in cosolvent starting either from the apo- or holocrystal structure of the Bcl-xL yield binding-site conformations similar to that found in the cocrystal structures of Bcl-xL. Hydrophobic binding hot spots identified using the conformations from the cosolvent simulations are in excellent agreement with experimental structural data of known inhibitors. Importantly, cosolvent simulations revealed the highly dynamic nature of the hydrophobic binding pockets in Bcl-xL and yielded new structural insights for the design of novel Bcl-xL small-molecule inhibitors.

Effects of water on the structure and low/high temperature stability of confined proteins

In this study well-characterized model proteins were confined in silica nanoporous matrices. Confinement of the proteins in silica matrices allowed us to explore the role of water hydrogen bonding on the structures of the proteins in a broad range of temperatures (-120 degrees C to 95 degrees C). At low temperatures confinement suppressed freezing of water, which remained in the liquid state. We obtained direct evidence that the changes in the hydrogen bonding of water induced changes in the structure of confined proteins. At high temperatures, a reduction of hydrogen bonding of water facilitated protein-silica interactions and the confined proteins underwent denaturation. However, the incorporation of the osmolyte, trehalose, reduced protein-silica interactions, and altered the hydrogen bonding of water. As a result, the high temperature thermal stability of the confined proteins was greatly improved.

Effect of solvation-related interaction on the low-temperature dynamics of proteins

The effect of solvation-related interaction on the low-temperature dynamics of proteins is studied by taking into account the desolvation barriers in the interactions of native contacts. It is found out that about the folding transition temperature, the protein folds in a cooperative manner, and the water molecules are expelled from the hydrophobic core at the final stage in the folding process. At low temperature, however, the protein would generally be trapped in many metastable conformations with some water molecules frozen inside the protein. The desolvation takes an important role in these processes. The number of frozen water molecules and that of frozen states of proteins are further analyzed with the methods based on principal component analysis (PCA) and the clustering of conformations. It is found out that both the numbers of frozen water molecules and the frozen states of the protein increase quickly below a certain temperature. Especially, the number of frozen states of the protein increases exponentially following the decrease in the temperature, which resembles the basic features of glassy dynamics. Interestingly, it is observed that the freezing of water molecules and that of protein conformations happen at almost the same temperature. This suggests that the solvation-related interaction performs an important role for the low-temperature dynamics of the model protein.

Influence of the "host-guest" interactions on the mobility of genistein/beta-cyclodextrin inclusion complex

Inclusion complexes of cyclodextrins with nonpolar drugs are a topic of current interest in pharmaceutical science, because they increase the aqueous solubility, chemical stability and bioavailability of poorly water-soluble drugs. By means of elastic incoherent neutron scattering (EINS) technique on the backscattering spectrometer IN13, we measured, in the 150-340 K temperature range, the mean square displacements (MSD) of hydrogen atoms derived from elastic spectra of inclusion complex of beta-cyclodextrin (beta-CyD) with genistein (Gen), a phytoestrogen of great interest for its antioxidant, anticarcinogenic and antiosteoporotic activities. The mobility of the complex has been compared with the single components and the physical mixture. The elastic intensity has been interpreted in terms of the double-well model in the whole temperature range. In the case of Gen, a mainly vibrational dynamics is revealed, while for the other samples, the elastic intensity becomes non-Gaussian above a temperature T(d) congruent with 220-230 K. This dynamical activation, which is represented by an Arrhenius trend, seems to be promoted by the crystallization water molecules. The dynamics of the Gen/beta-CyD inclusion complex is restricted with respect to the physical mixture, due to the action of specific "host-guest" interactions upon complexation. Finally, the dynamical response of our systems to temperature is put in relationship with their thermal stability.

Study of solvent-protein coupling effects by neutron scattering

The present work aims to characterize the dynamical behavior of proteins immersed in bio-preserving liquids and glasses. For this purpose, the protein dUTPase was chosen, while the selected solvents were glycerol, a triol, and some homologous disaccharides, i.e., trehalose, maltose, and sucrose, which are known to be very effective bio-preserving agents. The results highlight that the disaccharides show a slowing down effect on the water dynamics, which is stronger for trehalose than in the case of the other disaccharides. Furthermore, a characterization of the medium which hosts the protein is performed by using an operative definition of fragility based on the mean square displacement extracted by elastic incoherent neutron scattering, which is directly connected to Angell's kinetic fragility based on the viscosity. Finally, a study of the dynamics of the protein sequestered within the solvents is performed. The result shows that the protein dynamics is coupled with that of the surrounding matrix.

Frontiers in biotransport: water transport and hydration

Biotransport, by its nature, is concerned with the motions of molecules in biological systems while water remains as the most important and the most commonly studied molecule across all disciplines. In this review, we focus on biopreservation and thermal therapies from the perspective of water, exploring how its molecular motions, properties, kinetic, and thermodynamic transitions govern biotransport phenomena and enable preservation or controlled destruction of biological systems.

Frontiers in Biotransport: Water Transport and Hydration

Biotransport, by its nature, is concerned with the motions of molecules in biological systems while water remains as the most important and the most commonly studied molecule across all disciplines. In this review, we focus on biopreservation and thermal therapies from the perspective of water, exploring how its molecular motions, properties, kinetic, and thermodynamic transitions govern biotransport phenomena and enable preservation or controlled destruction of biological systems.

From shell to cell: neutron scattering studies of biological water dynamics and coupling to activity

An integrated picture of hydration shell dynamics and of its coupling to functional macromolecular motions is proposed from studies on a soluble protein, on a membrane protein in its natural lipid environment, and on the intracellular environment in bacteria and red blood cells. Water dynamics in multimolar salt solutions was also examined, in the context of the very slow water component previously discovered in the cytoplasm of extreme halophilic archaea. The data were obtained from neutron scattering by using deuterium labelling to focus on the dynamics of different parts of the complex systems examined.

Instantaneous normal modes and the protein glass transition

In the instantaneous normal mode method, normal mode analysis is performed at instantaneous configurations of a condensed-phase system, leading to modes with negative eigenvalues. These negative modes provide a means of characterizing local anharmonicities of the potential energy surface. Here, we apply instantaneous normal mode to analyze temperature-dependent diffusive dynamics in molecular dynamics simulations of a small protein (a scorpion toxin). Those characteristics of the negative modes are determined that correlate with the dynamical (or glass) transition behavior of the protein, as manifested as an increase in the gradient with T of the average atomic mean-square displacement at approximately 220 K. The number of negative eigenvalues shows no transition with temperature. Further, although filtering the negative modes to retain only those with eigenvectors corresponding to double-well potentials does reveal a transition in the hydration water, again, no transition in the protein is seen. However, additional filtering of the protein double-well modes, so as to retain only those that, on energy minimization, escape to different regions of configurational space, finally leads to clear protein dynamical transition behavior. Partial minimization of instantaneous configurations is also found to remove nondiffusive imaginary modes. In summary, examination of the form of negative instantaneous normal modes is shown to furnish a physical picture of local diffusive dynamics accompanying the protein glass transition.

Differences in internal dynamics of actin under different structural states detected by neutron scattering

F-actin, a helical polymer formed by polymerization of the monomers (G-actin), plays crucial roles in various aspects of cell motility. Flexibility of F-actin has been suggested to be important for such a variety of functions. Understanding the flexibility of F-actin requires characterization of a hierarchy of dynamical properties, from internal dynamics of the actin monomers through domain motions within the monomers and relative motions between the monomers within F-actin to large-scale motions of F-actin as a whole. As a first step toward this ultimate purpose, we carried out elastic incoherent neutron scattering experiments on powders of F-actin and G-actin hydrated with D(2)O and characterized the internal dynamics of F-actin and G-actin. Well established techniques and analysis enabled the extraction of mean-square displacements and their temperature dependence in F-actin and in G-actin. An effective force constant analysis with a model consisting of three energy states showed that two dynamical transitions occur at approximately 150 K and approximately 245 K, the former of which corresponds to the onset of anharmonic motions and the latter of which couples with the transition of hydration water. It is shown that behavior of the mean-square displacements is different between G-actin and F-actin, such that G-actin is "softer" than F-actin. The differences in the internal dynamics are detected for the first time between the different structural states (the monomeric state and the polymerized state). The different behavior observed is ascribed to the differences in dynamical heterogeneity between F-actin and G-actin. Based on structural data, the assignment of the differences observed in the two samples to dynamics of specific loop regions involved in the polymerization of G-actin into F-actin is proposed.

Molecular processes in biological thermosensation

Since thermal gradients are almost everywhere, thermosensation could represent one of the oldest sensory transduction processes that evolved in organisms. There are many examples of temperature changes affecting the physiology of living cells. Almost all classes of biological macromolecules in a cell (nucleic acids, lipids, proteins) can present a target of the temperature-related stimuli. This review discusses some features of different classes of temperature-sensing molecules as well as molecular and biological processes that involve thermosensation. Biochemical, structural, and thermodynamic approaches are applied in the paper to organize the existing knowledge on molecular mechanisms of thermosensation. Special attention is paid to the fact that thermosensitive function cannot be assigned to any particular functional group or spatial structure but is rather of universal nature. For instance, the complex of thermodynamic, structural, and functional features of hemoglobin family proteins suggests their possible accessory role as “molecular thermometers”.

Hydration-dependent dynamical transition in protein: protein interactions at approximately 240 K

Interprotein motions in low and fully hydrated carboxymyoglobin crystals are investigated using molecular dynamics simulation. Below approximately 240 K, the calculated dynamic structure factor exhibits a peak arising from interprotein vibration. Above approximately 240 K, the intermolecular fluctuations of the fully hydrated crystal increase drastically, whereas the low-hydration model exhibits no transition. Autocorrelation function analysis shows the transition to be dominated by the activation of diffusive intermolecular motion. The potential of mean force for the interaction remains quasiharmonic. The results indicate useful experimental avenues on protein:protein interactions to be explored using next-generation neutron sources.