A view of dynamics changes in the molten globule-native folding step by quasielastic neutron scattering

Nanosecond structural dynamics of intrinsically disordered β-casein micelles by neutron spectroscopy

β-casein undergoes a reversible endothermic self-association, forming protein micelles of limited size. In its functional state, a single β-casein monomer is unfolded, which creates a high structural flexibility, which is supposed to play a major role in preventing the precipitation of calcium phosphate particles. We characterize the structural flexibility in terms of nanosecond molecular motions, depending on the temperature by quasielastic neutron scattering. Our major questions are: Does the self-association reduce the chain flexibility? How does the dynamic spectrum of disordered caseins differ from a compactly globular protein? How does the dynamic spectrum of β-casein in solution differ from that of a protein in hydrated powder states? We report on two relaxation processes on a nanosecond and a sub-nanosecond timescale for β-casein in solution. Both processes are analyzed by Brownian oscillator model, by which the spring constant can be defined in the isotropic parabolic potential. The slower process, which is analyzed by neutron spin echo, seems a characteristic feature of the unfolded structure. It requires bulk solvent and is not seen in hydrated protein powders. The faster process, which is analyzed by neutron backscattering, has a smaller amplitude and requires hydration water, which is also observed with folded proteins in the hydrated state. The self-association had no significant influence on internal relaxation, and thus, a β-casein protein monomer flexibility is preserved in the micelle. We derive spring constants of the faster and slower motions of β-caseins in solution and compared them with those of some proteins in various states (folded or hydrated powder).

Activated nanoscale actin-binding domain motion in the catenin-cadherin complex revealed by neutron spin echo spectroscopy

As the core component of the adherens junction in cell-cell adhesion, the cadherin-catenin complex transduces mechanical tension between neighboring cells. Structural studies have shown that the cadherin-catenin complex exists as an ensemble of flexible conformations, with the actin-binding domain (ABD) of α-catenin adopting a variety of configurations. Here, we have determined the nanoscale protein domain dynamics of the cadherin-catenin complex using neutron spin echo spectroscopy (NSE), selective deuteration, and theoretical physics analyses. NSE reveals that, in the cadherin-catenin complex, the motion of the entire ABD becomes activated on nanosecond to submicrosecond timescales. By contrast, in the α-catenin homodimer, only the smaller disordered C-terminal tail of ABD is moving. Molecular dynamics (MD) simulations also show increased mobility of ABD in the cadherin-catenin complex, compared to the α-catenin homodimer. Biased MD simulations further reveal that the applied external forces promote the transition of ABD in the cadherin-catenin complex from an ensemble of diverse conformational states to specific states that resemble the actin-bound structure. The activated motion and an ensemble of flexible configurations of the mechanosensory ABD suggest the formation of an entropic trap in the cadherin-catenin complex, serving as negative allosteric regulation that impedes the complex from binding to actin under zero force. Mechanical tension facilitates the reduction in dynamics and narrows the conformational ensemble of ABD to specific configurations that are well suited to bind F-actin. Our results provide a protein dynamics and entropic explanation for the observed force-sensitive binding behavior of a mechanosensitive protein complex.

How Useful can the Voigt Profile be in Protein Folding Processes?

The analytical expression for the Voigt profile, along with its simplified forms for the Gaussian and Lorentzian dominance, is presented. The applicability of the Voigt profile in the description of anomalous diffusion phenomena, ubiquitous in different fields of science including protein folding, is discussed. It is shown that the Voigt profile is a good descriptor of the processes occurring in protein folding and in the native state. The usefulness of the Voigt profile in deriving important information of the diffusive motions in proteins from a quasielastic incoherent neutron scattering experiments is illustrated.

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.

Comparative Study of the Collective Dynamics of Proteins and Inorganic Nanoparticles

