Spectroscopic and Electrochemical Characterization of Cytochrome P450st‐DDAB Films on a Plastic‐Formed Carbon Electrode

We have studied the characterization of thermophilic cytochrome P450 (P450st)‐didodecyldimethylammonium bromide (DDAB) films by using UV‐vis absorption, resonance Raman spectroscopy, and electrochemical methods. The observed Raman spectrum indicated near‐native conformation of the heme iron in DDAB film on the surface of a glass slide, while on the surface of a plastic‐formed carbon (PFC) electrode, the conformation of P450st‐DDAB was very similar to that of heme‐DDAB film, suggesting the release of heme from P450st in DDAB films on PFC electrodes. When NaBr was added as salt to the casting solution, the result of Raman spectrum indicated near‐native conformation of P450st in DDAB film even on the PFC electrode, but no redox potential of P450st which has near native structure was observed. This study suggests the essential experimental conditions when working with heme protein‐DDAB films as, in some cases, heme iron from proteins is released on the surface of the electrode.     

Reaction of cytochrome P450BM3 and peroxynitrite yields nitrosyl complex

Peroxynitrite has come into the spotlight in recent years. Its effects on proteins have been implicated in several diseases such as acute lung injury, rheumatoid arthritis, implant rejection, artherosclerosis, Parkinson's disease, and Alzheimer's disease. Peroxynitrite is thought to inactivate a variety of proteins including thiolate-ligated heme proteins such as cytochrome P450 2B1 and PGI2 synthase, through the nitration of tyrosine residues. In previous studies it was reported that thiolate-ligated heme enzymes react with peroxynitrite to form a ferryl intermediate. In an effort to spectroscopically characterize this species in P450BM3, we discovered that the peroxynitrite-generated intermediate is not an FeIVoxo, but rather an iron-nitrosyl [FeNO]6 complex. We present density functional calculations as well as M?ssbauer and stopped-flow spectroscopic characterizations of the peroxynitrite-generated intermediate in P450BM3.     

How do azoles inhibit cytochrome P450 enzymes? A density functional study

To examine how azole inhibitors interact with the heme active site of the cytochrome P450 enzymes, we have performed a series of density functional theory studies on azole binding. These are the first density functional studies on azole interactions with a heme center and give fundamental insight into how azoles inhibit the catalytic function of P450 enzymes. Since azoles come in many varieties, we tested three typical azole motifs representing a broad range of azole and azole-type inhibitors: methylimidazolate, methyltriazolate, and pyridine. These structural motifs represent typical azoles, such as econazole, fluconazole, and metyrapone. The calculations show that azole binding is a stepwise mechanism whereby first the water molecule from the resting state of P450 is released from the sixth binding site of the heme to create a pentacoordinated active site followed by coordination of the azole nitrogen to the heme iron. This process leads to the breaking of a hydrogen bond between the resting state water molecule and the approaching inhibitor molecule. Although, formally, the water molecule is released in the first step of the reaction mechanism and a pentacoordinated heme is created, this does not lead to an observed spin state crossing. Thus, we show that release of a water molecule from the resting state of P450 enzymes to create a pentacoordinated heme will lead to a doublet to quartet spin state crossing at an Fe-OH(2) distance of approximately 3.0 A, while the azole substitution process takes place at shorter distances. Azoles bind heme with significantly stronger binding energies than a water molecule, so that these inhibitors block the catalytic cycle of the enzyme and prevent oxygen binding and the catalysis of substrate oxidation. Perturbations within the active site (e.g., a polarized environment) have little effect on the relative energies of azole binding. Studies with an extra hydrogen-bonded ethanol molecule in the model, mimicking the active site of the CYP121 P450, show that the resting state and azole binding structures are close in energy, which may lead to chemical equilibrium between the two structures, as indeed observed with recent protein structural studies that have demonstrated two distinct azole binding mechanisms to P450 heme.     

Novel Conformational Transitions of Human Cytochrome P450 2C8 during Thermal and Acid-induced Unfolding

