Dynamics of different functional parts of bacteriorhodopsin: H-2H labeling and neutron scattering

Quantifying Bacteriorhodopsin Activity as a Function of its Local Environment with a Raman-Based Assay

Bacteriorhodopsin (bR) is a transmembrane protein that functions as a light-driven proton pump in halophilic archaea. The bR photocycle has been well-characterized; however, these measurements almost exclusively measured purified bR, outside of its native membrane. To investigate what effect the cellular environment has on the bR photocycle, we have developed a Raman-based assay that can monitor the activity of the bR in a variety of conditions, including in its native membrane. The assay uses two continuous-wave lasers, one to initiate photochemistry and one to monitor bR activity. The excitation leads to the steady-state depletion of ground-state bR, which directly relates to the population of photocycle intermediate states. We have used this assay to monitor bR activity both in vitro and in vivo. Our in vitro measurements confirm that our assay is sensitive to bulk environmental changes reported in the literature. Our in vivo measurements show a decrease in bR activity with increasing extracellular pH for bR in its native membrane. The difference in activity with increasing pH indicates that the native membrane environment affects the function of bR. This assay opens the door to future measurements into understanding how the local environment of this transmembrane protein affects function.

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure-function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Lastly, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

Computational analysis of protein conformational heterogeneity

In this paper, we applied the molecular dynamics (MD) simulations and used thermolysin as the system to study the overall protein dynamics and side chain dihedral angles across the Arrhenius break. Simulations were performed at a high temperature of 36 °C which is above the previously observed Arrhenius break, and the lower temperature of 20 °C which is below the Arrhenius break. We observed different protein dynamics and conformational heterogeneity of side chain dihedral angles of thermolysin at the two temperatures. Our results indicated that certain regions of thermolysin have a higher level of fluctuation at lower temperature. A temperature dependent dihedral angles were also observed at the two temperatures. The changes observed in the protein indicated key areas of temperature sensitivity within the protein.Communicated by Ramaswamy H. Sarma.

Dynamic Coupling of Tyrosine 185 with the Bacteriorhodopsin Photocycle, as Revealed by Chemical Shifts, Assisted AF-QM/MM Calculations and Molecular Dynamic Simulations

Aromatic residues are highly conserved in microbial photoreceptors and play crucial roles in the dynamic regulation of receptor functions. However, little is known about the dynamic mechanism of the functional role of those highly conserved aromatic residues during the receptor photocycle. Tyrosine 185 (Y185) is one of the highly conserved aromatic residues within the retinal binding pocket of bacteriorhodopsin (bR). In this study, we explored the molecular mechanism of its dynamic coupling with the bR photocycle by automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) calculations and molecular dynamic (MD) simulations based on chemical shifts obtained by 2D solid-state NMR correlation experiments. We observed that Y185 plays a significant role in regulating the retinal cis–trans thermal equilibrium, stabilizing the pentagonal H-bond network, participating in the orientation switch of Schiff Base (SB) nitrogen, and opening the F42 gate by interacting with the retinal and several key residues along the proton translocation channel. Our findings provide a detailed molecular mechanism of the dynamic couplings of Y185 and the bR photocycle from a structural perspective. The method used in this paper may be applied to the study of other microbial photoreceptors.

A Minireview on Temperature Dependent Protein Conformational Sampling

In this minireview we discuss the role of the more subtle conformational change-protein conformational sampling and connect it to the classic relationship of protein structure and function. The theory of pre-existing functional states of protein are discussed in context of alternate protein conformational sampling. Last, we discuss how temperature, ligand binding and mutations affect the protein conformational sampling mode which is linked to the protein function regulation. The review includes several protein systems that showed temperature dependent protein conformational sampling. We also specifically included two enzyme systems, thermophilic alcohol dehydrogenase (ht-ADH) and thermolysin which we previously studied when discussing temperature dependent protein conformational sampling.

Dynamics of Bacteriorhodopsin in the Dark-Adapted State from Solution Nuclear Magnetic Resonance Spectroscopy