Molecular dynamics simulations of ubiquitin in water/glycerol solutions are used to test the suggestion by Karplus and coworkers that proteins in their biologically active state should exhibit a dynamics similar to 'surface-melted' inorganic nanoparticles (NPs). Motivated by recent studies indicating that surface-melted inorganic NPs are in a 'glassy' state that is an intermediate dynamical state between a solid and liquid, we probe the validity and significance of this proposed analogy. In particular, atomistic simulations of ubiquitin in solution based on CHARMM36 force field and pre-melted Ni NPs (Voter-Chen Embedded Atom Method potential) indicate a common dynamic heterogeneity, along with other features of glass-forming (GF) liquids such as collective atomic motion in the form of string-like atomic displacements, potential energy fluctuations and particle displacements with long range correlations ('colored' or 'pink' noise), and particle displacement events having a power law scaling in magnitude, as found in earthquakes. On the other hand, we find the dynamics of ubiquitin to be even more like a polycrystalline material in which the α-helix and β-sheet regions of the protein are similar to crystal grains so that the string-like collective atomic motion is concentrated in regions between the α-helix and β-sheet domains.

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.

Dynamical Behavior of Human α-Synuclein Studied by Quasielastic Neutron Scattering

α-synuclein (αSyn) is a protein consisting of 140 amino acid residues and is abundant in the presynaptic nerve terminals in the brain. Although its precise function is unknown, the filamentous aggregates (amyloid fibrils) of αSyn have been shown to be involved in the pathogenesis of Parkinson's disease, which is a progressive neurodegenerative disorder. To understand the pathogenesis mechanism of this disease, the mechanism of the amyloid fibril formation of αSyn must be elucidated. Purified αSyn from bacterial expression is monomeric but intrinsically disordered in solution and forms amyloid fibrils under various conditions. As a first step toward elucidating the mechanism of the fibril formation of αSyn, we investigated dynamical behavior of the purified αSyn in the monomeric state and the fibril state using quasielastic neutron scattering (QENS). We prepared the solution sample of 9.5 mg/ml purified αSyn, and that of 46 mg/ml αSyn in the fibril state, both at pD 7.4 in D2O. The QENS experiments on these samples were performed using the near-backscattering spectrometer, BL02 (DNA), at the Materials and Life Science Facility at the Japan Accelerator Research Complex, Japan. Analysis of the QENS spectra obtained shows that diffusive global motions are observed in the monomeric state but largely suppressed in the fibril state. However, the amplitude of the side chain motion is shown to be larger in the fibril state than in the monomeric state. This implies that significant solvent space exists within the fibrils, which is attributed to the αSyn molecules within the fibrils having a distribution of conformations. The larger amplitude of the side chain motion in the fibril state than in the monomeric state implies that the fibril state is entropically favorable.

Nanoscale protein dynamics: a new frontier for neutron spin echo spectroscopy

Recent studies show that neutron spin echo spectroscopy (NSE) can reveal long-range protein domain motions on nanometer lengthscales and on nanosecond to microsecond timescales. This unique capability of NSE provides new opportunities to understand protein dynamics and functions, such as how binding signals are propagated in a protein to distal sites. Here we review our applications of NSE to the study of nanoscale protein domain motions in a set of cell signaling proteins. We summarize the theoretical framework we have developed, which allows one to interpret the NSE data (Biophys. J. 99, 3473 (2010) and Proc. Natl. Acad. Sci. USA 102, 17646 (2005)). Our theoretical framework uses simple concepts from nonequilibrium statistical mechanics, and does not require elaborate molecular dynamics simulations, complex fits to rotational motion, or elastic network models. It is thus more robust than multiparameter techniques that require untestable assumptions. We also demonstrate our experimental scheme involving deuterium labeling of a protein domain or a subunit in a protein complex. We show that our selective deuteration scheme can highlight and resolve specific domain dynamics from the abundant global translational and rotational motions in a protein. Our approach thus clears significant hurdles to the application of NSE for the study of protein dynamics in solution.

Dynamics of thermodynamically stable, kinetically trapped, and inhibitor-bound states of pepsin

