Conformational fluctuation of native-like and molten-globule-like cytochrome c observed by time-resolved hole burning

Nanosecond to millisecond conformational fluctuations of Zn-substituted cytochrome c (ZnCytc) have been studied by the time-resolved transient hole-burning method. The investigation of low-temperature dynamics has been made on the ZnCytc solution sample in a water-glycerol mixture. The conformational fluctuations in the native-like and the molten-globule (MG)-like states have been compared for the aqueous solution samples at room temperature. ZnCytc in the MG-like state has been prepared by adding 200 mM NaClO4 to the protein solution with a pH of 2.1, and the formation of the MG-like state has been confirmed by both the far-UV CD and the visible absorption spectra. The hole spectrum of ZnCytc has been found to consist of two nearly degenerate components, that is, the Qx and Qy bands. The temporal change in the Qx component hole spectrum has been extracted by fitting the observed hole spectrum to the three-Gaussian form. The experimental results for ZnCytc dissolved in a water-glycerol mixture have revealed that the conformational fluctuation of ZnCytc is suppressed around 200 K, which is nearly the same temperature as the glass-like transition point of Zn-substituted myoglobin (ZnMb) and also as the glass-transition point of the solvent. This supports the idea of the solvent-induced glass-like transition of a protein. It has been also found that at physiological temperatures the time scale of the conformational fluctuation of ZnCytc lies around a few tens of nanoseconds, which is 2-3 orders of magnitude faster than that of ZnMb. The experimental results for the aqueous solution samples have shown that the difference between the native-like and the MG-like states is not conspicuous. However, they are indicative of the appearance of the slower conformational fluctuation in the MG-like state.     

Real-time observation of conformational fluctuations in Zn-substituted myoglobin by time-resolved transient hole-burning spectroscopy.

Equilibrium fluctuations of the protein conformation have been studied in myoglobin by a novel method of time-resolved transient hole-burning spectroscopy over a temperature range of 180-300 K and a time range of 10 ns to 10 ms. The temporal shift of the hole spectrum has been observed in a wide temperature region of 200-300 K. It has been found that the time behavior of the peak position of the hole is highly nonexponential and can be expressed by a stretched exponential function with a beta value of 0.22. As compared with the results for a dye solution sample, the time scale of the fluctuation of the protein conformation is much more weakly dependent on temperature. The time scale of the observed conformational dynamics shows a temperature dependence similar to that associated with the ligand escape process of myoglobin.     

Viscosity scaling and protein dynamics

The rates of molecular motions in the interior of some proteins were found to scale with an inverse power of the external solvent viscosity. The data were explained by a flexible protein structure whose dynamics is partially controlled by the solvent. Reaction dynamics in the presence of structural fluctuations with finite lifetimes lead to a dynamic friction coefficient defined by a generalized Langevin equation and a fluctuation-dissipation theorem. A model for the dynamic friction is derived assuming that the fluctuation spectrum at the reaction site involves two components: solvent-independent diffusion of local structural defects in the protein matrix and global fluctuations coupled to the solvent. The theory is applied to the viscosity dependence of molecular oxygen-binding rates in sperm whale myoglobin.     

HPLC study on the 'history' dependence of gramicidin A conformation in phospholipid model membranes

A novel HPLC methodology for the study of gramicidin A reconstituted in model membranes has been tested in comparison with circular dichroism data. It is shown that this chromatographic technique not only corroborates most of the recent spectroscopic results but allows one to explain them in terms of mass fractions of different actual conformational species of GA in the phospholipid assemblies. In particular, the dependence of the inserted peptide configuration on the organic solvent and other parameters involved in the 'history' of the sample preparation and handling has been analyzed by HPLC in two phospholipid model systems: small unilamellar vesicles and micelles. Moreover, a slow conformational transition of GA towards a beta 6.3-helical configuration, accelerated by heat incubation, has been also chromatographically visualized and quantitatively interpreted.     