Transitions among various heme coordination/spin states, heme environments and protein conformations of human cytochrome P450 2C8 were investigated under different denaturing conditions by means of electronic absorption and circular dichroism spectroscopies. It is the first report of it's kind. Our results indicated that the thermal and acid‐induced denaturation could convert P450 2C8 to various P420 forms. In the thermal unfolding process, the ferric P420 thermal form emerged with weakened Fe‐S (thiolate) bond. An absorption band at ca. 425 nm of the ferrous P420 2C8 thermal form was observed, suggesting that the axial Cys435 was protonated or displaced by other ligand. Moreover, the new coordination bond was stabilized when the temperature was cooled down. When binding with CO, the ferrous P420 2C8 thermal form had the protonated thiol of Cys435 as the axial ligand. X‐ray structure of P450 2C8 suggested that the specific structure of the β‐bulge where the axial cysteine ligand located might be the reason of the formation of these P420 2C8 thermal forms. In the acid‐induced unfolding studies, we found that at pH 3.0 the heme could be irreversibly released from the heme pocket of ferric and ferrous P450 2C8. Interestingly, the released heme could form a new coordination bond with an unidentified ligand at the surface of partially unfolded protein when binding with CO at reduced state.     

Enhanced electron transfer and lauric acid hydroxylation by site-directed mutagenesis of CYP119

CYP119, a cytochrome P450 from a thermophilic organism for which a crystal structure is available, is shown here to hydroxylate lauric acid in a reaction supported by putidaredoxin and putidaredoxin reductase. This fatty acid hydroxylation activity is increased 15-fold by T214V and D77R mutations. The T214V mutation increases the rate by facilitating substrate binding and enhancing the associated spin state change, whereas the D77R mutation improves binding of the heterologous redox partner putidaredoxin to CYP119 and the rate of electron transfer from it to the heme group. A sequence alignment with P450(cam) can, therefore, be used to identify a part of the binding site for putidaredoxin on an unrelated P450 enzyme. This information can be used to engineer by mutagenesis an improved complementarity of the protein-protein interface that results in improved electron transfer from putidaredoxin to the P450 enzyme. As a result, the catalytic activity of the thermo- and barostable CYP119 has been incorporated into a catalytic system that hydroxylates fatty acids.     

Oxygen activation and electron transfer in flavocytochrome P450 BM3

In flavocytochrome P450 BM3, there is a conserved phenylalanine residue at position 393 (Phe393), close to Cys400, the thiolate ligand to the heme. Substitution of Phe393 by Ala, His, Tyr, and Trp has allowed us to modulate the reduction potential of the heme, while retaining the structural integrity of the enzyme's active site. Substrate binding triggers electron transfer in P450 BM3 by inducing a shift from a low- to high-spin ferric heme and a 140 mV increase in the heme reduction potential. Kinetic analysis of the mutants indicated that the spin-state shift alone accelerates the rate of heme reduction (the rate determining step for overall catalysis) by 200-fold and that the concomitant shift in reduction potential is only responsible for a modest 2-fold rate enhancement. The second step in the P450 catalytic cycle involves binding of dioxygen to the ferrous heme. The stabilities of the oxy-ferrous complexes in the mutant enzymes were also analyzed using stopped-flow kinetics. These were found to be surprisingly stable, decaying to superoxide and ferric heme at rates of 0.01-0.5 s(-)(1). The stability of the oxy-ferrous complexes was greater for mutants with higher reduction potentials, which had lower catalytic turnover rates but faster heme reduction rates. The catalytic rate-determining step of these enzymes can no longer be the initial heme reduction event but is likely to be either reduction of the stabilized oxy-ferrous complex, i.e., the second flavin to heme electron transfer or a subsequent protonation event. Modulating the reduction potential of P450 BM3 appears to tune the two steps in opposite directions; the potential of the wild-type enzyme appears to be optimized to maximize the overall rate of turnover. The dependence of the visible absorption spectrum of the oxy-ferrous complex on the heme reduction potential is also discussed.     

Lignin Biodegradation by a Cytochrome P450 Enzyme: A Computational Study into Syringol Activation by GcoA

A recently characterized cytochrome P450 isozyme GcoA activates lignin components through a selective O-demethylation or alternatively an acetal formation reaction. These are important reactions in biotechnology and, because lignin is readily available; it being the main component in plant cell walls. In this work we present a density functional theory study on a large active site model of GcoA to investigate syringol activation by an iron(IV)-oxo heme cation radical oxidant (Compound I) leading to hemiacetal and acetal products. Several substrate-binding positions were tested and full energy landscapes calculated. The study shows that substrate positioning determines the product distributions. Thus, with the phenol group pointing away from the heme, an O-demethylation is predicted, whereas an initial hydrogen-atom abstraction of the weak phenolic O-H group would trigger a pathway leading to ring-closure to form acetal products. Predictions on how to engineer P450 GcoA to get more selective product distributions are given.     