The pepsin folding mechanism involves a prosegment (PS) domain that catalyzes folding, which is then removed, resulting in a kinetically trapped native state. Although native pepsin (Np) is kinetically stable, it is irreversibly denatured due to a large folding barrier, and in the absence of the PS it folds to a more thermodynamically stable denatured state, termed refolded pepsin (Rp). This system serves as a model to understand the nature of kinetic barriers and folding transitions between compact states. Quasielastic neutron scattering (QENS) was used to characterize and compare the flexibility of Np, as a kinetically trapped state, with that of Rp, as a thermodynamically stable fold. Additionally, the dynamics of Np were compared with those of a partially unfolded form and a thermally stabilized, inhibitor-bound form. QENS revealed length-scale-dependent differences between Np and Rp on a picosecond timescale and indicated greater flexibility in Np, leading to the conclusion that kinetic stabilization likely does not correspond to reduced internal dynamics. Furthermore, large differences were observed upon inhibition, indicating that QENS of proteins in solution may prove useful for examining the role of conformational entropy changes in ligand binding.

Activation of nanoscale allosteric protein domain motion revealed by neutron spin echo spectroscopy

NHERF1 is a multidomain scaffolding protein that assembles signaling complexes, and regulates the cell surface expression and endocytic recycling of a variety of membrane proteins. The ability of the two PDZ domains in NHERF1 to assemble protein complexes is allosterically modulated by the membrane-cytoskeleton linker protein ezrin, whose binding site is located as far as 110 Ångstroms away from the PDZ domains. Here, using neutron spin echo (NSE) spectroscopy, selective deuterium labeling, and theoretical analyses, we reveal the activation of interdomain motion in NHERF1 on nanometer length-scales and on submicrosecond timescales upon forming a complex with ezrin. We show that a much-simplified coarse-grained model suffices to describe interdomain motion of a multidomain protein or protein complex. We expect that future NSE experiments will benefit by exploiting our approach of selective deuteration to resolve the specific domain motions of interest from a plethora of global translational and rotational motions. Our results demonstrate that the dynamic propagation of allosteric signals to distal sites involves changes in long-range coupled domain motions on submicrosecond timescales, and that these coupled motions can be distinguished and characterized by NSE.

Proteins move! Protein dynamics and long-range allostery in cell signaling

An emerging point of view in protein chemistry is that proteins are not the static objects that are displayed in textbooks but are instead dynamic actors. Protein dynamics plays a fundamental role in many diseases, and spans a large hierarchy of timescales, from picoseconds to milliseconds or even longer. Nanoscale protein domain motion on length scales comparable to protein dimensions is key to understanding how signals are relayed through multiple protein-protein interactions. A canonical example is how the scaffolding proteins NHERF1 and ezrin work in coordination to assemble crucial membrane complexes. As membrane-cytoskeleton scaffolding proteins, these provide excellent prototypes for understanding how regulatory signals are relayed through protein-protein interactions between the membrane and the cytoskeleton. Here, we review recent progress in understanding the structure and dynamics of the interaction. We describe recent novel applications of neutron spin echo spectroscopy to reveal the dynamic propagation of allosteric signals by nanoscale protein motion, and present a guide to the future study of dynamics and its application to the cure of disease.

Ezrin induces long-range interdomain allostery in the scaffolding protein NHERF1

Scaffolding proteins are molecular switches that control diverse signaling events. The scaffolding protein Na(+)/H(+) exchanger regulatory factor 1 (NHERF1) assembles macromolecular signaling complexes and regulates the macromolecular assembly, localization, and intracellular trafficking of a number of membrane ion transport proteins, receptors, and adhesion/antiadhesion proteins. NHERF1 begins with two modular protein-protein interaction domains-PDZ1 and PDZ2-and ends with a C-terminal (CT) domain. This CT domain binds to ezrin, which, in turn, interacts with cytosekeletal actin. Remarkably, ezrin binding to NHERF1 increases the binding capabilities of both PDZ domains. Here, we use deuterium labeling and contrast variation neutron-scattering experiments to determine the conformational changes in NHERF1 when it forms a complex with ezrin. Upon binding to ezrin, NHERF1 undergoes significant conformational changes in the region linking PDZ2 and its CT ezrin-binding domain, as well as in the region linking PDZ1 and PDZ2, involving very long range interactions over 120 A. The results provide a structural explanation, at mesoscopic scales, of the allosteric control of NHERF1 by ezrin as it assembles protein complexes. Because of the essential roles of NHERF1 and ezrin in intracellular trafficking in epithelial cells, we hypothesize that this long-range allosteric regulation of NHERF1 by ezrin enables the membrane-cytoskeleton to assemble protein complexes that control cross-talk and regulate the strength and duration of signaling.

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.

