Radially softening diffusive motions in a globular protein

Dynamics of side chains in poly(quinoxaline-2,3-diyl)s studied via quasielastic neutron scattering

The side chain dynamics of poly(quinoxaline-2,3-diyl)s (PQXs) are expected to influence their conformation. To investigate these dynamics experimentally, quasielastic neutron scattering (QENS) was performed for PQXs in deuterated tetrahydrofuran (THF-d8) and deuterated 1,1,2-trichloroethane/THF (1,1,2-TCE-d3/THF-d8), in which they formed right-handed and left-handed helical structures, respectively. The mean-square displacement of the PQX side chains in 1,1,2-TCE-d3/THF-d8 was lower than that in THF-d8. Furthermore, QENS complementary studies and molecular dynamics simulations unraveled a coupling between the main-chain and side chain dynamics of PQXs, suggesting the possibility of controlling the main-chain helical chirality through the dynamics of chiral side chains. These insights present a novel strategy for the design of synthetic helical macromolecules with precise chirality control.

Molecular Dynamics of Lysozyme Amyloid Polymorphs Studied by Incoherent Neutron Scattering

Lysozyme amyloidosis is a hereditary disease, which is characterized by the deposition of lysozyme amyloid fibrils in various internal organs. It is known that lysozyme fibrils show polymorphism and that polymorphs formed at near-neutral pH have the ability to promote more monomer binding than those formed at acidic pH, indicating that only specific polymorphs become dominant species in a given environment. This is likely due to the polymorph-specific configurational diffusion. Understanding the possible differences in dynamical behavior between the polymorphs is thus crucial to deepen our knowledge of amyloid polymorphism and eventually elucidate the molecular mechanism of lysozyme amyloidosis. In this study, molecular dynamics at sub-nanosecond timescale of two kinds of polymorphic fibrils of hen egg white lysozyme, which has long been used as a model of human lysozyme, formed at pH 2.7 (LP27) and pH 6.0 (LP60) was investigated using elastic incoherent neutron scattering (EINS) and quasi-elastic neutron scattering (QENS). Analysis of the EINS data showed that whereas the mean square displacement of atomic motions is similar for both LP27 and LP60, LP60 contains a larger fraction of atoms moving with larger amplitudes than LP27, indicating that the dynamical difference between the two polymorphs lies not in the averaged amplitude, but in the distribution of the amplitudes. Furthermore, analysis of the QENS data showed that the jump diffusion coefficient of atoms is larger for LP60, suggesting that the atoms of LP60 undergo faster diffusive motions than those of LP27. This study thus characterizes the dynamics of the two lysozyme polymorphs and reveals that the molecular dynamics of LP60 is enhanced compared with that of LP27. The higher molecular flexibility of the polymorph would permit to adjust its conformation more quickly than its counterpart, facilitating monomer binding.

Mechanical Unfolding of Spectrin Repeats Induces Water-Molecule Ordering

Mechanical processes are involved at many stages of the development of living cells, and often external forces applied to a biomolecule result in its unfolding. Although our knowledge of the unfolding mechanisms and the magnitude of the forces involved has evolved, the role that water molecules play in the mechanical unfolding of biomolecules has not yet been fully elucidated. To this end, we investigated with steered molecular dynamics simulations the mechanical unfolding of dystrophin's spectrin repeat 1 and related the changes in the protein's structure to the ordering of the surrounding water molecules. Our results indicate that upon mechanically induced unfolding of the protein, the solvent molecules become more ordered and increase their average number of hydrogen bonds. In addition, the unfolded structures originating from mechanical pulling expose an increasing amount of the hydrophobic residues to the solvent molecules, and the uncoiled regions adapt a convex surface with a small radius of curvature. As a result, the solvent molecules reorganize around the protein's small protrusions in structurally ordered waters that are characteristic of the so-called "small-molecule regime," which allows water to maintain a high hydrogen bond count at the expense of an increased structural order. We also determined that the response of water to structural changes in the protein is localized to the specific regions of the protein that undergo unfolding. These results indicate that water plays an important role in the mechanically induced unfolding of biomolecules. Our findings may prove relevant to the ever-growing interest in understanding macromolecular crowding in living cells and their effects on protein folding, and suggest that the hydration layer may be exploited as a means for short-range allosteric communication.

Strong Adverse Contribution of Conformational Dynamics to Streptavidin-Biotin Binding