Electrochemistry of cytochrome P450 BM3 in sodium dodecyl sulfate films

Direct electrochemistry of the cytochrome P450 BM3 heme domain (BM3) was achieved by confining the protein within sodium dodecyl sulfate (SDS) films on the surface of basal-plane graphite (BPG) electrodes. Cyclic voltammetry revealed the heme FeIII/II redox couple at -330 mV (vs Ag/AgCl, pH 7.4). Up to 10 V/s, the peak current was linear with the scan rate, allowing us to treat the system as surface-confined within this regime. The standard heterogeneous rate constant determined at 10 V/s was estimated to be 10 s-1. Voltammograms obtained for the BM3-SDS-BPG system in the presence of dioxygen exhibited catalytic waves at the onset of FeIII reduction. The altered heme reduction potential of the BM3-SDS-graphite system indicates that SDS is likely bound in the enzyme active-site region. Compared to other P450-surfactant systems, we find redox potentials and electron-transfer rates that differ by approximately 100 mV and >10-fold, respectively, indicating that the nature of the surfactant environment has a significant effect on the observed heme redox properties.     

How Is a Metabolic Intermediate Formed in the Mechanism‐Based Inactivation of Cytochrome P450 by Using 1,1‐Dimethylhydrazine: Hydrogen Abstraction or Nitrogen Oxidation?

A precise understanding of the mechanism‐based inactivation of cytochrome P450 enzymes (P450s) at the quantum mechanical level should allow more reliable predictions of drug–drug interactions than those currently available. Hydrazines are among the molecules that act as mechanism‐based inactivators to terminate the function of P450s, which are essential heme enzymes responsible for drug metabolism in the human body. Despite its importance, the mechanism explaining how a metabolic intermediate (MI) is formed from hydrazine is not fully understood. We used density functional theory (DFT) calculations to compare four possible mechanisms underlying the reaction between 1,1‐dimethylhydrazine (or unsymmetrical dimethylhydrazine, UDMH) and the reactive compound I (Cpd I) intermediate of P450. Our DFT calculations provided a clear view on how an aminonitrene‐type MI is formed from UDMH. In the most favorable pathway, hydrogen is spontaneously abstracted from the N2 atom of UDMH by Cpd I, followed by a second hydrogen abstraction from the N2 atom by Cpd II. Nitrogen oxidation of nitrogen atoms and hydrogen abstraction from the C? H bond of the methyl group were found to be less favorable than the hydrogen abstraction from the N? H bond. We also found that the reaction of protonated UDMH with Cpd I is rather sluggish. The aminonitrene‐type MI binds to the ferric heme more strongly than a water molecule. This is consistent with the notion that the catalytic cycle of P450 is impeded when such an MI is produced through the P450‐catalyzed reaction.     

Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle

Molecular mechanics (MM) methods are computationally affordable tools for screening chemical libraries of novel compounds for sites of P450 metabolism. One challenge for MM methods has been the absence of a consistent and transferable set of parameters for the heme within the P450 active site. Experimental data indicate that mammalian P450 enzymes vary greatly in the size, architecture, and plasticity of their active sites. Thus, obtaining X-ray-based geometries for the development of accurate MM parameters for the major classes of hepatic P450 remains a daunting task. Our previous work with preliminary gas-phase quantum mechanics (QM)-derived atomic partial charges greatly improved the accuracy of docking studies of raloxifene to CYP3A4. We have therefore developed and tested a consistent set of transferable MM parameters based on gas-phase QM calculations of two model systems of the heme-a truncated (T-HM) and a full (F-HM) for four states of the P450 catalytic cycle. Our results indicate that the use of the atomic partial charges from the F-HM further improves the accuracy of docked predictions for raloxifene to CYP3A4. Different patterns for substrate docking are also observed depending on the choice of heme model and state. Newly parameterized heme models are tested in implicit and explicitly solvated MD simulations in the absence and presence of enzyme structures, for CYP3A4, and appear to be stable on the nanosecond simulation timescale. The new force field for the various heme states may aid the community for simulations of P450 enzymes and other heme-containing enzymes.     