Dynamic aspect of kinetics of reaction catalyzed by lactate dehydrogenase

The influence of solvent viscosity on the kinetic parameters of the pyruvate reduction reaction catalyzed by lactate dehydrogenase has been investigated. The viscosity was adjusted by sucrose and glycerol solutions at concentrations from 0 to 44% and from 0 to 63%, respectively. The reaction rate decreased abruptly with an increase in viscosity. The study of different reaction stages (enzyme-substrate complex formation, catalysis, inhibitory complex decomposition, competitive inhibition by chlorine ions) revealed that the catalysis (and the related conformational changes) is the only stage (of the above mentioned) that depends markedly on the solvent viscosity. The reaction is insensitive to the changes in the dielectric properties of the solution induced by the addition of alcohols and dioxane. The observed power dependence of the rate constant on viscosity is explained in terms of Kramer's theory which considers the proton transition through the activation barrier to be a diffusion in the field of random forces. The influence of solvent viscosity on enzymic kinetics indicates a direct relation between solvent dynamics and relevant protein conformational movements.     

A diffusion Michaelis-Menten mechanism: continuous conformational change in enzymatic kinetics

We present a simple model which extends the Michaelis-Menten mechanism by incorporating a continuous protein conformational change in enzymatic catalysis. This model can represent a quantitative version for "rack" or "induced fit" mechanisms. In the steady-state it leads to an equation of the Michaelis-Menten form, but with the catalytic step at the active site showing strong dependence on solvent viscosity. We suggest that a careful examination of solvent viscosity effects on enzymatic activity may serve as a test for the conformational change hypothesis.     

The membrane as an environment of minimal interconversion. A circular dichroism study on the solvent dependence of the conformational behavior of gramicidin in diacylphosphatidylcholine model membranes

The conformation of gramicidin in diacylphosphatidylcholine model membranes was investigated as a function of the solvent in which peptide and lipid are initially codissolved. By use of circular dichroism it is demonstrated that, upon removal of the solvent and hydration of the mixed gramicidin/lipid film, it is the conformational behavior of the peptide in the organic solvent that determines its final conformation in dimyristoylphosphatidylcholine model membranes. As a consequence, parameters that influence the conformation of the peptide in the solvent also play an essential role, such as the gramicidin concentration and the rate of interconversion between different conformations. Of the various solvents investigated, only with trifluoroethanol is it possible directly to incorporate gramicidin entirely in the beta 6.3-helical (channel) configuration. It is also shown that the conformation of gramicidin in the membrane varies with the peptide/lipid ratio, most likely as a result of intermolecular gramicidin-gramicidin interactions at higher peptide/lipid ratios, and that heat incubation leads to a conformational change in the direction of the beta 6.3-helical conformation. Using lipids with an acyl chain length varying from 12 carbon atoms in dilauroylphosphatidylcholine to 22 carbon atoms in dierucoylphosphatidylcholine, it was possible to investigate the acyl chain length dependence of the gramicidin conformation in model membranes prepared from these lipids with the use of different solvent systems. It is demonstrated for each solvent system that the distribution between different conformations is relatively independent of the acyl chain length but that the rate at which the conformation converts toward the beta 6.3-helical configuration upon heating of the samples is affected by the length of the acyl chain.(ABSTRACT TRUNCATED AT 250 WORDS)     

Intermolecular interactions of oxygenated sickle hemoglobin molecules in cells and cell-free solutions.