To achieve efficient proton pumping in the light-driven proton pump bacteriorhodopsin (bR), the protein must be tightly coupled to the retinal to rapidly convert retinal isomerization into protein structural rearrangements. Methyl group dynamics of bR embedded in lipid nanodiscs were determined in the dark-adapted state, and were found to be mostly well ordered at the cytosolic side. Methyl groups in the M145A mutant of bR, which displays only 10 % residual proton pumping activity, are less well ordered, suggesting a link between side-chain dynamics on the cytosolic side of the bR cavity and proton pumping activity. In addition, slow conformational exchange, attributed to low frequency motions of aromatic rings, was indirectly observed for residues on the extracellular side of the bR cavity. This may be related to reorganization of the water network. These observations provide a detailed picture of previously undescribed equilibrium dynamics on different time scales for ground-state bR.

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.

Multiscale nuclear magnetic relaxation dispersion of complex liquids in bulk and confinement

The nuclear magnetic relaxation dispersion (NMRD) technique consists of measurement of the magnetic-field dependence of the longitudinal nuclear-spin-lattice relaxation rate 1/T₁. Usually, the acquisition of the NMRD profiles is made using a fast field cycling (FFC) NMR technique that varies the magnetic field and explores a very large range of Larmor frequencies (10 kHz < ω₀/(2π) <40 MHz). This allows extensive explorations of the fluctuations to which nuclear spin relaxation is sensitive. The FFC technique thus offers opportunities on multiple scales of both time and distance for characterizing the molecular dynamics and transport properties of complex liquids in bulk or embedded in confined environments. This review presents the principles, theories and applications of NMRD for characterizing fundamental properties such as surface correlation times, diffusion coefficients and dynamical surface affinity (NMR wettability) for various confined liquids. The basic longitudinal and transverse relaxation equations are outlined for bulk liquids. The nuclear relaxation of a liquid confined in pores is considered in detail in order to find the biphasic fast exchange relations for a liquid at proximity of a solid surface. The physical-chemistry of liquids at solid surfaces induces striking differences between NMRD profiles of aprotic and protic (water) liquids embedded in calibrated porous disordered materials. A particular emphasis of this review concerns the extension of FFC NMR relaxation to industrial applications. For instance, it is shown that the FFC technique is sufficiently rapid for following the progressive setting of cement-based materials (plasters, cement pastes, concretes). The technique also allows studies of the dynamics of hydrocarbons in proximity of asphaltene nano-aggregates and macro-aggregates in heavy crude oils as a function of the concentration of asphaltenes. It also gives new information on the wettability of petroleum fluids (brine and oil) embedded in shale oil rocks. It is useful for understanding the relations and correlations between NMR relaxation times T₁ and T₂, diffusion coefficients D, and viscosity η of heavy crude oils. This is of particular importance for interpreting T₁, T₂, 2D T₁-T₂ and D-T₂ correlation spectra that could be obtained down-hole, thus giving a valuable tool for investigating in situ the molecular dynamics of petroleum fluids. Another domain of interest concerns biological applications. This is of particular importance for studying the complex dynamical spectrum of a folded polymeric structure that may span many decades in frequency or time. A direct NMRD characterization of water diffusional dynamics is presented at the protein interface. NMR experiments using a shuttle technique give results well above the frequency range accessible via the FFC technique; examples of this show protein dynamics over a range from fast and localized motions to slow and delocalized collective motions involving the whole protein. This review ends by an interpretation of the origin of the proton magnetic field dependence of T₁ for different biological tissues of animals; this includes a proposal for interpreting in vivo MRI data from human brain at variable magnetic fields, where the FFC relaxation analysis suggests that brain white-matter is distinct from grey-matter, in agreement with diffusion-weighted MRI imaging.

The fluctuating ribosome: thermal molecular dynamics characterized by neutron scattering

Conformational changes associated with ribosome function have been identified by X-ray crystallography and cryo-electron microscopy. These methods, however, inform poorly on timescales. Neutron scattering is well adapted for direct measurements of thermal molecular dynamics, the ‘lubricant’ for the conformational fluctuations required for biological activity. The method was applied to compare water dynamics and conformational fluctuations in the 30 S and 50 S ribosomal subunits from Haloarcula marismortui, under high salt, stable conditions. Similar free and hydration water diffusion parameters are found for both subunits. With respect to the 50 S subunit, the 30 S is characterized by a softer force constant and larger mean square displacements (MSD), which would facilitate conformational adjustments required for messenger and transfer RNA binding. It has been shown previously that systems from mesophiles and extremophiles are adapted to have similar MSD under their respective physiological conditions. This suggests that the results presented are not specific to halophiles in high salt but a general property of ribosome dynamics under corresponding, active conditions. The current study opens new perspectives for neutron scattering characterization of component functional molecular dynamics within the ribosome.

