Inhomogeneous broadening in spectral bands of carbonmonoxymyoglobin. The connection between spectral and functional heterogeneity

Imaging protein statistical substate occupancy in a spectrum-function phase space

Hemeprotein ligand rebinding studies reveal varying IR absorbance and rebinding functions across a cryogenic ensemble. Since IR-active vibrations and rebinding barriers couple to structural coordinates, spectral and functional heterogeneity arise from conformational heterogeneity. Modeling rebinding data as a spectrally resolved superposition of first-order rate processes and employing maximum entropy regularization, protein heterogeneity is imaged as an ensemble occupancy of a spectrum-function phase space. Results from myoglobin rebinding carbon monoxide are discussed.

Kinetic model of primary energy transfer and trapping in photosynthetic membranes

The picosecond time-domain incoherent singlet excitation transfer and trapping kinetics in core antenna of photosynthetic bacteria are studied in case of low excitation intensities by numerical integration of the appropriate master equation in a wide temperature range of 4-300 K. The essential features of our two-dimensional-lattice model are as follows: Förster excitation transfer theory, spectral heterogeneity of both the light-harvesting antenna and the reaction center, treatment of temperature effects through temperature dependence of spectral bands, inclusion of inner structure of the trap, and transition dipole moment orientation. The fluorescence kinetics is analyzed in terms of distributions of various kinetic components, and the influence of different inhomogeneities (orientational, spectral) is studied.A reasonably good agreement between theoretical and experimental fluorescence decay kinetics for purple photosynthetic bacterium Rhodospirillum rubrum is achieved at high temperatures by assuming relatively large antenna spectral inhomogeneity: 20 nm at the whole bandwidth of 40 nm. The mean residence time in the antenna lattice site (it is assumed to be the aggregate of four bacteriochlorophyll a molecules bound to proteins) is estimated to be approximately 12 ps. At 4 K only qualitative agreement between model and experiment is gained. The failure of quantitative fitting is perhaps due to the lack of knowledge about the real structure of antenna or local heating and cooling effects not taken into account.

Conformational substates of ferricytochrome c revealed by combined optical absorption and electronic circular dichroism spectroscopy at cryogenic temperature

We have investigated the heterogeneity of the Fe(III)-Met80 linkage of horse heart ferricytochrome c by probing the 695nm charge transfer band with absorption and electronic circular dichroism (ECD) spectroscopy. In order to verify the connection between conformational substates of the Fe(III)-Met80 linkage and the 695nm band spectral heterogeneity, we have performed experiments as a function of pH (neutral and acidic) and temperature (room and 20K). At room temperature, the ECD spectrum is blue shifted with respect to the absorption one; the shift is more pronounced at acidic pH and is compatible with the presence of sub-bands. ECD measurements at 20K highlighted the heterogeneous nature of the 695nm band and provided direct experimental evidence for the presence of sub-bands. Indeed, while the absorption spectra remained deceivingly unstructured, the ECD spectra showed well resolved peaks and shoulders. A consistent fit of the 20K absorption and ECD spectra showed that five Gaussians (each centered at the same frequency in the absorption and ECD spectrum) are able to reproduce the observed lineshapes. A careful analysis of frequency shifts and intensity ratios of these sub-bands enabled us to identify at least three distinct sub-bands arising from taxonomic conformational substates of the Fe(III)-Met80 linkage. In view of the major influence of the Fe(III)-Met80 linkage on the redox potential of ferricytochrome c, we speculate that these spectrally distinguishable substates may have different functional roles.

Charge density-dependent modifications of hydration shell waters by Hofmeister ions

Gadolinium (Gd(3+)) vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe, as a function of added Hofmeister series salts, changes in the OH stretching frequency derived from first-shell waters of aqueous Gd(3+) and of Gd(3+) coordinated to three different types of molecules: (i) a chelate (EDTA), (ii) structured peptides (mSE3/SE2) of the lanthanide-binding tags (LBTs) family with a single high-affinity binding site, and (iii) a calcium-binding protein (calmodulin) with four binding sites. The vibronic sideband (VSB) corresponding to the OH stretching mode of waters coordinated to Gd(3+), whose frequency is inversely correlated with the strength of the hydrogen bonding to neighboring waters, exhibits an increase in frequency when Gd(3+) becomes coordinated to either EDTA, calmodulin, or mSE3 peptide. In all of these cases, the addition of cation chloride or acetate salts to the solution increases the frequency of the vibronic band originating from the OH stretching mode of the coordinated waters in a cation- and concentration-dependent fashion. The cation dependence of the frequency increase scales with charge density of the cations, giving rise to an ordering consistent with the Hofmeister ordering. On the other hand, water Raman spectroscopy shows no significant change upon addition of these salts. Additionally, it is shown that the cation effect is modulated by the specific anion used. The results indicate a mechanism of action for Hofmeister series ions in which hydrogen bonding among hydration shell waters is modulated by several factors. High charge density cations sequester waters in a configuration that precludes strong hydrogen bonding to neighboring waters. Under such conditions, anion effects emerge as anions compete for hydrogen-bonding sites with the remaining free waters on the surface of the hydration shell. The magnitude of the anion effect is both cation and Gd(3+)-binding site specific.

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.

Spectroscopic study of substrate binding to the carbonmonoxy form of dehaloperoxidase from Amphitrite ornata