Deconstructing the native state: energy landscapes, function, and dynamics of globular proteins

Proteins are highly complex molecules with features exquisitely selected by nature to carry out essential biological functions. Physical chemistry and polymer physics provide us with the tools needed to make sense of this complexity. Upon translation, many proteins fold to a thermodynamically stable form known as the native state. The native state is not static, but consists of a hierarchy of conformations, that are continuously explored through dynamics. In this review we provide a brief introduction to some of the core concepts required in the discussion of the protein native dynamics using energy landscapes ideas. We first discuss recent works which have challenged the structure-function paradigm by demonstrating function in disordered proteins. Next we examine the hierarchical organization in the energy landscapes using atomistic molecular dynamics simulations and principal component analysis. In particular, the role of direct and water-mediated contacts in sculpting the landscape is elaborated. Another approach to studying the native state ensemble is based on choosing high-resolution order parameters for computing one- or two-dimensional free energy surfaces. We demonstrate that 2D free energy surfaces provide rich thermodynamic and kinetic information about the native state ensemble. Brownian dynamics simulations on such a surface indicate that protein conformational dynamics is weakly activated. Finally, we briefly discuss implicit and coarse-grained protein models and emphasize the solvent role in determining native state structure and dynamics.

Native and unfolded cytochrome c--comparison of dynamics using 2D-IR vibrational echo spectroscopy

Unfolded vs native CO-coordinated horse heart cytochrome c (h-cyt c) and a heme axial methionine mutant cyt c552 from Hydrogenobacter thermophilus ( Ht-M61A) are studied by IR absorption spectroscopy and ultrafast 2D-IR vibrational echo spectroscopy of the CO stretching mode. The unfolding is induced by guanidinium hydrochloride (GuHCl). The CO IR absorption spectra for both h-cyt c and Ht-M61A shift to the red as the GuHCl concentration is increased through the concentration region over which unfolding occurs. The spectra for the unfolded state are substantially broader than the spectra for the native proteins. A plot of the CO peak position vs GuHCl concentration produces a sigmoidal curve that overlays the concentration-dependent circular dichroism (CD) data of the CO-coordinated forms of both Ht-M61A and h-cyt c within experimental error. The coincidence of the CO peak shift curve with the CD curves demonstrates that the CO vibrational frequency is sensitive to the structural changes induced by the denaturant. 2D-IR vibrational echo experiments are performed on native Ht-M61A and on the protein in low- and high-concentration GuHCl solutions. The 2D-IR vibrational echo is sensitive to the global protein structural dynamics on time scales from subpicosecond to greater than 100 ps through the change in the shape of the 2D spectrum with time (spectral diffusion). At the high GuHCl concentration (5.1 M), at which Ht-M61A is essentially fully denatured as judged by CD, a very large reduction in dynamics is observed compared to the native protein within the approximately 100 ps time window of the experiment. The results suggest the denatured protein may be in a glassy-like state involving hydrophobic collapse around the heme.

Direct observation of fast protein conformational switching

Folded proteins can exist in multiple conformational substates. Each substate reflects a local minimum on the free-energy landscape with a distinct structure. By using ultrafast 2D-IR vibrational echo chemical-exchange spectroscopy, conformational switching between two well defined substates of a myoglobin mutant is observed on the approximately 50-ps time scale. The conformational dynamics are directly measured through the growth of cross peaks in the 2D-IR spectra of CO bound to the heme active site. The conformational switching involves motion of the distal histidine/E helix that changes the location of the imidazole side group of the histidine. The exchange between substates changes the frequency of the CO, which is detected by the time dependence of the 2D-IR vibrational echo spectrum. These results demonstrate that interconversion between protein conformational substates can occur on very fast time scales. The implications for larger structural changes that occur on much longer time scales are discussed.

