Temperature-dependent heme kinetics with nonexponential binding and barrier relaxation in the absence of protein conformational substates

Two-Dimensional Fluctuation Correlation Spectroscopy (2D-FlucCS): A Method to Determine the Origin of Relaxation Rate Dispersion

Relaxation rate dispersion, i.e., nonexponential or multicomponent kinetics, is observed in complex systems when measuring relaxation kinetics. Often, the origin of rate dispersion is associated with the heterogeneity in the system. However, both homogeneous (where all molecules experience the same rate but inherently nonexponential) and heterogeneous (where all molecules experience different rates) systems can exhibit rate dispersion. A multidimensional correlation analysis method has been demonstrated to detect and quantify rate dispersion observed in molecular rotation, diffusion, solvation, and reaction kinetics. One-dimensional (1D) autocorrelation function detects rate dispersion and measures its extent. Two-dimensional (2D) autocorrelation function measures the origin of rate dispersion and distinguishes homogeneous from heterogeneous. In a heterogeneous system, implicitly there exist subensembles of molecules experiencing different rates. A three-dimensional (3D) autocorrelation function measures subensemble exchange if present and reveals if the system possesses static or dynamic heterogeneity. This perspective discusses the principles, applications, and potential and also presents a future outlook of two-dimensional fluctuation correlation spectroscopy (2D-FlucCS). The method is applicable to any experiment or simulation where a time series of fluctuation in an observable (emission, scattering, current, etc.) around a mean value can be obtained in steady state (equilibrium or nonequilibrium), provided the system is ergodic.

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.

Solvent Composition Drives the Rebinding Kinetics of Nitric Oxide to Microperoxidase

The rebinding kinetics of NO after photodissociation from microperoxidase (Mp-9) is studied in different solvent environments. In mixed glycerol/water (G/W) mixtures the dissociating ligand rebinds with a yield close to 1 due to the cavities formed by the solvent whereas in pure water the ligand can diffuse into the solvent after photodissociation. In the G/W mixture, only geminate rebinding on the sub-picosecond and 5 ps time scales was found and the rebinding fraction is unity which compares well with available experiments. Contrary to that, simulations in pure water find two time scales – ~10 ps and ~200 ps - indicating that both, geminate rebinding and rebinding after diffusion of NO in the surrounding water contribute. The rebinding fraction is around 0.63 within 1 ns which is in stark contrast with experiment. Including ions (Na and Cl) at 0.15 M concentration in water leads to rebinding kinetics tending to that in the glycerol/water mixture and yields agreement with experiments. The effect of temperature is also probed and found to be non-negligible. The present simulations suggest that NO rebinding in Mp is primarily driven by thermal fluctuations which is consistent with recent resonance Raman spectroscopy experiments and simulations on MbNO.

A mechanism for CO regulation of ion channels

Despite being highly toxic, carbon monoxide (CO) is also an essential intracellular signalling molecule. The mechanisms of CO-dependent cell signalling are poorly defined, but are likely to involve interactions with heme proteins. One such role for CO is in ion channel regulation. Here, we examine the interaction of CO with K_ATP channels. We find that CO activates K_ATP channels and that heme binding to a CXXHX₁₆H motif on the SUR2A receptor is required for the CO-dependent increase in channel activity. Spectroscopic and kinetic data were used to quantify the interaction of CO with the ferrous heme-SUR2A complex. The results are significant because they directly connect CO-dependent regulation to a heme-binding event on the channel. We use this information to present molecular-level insight into the dynamic processes that control the interactions of CO with a heme-regulated channel protein, and we present a structural framework for understanding the complex interplay between heme and CO in ion channel regulation.

Ultrafast CO Kinetics in Heme Proteins: Adiabatic Ligand Binding and Heavy Atom Tunneling

