Hydroxylation by the hydroperoxy-iron species in cytochrome P450 enzymes

Intramolecular and intermolecular kinetic isotope effects (KIEs) were determined for hydroxylation of the enantiomers of trans-2-(p-trifluoromethylphenyl)cyclopropylmethane (1) by hepatic cytochrome P450 enzymes, P450s 2B1, Delta2B4, Delta2B4 T302A, Delta2E1, and Delta2E1 T303A. Two products from oxidation of the methyl group were obtained, unrearranged trans-2-(p-trifluoromethylphenyl)cyclopropylmethanol (2) and rearranged 1-(p-trifluoromethylphenyl)but-3-en-1-ol (3). In intramolecular KIE studies with dideuteriomethyl substrates (1-d(2)) and in intermolecular KIE studies with mixtures of undeuterated (1-d(0)) and trideuteriomethyl (1-d(3)) substrates, the apparent KIE for product 2 was consistently larger than the apparent KIE for product 3 by a factor of ca. 1.2. Large intramolecular KIEs found with 1-d(2) (k(H)/k(D) = 9-11 at 10 degrees C) were shown not to be complicated by tunneling effects by variable temperature studies with two P450 enzymes. The results require two independent isotope-sensitive processes in the overall hydroxylation reactions that are either competitive or sequential. Intermolecular KIEs were partially masked in all cases and largely masked for some P450s. The intra- and intermolecular KIE results were combined to determine the relative rate constants for the unmasking and hydroxylation reactions, and a qualitative correlation was found for the unmasking reaction and release of hydrogen peroxide from four of the P450 enzymes in the absence of substrate. The results are consistent with the two-oxidants model for P450 (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555), which postulates that a hydroperoxy-iron species (or a protonated analogue of this species) is a viable electrophilic oxidant in addition to the consensus oxidant, iron-oxo.     

细胞色素P450的电化学研究进展

细胞色素P450的电化学研究从一个侧面反映了为使细胞色素P450达到工业催化剂的最终目的人们所作的不懈努力。本文从细胞色素P450在电极上的电子转移研究,隧道扫描显微镜的微观成像研究和使用电极作为细胞色素P450的电子给体从而实现细胞色素P450底物转化三方面,评述了近年来细胞色素P450的电化学研究进展。     

Monooxygenation mechanism by cytochrome p-450

The substrate oxygenation mechanism by an ultimate species in monooxygenation by cytochrome P-450 (compound I) was investigated by the density functional theory method. An initial model compound was constructed from a structure obtained by 300-ps molecular dynamics simulation of compound I-formed P-450cam under physiologic conditions, and it consisted of porphine for protoporphyrin IX, S(-)-CH(3) for the side chain of Cys357 of the fifth ligand of heme, a methane molecule for the substrate, a heme iron, and an oxygen atom of the sixth ligand of heme. The results of the calculation revealed that the substrate oxygenation mechanism had four elementary processes, i.e., (1) formation of [FeOH](3+) and a substrate radical by hydrogen atom abstraction from the substrate caused by [FeO](3+), (2) rotation of the OH group of the sixth ligand of [FeOH](3+) produced by process 1, (3) substrate radical binding with the [FeOH](3+), and (4) elimination of the oxygenated substrate formed at the sixth ligand binding site. The rate-determining step is process 1, hydrogen atom abstraction from the substrate, and the activation energy was determined to be about 15 kcal/mol. For this reason, it is thought that this reaction occurs in vivo.     

Mechanisms of reaction in cytochrome P450: Hydroxylation of camphor in P450cam

The fundamental nature of reactivity in cytochrome P450 enzymes is currently controversial. Modelling of bacterial P450cam has suggested an important role for the haem propionates in the catalysis, though this finding has been questioned. Understanding the mechanisms of this enzyme family is important both in terms of basic biochemistry and potentially in the prediction of drug metabolism. We have modelled the hydroxylation of camphor by P450cam, using combined quantum mechanics/molecular mechanics (QM/MM) methods. A set of reaction pathways in the enzyme was determined. We were able to pinpoint the source of the discrepancies in the previous results. We show that when a correct ionization state is assigned to Asp297, no spin density appears on the haem propionates and the protein structure in this region remains preserved. These results indicate that the haem propionates are not involved in catalysis.     

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.     

