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

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

Ultrafast XANES Monitors Femtosecond Sequential Structural Evolution in Photoexcited Coenzyme B12

Polarized X-ray absorption near-edge structure (XANES) at the Co K-edge and broadband UV-vis transient absorption are used to monitor the sequential evolution of the excited-state structure of coenzyme B₁₂ (adenosylcobalamin) over the first picosecond following excitation. The initial state is characterized by sub-100 fs sequential changes around the central cobalt. These are polarized first in the y-direction orthogonal to the transition dipole and 50 fs later in the x-direction along the transition dipole. Expansion of the axial bonds follows on a ca. 200 fs time scale as the molecule moves out of the Franck-Condon active region of the potential energy surface. On the same 200 fs time scale there are electronic changes that result in the loss of stimulated emission and the appearance of a strong absorption at 340 nm. These measurements provide a cobalt-centered movie of the excited molecule as it evolves to the local excited-state minimum.

Ultrafast Structural Changes Decomposed from Serial Crystallographic Data

Direct visualization of electronic and molecular events during biochemical reactions is essential to mechanistic insights. This Letter presents an in-depth analysis of the serial crystallographic data sets collected by Barends and Schlichting et al. ( Science 2015 , 350 , 445 ) that probe the ligand photodissociation in carbonmonoxy myoglobin. This analysis reveals electron density changes caused by the formation of high-spin 3d atomic orbitals of the heme iron upon photolysis and their dynamic behaviors within the first few picoseconds. The heme iron is found popping out of and recoiling back into the heme plane in succession. These findings provide long-awaited visual validations for previous works using ultrafast spectroscopy and molecular dynamics simulations. Electron density variations are also found largely in the solvent during the first period of a low-frequency oscillation. This work demonstrates the importance of the analytical methods in detecting and isolating weak, transient signals of electronic changes arising from chemical reactions.

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.

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.

Moving in the Right Direction: Protein Vibrations Steering Function

Nearly all protein functions require structural change, such as enzymes clamping onto substrates, and ion channels opening and closing. These motions are a target for possible new therapies; however, the control mechanisms are under debate. Calculations have indicated protein vibrations enable structural change. However, previous measurements found these vibrations only weakly depend on the functional state. By using the novel technique of anisotropic terahertz microscopy, we find that there is a dramatic change to the vibrational directionality with inhibitor binding to lysozyme, whereas the vibrational energy distribution, as measured by neutron inelastic scattering, is only slightly altered. The anisotropic terahertz measurements provide unique access to the directionality of the intramolecular vibrations, and immediately resolve the inconsistency between calculations and previous measurements, which were only sensitive to the energy distribution. The biological importance of the vibrational directions versus the energy distribution is revealed by our calculations comparing wild-type lysozyme with a higher catalytic rate double deletion mutant. The vibrational energy distribution is identical, but the more efficient mutant shows an obvious reorientation of motions. These results show that it is essential to characterize the directionality of motion to understand and control protein dynamics to optimize or inhibit function.

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.

Investigations of the low-frequency spectral density of cytochrome c upon equilibrium unfolding

The equilibrium unfolding process of ferric horse heart cytochrome c (cyt c), induced by guanidinium hydrochloride (GdHCl), was studied using UV-vis absorption spectroscopy, resonance Raman spectroscopy, and vibrational coherence spectroscopy (VCS). The unfolding process was successfully fit using a three-state model which included the fully folded (N) and unfolded (U) states, along with an intermediate (I) assigned to a Lys bound heme. The VCS spectra revealed for the first time several low-frequency heme modes that are sensitive to cyt c unfolding: γ(a) (~50 cm(-1)), γ(b) (~80 cm(-1)), γ(c) (~100 cm(-1)), and ν(s)(His-Fe-His) at 205 cm(-1). These out-of-plane modes have potential functional relevance and are activated by protein-induced heme distortions. The free energies for the N-I and the I-U transitions at pH 7.0 and 20 °C were found to be 4.6 kcal/M and 11.6 kcal/M, respectively. Imidazole was also introduced to replace the methionine ligand so the unfolding can be modeled as a two-state system. The intensity of the mode γ(b)~80 cm(-1) remains nearly constant during the unfolding process, while the amplitudes of the other low frequency modes track with spectral changes observed at higher frequency. This confirms that the heme deformation changes are coupled to the protein tertiary structural changes that take place upon unfolding. These studies also reveal that damping of the coherent oscillations depends sensitively on the coupling between heme and the surrounding water solvent.

