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

Conformational heterogeneity suggests multiple substrate binding modes in CYP106A2

CYP106A2 (cytochrome P450meg) is a bacterial enzyme originally isolated from B. megaterium, and has been shown to hydroxylate a wide variety of substrates, including steroids. The regio- and stereochemistry of CYP106A2 hydroxylation has been shown to be dependent on a variety of factors, and hydroxylation often occurs at more than one site and/or with lack of stereospecificity for some substrates. Comprehensive backbone 15N, 1H and 13C resonance assignments based on multidimensional nuclear magnetic resonance (NMR) experiments performed with uniform and selective isotopically labeled CYP106A2 samples are reported herein, and broadening and splitting of resonances assigned to regions of the enzyme shown to be affected by substrate binding in other P450 enzymes indicate that substrate binding does not reduce structural heterogeneity as has been observed previously in P450 enzymes CYP101A1 and MycG. Paramagnetic relaxation enhancement (PRE) due to proximity between substrate protons and the heme iron were measured for three different substrates, and the relatively uniform nature of the PREs support the proposal that multiple substrate binding modes are occupied at saturating substrate concentrations.

Investigating the Active Oxidants Involved in Cytochrome P450 Catalyzed Sulfoxidation Reactions

Cytochrome P450 (CYP) heme-thiolate monooxygenases catalyze the hydroxylation of the C-H bonds of organic molecules. This reaction is initiated by a ferryl-oxo heme radical cation (Cpd I). These enzymes can also catalyze sulfoxidation reactions and the ferric-hydroperoxy complex (Cpd 0) and the Fe(III)-H₂ O₂ complex have been proposed as alternative oxidants for this transformation. To investigate this, the oxidation of 4-alkylthiobenzoic acids and 4-methoxybenzoic acid by the CYP199A4 enzyme from Rhodopseudomonas palustris HaA2 was compared using both monooxygenase and peroxygenase pathways. By examining mutants at the mechanistically important, conserved acid alcohol-pair (D251N, T252A and T252E) the relative amounts of the reactive intermediates that would form in these reactions were disturbed. Substrate binding and X-ray crystal structures helped to understand changes in the activity and enabled an attempt to evaluate whether multiple oxidants can participate in these reactions. In peroxygenase reactions the T252E mutant had higher activity towards sulfoxidation than O-demethylation but in the monooxygenase reactions with the WT enzyme the activity of both reactions was similar. The peroxygenase activity of the T252A mutant was greater for sulfoxidation reactions than the WT enzyme, which is the reverse of the activity changes observed for O-demethylation. The monooxygenase activity and coupling efficiency of sulfoxidation and oxidative demethylation were reduced by similar degrees with the T252A mutant. These observations infer that while Cpd I is required for O-dealkylation, another oxidant may contribute to sulfoxidation. Based on the activity of the CYP199A4 mutants it is proposed that this is the Fe(III)-H₂ O₂ complex which would be more abundant in the peroxide-driven reactions.

Hydroxylation Regiochemistry Is Robust to Active Site Mutations in Cytochrome P450cam (CYP101A1)

Cytochrome P450_cam (CYP101A1) catalyzes the hydroxylation of d-camphor by molecular oxygen. The enzyme-catalyzed hydroxylation exhibits a high degree of regioselectivity and stereoselectivity, with a single major product, d-5-exo-hydroxycamphor, suggesting that the substrate is oriented to facilitate this specificity. In previous work, we used an elastic network model and perturbation response scanning to show that normal deformation modes of the enzyme structure are highly responsive not only to the presence of a substrate but also to the substrate orientation. This work examines the effects of mutations near the active site on substrate localization and orientation. The investigated mutations were designed to promote a change in substrate orientation and/or location that might give rise to different hydroxylation products, while maintaining the same carbon and oxygen atom balances as in the wild type (WT) enzyme. Computational experiments and parallel in vitro site-directed mutations of CYP101A1 were used to examine reaction products and enzyme activity. ¹H-¹⁵N TROSY-HSQC correlation maps were used to compare the computational results with detectable perturbations in the enzyme structure and dynamics. We found that all of the mutant enzymes retained the same regio- and stereospecificity of hydroxylation as the WT enzyme, with varying degrees of efficiency, which suggests that large portions of the enzyme have been subjected to evolutionary pressure to arrive at the appropriate sequence-structure combination for efficient 5-exo hydroxylation of camphor.

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.