Identification of Specific Effect of Chloride on the Spectral Properties and Structural Stability of Multiple Extracellular Glutamic Acid Mutants of Bacteriorhodopsin

In the present work we combine spectroscopic, DSC and computational approaches to examine the multiple extracellular Glu mutants E204Q/E194Q, E204Q/E194Q/E9Q and E204Q/E194Q/E9Q/E74Q of bacteriorhodopsin by varying solvent ionic strength and composition. Absorption spectroscopy data reveal that the absorption maxima of multiple EC Glu mutants can be tuned by the chloride concentration in the solution. Visible Circular dichroism spectra imply that the specific binding of Cl- can modulate weakened exciton chromophore coupling and reestablish wild type-like bilobe spectral features of the mutants. The DSC data display reappearance of the reversible thermal transition, higher Tm of denaturation and an increase in the enthalpy of unfolding of the mutants in 1 M KCl solutions. Molecular dynamics simulations indicate high affinity binding of Cl- to Arg82 and to Gln204 and Gln194 residues in the mutants. Analysis of the experimental data suggests that simultaneous elimination of the negatively charged side chain of Glu194 and Glu204 is the major cause for mutants' alterations. Specific Cl- binding efficiently coordinates distorted hydrogen bonding interactions of the EC region and reconstitutes the conformation and structure stability of mutated bR in WT-like fashion.

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

BACKGROUND

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

SCOPE OF REVIEW

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

MAJOR CONCLUSIONS

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

GENERAL SIGNIFICANCE

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

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.

Dynamic footprint of sequestration in the molecular fluctuations of osteopontin

The sequestration of calcium phosphate by unfolded proteins is fundamental to the stabilization of biofluids supersaturated with respect to hydroxyapatite, such as milk, blood or urine. The unfolded state of osteopontin (OPN) is thought to be a prerequisite for this activity, which leads to the formation of core-shell calcium phosphate nanoclusters. We report on the structures and dynamics of a native OPN peptide from bovine milk, studied by neutron spectroscopy and small-angle X-ray and neutron scattering. The effects of sequestration are quantified on the nanosecond- ångström resolution by elastic incoherent neutron scattering. The molecular fluctuations of the free phosphopeptide are in agreement with a highly flexible protein. An increased resilience to diffusive motions of OPN is corroborated by molecular fluctuations similar to those observed for globular proteins, yet retaining conformational flexibilities. The results bring insight into the modulation of the activity of OPN and phosphopeptides with a role in the control of biomineralization. The quantification of such effects provides an important handle for the future design of new peptides based on the dynamics-activity relationship.

Specificity and affinity quantification of flexible recognition from underlying energy landscape topography

Flexibility in biomolecular recognition is essential and critical for many cellular activities. Flexible recognition often leads to moderate affinity but high specificity, in contradiction with the conventional wisdom that high affinity and high specificity are coupled. Furthermore, quantitative understanding of the role of flexibility in biomolecular recognition is still challenging. Here, we meet the challenge by quantifying the intrinsic biomolecular recognition energy landscapes with and without flexibility through the underlying density of states. We quantified the thermodynamic intrinsic specificity by the topography of the intrinsic binding energy landscape and the kinetic specificity by association rate. We found that the thermodynamic and kinetic specificity are strongly correlated. Furthermore, we found that flexibility decreases binding affinity on one hand, but increases binding specificity on the other hand, and the decreasing or increasing proportion of affinity and specificity are strongly correlated with the degree of flexibility. This shows more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Our work provides a theoretical foundation and quantitative explanation of the previous qualitative studies on the relationship among flexibility, affinity and specificity. In addition, we found that the folding energy landscapes are more funneled with binding, indicating that binding helps folding during the recognition. Finally, we demonstrated that the whole binding-folding energy landscapes can be integrated by the rigid binding and isolated folding energy landscapes under weak flexibility. Our results provide a novel way to quantify the affinity and specificity in flexible biomolecular recognition.

