Thermal motions in bacteriorhodopsin at different hydration levels studied by neutron scattering: correlation with kinetics and light-induced conformational changes

Temperature dependence of protein dynamics: computer simulation analysis of neutron scattering properties

The temperature dependence of the internal dynamics of an isolated protein, bovine pancreatic trypsin inhibitor, is examined using normal mode analysis and molecular dynamics (MD) simulation. It is found that the protein exhibits marked anharmonic dynamics at temperatures of approximately 100-120 K, as evidenced by departure of the MD-derived average mean square displacement from that of the harmonic model. This activation of anharmonic dynamics is at lower temperatures than previously detected in proteins and is found in the absence of solvent molecules. The simulation data are also used to investigate neutron scattering properties. The effects are determined of instrumental energy resolution and of approximations commonly used to extract mean square displacement data from elastic scattering experiments. Both the presence of a distribution of mean square displacements in the protein and the use of the Gaussian approximation to the dynamic structure factor lead to quantified underestimation of the mean square displacement obtained.

Protein flexibility from the dynamical transition: a force constant analysis

A standard analysis of the scattered neutron incoherent elastic intensity measured with very good energy resolution yields elastic scans, i.e., mean-square displacements of atomic motions (in a pico to nanosecond time scale) in a sample as a function of temperature. This provides a quick way for characterizing the dynamical behavior of biological macromolecules, such behavior being correlated with biological function and activity. Elastic scans of proteins exhibit a dynamical transition at approximately 200 K, marking a cross-over in molecular fluctuations between harmonic and nonharmonic dynamical regimes. This paper presents an approach allowing analysis of the elastic scan in terms of force constants and related parameters, such as the free energy barrier DeltaG at the transition. We find that the increased protein flexibility beyond the dynamical transition is associated with DeltaG approximately equals RT and effective force constants of the order of 0.1-3 N/m. The analysis provides a set of parameters for characterizing molecular resilience and exploring relations among dynamics, function, and activity in proteins.

Large-scale shape changes in proteins and macromolecular complexes

Proteins and RNA undergo intricate motions as they carry out functions in biological systems. These motions frequently entail large-scale conformational changes that induce changes in the surface structure, or shape, of a molecule. This review describes the experimental characterization of large-scale shape changes in proteins and macromolecular complexes and the effects of such changes on macromolecular behavior. We describe several important results that have been obtained by using small-angle scattering, which is emerging as a powerful technique for determining macromolecular shapes and elucidating the quaternary structure of macromolecular assemblies.

The dynamics of protein hydration water: a quantitative comparison of molecular dynamics simulations and neutron-scattering experiments

We present results from an extensive molecular dynamics simulation study of water hydrating the protein Ribonuclease A, at a series of temperatures in cluster, crystal, and powder environments. The dynamics of protein hydration water appear to be very similar in crystal and powder environments at moderate to high hydration levels. Thus, we contend that experiments performed on powder samples are appropriate for discussing hydration water dynamics in native protein environments. Our analysis reveals that simulations performed on cluster models consisting of proteins surrounded by a finite water shell with free boundaries are not appropriate for the study of the solvent dynamics. Detailed comparison to available x-ray diffraction and inelastic neutron-scattering data shows that current generation force fields are capable of accurately reproducing the structural and dynamical observables. On the time scale of tens of picoseconds, at room temperature and high hydration, significant water translational diffusion and rotational motion occur. At low hydration, the water molecules are translationally confined but display appreciable rotational motion. Below the protein dynamical transition temperature, both translational and rotational motions of the water molecules are essentially arrested. Taken together, these results suggest that water translational motion is necessary for the structural relaxation that permits anharmonic and diffusive motions in proteins. Furthermore, it appears that the exchange of protein-water hydrogen bonds by water rotational/librational motion is not sufficient to permit protein structural relaxation. Rather, the complete exchange of protein-bound water molecules by translational displacement seems to be required.

Moist and soft, dry and stiff: a review of neutron experiments on hydration-dynamics-activity relations in the purple membrane of Halobacterium salinarum

Twenty-five years of neutron experiments on hydration and thermal dynamics in purple membranes of Halobacterium salinarum are reviewed. Neutron diffraction, elastic and quasielastic scattering, allowed to map the distribution of water and lipids and to measure thermal fluctuations and correlation times in the membranes, under various conditions of temperature, hydration and lipid environment. Strong correlations were established between dynamics parameters and the activity of bacteriorhodopsin (the purple membrane protein), as a light driven proton pump supporting the hypothesis that the influence of hydration on activity is in fact due to its effects on membrane thermal dynamics. Hydrogen-deuterium labelling experiments highlighted stiffer and softer parts in the bacteriorhodopsin structure. The soft parts would allow the conformational changes involved in activity, while the stiffer ones may control a valve-like function in vectorial proton transfer.

