Roles of the axial push effect in cytochrome P450cam studied with the site-directed mutagenesis at the heme proximal site

Enhancing ferryl accumulation in H2O2-dependent cytochrome P450s

A facile strategy is presented to enhance the accumulation of ferryl (iron(IV)-oxo) species in H2O2 dependent cytochrome P450s (CYPs) of the CYP152 family. We report the characterization of a highly chemoselective CYP decarboxylase from Staphylococcus aureus (OleTSA) that is soluble at high concentrations. Examination of OleTSA Compound I (CpdI) accumulation with a variety of fatty acid substrates reveals a dependence on resting spin-state equilibrium. Alteration of this equilibrium through targeted mutagenesis of the proximal pocket favors the high-spin form, and as a result, enhances Cpd-I accumulation to nearly stoichiometric yields.

Agar plate-based activity assay for easy and fast screening of recombinant Pichia pastoris expressing unspecific peroxygenases

Unspecific peroxygenases (UPOs) are promising biocatalysts that catalyze oxyfunctionalization reactions without the need for costly cofactors. Pichia pastoris (reclassified as Komagataella phaffii) is considered an attractive host for heterologous expression of UPOs. However, integration of UPO-expression cassettes into the genome via a single cross-over yields recombinant Pichia transformants with different UPO gene copy numbers resulting in different expression levels. Selection of the most productive Pichia transformants by a commonly used screening in liquid medium in 96-well plates is laborious and lasts up to 5 days. In this work, we developed a simple two-step agar plate-based assay to screen P. pastoris transformants for UPO activity with less effort, within shorter time, and without automated screening devices. After cell growth and protein expression on agar plates supplemented with methanol and 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS), an additional top agar layer supplemented with ABTS and peroxide is added. UPO activity is visualized within 15 min by formation of green zones around UPO-secreting P. pastoris transformants. The assay was validated with two UPOs, AbrUPO from Aspergillus brasiliensis and evolved PaDa-I from Agrocybe aegerita. The assay results were confirmed in a quantitative 96-deep well plate screening in liquid medium.

Bidirectional Charge Transfer at the Heme Iron in Reversible and Quasi-irreversible Cytochrome P450 Inhibition

The coordination bonding between inhibitor ligands and heme iron plays a critical role in disrupting the essential catalytic functions of cytochrome P450 enzymes (P450s). Despite its intrinsic importance and consequential implications for human health, our current understanding of coordination bonding in P450 inhibition remains limited. To address this knowledge gap, we conducted a systematic theoretical analysis of the complexes between a ferric or a ferrous heme model and representative inhibitor ligands. Specifically, we evaluated the charge-transfer (CT) effect within these complexes by employing a series of theoretical methods based on density functional theory (DFT). Through a comprehensive analysis, we unveiled the relative significance of ligand-to-heme forward CT in the ferric and ferrous complexes of reversible inhibitors. In contrast, backward CT dominates over forward CT in the ferrous heme complexes of quasi-irreversible inhibitors. Further analysis using the compact frontier orbital method underscores the elevated electron-accepting abilities of quasi-irreversible inhibitors for π backdonation, which greatly amplifies their binding affinity for the ferrous heme. This study sheds light on the intricate mechanisms underlying P450 inhibition and provides valuable insights for future inhibitor design and development.

Rejigging Electron and Proton Transfer to Transition between Dioxygenase, Monooxygenase, Peroxygenase, and Oxygen Reduction Activity: Insights from Bioinspired Constructs of Heme Enzymes

Nature has employed heme proteins to execute a diverse set of vital life processes. Years of research have been devoted to understanding the factors which bias these heme enzymes, with all having a heme cofactor, toward distinct catalytic activity. Among them, axial ligation, distal super structure, and substrate binding pockets are few very vividly recognized ones. Detailed mechanistic investigation of these heme enzymes suggested that several of these enzymes, while functionally divergent, use similar intermediates. Furthermore, the formation and decay of these intermediates depend on proton and electron transfer processes in the enzyme active site. Over the past decade, work in this group, using in situ surface enhanced resonance Raman spectroscopy of synthetic and biosynthetic analogues of heme enzymes, a general idea of how proton and electron transfer rates relate to the lifetime of different O2 derived intermediates has been developed. These findings suggest that the enzymatic activities of all these heme enzymes can be integrated into one general cycle which can be branched out to different catalytic pathways by regulating the lifetime and population of each of these intermediates. This regulation can further be achieved by tuning the electron and proton transfer steps. By strategically populating one of these intermediates during oxygen reduction, one can navigate through different catalytic processes to a desired direction by altering proton and electron transfer steps.

Active Site Hydrogen Bonding Induced in Cytochrome P450cam by Effector Putidaredoxin

Cytochrome P450s are diverse and powerful catalysts that can activate molecular oxygen to oxidize a wide variety of substrates. Catalysis relies on effective uptake of two electrons and two protons. For cytochrome P450cam, an archetypal member of the superfamily, the second electron must be supplied by the redox partner putidaredoxin (Pdx). Pdx also plays an effector role beyond electron transfer, but after decades the mechanism remains under investigation. We applied infrared spectroscopy to heme-ligated CN^- to examine the influence of Pdx binding. The results indicate that Pdx induces the population of a conformation wherein the CN^- ligand forms a strong hydrogen bond to a solvent water molecule, experimentally corroborating the formation of a proposed proton delivery network. Further, characterization of T252A P450cam implicates the side chain of Thr252 in regulating the population equilibrium of hydrogen-bonded states within the P450cam/Pdx complex, which could underlie its role in directing activated oxygen toward product formation and preventing reaction uncoupling through peroxide release.

Navigating the Unnatural Reaction Space: Directed Evolution of Heme Proteins for Selective Carbene and Nitrene Transfer