Aromatic hydroxylation by cytochrome P450: model calculations of mechanism and substituent effects

The mechanism and selectivity of aromatic hydroxylation by cytochrome P450 enzymes is explored using new B3LYP density functional theory computations. The calculations, using a realistic porphyrin model system, show that rate-determining addition of compound I to an aromatic carbon atom proceeds via a transition state with partial radical and cationic character. Reactivity is shown to depend strongly on ring substituents, with both electron-withdrawing and -donating groups strongly decreasing the addition barrier in the para position, and it is shown that the calculated barrier heights can be reproduced by a new dual-parameter equation based on radical and cationic Hammett sigma parameters.     

Stochastic simulations of the cytochrome P450 catalytic cycle

Cytochrome P450 enzymes involve a complex reaction cycle which has been described, for the first time, by a stochastic simulation of the system in the present work. A series of models are developed for a basic catalytic cycle, employing a set of microscopic rate constants for the oxidation of p-alkoxyacylanilides catalyzed by the cytochrome P450 1A2. By analyzing the effects of low concentrations of enzyme and substrate on the system, and the dependence of the system on several rate constants, it is discovered that the system evolves along relatively stable patterns from its initial state, as indicated from different runs of simulations. Strong fluctuations appear at the entrance and exit of the pathway, with very weak fluctuations in the middle sections of the cycle. Although noises are apparent when the reactant populations are very low, basically, the fundamental feature of the P450 cycle based on a microscopic view is that it is deterministic in nature. Meanwhile, the mathematical models we developed are qualitatively validated by a comparison with those experimental results of the P450 cycle. The findings of this work will be helpful for a further deeper understanding of the catalytic mechanism of cytochrome P450 enzymes.     

Biotransformation of the sesquiterpene (+)-valencene by cytochrome P450cam and P450BM-3

The sesquiterpenoids are a large class of naturally occurring compounds with biological functions and desirable properties. Oxidation of the sesquiterpene (+)-valencene by wild type and mutants of P450cam from Pseudomonas putida, and of P450BM-3 from Bacillus megaterium, have been investigated as a potential route to (+)-nootkatone, a fine fragrance. Wild type P450cam did not oxidise (+)-valencene but the mutants showed activities up to 9.8 nmol (nmol P450)(-1) min(-1), with (+)-trans-nootkatol and (+)-nootkatone constituting >85% of the products. Wild type P450BM-3 and mutants had higher activities (up to 43 min(-1)) than P450cam but were much less selective. Of the many products, cis- and trans-(+)-nootkatol, (+)-nootkatone, cis-(+)-valencene-1,10-epoxide, trans-(+)-nootkaton-9-ol, and (+)-nootkatone-13S,14-epoxide were isolated from whole-cell reactions and characterised. The selectivity patterns suggest that (+)-valencene has one binding orientation in P450cam but multiple orientations in P450BM-3.     

An evaluation of molecular models of the cytochrome P450 Streptomyces griseolus enzymes P450SU1 and P450SU2

Summary P450SU1 and P450SU2 are herbicide-inducible bacterial cytochrome P450 enzymes from Streptomyces griseolus. They have two of the highest sequence identities to camphor hydroxylase (P450cam from Pseudomonas putida), the cytochrome P450 with the first known crystal structure. We have built several models of these two proteins to investigate the variability in the structures that can occur from using different modeling protocols. We looked at variability due to alignment methods, backbone loop conformations and refinement methods. We have constructed two models for each protein using two alignment algorithms, and then an additional model using an identical alignment but different loop conformations for both buried and surface loops. The alignments used to build the models were created using the Needleman-Wunsch method, adapted for multiple sequences, and a manual method that utilized both a dotmatrix search matrix and the Needleman-Wunsch method. After constructing the initial models, several energy minimization methods were used to explore the variability in the final models caused by the choice of minimization techniques. Features of cytochrome P450cam and the cytochrome P450 superfamily, such as the ferredoxin binding site, the heme binding site and the substrate binding site were used to evaluate the validity of the models. Although the final structures were very similar between the models with different alignments, active-site residues were found to be dependent on the conformations of buried loops and early stages of energy minimization. We show which regions of the active site are the most dependent on the particular methods used, and which parts of the structures seem to be independent of the methods.     