Flexibility in biomolecular recognition is crucial for the function. Flexibility often leads to moderate binding affinity but high binding specificity, challenging the conventional wisdom that high specificity is guaranteed by high affinity. Currently, understanding of the relationship between affinity and specificity in flexible biomolecular recognition is still obscure, even in a qualitative way. By exploring the intrinsic biomolecular recognition energy landscapes, we provided a novel way to quantify the thermodynamic intrinsic specificity by energy landscape topography and kinetic specificity by association rate. We show quantitatively that flexibility decreases binding affinity while increases binding specificity, and the relative changes in affinity and specificity are strongly correlated with the degree of flexibility. Our results show that more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Importantly, we demonstrated that flexibility modulates affinity and specificity through the underlying energy landscape. Our study establishes the quantitative relationship among flexibility, affinity and specificity, bridging the gap between theory and experiments.

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.

Zaccai neutron resilience and site-specific hydration dynamics in a globular protein

A discussion is presented of contributions of the Zaccai group to the understanding of flexibility in biological macromolecules using dynamic neutron scattering. The concept of resilience as introduced by Zaccai is discussed and investigated using molecular dynamics simulation on camphor-bound cytochrome P450. The resilience of hydrophilic residues is found to be more strongly affected by hydration than that of hydrophobic counterparts. The hydration-induced softening of protein propagates from the surface into the dry core. Moreover, buried hydrophilic residues behave more like those exposed on the protein surface, and are different from their hydrophobic counterparts.

Quantifying the topography of the intrinsic energy landscape of flexible biomolecular recognition

Biomolecular functions are determined by their interactions with other molecules. Biomolecular recognition is often flexible and associated with large conformational changes involving both binding and folding. However, the global and physical understanding for the process is still challenging. Here, we quantified the intrinsic energy landscapes of flexible biomolecular recognition in terms of binding-folding dynamics for 15 homodimers by exploring the underlying density of states, using a structure-based model both with and without considering energetic roughness. By quantifying three individual effective intrinsic energy landscapes (one for interfacial binding, two for monomeric folding), the association mechanisms for flexible recognition of 15 homodimers can be classified into two-state cooperative "coupled binding-folding" and three-state noncooperative "folding prior to binding" scenarios. We found that the association mechanism of flexible biomolecular recognition relies on the interplay between the underlying effective intrinsic binding and folding energy landscapes. By quantifying the whole global intrinsic binding-folding energy landscapes, we found strong correlations between the landscape topography measure Λ (dimensionless ratio of energy gap versus roughness modulated by the configurational entropy) and the ratio of the thermodynamic stable temperature versus trapping temperature, as well as between Λ and binding kinetics. Therefore, the global energy landscape topography determines the binding-folding thermodynamics and kinetics, crucial for the feasibility and efficiency of realizing biomolecular function. We also found "U-shape" temperature-dependent kinetic behavior and a dynamical cross-over temperature for dividing exponential and nonexponential kinetics for two-state homodimers. Our study provides a unique way to bridge the gap between theory and experiments.

Derivation of mean-square displacements for protein dynamics from elastic incoherent neutron scattering

The derivation of mean-square displacements from elastic incoherent neutron scattering (EINS) of proteins is examined, with the aid of experiments on camphor-bound cytochrome P450cam and complementary molecular dynamics simulations. It is shown that a q(4) correction to the elastic incoherent structure factor (where q is the scattering vector) can be simply used to reliably estimate from the experiment both the average mean-square atomic displacement, <Δr(2)> of the nonexchanged hydrogen atoms in the protein and its variance, σ(2). The molecular dynamics simulation results are in broad agreement with the experimentally derived <Δr(2)> and σ(2) derived from EINS on instruments at two different energy resolutions, corresponding to dynamics on the ∼100 ps and ∼1 ns time scales. Significant dynamical heterogeneity is found to arise from methyl-group rotations. The easy-to-apply q(4) correction extends the information extracted from elastic incoherent neutron scattering experiments and should be of wide applicability.