Dehaloperoxidase (DHP) is a globular heme enzyme found in the marine worm Amphitrite ornata that can catalyze the dehalogenation of halophenols to the corresponding quinones by using hydrogen peroxide as a cosubstrate. Its three-dimensional fold is surprisingly similar to that of the oxygen storage protein myoglobin (Mb). A key structural feature common to both DHP and Mb is the existence of multiple conformations of the distal histidine. In DHP, the conformational flexibility may be involved in promotion of substrate and cosubstrate entry and exit. Here we have explored the dynamics of substrate binding in DHP using Fourier transform infrared spectroscopy and flash photolysis. A number of discrete conformations at the active site were identified from the appearance of multiple CO absorbance bands in the infrared region of the spectrum. Upon photolysis at cryogenic temperatures, the CO molecules are trapped at docking sites within the protein matrix, as inferred from the appearance of several photoproduct bands characteristic of each site. Substrate binding stabilizes the protein by approximately 20 kJ/mol. The low yield of substrate-bound DHP at ambient temperature points toward a steric inhibition of substrate binding by carbon monoxide.

Bacteriorhodopsin photocycle at cryogenic temperatures reveals distributed barriers of conformational substates

The time course of thermal reactions after illumination of 100% humidified bacteriorhodopsin films was followed with FTIR spectroscopy between 125 and 195 K. We monitored the conversion of the initial photoproduct, K, to the next, L intermediate, and a shunt reaction of the L state directly back to the initial BR state. Both reactions can be described by either multiexponential kinetics, which would lead to apparent end-state mixtures that contain increasing amounts of the product, i.e., L or BR, with increasing temperature, or distributed kinetics. Conventional kinetic schemes that could account for the partial conversion require reversible reactions, branching, or parallel cycles. These possibilities were tested by producing K or L and monitoring their interconversion at a single temperature and by shifting the temperature upward or downward after an initial incubation and after their redistribution. The results are inconsistent with any conventional scheme. Instead, we attribute the partial conversions to the other alternative, distributed kinetics, observed previously in myoglobin, which arise from an ensemble of frozen conformational substates at the cryogenic temperatures. In this case, the time course of the reactions reflects the progressive depletion of distinct microscopic substates in the order of their increasing activation barriers, with a distribution width for K to L reaction of approximately 7 kJ/mol.

Different relaxations in myoglobin after photolysis

To clarify the interplay of kinetic hole-burning (KHB), structural relaxation, and ligand migration in myoglobin (Mb), we measured time-resolved absorption spectra in the Soret region after photolysis of carbon monoxide Mb (MbCO) in the temperature interval 120-260 K and in the time window 350 ns to 200 ms. The spectral contributions of both photolyzed (Mb*) and liganded Mb (MbCO) have been analyzed by taking into account homogeneous bandwidth, coupling to vibrational modes, and static conformational heterogeneity. We succeeded in separating the "time-dependent" spectral changes, and this work provides possibilities to identify the events in the process of ligand rebinding. KHB is dominant at T <190 K in both the Mb* and the MbCO components. For MbCO, conformational substates interconversion at higher temperatures tends to average out the KHB effect. At 230-260 K, whereas almost no shift is observed in the MbCO spectrum, a shift of the order of approximately 80 cm(-1) is observed in Mb*. We attribute this shift to protein relaxation coupled to ligand migration. The time dependence of the Mb* spectral shift is interpreted with a model that enables us to calculate the highly nonexponential relaxation kinetics. Fits of stretched exponentials to this kinetics yield Kohlrausch parameter values of 0.25, confirming the analogy between proteins and glasses.

Kinetic evidence for three photolyzable taxonomic conformational substates in oxymyoglobin

The kinetics of oxygen geminate binding with the taxonomic substates of MbO2 are reported. The maximum entropy method was used to analyze the rebinding kinetics of MbCO and MbO2 monitored in the Soret. The resulting rate distributions were found to consist of a small number of overlapping bands. A global parametric fit of a series of rate distributions recorded at several temperatures was performed using a Gaussian basis set to resolve the individual enthalpy distributions P(H). This approach was first validated by showing that the well-documented taxonomic substates of MbCO could be recovered. The method was then applied to MbO2. Three taxonomic substates were identified at pH 4.8, whereas only two of them contribute to oxygen geminate rebinding at pH 7.0. These findings show that, similarly to MbCO, MbO2 also exists as three photolyzable and kinetically different taxonomic substates and suggest reconsidering the issue of the photolysis quantum yield of MbO2.

The red-edge effects: 30 years of exploration

In 1970, three laboratories independently made a discovery that, for aromatic fluorophores embedded into different rigid and highly viscous media, the spectroscopic properties do not conform to classical rules. The fluorescence spectra can depend on excitation wavelength, and the excited-state energy transfer, if present, fails at the "red" excitation edge. These red-edge effects were related to the existence of excited-state distribution of fluorophores on their interaction energy with the environment and the slow rate of dielectric relaxation of this environment. In these conditions the site-selection can be provided by variation of the energy of illuminating light quanta, and the behaviour of selected species can be followed as a function of time and other variables. These observations found extensive application in different areas of research: colloid and polymer science, molecular biophysics, photochemistry and photobiology. In particular, they led to the development of very productive methods of studying the dynamics of dielectric relaxations in protein and membranes, using the tryptophan emission and the emission of a variety of probes. These studies were extended to the time domain with the observation of new site-selective effects in emission intensity and anisotropy decays. They stimulated the emergence and development of cryogenic energy-selective and single-molecular techniques that became valuable tools in their own right in chemistry and biophysics research. Site-selection effects were discovered for electron-transfer and proton-transfer reactions if they depended on the dynamics of the environment. This review is focused on the progress in the field of red-edge effects, their applications and prospects.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

