Ultrafast Dynamics and Vibrational Relaxation in Six-Coordinate Heme Proteins Revealed by Femtosecond Stimulated Raman Spectroscopy

Identifying the structural rearrangements during photoinduced reactions is a fundamental challenge for understanding from a microscopic perspective the dynamics underlying the functional mechanisms of heme proteins. Here, femtosecond stimulated Raman spectroscopy is applied to follow the ultrafast evolution of two different proteins, each bearing a six-coordinate heme with two amino acid axial ligands. By exploiting the sensitivity of Raman spectra to the structural configuration, we investigate the effects of photolysis and the binding of amino acid residues in cytochrome c and neuroglobin. By comparing the system response for different time delays and Raman pump resonances, we show how detailed properties of atomic motions and energy redistribution can be unveiled. In particular, we demonstrate substantially faster energy flow from the dissociated heme to the protein moiety in cytochrome c, which we assign to the presence of covalent heme–protein bonds.

Adiabatic Ligand Binding in Heme Proteins: Ultrafast Kinetics of Methionine Rebinding in Ferrous Cytochrome c

The dynamics of methionine geminate recombination following photodissociation in ferrous cytochrome c is investigated over a broad temperature range. The kinetic response, above the solvent glass transition ( T_g), is nearly monoexponential and displays a weak temperature dependence. Below T_g, the rebinding kinetics are nonexponential and can be explained using a quenched distribution of enthalpic rebinding barriers, arising from a relatively narrow distribution of heme out-of-plane displacements. The Arrhenius prefactor of this (Δ S = 2) reaction is ∼10¹¹ s^-1, which is similar to what has been found for the (Δ S = 1) NO binding reaction in heme proteins. This observation, along with other examples of ultrafast CO binding, provides strong evidence that ligand binding to heme is an adiabatic reaction with a spin-independent prefactor. In order to simultaneously account for the adiabatic nature of the reaction as well as the temperature dependence of both ultrafast CO and methionine geminate rebinding, it is proposed that a spin triplet state intersects and strongly couples to the reactant ( S = 2) and product ( S = 0) state surfaces in the transition state region along the reaction coordinate. It is also suggested that the nature of the intersecting triplet state and the reaction path may depend upon the proximity of the photolyzed ligand relative to the iron atom. At temperatures below ∼60 K, the kinetic data suggest that there is either an unexpected retardation of the heme photoproduct relaxation or that heavy atom quantum mechanical tunneling becomes significant.

Kinetic Analysis of a Globin-Coupled Histidine Kinase, AfGcHK: Effects of the Heme Iron Complex, Response Regulator, and Metal Cations on Autophosphorylation Activity

The globin-coupled histidine kinase, AfGcHK, is a part of the two-component signal transduction system from the soil bacterium Anaeromyxobacter sp. Fw109-5. Activation of its sensor domain significantly increases its autophosphorylation activity, which targets the His183 residue of its functional domain. The phosphate group of phosphorylated AfGcHK is then transferred to the cognate response regulator. We investigated the effects of selected variables on the autophosphorylation reaction's kinetics. The kcat values of the heme Fe(III)-OH(-), Fe(III)-cyanide, Fe(III)-imidazole, and Fe(II)-O2 bound active AfGcHK forms were 1.1-1.2 min(-1), and their Km(ATP) values were 18.9-35.4 μM. However, the active form bearing a CO-bound Fe(II) heme had a kcat of 1.0 min(-1) but a very high Km(ATP) value of 357 μM, suggesting that its active site structure differs strongly from the other active forms. The Fe(II) heme-bound inactive form had kcat and Km(ATP) values of 0.4 min(-1) and 78 μM, respectively, suggesting that its low activity reflects a low affinity for ATP relative to that of the Fe(III) form. The heme-free form exhibited low activity, with kcat and Km(ATP) values of 0.3 min(-1) and 33.6 μM, respectively, suggesting that the heme iron complex is essential for high catalytic activity. Overall, our results indicate that the coordination and oxidation state of the sensor domain heme iron profoundly affect the enzyme's catalytic activity because they modulate its ATP binding affinity and thus change its kcat/Km(ATP) value. The effects of the response regulator and different divalent metal cations on the autophosphorylation reaction are also discussed.