Molecular dynamics plays an important role for the biological function of proteins. For protein ligand interactions, changes of conformational entropy of protein and hydration layer are relevant for the binding process. Quasielastic neutron scattering (QENS) was used to investigate differences in protein dynamics and conformational entropy of ligand-bound and ligand-free streptavidin. Protein dynamics were probed both on the fast picosecond time scale using neutron time-of-flight spectroscopy and on the slower nanosecond time scale using high-resolution neutron backscattering spectroscopy. We found the internal equilibrium motions of streptavidin and the corresponding mean square displacements (MSDs) to be greatly reduced upon biotin binding. On the basis of the observed MSDs, we calculated the difference of conformational entropy ΔS_conf of the protein component between ligand-bound and ligand-free streptavidin. The rather large negative ΔS_conf value (-2 kJ mol^-1 K^-1 on the nanosecond time scale) obtained for the streptavidin tetramer seems to be counterintuitive, given the exceptionally high affinity of streptavidin-biotin binding. Literature data on the total entropy change ΔS observed upon biotin binding to streptavidin, which includes contributions from both the protein and the hydration water, suggest partial compensation of the unfavorable ΔS_conf by a large positive entropy gain of the surrounding hydration layer and water molecules that are displaced during ligand binding.

Homogeneous and heterogeneous dynamics in native and denatured bovine serum albumin

Quasielastic incoherent neutron spectroscopy experiments reveal that chemical denaturation significantly modifies the internal dynamics of bovine serum albumin.

Modulation of hemoglobin dynamics by an allosteric effector

Hemoglobin (Hb) is an extensively studied paradigm of proteins that alter their function in response to allosteric effectors. Models of its action have been used as prototypes for structure‐function relationships in many proteins, and models for the molecular basis of its function have been deeply studied and extensively argued. Recent reports suggest that dynamics may play an important role in its function. Relatively little is known about the slow, correlated motions of hemoglobin subunits in various structural states because experimental and computational strategies for their characterization are challenging. Allosteric effectors such as inositol hexaphosphate (IHP) bind to both deoxy‐Hb and HbCO, albeit at different sites, leading to a lowered oxygen affinity. The manner in which these effectors impact oxygen binding is unclear and may involve changes in structure, dynamics or both. Here we use neutron spin echo measurements accompanied by wide‐angle X‐ray scattering to show that binding of IHP to HbCO results in an increase in the rate of coordinated motions of Hb subunits relative to one another with little if any change in large scale structure. This increase of large‐scale dynamics seems to be coupled with a decrease in the average magnitude of higher frequency modes of individual residues. These observations indicate that enhanced dynamic motions contribute to the functional changes induced by IHP and suggest that they may be responsible for the lowered oxygen affinity triggered by these effectors.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Hierarchical molecular dynamics of bovine serum albumin in concentrated aqueous solution below and above thermal denaturation

The dynamics of proteins in solution is a complex and hierarchical process, affected by the aqueous environment as well as temperature. We present a comprehensive study on nanosecond time and nanometer length scales below, at, and above the denaturation temperature Td. Our experimental data evidence dynamical processes in protein solutions on three distinct time scales. We suggest a consistent physical picture of hierarchical protein dynamics: (i) self-diffusion of the entire protein molecule is confirmed to agree with colloid theory for all temperatures where the protein is in its native conformational state. At higher temperatures T > Td, the self-diffusion is strongly obstructed by cross-linking or entanglement. (ii) The amplitude of backbone fluctuations grows with increasing T, and a transition in its dynamics is observed above Td. (iii) The number of mobile side-chains increases sharply at Td while their average dynamics exhibits only little variations. The combination of quasi-elastic neutron scattering and the presented analytical framework provides a detailed microscopic picture of the protein molecular dynamics in solution, thereby reflecting the changes of macroscopic properties such as cluster formation and gelation.

Effects of pressure on the dynamics of an oligomeric protein from deep-sea hyperthermophile

Inorganic pyrophosphatase (IPPase) from Thermococcus thioreducens is a large oligomeric protein derived from a hyperthermophilic microorganism that is found near hydrothermal vents deep under the sea, where the pressure is up to 100 MPa (1 kbar). It has attracted great interest in biophysical research because of its high activity under extreme conditions in the seabed. In this study, we use the quasielastic neutron scattering (QENS) technique to investigate the effects of pressure on the conformational flexibility and relaxation dynamics of IPPase over a wide temperature range. The β-relaxation dynamics of proteins was studied in the time ranges from 2 to 25 ps, and from 100 ps to 2 ns, using two spectrometers. Our results indicate that, under a pressure of 100 MPa, close to that of the native environment deep under the sea, IPPase displays much faster relaxation dynamics than a mesophilic model protein, hen egg white lysozyme (HEWL), at all measured temperatures, opposite to what we observed previously under ambient pressure. This contradictory observation provides evidence that the protein energy landscape is distorted by high pressure, which is significantly different for hyperthermophilic (IPPase) and mesophilic (HEWL) proteins. We further derive from our observations a schematic denaturation phase diagram together with energy landscapes for the two very different proteins, which can be used as a general picture to understand the dynamical properties of thermophilic proteins under pressure.