This review is focused on the basic theory of diffusion-assisted reactions in microheterogeneous systems, from porous solids to self-organized colloids and biomolecules. Rich kinetic behaviors observed experimentally are explained in a unified fashion using simple concepts of competing distance and time scales of the reaction and the embedding structure. We mainly consider pseudo-first-order reactions, such as luminescence quenching, described by the Smoluchowski type of equation for the reactant pair distribution function with a sink term defined by the reaction mechanism. Microheterogeneity can affect the microscopic rate constant. It also enters the evolution equation through various spatial constraints leading to complicated boundary conditions and, possibly, to the reduction of dimensionality of the diffusion space. The reaction coordinate and diffusive motion along this coordinate are understood in a general way, depending on the problem at hand. Thus, the evolution operator can describe translational and rotational diffusion of molecules in a usual sense, it can be a discrete random walk operator when dealing with hopping of adsorbates in solids, or it can correspond to conformational fluctuations in proteins. Mathematical formulation is universal but physical consequences can be different. Understanding the principal features of reaction kinetics in microheterogeneous systems enables one to extract important structural and dynamical information about the host environments by analyzing suitably designed experiments, it helps building effective strategies for computer simulations, and ultimately opens possibilities for designing systems with controllable reactivity properties.

Thermal behavior of long wavelength absorption transitions in Spirulina platensis photosystem I trimers

In photosystem I trimers of Spirulina platensis a major long wavelength transition is irreversibly bleached by illumination with high-intensity white light. The photobleaching hole, identified by both absorption and circular dichroism spectroscopies, is interpreted as the inhomogeneously broadened Q(y) transition of a chlorophyll form that absorbs maximally near 709 nm at room temperature. Analysis of the mean square deviation of the photobleaching hole between 80 and 300 K, in the linear electron-phonon frame, indicates that the optical reorganization energy is 52 cm(-1), four times greater than that for the bulk, short-wavelength-absorbing chlorophylls, and the inhomogenous site distribution bandwidth is close to 150 cm(-1). The room temperature bandwidth, close to 18.5 nm, is dominated by thermal (homogeneous) broadening. Photobleaching induces correlated circular dichroism changes, of opposite sign, at 709 and 670 nm, which suggests that the long wavelength transition may be a low energy excitonic band, in agreement with its high reorganization energy. Clear identification of the 709-nm spectral form was used in developing a Gaussian description of the long wavelength absorption tail by analyzing the changing band shape during photobleaching using a global decomposition procedure. Additional absorption states near 720, 733, and 743 nm were thus identified. The lowest energy state at 743 nm is present in substoichiometric levels at room temperature and its presence was confirmed by fluorescence spectroscopy. This state displays an unusual increase in intensity upon lowering the temperature, which is successfully described by assuming the presence of low-lying, thermally populated states.

A thermal broadening study of the antenna chlorophylls in PSI-200, LHCI, and PSI core

The intact photosystem I of maize containing its full antenna complement (PSI-200) has been purified and fractionated into the core and outer antenna (LHCI) components. It is demonstrated by absorption and fluorescence spectroscopy that at least 80% of the long wavelength absorbing antenna pigments (red forms) are located in LHCI. Absorption spectra in the Qy region of all three preparations were measured between 72 and 300 K and subjected to a thermal broadening analysis. Data are interpreted in the linear electron-phonon coupling assumption, and the average optical reorganization energy (Snum) for the bulk pigment band and the red absorption tail determined. A marked asymmetry in Snum values across the absorption band is demonstrated. The bulk pigments in all three preparations have rather low values, in the range of 15-25 cm-1, suggesting that Stokes shifts for the absorption forms are in the 1. 5-3 nm range. On the other hand the red forms have markedly greater reorganization energies. While a direct thermal analysis of the red tail indicates minimum Snum values of around 60 cm-1, when the contribution of the red tail of the bulk pigments is corrected for in LHCI, the more reliable value of 110 cm-1 is obtained. These high Snum values for the red pigment forms suggest that they have unusually wide homogeneously broadened absorption bands and large Stokes shifts (6-11 nm).

Structural and dynamic properties of the homodimeric hemoglobin from Scapharca inaequivalvis Thr-72-->Ile mutant: molecular dynamics simulation, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy studies

Molecular dynamics simulations, low temperature visible absorption spectroscopy, and resonance Raman spectroscopy have been performed on a mutant of the Scapharca inaequivalvis homodimeric hemoglobin, where residue threonine 72, at the subunit interface, has been substituted by isoleucine. Molecular dynamics simulation indicates that in the Thr-72-->Ile mutant several residues that have been shown to play a role in ligand binding fluctuate around orientations and distances similar to those observed in the x-ray structure of the CO derivative of the native hemoglobin, although the overall structure remains in the T state. Visible absorption spectroscopy data indicate that in the deoxy form the Soret band is less asymmetric in the mutant than in the native protein, suggesting a more planar heme structure; moreover, these data suggest a similar heme-solvent interaction in both the liganded and unliganded states of the mutant protein, at variance with that observed in the native protein. The "conformation sensitive" band III of the deoxy mutant protein is shifted to lower energy by >100 cm-1 with respect to the native one, about one-half of that observed in the low temperature photoproducts of both proteins, indicating a less polar or more hydrophobic heme environment. Resonance Raman spectroscopy data show a slight shift of the iron-proximal histidine stretching mode of the deoxy mutant toward lower frequency with respect to the native protein, which can be interpreted in terms of either a change in packing of the phenyl ring of Phe-97, as also observed from the simulation, or a loss of water in the heme pocket. In line with this latter interpretation, the number of water molecules that dynamically enters the intersubunit interface, as calculated by the molecular dynamics simulation, is lower in the mutant than in the native protein. The 10-ns photoproduct for the carbonmonoxy mutant derivative has a higher iron-proximal histidine stretching frequency than does the native protein. This suggests a subnanosecond relaxation that is slowed in the mutant, consistent with a stabilization of the R structure. Taken together, the molecular dynamics and the spectroscopic data indicate that the higher oxygen affinity displayed by the Thr-72-->Ile mutant is mainly due to a local perturbation in the dimer interface that propagates to the heme region, perturbing the polarity of the heme environment and propionate interactions. These changes are consistent with a destabilization of the T state and a stabilization of the R state in the mutant relative to the native protein.

