High-level expression of recombinant rabbit cytochrome P450 2E1 in Escherichia coli C41 and its purification

Recombinant production of eukaryotic cytochrome P450s in microbial cell factories

Cytochrome P450s (P450s) comprise one of the largest known protein families. They occur in every kingdom of life and catalyze essential reactions, such as carbon source assimilation, synthesis of hormones and secondary metabolites, or degradation of xenobiotics. Due to their outstanding ability of specifically hydroxylating complex hydrocarbons, there is a great demand to use these enzymes for biocatalysis, including applications at an industrial scale. Thus, the recombinant production of these enzymes is intensively investigated. However, especially eukaryotic P450s are difficult to produce. Challenges are faced due to complex cofactor requirements and the availability of a redox-partner (cytochrome P450 reductase, CPR) can be a key element to get active P450s. Additionally, most eukaryotic P450s are membrane bound which complicates the recombinant production. This review describes current strategies for expression of P450s in the microbial cell factories Escherichia coli, Saccharomyces cerevisiae, and Pichia pastoris.

Heterologous expression of fungal cytochromes P450 (CYP5136A1 and CYP5136A3) from the white-rot basidiomycete Phanerochaete chrysosporium: Functionalization with cytochrome b5 in Escherichia coli

Cytochromes P450 from the white-rot basidiomycete Phanerochaete chrysosporium, CYP5136A1 and CYP5136A3, are capable of catalyzing oxygenation reactions of a wide variety of exogenous compounds, implying their significant roles in the metabolism of xenobiotics by the fungus. It is therefore interesting to explore their biochemistry to better understand fungal biology and to enable the use of fungal enzymes in the biotechnology sector. In the present study, we developed heterologous expression systems for CYP5136A1 and CYP5136A3 using the T7 RNA polymerase/promoter system in Escherichia coli. Expression levels of recombinant P450s were dramatically improved by modifications and optimization of their N-terminal amino acid sequences. A CYP5136A1 reaction system was reconstructed in E. coli whole cells by coexpression of CYP5136A1 and a redox partner, NADPH-dependent P450 reductase (CPR). The catalytic activity of CYP5136A1 was significantly increased when cytochrome b5 (Cyt-b5) was further coexpressed with CPR, indicating that Cyt-b5 supports electron transfer reactions from NAD(P)H to CYP5136A1. Notably, P450 reaction occurred in E. coli cells that harbored CYP5136A1 and Cyt-b5 but not CPR, implying that the reducing equivalents required for the P450 catalytic cycle were transferred via a CPR-independent pathway. Such an "alternative" electron transfer system in CYP5136A1 reaction was also demonstrated using purified enzymes in vitro. The fungal P450 reaction system may be associated with sophisticated electron transfer pathways.

Environmentally persistent free radical-containing particulate matter competitively inhibits metabolism by cytochrome P450 1A2

Combustion processes generate different types of particulate matter (PM) that can have deleterious effects on the pulmonary and cardiovascular systems. Environmentally persistent free radicals (EPFRs) represent a type of particulate matter that is generated after combustion of environmental wastes in the presence of redox-active metals and aromatic hydrocarbons. Cytochromes P450 (P450/CYP) are membrane-bound enzymes that are essential for the phase I metabolism of most lipophilic xenobiotics. The EPFR formed by chemisorption of 2-monochlorophenol to silica containing 5% copper oxide (MCP230) has been shown to generally inhibit the activities of different forms of P450s without affecting those of cytochrome P450 reductase and heme oxygenase-1. The mechanism of inhibition of rat liver microsomal CYP2D2 and purified rabbit CYP2B4 by MCP230 has been shown previously to be noncompetitive with respect to substrate. In this study, MCP230 was shown to competitively inhibit metabolism of 7-benzyl-4-trifluoromethylcoumarin and 7-ethoxyresorufin by the purified, reconstituted rabbit CYP1A2. MCP230 is at least 5- and 50-fold more potent as an inhibitor of CYP1A2 than silica containing 5% copper oxide and silica, respectively. Thus, even though PM generally inhibit multiple forms of P450, PM interacts differently with the forms of P450 resulting in different mechanisms of inhibition. P450s function as oligomeric complexes within the membrane. We also determined the mechanism by which PM inhibited metabolism by the mixed CYP1A2-CYP2B4 complex and found that the mechanism was purely competitive suggesting that the CYP2B4 is dramatically inhibited when bound to CYP1A2.

Inhibition of cytochrome P450 2B4 by environmentally persistent free radical-containing particulate matter