Compound I in heme thiolate enzymes: a comparative QM/MM study

This study directly compares the active species of heme enzymes, so-called Compound I (Cpd I), across the heme-thiolate enzyme family. Thus, sixty-four different Cpd I structures are calculated by hybrid quantum mechanical/molecular mechanical (QM/MM) methods using four different cysteine-ligated heme enzymes (P450(cam), the mutant P450(cam)-L358P, CPO and NOS) with varying QM region sizes in two multiplicities each. The overall result is that these Cpd I species are similar to each other with regard to many characteristic features. Hence, using the more stable CPO Cpd I as a model for P450 Cpd I in experiments should be a reasonable approach. However, systematic differences were also observed, and it is shown that NOS stands out in most comparisons. By analyzing the electrical field generated by the enzyme on the QM region, one can see that (a) the protein exerts a large influence and modifies all the Cpd I species compared with the gas-phase situation and (b) in NOS this field is approximately planar to the heme plane, whereas it is approximately perpendicular in the other enzymes, explaining the deviating results on NOS. The calculations on the P450(cam) mutant L358P show that the effects of removing the hydrogen bond between the heme sulfur and L358 are small at the Cpd I stage. Finally, Mossbauer parameters are calculated for the different Cpd I species, enabling future comparisons with experiments. These results are discussed in the broader context of recent findings of Cpd I species that exhibit large variations in the electronic structure due to the presence of the substrate.     

Tuning the regio- and stereoselectivity of C-H activation in n-octanes by cytochrome P450 BM-3 with fluorine substituents: evidence for interactions between a C-F bond and aromatic π systems

We employed the water-soluble cytochrome P450 BM-3 to study the activity and regiospecificity of oxidation of fluorinated n-octanes. Three mutations, A74G, F87V, and L188Q, were introduced into P450 BM-3 to allow the system to undergo n-octane oxidation. In addition, the alanine at residue 328 was replaced with a phenylalanine to introduce an aromatic residue into the hydrophobic pocket to examine whether or not van der Waals interactions between a C-F substituent in the substrate and the polarizable π system of the phenylalanine may be used to steer the positioning of the substrate within the active-site pocket of the enzyme and control the regioselectivity and stereoselectivity of hydroxylation. Interestingly, not only was the regioselectivity controlled when the fluorine substituent was judiciously positioned in the substrate, but the electron input into the iron-heme group became tightly coupled to the formation of product, essentially without abortive side reactions. Remarkable enhancement of the coupling efficiency between electron input and product formation was observed for a range of fluorinated octanes in the enzyme even without the A328F mutation, presumably because of interactions of the C-F substituent with the π system of the porphyrin macrocycle within the active-site pocket. Evidently, tightening the protein domain containing the heme pocket tunes the distribution of accessible enzyme conformations and the associated protein dynamics that activate the iron porphyrin for substrate hydroxylation to allow the reactions mediated by the high-valent Fe(IV)=O to become kinetically more commensurate with electron transfer from the flavin adenine dinucleotide (FAD)/flavin mononucleotide (FMN) reductase. These observations lend compelling evidence to support significant van der Waals interactions between the CF(2) group and aromatic π systems within the heme pocket when the fluorinated octane substrate is bound.     

Laser photoexcitation of NAD(P)H induces reduction of P450 BM3 heme domain on the microsecond time scale

We demonstrate that photoexcitation of NAD(P)H at 355 nm using a Nd:YAG laser leads to rapid reduction of the heme domain of the Bacillus megaterium fatty acid hydroxylase flavocytochrome P450 BM3. An aqueous electron derived from photoexcited NAD(P)H is rapidly transferred to the heme domain, enabling the formation of a carbon monoxy complex of the ferrous P450 (FeII-CO) on the microsecond time scale. Using this approach we have determined the limiting rate constant (1770 s-1 for substrate-free heme domain) for formation of the FeII-CO complex. We find no dependence of the observed rate of FeII-CO complex formation on NAD(P)H concentration but demonstrate a hyperbolic dependence on carbon monoxide concentration. The apparent dissociation constant for the complex of carbon monoxide bound noncovalently to the ferric form of the BM3 heme domain (and with NADH as reductant) is 323 microM. Binding of a P450 substrate (N-palmitoylglycine) weakened the complex between carbon monoxide and the ferric BM3 heme domain (Kd increased to 1404 microM) but enhanced the rate of formation of the FeII-CO complex (3036 s-1 for substrate-free heme domain). This study demonstrates the applicability of NAD(P)H photoexcitation as a method for rapid electron delivery to P450 enzymes and provides a new route to probing the P450 catalytic cycle and its transient intermediates.     