Vibrational probes and determinants of the S = 0 ⇌ S = 2 spin crossover in five-coordinate [Fe(TPP)(CN)]-

The low-frequency vibrational characterization of the spin-crossover complex, five-coordinate cyano(tetraphenylporphyrinato)iron(II), [Fe(TPP)(CN)](-), is reported. Nuclear resonance vibrational spectroscopy has been used to measure all low-frequency vibrations involving iron at several temperatures; this yields vibrational spectra of both the low- (S = 0) and high-spin (S = 2) states. Multitemperature oriented single-crystal measurements facilitate assignments of the vibrational character of all modes and are consistent with the DFT-predicted spectra. The availability of the entire iron vibrational spectrum allows for the complete correlation of the modes between the two spin states. These data demonstrate that not only do the frequencies of the vibrations shift to lower values for the high-spin species as would be expected owing to the weaker bonds in the high-spin state, but also the mixing of iron modes with ligand modes changes substantially. Diagrams illustrating the changing character of the modes and their correlation are given. The reduced iron-ligand frequencies are the primary factor in the entropic stabilization of the high-spin state responsible for the spin crossover.

Heme-protein vibrational couplings in cytochrome c provide a dynamic link that connects the heme-iron and the protein surface

The active site of cytochrome c (Cyt c) consists of a heme covalently linked to a pentapeptide segment (Cys-X-X-Cys-His), which provides a link between the heme and the protein surface, where the redox partners of Cyt c bind. To elucidate the vibrational properties of heme c, nuclear resonance vibrational spectroscopy (NRVS) measurements were performed on (57)Fe-labeled ferric Hydrogenobacter thermophilus cytochrome c(552), including (13)C(8)-heme-, (13)C(5)(15)N-Met-, and (13)C(15)N-polypeptide (pp)-labeled samples, revealing heme-based vibrational modes in the 200- to 450-cm(-1) spectral region. Simulations of the NRVS spectra of H. thermophilus cytochrome c(552) allowed for a complete assignment of the Fe vibrational spectrum of the protein-bound heme, as well as the quantitative determination of the amount of mixing between local heme vibrations and pp modes from the Cys-X-X-Cys-His motif. These results provide the basis to propose that heme-pp vibrational dynamic couplings play a role in electron transfer (ET) by coupling vibrations of the heme directly to vibrations of the pp at the protein-protein interface. This could allow for the direct transduction of the thermal (vibrational) energy from the protein surface to the heme that is released on protein/protein complex formation, or it could modulate the heme vibrations in the protein/protein complex to minimize reorganization energy. Both mechanisms lower energy barriers for ET. Notably, the conformation of the distal Met side chain is fine-tuned in the protein to localize heme-pp mixed vibrations within the 250- to 400-cm(-1) spectral region. These findings point to a particular orientation of the distal Met that maximizes ET.

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.

Hydrogen sulfide and hemeproteins: knowledge and mysteries

Historically, hydrogen sulfide (H(2)S) has been regarded as a poisonous gas, with a wide spectrum of toxic effects. However, like ·NO and CO, H(2)S is now referred to as a signaling gas involved in numerous physiological processes. The list of reports highlighting the physiological effects of H(2)S is rapidly expanding and several drug candidates are now being developed. As with ·NO and CO, not a single H(2)S target responsible for all the biological effects has been found till now. Nevertheless, it has been suggested that H(2)S can bind to hemeproteins, inducing different responses that can mediate its effects. For instance, the interaction of H(2)S with cytochrome c oxidase has been associated with the activation of the ATP-sensitive potassium channels, regulating muscle relaxation. Inhibition of cytochrome c oxidase by H(2)S has also been related to inducing a hibernation-like state. Although H(2)S might induce these effects by interacting with hemeproteins, the mechanisms underlying these interactions are obscure. Therefore, in this review we discuss the current state of knowledge about the interaction of H(2)S with vertebrate and invertebrate hemeproteins and postulate a generalized mechanism. Our goal is to stimulate further research aimed at evaluating plausible mechanisms that explain H(2)S reactivity with hemeproteins.

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.