The protein state of matter

Three ways are generally used to visualize proteins: (1) a static model in which the atomic positions are defined, (2) a dynamic model taking into account fluctuations, and (3) a reactive model that reflects the internal and external electric fields of the molecule. The properties of chromophoric prosthetic groups can be probed by optical spectroscopy, and when high resolution techniques are used, the results reveal information about the local electric fields in proteins, as influenced and determined by atomic positions and dynamics.

A synthetic picture of intramolecular dynamics of proteins. Towards a contemporary statistical theory of biochemical processes

An increasing body of experimental evidence indicates the slow character of internal dynamics of native proteins. The important consequence of this is that theories of chemical reactions, used hitherto, appear inadequate for description of most biochemical reactions. Construction of a contemporary, truly advanced statistical theory of biochemical processes will need simple but realistic models of microscopic dynamics of biomolecules. In this review, intended to be a contribution towards this direction, three topics are considered. First, an intentionally simplified picture of dynamics of native proteins which emerges from recent investigations is presented. Fast vibrational modes of motion, of periods varying from 10(-14) to 10(-11) s, are contrasted with purely stochastic conformational transitions. Significant evidence is adduced that the relaxation time spectrum of the latter spreads in the whole range from 10(-11) to 10(5) s or longer, and up to 10(-7) s it is practically quasi-continuous. Next, the essential ideas of the theory of reaction rates based on stochastic models of intramolecular dynamics are outlined. Special attention is paid to reactions involving molecules in the initial conformational substrates confirmed to the transition state, which is realized in actual experimental situations. And finally, the two best experimentally justified classes of models of conformational transition dynamics, symbolically referred to as "protein glass" and "protein machine", are described and applied to the interpretation of a few simple biochemical processes, perhaps the most important result reported is the demonstration of the possibility of predominance of the short initial condition-dependent stage of protein involved reactions over the main stage described by the standard kinetics. This initial stage, and not the latter, is expected to be responsible for the coupling of component reactions in the complete enzymatic cycles as well as more complex processes of biological free energy transduction.

Structural factors controlling ligand binding to myoglobin: a kinetic hole-burning study

Using temperature-derivative spectroscopy in the temperature range below 100 K, we have studied the dependence of the Soret band on the recombination barrier in sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation at 12 K. The spectra were separated into contributions from the photodissociated species, Mb*CO, and CO-bound myoglobin. The line shapes of the Soret bands of both photolyzed and liganded myoglobin were analyzed with a model that takes into account the homogeneous bandwidth, coupling of the electronic transition to vibrational modes, and static conformational heterogeneity. The analysis yields correlations between the activation enthalpy for rebinding and the model parameters that characterize the homogeneous subensembles within the conformationally heterogeneous ensemble. Such couplings between spectral and functional parameters arise when they both originate from a common structural coordinate. This effect is frequently denoted as "kinetic hole burning." The study of these correlations gives direct insights into the structure-function relationship in proteins. On the basis of earlier work that assigned spectral parameters to geometric properties of the heme, the connections with the heme geometry are discussed. We show that two separate structural coordinates influence the Soret line shape, but only one of the two is coupled to the enthalpy barrier for rebinding. We give evidence that this coordinate, contrary to widespread belief, is not the iron displacement from the mean heme plane.

Fourier transform infrared analysis of the interaction of azide with the active site of oxidized and reduced bovine Cu,Zn superoxide dismutase

Binding of azide to the native and arginine-modified bovine Cu,Zn superoxide dismutase in the oxidized and reduced form and to the copper-free derivative has been investigated by Fourier transform infrared spectroscopy. The antisymmetric stretching band of the azide is shifted to higher energy upon coordination to the copper atom of the oxidized form of the native enzyme. Similar spectral changes occur upon interaction of the anion with the Cu-diethylenetriamine model compound. On the other hand, interaction of azide with the native reduced form of the enzyme results in a band shift toward lower energy with respect to the free anion band. The same shift is observed after reaction of the azide with free lysine or arginine but not when it is reacted with other amino acid residues. The antisymmetric band of the azide is not perturbed by addition of the reduced arginine-modified enzyme; it is likely shifted toward higher energy upon addition of oxidized arginine-modified enzyme while it is again shifted toward lower energy in the presence of the copper-free derivative of the unmodified enzyme. It is concluded that azide does not directly coordinate to the copper in the reduced form of Cu,Zn superoxide dismutase but it remains in the active-site pocket in electrostatic interaction with the guanidinium group of Arg141, which is an invariant residue in this class of enzymes.

Structural characterization of the myoglobin active site using infrared crystallography

We use polarized IR absorption on single crystals to determine the orientation of carbon monoxide bound at the active site of myoglobin, and conclude that the C-O bond lies approximately 7 degrees from the normal to the mean plane of the heme. This result disagrees with much larger angular displacements reported in structural models derived from X-ray and neutron diffraction measurements. The insensitivity of the IR-derived orientation to changes in pH or crystal packing contrasts with the wide variations in CO orientation among diffraction-based models and suggests that the latter are in error. The small energies required to displace the C-O bond 7 degrees from its energetically preferred upright geometry suggest that distortion of the surrounding protein, rather than the relatively undeformable Fe-C-O unit, is the main steric mechanism inhibiting CO binding to myoglobin.