Combustion processes generate particulate matter (PM) that can affect human health. The presence of redox-active metals and aromatic hydrocarbons in the post-combustion regions results in the formation of air-stable, environmentally persistent free radicals (EPFRs) on entrained particles. Exposure to EPFRs has been shown to negatively influence pulmonary and cardiovascular functions. Cytochromes P450 (P450/CYP) are endoplasmic reticulum resident proteins that are responsible for the metabolism of foreign compounds. Previously, it was shown that model EPFRs, generated by exposure of silica containing 5% copper oxide (CuO-Si) to either dicholorobenzene (DCB230) or 2-monochlorophenol (MCP230) at ≥ 230 °C, inhibited six forms of P450 in rat liver microsomes (Toxicol. Appl. Pharmacol. (2014) 277:200-209). In this study, the inhibition of P450 by MCP230 was examined in more detail by measuring its effect on the rate of metabolism of 7-ethoxy-4-trifluoromethylcoumarin (7EFC) and 7-benzyloxyresorufin (7BRF) by the purified, reconstituted CYP2B4 system. MCP230 inhibited the CYP2B4-mediated metabolism of 7EFC at least 10-fold more potently than non-EPFR controls (CuO-Si, silica, and silica generated from heating silica and MCP at 50 °C, so that EPFRs were not formed (MCP50)). The inhibition by EPFRs was specific for the P450 and did not affect the ability of the redox partner, P450 reductase (CPR) from reducing cytochrome c. All of the PM inhibited CYP2B4-mediated metabolism noncompetitively with respect to substrate. When CYP2B4-mediated metabolism of 7EFC was measured as a function of the CPR concentration, the mechanism of inhibition was competitive. EPFRs likely inhibit CYP2B4-mediated substrate metabolism by physically disrupting the CPR·P450 complex.

Structure of pyrazole derivatives impact their affinity, stoichiometry, and cooperative interactions for CYP2E1 complexes

CYP2E1 plays a critical role in detoxification and carcinogenic activation of drugs, pollutants, and dietary compounds; however, these metabolic processes can involve poorly characterized cooperative interactions that compromise the ability to understand and predict CYP2E1 metabolism. Herein, we employed an array of ten azoles with an emphasis on pyrazoles to establish the selectivity of catalytic and cooperative CYP2E1 sites through binding and catalytic studies. Spectral binding studies for monocyclic azoles suggested two binding events, while bicyclic azoles suggested one. Pyrazole had moderate affinity toward the CYP2E1 catalytic site that improved when a methyl group was introduced at either position 3 or 4. The presence of methyl groups simultaneously at positions 3 and 5 blocked binding, and a phenyl group at position 3 did not improve binding affinity. In contrast, pyrazole fusion to a benzene or cyclohexane ring greatly increased affinity. The consequences of these binding events on CYP2E1 catalysis were studied through inhibition studies with 4-nitrophenol, a substrate known to bind both sites. Most pyrazoles shared a common mixed cooperative inhibition mechanism in which pyrazole binding rescued CYP2E1 from substrate inhibition. Overall, inhibitor affinities toward the CYP2E1 catalytic site were similar to those reported in binding studies, and the same trend was observed for binding at the cooperative site. Taken together, these studies identified key structural determinants in the affinity and stoichiometry of azole interactions with CYP2E1 and consequences on catalysis that further advance an understanding of the relationship between structure and function for this enzyme.

Effect of homomeric P450-P450 complexes on P450 function

Previous studies have shown that the presence of one P450 enzyme can affect the function of another. The goal of the present study was to determine if P450 enzymes are capable of forming homomeric complexes that affect P450 function. To address this problem, the catalytic activities of several P450s were examined in reconstituted systems containing NADPH-POR (cytochrome P450 reductase) and a single P450. CYP2B4 (cytochrome P450 2B4)-, CYP2E1 (cytochrome P450 2E1)- and CYP1A2 (cytochrome P450 1A2)-mediated activities were measured as a function of POR concentration using reconstituted systems containing different concentrations of P450. Although CYP2B4-dependent activities could be explained by a simple Michaelis-Menten interaction between POR and CYP2B4, both CYP2E1 and CYP1A2 activities generally produced a sigmoidal response as a function of [POR]. Interestingly, the non-Michaelis behaviour of CYP1A2 could be converted into a simple mass-action response by increasing the ionic strength of the buffer. Next, physical interactions between CYP1A2 enzymes were demonstrated in reconstituted systems by chemical cross-linking and in cellular systems by BRET (bioluminescence resonance energy transfer). Cross-linking data were consistent with the kinetic responses in that both were similarly modulated by increasing the ionic strength of the surrounding solution. Taken together, these results show that CYP1A2 forms CYP1A2-CYP1A2 complexes that exhibit altered catalytic activity.