Proton-transfer dynamics in the activation of cytochrome P450eryF

Molecular dynamics simulations are combined with quantum chemistry calculations of instantaneous proton-transfer energy profiles to investigate proton-transfer events in the transient pathway of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1), from the oxyferrous species to the catalytically active ferryl oxygen species (compound I). This reaction is one of the most fundamental unresolved aspects in the mechanism of oxidation that is common to all cytochrome P450s. We find that this process involves an ultrafast proton transfer from the crystallographic water molecule W519 to the distal oxygen bound to the heme group, and a subsequent proton-transfer event from W564 to W519. Both proton-transfer events are found to be endothermic in the oxyferrous state, suggesting that the oxyferrous reduction is mechanistically linked to the proton-transfer dynamics. These findings indicate that the hydrogen bond network, proximate to the O(2)-binding cleft, plays a crucial functional role in the enzymatic activation of P450s. Our results are consistent with the effect of mutations on the enzymatic efficacy.     

Does Hydrogen‐Bonding Donation to Manganese(IV)–Oxo and Iron(IV)–Oxo Oxidants Affect the Oxygen‐Atom Transfer Ability? A Computational Study

Iron(IV)–oxo intermediates are involved in oxidations catalyzed by heme and nonheme iron enzymes, including the cytochromes P450. At the distal site of the heme in P450 Compound I (Fe^IV–oxo bound to porphyrin radical), the oxo group is involved in several hydrogen‐bonding interactions with the protein, but their role in catalysis is currently unknown. In this work, we investigate the effects of hydrogen bonding on the reactivity of high‐valent metal–oxo moiety in a nonheme iron biomimetic model complex with trigonal bipyramidal symmetry that has three hydrogen‐bond donors directed toward a metal(IV)–oxo group. We show these interactions lower the oxidative power of the oxidant in reactions with dehydroanthracene and cyclohexadiene dramatically as they decrease the strength of the O? H bond (BDE_OH) in the resulting metal(III)–hydroxo complex. Furthermore, the distal hydrogen‐bonding effects cause stereochemical repulsions with the approaching substrate and force a sideways attack rather than a more favorable attack from the top. The calculations, therefore, give important new insights into distal hydrogen bonding, and show that in biomimetic, and, by extension, enzymatic systems, the hydrogen bond may be important for proton‐relay mechanisms involved in the formation of the metal–oxo intermediates, but the enzyme pays the price for this by reduced hydrogen atom abstraction ability of the intermediate. Indeed, in nonheme iron enzymes, where no proton relay takes place, there generally is no donating hydrogen bond to the iron(IV)–oxo moiety.     

Synthesis of a [2]rotaxane incorporating a "magic sulfur ring" by the thiol-ene click reaction

The mild and highly efficient thiol-ene click reaction has been used to construct a rotaxane incorporating dibenzo-24-crown-8 (DB24C8) and a dibenzylammonium-derived thread in high yield under the irradiation of UV light. A rotaxane containing a disulfide linkage in the macrocycle was also synthesized by the thiol-ene click reaction. It has been demonstrated that the formation of the [2]rotaxane with the disulfide bond in the macrocycle occurs by a mechanism that is different to the threading-followed-by-stoppering process. The successful construction of a rotaxane directly from its constituent components, the macrocycle containing a disulfide linkage and the dibenzylammonium hexafluorophosphate salt, suggests that the space within the macrocycle incorporating the disulfide linkage is smaller than the phenyl unit and a plausible reaction mechanism has been proposed as follows: A small amount of the initiator forms two radicals upon the absorption of UV irradiation; the radicals act as a "key" to "unlock" the disulfide bond in the macrocycle. The resulting crown ether like moiety in the macrocycle is clipped around the ammonium ion center in the dumb-bell-shaped compound. The [2]rotaxane is generated upon recombination of the disulfide linkage.     

Quantum mechanical/molecular mechanical study on the mechanisms of compound I formation in the catalytic cycle of chloroperoxidase: an overview on heme enzymes