Catalytic enhancement of the heme-based oxygen-sensing phosphodiesterase EcDOS by hydrogen sulfide is caused by changes in heme coordination structure

EcDOS is a heme-based O2-sensing phosphodiesterase in which O2 binding to the heme iron complex in the N-terminal domain substantially enhances catalysis toward cyclic-di-GMP, which occurs in the C-terminal domain. Here, we found that hydrogen sulfide enhances the catalytic activity of full-length EcDOS, possibly owing to the admixture of 6-coordinated heme Fe(III)-SH(-) and Fe(II)-O2 complexes generated during the reaction. Alanine substitution at Met95, the axial ligand for the heme Fe(II) complex, converted the heme Fe(III) complex into the heme Fe(III)-SH(-) complex, but the addition of Na2S did not further reduce it to the heme Fe(II) complex of the Met95Ala mutant, and no subsequent formation of the heme Fe(II)-O2 complex was observed. In contrast, a Met95His mutant formed a stable heme Fe(II)-O2 complex in response to the same treatment. An Arg97Glu mutant, containing a glutamate substitution at the amino acid that interacts with O2 in the heme Fe(II)-O2 complex, formed a stable heme Fe(II) complex in response to Na2S, but this complex failed to bind O2. Interestingly, the addition of Na2S promoted formation of verdoheme (oxygen-incorporated, modified protoporphyrin IX) in an Arg97Ile mutant. Catalytic enhancement by Na2S was similar for Met95 mutants and the wild type, but significantly lower for the Arg97 mutants. Thus, this study shows the first isolation of spectrometrically separated, stable heme Fe(III)-SH(-), heme Fe(II) and heme Fe(II)-O2 complexes of full-length EcDOS with Na2S, and confirms that external-ligand-bound, 6-coordinated heme Fe(III)-SH(-) or heme Fe(II)-O2 complexes critically contribute to the Na2S-induced catalytic enhancement of EcDOS.

Dynamics of the heme-binding bacterial gas-sensing dissimilative nitrate respiration regulator (DNR) and activation barriers for ligand binding and escape

DNR (dissimilative nitrate respiration regulator) is a heme-binding transcription factor that is involved in the regulation of denitrification in Pseudomonas aeruginosa. In the ferrous deoxy state, the heme is 6-coordinate; external NO and CO can replace an internal ligand. Using fluorescence anisotropy, we show that high-affinity sequence-specific DNA binding occurs only when the heme is nitrosylated, consistent with the proposed function of DNR as NO sensor and transcriptional activator. This role is moreover supported by the NO "trapping" properties revealed by ultrafast spectroscopy that are similar to those of other heme-based NO sensor proteins. Dissociated CO-heme pairs rebind in an essentially barrierless way. This process competes with migration out of the heme pocket. The latter process is thermally activated (Ea ∼ 7 kJ/mol). This result is compared with other heme proteins, including the homologous CO sensor/transcription factor CooA, variants of the 5-coordinate mycobacterial sensor DosT and the electron transfer protein cytochrome c. This comparison indicates that thermal activation of ligand escape from the heme pocket is specific for systems where an external ligand replaces an internal one. The origin of this finding and possible implications are discussed.

Heme-based globin-coupled oxygen sensors: linking oxygen binding to functional regulation of diguanylate cyclase, histidine kinase, and methyl-accepting chemotaxis

An emerging class of novel heme-based oxygen sensors containing a globin fold binds and senses environmental O2 via a heme iron complex. Structure-function relationships of oxygen sensors containing a heme-bound globin fold are different from those containing heme-bound PAS and GAF folds. It is thus worth reconsidering from an evolutionary perspective how heme-bound proteins with a globin fold similar to that of hemoglobin and myoglobin could act as O2 sensors. Here, we summarize the molecular mechanisms of heme-based oxygen sensors containing a globin fold in an effort to shed light on the O2-sensing properties and O2-stimulated catalytic enhancement observed for these proteins.