We have measured the intermolecular interactions of oxygenated sickle hemoglobin molecules in cells and in cell-free solutions, and have compared the results with similar data for liganded normal adult hemoglobin. The experiments involve the measurement of the spin-lattice relaxation time T1 of protons of solvent water molecules, as a function of an externally applied static magnetic field. From such data, one can derive a correlation time tauc, for each sample, which is a measure of the time taken for a hemoglobin molecule to randomize its orientation due to Brownian motion. Thus tauc is a measure of the freedom of rotational motion, on a molecular or microscopic level, of hemoglobin molecules. Intermolecular interactions will reduce this freedom of motion and lengthen tauc. We find that oxygenated sickle hemoglobin molecules have an additional intermolecular interaction not found for normal hemoglobin. This extra interaction is increased by the presence of either inorganic phosphate or diphosphoglycerate, and is greater for sickle hemoglobin within cells than in cell-free solutions. By comparing the present results with published data on the viscosity of oxygenated sickle and normal hemoglobin, we conclude that, at concentrations comparable to intracellular values, oxygenated sickle hemoglobin molecules form aggregates several tetramers in size. The possibility exists that these aggregates are the earliest stage of fiber formation itself, the physical basis of the sickling phenomena.     

Long time dynamics of Met-enkephalin: comparison of explicit and implicit solvent models

Met-enkephalin is one of the smallest opiate peptides. Yet, its dynamical structure and receptor docking mechanism are still not well understood. The conformational dynamics of this neuron peptide in liquid water are studied here by using all-atom molecular dynamics (MD) and implicit water Langevin dynamics (LD) simulations with AMBER potential functions and the three-site transferable intermolecular potential (TIP3P) model for water. To achieve the same simulation length in physical time, the full MD simulations require 200 times as much CPU time as the implicit water LD simulations. The solvent hydrophobicity and dielectric behavior are treated in the implicit solvent LD simulations by using a macroscopic solvation potential, a single dielectric constant, and atomic friction coefficients computed using the accessible surface area method with the TIP3P model water viscosity as determined here from MD simulations for pure TIP3P water. Both the local and the global dynamics obtained from the implicit solvent LD simulations agree very well with those from the explicit solvent MD simulations. The simulations provide insights into the conformational restrictions that are associated with the bioactivity of the opiate peptide dermorphin for the delta-receptor.     

Viscosity dependence of protein dynamics

The influence of solvent viscosity on protein dynamics was investigated with molecular dynamics simulations of factor Xa in two solvents differing only in viscosity, by a factor of 10. We obtained this viscosity change by changing the masses of the solvent atoms by a factor of 100. Equilibrium properties of the protein, that is, the average structure, its fluctuations, and the secondary structure, show no significant dependence on the solvent viscosity. The dynamic properties of the protein, that is, the atom-positional correlation times and torsional angle transitions, however, depend on the solvent viscosity. The protein appears to be much more mobile in the solvent of lower viscosity. It feels the influence of the solvent not only on the surface but even in its core. With increasing solvent viscosity, the positional relaxation times of atoms in the protein core increase as much as those of atoms on the protein surface, and the relative increase in the core is even larger than on the surface.     

Quaternary relaxations in sol-gel encapsulated hemoglobin studied via NIR and UV spectroscopy

In this work, we study the kinetics of the R --> T transition in hemoglobin using a combination of near-infrared and near-ultraviolet spectroscopy. We use a sol-gel encapsulation protocol to decelerate the conformational transitions and to avoid spectral perturbations arising from ligand migration and recombination. We monitor two spectroscopic markers: band III in the near-IR, which is a fine probe of the heme pocket conformation, and the tryptophan band in the near-UV, which probes the formation of the Trpbeta37-Aspalpha94 hydrogen bond, characteristic of the T structure, at the critical alpha1beta2 subunit interface. The time evolution of these two bands is monitored after deoxygenation of encapsulated oxyhemoglobin, obtained by diffusion of a reducing agent into the porous silica matrix. Characteristic spectral shifts are observed: comparison with myoglobin enables us to assign them to quaternary structure relaxations. Band III spectral relaxation is clearly nonexponential, and analysis with the Maximum Entropy Method enables us to identify three processes. On the other hand, near-UV spectral relaxation follows an exponential decay with a time constant closely corresponding to the second process observed in the near IR. Very interestingly, the rates of all processes markedly depend on the viscosity of the co-encapsulated solvent, following a power law. Our results reveal correlations between heme pocket relaxations, induced by the R --> T transition, and structural event(s) occurring at the alpha1beta2 interface and highlight their solvent dependence. The power law viscosity dependence of relaxation rates suggests that the observed protein relaxations are "slaved" to the co-encapsulated solvent. The stepwise character of the quaternary transition is also evidenced.     