Correlation between supercoiling and conformational motions of the bacterial flagellar filament

The bacterial flagellar filament is a very large macromolecular assembly of a single protein, flagellin. Various supercoiled states of the filament exist, which are formed by two structurally different conformations of flagellin in different ratios. We investigated the correlation between supercoiling of the protofilaments and molecular dynamics in the flagellar filament using quasielastic and elastic incoherent neutron scattering on the picosecond and nanosecond timescales. Thermal fluctuations in the straight L- and R-type filaments were measured and compared to the resting state of the wild-type filament. Amplitudes of motion on the picosecond timescale were found to be similar in the different conformational states. Mean-square displacements and protein resilience on the 0.1 ns timescale demonstrate that the L-type state is more flexible and less resilient than the R-type, whereas the wild-type state lies in between. Our results provide strong support that supercoiling of the protofilaments in the flagellar filament is determined by the strength of molecular forces in and between the flagellin subunits.

Vanishing amplitude of backbone dynamics causes a true protein dynamical transition: 2H NMR studies on perdeuterated C-phycocyanin

Using a combination of H2 nuclear magnetic resonance (NMR) methods, we study internal rotational dynamics of the perdeuterated protein C-phycocyanin (CPC) in dry and hydrated states over broad temperature and dynamic ranges with high angular resolution. Separating H2 NMR signals from methyl deuterons, we show that basically all backbone deuterons exhibit highly restricted motion occurring on time scales faster than microseconds. The amplitude of this motion increases when a hydration shell exists, while it decreases upon cooling and vanishes near 175 K. We conclude that the vanishing of the highly restricted motion marks a dynamical transition, which is independent of the time window and of a fundamental importance. This conclusion is supported by results from experimental and computational studies of the proteins myoglobin and elastin. In particular, we argue based on findings in molecular dynamics simulations that the behavior of the highly restricted motion of proteins at the dynamical transition resembles that of a characteristic secondary relaxation of liquids at the glass transition, namely the nearly constant loss. Furthermore, H2 NMR studies on perdeuterated CPC reveal that, in addition to highly restricted motion, small fractions of backbone segments exhibit weakly restricted dynamics when temperature and hydration are sufficiently high.

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.

Functional domain motions in proteins on the ~1-100 ns timescale: comparison of neutron spin-echo spectroscopy of phosphoglycerate kinase with molecular-dynamics simulation

Protein function often requires large-scale domain motion. An exciting new development in the experimental characterization of domain motions in proteins is the application of neutron spin-echo spectroscopy (NSE). NSE directly probes coherent (i.e., pair correlated) scattering on the ~1-100 ns timescale. Here, we report on all-atom molecular-dynamics (MD) simulation of a protein, phosphoglycerate kinase, from which we calculate small-angle neutron scattering (SANS) and NSE scattering properties. The simulation-derived and experimental-solution SANS results are in excellent agreement. The contributions of translational and rotational whole-molecule diffusion to the simulation-derived NSE and potential problems in their estimation are examined. Principal component analysis identifies types of domain motion that dominate the internal motion's contribution to the NSE signal, with the largest being classic hinge bending. The associated free-energy profiles are quasiharmonic and the frictional properties correspond to highly overdamped motion. The amplitudes of the motions derived by MD are smaller than those derived from the experimental analysis, and possible reasons for this difference are discussed. The MD results confirm that a significant component of the NSE arises from internal dynamics. They also demonstrate that the combination of NSE with MD is potentially useful for determining the forms, potentials of mean force, and time dependence of functional domain motions in proteins.

Activity and dynamics of an enzyme, pig liver esterase, in near-anhydrous conditions

Water is widely assumed to be essential for life, although the exact molecular basis of this requirement is unclear. Water facilitates protein motions, and although enzyme activity has been demonstrated at low hydrations in organic solvents, such nonaqueous solvents may allow the necessary motions for catalysis. To examine enzyme function in the absence of solvation and bypass diffusional constraints we have tested the ability of an enzyme, pig liver esterase, to catalyze alcoholysis as an anhydrous powder, in a reaction system of defined water content and where the substrates and products are gaseous. At hydrations of 3 (±2) molecules of water per molecule of enzyme, activity is several orders-of-magnitude greater than nonenzymatic catalysis. Neutron spectroscopy indicates that the fast (≤nanosecond) global anharmonic dynamics of the anhydrous functional enzyme are suppressed. This indicates that neither hydration water nor fast anharmonic dynamics are required for catalysis by this enzyme, implying that one of the biological requirements of water may lie with its role as a diffusion medium rather than any of its more specific properties.