The Heme-Based Oxygen-Sensor Phosphodiesterase Ec DOS (DosP): Structure-Function Relationships

Escherichia coli Direct Oxygen Sensor (Ec DOS, also known as Ec DosP) is a heme-based O2-sensing phosphodiesterase from Escherichia coli that catalyzes the conversion of cyclic-di-GMP to linear di-GMP. Cyclic-di-GMP is an important second messenger in bacteria, highlighting the importance of understanding structure-function relationships of Ec DOS. Ec DOS is composed of an N-terminal heme-bound O2-sensing PAS domain and a C-terminal phosphodiesterase catalytic domain. Notably, its activity is markedly enhanced by O2 binding to the heme Fe(II) complex in the PAS sensor domain. X-ray crystal structures and spectroscopic and catalytic characterization of the wild-type and mutant proteins have provided important structural and functional clues to understanding the molecular mechanism of intramolecular catalytic regulation by O2 binding. This review summarizes the intriguing findings that have obtained for Ec DOS.

Primary processes in heme-based sensor proteins

A wide and still rapidly increasing range of heme-based sensor proteins has been discovered over the last two decades. At the molecular level, these proteins function as bistable switches in which the catalytic activity of an enzymatic domain is altered mostly by binding or dissociation of small gaseous ligands (O2, NO or CO) to the heme in a sensor domain. The initial "signal" at the heme level is subsequently transmitted within the protein to the catalytic site, ultimately leading to adapted expression levels of specific proteins. Making use of the photolability of the heme-ligand bond that mimics thermal dissociation, early processes in this intra-protein signaling pathway can be followed using ultrafast optical spectroscopic techniques; they also occur on timescales accessible to molecular dynamics simulations. Experimental studies performed over the last decade on proteins including the sensors FixL (O2), CooA (CO) and soluble guanylate cyclase (NO) are reviewed with an emphasis on emerging general mechanisms. After heme-ligand bond breaking, the ligand can escape from the heme pocket and eventually from the protein, or rebind directly to the heme. Remarkably, in all sensor proteins the rebinding, specifically of the sensed ligand, is highly efficient. This "ligand trap" property possibly provides means to smoothen the effects of fast environmental fluctuations on the switching frequency. For 6-coordinate proteins, where exchange between an internal heme-bound residue and external gaseous ligands occurs, the study of early processes starting from the unliganded form indicates that mobility of the internal ligand may facilitate signal transfer. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Effect of DNA binding on geminate CO recombination kinetics in CO-sensing transcription factor CooA

Carbon monoxide oxidation activator (CooA) proteins are heme-based CO-sensing transcription factors. Here we study the ultrafast dynamics of geminate CO rebinding in two CooA homologues, Rhodospirillum rubrum (RrCooA) and Carboxydothermus hydrogenoformans (ChCooA). The effects of DNA binding and the truncation of the DNA-binding domain on the CO geminate recombination kinetics were specifically investigated. The CO rebinding kinetics in these CooA complexes take place on ultrafast time scales but remain non-exponential over many decades in time. We show that this non-exponential kinetic response is due to a quenched enthalpic barrier distribution resulting from a distribution of heme geometries that is frozen or slowly evolving on the time scale of CO rebinding. We also show that, upon CO binding, the distal pocket of the heme in the CooA proteins relaxes to form a very efficient hydrophobic trap for CO. DNA binding further tightens the narrow distal pocket and slightly weakens the iron-proximal histidine bond. Comparison of the CO rebinding kinetics of RrCooA, truncated RrCooA, and DNA-bound RrCooA proteins reveals that the uncomplexed and inherently flexible DNA-binding domain adds additional structural heterogeneity to the heme doming coordinate. When CooA forms a complex with DNA, the flexibility of the DNA-binding domain decreases, and the distribution of the conformations available in the heme domain becomes restricted. The kinetic studies also offer insights into how the architecture of the heme environment can tune entropic barriers in order to control the geminate recombination of CO in heme proteins, whereas spin selection rules play a minor or non-existent role.