Dynamics of carbon monoxide binding with cytochromes P-450

The dynamics of CO rebinding with cytochromes P-450cam, P-450scc, and P-450LM2 after laser flash photolysis have been investigated from 293 to 77 K, and the distribution functions of the rate parameters P(k) and of the activation enthalpy P(H) were determined using the maximum entropy method. In a fluid solvent, geminate rebinding is nonexponential, presumably because of a spectral shift induced by protein relaxation on the same time scale. Substrate binding increases the yield of the bimolecular process and decreases the bimolecular rate by 1 or 2 orders of magnitude. The amplitude of these effects seems to correlate with substrate specificity. In a rigid environment at low temperature, cytochromes P-450 exhibit a bimodal distribution of activation enthalpy; P(H) consists of two distinct bands which are in a thermal equilibrium even at 77 K. The results lead to a scheme in which a common structural perturbation splits the conformational substates of cytochromes P-450 into pairs of "doublet" substates with different dynamic properties. The hierarchy of conformational substates of cytochromes P-450 thus contrasts with that of oxygen-binding hemoproteins such as myoglobin.

Spectral diffusion and the energy landscape of a protein

We present a novel type of spectral diffusion experiment in the millikelvin range to characterize the energy landscape of a protein as compared with that of a glass. We measure the time evolution of spectral holes for more than 300 hr after well-defined initial nonequilibrium conditions. We show that the model of noninteracting two-level systems can describe spectral diffusion in the glass, but fails for the protein. Our results further demonstrate that randomness in the energy landscape of a protein shows features of organization. There are "deep minimum" states separated by barriers, the heights of which we are able to estimate. The energy landscape of a glass is featureless by comparison.

Stark effect experiments in cytochrome c-type proteins: structural hierarchies

We performed hole-burning Stark effect experiments on cytochrome c in which the iron of the herne was either removed or replaced by Zn. According to the experiments, the free-base compound has an effective inversion center, even in the protein. The Zn compound, on the other hand, shows quite peculiar features: in the low-frequency range of the inhomogeneous band, it definitely has a dipole moment, as indicated by a splitting of the hole in the external field. However, in the maximum of the inhomogeneous band, a severe charge redistribution occurs, as the experiments show. In addition to the Stark experiments, we performed calculations of the electrostatic fields at the pyrrole rings and at the metal site of the heme group. We interpret our findings with a model based on structural hierarchies: the protein can exist in a few subconformations, which can be distinguished through the structure of the heme pocket. The different pocket structures support different structures of the chromophore, which, in turn, can be distinguished through their behavior in an external field. These distinct structures, in turn, correspond to a rather broad distribution of protein structures, which leave, however, the pocket structure largely unchanged. These structures show up in inhomogeneous broadening.

Heme geometry in the 10 K photoproduct from sperm whale carbonmonoxymyoglobin

We have measured the Soret band of the photoproduct obtained by complete photolysis of sperm whale carbonmonoxymyoglobin at 10 K. The experimental spectrum has been modeled with an analytical expression that takes into account the homogeneous bandwidth, the coupling of the electronic transition with both high and low frequency vibrational modes, and the effects of static conformational heterogeneity. The comparison with deoxymyoglobin at low temperature reveals three main differences. In the photoproduct, the Soret band is shifted to red. The band is less asymmetric, and an enhanced coupling to the heme vibrational mode at 674 cm-1 is observed. These differences reflect incomplete relaxation of the active site after ligand dissociation. The smaller band asymmetry of the photoproduct can be explained by a smaller displacement of the iron atom from the mean porphyrin plane, in quantitative agreement with the X-ray structure analysis. The enhanced vibrational coupling is attributed to a subtle heme distortion from the planar geometry that is barely detectable in the X-ray structure.

Ligand binding and protein dynamics in cupredoxins

Type 1 copper sites bind nitric oxide (NO) in a photolabile complex. We have studied the NO binding properties of the type 1 copper sites in two cupredoxins, azurin and halocyanin, by measuring the temperature dependence of the ligand binding equilibria and the kinetics of the association reaction after photodissociation over a wide range of temperature (80-280 K) and time (10(-6)-10(2) s). In both proteins, we find nonexponential kinetics below 200 K that do not depend on the NO concentration. Consequently, this process is interpreted as geminate recombination. In azurin, the rebinding can be modeled with the Arrhenius law using a single pre-exponential factor of 10(8.3) s-1 and a Gaussian distribution of enthalpy barriers centered at 22 kJ/mol with a width [full width at half-maximum (FWHM)] of 11 kJ/mol. In halocyanin, a more complex behavior is observed. About 97% of the rebinding population can also be characterized by a Gaussian distribution of enthalpy barriers at 12 kJ/mol with a width of 6.0 kJ/mol (FWHM). The pre-exponential of this population is 1.6 x 10(12) s-1 at 100 K. After the majority population has rebound, a power-law phase that can be modeled with a gamma-distribution of enthalpy barriers is observed. Between 120 and 180 K, an additional feature that can be interpreted as a relaxation of the barrier distribution toward higher barriers shows up in the kinetics. Above 200 K, a slower, exponential rebinding appears in both cupredoxins. Since the kinetics depend on the NO concentration, this process is identified as bimolecular rebinding.(ABSTRACT TRUNCATED AT 250 WORDS)

Low temperature optical absorption spectroscopy: an approach to the study of stereodynamic properties of hemeproteins