High-Precision Megahertz-to-Terahertz Dielectric Spectroscopy of Protein Collective Motions and Hydration Dynamics

The low-frequency collective vibrational modes in proteins as well as the protein-water interface have been suggested as dominant factors controlling the efficiency of biochemical reactions and biological energy transport. It is thus crucial to uncover the mystery of the hydration structure and dynamics as well as their coupling to collective motions of proteins in aqueous solutions. Here, we report dielectric properties of aqueous bovine serum albumin protein solutions as a model system using an extremely sensitive dielectric spectrometer with frequencies spanning from megahertz to terahertz. The dielectric relaxation spectra reveal several polarization mechanisms at the molecular level with different time constants and dielectric strengths, reflecting the complexity of protein-water interactions. Combining the effective-medium approximation and molecular dynamics simulations, we have determined collective vibrational modes at terahertz frequencies and the number of water molecules in the tightly bound and loosely bound hydration layers. High-precision measurements of the number of hydration water molecules indicate that the dynamical influence of proteins extends beyond the first solvation layer, to around 7 Å distance from the protein surface, with the largest slowdown arising from water molecules directly hydrogen-bonded to the protein. Our results reveal critical information of protein dynamics and protein-water interfaces, which determine biochemical functions and reactivity of proteins.

Spectroscopic studies of the cytochrome P450 reaction mechanisms

The cytochrome P450 monooxygenases (P450s) are thiolate heme proteins that can, often under physiological conditions, catalyze many distinct oxidative transformations on a wide variety of molecules, including relatively simple alkanes or fatty acids, as well as more complex compounds such as steroids and exogenous pollutants. They perform such impressive chemistry utilizing a sophisticated catalytic cycle that involves a series of consecutive chemical transformations of heme prosthetic group. Each of these steps provides a unique spectral signature that reflects changes in oxidation or spin states, deformation of the porphyrin ring or alteration of dioxygen moieties. For a long time, the focus of cytochrome P450 research was to understand the underlying reaction mechanism of each enzymatic step, with the biggest challenge being identification and characterization of the powerful oxidizing intermediates. Spectroscopic methods, such as electronic absorption (UV-Vis), electron paramagnetic resonance (EPR), nuclear magnetic resonance (NMR), electron nuclear double resonance (ENDOR), Mössbauer, X-ray absorption (XAS), and resonance Raman (rR), have been useful tools in providing multifaceted and detailed mechanistic insights into the biophysics and biochemistry of these fascinating enzymes. The combination of spectroscopic techniques with novel approaches, such as cryoreduction and Nanodisc technology, allowed for generation, trapping and characterizing long sought transient intermediates, a task that has been difficult to achieve using other methods. Results obtained from the UV-Vis, rR and EPR spectroscopies are the main focus of this review, while the remaining spectroscopic techniques are briefly summarized. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

NADH reduction of nitroaromatics as a probe for residual ferric form high-spin in a cytochrome P450

The existence of a substrate-sensitive equilibrium between high spin (S=5/2) and low spin (S=1/2) ferric iron is a well-established phenomenon in the cytochrome P450 (CYP) superfamily, although its origins are still a subject of discussion. A series of mutations that strongly perturb the spin state equilibrium in the camphor hydroxylase CYP101A1 were recently described (Colthart et al., Sci. Rep. 6, 22035 (2016)). Wild type CYP101A1 as well as some CYP101A1 mutants are herein shown to be capable of catalyzing the reduction of nitroacetophenones by NADH to the corresponding anilino compounds (nitroreductase or NRase activity). The distinguishing characteristic between those mutants that catalyze the reduction and those that cannot appears to be the extent to which residual high spin form exists in the absence of the native substrate d-camphor, with those showing the largest spin state shifts upon camphor binding also exhibiting NRase activity. Optical and EPR spectroscopy was used to further examine these phenomena. These results suggest that reduction of nitroaromatics may provide a useful probe of residual high spin states in the CYP superfamily. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