Ultrafast ligand dynamics in the heme-based GAF sensor domains of the histidine kinases DosS and DosT from Mycobacterium tuberculosis

The transcriptional regulator DosR from M. tuberculosis plays a crucial role in the virulence to dormancy transition of the pathogen. DosR can be activated by DosT and DosS, two histidine kinases with heme-containing sensor GAF domains, capable of diatomic ligand binding. To investigate the initial processes occurring upon ligand dissociation, we performed ultrafast time-resolved absorption spectroscopy of the isolated sensor domains ligated with O(2), NO, and CO. The results reveal a relatively closed heme pocket for both proteins. For DosT the yield of O(2) escape from the heme pocket on the picoseconds time scale upon photodissociation was found to be very low (1.5%), similar to other heme-based oxygen sensor proteins, implying that this sensor acts as an effective O(2) trap. Remarkably, this yield is an order of magnitude higher in DosS (18%). For CO, by contrast, the fraction of CO rebinding within the heme pocket is higher in DosS. Experiments with mutant DosT sensor domains and molecular dynamics simulations indicate an important role in ligand discrimination of the distal tyrosine, present in both proteins, which forms a hydrogen bond with heme-bound O(2). We conclude that despite their similarity, DosT and DosS display ligand-specific different primary dynamics during the initial phases of intraprotein signaling. The distal tyrosine, present in both proteins, plays an important role in these processes.

Vibrational coherence spectroscopy of the heme domain in the CO-sensing transcriptional activator CooA

Femtosecond vibrational coherence spectroscopy was used to investigate the low-frequency vibrational dynamics of the heme in the carbon monoxide oxidation activator protein (CooA) from the thermophilic anaerobic bacterium Carboxydothermus hydrogenoformans (Ch-CooA). Low frequency vibrational modes are important because they are excited by the ambient thermal bath (k(B)T = 200 cm(-1)) and participate in thermally activated barrier crossing events. However, such modes are nearly impossible to detect in the aqueous phase using traditional spectroscopic methods. Here, we present the low frequency coherence spectra of the ferric, ferrous, and CO-bound forms of Ch-CooA in order to compare the protein-induced heme distortions in its active and inactive states. Distortions take place predominantly along the coordinates of low-frequency modes because of their weak force constants, and such distortions are reflected in the intensity of the vibrational coherence signals. A strong mode near ~90 cm(-1) in the ferrous form of Ch-CooA is suggested to contain a large component of heme ruffling, consistent with the imidazole-bound ferrous heme crystal structure, which shows a significant protein-induced heme distortion along this coordinate. A mode observed at ~228 cm(-1) in the six-coordinate ferrous state is proposed to be the ν(Fe-His) stretching vibration. The observation of the Fe-His mode indicates that photolysis of the N-terminal α-amino axial ligand takes place. This is followed by a rapid (~8.5 ps) transient absorption recovery, analogous to methionine rebinding in photolyzed ferrous cytochrome c. We have also studied CO photolysis in CooA, which revealed very strong photoproduct state coherent oscillations. The observation of heme-CO photoproduct oscillations is unusual because most other heme systems have CO rebinding kinetics that are too slow to make the measurement possible. The low frequency coherence spectrum of the CO-bound form of Ch-CooA shows a strong vibration at ~230 cm(-1) that is broadened and up-shifted compared to the ν(Fe-His) of Rr-CooA (216 cm(-1)). We propose that the stronger Fe-His bond is related to the enhanced thermal stability of Ch-CooA and that there is a smaller (time dependent) tilt of the histidine ring with respect to the heme plane in Ch-CooA. The appearance of strong modes at ~48 cm(-1) in both the ferrous and CO-bound forms of Ch-CooA is consistent with coupling of the heme doming distortion to the photolysis reaction in both samples. Upon CO binding and protein activation, a heme mode near 112 ± 5 cm(-1) disappears, probably indicating a decreased heme saddling distortion. This reflects changes in the heme environment and geometry that must be associated with the conformational transition activating the DNA-binding domain. Protein-specific DNA binding to the CO-bound form of Ch-CooA was also investigated, and although the CO rebinding kinetics are significantly perturbed, there are negligible changes in the low-frequency vibrational spectrum of the heme.