Neutron spin-echo studies of hemoglobin and myoglobin: multiscale internal dynamics

Neutron spin-echo spectroscopy was used to study structural fluctuations that occur in hemoglobin (Hb) and myoglobin (Mb) in solution. Using neutron spin-echo data up to a very high momentum transfer q ( approximately 0.62 A(-)(1)), we characterized the internal dynamics of these proteins at the levels of dynamic pair correlation function and self-correlation function in the time range of several picoseconds to a few nanoseconds. In the same protein solution, data transition from pair correlation motion to self-correlation motion as the momentum transfer q increases. At low q, coherent scattering dominates; at high q, observations are largely due to incoherent scattering. The low q data were interpreted in terms of an effective diffusion coefficient; on the other hand, the high q data were interpreted in terms of mean square displacements. Comparison of data from the two homologous proteins collected at different temperatures and protein concentrations was used to assess the contributions made by translational and rotational diffusion and internal modes of motion to the data. The temperature dependence of decay times can be attributed to changes in the viscosity and temperature of the solvent, as predicted by the Stokes-Einstein relationship. This is true for contributions from both diffusive and internal modes of motion, indicating an intimate relationship between the internal dynamics of the proteins and the viscosity of the solvent. Viscosity change associated with protein concentration can account for changes in diffusion observed at different concentrations, but is apparently not the only factor involved in the changes in internal dynamics observed with change in protein concentration. Data collected at high q indicate that internal modes in Mb are generally faster than those in Hb, perhaps due to the greater surface-to-volume ratio of Mb and the fact that surface groups tend to exhibit faster motion than buried groups. Comparison of data from Hb and data from Mb at low q indicates an unexpectedly rapid motion of Hb alphabeta dimers relative to one another. Dynamic motion of subunits is increasingly perceived as important to the allosteric behavior of Hb. Our data demonstrate that this motion is highly sensitive to protein concentration, temperature, and solvent viscosity, indicating that great care needs to be exercised in interpreting its effect on protein function.

Effects of co-solvents on peptide hydration water structure and dynamics

We evaluate the molecular response of hydration water as a function of temperature and proximity to the surface of the peptide N-acetyl-leucine-methylamide (NALMA) when in the presence of the kosmotrope co-solvent glycerol or the chaotrope co-solvent dimethyl sulfoxide (DMSO), using molecular dynamics simulation with a polarizable force field. These detailed microscopic studies complement established thermodynamic analysis on the role of co-solvents in shifting the equilibrium for proteins away from or towards the native folded state. We find that the structure of the water at the peptide interfaces reflects an increase in hydration number in the glycerol solution and a decrease in hydration numbers in the DMSO solution. While the water dynamics around NALMA in the presence of both co-solvents is slower than that observed with the water solvent alone, in the DMSO mixture we no longer measure a separation in water motion time scales at low temperatures as is seen in the pure water solvent, but rather one single relaxation time. In the glycerol, however, we do observe a separation of time scales at low temperatures, supporting the hypothesis that hydration water near a hydrophobic solute evolves on a separate time scale than the extensive hydrogen-bonding network of more bulk-like water. Our simulation studies highlight the differences in the two co-solvent solutions due to the relative frequency of water contacts with the hydrophobic vs. hydrophilic peptide surface, and direct water interactions with the co-solvents.

From powder to solution: hydration dependence of human hemoglobin dynamics correlated to body temperature

A transition in hemoglobin (Hb), involving partial unfolding and aggregation, has been shown previously by various biophysical methods. The correlation between the transition temperature and body temperature for Hb from different species, suggested that it might be significant for biological function. To focus on such biologically relevant human Hb dynamics, we studied the protein internal picosecond motions as a response to hydration, by elastic and quasielastic neutron scattering. Rates of fast diffusive motions were found to be significantly enhanced with increasing hydration from fully hydrated powder to concentrated Hb solution. In concentrated protein solution, the data showed that amino acid side chains can explore larger volumes above body temperature than expected from normal temperature dependence. The body temperature transition in protein dynamics was absent in fully hydrated powder, indicating that picosecond protein dynamics responsible for the transition is activated only at a sufficient level of hydration. A collateral result from the study is that fully hydrated protein powder samples do not accurately describe all aspects of protein picosecond dynamics that might be necessary for biological function.

Atomic-scale dynamics inside living cells explored by neutron scattering