The ultrafast kinetics of CO rebinding to carbon monoxide oxidation activator protein (ChCooA) are measured over a wide temperature range and compared with the kinetics of CO binding in other heme systems such as myoglobin (Mb) and hemoglobin (Hb). The Arrhenius prefactor for CO binding to ChCooA and protoheme (∼10¹¹ s^-1) is similar to what is found for spin-allowed NO binding to heme proteins and is several orders of magnitude larger than the prefactor of Mb and Hb (∼10⁹ s^-1). This indicates that the CO binding reaction is adiabatic, in contrast to the commonly held view that it is nonadiabatic due to spin-forbidden (ΔS = 2) selection rules. Under the adiabatic condition, entropic factors, rather than spin-selection rules, are the source of the reduced Arrhenius prefactors associated with CO binding in Mb and Hb. The kinetic response of ChCooA-CO is nonexponential at all temperatures, including 298 K, and is described quantitatively using a distribution of enthalpic rebinding barriers associated with heterogeneity in the heme doming conformation. Above the solvent glass transition (T_g ∼ 180 K), the rebinding progress slows as temperature increases, and this is ascribed to an evolution of the distribution toward increased heme doming and larger enthalpic barriers. Between T_g and ∼60 K, the nonexponential rebinding slows down as the temperature is lowered and the survival fraction follows the predictions expected for a quenched barrier distribution. Below ∼60 K the rebinding kinetics do not follow these predictions unless quantum mechanical tunneling along the heme doming coordinate is also included as an active channel for CO binding.

The Dynamics Behind the Affinity: Controlling Heme-Gas Affinity via Geminate Recombination and Heme Propionate Conformation in the NO Carrier Cytochrome c'

Nitric oxide (NO) sensors are heme proteins which may also bind CO and O₂. Control of heme-gas affinity and their discrimination are achieved by the structural properties and reactivity of the heme and its distal and proximal environments, leading to several energy barriers. In the bacterial NO sensor cytochrome c' from Alcaligenes xylosoxidans (AXCP), the single Leu16Ala distal mutation boosts the affinity for gas ligands by a remarkable 10⁶-10⁸-fold, transforming AXCP from one of the lowest affinity gas binding proteins to one of the highest. Here, we report the dynamics of diatomics after photodissociation from wild type and L16A-AXCP over 12 orders of magnitude in time. For the L16A variant, the picosecond geminate rebinding of both CO and NO appears with an unprecedented 100% yield, and no exit of these ligands from protein to solvent could be observed. Molecular dynamic simulations saliently demonstrate that dissociated CO stays within 4 Å from Fe²⁺, in contrast to wild-type AXCP. The L16A mutation confers a heme propionate conformation and docking site which traps the diatomics, maximizing the probability of recombination and directly explaining the ultrahigh affinities for CO, NO, and O₂. Overall, our results point to a novel mechanism for modulating heme-gas affinities in proteins.

Rebinding dynamics of NO to microperoxidase-8 probed by time-resolved vibrational spectroscopy

Femtosecond vibrational spectroscopy was used to probe the rebinding kinetics of NO to microperoxidase-8 (Mp), an ideal model system for the active site of ligand-binding heme proteins, including myoglobin and hemoglobin, after the photodeligation of MpNO in glycerol/water (G/W) solutions at 294 K. The geminate rebinding (GR) of NO to Mp in viscous solutions was highly efficient and ultrafast and negligibly dependent on the solution viscosity, which was adjusted by changing the glycerol content from 65% to 90% by volume in G/W mixtures. The kinetics of the GR of NO to Mp in viscous solutions was well represented by an exponential function with a time constant of ca. 11 ps. Although the kinetic traces of the GR of NO to Mp in solutions with three different viscosities (18, 81, and 252 cP) almost overlap, they show a slight difference early in the decay process. The kinetic traces were also described by the diffusion-controlled reaction theory with a Coulomb potential. Since the ligand is deligated in a neutral form, an ionic pair of NO(-) and Mp(+) may be produced before forming the Mp-NO bond by an electron transfer from Mp to NO as the deligated NO is sufficiently near to the Fe atom of Mp. The strong reactivity between NO and ferrous heme may arise from the Coulomb interaction between the reacting pair, which is consistent with the harpooning mechanism for NO binding to heme.

Investigations of the low frequency modes of ferric cytochrome c using vibrational coherence spectroscopy