Investigations of low-frequency vibrational dynamics and ligand binding kinetics of cystathionine beta-synthase

Vibrational coherence spectroscopy is used to study the low frequency dynamics of the truncated dimer of human cystathionine beta-synthase (CBS). CBS is a pyridoxal-5'-phosphate-dependent heme enzyme with cysteine and histidine axial ligands that catalyzes the condensation of serine and homocysteine to form cystathionine. A strong correlation between the "detuned" coherence spectrum (which probes higher frequencies) and the Raman spectrum is demonstrated, and a rich pattern of modes below 200 cm(-1) is revealed. Normal coordinate structural decomposition (NSD) of the ferric CBS crystal structure predicts the enhancement of normal modes with significant heme "doming", "ruffling", and "saddling" content, and they are observed in the coherence spectra near approximately 40, approximately 60, and approximately 90 cm(-1). When pH is varied, the relative intensities and frequencies of the low frequency heme modes indicate the presence of a unique protein-induced heme structural perturbation near pH 7 that differs from what is observed at higher or lower pH. For ferric CBS, we observe a new mode near approximately 25 cm(-1), possibly involving the response of the protein, which exhibits a phase jump of approximately pi for excitation on the blue and red side of the Soret band maximum. The low frequency vibrational coherence spectrum of ferrous CBS is also presented, along with our efforts to probe its NO-bound complex. The CO geminate rebinding kinetics of CBS are similar to the CO-bound form of the gene activator protein CooA, but with the appearance of a significant additional kinetic inhomogeneity. Analysis of this inhomogeneity suggests that it arises from the two subunits of CBS and leads to a factor of approximately 20 for the ratio of the average CO geminate rebinding rates of the two subunits.

Heme ligand binding properties and intradimer interactions in the full-length sensor protein dos from Escherichia coli and its isolated heme domain

Dos from Escherichia coli is a bacterial gas sensor protein comprising a heme-containing gas sensor domain and a phosphodiesterase catalytic domain. Using a combination of static light scattering and gel filtration experiments, we established that, as are many other sensor proteins, the full-length protein is dimeric. The full-length dimer (association constant <10 nm) is more stable than the dimeric heme domain (association constant approximately 1 mum), and the dimer interface presumably includes both sensor and catalytic domains. Ultrafast spectroscopic studies showed little influence of the catalytic domain on kinetic processes in the direct vicinity of the heme. By contrast, the properties of ligand (CO and O(2)) binding to the heme in the sensor domain, occurring on a microsecond to second time scale, were found to be influenced by (i) the presence of the catalytic domain, (ii) the dimerization state, and in dimers, (iii) the ligation state of the other subunit. These results imply allosteric interactions within dimers. Steady-state titrations demonstrated marked cooperativity in oxygen binding to both the full-length protein and the isolated heme domain, a feature not reported to date for any dimeric sensor protein. Analysis of a variety of time-resolved experiments showed that Met-95 plays a major role in the intradimer interactions. The intrinsic binding and dissociation rates of Met-95 to the heme were modulated approximately 10-fold by intradimer and sensor-catalytic domain interactions. Dimerization effects were also observed for cyanide binding to the ferric heme domains, suggesting a similar role for Met-95 in ferric proteins.

Arg97 at the heme-distal side of the isolated heme-bound PAS domain of a heme-based oxygen sensor from Escherichia coli (Ec DOS) plays critical roles in autoxidation and binding to gases, particularly O2