What is the Functional Role of N-terminal Transmembrane Helices in the Metabolism Mediated by Liver Microsomal Cytochrome P450 and its Reductase?

We sought to clarify on the hitherto unresolved role of N-terminal transmembrane segments (TMS) of cytochrome P450 (CYP) and its' reductase (CPR) in protein interaction/catalysis. TMS analyses show little evolutionary conservation in CYPs. The conserved CPR's TMS poses limited scope for predictable/consistent hetero-recognition with the wide bevy of CYPs' TMS, as evident from preliminary analyses and TMhit server predictions for inter-helical binding. Further, experimentations with four different CPR preparations (preps) and two liver microsomal CYPs (2C9 and 2E1) shows that the hydroxylated product formation rate is not quantitatively correlated to the extent of integrity of the CPR N-terms. Incorporation of cytochrome b (5) in some reactions afforded similar rates while employing either fully intact or partially intact CPR. A survey of literature shows that liver microsomal CYPs function quite well even without the TMS or with significantly altered TMS. These observations negate the hypothesis that N-term TMS of CPR or CYP is obligatory for CYP-CPR interaction and catalysis. Also, in CYP2E1-mediated hydroxylation of para-nitrophenol, the extent of intactness or truncation did not significantly affect the CPR preps' catalytic role at very low or high substrate concentrations. To interpret these results, we draw support from recently published research on reduced nicotinamide adenide dinucleotide phosphate oxidase (Takac et al., J Biol Chem, 286:13304-13313, 2011) and from our pertinent earlier works. We infer that CPR' free TMS segment could alter the diffusible reactive oxygen species' dynamics in the microenvironment, thereby altering the reaction outcome. Based on the evidence, we conclude that TMS merely facilitates "interaction/catalysis" by anchoring the CYP and CPR in the lipid interface.

Heterologous expression and mechanistic investigation of a fungal cytochrome P450 (CYP5150A2): involvement of alternative redox partners

A fungal cytochrome P450 monooxygenase (CYP5150A2) from the white-rot basidiomycete Phanerochaete chrysosporium was heterologously expressed in Escherichia coli and purified as an active form. The purified CYP5150A2 was capable of hydroxylating 4-propylbenzoic acid (PBA) with NADPH-dependent cytochrome P450 oxidoreductase (CPR) as the single redox partner; the reaction efficiency was improved by the addition of electron transfer protein cytochrome b5 (Cyt-b5). Furthermore, CYP5150A2 exhibited substantial activity with redox partners Cyt-b5 and NADH-dependent Cyt-b5 reductase (CB5R) even in the absence of CPR. These results indicated that a combination of CB5R and Cyt-b5 may be capable of donating both the first and the second electrons required for the monooxygenation reaction. Under reaction conditions in which the redox system was associated with the CB5R-dependent Cyt-b5 reduction system, the exogenous addition of CPR and NADPH had no effect on the PBA hydroxylation rate or on coupling efficiency, indicating that the transfer of the second electron from Cyt-b5 was the rate-limiting step in the monooxygenase system. In addition, the rate of PBA hydroxylation was significantly dependent on Cyt-b5 concentration, exhibiting Michaelis-Menten kinetics. This study provides indubitable evidence that the combination of CB5R and Cyt-b5 is an alternative redox partner facilitating the monooxygenase reaction catalyzed by CYP5150A2.

An expression of an insect membrane-bound cytochrome P450 CYP6AA3 in the Escherichia coli in relation to insecticide resistance in a malarial vector

This laboratory investigation was carried out at the Faculty of Sciences, Mahidol University, Thailand during October 2007 to May 2009. The objectives of this study include: the search for heterologous expression of the cytochrome P450 CYP6AA3 enzyme of the Anopheles minimus mosquitoes in relation to Malaria disease and to provide some information on molecular mechanism of insects' pyrethroid resistance. The polymerase chain reaction aided by the Pfu DNA polymerase and some specific generated primers were used to modify the CYP6AA3 gene. The PCR product was ligated with a predigested pET-3a at the NdeI and BamHI restriction sites. The modified CYP6AA3 enzyme was expressed in the Escherichia coli BL21 (DE3) pLysS in order to achieve a high amount of soluble form of its expression. The results showed that the use of the isopropyl-beta-D-thiogalactopyranoside (IPTG) and incubation together with ferric chloride and delta-aminolevulinic acid did not increase any soluble form of the CYP6AA3 enzyme. A significant amount of soluble enzyme was produced upon the replacement of the 30 N-terminal residues with a short peptide where it gave Ldelta30CYP6AA3 protein and after purification process was taken place, it yielded up to 10.64 mg 10 L(-1) or approximately 1 mg L(-1) of the homogenous Ldelta30CYP6AA3. When this purified Ldelta30CYP6AA3 protein was used in a metabolizing process with the cypermethrin, deltamethrin and permethrin substrates, it gave their apparent Km values for cypermethrin and deltamethrin of 12.5 and 23.5 microM, respectively. The heterologous expression carried out with the use of the E. coli gave a high amount of soluble CYP6AA3 enzyme of the An. minimus mosquitoes hence the modified technique being used was successfully achieved.