Single-particle neutron spectroscopy has contributed important experimental data on molecular dynamics in biological systems. The technique provides information on atomic and molecular motions in macromolecules on the picosecond to the nanosecond time scale, which are essential to biological function. Here, we report on recent neutron measurements performed directly in living cells by using isotope labelling to explore the dynamics of specific cellular components. The paper proposes an integrated view of results on atomic-scale cell water dynamics, internal and global macromolecular motions and solvent isotope effect on macromolecular dynamics. The work established the specific usefulness of the neutron scattering technique to get insight into biologically relevant dynamical features, in particular through comparative measurements. The method developed can now be applied to look for dynamical signatures related to cell characteristics in many different cell types and organelles.

Hydration water dynamics near biological interfaces

We performed classical molecular dynamics simulations using both fixed-charge and polarizable water and protein force fields to contrast the hydration dynamics near hydrophilic and amphiphilic peptides as a function of temperature. The high peptide concentrations we use serve as a model for the surface of folded proteins where hydration layers around each residue overlap significantly. Through simulation we determine that there are notable differences in the water dynamics analyzed from the outer and inner hydration layer regions of the amphiphilic peptide solution that explains the experimentally observed presence of two translational relaxations, while the hydrophilic peptide solution shows only a single non-Arrhenius translational process with no distinction between hydration layers. Given that water dynamics for the amphiphilic peptide system reproduces all known rotational and translational hydration dynamical anomalies exhibited by hydration water near protein surfaces, our analysis provides strong evidence that dynamical signatures near biological interfaces arises because of frustration in the hydration dynamics induced by chemical heterogeneity, as opposed to just topological roughness, of the protein surface.

Hemoglobin dynamics in red blood cells: correlation to body temperature

A transition in hemoglobin behavior at close to body temperature has been discovered recently by micropipette aspiration experiments on single red blood cells (RBCs) and circular dichroism spectroscopy on hemoglobin solutions. The transition temperature was directly correlated to the body temperatures of a variety of species. In an exploration of the molecular basis for the transition, we present neutron scattering measurements of the temperature dependence of hemoglobin dynamics in whole human RBCs in vivo. The data reveal a change in the geometry of internal protein motions at 36.9 degrees C, at human body temperature. Above that temperature, amino acid side-chain motions occupy larger volumes than expected from normal temperature dependence, indicating partial unfolding of the protein. Global protein diffusion in RBCs was also measured and the findings compared favorably with theoretical predictions for short-time self-diffusion of noncharged hard-sphere colloids. The results demonstrated that changes in molecular dynamics in the picosecond time range and angstrom length scale might well be connected to a macroscopic effect on whole RBCs that occurs at body temperature.

Distinguishing thermodynamic and kinetic views of the preferential hydration of protein surfaces

Motivated by a quasi-chemical view of protein hydration, we define specific hydration sites on the surface of globular proteins in terms of the local water density at each site relative to bulk water density. The corresponding kinetic definition invokes the average residence time for a water molecule at each site and the average time that site remains unoccupied. Bound waters are identified by high site occupancies using either definition. In agreement with previous molecular dynamics simulation studies, we find only a weak correlation between local water densities and water residence times for hydration sites on the surface of two globular proteins, lysozyme and staphylococcal nuclease. However, a strong correlation is obtained when both the average residence and vacancy times are appropriately taken into account. In addition, two distinct kinetic regimes are observed for hydration sites with high occupancies: long residence times relative to vacancy times for a single water molecule, and short residence times with high turnover involving multiple water molecules. We also correlate water dynamics, characterized by average occupancy and vacancy times, with local heterogeneities in surface charge and surface roughness, and show that both features are necessary to obtain sites corresponding to kinetically bound waters.

In vivo measurement of internal and global macromolecular motions in Escherichia coli

We present direct quasielastic neutron scattering measurements, in vivo, of macromolecular dynamics in Escherichia coli. The experiments were performed on a wide range of timescales to cover the large panel of internal and self-diffusion motions. Three major internal processes were extracted at physiological temperature: a fast picosecond process that corresponded to restricted jump diffusion motions and two slower processes that resulted from reorientational motions occurring in approximately 40 ps and 90 ps, respectively. The analysis of the fast process revealed that the cellular environment leads to an appreciable increase in internal molecular flexibility and diffusive motion rates compared with those evaluated in fully hydrated powders. The result showed that the amount of cell water plays a decisive role in internal molecular dynamics. Macromolecular interactions and confinement, however, attenuate slightly the lubricating effect of water, as revealed by the decrease of the in vivo parameters compared with those measured in solution. The study demonstrated that standard sample preparations do not mimic accurately the physiological environment and suggested that intracellular complexity participates in functional dynamics necessary for biological activity. Furthermore, the method allowed the extraction of the self-diffusion of E. coli macromolecules, which presented similar parameters as those extracted for hemoglobin in red blood cells.

Biophysical study of thermal denaturation of apo-calmodulin: dynamics of native and unfolded states