Engineering and assaying of cytochrome P450 biocatalysts

Cytochrome P450s constitute a highly fascinating superfamily of enzymes which catalyze a broad range of reactions. They are essential for drug metabolism and promise industrial applications in biotechnology and biosensing. The constant search for cytochrome P450 enzymes with enhanced catalytic performances has generated a large body of research. This review will concentrate on two key aspects related to the identification and improvement of cytochrome P450 biocatalysts, namely the engineering and assaying of these enzymes. To this end, recent advances in cytochrome P450 development are reported and commonly used screening methods are surveyed.     

Epoxidation of olefins by hydroperoxo-ferric cytochrome P450

The T252A mutant of cytochrome P450cam is unable to form the oxoferryl "active oxygen" intermediate, as judged by its inability to hydroxylate its normal substrate, camphor. In the present study, we demonstrate that T252A P450cam is nonetheless able to epoxidize olefins, due to the action of a second oxidant. However, as shown in earlier radiolytic studies and by the ability of T252A to reduce dioxygen to hydrogen peroxide, the mutant retains the ability to form the hydroperoxo-ferric reaction cycle intermediate. The present results provide strong evidence that hydroperoxo-ferric P450 can serve as a second electrophilic oxidant capable of olefin epoxidation.     

Role of the cytochrome P450 NocL in nocardicin A biosynthesis

The nocardicins are members of a class of monocyclic beta-lactam antibiotics produced by the actinomycete Nocardia uniformis subsp. tsuyamanesis (ATCC 21806). The oxime moiety in nocardicin A is required for full biological activity and is rare among natural products. A putative gene cluster encoding the enzymes of the nocardicin biosynthetic pathway has been identified. Among these, NocL shows similarity to cytochromes P450. We describe its heterologous expression and demonstrate its role in the formation of this unusual structural feature in nocardicin A. This is the first gene from this pathway to be expressed and its function established, and it is the first example of oxime formation from a primary amine in bacteria unambiguously linked to a cytochrome P450.     

QM/MM modeling of compound I active species in cytochrome P450, cytochrome C peroxidase, and ascorbate peroxidase

QM/MM calculations provide a means for predicting the electronic structure of the metal center in metalloproteins. Two heme peroxidases, Cytochrome c Peroxidase (CcP) and Ascorbate Peroxidase (APX), have a structurally very similar active site, yet have active intermediates with very different electronic structures. We review our recent QM/MM calculations on these systems, and present new computational data. Our results are in good agreement with experiment, and suggest that the difference in electronic structure is due to a large number of small differences in structure from one protein to another. We also discuss recent QM/MM calculations on the active species of cytochrome P450, in which a similar sensitivity of the electronic structure to the environment is found. However, this does not appear to explain different catalytic profiles of the different drug-metabolizing isoforms of this class of enzyme.     

Solvation of the active site of cytochrome P450-cam

Summary Energetically favorable water binding sites in the substrate pocket of cytochrome P450-cam have been predicted by a molecular mechanics method. Binding sites corresponding to all the experimentally observed water sites in this region of the enzyme were located. The calculations also indicate the presence of two further water binding sites. One of these is located in a hydrophobic region of the protein where a water molecule would not bind tightly to the substrate-free enzyme. However, in the substrate-bound enzyme, a water molecule in this region could donate a hydrogen bond of optimum geometry to the carbonyl oxygen atom of the camphor substrate and could therefore contribute to the correct positioning of the comphor substrate for 5-exo-hydroxylation. These calculations also suggest that a steric analogue of camphor, containing an alkyl group which could prevent a water molecule from binding in this region, might inhibit cytochrome P450-cam by forming a more stable enzyme-ligand complex than camphor itself.     

Active species of horseradish peroxidase (HRP) and cytochrome P450: two electronic chameleons