In this short review we show how suitable analysis of the temperature dependence of the optical absorption spectra of metalloproteins can give insight into their stereodynamic properties in the region of the chromophore. To this end, the theory of coupling between an intense allowed electronic transition of a chromophore and Franck-Condon active vibrations of the nearby atoms is applied to the Soret band of hemeproteins to obtain an analytical expression suitable for fitting the spectral profile at various temperatures. The reported approach enables one to separate the various contributions to the overall bandwidth together with the parameters that characterize the vibrational coupling. The thermal behavior of these quantities gives information on the dynamic properties of the active site and on their dependence upon protein structure and ligation state. The Soret band of hemeproteins appears to be coupled to high frequency vibrational modes of the heme group (as already shown by resonance Raman spectroscopy) and to a "bath" of low frequency modes most likely deriving from the bulk of the protein. For the deoxy derivatives inhomogeneous broadening arising from conformational heterogeneity appears to contribute substantially to the linewidth. The data indicate the onset, at temperatures near 180 K, of large scale anharmonic motions that can be attributed to jumping among different conformational substates of the protein.

Thermal broadening of the Soret band in heme complexes and in heme-proteins: role of iron dynamics

We report the thermal broadening of the Soret band in heme-CO, heme-OH and protoporphyrin IX in the temperature range 300-20 K. For protoporphyrin IX the temperature dependent Gaussian line broadening follows the behavior predicted by the harmonic approximation in the entire temperature range investigated. In contrast, for heme-CO and heme-OH the harmonic behavior is obeyed only up to about 180 K and an anomalous line broadening increase is observed at higher temperatures. This effect is attributed to the onset of anharmonic motions of the iron atom with respect to the porphyrin plane. Comparison with previously reported analogous data for heme proteins enables us to suggest that the onset of substrate interconversions in these latter systems can be reflected in motions of the iron atom with respect to the porphyrin plane.

Ligand binding to heme proteins: the effect of light on ligand binding in myoglobin

Extended illumination slows the rebinding of CO to myoglobin after photodissociation at cryogenic temperatures. Two types of models have been put forward to explain the effect: motions of the CO within the heme pocket or conformational transitions of the protein. To resolve this ambiguity, we have studied the effect of extended illumination on ligand binding to horse and sperm whale myoglobin (hMb and swMb) with temperature-derivative spectroscopy, monitoring the reaction in the CO stretch bands in the infrared and the conformation-sensitive band III near 760 nm. The experiments show that the stretch frequency of the photodissociated CO does not change upon illumination, implying that the slowing of the CO rebinding is caused by conformational relaxation of Mb from the bound state toward the deoxy structure. The light-induced relaxation (LIR) depends on the number of photons absorbed but not on the light intensity or duration separately. LIR occurs on photon absorption in either the bound or photodissociated state and depends on the temperature at which the MbCO is illuminated. The LIR proceeds in jumps through a small number of conformational substates. The effective barrier for rebinding increases with each step. The substates populated are similar to those found in the thermally-induced relaxation (TIR) that is observed above 160 K. LIR depends markedly on the structural details; it differs for swMbCO and hMbCO and even for the three A substates of swMbCO. Pronounced differences exist between the effects in MbCO and MbO2. The similarity of LIR and TIR leads to a revised model for ligand binding to swMbCO and hMbCO, in which the relaxation is crucial for the escape of the ligand from the pocket, as was first suggested by Friedman [Friedman, J. M. (1985) Science 228, 1273-1280].

Stereodynamic properties of the cooperative homodimeric Scapharca inaequivalvis hemoglobin studied through optical absorption spectroscopy and ligand rebinding kinetics

The study of the thermal evolution of the Soret band in heme proteins has proved to be a useful tool to understand their stereodynamic properties; moreover, it enables one to relate protein matrix fluctuations and functional behavior when carried out in combination with kinetic experiments on carbon monoxide rebinding after flash photolysis. In this work, we report the thermal evolution of the Soret band of deoxy, carbonmonoxy, and nitric oxide derivatives of the cooperative homodimeric Scapharca inaequivalvis hemoglobin in the temperature range 10-300 K and the carbon monoxide rebinding kinetics after flash photolysis in the temperature range 60-200 K. The two sets of results indicate that Scapharca hemoglobin has a very rigid protein structure compared with other hemeproteins. This feature is brought out i) by the absence of nonharmonic contributions to the soft modes coupled to the Soret band in the liganded derivatives, and ii) by the almost "in plane" position of the iron atom in the photoproduct obtained approximately 10(-8) s after dissociating the bound carbon monoxide molecule at 15 K.

Conformational relaxation and ligand binding in myoglobin