Bending elasticity of anti-parallel beta-sheets

Using a coarse-grained elastic model, we examine the bending properties of anti-parallel beta-sheets comprised of uniform amino-acid residues in vacuum as well as in explicit solvent. By comparing the conformational probability of the beta-sheet from molecular dynamics simulations with the same quantities obtained from the coarse-grained model, we compute the elastic bending constant, kappa. Equilibrium fluctuations of the beta-sheet and its response to external forces are well reproduced by a model with a uniform isotropic bending constant. An anisotropic bending model is also investigated, although the computed anisotropy is relatively weak and most of the observed properties are well described by an isotropic model. The presence of explicit solvent also lowers the bending constant. The sequence dependence of our result and its implications in protein conformational dynamics are discussed.     

Hydration dependence of backbone and side chain polylysine dynamics: a 13C solid-state NMR and IR spectroscopy study

The molecular dynamics of solid poly-L-lysine has been studied by the following natural abundance (13)C-NMR relaxation methods: measurements of the relaxation times T(1) at two resonance frequencies, off-resonance T(1rho) at two spin-lock frequencies, and proton-decoupled T(1rho). Experiments were performed at different temperatures and hydration levels (up to 17% H(2)O by weight). The natural abundance (13)C-CPMAS spectrum of polylysine provides spectral resolution of all types of backbone and side chain carbons and thus, dynamic parameters could be determined separately for each of them. At the same time, the conformational properties of polylysine were investigated by Fourier transform infrared spectroscopy. The data obtained from the different NMR experiments were simultaneously analyzed using the correlation function formalism and model-free approach. The results indicate that in dry polylysine both backbone and side chains take part in two low amplitude motions with correlation times of the order of 10(-4) s and 10(-9) s. Upon hydration, the dynamic parameters of the backbone remain almost constant except for the amplitude of the slower process that increases moderately. The side chain dynamics reveals a much stronger hydration response: the amplitudes of both slow and fast motions increase significantly and the correlation time of the slow motion shortens by about five orders of magnitude, and at hydration levels of more than 10% H(2)O fast and slow side chain motions are experimentally indistinguishable. These changes in the molecular dynamics cannot be ascribed to any hydration-dependent conformational transitions of polylysine because IR spectra reveal almost no hydration dependence in either backbone or side chain absorption domains. The physical nature of the fast and slow motions, their correlation time distributions, and hydration dependence of microdynamic parameters are discussed.     

Conformational states and fluctuations in endothelial nitric oxide synthase under calmodulin regulation

Mechanisms that regulate nitric oxide synthase enzymes (NOS) are of interest in biology and medicine. Although NOS catalysis relies on domain motions and is activated by calmodulin (CaM) binding, the relationships are unclear. We used single-molecule fluorescence resonance energy transfer (FRET) spectroscopy to elucidate the conformational states distribution and associated conformational fluctuation dynamics of the two NOS electron transfer domains in an FRET dye-labeled endothelial NOS reductase domain (eNOSr) and to understand how CaM affects the dynamics to regulate catalysis by shaping the spatial and temporal conformational behaviors of eNOSr. In addition, we developed and applied a new imaging approach capable of recording three-dimensional FRET efficiency versus time images to characterize the impact on dynamic conformal states of the eNOSr enzyme by the binding of CaM, which identifies clearly that CaM binding generates an extra new open state of eNOSr, resolving more detailed NOS conformational states and their fluctuation dynamics. We identified a new output state that has an extra open conformation that is only populated in the CaM-bound eNOSr. This may reveal the critical role of CaM in triggering NOS activity as it gives conformational flexibility for eNOSr to assume the electron transfer output FMN-heme state. Our results provide a dynamic link to recently reported EM static structure analyses and demonstrate a capable approach in probing and simultaneously analyzing all of the conformational states, their fluctuations, and the fluctuation dynamics for understanding the mechanism of NOS electron transfer, involving electron transfer among FAD, FMN, and heme domains, during nitric oxide synthesis.     