Terahertz spectroscopy of bacteriorhodopsin and rhodopsin: similarities and differences

We studied the low-frequency terahertz spectroscopy of two photoactive protein systems, rhodopsin and bacteriorhodopsin, as a means to characterize collective low-frequency motions in helical transmembrane proteins. From this work, we found that the nature of the vibrational motions activated by terahertz radiation is surprisingly similar between these two structurally similar proteins. Specifically, at the lowest frequencies probed, the cytoplasmic loop regions of the proteins are highly active; and at the higher terahertz frequencies studied, the extracellular loop regions of the protein systems become vibrationally activated. In the case of bacteriorhodopsin, the calculated terahertz spectra are compared with the experimental terahertz signature. This work illustrates the importance of terahertz spectroscopy to identify vibrational degrees of freedom which correlate to known conformational changes in these proteins.

Fundamental and biotechnological applications of neutron scattering measurements for macromolecular dynamics

To explore macromolecular dynamics on the picosecond timescale, we used neutron spectroscopy. First, molecular dynamics were analyzed for the hyperthermophile malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, the lactate dehydrogenase from Oryctolagus cunniculus muscle. Hyperthermophiles have elaborate molecular mechanisms of adaptation to extremely high temperature. Using a novel elastic neutron scattering approach that provides independent measurements of the global flexibility and of the structural resilience (rigidity), we have demonstrated that macromolecular dynamics represents one of these molecular mechanisms of thermoadaptation. The flexibilities were found to be similar for both enzymes at their optimal activity temperature and the resilience is higher for the hyperthermophilic protein. Secondly, macromolecular motions were examined in a native and immobilized dihydrofolate reductase (DHFR) from Escherichia coli. The immobilized mesophilic enzyme has increased stability and decreased activity, so that its properties are changed to resemble those of the thermophilic enzyme. Are these changes reflected in dynamical behavior? For this study, we performed quasielastic neutron scattering measurements to probe the protein motions. The residence time is 7.95 ps for the native DHFR and 20.36 ps for the immobilized DHFR. The average height of the potential barrier to local motions is therefore increased in the immobilized DHFR, with a difference in activation energy equal to 0.54 kcal/mol, which is, using the theoretical rate equation, of the same order than expected from calculation.

Dynamics of immobilized and native Escherichia coli dihydrofolate reductase by quasielastic neutron scattering

The internal dynamics of native and immobilized Escherichia coli dihydrofolate reductase (DHFR) have been examined using incoherent quasielastic neutron scattering. These results reveal no difference between the high frequency vibration mean-square displacement of the native and the immobilized E. coli DHFR. However, length-scale-dependent, picosecond dynamical changes are found. On longer length scales, the dynamics are comparable for both DHFR samples. On shorter length scales, the dynamics is dominated by local jump motions over potential barriers. The residence time for the protons to stay in a potential well is tau = 7.95 +/- 1.02 ps for the native DHFR and tau = 20.36 +/- 1.80 ps for the immobilized DHFR. The average height of the potential barrier to the local motions is increased in the immobilized DHFR, and may increase the activation energy for the activity reaction, decreasing the rate as observed experimentally. These results suggest that the local motions on the picosecond timescale may act as a lubricant for those associated with DHFR activity occurring on a slower millisecond timescale. Experiments indicate a significantly slower catalytic reaction rate for the immobilized E. coli DHFR. However, the immobilization of the DHFR is on the exterior of the enzyme and essentially distal to the active site, thus this phenomenon has broad implications for the action of drugs distal to the active site.

Temperature-dependent structural and functional features of a hyperthermostable enzyme using elastic neutron scattering