Dynamic properties of photosystem II membranes at physiological temperatures characterized by elastic incoherent neutron scattering. Increased flexibility associated with the inactivation of the oxygen evolving complex

Elastic incoherent neutron scattering (EINS), a non-invasive technique which is capable of measuring the mean square displacement of atoms in the sample, has been widely used in biology for exploring the dynamics of proteins and lipid membranes but studies on photosynthetic systems are scarce. In this study we investigated the dynamic characteristics of Photosystem II (PSII) membrane fragments between 280 and 340 K, i.e., in the physiological temperature range and in the range of thermal denaturation of some of the protein complexes. The mean square displacement values revealed the presence of a hydration-sensitive transition in the sample between 310 and 320 K, suggesting that the oxygen evolving complex (OEC) plays an important role in the transition. Indeed, in samples in which the OEC had been removed by TRIS- or heat-treatments (323 and 333 K) no such transition was found. Further support on the main role of OEC in these reorganizations is provided by data obtained from differential scanning calorimetry experiments, showing marked differences between the untreated and TRIS-treated samples. In contrast, circular dichroism spectra exhibited only minor changes in the excitonic interactions below 323 K, showing that the molecular organization of the pigment-protein complexes remains essentially unaffected. Our data, along with earlier incoherent neutron scattering data on PSII membranes at cryogenic temperatures (Pieper et al., Biochemistry 46:11398-11409, 2007), demonstrate that this technique can be applied to characterize the dynamic features of PSII membranes, and can be used to investigate photosynthetic membranes under physiologically relevant experimental conditions.

Sequence composition and environment effects on residue fluctuations in protein structures

Structure fluctuations in proteins affect a broad range of cell phenomena, including stability of proteins and their fragments, allosteric transitions, and energy transfer. This study presents a statistical-thermodynamic analysis of relationship between the sequence composition and the distribution of residue fluctuations in protein-protein complexes. A one-node-per-residue elastic network model accounting for the nonhomogeneous protein mass distribution and the interatomic interactions through the renormalized inter-residue potential is developed. Two factors, a protein mass distribution and a residue environment, were found to determine the scale of residue fluctuations. Surface residues undergo larger fluctuations than core residues in agreement with experimental observations. Ranking residues over the normalized scale of fluctuations yields a distinct classification of amino acids into three groups: (i) highly fluctuating-Gly, Ala, Ser, Pro, and Asp, (ii) moderately fluctuating-Thr, Asn, Gln, Lys, Glu, Arg, Val, and Cys, and (iii) weakly fluctuating-Ile, Leu, Met, Phe, Tyr, Trp, and His. The structural instability in proteins possibly relates to the high content of the highly fluctuating residues and a deficiency of the weakly fluctuating residues in irregular secondary structure elements (loops), chameleon sequences, and disordered proteins. Strong correlation between residue fluctuations and the sequence composition of protein loops supports this hypothesis. Comparing fluctuations of binding site residues (interface residues) with other surface residues shows that, on average, the interface is more rigid than the rest of the protein surface and Gly, Ala, Ser, Cys, Leu, and Trp have a propensity to form more stable docking patches on the interface. The findings have broad implications for understanding mechanisms of protein association and stability of protein structures.

The low-temperature inflection observed in neutron scattering measurements of proteins is due to methyl rotation: direct evidence using isotope labeling and molecular dynamics simulations

There is increasing interest in the contribution of methyl groups to the overall dynamics measured by neutron scattering experiments of proteins. In particular an inflection observed in atomic mean square displacements measured as a function of temperature on high resolution spectrometers (approximately 1 microeV) was explained by the onset of methyl group rotations. By specifically labeling a non-methyl-containing side-chain in a native protein system, the purple membrane, and performing neutron scattering measurements, we here provide direct experimental evidence that the observed inflection is indeed due to methyl group rotations. Molecular dynamics simulations reproduce the experimental data, and their analysis suggests that the apparent transition is due to methyl group rotation entering the finite instrumental resolution of the spectrometer. Methyl group correlation times measured by solid state NMR in the purple membrane, taken from previous work, support the interpretation.

Study of solvent-protein coupling effects by neutron scattering

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

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.

Consistent improvement of cross-docking results using binding site ensembles generated with elastic network normal modes