Dynamical transition of myoglobin in a crystal: comparative studies of X-ray crystallography and Mössbauer spectroscopy

The crystallographic normal mode refinements of myoglobin at a wide range of temperature from 40 K to 300 K were carried out to study the temperature dependence of the internal atomic fluctuations. The refinement method decomposes the mean square displacement from the average position, (deltar2), into the contributions from the internal degrees of freedom and those from the external degrees of freedom. The internal displacements show linear temperature dependence as (deltar2)=alphaT+beta, throughout the temperature range measured here, and exhibit no obvious change in the slope alpha at the dynamical transition temperature (Tc=ca. 180 K). The slope alpha is practically the same as the value predicted theoretically by normal mode analysis. Such linear dependence is considered to be due to the following reason. The crystallographic Debye-Waller factor represents the static distribution caused by convolution of temperature-dependent normal mode motions and a temperature-independent set of the conformational substates. In contrast, M?ssbauer absorption spectroscopy shows a clear increase in the gradient alpha at Tc. This difference from X-ray diffraction originates from the incoherent nature of the M?ssbauer effect together with its high-energy resolution, which yields the self-correlation, and the temporal behavior of individual Fe atoms in the myoglobin crystal.     

The influence of a protein on water dynamics in its vicinity investigated by molecular dynamics simulation

A system containing the globular protein ubiquitin and 4,197 water molecules has been used for the analysis of the influence exerted by a protein on solvent dynamics in its vicinity. Using Voronoi polyhedra, the solvent has been divided into three subsets, i.e., the first and second hydration shell, and the remaining bulk, which is hardly affected by the protein. Translational motion in the first shell is retarded by a factor of 3 in comparison to bulk. Several molecules in the first shell do not reach the diffusive regime within 100 ps. Shell-averaged orientational autocorrelation functions, which are also subject to a retardation effect, cannot be modeled by a single exponential time law, but are instead well-described by a Kohlrausch-Williams-Watts (KWW) function. The underlying distribution of single-molecule rotational correlation times is both obtained directly from the simulation and derived theoretically. The temperature dependence of reorientation is characterized by a strongly varying correlation time, but a virtually temperature-independent KWW exponent. Thus, the coupling of water structure relaxation with the respective environment, which is characteristic of each solvation shell, is hardly affected by temperature. In other words, the functional form of the distributions of single-molecule rotational correlation times is not subject to a temperature effect. On average, a correlation between reorientation and lifetimes of neighborhood relations is observed. © 1996 Wiley-Liss, Inc.     

Viscous cosolvent effect on the ultrasonic absorption of bovine serum albumin.

Protein-ligand binding and enzyme activity have been shown to be regulated by solvent viscosity, induced by the addition of viscous cosolvents. This was indirectly interpreted as an effect on protein dynamics. However, viscous cosolvents might affect dynamic, e.g., viscosity, as well as thermodynamic properties of the solution, e.g., activity of solution components. This work was undertaken to examine the effect of viscous cosolvent on the structural dynamics of proteins and its correlation with dynamic and thermodynamic solution properties. For this purpose we studied the effect of viscous cosolvent on the specific ultrasonic absorption, delta mu, of bovine serum albumin, at pH = 7.0 and at 21 degrees C, and frequency range of 3-4 MHz. Ultrasonic absorption (UA) directly probes protein dynamics related to energy dissipation processes. It was found that the addition of sucrose, glycerol, or ethylene glycol increased the BSA delta mu. This increase correlates well with the solvent viscosity, but not with the cosolvent mass concentration, activity of the solvent components, dielectric constant, or the hydration of charged groups. On the grounds of these results and previously reported findings, as well as theoretical considerations, we propose the following mechanism for the solvent viscosity effect on the protein structural fluctuations, reflected in the UA: increased solvent viscosity alters the frequency spectrum of the polypeptide chain movements; attenuating the fast (small amplitude) movements, and enhancing the slow (large amplitude) ones. This modulates the interaction strength between the polypeptide and water species that "lubricates" the chain's movements, leading to larger protein-volume fluctuation and higher ultrasonic absorption. This study demonstrates that solvent viscosity is a regulator of protein structural fluctuations.     