Despite the astonishing diversity of naturally occurring biocatalytic processes, enzymes do not catalyze many of the transformations favored by synthetic chemists. Either nature does not care about the specific products, or if she does, she has adopted a different synthetic strategy. In many cases, the appropriate reagents used by synthetic chemists are not readily accessible to biological systems. Here, we discuss our efforts to expand the catalytic repertoire of enzymes to encompass powerful reactions previously known only in small-molecule catalysis: formation and transfer of reactive carbene and nitrene intermediates leading to a broad range of products, including products with bonds not known in biology. In light of the structural similarity of iron carbene (Fe=C(R¹)(R²)) and iron nitrene (Fe=NR) to the iron oxo (Fe=O) intermediate involved in cytochrome P450-catalyzed oxidation, we have used synthetic carbene and nitrene precursors that biological systems have not encountered and repurposed P450s to catalyze reactions that are not known in the natural world. The resulting protein catalysts are fully genetically encoded and function in intact microbial cells or cell-free lysates, where their performance can be improved and optimized by directed evolution. By leveraging the catalytic promiscuity of P450 enzymes, we evolved a range of carbene and nitrene transferases exhibiting excellent activity toward these new-to-nature reactions. Since our initial report in 2012, a number of other heme proteins including myoglobins, protoglobins, and cytochromes c have also been found and engineered to promote unnatural carbene and nitrene transfer. Due to the altered active-site environments, these heme proteins often displayed complementary activities and selectivities to P450s.

Using wild-type and engineered heme proteins, we and others have described a range of selective carbene transfer reactions, including cyclopropanation, cyclopropenation, Si–H insertion, B–H insertion, and C–H insertion. Similarly, a variety of asymmetric nitrene transfer processes including aziridination, sulfide imidation, C–H amidation, and, most recently, C–H amination have been demonstrated. The scopes of these biocatalytic carbene and nitrene transfer reactions are often complementary to the state-of-the-art processes based on small-molecule transition-metal catalysts, making engineered biocatalysts a valuable addition to the synthetic chemist’s toolbox. Moreover, enabled by the exquisite regio- and stereocontrol imposed by the enzyme catalyst, this biocatalytic platform provides an exciting opportunity to address challenging problems in modern synthetic chemistry and selective catalysis, including ones that have eluded synthetic chemists for decades.

Conformational Change Induced by Putidaredoxin Binding to Ferrous CO-ligated Cytochrome P450cam Characterized by 2D IR Spectroscopy

The importance of conformational dynamics to protein function is now well-appreciated. An outstanding question is whether they are involved in the effector role played by putidaredoxin (Pdx) in its reduction of the O₂ complex of cytochrome P450cam (P450cam), an archetypical member of the cytochrome P450 superfamily. Recent studies have reported that binding of Pdx induces a conformational change from a closed to an open state of ferric P450cam, but a similar conformational change does not appear to occur for the ferrous, CO-ligated enzyme. To better understand the effector role of Pdx when binding the ferrous, CO-ligated P450cam, we applied 2D IR spectroscopy to compare the conformations and dynamics of the wild-type (wt) enzyme in the absence and presence of Pdx, as well as of L358P P450cam (L358P), which has served as a putative model for the Pdx complex. The CO vibrations of the Pdx complex and L358P report population of two conformational states in which the CO experiences distinct environments. The dynamics among the CO frequencies indicate that the energy landscape of substates within one conformation are reflective of the closed state of P450cam, and for the other conformation, differ from the free wt enzyme, but are equivalent between the Pdx complex and L358P. The two states co-populated by the Pdx complex are postulated to reflect a loosely bound encounter complex and a more tightly bound state, as is commonly observed for the dynamic complexes of redox partners. Significantly, this study shows that the binding of Pdx to ferrous, CO-ligated P450cam does perturb the conformational ensemble in a way that might underlie the effector role of Pdx.

Proximal Pocket Controls Alkene Oxidation Selectivity of Cytochrome P450 and Chloroperoxidase toward Small, Nonpolar Substrates

This paper examines the influence of the proximal pockets of cytochrome P450_CAM and chloroperoxidase (CPO) on the relative favorability of catalytic epoxidation and allylic hydroxylation of olefins, a type of alkene oxidation selectivity. The study employs quantum mechanical models of the active site to isolate the proximal pocket's influence on the barrier for the selectivity-determining step for each reaction, using cyclohexene and cis-β-methylstyrene as substrates. The proximal pocket is found to preference epoxidation by 2-5 kcal/mol, the largest value being for CPO, converting the active heme-thiolate moiety from being intrinsically hydroxylation-selective to being intrinsically epoxidation-selective. This theoretical study, the first to correctly predict these enzymes' preference for epoxidation of allylic substrates, strongly suggests that the proximal pocket is the key determinant of alkene oxidation selectivity. The selectivity for epoxidation can be rationalized in terms of the proximal pocket's modulation of the thiolate's electron "push" and consequent influence on the heme redox potential and the basicity of the trans ligand.

Enzymatic Activation of the Emerging Drug Resveratrol