Apo-calmodulin, a small, mainly alpha, soluble protein is a calcium-dependent protein activator. This article presents a study of internal dynamics of native and thermal unfolded apo-calmodulin, using quasi-elastic neutron scattering. This technique can probe protein internal dynamics in the picosecond timescale and in the nanometer length-scale. It appears that a dynamical transition is associated with thermal denaturation of apo-calmodulin. This dynamical transition goes together with a decrease of the confinement of hydrogen atoms, a decrease of immobile protons proportion and an increase of dynamical heterogeneity. The comparison of native and unfolded states dynamics suggests that the dynamics of protein atoms is more influenced by their distance to the backbone than by their solvent exposure.

A molecular view of melting in anhydrous phospholipidic membranes

A high-flux backscattering spectrometer and a time-of-flight disk chopper spectrometer are used to probe the molecular mobility of model freeze-dried phospholipid liposomes at a range of temperatures surrounding the main melting transition. Using specific deuteration, quasielastic neutron scattering provides evidence that, in contrast to the hydrocarbon chains, the headgroups of the phospholipid molecules do not exhibit a sharp melting transition. The onset of motion in the tails is located at temperatures far below the calorimetric transition. Long-range motion is achieved through the onset of whole-lipid translation at the melting temperature. Atomistic simulations are performed on a multibilayer model at conditions corresponding to the scattering experiments. The model provides a good description of the dynamics of the system, with predictions of the scattering functions that agree with experimental results. The analysis of both experimental data and results of simulations supports a picture of a gradual melting of the heterogeneous hydrophobic domain, with part of the chains spanning increasingly larger volumes and part of them remaining effectively immobile until the thermodynamic phase transition occurs.

C-phycocyanin hydration water dynamics in the presence of trehalose: an incoherent elastic neutron scattering study at different energy resolutions

We present a study of C-phycocyanin hydration water dynamics in the presence of trehalose by incoherent elastic neutron scattering. By combining data from two backscattering spectrometers with a 10-fold difference in energy resolution we extract a scattering law S(Q,omega) from the Q-dependence of the elastic intensities without sampling the quasielastic range. The hydration water is described by two dynamically different populations--one diffusing inside a sphere and the other diffusing quasifreely--with a population ratio that depends on temperature. The scattering law derived describes the experimental data from both instruments excellently over a large temperature range (235-320 K). The effective diffusion coefficient extracted is reduced by a factor of 10-15 with respect to bulk water at corresponding temperatures. Our approach demonstrates the benefits and the efficiency of using different energy resolutions in incoherent elastic neutron scattering over a large angular range for the study of biological macromolecules and hydration water.

Picosecond dynamics of T and R forms of aspartate transcarbamylase: A neutron scattering study