Coupling of protein and environment fluctuations

We review the concepts of protein dynamics developed over the last 35years and extend applications of the unified model of protein dynamics to heat flow and spatial fluctuations in hydrated myoglobin (Mb) powders. Differential scanning calorimetry (DSC) and incoherent neutron scattering (INS) data on hydration Mb powders are explained by the temperature-dependence of the hydration-shell β(h) process measured by dielectric relaxation spectroscopy (DRS). The unified model explains the temperature dependence of DSC and INS data as a kinetic effect due to a fixed experimental time window and a broad distribution of hydration-shell β(h) fluctuation rates. We review the slaving of large scale protein motions to the bulk solvent α process, and the metastability of Mb molecules in glass forming bulk solvent at low temperatures. This article is part of a Special Issue entitled: "Protein Dynamics: Experimental and Computational Approaches".     

Diffusional barrier in the unfolding of a small protein

To determine how the dynamics of the polypeptide chain in a protein molecule are coupled to the bulk solvent viscosity, the unfolding by urea of the small protein barstar was studied in the presence of two viscogens, xylose and glycerol. Thermodynamic studies of unfolding show that both viscogens stabilize barstar by a preferential hydration mechanism, and that viscogen and urea act independently on protein stability. Kinetic studies of unfolding show that while the rate-limiting conformational change during unfolding is dependent on the bulk solvent viscosity, eta, its rate does not show an inverse dependence on eta, as expected by Kramers' theory. Instead, the rate is found to be inversely proportional to an effective viscosity, eta + xi, where xi is an adjustable parameter which needs to be included in the rate equation. xi is found to have a value of -0.7 cP in xylose and -0.5 cP in glycerol, in the case of unfolding, at constant urea concentration as well as under isostability conditions. Hence, the unfolding protein chain does not experience the bulk solvent viscosity, but instead an effective solvent viscosity, which is lower than the bulk solvent viscosity by either 0.7 cP or 0.5 cP. A second important result is the validation of the isostability assumption, commonly used in protein folding studies but hitherto untested, according to which if a certain concentration of urea can nullify the effect of a certain concentration of viscogen on stability, then the same concentrations of urea and viscogen will also not perturb the free energy of activation of the unfolding of the protein.     

Internal friction controls the speed of protein folding from a compact configuration

Several studies have found millisecond protein folding reactions to be controlled by the viscosity of the solvent: Reducing the viscosity allows folding to accelerate. In the limit of very low solvent viscosity, however, one expects a different behavior. Internal interactions, occurring within the solvent-excluded interior of a compact molecule, should impose a solvent-independent upper limit to folding speed once the bulk diffusional motions become sufficiently rapid. Why has this not been observed? We have studied the effect of solvent viscosity on the folding of cytochrome c from a highly compact, late-stage intermediate configuration. Although the folding rate accelerates as the viscosity declines, it tends toward a finite limiting value approximately 10(5) s(-1) as the viscosity tends toward zero. This limiting rate is independent of the cosolutes used to adjust solvent friction. Therefore, interactions within the interior of a compact denatured polypeptide can limit the folding rate, but the limiting time scale is very fast. It is only observable when the solvent-controlled stages of folding are exceedingly rapid or else absent. Interestingly, we find a very strong temperature dependence in these "internal friction"-controlled dynamics, indicating a large energy scale for the interactions that govern reconfiguration within compact, near-native states of a protein.