Femtosecond vibrational coherence spectroscopy is used to investigate the low frequency vibrational dynamics of the electron transfer heme protein, cytochrome c (cyt c). The vibrational coherence spectra of ferric cyt c have been measured as a function of excitation wavelength within the Soret band. Vibrational coherence spectra obtained with excitation between 412 and 421 nm display a strong mode at ∼44 cm^–1 that has been assigned to have a significant contribution from heme ruffling motion in the electronic ground state. This assignment is based partially on the presence of a large heme ruffling distortion in the normal coordinate structural decomposition (NSD) analysis of the X-ray crystal structures. When the excitation wavelength is moved into the ∼421–435 nm region, the transient absorption increases along with the relative intensity of two modes near ∼55 and 30 cm^–1. The intensity of the mode near 44 cm^–1 appears to minimize in this region and then recover (but with an opposite phase compared to the blue excitation) when the laser is tuned to 443 nm. These observations are consistent with the superposition of both ground and excited state coherence in the 421–435 nm region due to the excitation of a weak porphyrin-to-iron charge transfer (CT) state, which has a lifetime long enough to observe vibrational coherence. The mode near 55 cm^–1 is suggested to arise from ruffling in a transient CT state that has a less ruffled heme due to its iron d⁶ configuration.

Computational and theoretical methods for protein folding

A computational approach is essential whenever the complexity of the process under study is such that direct theoretical or experimental approaches are not viable. This is the case for protein folding, for which a significant amount of data are being collected. This paper reports on the essential role of in silico methods and the unprecedented interplay of computational and theoretical approaches, which is a defining point of the interdisciplinary investigations of the protein folding process. Besides giving an overview of the available computational methods and tools, we argue that computation plays not merely an ancillary role but has a more constructive function in that computational work may precede theory and experiments. More precisely, computation can provide the primary conceptual clues to inspire subsequent theoretical and experimental work even in a case where no preexisting evidence or theoretical frameworks are available. This is cogently manifested in the application of machine learning methods to come to grips with the folding dynamics. These close relationships suggested complementing the review of computational methods within the appropriate theoretical context to provide a self-contained outlook of the basic concepts that have converged into a unified description of folding and have grown in a synergic relationship with their computational counterpart. Finally, the advantages and limitations of current computational methodologies are discussed to show how the smart analysis of large amounts of data and the development of more effective algorithms can improve our understanding of protein folding.

Investigations of heme ligation and ligand switching in cytochromes p450 and p420

It is generally accepted that the inactive P420 form of cytochrome P450 (CYP) involves the protonation of the native cysteine thiolate to form a neutral thiol heme ligand. On the other hand, it has also been suggested that recruitment of a histidine to replace the native cysteine thiolate ligand might underlie the P450 → P420 transition. Here, we discuss resonance Raman investigations of the H93G myoglobin (Mb) mutant in the presence of tetrahydrothiophene (THT) or cyclopentathiol (CPSH), and on pressure-induced cytochrome P420cam (CYP101), that show a histidine becomes the heme ligand upon CO binding. The Raman mode near 220 cm⁻¹, normally associated with the Fe-histidine vibration in heme proteins, is not observed in either reduced P420cam or the reduced H93G Mb samples, indicating that histidine is not the ligand in the reduced state. The absence of a mode near 220 cm⁻¹ is also inconsistent with a generalization of the suggestion that the 221 cm⁻¹ Raman mode, observed in the P420-CO photoproduct of inducible nitric oxide synthase (iNOS), arises from a thiol-bound ferrous heme. This leads us to assign the 218 cm⁻¹ mode observed in the 10 ns P420cam-CO photoproduct Raman spectrum to a Fe-histidine vibration, in analogy to many other histidine-bound heme systems. Additionally, the inverse correlation plots of the νFe-His and νCO frequencies for the CO adducts of P420cam and the H93G analogs provide supporting evidence that histidine is the heme ligand in the P420-CO-bound state. We conclude that, when CO binds to the ferrous P420 state, a histidine ligand is recruited as the heme ligand. The common existence of an HXC-Fe motif in many CYP systems allows the C → H ligand switch to occur with only minor conformational changes. One suggested conformation of P420-CO involves the addition of another turn in the proximal L helix so that, when the protonated Cys ligand is dissociated from the heme, it can become part of the helix, and the heme is ligated by the His residue from the adjoining loop region. In other systems, such as iNOS and CYP3A4 (where the HXC-Fe motif is not found), a somewhat larger conformational change would be necessary to recuit a nearby histidine.

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.

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.