The representation of protein flexibility is still a challenge for the state-of-the-art flexible ligand docking protocols. In this article, we use a large and diverse benchmark to prove that is possible to improve consistently the cross-docking performance against a single receptor conformation, using an equilibrium ensemble of binding site conformers. The benchmark contained 28 proteins, and our method predicted the top-ranked near native ligand poses 20% more efficiently than using a single receptor. The multiple conformations were derived from the collective variable space defined by all heavy-atom elastic network normal modes, including backbone and side chains. We have found that the binding site displacements for best positioning of the ligand seem rather independent from the global collective motions of the protein. We also found that the number of binding site conformations needed to represent nonredundant flexibility was < 100. The ensemble of receptor conformations can be generated at our Web site at http://abagyan.scripps.edu/MRC.

Instantaneous normal modes and the protein glass transition

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

The bacteriorhodopsin carboxyl-terminus contributes to proton recruitment and protein stability

We examined functional and structural roles for the bacteriorhodopsin (bR) carboxyl-terminus. The extramembranous and intracellular carboxyl-terminus was deleted by insertion of premature translation stop codons. Deletion of the carboxyl-terminus had no effect on purple membrane (PM) lattice dimensions, sheet size, or the electrogenic environment of the ground-state chromophore. Removal of the distal half of the carboxyl-terminus had no effect on light-activated proton pumping, however, truncation of the entire carboxyl-terminus accelerated the rates of M-state decay and proton uptake approximately 3.7-fold and severely distorted the kinetics of proton uptake. Differential scanning calorimetry (DSC) and SDS denaturation demonstrated that removal of the carboxyl-terminus decreased protein stability. The DSC melting temperature was lowered by 6 degrees C and the calorimetric enthalpy reduced by 50% following removal of the carboxyl-terminus. Over the time range of milliseconds to hours at least 3 phases were required to describe the SDS denaturation kinetics for each bR construction. The fastest phases were indistinguishable for all bR's, and reflected PM solubilization. At pH 7.4, 20 degrees C, and in 0.3% SDS (w/v) the half-times of bR denaturation were 19.2 min for the wild-type, 12.0 min for the half-truncation and 3.6 min for the full-truncation. Taken together the results of this study suggest that the bR ground state exhibits two "domains" of stability: (1) a core chromophore binding pocket domain that is insensitive to carboxyl-terminal interactions and (2) the surrounding helical bundle whose contributions to protein stability and proton pumping are influenced by long-range interactions with the extramembranous carboxyl-terminus.

New sources and instrumentation for neutrons in biology

Neutron radiation offers significant advantages for the study of biological molecular structure and dynamics. A broad and significant effort towards instrumental and methodological development to facilitate biology experiments at neutron sources worldwide is reviewed.

Deuterium labeling for neutron structure-function-dynamics analysis

Neutron scattering and diffraction provide detailed information on the structure and dynamics of biological materials across time and length scales that range from picoseconds to nanoseconds and from 1 to 10,000 A, respectively. The particular sensitivity of neutrons to the isotopes of hydrogen makes selective deuterium labeling of biological systems an essential tool for maximizing the return from neutron scattering experiments. In neutron protein crystallography, the use of fully deuterated protein crystals improves the signal-to-noise ratio of the data by an order of magnitude and enhances the visibi-lity of the molecular structure (Proc Natl Acad Sci U S A 97:3872-3877, 2000; Acta Crystallogr D Biol Crystallogr 61:1413-1417, 2005; Acta Crystallogr D Biol Crystallogr 61:539-544, 2005). In solution and surface scattering experiments, the incorporation of deuterium-labeled subunits or components into complex assemblies or structures makes it possible to deconvolute the scattering of the labeled and unlabeled subunits and to determine their relative dispositions within the complex (J Mol Biol 93:255-265, 1975). With multiple labeling patterns, it is also possible to reconstruct the locations of multiple subunits in ternary and higher-order complexes (Science 238:1403-1406, 1987; J Mol Biol 271:588-601, 1997; J Biol Chem 275:14432-14439, 2000; Biochemistry 42:7790-7800, 2003). In inelastic neutron scattering experiments, which probe hydrogen dynamics in biological materials, the application of site, residue, or region-specific hydrogen-deuterium-labeling patterns can be used to distinguish and highlight the specific dynamics within a system (Proc Natl Acad Sci U S A 95:4970-4975, 1998).Partial, selective, or fully deuterated proteins can be readily produced by endogenous expression of recombinant proteins in bacterial systems that are adapted to growth in D(2)O solution and using selectively deuterated carbon sources. Adaptation can be achieved either by gradual step-wise increase in D(2)O concentration or, more directly, by plating cells on media of choice and selecting colonies that perform best for subsequent culture and inoculation. Scale-up growth and expression is typically performed in standard shaker flasks using either commercial or "home-grown" rich media (derived, for example, from cell lysates produced from algae grown in D(2)O) or under more controlled conditions in defined minimal media. Cell growth is typically slower in deuterated media (>5 times slower) and yields are correspondingly lower. Once the target protein has been expressed, purification proceeds by the protocols developed for the hydrogenated protein. The deuteration levels of the final product are determined by mass spectrometry.

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.