The dynamic behavior of an endoglucanase from the hyperthermophilic microorganism Pyrococcus furiosus was investigated using elastic neutron scattering. The temperature dependence of the atomic motions was correlated with conformational and functional characteristics of the enzyme. The onset of biological function at temperatures higher than approximately 25 degrees C (the hyperthermostable enzyme is essentially inactive at room temperature) was associated with a dynamical transition in the anharmonic motions domain. This transition from the nonactive to the enzymatically active conformation involved structurally similar conformational substates in the energy landscape. From the mean-square displacement of the protein atoms, the molecular flexibility and the effective force constants were calculated at different temperature zones. The results showed that the activity increases at higher temperatures where the intramolecular bonds are weakened and the overall rigidity of the protein is decreased. Further temperature increase resulted in significantly increased atomic fluctuations featuring heat denaturation of the protein.

Molecular mechanism underlying the thermal stability and pH-induced unfolding of CHABII

The 37-residue alpha/beta protein CHABII was previously demonstrated to undergo a gradual pH-induced unfolding. It has been shown that even at pH 4.0 CHABII still retained a highly native-like secondary structure and tertiary topology although its tight side-chain packing was severely disrupted, typical of the molten globule state. Here, we have expressed and refolded the recombinant proteins of CHABII and its mutant [Phe21]-CHABII, and subsequently conducted extensive CD and NMR characterizations. The results indicated: (1) replacement of His21 by Phe in [Phe21]-CHABII eliminated the pH-induced unfolding from pH 6.5 to 4.0, indicating that His21 was responsible for the observed pH-induced unfolding of CHABII. Further examinations revealed that although the pH-induced unfolding of CHABII was also triggered by the protonation of the His residue as previously uncovered for apomyoglobin, their molecular mechanisms are different. (2) Monitoring the pH-induced unfolding by 1H-15N HSQC spectroscopy allowed us to visualize the gradual development of the CHABII molten globule. At pH 4.0, the HSQC spectrum of CHABII was poorly dispersed with dispersions of approximately 1 ppm over proton dimension and 10 ppm over 15N dimension, characteristic of severely or even "completely unfolded" proteins. One the other hand, unambiguous assignments of the NOESY spectra of CHABII led to the identification of the persistent medium and long-range NOEs at pH 4.0, which define a highly native-like secondary structure and tertiary packing. This implies that the degree of the native-like topology might be underestimated in the previous characterization of partially folded and even completely unfolded proteins. (3) Replacement of His21 by Phe with higher side-chain hydrophobicity only caused a minor structural rearrangement but considerably enhanced the packing interaction of the hydrophobic core, as evident from a dramatic increase in NOE contacts in [Phe21]-CHABII. The enhancement led to an increase of the thermal stability of [Phe21]-CHABII by approximately 17 deg. C.

Cis/trans isomerization in HIV-1 capsid protein catalyzed by cyclophilin A: insights from computational and theoretical studies

A network of protein vibrations has recently been identified in the enzyme cyclophilin A (CypA) that is associated with its peptidyl-prolyl cis/trans isomerization activity of small peptide substrates. It has been suggested that this network may have a role in promoting the catalytic step during the isomerization reaction. This work presents the results from the characterization of this network during the isomerization of the Gly89-Pro90 peptide bond in the N-terminal domain of the capsid protein (CA(N)) from human immunodeficiency virus type 1 (HIV-1), which is a naturally occurring, biologically relevant protein substrate for CypA. A variety of computational and theoretical studies are utilized to investigate the protein dynamics of the CypA-CA(N) complex, at multiple time scales, during the isomerization step. The results provide insights into the detailed mechanism of isomerization and confirm the presence of previously reported network of protein vibrations coupled to the reaction. Conserved CypA residues at the complex interface and at positions distal to the interface form parts of this network. There is HIV-1 related medical interest in CypA; incorporation of CypA, complexed with the capsid protein, into the virion is required for the infectious activity of HIV-1. Interaction energy and dynamical cross-correlation calculations are used for a detailed investigation of the protein-protein interactions in the CypA-CA(N) complex. The results show that CA(N) residues His87-Ala-Gly-Pro-Ile-Ala92 form the majority of the interactions with CypA residues. New protein-protein interactions distal to the active site (CypA Arg148-CA(N) Gln95 and CypA Arg148-CA(N) Asn121) are also identified.

Adaptation to extreme environments: macromolecular dynamics in bacteria compared in vivo by neutron scattering