Spectroscopic analysis of protein Fe–NO complexes

The toxic free radical NO (nitric oxide) has diverse biological roles in eukaryotes and bacteria, being involved in signalling, vasodilation, blood clotting and immunity, and as an intermediate in microbial denitrification. The predominant biological mechanism of detecting NO is through the formation of iron nitrosyl complexes, although this is a deleterious process for other iron-containing enzymes. We have previously applied techniques such as UV–visible and EPR spectroscopy to the analysis of protein Fe–NO complex formation in order to study how NO controls the activity of the bacterial transcriptional regulators NorR and NsrR. These studies have analysed NO-dependent biological activity both in vitro and in vivo using diverse biochemical, molecular and spectroscopic methods. Recently, we have applied ultrafast 2D-IR (two-dimensional IR) spectroscopy to the analysis of NO–protein interactions using Mb (myoglobin) and Cc (cytochrome c) as model haem proteins. The ultrafast fluctuations of Cc and Mb show marked differences, indicating altered flexibility of the haem pockets. We have extended this analysis to bacterial catalase enzymes that are known to play a role in the nitrosative stress response by detoxifying peroxynitrite. The first 2D-IR analysis of haem nitrosylation and perspectives for the future are discussed.

Investigation of the low frequency dynamics of heme proteins: native and mutant cytochrome P450(cam) and redox partner complexes

Vibrational coherence spectroscopy (VCS) is used to investigate the low-frequency dynamics of camphor-free and camphor-bound cytochrome P450(cam) (CYP 101) and its L358P mutant. The low-frequency heme vibrations are found to be perturbed upon binding to the electron transfer partner putidaredoxin (Pdx). A strong correlation between the "detuned" vibrational coherence spectrum, which monitors frequencies between 100 and 400 cm(-1), and the lower frequency part of the Raman spectrum is also demonstrated. The very low frequency region ≤200 cm(-1), uniquely accessed by open-band VCS measurements, reveals a mode near 103 cm(-1) in P450(cam) when camphor is not present in the distal pocket. This reflects the presence of a specific heme distortion, such as saddling or ruffling, in the substrate-free state where water is coordinated to the low-spin iron atom. Such distortions are likely to retard the rate of electron transfer to the substrate-free protein. The presence of strong mode near ∼33 cm(-1) in the camphor-bound form suggests a significant heme-doming distortion, which is supported by analysis using normal coordinate structural decomposition. Pdx also displays a strong coherent vibration near 30 cm(-1) that in principle could be involved in vibrational resonance with its electron transfer target. A splitting of the 33 cm(-1) feature and intensification of a mode near 78 cm(-1) appear when the P450(cam)/Pdx complex is formed. These observations are consistent with vibrational mixing and heme geometric distortions upon Pdx binding that are coincident with the increased thiolate electron donation to the heme. The appearance of a mode near 65 cm(-1) in the coherence spectra of the L358P mutant is comparable to the mode at 78 cm(-1) seen in the P450(cam)/Pdx complex and is consistent with the view that the heme and its environment in the L358P mutant are similar to the Pdx-bound native protein. Resonance Raman spectra are presented for both P450(cam) and the L358P mutant and the changes are correlated with an increased amount of thiolate electron donation to the heme in the mutant sample.

Relationship between protein structural fluctuations and rebinding dynamics in ferric haem nitrosyls

The interaction of nitric oxide (NO) with haem proteins is widespread in biology. In the current paper, we present the first ultrafast 2D-IR (two-dimensional infrared) spectroscopic analysis of haem nitrosylation, which has been combined with time-resolved IR pump–probe studies to investigate the relationship between equilibrium vibrational dynamics of the haem environment and ligand rebinding behaviour following photolysis of NO from the Fe(III)–NO site. Studies of two haem proteins, Mb (myoglobin) and Cc (cytochrome c), which play different physiological roles, reveal marked contrasts in the ultrafast fluctuations of the protein pockets containing the haem, showing that the Mb pocket is somewhat more flexible than that of Cc. This correlates strongly with slower observed photolysis rebinding kinetics of Mb–NO compared with Cc–NO, and indicates a direct link between ultrafast fluctuations and biological functionality. Furthermore, this indicates the validity of linear response theories in relation to protein ligand binding. Finally, 2D-IR shows that Cc–NO displays two distinct structural sub-sites at room temperature that do not exchange on the timescales accessible via the NO vibrational lifetime.

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.

Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).