Single mutations change CYP2F3 from a dehydrogenase of 3-methylindole to an oxygenase

Pulmonary cytochrome P450 2F3 (CYP2F3) catalyzes the dehydrogenation of the pneumotoxin 3-methylindole (3MI) to an electrophilic intermediate, 3-methyleneindolenine, which is responsible for the toxicity of the parent compound. Members of the CYP2F subfamily are the only enzymes known to exclusively dehydrogenate 3MI, without detectable formation of oxygenation products. Thus, CYP2F3 is an attractive model to study dehydrogenation mechanisms. The purpose of this study was to identify specific residues that could facilitate 3MI dehydrogenation. Both single and double mutations were constructed to study the molecular mechanisms that direct dehydrogenation. Double mutations in substrate recognition sites (SRS) 1 produced an inactive enzyme, while double mutants in SRS 4 did not alter 3MI metabolism. However, double mutations in SRS 5 and SRS 6 successfully introduced oxygenase activity to CYP2F3. Single mutations in SRS 5, SRS 6 and near SRS 2 also introduced 3MI oxygenase activity. Mutants S474H and D361T oxygenated 3MI but also increased dehydrogenation rates, while G214L, E215Q and S475I catalyzed 3MI oxygenation exclusively. A homology model of CYP2F3 was precisely consistent with specific dehydrogenation of 3MI via initial hydrogen atom abstraction from the methyl group. In addition, intramolecular kinetic deuterium isotope studies demonstrated an isotope effect ( K H/ K D) of 6.8. This relatively high intramolecular deuterium isotope effect confirmed the initial hydrogen abstraction step; a mutant (D361T) that retained the dehydrogenation reaction exhibited the same deuterium isotope effect. The results showed that a single alteration, such as a serine to isoleucine change at residue 475, dramatically switched catalytic preference from dehydrogenation to oxygenation.

Heteromeric complex formation between CYP2E1 and CYP1A2: evidence for the involvement of electrostatic interactions

Mixed reconstituted systems containing CYP2B4, CYP1A2, and NADPH-cytochrome P450 reductase were previously shown to exhibit a dramatic inhibition of 7-pentoxyresorufin O-dealkylation (PROD) when compared to simple reconstituted systems containing reductase and a single P450 enzyme, results consistent with the formation of CYP1A2-CYP2B4 complexes where the reductase binds with high affinity to the CYP1A2 moiety of the complex. In this report, we provide evidence for an interaction between CYP1A2 and CYP2E1. Synergism of 7-ethoxyresorufin O-deethylation (EROD) and PROD was observed when these P450s were combined in mixed reconstituted systems at subsaturating reductase concentrations. Higher ionic strength attenuated the synergistic stimulation of both PROD and EROD in mixed reconstituted systems, consistent with disruption of heteromeric CYP2E1-CYP1A2 complexes. The effect of ionic strength was further examined as a function of reductase concentration. At lower ionic strength, there was a significant synergistic stimulation of EROD. This synergistic stimulation diminished with increasing reductase concentration, resulting in an additive response as reductase became saturating. Interestingly, at high ionic strength, the synergism of EROD in the mixed reconstituted system was not observed. In contrast, mixed reconstituted systems containing CYP2E1 and CYP2B4 did not provide evidence for the formation of these heteromeric P450-P450 complexes. The synergistic stimulation observed with the reductase-CYP1A2-CYP2E1 mixed reconstituted system is consistent with the formation of a CYP1A2-CYP2E1 complex. Taken together with the lack of a kinetically detectable interaction between CYP2B4 and CYP2E1, and the previously reported CYP1A2-CYP2B4 interaction, these results suggest that CYP1A2 may facilitate the formation of complexes with other P450 enzymes.

CYP2E1 substrate inhibition. Mechanistic interpretation through an effector site for monocyclic compounds