The catalytic activity of heme-regulated phosphodiesterase from Escherichia coli (Ec DOS) on cyclic di-GMP is markedly enhanced upon binding of gas molecules, such as O2 and CO, to the heme iron complex in the sensor domain. Arg97 interacts directly with O2 bound to Fe(II) heme in the crystal structure of the isolated heme-bound sensor domain with the PAS structure (Ec DOS-PAS) and may thus be critical in ligand recognition. To establish the specific role of Arg97, we generated Arg97Ala, Arg97Glu, and Arg97Ile mutant Ec DOS-PAS proteins and examined binding to O2, CO, and cyanide, as well as redox potentials. The autoxidation rates of the Arg97Ala and Arg97Glu mutant proteins were up to 2000-fold higher, while the O2 dissociation rate constant for dissociation from the Fe(II)-O2 heme complex of the Arg97Ile mutant was 100-fold higher than that of the wild-type protein. In contrast, the redox potential values of the mutant proteins were only slightly different from that of the wild type (within 10 mV). Accordingly, we propose that Arg97 plays critical roles in recognition of the O2 molecule and redox switching by stabilizing the Fe(II)-O2 complex, thereby anchoring O2 to the heme iron and lowering the autoxidation rate to prevent formation of Fe(III) hemin species not regulated by gas molecules. Arg97 mutations significantly influenced interactions with the internal ligand Met95, during CO binding to the Fe(II) complex. Moreover, the binding behavior of cyanide to the Fe(III) complexes of the Arg mutant proteins was similar to that of O2, which is evident from the Kd values, suggestive of electrostatic interactions between cyanide and Arg97.

Ultrafast heme-residue bond formation in six-coordinate heme proteins: implications for functional ligand exchange

A survey is presented of picosecond kinetics of heme-residue bond formation after photolysis of histidine, methionine, or cysteine, in a broad range of ferrous six-coordinate heme proteins. These include human neuroglobin, a bacterial heme-binding superoxide dismutase (SOD), plant cytochrome b 559, the insect nuclear receptor E75, horse heart cytochrome c and the heme domain of the bacterial sensor protein Dos. We demonstrate that the fastest and dominant phase of binding of amino acid residues to domed heme invariably takes place with a time constant in the narrow range of 5-7 ps. Remarkably, this is also the case in the heme-binding SOD, where the heme is solvent-exposed. We reason that this fast phase corresponds to barrierless formation of the heme-residue bond from a configuration close to the bound state. Only in proteins where functional ligand exchange occurs, additional slower rebinding takes place on the time scale of tens of picoseconds after residue dissociation. We propose that the presence of these slower phases reflects flexibility in the heme environment that allows external ligands (O2, CO, NO, . . .) to functionally replace the internal residue after thermal dissociation of the heme-residue bond.

Low-frequency dynamics of Caldariomyces fumago chloroperoxidase probed by femtosecond coherence spectroscopy

Ultrafast laser spectroscopy techniques are used to measure the low-frequency vibrational coherence spectra and nitric oxide rebinding kinetics of Caldariomyces fumago chloroperoxidase (CPO). Comparisons of the CPO coherence spectra with those of other heme species are made to gauge the protein-specific nature of the low-frequency spectra. The coherence spectrum of native CPO is dominated by a mode that appears near 32-33 cm(-1) at all excitation wavelengths, with a phase that is consistent with a ground-state Raman-excited vibrational wavepacket. On the basis of a normal coordinate structural decomposition (NSD) analysis, we assign this feature to the thiolate-bound heme doming mode. Spectral resolution of the probe pulse ("detuned" detection) reveals a mode at 349 cm(-1), which has been previously assigned using Raman spectroscopy to the Fe-S stretching mode of native CPO. The ferrous species displays a larger degree of spectral inhomogeneity than the ferric species, as reflected by multiple shoulders in the optical absorption spectra. The inhomogeneities are revealed by changes in the coherence spectra at different excitation wavelengths. The appearance of a mode close to 220 cm(-1) in the coherence spectrum of reduced CPO excited at 440 nm suggests that a subpopulation of five coordinated histidine-ligated hemes is present in the ferrous state at a physiologically relevant pH. A significant increase in the amplitude of the coherence signal is observed for the resonance with the 440 nm subpopulation. Kinetics measurements reveal that nitric oxide binding to ferric and ferrous CPO can be described as a single-exponential process, with rebinding time constants of 29.4 +/- 1 and 9.3 +/- 1 ps, respectively. This is very similar to results previously reported for nitric oxide binding to horseradish peroxidase.