The active site of HRP Compound I (Cpd I) is modeled using hybrid density functional theory (UB3LYP). The effects of neighboring amino acids and of environmental polarity are included. The low-lying states have porphyrin radical cationic species (Por(*)(+)). However, since the Por(*)(+) species is a very good electron acceptor, other species, which can be either the ligand or side chain amino acid residues, may participate in electron donation to the Por(*)(+) moiety, thereby making Cpd I behave like a chemical chameleon. Thus, this behavior that was noted before for Cpd I of P450 is apparently much more wide ranging than initially appreciated. Since chemical chameleonic behavior property was found to be expressed not only in the properties of Cpd I itself, but also in its reactivity, the roots of this phenomenon are generalized. A comparative discussion of Cpd I species follows for the enzymes HRP, CcP, APX, CAT (catalase), and P450.     

DNA-mediated assembly of cytochrome P450 BM3 subdomains

Cytochrome P450 BM3 is a versatile enzyme, which holds great promise for applications in biocatalysis and biomedicine. We here report on the generation of a hybrid DNA-protein device based on the two subdomains of BM3, the reductase domain BMR and the porphyrin domain BMP. Both subdomains were fused genetically to the HaloTag protein, a self-labeling enzyme, allowing for the bioconjugation with chloroalkane-modified oligonucleotides. The subdomain-DNA-chimeras could be reassembled by complementary oligonucleotides, thus leading to reconstitution of the monooxygenase activity of BM3 holoenzyme, as demonstrated by conversion of the reporter substrate 12-pNCA. Arrangement of the two chimeras on a switchable DNA scaffold allowed one to control the distance between both subdomains, as indicated by the DNA-dependent activity of the holoenzyme. Furthermore, a switchable chimeric device was constructed, in which monooxygenase activity could be turned off by DNA strand displacement. This study demonstrates that P450 BM3 engineering and strategies of DNA nanotechnology can be merged to open up novel ways for the development of novel screening systems or responsive catalysts with potential applications in drug delivery.     

Structural fine-tuning of a multifunctional cytochrome P450 monooxygenase

AurH is a unique cytochrome P450 monooxygenase catalyzing the stepwise formation of a homochiral oxygen heterocycle, a key structural and pharmacophoric component of the antibiotic aureothin. The exceptional enzymatic reaction involves a tandem oxygenation process including a regio- and stereospecific hydroxylation, followed by heterocyclization. For the structural and biochemical basis of this unparalleled sequence, four crystal structures of AurH variants in different conformational states and in complex with the P450 inhibitor ancymidol were solved, which represent the first structures of the CYP151A group. Structural data in conjunction with computational docking, site-directed mutagenesis, and chemical analyses unveiled a switch function when recognizing the two substrates, deoxyaureothin and the hydroxylated intermediate, thus allowing the second oxygenation-heterocyclization step. Furthermore, we were able to modify the chemo- and regioselectivity of AurH, yielding mutants that catalyze the regioselective six-electron transfer of a nonactivated methyl group to a carboxylic acid via hydroxyl and aldehyde intermediates.     

Molecular modeling of cytochrome P450 3A4

The three-dimensional structure of human cytochrome P450 3A4 was modeled based on crystallographic coordinates of four bacterial P450s: P450 BM-3, P450cam, P450terp, and P450eryF. The P450 3A4 sequence was aligned to those of the known proteins using a structure-based alignment of P450 BM-3, P450cam, P450terp, and P450eryF. The coordinates of the model were then calculated using a consensus strategy, and the final structure was optimized in the presence of water. The P450 3A4 model resembles P450 BM-3 the most, but the B helix is similar to that of P450eryF, which leads to an enlarged active site when compared with P450 BM-3, P450cam, and P450terp. The 3A4 residues equivalent to known substrate contact residues of the bacterial proteins and key residues of rat P450 2B1 are located in the active site or the substrate access channel. Docking of progesterone into the P450 3A4 model demonstrated that the substrate bound in a 6-orientation can interact with a number of active site residues, such as 114, 119, 301, 304, 305, 309, 370, 373, and 479, through hydrophobic interactions. The active site of the enzyme can also accommodate erythromycin, which, in addition to the residues listed for progesterone, also contacts residues 101, 104, 105, 214, 215, 217, 218, 374, and 478. The majority of 3A4 residues which interact with progesterone and/or erythromycin possess their equivalents in key residues of P450 2B enzymes, except for residues 297, 480 and 482, which do not contact either substrate in P450 3A4. The results from docking of progesterone and erythromycin into the enzyme model make it possible to pinpoint residues which may be important for 3A4 function and to target them for site-directed mutagenesis.