Dynamic fluctuation of proteins watched in real time

The dynamic nature of protein function is a fundamental concept in the physics of proteins. Although the basic general ideas are well accepted most experimental evidence has an indirect nature. The detailed characterization of the dynamics is necessary for the understanding in detail. The dynamic fluctuations thought crucial for the function span an extremely broad time, starting from the picosecond regime. Recently, a few new experimental techniques emerged that permit the observation of dynamical phenomena directly. Notably, pulsed infrared (IR) spectroscopy has been applied with great success to observe structural changes with picosecond time resolution. Using two-dimensional-IR vibrational echo chemical exchange spectroscopy Ishikawa and co-workers [Ishikawa et al. (2008), Proc. Natl. Acad. Sci. U.S.A. 101, 14402-14407] managed to observe the transition between well defined conformational substrates of carbonmonoxy myoglobin directly. This is an important step in improving our insight into the details of protein function.

Picometer-scale conformational heterogeneity separates functional from nonfunctional states of a photoreceptor protein

Protein structural fluctuations occur over a wide spatial scale, ranging from minute, picometer-scale displacements, to large, interdomain motions and partial unfolding. While large-scale protein structural changes and their effects on protein function have been the focus of much recent attention, small-scale fluctuations have been less well studied, and are generally assumed to have proportionally smaller effects. Here we use the bacterial photoreceptor photoactive yellow protein (PYP) to test if subtle structural changes do, indeed, imply equally subtle functional effects. We flash froze crystals of PYP to trap the protein's conformational ensemble, and probed the molecules in this ensemble for their ability to facilitate PYP's biological function (i.e., light-driven isomerization of its chromophore). Our results indicate that the apparently homogeneous structural state observed in a 0.82 A crystal structure in fact comprises an ensemble of conformational states, in which subpopulations with nearly identical structures display dramatically different functional properties.

Large scale motions in a biosensor protein glucose oxidase: a combined approach by QENS, normal mode analysis, and molecular dynamics studies

The characteristics of the glucose oxidase were studied using a combination of experimental and theoretical techniques. Quasi elastic neutron scattering experiments were used to obtain the vibrational frequencies of the protein. These were compared to theoretical results obtained by normal mode analysis. Results indicate a good match between the experimental and theoretical values. Molecular dynamic simulation with covariant analysis was used to study the structure and dynamics of glucose oxidase. Various parameters like the radius of gyration, root mean square fluctuations, solvent accessibility were studied for evaluating the structural stability of the protein. The frequency of vibration calculated from the three methods is used to derive the large scale motions. Theses studies were used to predict the suitable lysine residues for linkage with carbon nanotubes.

Role of extracellular glutamic acids in the stability and energy landscape of bacteriorhodopsin