Novel mechanism for temperature-independent transitions in flexible molecules: role of thermodynamic fluctuations

A novel physical mechanism is proposed to explain the temperature-independent transition reactions in molecular systems. The mechanism becomes effective in the case of conformation transitions between quasi-isoenergetic molecular states. It is shown that at room temperatures, stochastic broadening of molecular energy levels predominates the energy of low-frequency vibrations accompanying the transition. This leads to a cancellation of temperature dependence in the stochastically averaged rate constants. As an example, a physical interpretation of temperature-independent onset of P2X{3} receptor desensitization in neuronal membranes is provided.

Measurements of heme relaxation and ligand recombination in strong magnetic fields

Heme cooling signals and diatomic ligand recombination kinetics are measured in strong magnetic fields (up to 10 T). We examined diatomic ligand recombination to heme model compounds (NO and CO), myoglobin (NO and O(2)), and horseradish peroxidase (NO). No magnetic field induced rate changes in any of the samples were observed within the experimental detection limit. However, in the case of CO binding to heme in glycerol and O(2) binding to myoglobin, we observe a small magnetic field dependent change in the early time amplitude of the optical response that is assigned to heme cooling. One possibility, consistent with this observation, is that there is a weak magnetic field dependence of the nonradiative branching ratio into the vibrationally hot electronic ground state during CO photolysis. Ancillary studies of the "spin-forbidden" CO binding reaction in a variety of heme compounds in the absence of magnetic field demonstrate a surprisingly wide range for the Arrhenius prefactor. We conclude that CO binding to heme is not always retarded by unfavorable spin selection rules involving a double spin-flip superexchange mechanism. In fact, it appears that the small prefactor ( approximately 10(9) s(-1)) found for CO rebinding to Mb may be anomalous, rather than the general rule for heme-CO rebinding. These results point to unresolved fundamental issues that underlie the theory of heme-ligand photolysis and rebinding.

Low-frequency mode activity of heme: femtosecond coherence spectroscopy of iron porphine halides and nitrophorin

The low-frequency mode activity of metalloporphyrins has been studied for iron porphine-halides (Fe(P)(X), X = Cl, Br) and nitrophorin 4 (NP4) using femtosecond coherence spectroscopy (FCS) in combination with polarized resonance Raman spectroscopy and density functional theory (DFT). It is confirmed that the mode symmetry selection rules for FCS are the same as for Raman scattering and that both Franck-Condon and Jahn-Teller mode activities are observed for Fe(P)(X) under Soret resonance conditions. The DFT-calculated low-frequency (20-400 cm (-1)) modes, and their frequency shifts upon halide substitution, are in good agreement with experimental Raman and coherence data, so that mode assignments can be made. The doming mode is located at approximately 80 cm (-1) for Fe(P)(Cl) and at approximately 60 cm (-1) for Fe(P)(Br). NP4 is also studied with coherence techniques, and the NO-bound species of ferric and ferrous NP4 display a mode at approximately 30-40 cm (-1) that is associated with transient heme doming motion following NO photolysis. The coherence spectra of three ferric derivatives of NP4 with different degrees of heme ruffling distortion are also investigated. We find a mode at approximately 60 cm (-1) whose relative intensity in the coherence spectra depends quadratically on the magnitude of the ruffling distortion. To quantitatively account for this correlation, a new "distortion-induced" Raman enhancement mechanism is presented. This mechanism is unique to low-frequency "soft modes" of the molecular framework that can be distorted by environmental forces. These results demonstrate the potential of FCS as a sensitive probe of dynamic and functionally important nonplanar heme vibrational excitations that are induced by the protein environmental forces or by the chemical reactions in the aqueous phase.

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.

Coherence spectroscopy investigations of the low-frequency vibrations of heme: effects of protein-specific perturbations

