- Open Access
Biochemical characterization of three new α-olefin-producing P450 fatty acid decarboxylases with a halophilic property
Biotechnology for Biofuels volume 12, Article number: 79 (2019)
The CYP152 family member OleTJE from Jeotgalicoccus sp. ATCC 8456 has been well-known to catalyze the unusual one-step decarboxylation of free fatty acids towards the formation of terminal alkenes. Efforts to tune up its decarboxylation activity for better production of biological alkenes have been extensively explored via approaches such as site-directed mutagenesis and electron source engineering, but with limited success. To gain more insights into the decarboxylation mechanism and reaction bifurcation (decarboxylation versus hydroxylation), we turned to an alternative approach to explore the natural CYP152 resources for a better variety of enzyme candidates.
We biochemically characterized three new P450 fatty acid decarboxylases including OleTJH, OleTSQ and OleTSA, with respect to their substrate specificity, steady-state kinetics, and salt effects. These enzymes all act as an OleTJE-like fatty acid decarboxylase being able to decarboxylate a range of straight-chain saturated fatty acids (C8–C20) to various degrees. Site-directed mutagenesis analysis to the lower activity P450 enzyme OleTSA revealed a number of key amino acid residues within the substrate-binding pocket (T47F, I177L, V319A and L405I) that are important for delicate substrate positioning of different chain-length fatty acids and thus the decarboxylation versus hydroxylation chemoselectivity, in particular for the mid-chain fatty acids (C8–C12). In addition, the three new decarboxylases exhibited optimal catalytic activity and stability at a salt concentration of 0.5 M, and were thus classified as moderate halophilic enzymes.
The P450 fatty acid decarboxylases OleTJE, OleTJH, OleTSQ and OleTSA belong to a novel group of moderate halophilic P450 enzymes. OleTJH from Jeotgalicoccus halophilus shows the decarboxylation activity, kinetic parameters, as well as salt tolerance and stability that are comparable to OleTJE. Site-directed mutagenesis of several key amino acid residues near substrate-binding pocket provides important guidance for further engineering of these P450 fatty acid decarboxylases that hold promising application potential for production of α-olefin biohydrocarbons.
Development of renewable, sustainable, and cost-effective biofuels has been driven by the shortage of fossil fuels, serious environmental problems, and ever-changing geopolitical factors [1, 2]. Among different types of biofuels, biohydrocarbons have attracted much attention as an ideal alternative to petroleum-based fuels due to their high energy content, low hygroscopicity, and compatibility with existing engine and distribution systems [3,4,5,6,7].
Cytochrome P450 (CYP) enzymes are a superfamily of proteins with a thiolate-heme prosthetic group, which are renowned for their versatile catalytic activities and exceptional capability to accept a vast variety of substrates [8,9,10]. According to the catalytic properties, P450 enzymes can be classified into monooxygenases, peroxidases, and peroxygenases . The biofuel-related P450 fatty acid decarboxylases (FADCs), such as OleTJE from Jeotgalicoccus sp. ATCC 8456  and CYP-Sm46 from Staphylococcus massiliensis S46 , belong to the CYP152 peroxygenase family. This family of P450 enzymes utilizes H2O2, instead of O2 which is employed by most of P450 monooxygenases, as the oxidant to support the unique oxidative decarboxylation reactions (Fig. 1) that convert the Cn (n = 4–22) chain length free fatty acids (FFAs) into Cn−1 chain length 1-alkenes (i.e., α-olefins) [13, 14]. Since α-olefins are both excellent biofuel molecules and useful precursors of lubricants, detergents and other chemicals [15, 16], P450 FADCs hold promising application potential for production of biological α-olefins.
A majority of CYP152 peroxygenases catalyze FFA hydroxylation and decarboxylation reactions simultaneously. According to the main catalytic reaction types, this family of P450 enzymes can be classified into P450BSβ-like FFA hydroxylases (i.e., those with FFA hydroxylation as major chemistry, such as P450BSβ mainly produces similar amounts of α- and β-hydroxyl fatty acids, P450SPα and CYP-Aa162 mainly catalyze the Cα hydroxylation of FFAs, CYP-MP introduces the hydroxyl group at various carbon positions, and OleTMC prefers to hydroxylate the long-chain fatty acids) and OleTJE-like FFA decarboxylases (i.e., those preferentially catalyze FFA decarboxylation reaction, like OleTJE and CYP-Sm46Δ29) [11, 12, 14, 17,18,19]. To understand the unusual decarboxylation mechanism for activity and selectivity optimization, a growing number of studies on the prototypic P450 FADC OleTJE have been carried out [20,21,22,23,24,25,26,27,28,29]. For example, Munro et al. resolved the crystal structures of wild-type OleTJE  and a number of mutants such as R245L/E, F79A/Y/W and H85Q , revealing the important roles of Arg245, His85, and Phe79 in both catalytic activity and decarboxylation/hydroxylation bifurcation, as well as a group of active site residues responsible for productive fatty acid substrate binding. Furthermore, systematic mutagenesis analyses of a select group of active site residues including Arg245, His85 and Ile170 suggested the accurate substrate positioning is essential for decarboxylation activity [19, 27, 30, 31].
Interestingly, OleTJE, the first identified P450 FADC remains the best biocatalyst for α-olefin production in terms of both catalytic efficiency and chemoselectivity (i.e., decarboxylation versus hydroxylation) when compared to other biochemically characterized P450 FADCs including P450BSβ, CYP-MP, OleTMC, CYP-Aa162, and CYP-Sm46Δ29 [12, 19, 32]. Of note, some approaches such as redox partner engineering  and development of photocatalytic systems [33, 34] were also unsuccessful to improve the decarboxylation activity of OleTJE. Thus, it is highly expected to discover or engineer a novel CYP152 biocatalyst that can convert FFAs to 1-alkenes more efficiently and selectively.
In this work, we in vitro characterized three novel CYP152 FADCs including OleTJH from Jeotgalicoccus halophilus (CYP152L1_ortholog, GenBank accession number: WP_092595307), OleTSQ from Salinicoccus qingdaonensis (CYP152L8, WP_092983663), and OleTSA from Staphylococcus aureus (CYP152L7, WP_049319149). We determined the substrate preference, kinetic parameters and salt tolerance of these enzymes for the first time. Moreover, site-directed mutagenesis on OleTSA, the OleT with relatively lower activity, was performed in order to further understand the function of active site residues in the catalytic activity and chemoselectivity of P450 FADCs.
Genome mining of OleTJE-like P450 fatty acid decarboxylases
The number of biochemically characterized P450 FADCs is far less than that of existing CYP152 sequences with potential FFA oxidation activities in GenBank. To explore more P450 FADCs that may possess greater decarboxylation activity and/or selectivity, we built a phylogenetic tree (Fig. 2) based on the protein sequences of the two known P450 FADCs OleTJE (CYP152L1) and CYP-Sm46Δ29 (CYP152L2), and their homologous sequences with a sequence identity higher than 60%. Interestingly, these sequences mostly originate from the genera of Jeotgalicoccus, Staphylococcus, and Salinicoccus, which are well-known microorganisms associated with halophilicity or salt tolerance [35,36,37,38]. This finding suggests a possibility of identifying more 1-alkene producers and FADCs with higher decarboxylation activities from these groups of microorganisms.
Next, based on the conserved residues that are believed to be essential for the decarboxylation activity, including Phe79, His85, Ile170 and Arg245 (OleTJE numbering, Additional file 1: Figure S1), we from each major branch selected OleTJH from Jeotgalicoccus halophilus (CYP152L1_ortholog, GenBank accession number: WP_092595307), OleTSQ from Salinicoccus qingdaonensis (CYP152L8, WP_092983663), and OleTSA from Staphylococcus aureus (CYP152L7, WP_049319149) for our following biochemical characterization. According to the protein sequence alignment (Additional file 1: Table S1 and Additional file 2), OleTJH, OleTSQ, and OleTSA show 93% (62%), 76% (63%), and 69% (64%) amino acid sequence identity to OleTJE (CYP-Sm46Δ29), respectively.
Substrate specificity and chemoselectivity
The codon-optimized genes that encode OleTJH, OleTSQ, and OleTSA were individually expressed in Escherichia coli BL21(DE3). The resultant N-terminal His6-tagged recombinant proteins were purified to homogeneity using nickel affinity chromatography (Additional file 1: Figure S2). As expected, the three enzymes showed characteristic CO-bound reduced difference spectra (Additional file 1: Figure S3), indicative of their functional expression.
Using purified OleT enzymes, we determined their activities towards a range of straight-chain saturated fatty acids (C8–C20) with H2O2 as cofactor. As results, all the four enzymes exhibited similar substrate preference profiles with decanoic acid (C10) or lauric acid (C12) as their optimal substrate (Fig. 3). The lauric acid conversion ratios were 93.8 ± 6.1%, 98.6 ± 0.6%, 99.1 ± 0.2%, and 86.2 ± 1.9% for OleTJE, OleTJH, OleTSQ, and OleTSA, respectively. With respect to the productivity of 1-alkenes, lauric acid was the best substrate unanimously. Of note, OleTJE, OleTJH, and OleTSQ exhibited higher decarboxylation (DC) activities than hydroxylation (HD) activities with the DC/HD values up to 38.3 when C10–C14 FFAs were used as substrates. However, OleTSA only showed moderately higher alkene production from its optimal substrates lauric acid (DC/HD = 2.6) and myristic acid (DC/HD = 4.0). By contrast, such chemoselectivity was entirely reversed when the long-chain FFA arachidic acid (C20) was the substrate, suggesting that substrates with different carbon-chain length could adopt distinct substrate-binding modes, thereby impacting the preference for decarboxylation or hydroxylation.
Subsequently, we elected to perform more detailed analyses of the myristic acid conversions catalyzed by all four OleT enzymes for two reasons: (1) myristic acid gave the highest DC/HD values among the tested FFAs, thus being the optimal decarboxylation substrate; (2) the commercial availability of α- and β-hydroxy myristic acid would enable quantification of the hydroxylation products. As expected, OleTJE, OleTJH, OleTSQ, and OleTSA all produced 1-tridecene as the dominant product. As for hydroxylation products, β-OH-myristic acid was the major hydroxylation product (6.8–15.4% of the total products), α-OH-myristic acid only accounted for less than 1% of the total products, and no γ-, δ-, or ε-hydroxylated products were detected (Table 1). Taken together, OleTJH, OleTSQ, and OleTSA are three new OleTJE-like FADCs, which can sometimes outperform the decarboxylation activity of OleTJE, depending on the FFA substrates to be decarboxylated.
To further quantitatively evaluate the catalytic efficiencies of OleTJH, OleTSQ, and OleTSA, we determined their steady-state kinetic parameters (Table 2) towards the optimal substrate lauric acid (C12) with OleTJE as a control (Additional file 1: Figure S4). By measuring the initial substrate consumption rates with gas chromatography, we found that the kinetic parameters of OleTJH were mostly similar to those of the control FADC OleTJE with a slightly higher catalytic activity (kcat). In comparison, OleTSQ and OleTSA displayed relatively lower catalytic efficiency (kcat/Km) due to their attenuated kcat values, albeit a slightly higher lauric acid substrate-binding affinity reflected by their lower Km. This may seem inconsistent with the previous results that OleTSQ showed higher conversion rate (99.1 ± 0.2%) towards lauric acid than OleTJE (93.8 ± 6.1%). We reason this contradiction might be that OleTSQ could be more likely prone to H2O2 inactivation when used at low enzyme concentrations as in the kinetic studies.
Next, we evaluated the reaction turnover numbers (TONs) of the four FADCs at a higher substrate concentration (1 mM) for longer reaction time (12 h). A previous study showed that the fed-batch addition of H2O2 reduced enzyme inactivation and improved catalytic conversion . We thus performed the assay with one-off and fed-batch supplementation of H2O2, respectively. As expected, the fed-batch addition of H2O2 resulted in higher TONs in general (Table 2). The maximum TON of 185 was obtained from OleTJE while OleTJH showed the best TON in 1 mM H2O2. OleTJE and OleTJH seemed to be relatively more resistant to H2O2 inactivation than OleTSA and OleTSQ, as more than twofold (2.5 and 2.7) TON reductions were observed in H2O2 batch addition assays of the latter two enzymes compared with their fed-batch addition assays. Nevertheless, the three new FADCs exhibited strong decarboxylation potentials and similar kinetic parameters to OleTJE. Future reaction process optimization such as in situ generation of H2O2 [33, 34] or utilization of a redox cascade system [22, 26, 31] could be considered to attenuate the catalyst inactivation by high concentration H2O2.
Effects of salt concentration on OleT enzymes
Considering the facts that precipitation of OleTJE occurs in low salt buffer  and that the OleT enzyme host strains including Jeotgalicoccus, Salinicoccus, and Staphylococcus species are often associated with the property of halophilicity or salt tolerance [35,36,37,38], we investigated the effects of salt concentration on the activity and stability of the OleT enzymes. Experimentally, the activities of OleTJE, OleTJH, OleTSQ, and OleTSA against lauric acid were evaluated under a range of NaCl concentrations. The overall and decarboxylation activities of all four enzymes remained almost unchanged (> 90%) within up to 2 M salt concentration (Fig. 4). However, when NaCl concentration was 3 M or higher, the activities of all enzymes decreased gradually to different extents, with the only exception of OleTJE at 3 M NaCl. However, they all remained at least 22.7% activity in saturated salt concentration, even up to 66.8% for OleTJE.
Next, we examined the salt stability of these OleT enzymes. Upon a 12-h treatment with different salt solution buffers. It is evident that OleTJE and OleTJH showed much better stabilities than OleTSQ and OleTSA under all tested NaCl concentrations. Again, OleTJE and OleTJH were even able to retain a majority of overall and decarboxylation activities at the saturated concentration (Table 3). Overall, the majority of activity of all four enzymes was maintained in buffer with a NaCl concentration of 0.5–2 M. According to the most widely used standards for halophilicity by Kushner, i.e., extreme halophiles (best in 2.5–5.2 M salt), borderline extreme halophiles (best in 1.5–4.0 M salt), and moderate halophiles (best in 0.5–2.5 M salt) [39,40,41,42], which can also be applied to halophilic enzymes [43, 44], these P450 FADCs belong to moderate halophilic enzymes.
Site-directed mutagenesis of OleTSA fatty acid decarboxylase
Despite significant advances in understanding of the OleTJE decarboxylation mechanism [20, 23, 28, 29, 31, 45, 46], an effective strategy for engineering a better P450 FADC remains unclear. In this study, the four analogous OleT enzymes with characterized substrate specificity, product distribution and kinetic parameters provided an outstanding opportunity to further dissect the residues important for the decarboxylation activity. Specifically, we discovered, by comparing the crystal structure of OleTJE (PDB ID #: 4L40) and the three modeled structures of the new FADCs, that OleTJE, OleTJH, OleTSQ, and OleTSA have highly similar substrate-binding pockets with only five out of 19 residues that are different from one another. Comparatively, OleTJE and OleTSA have the most different substrate-binding pockets (Fig. 5 and Additional file 1: Table S2), which is consistent with their largest difference in catalytic activities among the four tested OleT enzymes towards mid-chain FFA substrates (C8–C12) (Fig. 3 and Table 2).
Following these analyses, we were drawn to mutate Thr47, Ile177, Val319, and Leu405 of OleTSA into their counterparts in OleTJE in consideration of the highest catalytic efficiency of OleTJE and the most different catalytic activities and substrate-binding pockets between these two FADCs (Additional file 1: Table S2). Specifically, we prepared four single mutants including OleTSA–T47F (M1), OleTSA–I177L (M2), OleTSA–V319A (M3) and OleTSA–L405I (M4), and one quadruple mutant OleTSA–T47F-I177L–V319A–L405I (M5). Measurements of the in vitro activities of these mutants towards C8–C12 FFAs (Fig. 6 and Additional file 1: Figure S7) showed that (1) for caprylic acid (C8), M1, M2, and M4 displayed enhanced overall and decarboxylation activities, while M3 and M5 exhibited decreased overall and decarboxylation activities when compared to the wild-type OleTSA; (2) for decanoic acid (C10), all four single mutants showed improved overall and decarboxylation activities relative to their parental enzyme, but the activity of the quadruple mutant decreased; and (3) for lauric acid (C12), all mutants demonstrated > 20% improvement of substrate conversion ratios compared to the wild-type OleTSA. However, except for M2, other mutations did not lead to higher levels of 1-undecene (C11) production.
Collectively, all four tested positions could become the “hot spots” for the future iterative saturation mutagenesis (ISM) and combinatorial active-site saturation testing (CASTing) analyses [47, 48]. Here, the attenuated activity of OleTSA–T47F–I177L–V319A–L405I suggests there might not be an expected synergistic effect between these four residues.
In this study, we have expressed, purified and characterized three new P450 FADCs in vitro, namely, OleTJH, OleTSQ, and OleTSA that can efficiently catalyze the conversions from FFAs to 1-alkenes. Among all substrates (C8–C20 FFAs), these three new enzymes exhibit higher decarboxylation activity towards lauric acid and myristic acid, with the lower decarboxylation activity for eicosanoic acid. The strict carbon chain length dependence for decarboxylative activity is consistent with the well-characterized FADCs OleTJE  and CYP-Sm46∆29 , which further emphasize the importance of substrate identity in achieving the desirable reaction type. Our finding that β-hydroxyl product rather than α-hydroxyl product is the major by-product in all alkene-generating catalysis here supports the notion that Cβ–H abstraction may be one of the key factors that affect the production of 1-alkenes, which is well in line with the proposed mechanism by Makris et al. that Compound I (an iron(IV)-oxo cation radical) abstracts the Cβ–H atom to initiate the decarboxylation of fatty acids .
In general, salt solutions could impact protein folding mode, and hence the conformational stability of proteins [49,50,51]. We found that OleTJH, OleTSQ, OleTSA and OleTJE exhibit optimal activities as well as stabilities in salt solutions with the NaCl concentration ranging from 0.5 to 2 M (2.9–11.7%), thus classifying these FADCs as moderate halophilic enzymes. This seems to be consistent with the halophilic property of their native host microorganisms which have shown optimal growth at the NaCl strength of 2–3% (Jeotgalicoccus halophilus) , 3% (Salinicoccus qingdaonensis)  and 2.1–7% (Staphylococcus aureus) , respectively, thus being classified as moderate halophilic bacterial strains.
However, not all halophiles bear halophilic (P450) proteins. Two other P450 enzymes CYP102A26  and CYP116B62  from the halophilic bacteria Pontibacillus halophilus and Halomonas sp. NCIMB 172 exhibited decreased enzymatic activity in the presence of NaCl at 0.025 M and 0.05 M, respectively. This difference may be related to the different strategies how halophiles adapt to high salinity (“salt in” or “salt out”). Most halophilic bacteria use a “salt out” strategy to exclude salt to protect their non-halophilic enzymes . While in a ‘salt in’ strategy that is used primarily by haloarchaea to accumulate high concentrations of salt, their enzymes have been examined to tolerate or even require 4–5 M salt .
Previous studies suggested that hydrogen bonds between negatively charged side chains and water molecules are critical for halophilic proteins to maintain a stable hydration shell under water-limited conditions [55, 56]. In a detailed study comparing 15 pairs of homologous non-halophilic and halophilic proteins (γ-glutamyltranspeptidases), it was shown that halophilic proteins tend to have increased number of acidic amino acids in their total surface residues, and within the whole protein sequences . Thus, the ratio of acidic amino acids (Glu/Asp) in the total protein sequence of a halophilic protein has been considered as an important parameter to distinguish it from non-halophilic proteins from the same class of enzymes [39,40,41]. Notably, this is the first time that ‘halophilic P450 enzymes’ has been proposed for the need of salt to maintain the enzyme activity and stability. We subsequently calculated the ratios of acidic amino acids (Glu/Asp) in the total protein sequences for the four P450 halophilic FADCs. An average ratio of 14% (OleTJE: 14.5%, OleTJH: 14.2%, OleTSQ: 14.6%, OleTSA: 14.5%) (Figure S8) was obtained, which is comparable to a halophilic alkaline phosphatase (TAP: 14.0%) from Halomonas sp. 593  and significantly higher than those non-halophilic P450 enzymes such as P450BSβ (11.7%) and P450SPα (11.8%) [14, 17]. However, detailed salt-tolerance mechanisms used by these P450s have yet to be elucidated. Furthermore, we suggest that a halophilic bacterial host for P450 FADCs could be engineered for efficient production of 1-alkenes based on the halophilic property of these CYP152 peroxygenases.
In our mutagenesis analyses, according to the catalytic activity profiles of the OleTSA mutants, T47F mutation seemed to show more impact on the conversion of longer-chain FFAs (a 1.3-fold increase toward C12) than on that of the shorter chain C8 fatty acid (1.2-fold increase). Considering this residue is located at the top (the alkyl end) of the binding pocket, the bulky nature of the Phe side chain may help push an improved docking and binding of the longer chain fatty acid substrates. However, it is too far to the active site to make an impact on the chemo-selectivity of the enzymes (Fig. 5 and Additional file 1: Figure S7). Another mutation I177L is located in a structurally disordered region, the F–G loop. This distal loop area has been recently studied by Makris et al. suggesting its necessity for both substrate positioning and product release . A single amino acid change at this site (L176G) switched OleTJE from a decarboxylase to a hydroxylase . Our results with the mutant OleTSA-I177L displayed improved catalytic conversion as well as the DC over HD selectivity towards all tested substrates, in particular the mid-chain fatty acid C12 (Additional file 1: Figure S7), demonstrating the crucial role of this leucine residue in OleT enzymes for the decarboxylation reaction. After all, the Ile to Leu replacement is only a minor structural change; the data suggest a delicate and accurate substrate positioning is required for regio- and chemo-selectivity of these P450 FADCs. As for the mutation L405I, since it is more into the middle and towards the bottom (the acyl end) of the binding pocket, the mutant OleTSA-L405I exhibited more prominent effect on the DC/HD selectivity towards the shorter-chain C8 fatty acid than by the other mutants (Additional file 1: Figure S7). However, despite the overall modulation of these mutants on the catalytic conversion and reaction selectivity of OleTSA, none of them are superior to OleTJE in catalytic activity or decarboxylation selectivity. It was also difficult to see a synergistic effect from the quadruple mutant (Fig. 6).
These results together with previous studies [12, 24,25,26, 31], suggest that it remains challenging to rationally engineer a better substrate-binding pocket in OleT enzymes for better overall and decarboxylation activities. The residues that are distant from the substrate-binding site may be worth more attention in the future, which apparently requires a high-throughput screening (HTS) assay to enable the directed evolution [59,60,61] efforts. The development of such a HTS method is currently ongoing in our laboratory. Better unnatural FADCs would not only benefit the mechanistic understanding, but also increase the opportunity for industrialization of the intriguing α-olefin-producing enzymes.
In this study, we biochemically characterized three new P450 FADCs that are able to efficiently decarboxylate a range of saturated fatty acids (C8–C20), with OleTJH showing the decarboxylation activities, kinetic parameters, as well as salt tolerance and stability that are comparable to OleTJE. All four tested P450 FADCs exhibited moderate halophilicity unanimously. Further mutagenesis analysis based on protein sequence and activity comparisons of OleTJE and OleTSA provides more insights into the unique catalytic mechanism of P450 FADCs. The results lay out an advanced foundation for the future engineering efforts on this class of moderate halophilic P450 enzymes in order to produce desirable carbon chain length 1-alkenes efficiently.
The strains of E. coli DH5α and BL21(DE3) were preserved by our laboratory. All chemicals and antibiotics were obtained from TCI (Shanghai, China), Solarbio (Beijing, China), Sigma Aldrich (St. Louis, MO, USA) or Thermo Scientific (Shanghai, China). I-5™ 2× High-Fidelity Master Mix and Trelief SoSoo Cloning Kit were obtained from TsingKe (Beijing, China). Plasmid Miniprep Kit from TsingKe was used to prepare plasmid DNA from E. coli DH5α. ClonExpress II One-Step Cloning Kit was purchased from Vazyme (Nanjing, China). The 10× QuickRun™ Fast Running Buffer and FlexiRun™ premixed gel solution for SDS-PAGE were obtained from MDBio (Xinbei, China). Ni–NTA resin used for protein purification was supplied by Sangon Biotech (Shanghai, China). PD-10 desalting columns were purchased from GE Healthcare (Piscataway, NJ, USA). Millipore Amicon Ultra centrifugal filters were obtained from Millipore (Billerica, MA, USA).
The CYP152 family protein sequences were BLAST searched against the protein database on National Center for Biotechnology Information (NCBI) based on the query sequences of OleTJE (CYP152L1) and CYP-Sm46Δ29 (CYP152L2). Proteins with a sequence identity greater than 60% were selected and the sequence alignments were performed using ClustalW (https://www.genome.jp/tools-bin/clustalw). The phylogenetic tree was built using the neighbor-joining (NJ) method in MEGA 7.0 package. Bootstrap values shown next to the branches were computed from 1000 bootstrap tests.
Molecular cloning and protein purification
The gene sequences encoding OleTJH from Jeotgalicoccus halophilus (CYP152L1_ortholog, GenBank accession number: WP_092595307), OleTSQ from Salinicoccus qingdaonensis (CYP152L8, GenBank accession number: WP_092983663), and OleTSA from Staphylococcus aureus (CYP152L7, GenBank accession number: WP_049319149) were codon-optimized and synthesized by Qinglan (Yixing, China), and then cloned into the vector pET28b via the NdeI/XhoI restriction sites for expression of the N-terminal His6-tagged proteins. The generation of single mutation gene constructs of OleTSA–T47F, OleTSA–I177L, OleTSA–V319A, OleTSA–L405I was achieved by site-directed mutagenesis via overlap extension PCR . The sequences of primers used in this study are listed in Additional file 1: Table S1). All cloned sequences were confirmed by DNA sequencing at Sangon Biotech (Shanghai, China), and then used to transform E. coli BL21 (DE3) for protein expression and purification.
The E. coli BL21 (DE3) cells carrying the recombinant expression vector were grown 12 h at 37 °C with shaking at 220 rpm and then used as seed cultures to inoculate (1:100 ratio) a modified Terrific Broth medium containing 4% glycerol, a rare salt solution  and 1 mM thiamine. Cells were grown at 37 °C for 3–4 h until the optical density at 600 nm (OD600) reached 0.8–1.0, to which δ-aminolevulinic acid (5-ALA, 0.5 mM) and isopropyl β-d-1-thiogalactopyranoside (IPTG, 0.2 mM) were added, and followed by 24 h of cultivation at 18 °C. The cells were harvested (6000 ×g, 4 °C, 10 min) and then stored at − 80 °C.
Purification of His-tagged OleT enzymes was carried out as described by Liu et al.  with slight modifications. The cell pellets stored at − 80 °C were taken out to melt at room temperature, then all the following steps were performed at 4 °C. Briefly, the cell pellets were re-suspended in 50 mL lysis buffer (pH 8.0, 50 mM NaH2PO4, 300 mM NaCl, 10% glycerol, 10 mM imidazole) through vortexing and were then disrupted by ultrasonication. Cell-free lysate was obtained by centrifuging at 10,000 ×g for 1 h at 4 °C, to which 1–2 mL Ni–NTA resin slurry was added and mixed gently for 2 h. The mixture was loaded onto an empty column and washed with approximately 200 mL wash buffer (pH 8.0, 50 mM NaH2PO4, 300 mM NaCl, 10% glycerol, 20 mM imidazole) until no protein was eluted in flow-through. His-tagged proteins bound to Ni–NTA resin were eluted with 10 mL elution buffer (pH 8.0, 50 mM NaH2PO4, 500 mM NaCl, 10% glycerol, 250 mM imidazole). The eluents were concentrated with an Amicon Ultra centrifugal filter (30 kDa cutoff) and were buffer-exchanged into storage buffer (pH 7.4, 50 mM NaH2PO4, 500 mM NaCl, 10% glycerol). The final purified proteins were flash frozen by liquid nitrogen and stored at − 80 °C for later use.
UV–visible spectroscopic characterization and determination of enzyme concentration
The analysis of the UV–visible spectroscopic properties was carried out as described by Xu et al. with minor modifications . In general, the UV–visible spectroscopic properties of the His-tagged proteins were performed on a Cary 60 UV–visible spectrophotometer (Varian, UK). For preparation of the dithionite-reduced ferrous-CO complex of each enzyme, the purified ferric enzymes were diluted in Storage buffer (pH 7.4, 50 mM NaH2PO4, 500 mM NaCl, 10% glycerol) and subjected to slow CO bubbling for 40–50 s for initial scan at 350–600 nm to record the spectrum of the resting-state CO-bound proteins, and then followed by sufficient reduction of the protein by sodium dithionate to obtain the CO-bound reduced difference spectrum. The protein concentration was calculated by the reduced differential extinction coefficient ε450–490 nm of 91,000 M−1 cm−1 for the functional P450 concentration .
In vitro enzymatic assay
Typical assays containing 1 μM of each enzyme (OleTJE, OleTJH, OleTSQ, OleTSA or its individual mutant), 200 μM fatty acid substrate (C8–C20 FFAs prepared from stock solution in DMSO (20 mM)), and 220 μM H2O2 in 200 μL of buffer (pH 7.4, 50 mM NaH2PO4, 500 mM NaCl, 10% glycerol) were carried out at 30 °C for 2 h. This amount of hydrogen peroxide was proven to be sufficient (Additional file 1: Figure S9). Reactions were quenched by adding 20 μL 10 M HCl, then heptadecanoic acid (C17) was added as internal standard and the mixture was extracted by 150 μL ethyl acetate. The organic phase was analyzed by gas chromatography (GC) as described below.
For detection of 1-heptene (C7) product generated from caprylic acid (C8) fatty acid decarboxylation, 1.5 mL polytetrafluorethylene (PTFE) septum-sealed glass bottles were used for 500 μL reaction systems containing 1 μM enzyme, 200 μM C8 fatty acid substrate, 220 μM H2O2, and 200 μM 1-nonene (C9) as the internal standard. The reactions were incubated at 30 °C for 2 h with shaking at 100 rpm. Then the reactions were placed at 4 °C for 12 h to stop reactions prior to heating at 40 °C for 20 min for headspace sampling using a gas-tight Hamilton syringe for GC–MS analysis. Different concentrations of the authentic 1-heptene standard incubated under the same conditions as reactions were analyzed using the same GC–MS method to obtain the standard curve. After headspace sampling, the remaining C8 substrate in reactions was extracted as above described for further GC measurement and analysis.
Effects of salt concentration
Effect of salt (NaCl) was tested by measuring the decarboxylation activity in the reaction buffer (pH 7.4, 50 mM NaH2PO4, 10% glycerol) containing different concentrations of NaCl ranging from 0 M, 0.5 M, 1 M, 2 M, 3 M to saturated salinity (> 5 M) at 30 °C for 2 h. To determine their halostability, these enzymes were pre-incubated in buffer (pH 7.4, 50 mM NaH2PO4, 10% glycerol) containing different concentrations of NaCl (0 M, 0.5 M, 2 M, and saturated salt solution) at 30 °C for 12 h with non-treated enzyme as the control. The residual activities were then measured as in typical assays. For all reactions in salt effect experiments, 1–3 μL of enzymes stored at − 80 °C in storage buffer (pH 7.4, 50 mM NaH2PO4, 500 mM NaCl, 10% glycerol) were used in 200 μL reaction systems containing different concentrations of salt.
Steady-state kinetic analysis
Determination of the steady-stated kinetic parameters was carried out as described by Xu et al. with minor modifications . Briefly, 15–40 nM enzyme was used with a range of different substrate (C12) concentrations in a 1 mL reaction system (pH 7.4, 50 mM NaH2PO4, 500 mM NaCl, 10% glycerol). The reaction was initiated by adding an excess amount of H2O2 (330 μM) at 30 °C. Aliquots (200 μL) of reactions were removed and quenched at fixed time points (0, 0.5, 1, 2 min) by adding 20 μL 10 M HCl for reaction termination. Sample extraction and analysis were performed the same as above for GC analysis. Initial rates were calculated in terms of the substrate consumption by each enzyme. Kinetic analyses were performed using OriginPro 8.5 program.
The 3D structures of OleTJH, OleTSQ and OleTSA were generated by using the homology modeling function of Phyre2 (Protein Homology/analogy Recognition Engine V2.0) program. The structure of OleTJE (PDB ID #: 4L40) and the modeled structures were analyzed with PyMOL (V2.2.0).
The hydrocarbon and fatty acid samples were analyzed by the methods from Guan et al. . The Agilent 7890B gas chromatograph equipped with a capillary column HP-5 (Agilent Technologies, Santa Clara, CA, USA; cross-linked polyethylene glycerol, i.d. 0.25 μm film thickness, 30 m by 0.32 mm) was used for analyses. The flow rate of helium was set to 1 mL per min. The oven program was set initially at 40 °C for 4 min, then increased to 280 °C by the rate of 10 °C per min and held for 5 min. The injecting temperature was set to 280 °C under splitless injection conditions with 1 μL injection volume. Under these conditions, the retention times and signal intensity of fatty acids and terminal alkenes products were determined and compared with corresponding authentic standards (FFAs: C8–C20, 1-alkenes: C7–C19) and the internal standard (1-heptadecanoic acid (C17)). For analyses of GC–MS, the gas chromatography was equipped with an Agilent 5975C MSD single quadrupole mass spectrometer operated under electron ionization mode at 70 eV in the scan range of 50–500 m/z. For detection of α- and β-hydroxyl products of myristic acid, the samples extracted from the myristic acid (C14) reactions were derivatized with an equal volume of N,O-bis (trimethylsilyl) trifluoroacetamide (BSTFA) with 1% trimethylchlorosilane at 72 °C for 15 min before GC–MS analysis. The GC–MS analysis used the previous protocol adapted from Rude et al. with the Agilent J&W DB-5MS column (i.d. 0.25 μm film thickness, 50 m by 0.25 mm). Peak identity was determined by comparisons of the retention time and fragmentation pattern with the authentic standard compounds where available and to National Institute of Standards and Technology, USA mass spectral database. We found that the sum of all products (α-, β-hydroxy myristic acids and 1-tridecene) almost accounted for 99% of the substrate consumption, so we subtracted the 1-alkene production from the total substrate consumption to quantify the percentage of hydroxylated products for all substrates unless otherwise stated. For detection of 1-heptene (C7) product produced by caprylic acid (C8) decarboxylation, 500 μL of the reaction headspace sample was injected into GC–MS system with a Hamilton needle syringe and analyzed by authentic standard curves as well as reaction controls. The oven temperature procedure was set at follows: 40 °C for 2 min, then increased to 100 °C by 5 °C min−1 and held for 2 min.
- CYP or P450:
cytochrome P450 enzyme
free fatty acid
fatty acid decarboxylases
polymerase chain reaction
gas chromatography–mass spectrometry
Stephanopoulos G. Challenges in engineering microbes for biofuels production. Science. 2007;315:801–4.
Shafiee S, Topal E. When will fossil fuel reserves be diminished? Energy Policy. 2009;37:181–9.
Liao W, Heijungs R, Huppes G. Natural resource demand of global biofuels in the Anthropocene: a review. Renew Sustain Energy Rev. 2012;16:996–1003.
Dutta K, Daverey A, Lin J-G. Evolution retrospective for alternative fuels: first to fourth generation. Renew Energy. 2014;69:114–22.
Schirmer A, Rude MA, Li X, Popova E, del Cardayre SB. Microbial biosynthesis of alkanes. Science. 2010;329:559–62.
Herman NA, Zhang W. Enzymes for fatty acid-based hydrocarbon biosynthesis. Curr Opin Chem Biol. 2016;35:22–8.
Lin FM, Marsh ENG, Lin XN. Recent progress in hydrocarbon biofuel synthesis: pathways and enzymes. Chin Chem Lett. 2015;26:431–4.
Denisov IG, Makris TM, Sligar SG, Ilme S. Structure and chemistry of cytochrome P450. Chem Rev. 2005;105:2253–77.
Coon MJ. Cytochrome P450: nature’s most versatile biological catalyst. Annu Rev Pharmacol Toxicol. 2005;45:1–25.
Hrycay EG, Bandiera SM. The monooxygenase, peroxidase, and peroxygenase properties of cytochrome P450. Arch Biochem Biophys. 2012;522:71–89.
Rude MA, Baron TS, Shane B, Murtaza A, Cardayre SB, Del Andreas S. Terminal olefin (1-alkene) biosynthesis by a novel p450 fatty acid decarboxylase from Jeotgalicoccus species. Appl Environ Microbiol. 2011;77:1718.
Xu H, Ning L, Yang W, Fang B, Wang C, Wang Y, Xu J, Collin S, Laeuffer F, Fourage L, Li S. In vitro oxidative decarboxylation of free fatty acids to terminal alkenes by two new P450 peroxygenases. Biotechnol Biofuels. 2017;10:208.
Hrycay EG, Bandiera SM. Monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 enzymes. Oxygen Transport Tissue XXXIII. 2016;851:1–61.
Matsunaga I, Sumimoto T, Ueda A, Kusunose E, Ichihara K. Fatty acid-specific, regiospecific, and stereospecific hydroxylation by cytochrome P450 (CYP152B1) from Sphingomonas paucimobilis: substrate structure required for alpha-hydroxylation. Lipids. 2000;35:365–71.
Rui Z, Li X, Zhu X, Liu J, Domigan B, Barr I, Cate JH, Zhang W. Microbial biosynthesis of medium-chain 1-alkenes by a nonheme iron oxidase. Proc Natl Acad Sci USA. 2014;111:18237–42.
Rui Z, Harris NC, Zhu X, Huang W, Zhang W. Discovery of a family of desaturase-like enzymes for 1-alkene biosynthesis. Acs Catal. 2015;5:7091–4.
Shoji O, Fujishiro T, Nakajima H, Kim M, Nagano S, Shiro Y, Watanabe Y. Hydrogen peroxide dependent monooxygenations by tricking the substrate recognition of cytochrome P450BSβ. Angew Chem Int Ed Engl. 2007;119:3730–3.
Girhard M, Kunigk E, Tihovsky S, Shumyantseva VV, Urlacher VB. Light-driven biocatalysis with cytochrome P450 peroxygenases. Biotechnol Appl Biochem. 2013;60:111–8.
Lee JW, Niraula NP, Trinh CT. Harnessing a P450 fatty acid decarboxylase from Macrococcus caseolyticus for microbial biosynthesis of odd chain terminal alkenes. Metab Eng Commun. 2018;7:e00076.
Grant JL, Hsieh CH, Makris TM. Decarboxylation of fatty acids to terminal alkenes by cytochrome P450 compound I. J Am Chem Soc. 2015;137:4940–3.
Grant JL, Mitchell ME, Makris TM. Catalytic strategy for carbon–carbon bond scission by the cytochrome P450 OleT. Proc Natl Acad Sci USA. 2016;113:10049.
Dennig A, Kuhn M, Tassoti S, Thiessenhusen A, Gilch S, Bulter T, Haas T, Hall M, Faber K. Oxidative decarboxylation of short-chain fatty acids to 1-alkenes. Angew Chem Int Ed Engl. 2015;54:8819–22.
Hsieh CH, Huang X, Amaya JA, Rutland CD, Keys CL, Groves JT, Austin RN, Makris TM. The enigmatic P450 decarboxylase OleT is capable of, but evolved to frustrate, oxygen rebound chemistry. Biochemistry. 2017;56:3347–57.
Wang J, Lonsdale R, Reetz MT. Exploring substrate scope and stereoselectivity of P450 peroxygenase OleTJE in olefin-forming oxidative decarboxylation. Chem Commun. 2016;52:8131.
Hsieh CH, Makris TM. Expanding the substrate scope and reactivity of cytochrome P450 OleT. Biochem Biophys Res Commun. 2016;476:462–6.
Liu Y, Wang C, Yan J, Zhang W, Guan W, Lu X, Li S. Hydrogen peroxide-independent production of α -alkenes by OleT JE P450 fatty acid decarboxylase. Biotechnol Biofuels. 2014;7:28.
Faponle AS, Quesne MG, de Visser SP. Origin of the regioselective fatty acid hydroxylation versus decarboxylation by a cytochrome P450 peroxygenase: what drives the reaction to biofuel production? Chemistry. 2016;22:5478–83.
Belcher J, Mclean KJ, Matthews S, Woodward LS, Fisher K, Rigby SEJ, Nelson DR, Potts D, Baynham MT, Parker DA. Structure and biochemical properties of the alkene producing cytochrome P450 OleTJE (CYP152L1) from the Jeotgalicoccus sp. 8456 bacterium. J Biol Chem. 2014;289:6535–50.
Du J, Liu L, Guo L, Yao X, Yang J. Molecular basis of P450 OleTJE: an investigation of substrate binding mechanism and major pathways. J Comput Aid Mol Des. 2017;31:483–95.
Matthews S, Belcher JD, Tee KL, Girvan HM, McLean KJ, Rigby SE, Levy CW, Leys D, Parker DA, Blankley RT, Munro AW. Catalytic determinants of alkene production by the cytochrome P450 peroxygenase OleTJE. J Biol Chem. 2017;292:5128–43.
Fang B, Xu H, Liu Y, Qi F, Zhang W, Chen H, Wang C, Wang Y, Yang W, Li S. Mutagenesis and redox partners analysis of the P450 fatty acid decarboxylase OleTJE. Sci Rep. 2017;7:44258.
Amaya JA, Rutland CD, Makris TM. Mixed regiospecificity compromises alkene synthesis by a cytochrome P450 peroxygenase from Methylobacterium populi. J Inorg Biochem. 2016;158:11–6.
Zachos I, Gassmeyer SK, Bauer D, Sieber V, Hollmann F, Kourist R. Photobiocatalytic decarboxylation for olefin synthesis. Chem Commun. 2015;51:1918–21.
Köninger K, Grote M, Zachos I, Hollmann F, Kourist R. Light-driven enzymatic decarboxylation. J Vis Exp Jove. 2016;111:53439.
Koujima I. Salt-tolerance mechanism of Staphylococcus aureus: role of proline and water in osmoregulation of S. aureus (author’s transl). Nippon Saikingaku Zasshi. 1978;33:643–9.
Chen Y, Cui X, Pukall R, Li H, Yang Y, Xu L, Wen M, Peng Q, Jiang C. Salinicoccus kunmingensis sp. nov., a moderately halophilic bacterium isolated from a salt mine in Yunnan, south-west China. Int J Syst Evol Microbiol. 2007;57:2327–32.
Qu Z, Li Z, Zhang X, Zhang X. Salinicoccus qingdaonensis sp. nov., isolated from coastal seawater during a bloom of green algae. Int J Syst Evol Micr. 2012;62:545.
Liu W, Jiang L, Guo C, Yang S. Jeotgalicoccus halophilus sp. nov., isolated from salt lakes. Int J Syst Evol Microbiol. 2011;61:1720–4.
Ma Y, Galinski EA, Grant WD, Oren A, Ventosa A. Halophiles 2010: life in saline environments. Appl Environ Microb. 2010;76:6971–81.
Madern D, Ebel C, Zaccai G. Halophilic adaptation of enzymes. Extremophiles. 2000;4:91–8.
DasSarma S, DasSarma P. Halophiles and their enzymes: negativity put to good use. Curr Opin Microbiol. 2015;25:120–6.
Ventosa A, Nieto JJ, Oren A. Biology of moderately halophilic aerobic bacteria. Microbiol Mol Biol Rev. 1998;62:504–44.
Khire JM. Production of moderately halophilic amylase by newly isolated Micrococcus sp. 4 from a salt-pan. Lett Appl Microbiol. 1994;19:210–2.
Miyashita Y, Ohmae E, Ikura T, Nakasone K, Katayanagi K. Halophilic mechanism of the enzymatic function of a moderately halophilic dihydrofolate reductase from Haloarcula japonica strain TR-1. Extremophiles. 2017;21:591–602.
Amaya JA, Rutland CD, Leschinsky N, Makris TM. A distal loop controls product release and chemo- and regioselectivity in cytochrome P450 decarboxylases. Biochemistry. 2017;57:344–53.
Pickl M, Kurakin S, Cantú Reinhard FG, Schmid P, Pöcheim A, Winkler CK, Kroutil W, Visser SP, Faber K. Mechanistic studies of fatty acid activation by CYP152 peroxygenases reveal unexpected desaturase activity. ACS Catal. 2018;9:565–77.
Reetz MT, José Daniel C. Iterative saturation mutagenesis (ISM) for rapid directed evolution of functional enzymes. Nat Protoc. 2007;2:891–903.
Reetz MT, Bocola M, Carballeira JD, Zha D, Vogel A. Expanding the range of substrate acceptance of enzymes: combinatorial active-site saturation test. Angew Chem Int Ed Engl. 2005;117:4192–6.
Christian JHB, Waltho JA. Solute concentrations within cells of halophilic and non-halophilic bacteria. BBA Bioenerg. 1962;65:506–8.
Bonneté F, Madern D, Zaccai G. Stability against denaturation mechanisms in halophilic malate dehydrogenase “adapt” to solvent conditions. J Mol Biol. 1994;244:436–47.
Ebel C, Faou P, Kernel B, Zaccai G. Relative role of anions and cations in the stabilization of halophilic malate dehydrogenase. Biochemistry. 1999;38:9039–47.
Thomas LV, Wimpenny JW. Investigation of the effect of combined variations in temperature, pH, and NaCl concentration on nisin inhibition of Listeria monocytogenes and Staphylococcus aureus. Appl Environ Microb. 1996;62:2006–12.
Porter JL, Manning J, Sabatini S, Tavanti M, Turner NJ, Flitsch SL. Characterisation of CYP102A25 from B. marmarensis and CYP102A26 from P. halophilus: P450 homologues of BM3 with preference towards hydroxylation of medium chain fatty acids. Chembiochem. 2017;19:513–20.
Porter JL, Sabatini S, Manning J, Tavanti M, Galman JL, Turner NJ, Flitsch SL. Cloning, expression and characterisation of P450-Hal1 (CYP116B62) from Halomonas sp. NCIMB 172: a self-sufficient P450 with high expression and diverse substrate scope. Enzyme Microb Tech. 2018;113:1–8.
Dym O, Mevarech M, Sussman JL. Structural features that stabilize halophilic malate dehydrogenase from an archaebacterium. Science. 1995;267:1344–6.
Britton KL, Baker PJ, Fisher M, Ruzneinikov S, Grlmour DJ, Bonete MJ, Ferrer J, Pire C, Esclapez J, Rice DW. Analysis of protein solvent interactions in glucose dehydrogenase from the extreme halophile Haloferax mediterranei. Proc Natl Acad Sci USA. 2006;103:4846–51.
Pica A, Krauss IR, Castellano I, Cara FL, Graziano G, Sica F, Merlino A. Effect of NaCl on the conformational stability of the thermophilic γ-glutamyltranspeptidase from Geobacillus thermodenitrificans: implication for globular protein halotolerance. BBA Proteins Proteom. 2013;1834:149–57.
Ishibashi M, Yamashita S, Tokunaga M. Characterization of halophilic alkaline phosphatase from Halomonas sp. 593, a moderately halophilic bacterium. Biosci biotechnol biochem. 2005;69:1213–6.
Packer MS, Liu DR. Methods for the directed evolution of proteins. Nat Rev Genet. 2015;16:379–94.
Lane MD, Seelig B. Advances in the directed evolution of proteins. Curr Opin Chem Biol. 2014;22:129–36.
Arnold FH, Volkov AA. Directed evolution of biocatalysts. Curr Opin Chem Biol. 1999;3:54–9.
Ho SN, Hunt HD, Horton RM, Pullen JK, Pease LR. Site-directed mutagenesis by overlap extension using the polymerase chain-reaction. Gene. 1989;77:51–9.
Omura T, Sato R. The carbon monoxide-binding pigment of liver microsomes. I. Evidence for its hemoprotein nature. J Biol Chem. 1964;239:2370–8.
Guan W, Zhao H, Lu X, Wang C, Yang M, Bai F. Quantitative analysis of fatty-acid-based biofuels produced by wild-type and genetically engineered cyanobacteria by gas chromatography-mass spectrometry. J Chromatogra A. 2011;1218:8289–93.
SL, HX and YJ conceived and designed the study. YJ, ZL and CW performed the experiments and analyzed the results. YJ, HX and SL wrote the manuscript. YJZ discussed the results and revised the manuscript. All authors read and approved the final manuscript.
We thank Prof. David Nelson for his kind help for nomenclature of the three new P450 enzymes.
The authors declare that they have no competing interests.
Availability of supporting data
All data generated or analyzed during this study are included in this manuscript (and its additional files).
Consent for publication
Ethics approval and consent to participate
This work was supported by the National Science Foundation of China (NSFC 31770844), the DICP&QIBEBT program (Grant: DICP & QIBEBT UN201706) and the Dalian National Laboratory for Clean Energy (DNL), CAS, the Chinese Academy of Sciences (QYZDB-SSW-SMC042), and the Shandong Provincial Natural Science Foundation (ZR2017ZB0207).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.