Photoacid Dynamics in the Green Fluorescent Protein.

The photoacid dynamics of fluorescent proteins include both electronic excited- and ground-state mechanisms of proton transfer. The associated characteristic timescales of these reactions range over many orders of magnitude, and the tunneling, barrier crossing, and relevant thermodynamics have in certain cases been linked to coherent nuclear motion. We review the literature and summarize the experiments and theory that demonstrate proton tunneling in the electronic ground state of the green fluorescent protein (GFP). We also discuss the excited-state proton-transfer reaction of GFP that takes place on the picosecond timescale. Although this reaction has been investigated using several vibrational spectroscopic methods, the interpretation remains unsettled. We discuss recent advances as well as remaining questions, in particular those related to the vibrational mode couplings that involve low-frequency modulations of chromophore vibrations on the timescale of proton transfer.

Electrical unfolding of cytochrome c during translocation through a nanopore constriction.

Many small proteins move across cellular compartments through narrow pores. In order to thread a protein through a constriction, free energy must be overcome to either deform or completely unfold the protein. In principle, the diameter of the pore, along with the effective driving force for unfolding the protein, as well as its barrier to translocation, should be critical factors that govern whether the process proceeds via squeezing, unfolding/threading, or both. To probe this for a well-established protein system, we studied the electric-field-driven translocation behavior of cytochrome c (cyt c) through ultrathin silicon nitride (SiNx) solid-state nanopores of diameters ranging from 1.5 to 5.5 nm. For a 2.5-nm-diameter pore, we find that, in a threshold electric-field regime of ∼30 to 100 MV/m, cyt c is able to squeeze through the pore. As electric fields inside the pore are increased, the unfolded state of cyt c is thermodynamically stabilized, facilitating its translocation. In contrast, for 1.5- and 2.0-nm-diameter pores, translocation occurs only by threading of the fully unfolded protein after it transitions through a higher energy unfolding intermediate state at the mouth of the pore. The relative energies between the metastable, intermediate, and unfolded protein states are extracted using a simple thermodynamic model that is dictated by the relatively slow (∼ms) protein translocation times for passing through the nanopore. These experiments map the various modes of protein translocation through a constriction, which opens avenues for exploring protein folding structures, internal contacts, and electric-field-induced deformability.

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 (∼1011 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 (∼109 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 (Tg ∼ 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 Tg 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.

Proton-Coupled Electron Transfer and the "Linear Approximation" for Coupling to the Donor-Acceptor Distance Fluctuations.

The often-used "linear approximation" for treating the coupling of proton donor-acceptor (D-A) distance fluctuations to proton-coupled electron transfer tunneling reactions is systematically examined. The accuracy of this approximation is found to depend on the potential energy surfaces that are used to describe both the tunneling particle vibrations and the D-A coordinate probability distribution. Harmonic treatment of both the tunneling particle and the D-A coordinates results in a significant breakdown of the linear approximation when the width of the D-A distribution exceeds ∼0.05 Å. When a symmetric back-to-back Morse potential is used to describe the tunneling particle vibrations in the reactant and product states, the D-A distribution width can be expanded up to ∼0.09 Å before the rates calculated using the linear approximation exceed the exact result by an order of magnitude. Incorporation of a more realistic anharmonic D-A potential, based on quantum calculations, includes the important electronic D-A repulsion energy so that the sampling of short D-A distances is reduced. This approach improves the linear approximation as long as the D-A distribution has a width less than ∼0.1 Å. The conditions for the validity of the linear approximation are found to be more fragile when the calculated kinetic isotope effect (KIE) and its temperature dependence are also taken into account.

Investigations of heme distortion, low-frequency vibrational excitations, and electron transfer in cytochrome c.

Cytochrome (cyt) c is an important electron transfer protein. The ruffling deformation of its heme cofactor has been suggested to relate to its electron transfer rate. However, there is no direct experimental evidence demonstrating this correlation. In this work, we studied Pseudomonas aeruginosa cytochrome c551 and its F7A mutant. These two proteins, although similar in their X-ray crystal structure, display a significant difference in their heme out-of-plane deformations, mainly along the ruffling coordinate. Resonance Raman and vibrational coherence measurements also indicate significant differences in ruffling-sensitive modes, particularly the low-frequency γa mode found between ∼50-60 cm(-1). This supports previous assignments of γa as having a large ruffling content. Measurement of the photoreduction kinetics finds an order of magnitude decrease of the photoreduction cross-section in the F7A mutant, which has nearly twice the ruffling deformation as the WT. Additional measurements on cytochrome c demonstrate that heme ruffling is correlated exponentially with the electron transfer rates and suggest that ruffling could play an important role in redox control. A major relaxation of heme ruffling in cytochrome c, upon binding to the mitochondrial membrane, is discussed in this context.

Deep proton tunneling in the electronically adiabatic and non-adiabatic limits: comparison of the quantum and classical treatment of donor-acceptor motion in a protein environment.

Analytical models describing the temperature dependence of the deep tunneling rate, useful for proton, hydrogen, or hydride transfer in proteins, are developed and compared. Electronically adiabatic and non-adiabatic expressions are presented where the donor-acceptor (D-A) motion is treated either as a quantized vibration or as a classical "gating" distribution. We stress the importance of fitting experimental data on an absolute scale in the electronically adiabatic limit, which normally applies to these reactions, and find that vibrationally enhanced deep tunneling takes place on sub-ns timescales at room temperature for typical H-bonding distances. As noted previously, a small room temperature kinetic isotope effect (KIE) does not eliminate deep tunneling as a major transport channel. The quantum approach focuses on the vibrational sub-space composed of the D-A and hydrogen atom motions, where hydrogen bonding and protein restoring forces quantize the D-A vibration. A Duschinsky rotation is mandated between the normal modes of the reactant and product states and the rotation angle depends on the tunneling particle mass. This tunnel-mass dependent rotation contributes substantially to the KIE and its temperature dependence. The effect of the Duschinsky rotation is solved exactly to find the rate in the electronically non-adiabatic limit and compared to the Born-Oppenheimer (B-O) approximation approach. The B-O approximation is employed to find the rate in the electronically adiabatic limit, where we explore both harmonic and quartic double-well potentials for the hydrogen atom bound states. Both the electronically adiabatic and non-adiabatic rates are found to diverge at high temperature unless the proton coupling includes the often neglected quadratic term in the D-A displacement from equilibrium. A new expression is presented for the electronically adiabatic tunnel rate in the classical limit for D-A motion that should be useful to experimentalists working near room temperature. This expression also holds when a broad protein conformational distribution of D-A equilibrium distances dominates the spread of the D-A vibrational wavefunction.

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 identification of reactive porphyrin motions.

Nuclear resonance vibrational spectroscopy (NRVS) reveals the vibrational dynamics of a Mössbauer probe nucleus. Here, (57)Fe NRVS measurements yield the complete spectrum of Fe vibrations in halide complexes of iron porphyrins. Iron porphine serves as a useful symmetric model for the more complex spectrum of asymmetric heme molecules that contribute to numerous essential biological processes. Quantitative comparison with the vibrational density of states (VDOS) predicted for the Fe atom by density functional theory calculations unambiguously identifies the correct sextet ground state in each case. These experimentally authenticated calculations then provide detailed normal mode descriptions for each observed vibration. All Fe-ligand vibrations are clearly identified despite the high symmetry of the Fe environment. Low frequency molecular distortions and acoustic lattice modes also contribute to the experimental signal. Correlation matrices compare vibrations between different molecules and yield a detailed picture of how heme vibrations evolve in response to (a) halide binding and (b) asymmetric placement of porphyrin side chains. The side chains strongly influence the energetics of heme doming motions that control Fe reactivity, which are easily observed in the experimental signal.

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).

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.

Oxygen-Induced In Situ Manipulation of the Interlayer Coupling and Exciton Recombination in Bi2Se3/MoS2 2D Heterostructures.

Two-dimensional (2D) heterostructures are more than a sum of the parent 2D materials, but are also a product of the interlayer coupling, which can induce new properties. In this paper, we present a method to tune the interlayer coupling in Bi2Se3/MoS2 2D heterostructures by regulating the oxygen presence in the atmosphere, while applying laser or thermal energy. Our data suggest that the interlayer coupling is tuned through the diffusive intercalation and deintercalation of oxygen molecules. When one layer of Bi2Se3 is grown on monolayer MoS2, an influential interlayer coupling is formed, which quenches the signature photoluminescence (PL) peaks. However, thermally treating in the presence of oxygen disrupts the interlayer coupling, facilitating the emergence of the MoS2 PL peak. Our density functional theory calculations predict that intercalated oxygen increases the interlayer separation ∼17%, disrupting the interlayer coupling and inducing the layers to behave more electronically independent. The interlayer coupling can then be restored by thermally treating in N2 or Ar, where the peaks will requench. Hence, this is an interesting oxygen-induced switching between "non-radiative" and "radiative" exciton recombination. This switching can also be accomplished locally, controllably, and reversibly using a low-power focused laser, while changing the environment from pure N2 to air. This allows for the interlayer coupling to be precisely manipulated with submicron spatial resolution, facilitating site-programmable 2D light-emitting pixels whose emission intensity could be precisely varied by a factor exceeding 200×. Our results show that these atomically thin 2D heterostructures may be excellent candidates for oxygen sensing.

Environmental filtering of native and non-native stream macrophyte assemblages by habitat disturbances in an agricultural landscape.

Understanding how inter-specific variation in functional traits affects native and non-native species responses to stream disturbances, is necessary to inform management strategies, providing tools for biomonitoring, conservation and restoration. This study used a functional trait approach to characterise the responses of macrophyte assemblages to reach-scale disturbances (measured by lack of riparian shading, altered hydromorphology and eutrophication), from 97 wadeable stream sites in an agriculturally impacted region of New Zealand. To determine whether macrophyte assemblages differed due to disturbances, we examined multidimensional assemblage functional structure in relation to eleven functional traits and further related two functional diversity indices (entropy and originality) to disturbances. Macrophyte assemblages showed distinct patterns in response to disturbances, with riparian shading and hydromorphological conditions being the strongest variables shaping macrophyte functional structure. In the multidimensional space, most of the non-native species were associated with disturbed conditions. These species had traits allowing faster colonisation rates (higher number of reproductive organs and larger root-rhizome system) and superior competitive abilities for resources (tall and dense canopy, heterophylly and greater preferences for light and nitrogen). In addition, lack of riparian shading increased the abundance of functionally distinct species (i.e. entropy), and eutrophication resulted in the growth of functionally unique species (i.e. originality). We demonstrated that stream reach-scale habitat disturbances were associated to a dominance of more productive species, equating to a greater abundance of non-native species. This, can result in a displacement of native species, habitat alterations, and changes to higher trophic level assemblages. Our results suggests that reach-scale management efforts such as the conservation and restoration of riparian vegetation that provides substantial shading and hydromorphologically diverse in-stream habitat, would have beneficial direct and indirect effects on ecosystem functioning, and contribute to the mitigation of land-use impacts.

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 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.

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.

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 ( Tg), is nearly monoexponential and displays a weak temperature dependence. Below Tg, 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 ∼1011 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.

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.