The plant originated stilbene "resveratrol" (3,4',5-trans-trihydroxystilbene) is well known for its diverse health benefits including anti-tumor, anti-inflammatory, anti-microbial, and anti-oxidant properties. Besides a significant amount of reports on different aspects of its application as prodrug in the last 50 years, still, a strategy leading to the production of the active drug is missing. The aim of this work was to evaluate the enzymatic activation of prodrug resveratrol to the effective drug piceatannol, without engaging expensive cofactors. Five different heme proteins were analyzed for the transformation of resveratrol. Kinetic parameters of resveratrol transformation and analysis of the transformed products were conducted through HPLC and GC-MS. Effect of pH and organic solvent on the transformation process had also been evaluated. Among all tested heme proteins, only a variant of cytochrome P450_BM3 from Bacillus megaterium (CYP_BM3F87A) was found suitable for piceatannol production. The most suitable pH for the reaction conditions was 8.5, while organic solvents did not show any effect on transformation. For resveratrol transformation, the turnover rate (k _cat) was 21.7 (± 0.6) min^-1, the affinity constant (K _M) showed a value of 55.7 (± 16.7) μM for a catalytic efficiency (k _cat/K _M) of 389 min^-1 mM^-1. GC-MS analysis showed that the only product from resveratrol transformation by cytochrome P450_BM3 is the biologically active piceatannol. The enzymatic transformation of resveratrol, an emerging compound with medical interest, to active product piceatannol by a variant of cytochrome P450_BM3 in the absence of expensive NADPH cofactor is demonstrated. This enzymatic process is economically attractive and can be scaled up to cover the increasing medical demand for piceatannol.

Effect of Redox Partner Binding on Cytochrome P450 Conformational Dynamics

Previous crystal structures of cytochrome P450cam complexed with its redox partner, putidaredoxin (Pdx), shows that P450cam adopts the open conformation. It has been hypothesized that the Pdx-induced shift toward the open state frees the essential Asp251 from salt bridges with Arg186 and Lys178 so that Asp251 can participate in a proton relay network required for O₂ activation. This in part explains why P450cam has such a strict requirement for Pdx. One problem with this view is that looser substrate-protein interactions in the open state may not be compatible with the observed regio- and stereoselective hydroxylation. In the present study, molecular dynamics simulations show that Pdx binding favors a conformation that stabilizes the active site and decreases camphor mobility yet retains a partially open conformation compatible with the required proton relay network. The R186A mutant which frees Asp251 in the absence of Pdx retains good enzyme activity, and the crystal structure shows that product, 5-exo-hydroxycamphor, is bound. This indicates that rupture of the Asp251-Arg186 relaxes selectivity with respect to source of electrons and enables X-ray generated reducing equivalents to support substrate hydroxylation. These combined computational and experimental results are consistent with the proposed role of Pdx in assisting the release of Asp251 from ion pairs so that it can participate in proton-coupled electron transfer.

Identification of productive and futile encounters in an electron transfer protein complex

Well-defined, stereospecific states in protein complexes are often in exchange with an ensemble of more dynamic orientations: the encounter states. The structure of the stereospecific complex between cytochrome P450cam and putidaredoxin was solved recently by X-ray diffraction as well as paramagnetic NMR spectroscopy. Other than the stereospecific complex, the NMR data clearly show the presence of additional states in the complex in solution. In these encounter states, populated for a small percentage of the time, putidaredoxin assumes multiple orientations and samples a large part of the surface of cytochrome P450cam. To characterize the nature of the encounter states, an extensive paramagnetic NMR dataset has been analyzed using the Maximum Occurrence of Regions methodology. The analysis reveals the location and maximal spatial extent of the additional states needed to fully explain the NMR data. Under the assumption of sparsity of the size of the conformational ensemble, several minor states can be located quite precisely. The distribution of these minor states correlates with the electrostatic potential map around cytochrome P450cam. Whereas some minor states are on isolated positively charged patches, others are connected to the stereospecific site via positively charged paths. The existence of electrostatically favorable pathways between the stereospecific interaction site and the different minor states or lack thereof suggests a means to discriminate between productive and futile encounter states.

How the Proximal Pocket May Influence the Enantiospecificities of Chloroperoxidase-Catalyzed Epoxidations of Olefins

Chloroperoxidase-catalyzed enantiospecific epoxidations of olefins are of significant biotechnological interest. Typical enantiomeric excesses are in the range of 66%-97% and translate into free energy differences on the order of 1 kcal/mol. These differences are generally attributed to the effect of the distal pocket. In this paper, we show that the influence of the proximal pocket on the electron transfer mechanism in the rate-limiting event may be just as significant for a quantitatively accurate account of the experimentally-measured enantiospecificities.

Chloroperoxidase-Catalyzed Epoxidation of Cis-β-Methylstyrene: NH-S Hydrogen Bonds and Proximal Helix Dipole Change the Catalytic Mechanism and Significantly Lower the Reaction Barrier

Proximal hydrogen bonding of the axial sulfur with the backbone amides (NH-S) is a conserved feature of heme-thiolate enzymes such as chloroperoxidase (CPO) and cytochrome P450 (P450). In CPO, the effect of NH-S bonds is amplified by the dipole moment of the proximal helix. Our gas-phase DFT studies show that the proximal pocket effect significantly enhances CPO's reactivity toward the epoxidation of olefinic substrates. Comparison of models with and without proximal pocket residues shows that with them, the barrier for Cβ-O bond formation is lowered by about ∼4.6 kcal/mol, while Cα-O-Cβ ring closure becomes barrierless. The dipole moment of the proximal helix was estimated to contribute 1/3 of the decrease, while the rest is attributed to the effect of NH-S bonds. The decrease of the reaction barrier correlates with increased electron density transfer to residues of the proximal pocket. The effect is most pronounced on the doublet spin surface and involves a change in the electron-transfer mechanism. A full enzyme QMMM study on the doublet spin surface gives about the same barrier as the gas-phase DFT study. The free-energy barrier was estimated to be in agreement with the experimental results for the CPO-catalyzed epoxidation of styrene.

Dissecting structural and electronic effects in inducible nitric oxide synthase

Nitric oxide synthases (NOSs) are haem-thiolate enzymes that catalyse the conversion of L-arginine (L-Arg) into NO and citrulline. Inducible NOS (iNOS) is responsible for delivery of NO in response to stressors during inflammation. The catalytic performance of iNOS is proposed to rely mainly on the haem midpoint potential and the ability of the substrate L-Arg to provide a hydrogen bond for oxygen activation (O-O scission). We present a study of native iNOS compared with iNOS-mesohaem, and investigate the formation of a low-spin ferric haem-aquo or -hydroxo species (P) in iNOS mutant W188H substituted with mesohaem. iNOS-mesohaem and W188H-mesohaem were stable and dimeric, and presented substrate-binding affinities comparable to those of their native counterparts. Single turnover reactions catalysed by iNOSoxy with L-Arg (first reaction step) or N-hydroxy-L-arginine (second reaction step) showed that mesohaem substitution triggered higher rates of Fe(II)O₂ conversion and altered other key kinetic parameters. We elucidated the first crystal structure of a NOS substituted with mesohaem and found essentially identical features compared with the structure of iNOS carrying native haem. This facilitated the dissection of structural and electronic effects. Mesohaem substitution substantially reduced the build-up of species P in W188H iNOS during catalysis, thus increasing its proficiency towards NO synthesis. The marked structural similarities of iNOSoxy containing native haem or mesohaem indicate that the kinetic behaviour observed in mesohaem-substituted iNOS is most heavily influenced by electronic effects rather than structural alterations.

The crystal structure of the versatile cytochrome P450 enzyme CYP109B1 from Bacillus subtilis

The crystal structure of the versatile CYP109B1 enzyme from Bacillus subtilis has been solved at 1.8 Å resolution. This is the first structure of an enzyme from this CYP family, whose members are prevalent across diverse species of bacteria. In the crystal structure the enzyme has an open conformation with an access channel leading from the heme to the surface. The substrate-free structure reveals the location of the key residues in the active site that are responsible for binding the substrate in the correct orientation for regioselective oxidation. Importantly, there are significant differences among these residues in members of the CYP109 and closely related CYP106 families and these likely account for the variations in substrate binding and oxidation profiles observed with these enzymes. A whole-cell oxidation biosystem was developed, which contains CYP109B1 and a phthalate family oxygenase reductase (PFOR), from Pseudomonas putida KT24440, as the electron transfer partner. This electron transfer system is able to support CYP109B1 activity resulting in the regioselective hydroxylation of both α- and β-ionone in vivo and in vitro. The PFOR is therefore a versatile electron transfer partner that is able to support the activity of CYP enzymes from other bacterium. The crystal structure of CYP109B1 has a positively charged proximal face and this explains why it can interact with PFOR and adrenodoxin which are predominantly negatively charged around their [2Fe-2S] clusters.

Monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 enzymes

This review examines the monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 (CYP) enzymes in bacterial, archaeal and mammalian systems. CYP enzymes catalyze monooxygenation reactions by inserting one oxygen atom from O2 into an enormous number and variety of substrates. The catalytic versatility of CYP stems from its ability to functionalize unactivated carbon-hydrogen (C-H) bonds of substrates through monooxygenation. The oxidative prowess of CYP in catalyzing monooxygenation reactions is attributed primarily to a porphyrin π radical ferryl intermediate known as Compound I (CpdI) (Por•+FeIV=O), or its ferryl radical resonance form (FeIV-O•). CYP-mediated hydroxylations occur via a consensus H atom abstraction/oxygen rebound mechanism involving an initial abstraction by CpdI of a H atom from the substrate, generating a highly-reactive protonated Compound II (CpdII) intermediate (FeIV-OH) and a carbon-centered alkyl radical that rebounds onto the ferryl hydroxyl moiety to yield the hydroxylated substrate. CYP enzymes utilize hydroperoxides, peracids, perborate, percarbonate, periodate, chlorite, iodosobenzene and N-oxides as surrogate oxygen atom donors to oxygenate substrates via the shunt pathway in the absence of NAD(P)H/O2 and reduction-oxidation (redox) auxiliary proteins. It has been difficult to isolate the historically elusive CpdI intermediate in the native NAD(P)H/O2-supported monooxygenase pathway and to determine its precise electronic structure and kinetic and physicochemical properties because of its high reactivity, unstable nature (t½~2 ms) and short life cycle, prompting suggestions for participation in monooxygenation reactions of alternative CYP iron-oxygen intermediates such as the ferric-peroxo anion species (FeIII-OO-), ferric-hydroperoxo species (FeIII-OOH) and FeIII-(H2O2) complex.

Fixed distance photoinduced electron transfer between Fe and Zn porphyrins encapsulated within the Zn HKUST-1 metal organic framework

An attractive strategy for the development of photocatalytic metal organic framework (MOF) materials is to co-encapsulate a photoactive electron donor with a catalytic electron acceptor within the MOF.

A comprehensive test set of epoxidation rate constants for iron(iv)-oxo porphyrin cation radical complexes

Trends in oxygen atom transfer to Compound I of the P450 models with an extensive test set have been studied and show a preferred regioselectivity of epoxidation over hydroxylation in the gas-phase for the first time.

Cytochrome P450 enzymes are heme based monoxygenases that catalyse a range of oxygen atom transfer reactions with various substrates, including aliphatic and aromatic hydroxylation as well as epoxidation reactions. The active species is short-lived and difficult to trap and characterize experimentally, moreover, it reacts in a regioselective manner with substrates leading to aliphatic hydroxylation and epoxidation products, but the origin of this regioselectivity is poorly understood. We have synthesized a model complex and studied it with low-pressure Fourier transform-ion cyclotron resonance (FT-ICR) mass spectrometry (MS). A novel approach was devised using the reaction of [Fe^III(TPFPP)]⁺ (TPFPP = meso-tetrakis(pentafluorophenyl)porphinato dianion) with iodosylbenzene as a terminal oxidant which leads to the production of ions corresponding to [Fe^IV(O)(TPFPP⁺˙)]⁺. This species was isolated in the gas-phase and studied in its reactivity with a variety of olefins. Product patterns and rate constants under Ideal Gas conditions were determined by FT-ICR MS. All substrates react with [Fe^IV(O)(TPFPP⁺˙)]⁺ by a more or less efficient oxygen atom transfer process. In addition, substrates with low ionization energies react by a charge-transfer channel, which enabled us to determine the electron affinity of [Fe^IV(O)(TPFPP⁺˙)]⁺ for the first time. Interestingly, no hydrogen atom abstraction pathways are observed for the reaction of [Fe^IV(O)(TPFPP⁺˙)]⁺ with prototypical olefins such as propene, cyclohexene and cyclohexadiene and also no kinetic isotope effect in the reaction rate is found, which suggests that the competition between epoxidation and hydroxylation – in the gas-phase – is in favour of substrate epoxidation. This notion further implies that P450 enzymes will need to adapt their substrate binding pocket, in order to enable favourable aliphatic hydroxylation over double bond epoxidation pathways. The MS studies yield a large test-set of experimental reaction rates of iron(iv)–oxo porphyrin cation radical complexes, so far unprecedented in the gas-phase, providing a benchmark for calibration studies using computational techniques. Preliminary computational results presented here confirm the observed trends excellently and rationalize the reactivities within the framework of thermochemical considerations and valence bond schemes.

Trends in predicted chemoselectivity of cytochrome P450 oxidation: B3LYP barrier heights for epoxidation and hydroxylation reactions

Prediction of epoxide formation in drug metabolism is a difficult but important task, as epoxide formation is linked to drug toxicity. A comparison of the energy barriers for cytochrome P450 mediated epoxidation of alkenes to the barriers for the hydroxylation of an aliphatic carbon atom next to a double bond has been performed using B3LYP and B3LYP-D3. Relevant experimental data on oxidation selectivity has also been assessed. The results show that density functional theory, when using B3LYP-D3, does well in reproducing the experimental trends. Considering that the comparison involves chemical steps with quite different features this is remarkable. We also find that B3LYP consistently underestimates the hydrogen abstraction barriers relative to the epoxidation barriers, and that including a dispersion correction reduces this problem.

Electrochemistry of mammalian cytochrome P450 2B4 indicates tunable thermodynamic parameters in surfactant films

Electrochemical methods continue to present an attractive means for achieving in vitro biocatalysis with cytochromes P450; however understanding fully the nature of electrode-bound P450 remains elusive. Herein we report thermodynamic parameters using electrochemical analysis of full-length mammalian microsomal cytochrome P450 2B4 (CYP 2B4) in didodecyldimethylammonium bromide (DDAB) surfactant films. Electronic absorption spectra of CYP 2B4-DDAB films on silica slides reveal an absorption maximum at 418nm, characteristic of low-spin, six-coordinate, water-ligated Fe(III) heme in P450. The Fe(III/II) and Fe(II/I) redox couples (E1/2) of substrate-free CYP 2B4 measured by cyclic voltammetry are -0.23V and -1.02V (vs. SCE, or 14mV and -776mV vs. NHE) at 21°C. The standard heterogeneous rate constant for electron transfer from the electrode to the heme for the Fe(III/II) couple was estimated at 170s(-1). Experiments indicate that the system is capable of catalytic reduction of dioxygen, however substrate oxidation was not observed. From the variation of E1/2 with temperature (18-40°C), we have measured entropy and enthalpy changes that accompany heme reduction, -151Jmol(-1)K(-1) and -46kJmol(-1), respectfully. The corresponding entropy and enthalpy values are less for the six-coordinate low-spin, imidazole-ligated enzyme (-59Jmol(-1)K(-1) and -18kJmol(-1)), consistent with limited conformational changes upon reduction. These thermodynamic parameters are comparable to those measured for bacterial P450 from Bacillus megaterium (CYP BM3), confirming our prior reports that the surfactant environment exerts a strong influence on the redox properties of the heme.

Synergistic effects of mutations in cytochrome P450cam designed to mimic CYP101D1

A close orthologue to cytochrome P450cam (CYP101A1) that catalyzes the same hydroxylation of camphor to 5-exo-hydroxycamphor is CYP101D1. There are potentially important differences in and around the active site that could contribute to subtle functional differences. Adjacent to the heme iron ligand, Cys357, is Leu358 in P450cam, whereas this residue is Ala in CYP101D1. Leu358 plays a role in binding of the P450cam redox partner, putidaredoxin (Pdx). On the opposite side of the heme, about 15-20 Å away, Asp251 in P450cam plays a critical role in a proton relay network required for O2 activation but forms strong ion pairs with Arg186 and Lys178. In CYP101D1 Gly replaces Lys178. Thus, the local electrostatic environment and ion pairing are substantially different in CYP101D1. These sites have been systematically mutated in P450cam to the corresponding residues in CYP101D1 and the mutants analyzed by crystallography, kinetics, and UV-vis spectroscopy. Individually, the mutants have little effect on activity or structure, but in combination there is a major drop in enzyme activity. This loss in activity is due to the mutants being locked in the low-spin state, which prevents electron transfer from the P450cam redox partner, Pdx. These studies illustrate the strong synergistic effects on well-separated parts of the structure in controlling the equilibrium between the open (low-spin) and closed (high-spin) conformational states.

Catalytic intermediates of inducible nitric-oxide synthase stabilized by the W188H mutation

Nitric-oxide synthase (NOS) catalyzes nitric oxide (NO) synthesis via a two-step process: L-arginine (L-Arg) → N-hydroxy-L-arginine → citrulline + NO. In the active site the heme is coordinated by a thiolate ligand, which accepts a H-bond from a nearby tryptophan residue, Trp-188. Mutation of Trp-188 to histidine in murine inducible NOS was shown to retard NO synthesis and allow for transient accumulation of a new intermediate with a Soret maximum at 420 nm during the L-Arg hydroxylation reaction (Tejero, J., Biswas, A., Wang, Z. Q., Page, R. C., Haque, M. M., Hemann, C., Zweier, J. L., Misra, S., and Stuehr, D. J. (2008) J. Biol. Chem. 283, 33498-33507). However, crystallographic data showed that the mutation did not perturb the overall structure of the enzyme. To understand how the proximal mutation affects the oxygen chemistry, we carried out biophysical studies of the W188H mutant. Our stopped-flow data showed that the 420-nm intermediate was not only populated during the L-Arg reaction but also during the N-hydroxy-L-arginine reaction. Spectroscopic data and structural analysis demonstrated that the 420-nm intermediate is a hydroxide-bound ferric heme species that is stabilized by an out-of-plane distortion of the heme macrocycle and a cation radical centered on the tetrahydrobiopterin cofactor. The current data add important new insights into the previously proposed catalytic mechanism of NOS (Li, D., Kabir, M., Stuehr, D. J., Rousseau, D. L., and Yeh, S. R. (2007) J. Am. Chem. Soc. 129, 6943-6951).

NADPH-cytochrome P450 oxidoreductase: prototypic member of the diflavin reductase family

NADPH-cytochrome P450 oxidoreductase (CYPOR) and nitric oxide synthase (NOS), two members of the diflavin oxidoreductase family, are multi-domain enzymes containing distinct FAD and FMN domains connected by a flexible hinge. FAD accepts a hydride ion from NADPH, and reduced FAD donates electrons to FMN, which in turn transfers electrons to the heme center of cytochrome P450 or NOS oxygenase domain. Structural analysis of CYPOR, the prototype of this enzyme family, has revealed the exact nature of the domain arrangement and the role of residues involved in cofactor binding. Recent structural and biophysical studies of CYPOR have shown that the two flavin domains undergo large domain movements during catalysis. NOS isoforms contain additional regulatory elements within the reductase domain that control electron transfer through Ca(2+)-dependent calmodulin (CaM) binding. The recent crystal structure of an iNOS Ca(2+)/CaM-FMN construct, containing the FMN domain in complex with Ca(2+)/CaM, provided structural information on the linkage between the reductase and oxgenase domains of NOS, making it possible to model the holo iNOS structure. This review summarizes recent advances in our understanding of the dynamics of domain movements during CYPOR catalysis and the role of the NOS diflavin reductase domain in the regulation of NOS isozyme activities.

Effects of Dispersion in Density Functional Based Quantum Mechanical/Molecular Mechanical Calculations on Cytochrome P450 Catalyzed Reactions

Density functional theory (DFT) based quantum mechanical/molecular mechanical (QM/MM) calculations have provided valuable insight into the reactivity of the cytochrome P450 family of enzymes (P450s). A failure of commonly used DFT methods, such as B3LYP, is the neglect of dispersion interactions. An empirical dispersion correction has been shown to improve the accuracy of gas phase DFT calculations of P450s. The current work examines the effect of the dispersion correction in QM/MM calculations on P450s. The hydrogen abstraction from camphor, and hydrogen abstraction and C-O addition of cyclohexene and propene by P450cam have been modeled, along with the addition of benzene to Compound I in CYP2C9, at the B3LYP-D2/CHARMM27 level of theory. Single point energy calculations were also performed at the B3LYP-D3//B3LYP-D2/CHARMM27 level. The dispersion corrections lower activation energy barriers significantly (by ∼5 kcal/mol), as seen for gas phase calculations, but has a small effect on optimized geometries.These effects are likely to be important in modeling reactions catalyzed by other enzymes also. Given the low computational cost of including such dispersion corrections, we recommend doing so in all B3LYP based QM/MM calculations.

Crystal structures of substrate-free and nitrosyl cytochrome P450cin: implications for O(2) activation

The crystal structure of the P450cin substrate-bound nitric oxide complex and the substrate-free form have been determined revealing a substrate-free structure that adopts an open conformation relative to the substrate-bound structure. The region of the I helix that forms part of the O(2) binding pocket shifts from an α helix in the substrate-free form to a π helix in the substrate-bound form. Unique to P450cin is an active site residue, Asn242, in the I helix that H-bonds with the substrate. In most other P450s this residue is a Thr and plays an important role in O(2) activation by participating in an H-bonding network required for O(2) activation. The π/α I helix transition results in the carbonyl O atom of Gly238 moving in to form an H-bond with the water/hydroxide ligand in the substrate-free form. The corresponding residue, Gly248, in the substrate-free P450cam structure experiences a similar motion. Most significantly, in the oxy-P450cam complex Gly248 adopts a position midway between the substrate-free and -bound states. A comparison between these P450cam and the new P450cin structures provides insights into differences in how the two P450s activate O(2). The structure of P450cin complexed with nitric oxide, a close mimic of the O(2) complex, shows that Gly238 is likely to form tighter interactions with ligands than the corresponding Gly248 in P450cam. Having a close interaction between an H-bond acceptor, the Gly238 carbonyl O atom, and the distal oxygen atom of O(2) will promote protonation and hence further reduction of the oxy complex to the hydroperoxy intermediate resulting in heterolytic cleavage of the peroxide O-O bond and formation of the active ferryl intermediate required for substrate hydroxylation.

Effect of the disease-causing R266K mutation on the heme and PLP environments of human cystathionine β-synthase

Cystathionine β-synthase (CBS) is an essential pyridoxal 5'-phosphate (PLP)-dependent enzyme of the transsulfuration pathway that condenses serine with homocysteine to form cystathionine; intriguingly, human CBS also contains a heme b cofactor of unknown function. Herein we describe the enzymatic and spectroscopic properties of a disease-associated R266K hCBS variant, which has an altered hydrogen-bonding environment. The R266K hCBS contains a low-spin, six-coordinate Fe(III) heme bearing a His/Cys ligation motif, like that of WT hCBS; however, there is a geometric distortion that exists at the R266K heme. Using rR spectroscopy, we show that the Fe(III)-Cys(thiolate) bond is longer and weaker in R266K, as evidenced by an 8 cm(-1) downshift in the ν(Fe-S) resonance. Presence of this longer and weaker Fe(III)-Cys(thiolate) bond is correlated with alteration of the fluorescence spectrum of the active PLP ketoenamine tautomer. Activity data demonstrate that, relative to WT, the R266K variant is more impaired in the alternative cysteine-synthesis reaction than in the canonical cystathionine-synthesis reaction. This diminished cysteine synthesis activity and a greater sensitivity to exogenous PLP correlate with the change in PLP environment. Fe-S(Cys) bond weakening causes a nearly 300-fold increase in the rate of ligand switching upon reduction of the R266K heme. Combined, these data demonstrate cross talk between the heme and PLP active sites, consistent with previous proposals, revealing that alteration of the Arg(266)-Cys(52) interaction affects PLP-dependent activity and dramatically destabilizes the ferrous thiolate-ligated heme complex, underscoring the importance of this hydrogen-bonding residue pair.

A single-site mutation (F429H) converts the enzyme CYP 2B4 into a heme oxygenase: a QM/MM study

The intriguing deactivation of the cytochrome P450 (CYP) 2B4 enzyme induced by mutation of a single residue, Phe429 to His, is explored by quantum mechanical/molecular mechanical calculations of the O-OH bond activation of the (Fe(3+)OOH)(-) intermediate. It is found that the F429H mutant of CYP 2B4 undergoes homolytic instead of heterolytic O-OH bond cleavage. Thus, the mutant acquires the following characteristics of a heme oxygenase enzyme: (a) donation by His429 of an additional NH---S H-bond to the cysteine ligand combined with the presence of the substrate retards the heterolytic cleavage and gives rise to homolytic O-OH cleavage, and (b) the Thr302/water cluster orients nascent OH(•) and ensures efficient meso hydroxylation.

Influence of heme-thiolate in shaping the catalytic properties of a bacterial nitric-oxide synthase

Nitric-oxide synthases (NOS) are heme-thiolate enzymes that generate nitric oxide (NO) from L-arginine. Mammalian and bacterial NOSs contain a conserved tryptophan (Trp) that hydrogen bonds with the heme-thiolate ligand. We mutated Trp(66) to His and Phe (W66H, W66F) in B. subtilis NOS to investigate how heme-thiolate electronic properties control enzyme catalysis. The mutations had opposite effects on heme midpoint potential (-302, -361, and -427 mV for W66H, wild-type (WT), and W66F, respectively). These changes were associated with rank order (W66H < WT < W66F) changes in the rates of oxygen activation and product formation in Arg hydroxylation and N-hydroxyarginine (NOHA) oxidation single turnover reactions, and in the O(2) reactivity of the ferrous heme-NO product complex. However, enzyme ferrous heme-O(2) autoxidation showed an opposite rank order. Tetrahydrofolate supported NO synthesis by WT and the mutant NOS. All three proteins showed similar extents of product formation (L-Arg → NOHA or NOHA → citrulline) in single turnover studies, but the W66F mutant showed a 2.5 times lower activity when the reactions were supported by flavoproteins and NADPH. We conclude that Trp(66) controls several catalytic parameters by tuning the electron density of the heme-thiolate bond. A greater electron density (as in W66F) improves oxygen activation and reactivity toward substrate, but decreases heme-dioxy stability and lowers the driving force for heme reduction. In the WT enzyme the Trp(66) residue balances these opposing effects for optimal catalysis.

The conserved Trp-Cys hydrogen bond dampens the "push effect" of the heme cysteinate proximal ligand during the first catalytic cycle of nitric oxide synthase

Residues surrounding and interacting with the heme proximal ligand are important for efficient catalysis by heme proteins. The nitric oxide synthases (NOSs) are thiolate-coordinated enzymes that catalyze the hydroxylation of l-Arg in the first of the two catalytic cycles needed to synthesize nitric oxide. In NOSs, the indole NH group of a conserved tryptophan [W56 of the bacterial NOS-like protein from Staphylococcus aureus (saNOS)] forms a hydrogen bond with the heme proximal cysteinate ligand. The purpose of this study was to determine the impact of increasing (W56F and W56Y variants) or decreasing (W56H variant) the electron density of the proximal cysteinate ligand on molecular oxygen (O(2)) activation using saNOS as a model. We show that the removal of the indole NH···S(-) bond for W56F and W56Y caused an increase in the electron density of the cysteinate. This was probed by the decrease of the midpoint reduction potential (E(1/2)) along with weakened σ-bonding and strengthened π-backbonding with distal ligands (CO and O(2)). On the other hand, the W56H variant showed stronger Fe-OO and Fe-CO bonds (strengthened σ-bonding) along with an elevated E(1/2), which is consistent with the formation of a strong NH···S(-) hydrogen bond from H56. We also show here that changing the electron density of the proximal thiolate controls its "push effect"; whereas the rates of both O(2) activation and autoxidation of the Fe(II)O(2) complex increase with the stronger push effect created by removing the indole NH···S(-) hydrogen bond (W56F and W56Y variants), the W56H variant showed an increased stability of the complex against autoxidation and a slower rate of O(2) activation. These results are discussed with regard to the roles played by the conserved tryptophan-cysteinate interaction in the first catalytic cycle of NOS.

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.

Mesohaem substitution reveals how haem electronic properties can influence the kinetic and catalytic parameters of neuronal NO synthase

NOSs (NO synthases, EC 1.14.13.39) are haem-thiolate enzymes that catalyse a two-step oxidation of L-arginine to generate NO. The structural and electronic features that regulate their NO synthesis activity are incompletely understood. To investigate how haem electronics govern the catalytic properties of NOS, we utilized a bacterial haem transporter protein to overexpress a mesohaem-containing nNOS (neuronal NOS) and characterized the enzyme using a variety of techniques. Mesohaem-nNOS catalysed NO synthesis and retained a coupled NADPH consumption much like the wild-type enzyme. However, mesohaem-nNOS had a decreased rate of Fe(III) haem reduction and had increased rates for haem-dioxy transformation, Fe(III) haem-NO dissociation and Fe(II) haem-NO reaction with O2. These changes are largely related to the 48 mV decrease in haem midpoint potential that we measured for the bound mesohaem cofactor. Mesohaem nNOS displayed a significantly lower Vmax and KmO2 value for its NO synthesis activity compared with wild-type nNOS. Computer simulation showed that these altered catalytic behaviours of mesohaem-nNOS are consistent with the changes in the kinetic parameters. Taken together, the results of the present study reveal that several key kinetic parameters are sensitive to changes in haem electronics in nNOS, and show how these changes combine to alter its catalytic behaviour.

Molecular basis for the inability of an oxygen atom donor ligand to replace the natural sulfur donor heme axial ligand in cytochrome P450 catalysis. Spectroscopic characterization of the Cys436Ser CYP2B4 mutant

All cytochrome P450s (CYPs) contain a cysteinate heme iron proximal ligand that plays a crucial role in their mechanism of action. Conversion of the proximal Cys436 to Ser in NH(2)-truncated microsomal CYP2B4 (ΔCYP2B4) transforms the enzyme into a two-electron NADPH oxidase producing H(2)O(2) without monooxygenase activity [K.P. Vatsis, H.M. Peng, M.J. Coon, J. Inorg. Biochem. 91 (2002) 542-553]. To examine the effects of this ligation change on the heme iron spin-state and coordination structure of ΔC436S CYP2B4, the magnetic circular dichroism and electronic absorption spectra of several oxidation/ligation states of the variant have been measured and compared with those of structurally defined heme complexes. The spectra of the substrate-free ferric mutant are indicative of a high-spin five-coordinate structure ligated by anionic serinate. The spectroscopic properties of the dithionite-reduced (deoxyferrous) protein are those of a five-coordinate (high-spin) state, and it is concluded that the proximal ligand has been protonated to yield neutral serine (ROH-donor). Low-spin six-coordinate ferrous complexes of the mutant with neutral sixth ligands (NO, CO, and O(2)) examined are also likely ligated by neutral serine, as would be expected for ferric complexes with anionic sixth ligands such as the hydroperoxo-ferric catalytic intermediate. Ligation of the heme iron by neutral serine vs. deprotonated cysteine is likely the result of the large difference in their acidity. Thus, without the necessary proximal ligand push of the cysteinate, although the ΔC436S mutant can accept two electrons and two protons, it is unable to heterolytically cleave the O-O bond of the hydroperoxo-ferric species to generate Compound I and hydroxylate the substrate.

Effects of substrate, protein environment, and proximal ligand mutation on compound I and compound 0 of chloroperoxidase

This paper investigates the enzyme chloroperoxidase (CPO) by means of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. The effects of anionic substrate, protein environment, and proximal ligand mutation on the high-valent iron-oxo species, compound I (Cpd I), and the ferric hydroperoxide complex, compound 0 (Cpd 0), are studied. The results indicate that the presence of an anionic substrate (acetate) and the protonation state of one critical residue (Glu104) have a considerable impact on the relative stabilities of Cpd I and Cpd 0. In the absence of the substrate or when the substrate is protonated, Cpd I is considerably more stable, and its formation barrier is smaller than in the case where the substrate is in its anionic state and when Glu104 is deprotonated. This trend, which is shown to be a simple manifestation of the Hammond principle, reproduces the experimental observation that the working pH of the enzyme is acidic. Furthermore, in the absence of substrate (or when it is protonated), the relative Cpd 0/Cpd I energies are found to be a good index of Cpd I stability in heme enzymes and to follow the experimental order: horseradish peroxidase (HRP) > CPO > P450. In silico mutation of the proximal ligand from cysteine to selenocysteine was found to have no effect at all on the properties of Cpd I (e.g., spin density on the chalcogen, Mössbauer parameters, etc.) and its relative stability to Cpd 0 or on the corresponding barrier for formation. This surprising finding shows that the polar CPO pocket applies a leveling effect that stabilizes the anionic forms of the proximal ligands (CysS(-) and CysSe(-)). This in turn means that the Se-Cpd I of the mutant CPO is observable.

Heme-coordinating inhibitors of neuronal nitric oxide synthase. Iron-thioether coordination is stabilized by hydrophobic contacts without increased inhibitor potency

The heme-thioether ligand interaction often occurs between heme iron and native methionine ligands, but thioether-based heme-coordinating (type II) inhibitors are uncommon due to the difficulty in stabilizing the Fe-S bond. Here, a thioether-based inhibitor (3) of neuronal nitric oxide synthase (nNOS) was designed, and its binding was characterized by spectrophotometry and crystallography. A crystal structure of inhibitor 3 coordinated to heme iron was obtained, representing, to our knowledge, the first crystal structure of a thioether inhibitor complexed to any heme enzyme. A series of related potential inhibitors (4-8) also were evaluated. Compounds 4-8 were all found to be type I (non-heme-coordinating) inhibitors of ferric nNOS, but 4 and 6-8 were found to switch to type II upon heme reduction to the ferrous state, reflecting the higher affinity of thioethers for ferrous heme than for ferric heme. Contrary to what has been widely thought, thioether-heme ligation was found not to increase inhibitor potency, illustrating the intrinsic weakness of the thioether-ferric heme linkage. Subtle changes in the alkyl groups attached to the thioether sulfur caused drastic changes in the binding conformation, indicating that hydrophobic contacts play a crucial role in stabilizing the thioether-heme coordination.

Probing the role of the proximal heme ligand in cytochrome P450cam by recombinant incorporation of selenocysteine

The unique monooxygenase activity of cytochrome P450cam has been attributed to coordination of a cysteine thiolate to the heme cofactor. To investigate this interaction, we replaced cysteine with the more electron-donating selenocysteine. Good yields of the selenoenzyme were obtained by bacterial expression of an engineered gene containing the requisite UGA codon for selenocysteine and a simplified yet functional selenocysteine insertion sequence (SECIS). The sulfur-to-selenium substitution subtly modulates the structural, electronic, and catalytic properties of the enzyme. Catalytic activity decreases only 2-fold, whereas substrate oxidation becomes partially uncoupled from electron transfer, implying a more complex role for the axial ligand than generally assumed.