How soft is a protein? A protein dynamics force constant measured by neutron scattering

An effective environmental force constant is introduced to quantify the molecular resilience (or its opposite, "softness") of a protein structure and relate it to biological function and activity. Specific resilience-function relations were found in neutron-scattering experiments on purple membranes containing bacteriorhodopsin, the light-activated proton pump of halobacteria; the connection between resilience and stability is illustrated by a study of myoglobin in different environments. Important advantages of the neutron method are that it can characterize the dynamics of any type of biological sample-which need not be crystalline or monodisperse-and that it enables researchers to focus on the dynamics of specific parts of a complex structure with deuterium labeling.

Enzyme dynamics and activity: time-scale dependence of dynamical transitions in glutamate dehydrogenase solution

We have examined the temperature dependence of motions in a cryosolution of the enzyme glutamate dehydrogenase (GDH) and compared these with activity. Dynamic neutron scattering was performed with two instruments of different energy resolution, permitting the separate determination of the average dynamical mean square displacements on the sub-approximately 100 ps and sub-approximately 5 ns time scales. The results demonstrate a marked dependence on the time scale of the temperature profile of the mean square displacement. The lowest temperature at which anharmonic motion is observed is heavily dependent on the time window of the instrument used to observe the dynamics. Several dynamical transitions (inflexions of the mean squared displacement) are observed in the slower dynamics. Comparison with the temperature profile of the activity of the enzyme in the same solvent reveals dynamical transitions that have no effect on GDH function.

Protein, lipid and water organization in bacteriorhodopsin crystals: a molecular view of the purple membrane at 1.9 A resolution

BACKGROUND

Bacteriorhodopsin (bR) from Halobacterium salinarum is a proton pump that converts the energy of light into a proton gradient that drives ATP synthesis. The protein comprises seven transmembrane helices and in vivo is organized into purple patches, in which bR and lipids form a crystalline two-dimensional array. Upon absorption of a photon, retinal, which is covalently bound to Lys216 via a Schiff base, is isomerized to a 13-cis,15-anti configuration. This initiates a sequence of events - the photocycle - during which a proton is transferred from the Schiff base to Asp85, followed by proton release into the extracellular medium and reprotonation from the cytoplasmic side.

RESULTS

The structure of bR in the ground state was solved to 1.9 A resolution from non-twinned crystals grown in a lipidic cubic phase. The structure reveals eight well-ordered water molecules in the extracellular half of the putative proton translocation pathway. The water molecules form a continuous hydrogen-bond network from the Schiff-base nitrogen (Lys216) to Glu194 and Glu204 and includes residues Asp85, Asp212 and Arg82. This network is involved both in proton translocation occurring during the photocycle, as well as in stabilizing the structure of the ground state. Nine lipid phytanyl moieties could be modeled into the electron-density maps. Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) analysis of single crystals demonstrated the presence of four different charged lipid species.

CONCLUSIONS

The structure of protein, lipid and water molecules in the crystals represents the functional entity of bR in the purple membrane of the bacteria at atomic resolution. Proton translocation from the Schiff base to the extracellular medium is mediated by a hydrogen-bond network that involves charged residues and water molecules.

Simulation analysis of the retinal conformational equilibrium in dark-adapted bacteriorhodopsin

In dark-adapted bacteriorhodopsin (bR) the retinal moiety populates two conformers: all-trans and (13,15)cis. Here we examine factors influencing the thermodynamic equilibrium and conformational transition between the two forms, using molecular mechanics and dynamics calculations. Adiabatic potential energy mapping indicates that whereas the twofold intrinsic torsional potentials of the C13==C14 and C15==N16 double bonds favor a sequential torsional pathway, the protein environment favors a concerted, bicycle-pedal mechanism. Which of these two pathways will actually occur in bR depends on the as yet unknown relative weight of the intrinsic and environmental effects. The free energy difference between the conformers was computed for wild-type and modified bR, using molecular dynamics simulation. In the wild-type protein the free energy of the (13,15)cis retinal form is calculated to be 1.1 kcal/mol lower than the all-trans retinal form, a value within approximately kBT of experiment. In contrast, in isolated retinal the free energy of the all-trans state is calculated to be 2.1 kcal/mol lower than (13,15)cis. The free energy differences are similar to the adiabatic potential energy differences in the various systems examined, consistent with an essentially enthalpic origin. The stabilization of the (13,15)cis form in bR relative to the isolated retinal molecule is found to originate from improved protein-protein interactions. Removing internal water molecules near the Schiff base strongly stabilizes the (13,15)cis form, whereas a double mutation that removes negative charges in the retinal pocket (Asp85 to Ala; Asp212 to Ala) has the opposite effect.