Absorption spectroscopy with nanosecond time resolution shows that myoglobin undergoes conformational relaxation on the same time scale as geminate rebinding of carbon monoxide. Ligand rebinding following photodissociation of the heme-CO complex was measured from the amplitude of the average difference spectrum, while conformational changes were measured from changes in the detailed shape of the Soret spectra of the deoxyhemes. Experiments in which the solvent viscosity was varied between 1 and 300 cP and the temperature between 268 and 308 K were analyzed by fitting the multiwavelength kinetic data with both empirical and molecular models. Novel numerical techniques were employed in fitting the data, including the use of singular value decomposition to remove the effects of temperature and solvent on the spectra and of a Monte Carlo method to overcome the multiple minimum problem in searching parameter space. The molecular model is the minimal model that incorporates all of the major features of myoglobin kinetics at ambient temperatures, including a fast and slow rebinding conformation and two geminate states for each conformation. The results of fitting the kinetic data with this model indicate that the geminate-rebinding rates for the two conformations differ by at least a factor of 100. The differences between the spectra of the two conformations generated from the fits are similar to the differences between those of the R and T conformations of hemoglobin. In modeling the data, the dependence of the rates on temperature and viscosity was parametrized using a modification of Kramers theory which includes the contributions of both protein and solvent to the friction. The rate of the transition from the fast to the slow rebinding conformation is found to be inversely proportional to the viscosity when the viscosity exceeds about 30 cP and nearly viscosity independent at low viscosity. The viscosity dependence at high viscosities suggests that the two conformations differ by the global displacement of protein atoms on the proximal side of the heme observed by X-ray crystallography. We suggest that the conformational change observed in our experiments corresponds to the final portion of the nonexponential conformational relaxation recently observed by Anfinrud and co-workers, which begins on a picosecond time scale. Furthermore, extrapolation of our data to temperatures near that of the solvent glass transition suggests that this conformational relaxation may very well be the one postulated by Frauenfelder and co-workers to explain the decrease in the rate of geminate rebinding with increasing temperature above 180 K.

The transition from inhomogeneous to homogeneous kinetics in CO binding to myoglobin

Heme proteins react inhomogeneously with ligands at cryogenic temperatures and homogeneously at room temperature. We have identified and characterized a transition from inhomogeneous to homogeneous behavior at intermediate temperatures in the time dependence of CO binding to horse myoglobin. The turnover is attributed to a functionally important tertiary protein relaxation process during which the barrier increases dynamically. This is verified by a combination of theory and multipulse measurements. A likely biological significance of this effect is in the autocatalysis of the ligand release process.

Energy transfer in the inhomogeneously broadened core antenna of purple bacteria: a simultaneous fit of low-intensity picosecond absorption and fluorescence kinetics

The excited state decay kinetics of chromatophores of the purple photosynthetic bacterium Rhodospirillum rubrum have been recorded at 77 K using picosecond absorption difference spectroscopy under strict annihilation free conditions. The kinetics are shown to be strongly detection wavelength dependent. A simultaneous kinetic modeling of these experiments together with earlier fluorescence kinetics by numerical integration of the appropriate master equation is performed. This model, which accounts for the spectral inhomogeneity of the core light-harvesting antenna of photosynthetic purple bacteria, reveals three qualitatively distinct stages of excitation transfer with different time scales. At first a fast transfer to a local energy minimum takes place (approximately 1 ps). This is followed by a much slower transfer between different energy minima (10-30 ps). The third component corresponds to the excitation transfer to the reaction center, which depends on its state (60 and 200 ps for open and closed, respectively) and seems also to be the bottleneck in the overall trapping time. An acceptable correspondence between theoretical and experimental decay kinetics is achieved at 77 K and at room temperature by assuming that the width of the inhomogeneous broadening is 10-15 nm and the mean residence time of the excitation in the antenna lattice site is 2-3 ps.

Structural heterogeneity of the Fe(2+)-N epsilon (HisF8) bond in various hemoglobin and myoglobin derivatives probed by the Raman-active iron histidine stretching mode

