Emma A. Hall and
Stephen G. Bell*
School of Chemistry and Physics, University of Adelaide, SA 5005, Australia. E-mail: stephen.bell@adelaide.edu.au; Fax: +61 8 83034380; Tel: +61 8 83134822
First published on 12th December 2014
CYP101B1 from Novosphingobium aromaticivorans DSM12444 is a homologue of CYP101A1 (P450cam) from Pseudomonas putida and the CYP101D1, CYP101D2 and CYP101C1 enzymes from the same bacterium. CYP101B1 binds norisoprenoids more tightly than camphor and efficiently hydroxylates substrates in combination with ferredoxin reductase, ArR, and [2Fe–2S] ferredoxin, Arx, electron transfer partners. The norisoprenoids, α-ionone and β-damascone are both oxidised by CYP101B1 with high product formation activity, >500 min−1. α-Ionone oxidation occurred regioselectively at the allylic C3 position while β-damascone was hydroxylated predominantly at C3, 86%, with the main competing minor product arising from oxidation at the allylic C4 position (11%). When incorporated into a whole-cell oxidation system, with ArR and Arx, CYP101B1 is also capable of oxidising the aromatic compound indole. Other aromatic molecules including phenylcyclohexane and p-cymene were tested and both were hydroxylated by CYP101B1. Phenylcyclohexane was selectively oxidised to trans-4-phenylcyclohexanol while p-cymene was hydroxylated at the benzylic carbons to yield a mixture of isopropylbenzyl alcohol and p-α,α-trimethylbenzylalcohol. Trans-4-Phenylcyclohexanol was formed with a product formation rate of 141 min−1 and was five times more active than the oxidation of p-cymene.
We have reported the heterologous production and purification of the majority of the sixteen P450 enzymes from Novosphingobium aromaticivorans DSM12444.5,9 Potential substrates for many of these enzymes have been identified and these span a broad range of organic molecules including terpenoids, linear alkanes and polyaromatic hydrocarbons.5,10,11 The CYP enzymes of this bacterium are presumably involved in metabolic processes which enable it to survive in the oligotrophic environments where it lives.12,13 A class I electron transfer system, consisting of a flavin-dependent ferredoxin reductase, ArR, and a [2Fe–2S] ferredoxin, Arx, has been identified which is able to reconstitute the monooxygenase activity of at least five of these enzymes (CYP101B1, CYP101C1, CYP101D1, CYP101D2 and CYP111A2).5,9,14 Whole-cell systems capable of product formation on the gram-per-litre scale in shake flasks have also been constructed.9 Like P450cam (CYP101A1), CYP101D1 and CYP101D2 both bind and oxidise camphor (yielding ≥98% 5-exo-hydroxycamphor). While CYP101B1 catalyses the unselective oxidation of camphor generating four products, CYP101C1 does not. CYP101C1 and CYP101B1 are able to bind and oxidise β-ionone. Crystal structures of CYP101D1, CYP101D2 and CYP101C1 and the ArR and Arx electron transfer partners reveal striking similarities between these CYP101 family enzymes and CYP101A1.15–18 The enzymes also provide templates in order to understand how variations in the amino acid sequence of closely related P450 enzymes can lead to different substrate binding preferences, altered product selectivity and a predilection for specific electron transfer partners.15–17
CYP101B1 when combined with ArR and Arx in a whole-cell oxidation system was found to turn the culture medium blue indicating that it is able to hydroxylate indole to 3-hydroxyindole which can then dimerise to form indigo.9,19 Novosphingobium species are able to degrade a wide variety of aromatic hydrocarbons and therefore CYP101B1 has the potential to be a biocatalyst for the oxidation of these hydrophobic substrates.9,12,13 Here we report our investigations into the substrate range of CYP101B1 and show it is a highly active and regioselective biocatalyst for the oxidation of α-ionone and β-damascone. We also investigate the ability of CYP101B1 to oxidise hydrophobic aromatic ring containing substrates such as phenylcyclohexane and p-cymene.
Despite the relatively low sequence similarity with CYP101B2 from N. tardaugens a comparison of the residues that are likely to make up the active site and substrate access channel revealed the almost complete conservation of amino acids (Fig. 1 and Table S2†). The only difference was the conservative replacement of Ile236 in CYP101B1 with Val244 in CYP101B2. This is the equivalent position to Val247 in CYP101A1 a residue located in the active site which is important for substrate binding.24 This comparison also revealed important differences in the active sites of the two CYP101B subfamily enzymes when compared to the other members of the CYP101 family. For example sequence alignment reveals that, Val91 aligns with Phe and Trp residues (Phe87; CYP101A1 numbering given) and His100 with Tyr residues (Tyr96). However there are similarities with Thr173 aligning with Thr residues (Thr185) and Val284 and Val385 with Val residues (Val295 and Val396; Table S2†).25,26 Additionally residues such as Thr85 and Ile384 (Thr101 and Ile395) align with other Thr and Ile residues in all the CYP101 family members except CYP101C1.16
In CYP101B1 there are also changes to the residues that comprise the potassium binding site and salt bridge network of CYP101A1. Glu84 and Tyr96 of the CYP1011A1 potassium binding site are replaced with Arg72 and His80, respectively in CYP101B1. In the salt bridge network most of the residues are conserved but there are modifications with Asp97 being replaced by Gly81 and Lys178 by Arg166. These differences when taken together with those in active site are likely to account for the lower affinity for camphor and the altered substrate profile and C–H bond oxidation selectivity of CYP101B1.5 It is of note that the number of different active site residues in CYP101B1 compared to CYP101A1 is greater than in both CYP101D enzymes but fewer when compared to CYP101C1 (Table S2† and Fig. 1 and 2).16
The acid alcohol pair of the I-helix, which is involved in proton delivery and oxygen activation, and the glutamate and arginine pair (EXXR), which is proposed to be involved in heme incorporation, are conserved in CYP101B1.27 The normally conserved phenylalanine amino acid seven residues towards the N-terminus from the heme iron ligating cysteine, Cys346, is also present (Fig. 1).28 One potentially interesting point of difference in the CYP101B1 sequence is the presence of a proline residue, Pro347, one residue to the C-terminal side of the proximal cysteine. Proline occurs naturally at this position in chloroperoxidase and P450s such as CYP7A1, and CYP121A1.27,29,30 The introduction of proline residues at this position is known to have distinctive effects in other CYP enzymes, with the Leu358Pro mutant of CYP101A1 exhibiting the same spectral perturbations induced by putidaredoxin binding.31,32 The Ile401Pro mutant of CYP102A1 (P450Bm-3) displays enhanced oxidation rates with non-natural substrates and exhibits unusual electronic properties and structural changes.33,34
![]() | ||
Fig. 3 UV/Vis analysis of the CYP101B1 spin state shift with (a) α-ionone (b) β-damascone (c) phenylcyclohexane and (d) p-cymene. |
![]() | ||
Fig. 4 Substrate dissociation constant analysis of (a) α-ionone, (b) β-damascone, (c) phenylcyclohexane and (d) p-cymene with CYP101B1. |
CYP101B1/substrate | %HS heme | Kd (μM) | NADH consumption frequency (min−1) | Product formation rate (min−1) | Coupling (%) |
---|---|---|---|---|---|
β-Ionone | ≥95% | 0.23 ± 0.1 | 1600 ± 100 | 1010 ± 60 | 63 |
α-Ionone | ≥95% | 0.26 ± 0.04 | 1380 ± 140 | 663 ± 56 | 48 |
β-Damascone | 80% | 8.3 ± 0.9 | 931 ± 13 | 562 ± 12 | 60 |
Phenylcyclohexane | 20% | 7.8 ± 0.9 | 293 ± 9.0 | 141 ± 17 | 48 |
p-Cymene | 5% | 23.0 ± 3.0 | 197 ± 27 | 25.3 ± 6.0 | 13 |
The oxidation of β-damascone generated one major product (∼86%) and two minor products (Scheme 1 and Fig. 5). GC analysis revealed two major products from α-ionone oxidation in 63:
37% ratio (Fig. 5). The product formation rate of α-ionone and β-damascone were similar, 663 ± 56 min−1 and 562 ± 12 min−1, with the coupling of NADH to product formation for β-damascone (60%) compensating for its lower NADH turnover frequency when compared to α-ionone (Table 1).
Using a whole-cell oxidation system which produced CYP101B1, Arx and ArR in E. coli we were able to generate 100–200 mg of a mixture of the products formed from each norisoprenoid. The whole-cell biotransformation of α-ionone generates three products all in roughly equal amounts (30:
37
:
33% ratio, Fig. 5). The first two of these correspond to the two products observed in vitro and have masses indicating they are hydroxylated monooxygenase products (Mass 208.2 AMU). The third product, which was not observed in the in vitro turnovers, has a mass indicating it arises from further oxidation of an alcohol product to a ketone (Mass 206.1 AMU). The whole-cell oxidation of β-damascone and β-ionone yielded predominantly one main product (Fig. 5) which corresponded to the products of the in vitro turnover systems. β-Damascone oxidation produced two minor products in low yields. Purification and isolation of the individual products was carried out using silica gel chromatography with only the fractions containing predominantly a single compound being taken forward for characterisation. We isolated between 2–10 mg pure product from each turnover of β-ionone, α-ionone and β-damascone after column chromatography, with the remaining product being found in mixed fractions. These pure samples were used to identify the products using NMR and mass spectrometry (Experimental section and ESI†). The products of α-ionone oxidation were identified as cis- and trans-3-hydroxy-α-ionone and 3-oxo-α-ionone (Scheme 1). The lack of significant quantities of the further oxidation product 3-oxo-α-ionone in the in vitro enzyme assays suggests it is only produced when the amount of α-ionone substrate is limited relative to the two hydroxylated products or that the increased solubility of the alcohols allow them to enter the cells more readily than the less soluble substrate. Furthermore the proportion of the cis and trans isomers were different in the whole-cell oxidation turnover suggesting that the trans isomer is more readily converted to 3-oxo-α-ionone. The regioselective oxidation of α-ionone at C3 and the absence of any evidence of the epoxidation of the double bond suggests that the substrate must be held in the active site with both the C3 pseudoequatorial and pseudoaxial C–H bonds held close to the heme iron.
β-Damascone oxidation by CYP101B1 yielded 86% of the majority product, 3-hydroxy-β-damascone, and 11% of the more abundant minor product, 4-hydroxy-β-damascone (Scheme 1). The other minor product that was produced (≤3% of an in vitro reaction) was not generated in a sufficiently pure enough form in high enough yields for NMR characterisation. GC-MS analysis indicated that it is likely to be a further oxidation product of one of the alcohols (Mass 206.2 AMU). These results are similar to those obtained previously with β-ionone. The major product from β-ionone was confirmed by NMR as 3-hydroxy-β-ionone and the minor product co-eluted with 4-hydroxy-β-ionone produced by CYP101C1 and CYP102A1.9,16,35 The ratio of 3-hydroxy- to 4-hydroxy-β-damascone is similar to that observed for the equivalent products of β-ionone suggesting these two substrates are bound in a similar orientation in the active site with the C3 carbon being closest to the heme iron followed by C4.
When using the CYP101B1 whole-cell oxidation system we previously observed the formation of a blue coloured pigment in cells after induction, which can be the result of indigo formation from indole oxidation.9 Here we show that the addition of indole to this system resulted in increased levels of blue pigment (Fig. S1†). This indicates that indole and related aromatic molecules may be viable substrates for oxidation by CYP101B1. Therefore we tested a small selection of aromatic substrates for their ability to bind to and be oxidised by CYP101B1 (indole, p-cymene and phenylcyclohexane). p-Cymene resulted in virtually no shift to the high spin form compared to indole which induced a small spin state shift (15%). Phenylcyclohexane induced a slightly larger shift of 20% but bound to CYP101B1 with comparable affinity to β-damascone (Table 1, Fig. 3 and 4).
Both phenylcyclohexane and p-cymene induced moderate activity with CYP101B1 as measured by the NADH turnover frequency, but both are lower when compared to the norisoprenoids. Phenylcyclohexane was more active, 293 ± 9.0 min−1, than p-cymene, 197 ± 6.0 min−1. Indole did not induce any enzyme activity above that of the leak rate. HPLC and GC analysis of phenylcyclohexane turnover indicated that a single product was formed (Fig. 6). This was isolated using the whole-cell oxidation system, purified (∼10 mg) and characterised using NMR as trans-4-phenylcyclohexanol. Coupling of the NADH consumption to product formation was calculated to be 48% resulting in a product formation rate of 141 ± 17 min−1. Two products were observed in small quantities in the turnover of p-cymene. These were identified via GC and HPLC co-elution experiments as p-α,α-trimethylbenzyl alcohol (25%) and isopropylbenzyl alcohol (75%). The coupling of NADH to product formation and the product formation rates are low (13% and 25.3 ± 6.0 min−1).
Overall the data suggests certain aromatic substrates such as phenylcyclohexane are good substrates for CYP101B1 resulting in moderate activity and regioselective hydroxylations. However the norisoprenoid type substrates, ionone and damascone, show high affinities and activities for CYP101B1 and are efficiently oxidised with high regioselectivity and rates of product formation. α-Ionone is solely oxidised at the allylic C3 C–H compared to β-ionone and β-damascone where C3 oxidation dominates but with a contribution from allylic C4 hydroxylation.
Potassium ions have been shown to increase the strength of camphor binding to CYP101A1. The behaviour of CYP101 family enzyme from N. aromaticivorans is different due to changes in the residues that comprise the camphor binding site between the CYP101 family members. The K+-binding site in CYP101A1 is octahedral, with the backbone carbonyl oxygens of Glu84, Gly93, Glu94 and Tyr96, and two water molecules, as ligands.36 These residues are conserved in CYP101D2, but not in CYP101D1, CY101C1 and CYP101B1. As well as the change at Tyr96, Glu84, is replaced with Arg residues in CYP101B1 and CYP101D1 and Gln in CYP101C1 which significantly reduces the effect of potassium on substrate binding.14 We have previously reported that, in common with these other CYP101 enzymes, CYP101B1 shows decreased activity at higher ionic strength predominantly due to an increasing KM.14 The proximal face of CYP101B1 is likely to consist of more positively charged than negatively charged residues, similar to CYP101C1, CYP101D1 and CYP101D2, allowing for strong interactions with the negatively charged surface potential of Arx in the vicinity of the [2Fe–2S] cluster.14
The side-chains of Lys178, Asp182, Thr185, Arg186, Asp251, and Asn255 are involved in the extensive salt bridge interactions in camphor-bound CYP101A1.37,38 In addition the Asp97 (B′ helix) and Lys197 (G helix) interaction has been proposed to play a role in camphor entry and product exit. Structural analysis of CYP101D1, CYP101D2 and CYP101C1 revealed that this salt bridge network was less extensive in these CYP101 family members. In CYP101B1 the changes of Lys178 to an arginine and the more dramatic alteration of Asp97 to a glycine would also be expected to alter this salt bridge network. The modifications in these regions along with those in the active site will alter the substrate range and C–H bond oxidation selectively of CYP101B1, relative to the other CYP101 subfamily members. Overall these changes in CYP101B1, when compared with CYP101A1, are fewer than those of CYP101C1 which agrees with the observation that CYP101B1 retains some binding affinity for camphor while binding of this substrate is virtually abolished in CYP101C1.5,9
CYP101B1 is able to bind norisoprenoids and camphor as well as selected substituted aromatics. It is tempting to speculate that the proline residue one position to the C-terminus of the proximal iron binding cysteinate residue contributes to the seemingly broader substrate range of this enzyme, though this remains to be confirmed. This would agree with data from the Ile401Pro variant of CYP102A1 but the CYP121A1 enzyme from Mycobacterium tuberculosis which has a proline in the same position has been reported to have a very specific substrate and reaction type, the C–C coupling of a cyclodipeptide.39
The physiological substrate of CYP101B1 is yet to be confirmed. However, the norisoprenoid class of compounds such as β- and α-ionone and β-damascone are efficiently oxidised, the binding affinity for these substrates is relatively strong and they induce large type I spin state shifts. Norisoprenoids can be derived from the breakdown of carotenoids and N. aromaticivorans contains at least one carotenoid dioxygenase, though this was identified as a 13,14-dioxygenase.40 The relative location of the ketone and alkene functionalities in the butenone side chain has more of an effect on the binding affinity than the location of the alkene double bond in the cyclohexyl ring with α- and β-ionone being better substrates than β-damascone. The location of the double bond in the cyclohexyl ring does have an effect on the selectivity of the hydroxylation. The allylic C3 C–H bonds of α-ionone being exclusively oxidised but for β-ionone and β-damascone allylic C4 hydroxylation is observed alongside C3 oxidation. The C–H bonds of the C3 in the norisoprenoid substrates, which are oxidised in the major products, are presumably held closer to the heme iron than any other C–H bond. The lower regioselectivity of β-ionone and β-damascone is presumably due to the more reactive hydrogens on the allylic C4 being in close enough proximity to the heme iron enabling them to compete with oxidation at C3. In α-ionone the C4 hydrogen is vinylic and the reactive alkene double bond must be too far away from the heme iron for epoxidation to occur.
Microbial oxidative biotransformations of α- and β-ionones have been reported using different bacteria including Streptomyces.21 P450 enzymes from S. griseolus and S. griseus were shown to oxidise α- and β-ionone at the allylic C3 and C4 positions, respectively.41 Variants of the self-sufficient cytochrome P450, CYP102A1, are able to selectively oxidise β-ionone at the allylic C4 position.35 Others have reported the stereoselective oxidation of the different enantiomers of α-ionone at the C3 position.42 The best reported product formation rates and coupling of these variants of the highly efficient and self-sufficient CYP102A1 enzyme with norisoprenoids are 280 min−1 and 35% for β-ionone, 152 min−1 and 39% for (R)-α-ionone and 98 and 44% for (S)-α-ionone, which are lower than those observed with the CYP101B1/ArR/Arx system.42
Other researchers have reported the binding and oxidation of ionones by the CYP264B1, CYP109B1 and CYP109D1 bacterial enzymes. The CYP109 enzymes were regioselective for oxidation at the allylic position while CYP264B1 oxidised both α- and β-ionones solely at C3.20,23,43 Non-physiological electron transfer partners were used for all three of these cytochrome P450 enzymes and the catalytic rates were slow compared to those reported for CYP102A1 and CYP101B1 (<5 min−1).
CYP101C1 from the same Novosphingobium bacterium is interesting in that it can rapidly oxidise β-ionone (PFR ∼2000 min−1) to a mixture of 4-hydroxy- (75%) and 3-hydroxy-β-ionone (25%) but the turnover with α-ionone is much slower (PFR ∼ 220 min−1). CYP101C1 oxidation of α-ionone is also regioselective for allylic C3 oxidation. It can also oxidise β-damascone to a mixture of 3- and 4-hydroxy products as well as further oxidation products. CYP101B1 is more active than all the other P450s tested so far for α-ionone oxidation. Whether this reaction is stereoselective with enantiomerically pure α-ionones remains to be determined. With β-ionone and β-damascone CYP101B1 is less active than CYP101C1 but alters the selectivity to favour oxidation predominantly at the C3 rather than the allylic C4 position. Of the other enzymes utilised to oxidise ionones, CYP102A1 mutants and the CYP109 enzymes family are less active but are more selective for allylic hydroxylation at C4 while CYP264B1 is very selective for oxidation at C3 but has very low activity.
The selective oxidation of the unrelated and more hydrophobic aromatic substrate phenylcyclohexane to trans-4-phenylcyclohexanol by CYP101B1 on the aliphatic ring at moderate activity is interesting. Tyrosine variants of CYP101A1 have been used to oxidise phenylcyclohexane at predominantly the 4-position though the activities were lower. Other mutations in the active site were also shown to alter the selectivity though none were completely regioselective.44–46 CYP108D1, also from N. aromaticivorans, has also been reported to tightly bind phenylcyclohexane and oxidise it to trans-4-phenylcyclohexanol, though the activity using non-physiological electron transfer partners is very low.47 The oxidation of p-cymene was significantly slower than phenylcyclohexane suggesting the smaller size of this substrates is detrimental to binding and productive and selective turnover. p-Cymene can bind with either methyl or isopropyl substituent closest to the heme iron. The combination of selective oxidations, high monooxygenase activities and the large total number of turnovers obtained using the whole-cell oxidation system make CYP101B1 ideally suited to the biocatalytic generation of hydroxylated products.
To determine the dissociation constant the CYP101B1 enzyme was prepared as above in a total volume 2.5 mL and used to baseline the spectrophotometer. Aliquots of substrate (0.5–2 μL) were added using a Hamilton syringe from a 1, 10 or 100 mM stock solution in ethanol. The solution was mixed and the peak-to-trough difference in absorbance recorded between 700 nm and 250 nm. Further aliquots of substrate were added until the peak-to-trough difference did not shift further. The apparent dissociation constants, Kd, were obtained by fitting the peak-to-trough difference against substrate concentration to a hyperbolic function (eqn (1)):
![]() | (1) |
Both ionone substrates exhibited tight binding, with Kd < 5[E], and in these instances the data were fitted to the tight binding quadratic equation (eqn (2)):
![]() | (2) |
NADH turnover assays were performed with mixtures (1.2 mL) containing 50 mM Tris, pH 7.4, 0.5 μM CYP101B1, 5 μM Arx, 1 μM ArR and 100 μg mL−1 bovine liver catalase. The mixtures were oxygenated and then equilibrated at 30 °C for 2 min. Substrates were added as 100 mM stock solutions in ethanol to a final concentration of 1 mM. NADH was added to ∼320 μM, final A340 = 2.00, and the absorbance at 340 nm was monitored. The rate of NADH turnover was calculated using ε340 = 6.22 mM−1 cm−1 (Table 1).
GC-MS analyses were carried out with the injector held at 250 °C and the interface at 280 °C. The oven temperature was varied as follows: for norisoprenoids and phenylcyclohexane, 120 °C for 3 minutes followed by a gradient of 10 °C min−1 for 10 min and held at 220 °C for a further 7 minutes; p-cymene 80 °C for 3 minutes followed by a gradient of 10 °C min−1 for 14 min and held at 220 °C for a further 3 minutes. The GC retention times were as follows: β-ionone, 7.5 min; 4-hydroxy-β-ionone, 9.9 min; 3-hydroxy-β-ionone, 10.3 min; α-ionone, 6.7 min; trans-3-hydroxy-α-ionone, 9.4 min; cis-3-hydroxy-α-ionone, 9.5 min; 3-oxo-α-ionone, 9.8 min; β-damascone, 6.5 min; 4-hydroxy-β-damascone, 9.1 min; 3-hydroxy-β-damascone, 9.3 min; phenylcyclohexane, 5.1 min; trans-4-hydroxyphenylcyclohexane, 7.7 min; p-cymene 4.4 min; isopropylbenzyl alcohol, 8.9 min (14.0 after derivatisation); p-α,α-trimethylbenzylalcohol, 7.2 min; thymol, 8.9 min (11.9 after derivatisation); carvacrol, 9.1 min.
To isolate and identify products where standards were not available a whole-cell oxidation system utilising the plasmids pETDuetArx/ArR and pRSFDuetArx/CYP101B1 was used to oxidise β-ionone, α-ionone, β-damascone and phenylcyclohexane, as described previously.9 The supernatant (200 mL) was extracted in ethylacetate (3 × 100 mL), washed with brine (100 mL) and dried with magnesium sulphate and the organic extracts were pooled and the solvent was removed by vacuum distillation and then under a stream of nitrogen. The products were purified using silica gel chromatography using a hexane–ethyl acetate stepwise gradient ranging from 80:
20 to 50
:
50 hexane to ethyl acetate using 2.5% increases every 100 mL. The composition of the fractions was assessed by TLC and GC-MS and those containing single products (≥95%) were combined for characterisation. The solvent was removed under reduced pressure. The purified product was dissolved in CDCl3 and the organics characterised by NMR spectroscopy and GC-MS. NMR spectra were acquired on a Varian Inova-600 and spectrometer operating at 600 MHz for 1H and 150 MHz for 13C.
A combination of 1H, 13C, COSY, HSQC and HMBC experiments was used to determine the structures of the products. The cis/trans configuration of 3-hydroxy-α-ionone diastereomers was attributed from the assignments of Yamazaki52 and the selective nOe difference measurements in ROESY experiments. The diagnostic nOe was between the H3 and H6 protons in the cis diastereomer with the C11 methyl 3H signal.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4ra14010a |
This journal is © The Royal Society of Chemistry 2015 |