A Comparative Analysis of the Effector Role of Redox Partner Binding in Bacterial P450s

The camphor monooxygenase, cytochrome P450cam, exhibits a strict requirement for its own redox partner, putidaredoxin (Pdx), a two-iron-sulfur ferredoxin. The closest homologue to P450cam, CYP101D1, is structurally very similar, uses a similar redox partner, and exhibits nearly identical enzymatic properties in the monooxygenation of camphor to give the same single 5-exo-hydroxy camphor product. However, CYP101D1 does not strictly require its own ferredoxin (Arx) for activity because Pdx can support CYP101D1 catalysis but Arx cannot support P450cam catalysis. We have further examined the differences between these two P450s by determining the effect of spin equilibrium, redox properties, and stability of oxygen complexes. We find that Arx shifts the spin state equilibrium toward high-spin, which is the opposite of the effect of Pdx on P450cam. In both P450s, redox partner binding destabilizes the oxy-P450 complex but this effect is much weaker with CYP101D1. In addition, resonance Raman data show that structural perturbations observed in P450cam upon addition of Pdx are absent in CYP101D1. These data indicate that Arx does not play the same effector role in catalysis as Pdx does with P450cam. The most relevant structural difference between these two P450s centers on a catalytically important Asp residue required for proton-coupled electron transfer. We postulate that with P450cam larger Pdx-assisted motions are required to free this Asp for catalysis while the smaller number of restrictions in CYP101D1 precludes the need for redox partner-assisted structural changes.

Detection of substrate-dependent conformational changes in the P450 fold by nuclear magnetic resonance

Cytochrome P450 monooxygenases typically catalyze the insertion of one atom of oxygen from O2 into unactivated carbon-hydrogen and carbon-carbon bonds, with concomitant reduction of the other oxygen atom to H2O by NAD(P)H. Comparison of the average structures of the camphor hydroxylase cytochrome P450(cam) (CYP101) obtained from residual dipolar coupling (RDC)-restrained molecular dynamics (MD) in the presence and absence of substrate camphor shows structural displacements resulting from the essential collapse of the active site upon substrate removal. This collapse has conformational consequences that extend across the protein structure, none of which were observed in analogous crystallographic structures. Mutations were made to test the involvement of the observed conformational changes in substrate binding and recognition. All of the mutations performed based upon the NMR-detected perturbations, even those remote from the active site, resulted in modified substrate selectivity, enzyme efficiency and/or haem iron spin state. The results demonstrate that solution NMR can provide insights into enzyme structure-function relationships that are difficult to obtain by other methods.

Coherent Nuclear Wave Packets in Q States by Ultrafast Internal Conversions in Free Base Tetraphenylporphyrin

Persistence of vibrational coherence in electronic transition has been noted especially in biochemical systems. Here, we report the dynamics between electronic excited states in free base tetraphenylporphyrin (H2TPP) by time-resolved fluorescence with high time resolution. Following the photoexcitation of the B state, ultrafast internal conversion occurs to the Qx state directly as well as via the Qy state. Unique and distinct coherent nuclear wave packet motions in the Qx and Qy states are observed through the modulation of the fluorescence intensity in time. The instant, serial internal conversions from the B to the Qy and Qx states generate the coherent wave packets. Theory and experiment show that the observed vibrational modes involve the out-of-plane vibrations of the porphyrin ring that are strongly coupled to the internal conversion of H2TPP.

Low frequency dynamics of the nitrogenase MoFe protein via femtosecond pump probe spectroscopy - Observation of a candidate promoting vibration