Mean macromolecular dynamics was quantified in vivo by neutron scattering in psychrophile, mesophile, thermophile and hyperthermophile bacteria. Root mean square atomic fluctuation amplitudes determining macromolecular flexibility were found to be similar for each organism at its physiological temperature ( approximately 1 A in the 0.1 ns timescale). Effective force constants determining the mean macromolecular resilience were found to increase with physiological temperature from 0.2 N/m for the psychrophiles, which grow at 4 degrees C, to 0.6 N/m for the hyperthermophiles (85 degrees C), indicating that the increase in stabilization free energy is dominated by enthalpic rather than entropic terms. Larger resilience allows macromolecular stability at high temperatures, while maintaining flexibility within acceptable limits for biological activity.

Hydration dynamics near a model protein surface

The evolution of water dynamics from dilute to very high concentration solutions of a prototypical hydrophobic amino acid with its polar backbone, N-acetyl-leucine-methylamide (NALMA), is studied by quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulation for both the completely deuterated and completely hydrogenated leucine monomer. The NALMA-water system and the QENS data together provide a unique study for characterizing the dynamics of different hydration layers near a prototypical hydrophobic side chain and the backbone of which it is attached. We observe several unexpected features in the dynamics of these biological solutions under ambient conditions. The NALMA dynamics shows evidence of de Gennes narrowing, an indication of coherent long timescale structural relaxation dynamics. The translational and rotational water dynamics at the highest solute concentrations are found to be highly suppressed as characterized by long residential time and slow diffusion coefficients. The analysis of the more dilute concentration solutions models the first hydration shell with the 2.0 M spectra. We find that for outer layer hydration dynamics that the translational diffusion dynamics is still suppressed, although the rotational relaxation time and residential time are converged to bulk-water values. Molecular dynamics analysis of the first hydration shell water dynamics shows spatially heterogeneous water dynamics, with fast water motions near the hydrophobic side chain, and much slower water motions near the hydrophilic backbone. We discuss the hydration dynamics results of this model protein system in the context of protein function and protein-protein recognition.

A measure of conformational entropy change during thermal protein unfolding using neutron spectroscopy

Thermal unfolding of proteins at high temperatures is caused by a strong increase of the entropy change which lowers Gibbs free energy change of the unfolding transition (DeltaG(unf) = DeltaH - TDeltaS). The main contributions to entropy are the conformational entropy of the polypeptide chain itself and ordering of water molecules around hydrophobic side chains of the protein. To elucidate the role of conformational entropy upon thermal unfolding in more detail, conformational dynamics in the time regime of picoseconds was investigated with neutron spectroscopy. Confined internal structural fluctuations were analyzed for alpha-amylase in the folded and the unfolded state as a function of temperature. A strong difference in structural fluctuations between the folded and the unfolded state was observed at 30 degrees C, which increased even more with rising temperatures. A simple analytical model was used to quantify the differences of the conformational space explored by the observed protein dynamics for the folded and unfolded state. Conformational entropy changes, calculated on the basis of the applied model, show a significant increase upon heating. In contrast to indirect estimates, which proposed a temperature independent conformational entropy change, the measurements presented here, demonstrated that the conformational entropy change increases with rising temperature and therefore contributes to thermal unfolding.

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.

Fast dynamics of halophilic malate dehydrogenase and BSA measured by neutron scattering under various solvent conditions influencing protein stability

Protein thermal dynamics was evaluated by neutron scattering for halophilic malate dehydrogenase from Haloarcula marismortui (HmMalDH) and BSA under different solvent conditions. As a measure of thermal stability in each case, loss of secondary structure temperatures were determined by CD. HmMalDH requires molar salt and has different stability behavior in H(2)O, D(2)O, and in NaCl and KCl solvents. BSA remains soluble in molar NaCl. The neutron experiments provided values of mean-squared atomic fluctuations at the 0.1 ns time scale. Effective force constants, characterizing the mean resilience of the protein structure, were calculated from the variation of the mean-squared fluctuation with temperature. For HmMalDH, resilience increased progressively with increasing stability, from molar NaCl in H(2)O, via molar KCl in D(2)O, to molar NaCl in D(2)O. Surprisingly, however, the opposite was observed for BSA; its resilience is higher in H(2)O where it is less stable than in D(2)O. These results confirmed the complexity of dynamics-stability relationships in different proteins. Softer dynamics for BSA in D(2)O showed that the higher thermostability is associated with entropic fluctuations. In the halophilic protein, higher stability is associated with increased resilience showing the dominance of enthalpic terms arising from bonded interactions. From previous data, it is suggested that these are associated with hydrated ion binding stabilizing the protein in the high-salt solvent.