In this study we offer a mechanistic interpretation of the previously known but unexplained substrate inhibition observed for CYP2E1. At low substrate concentrations, p-nitrophenol (pNP) was rapidly turned over (47 min(-1)) with relatively low K(m) (24 microM); nevertheless, at concentrations of >100 microM, the rate of pNP oxidation gradually decreased as a second molecule bound to CYP2E1 through an effector site (K(ss) = 260 microm), which inhibited activity at the catalytic site. 4-Methylpyrazole (4MP) was a potent inhibitor for both sites through a mixed inhibition mechanism. The K(i) for the catalytic site was 2.0 microM. Although we were unable to discriminate whether an EIS or ESI complex formed, the respective inhibition constants were far lower than K(ss). Bicyclic indazole (IND) inhibited catalysis through a single CYP2E1 site (K(i) = 0.12 microM). Similarly, 4MP and IND yielded type II binding spectra that reflected the association of either two 4MP or one IND molecule(s) to CYP2E1, respectively. Based on computational docking studies with a homology model for CYP2E1, the two sites for monocyclic molecules, pNP and 4MP, exist within a narrow channel connecting the active site to the surface of the enzyme. Because of the presence of the heme iron, one site supports catalysis, whereas the other more distal effector site binds molecules that can influence the binding orientation and egress of molecules for the catalytic site. Although IND did not bind these sites simultaneously, the presence of IND at the catalytic site blocked binding at the effector site.

Global analysis of protein-protein interactions reveals multiple CYP2E1-reductase complexes

Although a single binary functional complex between cytochrome P450 (P450 or CYP for a specific isoform) and cytochrome P450 reductase (CPR) has been generally accepted in the literature, this simple model failed to explain the experimentally observed catalytic activity of recombinant CYP2E1 in dependence on the total concentration of the added CPR-K56Q mutant. Our rejection of the simplest 1:1 binding model was based on two independent lines of experimental evidence. First, under the assumption of the 1:1 binding model, separate analyses of titration curves obtained while varying either P450 or CPR concentrations individually produced contradictory results. Second, an asymmetric Job plot suggested the existence of higher order molecular complexes. To identify the most probable complexation mechanism, we generated a comprehensive data set where the concentrations of both P450 and P450 were varied simultaneously, rather than one at a time. The resulting two-dimensional data were globally fit to 32 candidate mechanistic models, involving the formation of binary, ternary, and quaternary P450.CPR complexes, in the absence or presence or P450 and CPR homodimers. Of the 32 candidate models (mechanisms), two models were approximately equally successful in explaining our experimental data. The first plausible model involves the binary complex P450.CPR, the quaternary complex (P450)2.(CPR)2, and the homodimer (P450)2. The second plausible model additionally involves a weakly bound ternary complex (P450)2.CPR. Importantly, only the binary complex P450.CPR seems catalytically active in either of the two most probable mechanisms.

CYP2E1 active site residues in substrate recognition sequence 5 identified by photoaffinity labeling and homology modeling

Despite its biological importance, our knowledge of active site structure and relevance of critical amino acids in CYP2E1 catalytic processes remain limited. In this study, we identified CYP2E1 active site residues using photoaffinity labeling with 7-azido-4-methylcoumarin (AzMC) coupled with a CYP2E1 homology model. In the absence of light, AzMC was an effective competitor against substrate p-nitrophenol oxidation by CYP2E1. Photoactivation of AzMC led to a concentration-dependent loss in CYP2E1 activity and structural integrity resulting from the modification of both heme and protein. The photo-labeling reaction degraded heme and produced a possible heme adduct. Probe incorporation into the protein occurred at multiple sites within substrate recognition sequence 5 (SRS-5). Based on a CYP2E1 homology model, we hypothesize AzMC labels SRS-5 residues, Leu363, Val364, and Leu368, in the active site. In addition, we propose a series of phenylalanines, especially Phe106, mediate contacts with the coumarin.

Altering the regioselectivity of cytochrome P450 CYP102A3 of Bacillus subtilis by using a new versatile assay system

A novel monooxygenase (CYP102A3) has been discovered within the Bacillus subtilis genome that reveals a similarity of 76 % to the well-known cytochrome P450 BM-3 of B. megaterium (CYP102A1). Both enzymes are natural fusion proteins consisting of a heme domain and a FAD/FMN-reductase domain. Because of their high turnover rates, these biocatalysts are of special interest for industrial applications, but show only limited regioselectivity. In this work, the regioselectivity of CYP102A3 was changed by directed evolution and protein design to hydroxylate substrates not only in different subterminal, but also to a high extent, in terminal carbon chain positions. To enable a high-throughput screening procedure, a very versatile assay was developed that is capable of discriminating between terminal and subterminal hydroxylation of carbon chains. A double mutant of CYP102A3 was obtained that produces 48 % octan-1-ol as the main product of the reaction.