We have used femtosecond pump-probe spectroscopy (FPPS) to study the FeMo-cofactor within the nitrogenase (N2ase) MoFe protein from Azotobacter vinelandii. A sub-20-fs visible laser pulse was used to pump the sample to an excited electronic state, and a second sub-10-fs pulse was used to probe changes in transmission as a function of probe wavelength and delay time. The excited protein relaxes to the ground state with a ~1.2ps time constant. With the short laser pulse we coherently excited the vibrational modes associated with the FeMo-cofactor active site, which are then observed in the time domain. Superimposed on the relaxation dynamics, we distinguished a variety of oscillation frequencies with the strongest band peaks at ~84, 116, 189, and 226cm(-1). Comparison with data from nuclear resonance vibrational spectroscopy (NRVS) shows that the latter pair of signals comes predominantly from the FeMo-cofactor. The frequencies obtained from the FPPS experiment were interpreted with normal mode calculations using both an empirical force field (EFF) and density functional theory (DFT). The FPPS data were also compared with the first reported resonance Raman (RR) spectrum of the N2ase MoFe protein. This approach allows us to outline and assign vibrational modes having relevance to the catalytic activity of N2ase. In particular, the 226cm(-1) band is assigned as a potential 'promoting vibration' in the H-atom transfer (or proton-coupled electron transfer) processes that are an essential feature of N2ase catalysis. The results demonstrate that high-quality room-temperature solution data can be obtained on the MoFe protein by the FPPS technique and that these data provide added insight to the motions and possible operation of this protein and its catalytic prosthetic group.

Ultrafast infrared spectroscopy reveals water-mediated coherent dynamics in an enzyme active site

Ultrafast infrared spectroscopy provides insights into the dynamic nature of water in the active sites of catalase and peroxidase enzymes.

Understanding the impact of fast dynamics upon the chemical processes occurring within the active sites of proteins and enzymes is a key challenge that continues to attract significant interest, though direct experimental insight in the solution phase remains sparse. Similar gaps in our knowledge exist in understanding the role played by water, either as a solvent or as a structural/dynamic component of the active site. In order to investigate further the potential biological roles of water, we have employed ultrafast multidimensional infrared spectroscopy experiments that directly probe the structural and vibrational dynamics of NO bound to the ferric haem of the catalase enzyme from Corynebacterium glutamicum in both H₂O and D₂O. Despite catalases having what is believed to be a solvent-inaccessible active site, an isotopic dependence of the spectral diffusion and vibrational lifetime parameters of the NO stretching vibration are observed, indicating that water molecules interact directly with the haem ligand. Furthermore, IR pump–probe data feature oscillations originating from the preparation of a coherent superposition of low-frequency vibrational modes in the active site of catalase that are coupled to the haem ligand stretching vibration. Comparisons with an exemplar of the closely-related peroxidase enzyme family shows that they too exhibit solvent-dependent active-site dynamics, supporting the presence of interactions between the haem ligand and water molecules in the active sites of both catalases and peroxidases that may be linked to proton transfer events leading to the formation of the ferryl intermediate Compound I. In addition, a strong, water-mediated, hydrogen bonding structure is suggested to occur in catalase that is not replicated in peroxidase; an observation that may shed light on the origins of the different functions of the two enzymes.

Terahertz underdamped vibrational motion governs protein-ligand binding in solution

Low-frequency collective vibrational modes in proteins have been proposed as being responsible for efficiently directing biochemical reactions and biological energy transport. However, evidence of the existence of delocalized vibrational modes is scarce and proof of their involvement in biological function absent. Here we apply extremely sensitive femtosecond optical Kerr-effect spectroscopy to study the depolarized Raman spectra of lysozyme and its complex with the inhibitor triacetylchitotriose in solution. Underdamped delocalized vibrational modes in the terahertz frequency domain are identified and shown to blue-shift and strengthen upon inhibitor binding. This demonstrates that the ligand-binding coordinate in proteins is underdamped and not simply solvent-controlled as previously assumed. The presence of such underdamped delocalized modes in proteins may have significant implications for the understanding of the efficiency of ligand binding and protein-molecule interactions, and has wider implications for biochemical reactivity and biological 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.

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.