Dynamic regimes and correlated structural dynamics in native and denatured alpha-lactalbumin

Understanding the mechanisms of protein folding requires knowledge of both the energy landscape and the structural dynamics of a protein. We report a neutron-scattering study of the nanosecond and picosecond dynamics of native and the denatured alpha-lactalbumin. The quasielastic scattering intensity shows that there are alpha-helical structure and tertiary-like side-chain interactions fluctuating on sub-nanosecond time-scales under extremely denaturing conditions and even in the absence of disulfide bonds. Based on the length-scale dependence of the decay rate of the measured correlation functions, the nanosecond dynamics of the native and the variously denatured proteins have three dynamic regimes. When 0.05<Q<0.5 A(-1) (where the scattering vector, Q, is inversely proportional to the length-scale), the decay rate, Gamma, shows a power law relationship, Gamma proportional to Q(2.42+/-0.08), that is analogous to the dynamic behavior of a random coil. However, when 0.5<Q<1.0 A(-1), the decay rate exhibits a Gamma proportional to Q(1.0+/-0.2) relationship. The effective diffusion constant of the protein decreases with increasing Q, a striking dynamic behavior that is not found in any chain-like macromolecule. We suggest that this unusual dynamics is due to the presence of a strongly attractive force and collective conformational fluctuations in both the native and the denatured states of the protein. Above Q>1.0 A(-1) is a regime that displays the local dynamic behavior of individual residues, Gamma proportional to Q(1.8+/-0.3). The picosecond time-scale dynamics shows that the potential barrier to side-chain proton jump motion is reduced in the molten globule and in the denatured proteins when compared to that of the native protein. Our results provide a dynamic view of the native-like topology established in the early stages of protein folding.

Radially softening diffusive motions in a globular protein

Molecular dynamics simulation, quasielastic neutron scattering and analytical theory are combined to characterize diffusive motions in a hydrated protein, C-phycocyanin. The simulation-derived scattering function is in approximate agreement with experiment and is decomposed to determine the essential contributions. It is found that the geometry of the atomic motions can be modeled as diffusion in spheres with a distribution of radii. The time dependence of the dynamics follows stretched exponential behavior, reflecting a distribution of relaxation times. The average side chain and backbone dynamics are quantified and compared. The dynamical parameters are shown to present a smooth variation with distance from the core of the protein. Moving outward from the center of the protein there is a progressive increase of the mean sphere size, accompanied by a narrowing and shifting to shorter times of the relaxation time distribution. This smooth, "radially softening" dynamics may have important consequences for protein function. It also raises the possibility that the dynamical or "glass" transition with temperature observed experimentally in proteins might be depth dependent, involving, as the temperature decreases, progressive freezing out of the anharmonic dynamics with increasing distance from the center of the protein.

Nmr probes of molecular dynamics: overview and comparison with other techniques

NMR spin relaxation spectroscopy is a powerful approach for characterizing intramolecular and overall rotational motions in proteins. This review describes experimental methods for measuring laboratory frame spin relaxation rate constants by high-resolution solution-state NMR spectroscopy, together with theoretical approaches for interpreting spin relaxation data in order to quantify protein conformational dynamics on picosecond-nanosecond time scales. Recent applications of these techniques to proteins are surveyed, and investigations of the contribution of conformational chain entropy to protein function are highlighted. Insights into the dynamical properties of proteins obtained from NMR spin relaxation spectroscopy are compared with results derived from other experimental and theoretical techniques.