Femtosecond coherence spectroscopy is used to probe the low-frequency (20-200 cm(-1)) vibrational modes of heme proteins in solution. Horseradish peroxidase (HRP), myoglobin (Mb), and Campylobacter jejuni globin (Cgb) are compared and significant differences in the coherence spectra are revealed. It is concluded that hydrogen bonding and ligand charge do not strongly affect the low-frequency coherence spectra and that protein-specific deformations of the heme group lower its symmetry and control the relative spectral intensities. Such deformations potentially provide a means for proteins to tune heme reaction coordinates, so that they can perform a broad array of specific functions. Native HRP displays complex spectral behavior above approximately 50 cm(-1) and very weak activity below approximately 50 cm(-1). Binding of the substrate analog, benzhydroxamic acid, leads to distinct changes in the coherence and Raman spectra of HRP that are consistent with the stabilization of a heme water ligand. The CN derivatives of the three proteins are studied to make comparisons under conditions of uniform heme coordination and spin-state. MbCN is dominated by a doming mode near 40 cm(-1), while HRPCN displays a strong oscillation at higher frequency (96 cm(-1)) that can be correlated with the saddling distortion observed in the X-ray structure. In contrast, CgbCN displays low-frequency coherence spectra that contain strong modes near 30 and 80 cm(-1), probably associated with a combination of heme doming and ruffling. HRPNO displays a strong doming mode near 40 cm(-1) that is activated by photolysis. The damping of the coherent motions is significantly reduced when the heme is shielded from solvent fluctuations by the protein material and reduced still further when T approximately < 50 K, as pure dephasing processes due to the protein-solvent phonon bath are frozen out.

Investigations of vibrational coherence in the low-frequency region of ferric heme proteins

Femtosecond coherence spectroscopy is applied to a series of ferric heme protein samples. The low-frequency vibrational spectra that are revealed show dominant oscillations near 40 cm(-1). MbCN is taken as a typical example of a histidine-ligated, six-coordinate, ferric heme and a comprehensive spectroscopic analysis is carried out. The results of this analysis reveal a new heme photoproduct species, absorbing near 418 nm, which is consistent with the photolysis of the His(93) axial ligand. The photoproduct undergoes subsequent rebinding/recovery with a time constant of approximately 4 ps. The photoproduct lineshapes are consistent with a photolysis quantum yield of 75-100%, although the observation of a relatively strong six-coordinate heme coherence near 252 cm(-1) (assigned to nu(9) in the MbCN Raman spectrum) suggests that the 75% lower limit is much more likely. The phase and amplitude excitation profiles of the low-frequency mode at 40 cm(-1) suggest that this mode is strongly coupled to the MbCN photoproduct species and it is assigned to the doming mode of the transient penta-coordinated material. The absolute phase of the 40 cm(-1) mode is found to be pi/2 on the red side of 418 nm and it jumps to 3pi/2 as excitation is tuned to the blue side of 418 nm. The absolute phase of the 40 cm(-1) signal is not explained by the standard theory for resonant impulsive stimulated Raman scattering. New mechanisms that give a dominant momentum impulse to the resonant wavepacket, rather than a coordinate displacement, are discussed. The possibilities of heme iron atom recoil after photolysis, as well as ultrafast nonradiative decay, are explored as potential ways to generate the strong momentum impulse needed to understand the phase properties of the 40 cm(-1) mode.

Low frequency spectral density of ferrous heme: perturbations induced by axial ligation and protein insertion

Femtosecond coherence spectroscopy is used to probe low frequency (20-400 cm(-1)) modes of the ferrous heme group in solution, with and without 2-methyl imidazole (2MeIm) as an axial ligand. The results are compared to heme proteins (CPO, P450(cam), HRP, Mb) where insertion of the heme into the protein results in redistribution of the low frequency spectral density and in (approximately 60%) longer damping times for the coherent signals. The major effect of imidazole ligation to the ferrous heme is the "softening" of the low frequency force constants by a factor of approximately 0.6 +/- 0.1. The functional consequences of imidazole ligation are assessed and it is found that the enthalpic CO rebinding barrier is increased significantly when imidazole is bound. The force constant softening analysis, combined with the kinetics results, indicates that the iron is displaced by only approximately 0.2 A from the heme plane in the absence of the imidazole ligand, whereas it is displaced by approximately 0.4 A when imidazole (histidine) is present. This suggests that binding of imidazole (histidine) as an axial ligand, and the concomitant softening of the force constants, leads to an anharmonic distortion of the heme group that has significant functional consequences.