E. coli aspartate transcarbamylase (ATCase) is a 310 kDa allosteric enzyme which catalyses the first committed step in pyrimidine biosynthesis. The binding of its substrates, carbamylphosphate and aspartate, induces significant conformational changes. This enzyme shows homotropic cooperative interactions between the catalytic sites for the binding of aspartate. This property is explained by a quaternary structure transition from T state (aspartate low affinity) to R state (aspartate high affinity) accompanied by a 5% increase of radius of gyration of ATCase. The same quaternary structure change is observed upon binding of the bisubstrate analogue PALA (N-(phosphonacetyl)-L-aspartate. Owing to the large incoherent neutron scattering cross-section of the hydrogen atom and the abundance of this element in proteins, inelastic neutron scattering gives a global view of protein dynamics as sensed via the individual motions of its hydrogen atoms. We present neutron scattering results of the local dynamics (few angstroms), at short time (few tens of picoseconds), of ATCase in T and R forms. Compared to the T form, we observe an increased mobility of the protein in the R form that we associate to an increase of accessible surface area to the solvent. Beyond this specific result, this highlights the key role of the accessible surface area (ASA) in dynamic contribution to inelastic neutron data in the picosecond time scale. In particular, we want to stress out (i) that a difference at the picosecond time scale does not allow to conclude to a difference in the dynamics at a longer time scale and to address whether the T state is looser than the R state (ii) how challenging is, any comparison in terms of general dynamics (tense or relaxed) between dynamic values deduced from experimental neutron data on proteins with different sequences and therefore ASA. This caveat holds particularly when comparing dynamics of a mesophile with the corresponding extremophile.

How do side chains orient globally in protein structures?

An angle Omega is defined to serve as a metric for global side-chain orientations, which reflects the orientation of the side chain relative to the radial vector from the center of the protein to an amino acid. The side-chain orientations of buried residues exhibit characteristically different orientations than do exposed residues, in both monomeric and dimeric structures. Overall, buried side chains point mostly inward, whereas surface side chains tend to point outward from the surface. This difference in behavior also correlates well with the residue hydrophobicity; so a global side-chain orientation can be viewed as a direct structural manifestation of hydrophobicity. When various solvent-accessible layers are considered, the behavior is relatively continuous between centrally located and exposed residues. In the case of interfacial residues between subunits, there are statistically significant differences between exposed residues and interface residues for ALA, ARG, ASN, ASP, GLU, HIS, LYS, THR, VAL, MET, PRO, and overall the interface residues have an increased tendency to point inward. Presumably, these substantial differences in orientations of side chains may be a manifestation of hydrophobic forces.

Picosecond-time-scale fluctuations of proteins in glassy matrices: the role of viscosity

Through elastic neutron scattering we investigated the fast dynamics of lysozyme in hydrated powder form or embedded in glycerol-water and glucose-water matrices. We calculated the relaxational contribution to the mean square displacements of protein hydrogen atoms. We found that the inverse of this quantity is linearly proportional to the logarithm of the viscosity of the solvent glassy matrix. This relationship suggests a close connection between the picosecond-time-scale dynamics of protein side chains and the solvent structural relaxation.

Molecular view of water dynamics near model peptides

Incoherent quasi-elastic neutron scattering (QENS) has been used to measure the dynamics of water molecules in solutions of a model protein backbone, N-acetyl-glycine-methylamide (NAGMA), as a function of concentration, for comparison with results for water dynamics in aqueous solutions of the N-acetyl-leucine-methylamide (NALMA) hydrophobic peptide at comparable concentrations. From the analysis of the elastic incoherent structure factor, we find significant fractions of elastic intensity at high and low concentrations for both solutes, which corresponds to a greater population of protons with rotational time scales outside the experimental resolution (>13 ps). The higher-concentration solutions show a component of the elastic fraction that we propose is due to water motions that are strongly coupled to the solute motions, while for low-concentration solutions an additional component is activated due to dynamic coupling between inner and outer hydration layers. An important difference between the solute types at the highest concentration studied is found from stretched exponential fits to their experimental intermediate scattering functions, showing more pronounced anomalous diffusion signatures for NALMA, including a smaller stretched exponent beta and a longer structural relaxation time tau than those found for NAGMA. The more normal water diffusion exhibited near the hydrophilic NAGMA provides experimental support for an explanation of the origin of the anomalous diffusion behavior of NALMA as arising from frustrated interactions between water molecules when a chemical interface is formed upon addition of a hydrophobic side chain, inducing spatial heterogeneity in the hydration dynamics in the two types of regions of the NALMA peptide. We place our QENS measurements on model biological solutes in the context of other spectroscopic techniques and provide both confirming as well as complementary dynamic information that attempts to give a unifying molecular view of hydration dynamics signatures near peptides and proteins.

Quasielastic neutron scattering and relaxation processes in proteins: analytical and simulation-based models

The present article gives an overview of analytical and simulation approaches to describe the relaxation dynamics of proteins. Particularly emphasised are recent developments of theoretical models, such as fractional Brownian dynamics. The latter connects dynamical events seen on the pico- to nanosecond time scale, accessible to quasielastic neutron scattering, and functional dynamics of proteins on much longer time scales.

Protein dynamics in solution and powder measured by incoherent elastic neutron scattering: the influence of Q-range and energy resolution

Incoherent elastic neutron scattering (IENS) has been widely used to measure intramolecular atomic mean square displacements (MSDs) of proteins in powder and in solution. The instrumental energy resolution and the wave vector transfer (Q-range) determine, respectively, the time and length scales of observable motions. In order to investigate contributions of diffusive motions to MSDs measured by this method, we calculated the elastic intensity for several simple scattering functions. We showed that continuous translational diffusion contributes to MSDs in a Q-range where the energy width of the scattering function is of the order of the instrumental energy resolution. We discuss the choice of instruments adapted to focus on intramolecular motions in the presence of solvent or global macromolecular diffusion. The concepts developed are applied to interpret experimental data from H(2)O- and D(2)O-hydrated proteins. Finally, analogies between the Gaussian approximation in IENS and the Guinier approximation in small-angle scattering are discussed.

The role of dynamics in enzyme activity

Although protein function is thought to depend on flexibility, precisely how the dynamics of the molecule and its environment contribute to catalytic mechanisms is unclear. We review experimental and computational work relating to enzyme dynamics and function, including the role of solvent. The evidence suggests that fast motions on the 100 ps timescale, and any motions coupled to these, are not required for enzyme function. Proteins where the function is electron transfer, proton tunneling, or ligand binding may have different dynamical dependencies from those for enzymes, and enzymes with large turnover numbers may have different dynamical dependencies from those that turn over more slowly. The timescale differences between the fastest anharmonic fluctuations and the barrier-crossing rate point to the need to develop methods to resolve the range of motions present in enzymes on different time- and lengthscales.

Translational hydration water dynamics drives the protein glass transition

Experimental and computer simulation studies have revealed the presence of a glass-like transition in the internal dynamics of hydrated proteins at approximately 200 K involving an increase of the amplitude of anharmonic dynamics. This increase in flexibility has been correlated with the onset of protein activity. Here, we determine the driving force behind the protein transition by performing molecular dynamics simulations of myoglobin surrounded by a shell of water. A dual heat bath method is used with which, in any given simulation, the protein and solvent are held at different temperatures, and sets of simulations are performed varying the temperature of the components. The results show that the protein transition is driven by a dynamical transition in the hydration water that induces increased fluctuations primarily in side chains in the external regions of the protein. The water transition involves activation of translational diffusion and occurs even in simulations where the protein atoms are held fixed.

Kinetics of adsorption, desorption, and exchange of alpha-chymotrypsin and lysozyme on poly(ethyleneterephthalate) tracked film and track-etched membrane

Adsorption kinetics of (125)I-radiolabeled alpha-chymotrypsin at pH 8.6 was studied in a laminar regime between two walls of poly(ethyleneterephthalate) tracked films and membranes. Adsorption kinetics in the presence of solution (10 microg/mL), desorption by rinsing with buffer, and the following exchange of proteins by flowing unlabeled solution were measured. At pH 8.6, alpha-chymotrypsin is almost neutral and can be mostly removed from the film surface, contrary to positive lysozyme adsorbed at pH 7.4. Results suggest that alpha-chymotrypsin is irreversibly adsorbed in pores, while desorption and exchange occur on membrane flat faces. A method is proposed to determine adsorption kinetics in the pores. Kinetics of desorption and exchange of alpha-chymotrypsin from the film surface can be described by stretched exponential functions in the examined time domain with the same exponent, beta approximately 0.62, which does not depend also on the former adsorption duration. However, the mean residence time at the interface is about 2.5 times greater in the presence of only the buffer than that in the presence of solution. This effect could be explained by a fast exchange at the arrival of unlabeled solution for a part of the adsorbed population.

Energy resolution and dynamical heterogeneity effects on elastic incoherent neutron scattering from molecular systems

Incoherent neutron scattering is widely used to probe picosecond-nanosecond time scale dynamics of molecular systems. In systems of spatially confined atoms the relatively high intensity of elastic incoherent neutron scattering is often used to obtain a first estimate of the dynamics present. For many complex systems, however, experimental elastic scattering is difficult to interpret unambiguously using analytical dynamical models that go beyond the determination of an average mean-square displacement. To circumvent this problem a description of the scattering is derived here that encompasses a variety of analytical models in a common framework. The framework describes the time-converged part of the dynamic structure factor [the elastic incoherent scattering function (EISF)] and lends itself to practical use by explicitly incorporating effects due to the finite energy resolution of the instrument used. The dependence of the elastic scattering on wave vector is examined, and it is shown how heterogeneity in the distribution of mean-square displacements can be related to deviations of the scattering from Gaussian behavior. In this case, a correction to fourth order in the scattering vector can be used to extract the variance of the distribution of mean-square displacements. The formalism is used in a discussion of measurements on dynamics accompanying the glass transition in molecular systems. By fitting to experimental data obtained on a protein solution the present methodology is used to show how the existence of a temperature-dependent relaxation frequency can lead to a transition in the measured mean-square displacement in the absence of an EISF change.

Dynamic transition associated with the thermal denaturation of a small Beta protein

We studied the temperature dependence of the picosecond internal dynamics of an all-beta protein, neocarzinostatin, by incoherent quasielastic neutron scattering. Measurements were made between 20 degrees C and 71 degrees C in heavy water solution. At 20 degrees C, only 33% of the nonexchanged hydrogen atoms show detectable dynamics, a number very close to the fraction of protons involved in the side chains of random coil structures, therefore suggesting a rigid structure in which the only detectable diffusive movements are those involving the side chains of random coil structures. At 61.8 degrees C, although the protein structure is still native, slight dynamic changes are detected that could reflect enhanced backbone and beta-sheet side-chain motions at this higher temperature. Conversely, all internal dynamics parameters (amplitude of diffusive motions, fraction of immobile scatterers, mean-squared vibration amplitude) rapidly change during heat-induced unfolding, indicating a major loss of rigidity of the beta-sandwich structure. The number of protons with diffusive motion increases markedly, whereas the volume occupied by the diffusive motion of protons is reduced. At the half-transition temperature (T = 71 degrees C) most of backbone and beta-sheet side-chain hydrogen atoms are involved in picosecond dynamics.