The formation of Compound I (Cpd I), the active species of the enzyme chloroperoxidase (CPO), was studied using QM/MM calculation. Starting from the substrate complex with hydrogen peroxide, FeIII-HOOH, we examined two alternative mechanisms on the three lowest spin-state surfaces. The calculations showed that the preferred pathway involves heterolytic O-O cleavage that proceeds via the iron hydroperoxide species, i.e., Compound 0 (Cpd 0), on the doublet-state surface. This process is effectively concerted, with a barrier of 12.4 kcal/mol, and is catalyzed by protonation of the distal OH group of Cpd 0. By comparison, the path that involves a direct O-O cleavage from FeIII-HOOH is less favored. A proton coupled electron transfer (PCET) feature was found to play an important role in the mechanism nascent from Cpd 0. Initially, the O-O cleavage progresses in a homolytic sense, but as soon as the proton is transferred to the distal OH, it triggers an electron transfer from the heme-oxo moiety to form water and Cpd I. This study enables us to generalize the mechanisms of O-O activation, elucidated so far by QM/MM calculations, for other heme enzymes, e.g., cytochrome P450cam, horseradish peroxidase (HRP), nitric oxide synthase (NOS), and heme oxygenase (HO). Much like for CPO, in the cases of P450 and HRP, the PCET lowers the barrier below the purely homolytic cleavage alternative (in our case, the homolytic mechanism is calculated directly from FeIII-HOOH). By contrast, the absence of PCET in HO, along with the robust water cluster, prefers a homolytic cleavage mechanism.     

Cofactor‐Free Light‐Driven Whole‐Cell Cytochrome P450 Catalysis

Cytochromes P450 can catalyze various regioselective and stereospecific oxidation reactions of non‐functionalized hydrocarbons. Here, we have designed a novel light‐driven platform for cofactor‐free, whole‐cell P450 photo‐biocatalysis using eosin Y (EY) as a photosensitizer. EY can easily enter into the cytoplasm of Escherichia coli and bind specifically to the heme domain of P450. The catalytic turnover of P450 was mediated through the direct transfer of photoinduced electrons from the photosensitized EY to the P450 heme domain under visible light illumination. The photoactivation of the P450 catalytic cycle in the absence of cofactors and redox partners is successfully conducted using many bacterial P450s (variants of P450 BM3) and human P450s (CYPs 1A1, 1A2, 1B1, 2A6, 2E1, and 3A4) for the bioconversion of different substrates, including marketed drugs (simvastatin, lovastatin, and omeprazole) and a steroid (17β‐estradiol), to demonstrate the general applicability of the light‐driven, cofactor‐free system.     

Conformational landscape surfing induced by off-on pi-pi stacking in a porphyrin-quinone dyad

A covalently linked porphyrin-quinone dyad crystallizes with two orientations of the quinone, extended away from (off) and cofacial with the porphyrin macrocycle (on), which induce different conformations of the macrocycle and model the recently proposed structural effect of a nearby residue on the heme prosthetic group of a nitric oxide synthase.     

Hartree–Fock and density functional theory calculations for the reaction mechanism of nitric oxide reductase cytochrome P450nor from Fusarium oxysporum

A reaction mechanism of a nitric oxide reductase, cytochrome P450nor (P450nor) from Fusarium oxysporum, was clarified by using Density functional theory and Hartree–Fock calculations. In this reaction mechanism, molecular orbital (MO) analysis revealed that the NO ligand dissociates from the heme iron immediately after one-electron reduction by NADH, and MO energy analysis revealed that NADH acts as a one-electron reducer, not as a two-electron reducer, and that NADH has a pivotal role different from other one-electron reducers. The role of NADH is to act as a double one-electron donor (i.e. one-electron transfer occurring twice) and to combine with the NO⁻ molecule by charge recombination reaction. Our quantum chemical calculations indicated that all reactions occurring in the heme pocket are too fast to become rate-limiting. Therefore, the rate-limiting steps in the proposed reaction mechanism are the process of capturing NO and NADH into the heme pocket and the process of expelling a product generated in the heme pocket. Kinetics of these processes was discussed based on large-amplitude vibration, which helps capturing and expelling processes in a widely opened heme pocket of P450nor. The reaction mechanism proposed here well explains published experimental data.