We have examined the Fe(2+)-N epsilon (HisF8) complex in hemoglobin A (HbA) by measuring the band profile of its Raman-active nu Fe-His stretching mode at pH 6.4, 7.0, and 8.0 using the 441-nm line of a HeCd laser. A line shape analysis revealed that the band can be decomposed into five different sublines at omega 1 = 195 cm-1, omega 2 = 203 cm-1, omega 3 = 212 cm-1, omega 4 = 218 cm-1, and omega 5 = 226 cm-1. To identify these to the contributions from the different subunits we have reanalyzed the nu Fe-His band of the HbA hybrids alpha(Fe)2 beta(Co)2 and alpha(Co)2 beta(Fe)2 reported earlier by Rousseau and Friedman (D. Rousseau and J. M. Friedman. 1988. In Biological Application on Raman Spectroscopy. T. G. Spiro, editor, 133-216). Moreover we have reanalyzed other Raman bands from the literature, namely the nu Fe-His band of the isolated hemoglobin subunits alpha SH- and beta SH-HbA, various hemoglobin mutants (i.e., Hb(TyrC7 alpha-->Phe), Hb(TyrC7 alpha-->His), Hb M-Boston and Hb M-Iwate), N-ethylmaleimide-des(Arg141 alpha) hemoglobin (NES-des(Arg141 alpha)HbA) and photolyzed carbonmonoxide hemoglobin (Hb*CO) measured 25 ps and 10 ns after photolysis. These molecules are known to exist in different quaternary states. All bands can be decomposed into a set of sublines exhibiting frequencies which are nearly identical to those found for deoxyhemoglobin A. Additional sublines were found to contribute to the nu Fe-His band of NES-des(Arg141 alpha) HbA and the Hb*CO species. The peak frequencies of the bands are determined by the most intensive sublines. Moreover we have measured the nu Fe-His band of deoxyHbA at 10 K in an aqueous solution and in a 80% glycerol/water mixture. Its subline composition at this temperature depends on the solvent and parallels that of more R-like hemoglobin derivatives. We have also measured the optical charge transfer band III of deoxyHbA at room temperature and found, that at least three subbands are required to fit its asymmetric band shape. This corroborates the findings on the nu Fe-His band in that it is indicative of a heterogeneity of the Fe(2+)-N epsilon(HisF8) bond. Finally we measured the nu Fe-His band of horse heart deoxyMb at different temperatures and decomposed it into three different sublines. In accordance with what was obtained for HbA their intensities rather than their frequencies are temperature-dependent. By comparison with VFe-His bands of some Mb mutants (i.e., Mb(His E7.->Gly) and Mb(HisE7__*Met) we suggest that these sublines may be attributed to different conformations of the heme pocket. Our data show, that the V Fe-His band is governed by at least two different coordinates x and y determining its frequency and intensity, respectively. While the former can be assigned to the tilt angle theta between the Fe2+-NJ(HisF8) bond and the heme normal and/or to the displacement delta of the iron from the heme plane, variations in the intensity may be caused by changes of the azimuthal angle phi formed by the projection of the proximal imidazole and the N(l)-Fe2+-N(III) line of the heme. The sublines are therefore interpreted as resulting from different conformational substates of the Fe2+-N(HisFa) complex which differ in terms of x (theta and/or delta). Each of them may further be subdivided in sub-substates with respect to the coordinate y (theta). Quaternary and tertiary transitions of the protein alter the population of these substates thus giving rise to a redistribution of the VFe-HiS sublines which shifts the corresponding peak frequency to higher values.

Protein dynamics: conformational disorder, vibrational coupling and anharmonicity in deoxy-hemoglobin and myoglobin

In this work we study the temperature dependence of the Soret band lineshape of deoxymyoglobin and deoxyhemoglobin, in the range 300-20 K. To fit the measured spectra we use an approach originally proposed by Champion and coworkers (Srajer et al. 1986; Srajer and Champion 1991). The band profile is modelled as a Voigt function that accounts for the coupling with low frequency vibrational modes, whereas the coupling with high frequency modes is responsible for the vibronic structure of the spectra. Moreover, owing to the position of the iron atom out of the mean heme plane, inhomogeneous broadening brings about a non-Gaussian distribution of 0-0 electronic transition frequencies. The reported analysis enables us to isolate the various contributions to the overall bandwidth, and their temperature dependence points out the relevance of low frequency vibrations and of large scale anharmonic motions starting at temperatures higher than 170 K. Information on the mean iron-heme plane distance and on its temperature dependence, as well as on the heme pocket conformational disorder, is also obtained.

Vibronic energy map and excited state vibrational characteristics of magnesium myoglobin determined by energy-selective fluorescence

The vibrational frequencies of the singlet excited state of Mg-substituted myoglobin and relative absorption probabilities were determined by fluorescence line-narrowing spectroscopy. These spectra contain information on the structure of the excited state species, and the availability of vibrationally resolved spectra from excited state biomolecules should aid in elucidating their structure and reactivity.

Conformational substates in azurin

Azurin is a small blue copper protein in the electron transfer chain of denitrifying bacteria. It forms a photolabile complex with nitric oxide (NO) at low temperatures. We studied the temperature dependence of the ligand binding equilibrium and the kinetics of the association reaction after photodissociation over a wide range of temperature (80-280 K) and time (10(-6)-10(2) s). The nonexponential rebinding below 200 K is independent of the NO concentration and is interpreted as internal recombination. The rebinding can be modeled with the Arrhenius law by using a single preexponential factor of 6.3 x 10(8) s-1 and a Gaussian distribution of enthalpy barriers centered at 23 kJ/mol with a width of 11 kJ/mol. Above 200 K, a slower, exponential rebinding process appears. The dependence of the kinetics on the NO concentration characterizes this reaction as bimolecular rebinding. The binding kinetics of NO to azurin show impressive analogies to the binding of carbon monoxide to myoglobin. We conclude that conformational substates occur not only in heme proteins but also in proteins with different active sites and secondary structures.

Does biocatalysis involve inhomogeneous kinetics?

Biocatalytic reactions can occur according to two very different mechanisms: homogeneous, which is described by standard transition state theory (TST) and its modifications, and inhomogeneous (polychromatic), which is characteristic for some of the charge-transfer reactions in liquids and amorphous solids. While most data published on enzyme reactions are interpreted on the basis of homogeneous kinetics, the important recent findings suggest the involvement of inhomogeneous kinetics mechanism.

Protein dynamics. Vibrational coupling, spectral broadening mechanisms, and anharmonicity effects in carbonmonoxy heme proteins studied by the temperature dependence of the Soret band lineshape

In this work we study the temperature dependence of the Soret band lineshape of the carbonmonoxy derivatives of sperm whale myoglobin, human hemoglobin, and its isolated alpha and beta subunits. To fit the observed spectral profile we use an analytic expression derived for a system whereby a single electronic transition is coupled to Franck-Condon active vibrational modes, within the adiabatic and harmonic approximation. The vibronic structure of the spectra arises from the coupling with high frequency modes; these modes contribute to the total line shape through a series of Lorentzians with peak positions at vibrational overtones and half width related to the time constant of the population decay of the excited electronic state (homogeneous broadening); moreover, the coupling with low frequency modes broadens each Lorentzian to a Voigtian. Inhomogeneous broadening is modeled as a gaussian distribution of the 0-0 transition frequencies and is therefore added as a constant term to the previous gaussian width. This spectral deconvolution enables us to investigate the different contributions to line broadening and the parameters that characterize the vibrational coupling, as well as their dependence upon protein and solvent composition. The investigation is carried out as a function of temperature in the range 20-300 K; relevant information is obtained by comparing experimental results with theoretical predictions. This work supports a description of the investigated proteins as heterogeneous systems, whose heterogeneity depends on the particular protein and on the composition of the external matrix. The delocalized pi electron cloud of the porphyrin ring is coupled not only to the high frequency vibrational modes of the active site but also to a "bath" of lower frequency modes that involve the entire protein; moreover at suitable temperatures (approximately 200 K), anharmonic motions, which are an obvious prerequisite for the jumping among different conformational substates, become evident.