Bacteriorhodopsin (BR), a specialized nanomachine, converts light energy into a proton gradient to power Halobacterium salinarum. In this work, we analyze the mechanical stability of a BR triple mutant in which three key extracellular residues, Glu(9), Glu(194), and Glu(204), were mutated simultaneously to Gln. These three Glu residues are involved in a network of hydrogen bonds, in cation binding, and form part of the proton release pathway of BR. Changes in these features and the robust photocycle dynamics of wild-type (WT) BR are apparent when the three extracellular Glu residues are mutated to Gln. It is speculated that such functional changes of proteins go hand in hand with changes in their mechanical properties. Here, we apply single-molecule dynamic force spectroscopy to investigate how the Glu to Gln mutations change interactions, reaction pathways, and the energy barriers of the structural regions of WT BR. The altered heights and positions of individual energy barriers unravel the changes in the mechanical and the unfolding kinetic properties of the secondary structures of WT BR. These changes in the mechanical unfolding energy landscape cause the proton pump to choose unfolding pathways differently. We suggest that, in a similar manner, the changed mechanical properties of mutated BR alter the functional energy landscape favoring different reaction pathways in the light-induced proton pumping mechanism.

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

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

Dynamical heterogeneity of specific amino acids in bacteriorhodopsin

Components of biological macromolecules, complexes and membranes are animated by motions occurring over a wide range of time and length scales, the synergy of which is at the basis of biological activity. Understanding biological function thus requires a detailed analysis of the underlying dynamical heterogeneity. Neutron scattering, using specific isotope labeling, and molecular dynamics simulations were combined in order to study the dynamics of specific amino acid types in bacteriorhodopsin within the purple membrane (PM) of Halobacterium salinarum. Motions of leucine, isoleucine and tyrosine residues on the pico- to nanosecond time scale were examined separately as a function of temperature from 20 to 300 K. The dynamics of the three residue types displayed different temperature dependence: isoleucine residues have larger displacements compared to the global PM above 120 K; leucine residues have displacements similar to that of PM in the entire temperature range studied; and tyrosine residues have displacements smaller than that of the average membrane in an intermediate temperature range. Experimental features were mostly well reproduced by molecular dynamics simulations performed at five temperatures, which allowed the dynamical characterisation of the amino acids under study as a function of local environment. The resulting dynamical map of bacteriorhodopsin revealed that movements of a specific residue are determined by both its environment and its residue type.

Hydration dependence of active core fluctuations in bacteriorhodopsin

We used neutron scattering and specific hydrogen-deuterium labeling to investigate the thermal dynamics of isotope-labeled amino acids and retinal, predominantly in the active core and extracellular moiety of bacteriorhodopsin (BR) in the purple membrane and the dynamical response to hydration. Measurements on two neutron spectrometers allowed two populations of motions to be characterized. The lower amplitude motions were found to be the same for both the labeled amino acids and retinal of BR and the global membrane. The larger amplitude dynamics of the labeled part, however, were found to be more resilient than the average membrane, suggesting their functional importance. The response to hydration was characterized, showing that the labeled part of BR is not shielded from hydration effects. The results suggest that the inhibition of high-amplitude motions by lowering hydration may play a key role in the slowing down of the photocycle and the proton pumping activity of BR.

Dynamics of hydration water in deuterated purple membranes explored by neutron scattering

The function and dynamics of proteins depend on their direct environment, and much evidence has pointed to a strong coupling between water and protein motions. Recently however, neutron scattering measurements on deuterated and natural-abundance purple membrane (PM), hydrated in H2O and D2O, respectively, revealed that membrane and water motions on the ns-ps time scale are not directly coupled below 260 K (Wood et al. in Proc Natl Acad Sci USA 104:18049-18054, 2007). In the initial study, samples with a high level of hydration were measured. Here, we have measured the dynamics of PM and water separately, at a low-hydration level corresponding to the first layer of hydration water only. As in the case of the higher hydration samples previously studied, the dynamics of PM and water display different temperature dependencies, with a transition in the hydration water at 200 K not triggering a transition in the membrane at the same temperature. Furthermore, neutron diffraction experiments were carried out to monitor the lamellar spacing of a flash-cooled deuterated PM stack hydrated in H2O as a function of temperature. At 200 K, a sudden decrease in lamellar spacing indicated the onset of long-range translational water diffusion in the second hydration layer as has already been observed on flash-cooled natural-abundance PM stacks hydrated in D2O (Weik et al. in J Mol Biol 275:632-634, 2005), excluding thus a notable isotope effect. Our results reinforce the notion that membrane-protein dynamics may be less strongly coupled to hydration water motions than the dynamics of soluble proteins.