Two New Unspecific Peroxygenases from Heterologous Expression of Fungal Genes in Escherichia coli

UPOs catalyze regio- and stereoselective oxygenations of both aromatic and aliphatic compounds. Similar reactions were previously described for cytochrome P450 monooxygenases, but UPOs have the noteworthy biotechnological advantage of being stable enzymes requiring only H2O2 to be activated. Both characteristics are related to the extracellular nature of UPOs as secreted proteins. In the present study, the limited repertoire of UPO enzymes available for organic synthesis and other applications is expanded with the description of two new ascomycete UPOs obtained by Escherichia coli expression of the corresponding genes as soluble and active enzymes. Moreover, directed mutagenesis in E. coli, together with enzyme molecular modeling, provided relevant structure-function information on aromatic substrate oxidation by these two new biocatalysts.

often higher, although the specific activities were significantly lower. It is worth mentioning that activities at the purification process were measured under pH and substrate concentration conditions that are not the optimal conditions determined later with the purified enzymes, as described below.
Characterization of the new UPOs. The molecular masses of the reduced and denatured UPOs of C. virescens (ϳ24 kDa) and D. caldariorum (ϳ27 kDa) estimated by SDS-PAGE together with standard proteins ( Fig. 2A and B insets, respectively) were slightly lower than those calculated from their amino acid sequences (29,616 and 29,609 Da, respectively). Moreover, the masses of the native (not unfolded/denatured) UPOs, estimated by Sepharose 12 SEC (Fig. 2) calibrated with selected standards (see Fig. S5 in the supplemental material), corresponded to 52 kDa and 55 kDa (Fig. 2, insets) for the C. virescens and D. caldariorum UPOs, respectively. These results revealed that both UPOs most probably occur as dimers.
The amino acid sequences of the peptides from trypsin hydrolysis of the recombinant C. virescens and D. caldariorum UPOs obtained by nano-scale liquid chromatographic-tandem mass spectrometry (nLC-MS/MS), which covered 80% and 70% of the total protein, respectively (Fig. S2), were identical to the initial UPO sequences used for nucleotide sequence synthesis and optimization. This agreement confirmed that the purified proteins corresponded to the synthesized genes used for E. coli transformation and that no errors were made during the gene synthesis and expression and enzyme purification processes.
Proper folding of and cofactor incorporation into the enzymes produced in E. coli were confirmed by their UV-visible (UV-vis) spectra (at the resting state), showing the strong Soret band at 420 to 430 nm and the two minor peaks at ϳ540 and ϳ580 nm ( Fig. 3A and C). From these spectra, Reinheitszahl values (A 420 /A 280 ratios) of 1.6 and 1.3 were obtained for the C. virescens and D. caldariorum UPOs, respectively. Moreover, the presence of a heme-thiolate cofactor (with a cysteine ligand of heme iron) was confirmed by the Soret band displacement to 440 to 450 nm in the spectra of enzyme-CO complexes ( Fig. 3B and D), as reported for peroxygenases and cytochrome P450s (5,17). The heme content in the purified enzyme samples (ϳ0.5 mol of heme per mole of protein) was measured with the pyridine hemochrome method and used to calculate the extinction coefficients applied to estimate the concentration of active enzyme in kinetic and other analyses.
Catalytic properties. Four usual peroxidase/peroxygenase substrates, i.e., the simple aromatics veratryl alcohol, benzyl alcohol, and naphthalene and the aromatic azo dye ABTS, were used to investigate the kinetic properties of the new C. virescens and D. caldariorum UPOs. For proper comparison, the kinetic curves against different substrate concentrations were estimated at optimal pH values determined previously (see Fig. S6 in the supplemental material). Veratryl alcohol oxidation by the C. virescens  S7D) exhibits sigmoidal kinetics (fitted to the Hill equation). From the above-described kinetic curves, the appropriate equations (shown in Materials and Methods) were used to calculate the turnover numbers (k cat ), Michaelis constants (K m ), and catalytic efficiencies (k cat /K m ) of the two enzymes oxidizing each of the four substrates, together with the inhibition constant (K i ) and Hill coefficient (n H ) values, when applicable ( Table 2). Inhibition of the D. caldariorum UPO by veratryl alcohol occurs at high substrate concentrations, with a K i Ͼ280-fold higher than the K m , while the K i values were only ϳ50 and ϳ30-fold higher than the K m values for benzyl alcohol oxidation by the D. caldariorum UPO and ABTS oxidation by the C. virescens UPO, respectively. Finally, the positive cooperative degree (n H ) in the sigmoidal kinetics observed for the C. virescens UPO revealed that binding of one benzyl alcohol molecule improves the enzyme interaction with additional substrate molecules.
Regarding catalytic efficiency, apparent affinity, and maximal turnover values (Table  2), both similarities and differences between the two UPOs were observed. First, ABTS appears to be the best substrate for both of them, as shown by catalytic efficiency, which is notably high for the C. virescens UPO due to its high turnover number (k cat of ϳ160 s Ϫ1 ). Concerning the simple aromatic substrates assayed (veratryl and benzyl alcohols and naphthalene), it is interesting that all of them were better oxidized by the D. caldariorum UPO. This is due to the extremely low apparent affinity of the C. virescens UPO for both veratryl (with a K m of Ͼ2,900 M) and benzyl (with a K m of Ͼ7000 M) alcohols and its very low turnover number with naphthalene (k cat of ϳ1 s Ϫ1 ).
Finally, oleic acid was used to explore the fatty acid oxyfunctionalization ability of the two UPOs, analyzed by gas chromatography-mass spectrometry (GC-MS). As shown in Fig. 4A, 9,10-epoxyoctadecanoic acid is the main product of the oleic acid reaction with the C. virescens UPO, followed by 11-hydroxyoctadec-9-enoic acid [labeled (-7)-OH] and traces of hydroxy/keto derivatives of oleic acid. In contrast, no epoxides were found in the reaction with D. caldariorum UPO, which generated only the four 14-  hydroxy to 17-hydroxy (-4 to -1) derivatives and the two 14-keto and 15-keto (-4 and -3) derivatives of oleic acid, together with a small amount of a hydroxy/keto derivative, as shown in Fig. 4B. For comparative purposes, the results from oleic acid treatment with M. rotula UPO (12), also obtained by heterologous expression in E. coli, are included in Fig. 4C. Analysis of site-directed variants at the heme pocket. Enzyme molecular models were built, using available UPO crystal structures as templates, to explain the very different kinetic constants of the two new UPOs on simple aromatic substrates (and ABTS). Residues at the upper side of the heme pocket (opposite the lower side with a proximal cysteine acting as heme iron ligand) were compared in homology molecular models of the two new UPOs ( Fig. 5A and B) and crystal structures of the A. aegerita and M. rotula UPOs ( Fig. 5C and D). In addition to the proximal cysteine (Cys19/Cys15), three residues (His90/His86, Glu162/Glu156, and Tyr166/Tyr160) are conserved in the C. virescens/D. caldariorum UPOs. Both enzymes also share the presence of a phenylala- nine residue at this side of the heme cavity, being Phe88 in the C. virescens UPO and Phe152 in the D. caldariorum UPO (with Leu84 and Thr158 as their homologous residues in the second enzyme).
To investigate the functional implications of the above-described structural differences in active-site architecture, the F88L and T158F variants of the C. virescens UPO were obtained by site-directed mutagenesis (the double variant could not be expressed as a soluble active enzyme in E. coli) and purified to electrophoretic homogeneity (see Fig. S8 in the supplemental material) by the same protocol used for the wild-type enzymes, showing similar UV-vis spectra (see Fig. S9 in the supplemental material). Their optimal pH values (see Fig. S10 in the supplemental material) and kinetic curves (see Fig. S11 in the supplemental material) for the same aromatic substrates used with the wild-type enzymes were then obtained, and the corresponding kinetic constants were calculated.
Replacing Phe88 in the C. virescens UPO by a leucine residue, as found in the D. caldariorum UPO, strongly increased its catalytic efficiency in oxidizing veratryl alcohol (21-fold), benzyl alcohol (11-fold), and naphthalene (5-fold) ( Table 3). These differences  Table 3, compared with the same region in the crystal structures of the A. aegerita (C) and M. rotula (D) UPOs (PDB entries 2YP1 and 5FUJ, respectively).  Table 2 for the wild-type UPOs (at the optimal pH conditions shown in Fig. S10 in the supplemental material).
are on the same order of those found when the two wild-type enzymes were compared, with the D. caldariorum UPO showing 10-to 30-fold-higher catalytic efficiency in oxidizing these aromatic substrates than the C. virescens UPO (Table 2). In contrast, the F88L variation slightly decreased the efficiency of the C. virescens UPO in oxidizing ABTS, in agreement with the preference for this dye shown by this enzyme compared with the D. caldariorum UPO, which was more evident for the T158F variant described below.
The improvement in catalytic efficiency on simple aromatic substrates caused by the F88L mutation is due to a significant decrease of the K m , indicating that in the absence of this phenylalanine residue, a smaller amount of substrate is needed to attain maximal transformation. The additional turnover increase observed for veratryl alcohol oxidation by the F88L variant would be due to better (closer) positioning of the substrate in the Phe88-less variant.
Concerning the T158F mutation, the corresponding variant (Table 3) improved the oxidation only of benzyl alcohol but decreased the oxidation of veratryl alcohol and naphthalene compared with the wild-type enzyme ( Table 2). Most probably, the access of these two bulkier aromatic compounds (compared with benzyl alcohol) is hampered by the introduction of a phenylalanine residue at the active site of the C. virescens UPO. However, the effect of introducing a bulky side chain in the T158F mutant was especially evidenced by the 20-fold decrease of catalytic efficiency in oxidizing the large molecule ABTS, in agreement with the lower preference of the D. caldariorum UPO for ABTS than for simple aromatics.

DISCUSSION
UPOs are a promising enzyme family to be explored. The best known peroxideactivated enzymes (peroxidases), such as classic horseradish peroxidase (HRP), prokaryotic cytochrome c peroxidase (CcP), and fungal ligninolytic peroxidases (18,19), belong to the catalase-peroxidase superfamily (20), with the latter two being the first (21) and second/third (22)(23)(24) peroxidases with a known molecular structure. In recent years, new peroxide-activated enzymes have been described, including the so-called dyedecolorizing peroxidases (DyPs) in the CDE superfamily and the fungal UPO-type enzymes (including CPO) in the HTP superfamily (25)(26)(27)(28)(29). UPOs have been characterized from a variety of fungi, bacteria, and archaea (29,30), but amazingly, only the structure of the A. aegerita UPO has been published to date (31), together with the classic CPO (32). Thus, the structural-functional information and biotechnological potential of UPOs are strongly limited by the low number of isolated proteins, despite the facts that the first UPO was described 15 years ago (5) and a high number of putative UPO-encoding genes have been identified in sequenced fungal genomes and other databases (3).
The low number of wild (nonrecombinant) UPOs available is due to difficulties in (i) growing higher fungi in pure culture due to complex growth requirements, (ii) detecting in pure culture the UPOs identified in genomes/databases, and (iii) obtaining satisfactory amounts of enzyme. Moreover, the heterologous expression of UPO genes, which should overcome these difficulties, is far from being a straightforward process. Among the host systems used for expressing UPO-encoding genes, production using A. oryzae is a proprietary technology of the company Novozymes (Bagsvaerd, Denmark), and therefore, the strains, vectors, and growth conditions are not publicly available. Yeast expression, as far as reported, requires previous directed evolution introducing mutations that result in structural differences with respect to the wild-type enzymes (33). Overexpression in E. coli as inclusion bodies followed by in vitro activation is routinely used for fungal (class II) peroxidases (34)(35)(36)(37)(38)(39), but a protocol for in vitro activation of UPOs from inclusion bodies has not been developed yet.
Enlarging UPO diversity. Fortunately, successful expression of the previously described UPO of M. rotula (7) as a soluble and active recombinant enzyme (rMroUPO) in E. coli has been recently reported (12), covered by a patent application (16). This method was applied here to further explore UPO diversity. In this way, we were able to describe the eighth and ninth UPOs, from C. virescens and D. caldariorum, after those previously described from A. aegerita, C. radians, C. cinerea, C. globosum, H. insolens, M. rotula, and M. wettsteinii (5)(6)(7)(8)(9)11). The original nucleotide sequences of the C. virescens and D. caldariorum UPOs are not available. Only the amino acid sequences were provided by Lund et al. (15), with a putative peroxygenase of Podospora comata (accession no. VBB75842.1) and a hypothetical protein from the sequenced genome (40) of a Daldinia species (accession no. OTB17553.1) being the most related sequences in the NCBI database (76 to 77% amino acid identity), respectively. However, optimized nucleotide sequences for E. coli expression as active enzymes could be obtained, and they are now the second and third active recombinant UPOs heterologously expressed in E. coli. Interestingly, their production levels (7 and 3 mg of highly enriched enzyme per liter of bacterial culture, respectively) are higher than that obtained for the M. rotula UPO using the same expression system. However, the heterologous production of UPOs is not a straightforward process, as mentioned above, and only two of the eight UPO sequences assayed could be successfully transformed into soluble and active enzymes in E. coli.
The use of a lactose-containing self-induction medium (41) is an important aspect of E. coli expression of UPOs that, together with the low growth temperature, promotes proper protein folding as an active enzyme. Moreover, since UPOs are fungal proteins, the design and optimization of nucleotide sequences for E. coli expression is a key aspect for their heterologous production. With this aim, the Optimizer software (42) was used, providing both "fully optimized" sequences (i.e., those always using the most frequent codon for each amino acid in the host strain) and "randomly optimized" sequences (alternating several possible codons per amino acid). For UPO expression in E. coli, the first approach was not successful, most probably because of inefficient folding. Random optimization was then screened to identify sequences (see Fig. S1 in the supplemental material), enabling us to obtain soluble and active enzymes. Additional improvements in nucleotide sequence optimization for this E. coli system could further contribute to efficient expression of these (and other) UPO genes, providing an adequate speed balance between protein synthesis and folding (43).
The two new UPOs are ascomycete enzymes to be classified in the group of short UPOs (1), due to their shorter polypeptide chains of 240 to 270 residues (and molecular masses of around 25 kDa). The UPOs from C. virescens and D. caldariorum appear to be dimeric enzymes, as shown by the 2-fold-higher molecular mass estimated by SEC than by SDS-PAGE under reductive conditions. Cys241 of D. caldariorum UPO and one of the C. virescens UPO cysteine residues (Cys178, Cys230, or Cys235) could be involved in dimerization by formation of intermolecular disulfide bridges, as reported for the M. rotula UPO (44) based on its crystal structure (PDB no. 5FUJ and 5FUK). Another dimeric UPO has been reported from M. wettsteinii, based on combined SDS-PAGE and SEC results (8).
Catalytic properties of the new UPOs. In general terms, the optimal pH for catalysis by the two new enzymes (pH 3 to 6) is more acidic than reported for the previously characterized UPOs; e.g., pH optima near neutrality have been reported for veratryl alcohol oxidation by the A. aegerita and C. globosum UPOs, as well as for naphthalene oxidation by the latter enzyme and the M. rotula UPO (5,7,9). For the two new enzymes, ABTS is also the most efficiently transformed substrate of the A. aegerita and C. globosum UPOs and the second-better substrate of the M. rotula UPO after benzyl alcohol (5,7,9). In a similar way, naphthalene is less efficiently oxidized by the two new enzymes, in agreement with results reported for other UPOs (7,9). Among the two other aromatic compounds, benzyl alcohol is much better oxidized than veratryl alcohol, as also reported for the other characterized UPOs (5,7,9), a fact that can be related to its less bulky nature (due to the absence of methoxy substituents), facilitating its access to the buried heme cofactor.
In addition to the above-mentioned aromatics, UPOs are also active on aliphatic substrates. This includes unsaturated fatty acids, as recently shown for the M. rotula UPO (12,45). In a preliminary analysis, the oxyfunctionalization pattern for the latter enzyme was compared with those obtained for the two new UPOs using oleic acid as the substrate. Interestingly, the three enzymes showed different oxyfunctionalization patterns. While a mixture of epoxide and hydroxy-epoxides (with hydroxy groups at -7, -1, and other positions) had been reported for the wild (45) and recombinant (12) M. rotula UPO, nearly full conversion with the new UPOs yields as a main product either epoxide by the C. virescens UPO or a variety of subterminal hydroxy/keto derivatives by the D. caldariorum UPO. A complete study of unsaturated fatty acid epoxidation by the C. virescens UPO has been recently reported (54). While epoxides are investigated as reactive molecules for different applications, hydroxy fatty acids are of interest in the manufacture of ester homopolymers of fatty acids.
Structural-functional aspects of UPO activity on aromatic substrates. Both the reductive (on H 2 O 2 ) and oxidative (on aromatic and aliphatic substrates) UPO halfreactions (Fig. 1, right and left, respectively) take place at the upper side of the heme pocket (1). The main structural difference between the upper heme pockets of the UPOs of C. virescens and D. caldariorum concerns the location of a phenylalanine residue, being Phe88 and Phe158, respectively. A key role of phenylalanine residues at the heme pocket and access channel on UPO activity on aromatic and other substrates has been suggested in the first structural-functional studies on these enzymes, discussed below.
The first UPO, isolated from A. aegerita, has up to 27 phenylalanine residues in its mature protein structure (Fig. 6C). The involvement of several heme environment and channel residues (namely, Phe69, Phe76, Phe121, Phe191, and Phe199) in oxidation of polycyclic and other aromatic substrates by this UPO was suggested when its crystal structure was reported for the first time (31). Inversely, the preference of the M. rotula UPO for aliphatic compounds, compared with the A. aegerita UPO preferring aromatics (3,(45)(46)(47), could be related to the much lower number of phenylalanine residues at its heme pocket and access channel (only 9 phenylalanines in the whole protein structure) (Fig. 6D). An intermediate number of phenylalanines (17 and 16, respectively) are present in the new C. virescens and D. caldariorum UPOs (Fig. 6A and B, respectively). Interestingly, a docking study of the A. aegerita UPO evolved for heterologous expression in yeast (the so-called PaDa-I variant) identified three of the above-mentioned phenylalanine residues (Phe121, Phe199, and, to a lesser extent, Phe69) as being involved in interaction with all the (aromatic) substrates assayed (33).
As mentioned above, the two new UPOs exhibited differences in the location of an active-site phenylalanine ( Fig. 5A and B) and strongly different kinetic constants on aromatics that were investigated by directed mutagenesis of the C. virescens enzyme. Remarkably, replacing Phe88 by a leucine residue, as found in the D. caldariorum UPO, strongly increased the catalytic efficiency of the C. virescens UPO in oxidizing veratryl alcohol and naphthalene (22-fold and 5-fold, respectively), with parallel (8-fold and 2-fold) increases in k cat , although Phe88 occupies a position equivalent to the abovementioned Phe121 in the A. aegerita UPO (Fig. 5C). The opposite tendency was observed when Thr158 (similar to Thr192 in the A. aegerita UPO [Fig. 5C]) was replaced in the T158F mutation, which decreases (Ͼ20-fold) the UPO catalytic efficiency for the aromatic azo dye ABTS, with a simultaneous decrease of k cat (14-fold) and increase of K m (1.5-fold), with the catalytic activity of the C. virescens UPO again approaching that of the D. caldariorum UPO.
In conclusion, the structural-functional information on the two new UPOs supports the importance of active-site phenylalanines on UPO reactions with aromatic compounds (31,33). However, it clearly associates the higher activity on simple aromatics in these UPOs with the absence of a phenylalanine residue at the position equivalent to Phe121 in the A. aegerita UPO, while the opposite role has been proposed for this residue in the latter enzyme. In a similar way, the absence of a second phenylalanine in position 158 of the C. virescens UPO (which is also absent in the A. aegerita UPO) confers to the former enzyme its remarkable activity on ABTS compared with that of the D. caldariorum UPO.

MATERIALS AND METHODS
Enzyme production. The protein sequences of eight putative UPOs from the ascomycetes C. virescens, D. caldariorum, and Xylaria polymorpha, the basidiomycetes Laccaria bicolor and Rhodotorula graminis, the mucoromycete Rhizopus delemar, the glomeromycete Rhizophagus irregularis, and the chytrid Spizellomyces punctatus were obtained from NCBI database and reference 15. DNA sequences coding for these proteins were optimized for E. coli expression using the different options provided by the Optimizer software (42), synthesized by ATG:biosynthetics GmbH (Germany), and cloned in the pET23a and His-tag-containing pET28 plasmids under the control of the T7lac promoter.
The plasmids were transformed into competent E. coli C41 cells, which were grown for 4 to 5 days in lactose-containing autoinduction medium ZYM-5052 (41) supplemented with hemin at low temperature and agitation (16°C and 180 rpm) to obtain UPOs as soluble and active proteins (12,16). Cells were harvested (10 min at 8,000 rpm) and suspended in 50 to 70 ml of 10 mM Tris (pH 7.4), 5 mM dithiothreitol, and 20 mM EDTA for cultures with pET23-cloned genes or in 10 mM Tris (pH 8), 0.3 M NaCl, and 20 mM imidazole for cultures with the pET28-cloned genes. Lysozyme (2 mg · ml Ϫ1 ), DNase, and sonication were used to lyse the cells, and the soluble fractions containing the recombinant UPOs were recovered after debris removal (45 min at 15,000 rpm followed by 1 h at 36,000 rpm).
As a step before protein purification, an assay to detect the presence/absence of UPO in the E. coli extracts was performed based on the capacity of reduced (ferrous) heme-thiolate proteins to form a complex with CO, which shifts the main peak from 420 nm to 440 to 450 nm, as shown in difference spectra (48). For this purpose, the soluble fractions (200 l) were mixed with 200 l of phosphate (pH 8) and adjusted to 1 ml with Milli-Q water, and their spectra (300 to 700 nm) were recorded with a Cary 60 spectrophotometer before (Fe 3ϩ resting-state enzyme) and after (Fe 2ϩ reduced enzyme) the addition of Na 2 S 2 O 4 , followed by flushing with CO (CO-enzyme complex).
Enzyme purification. The C. virescens and D. caldariorum recombinant proteins showed the 440-to 450-nm maximum in the above-described assay and were therefore subjected to further purification using two different approaches, depending on whether the gene was cloned in the pET23a or in the pET28 plasmid.
The C. virescens recombinant enzyme could be successfully purified from the E. coli lysates with the gene cloned in the plasmid without a His tag (pET23) by using two chromatographic steps in an Äkta (GE Healthcare) fast protein liquid chromatography (FPLC) system. Since the C. virescens UPO has a theoretical isoelectric point of 7.7, cation-exchange chromatography was first applied using a HiTrap SPFF column (GE Healthcare, USA) in 10 mM Tris (pH 7.4) to retain positively charged proteins. The retained proteins were eluted with a gradient of the same buffer supplemented with 1 M NaCl and concentrated with an Amicon 3K instrument (Sigma-Aldrich, USA). The second purification step was SEC with a Superdex 75 column (10/300 GL; GE Healthcare, USA) in 10 mM Tris (pH 7.4) with 0.15 M NaCl. Protein elution was followed at 420 nm (heme proteins) and 280 nm (all proteins).
Purification of the D. caldariorum enzyme was facilitated by pET28 cloning, after unsuccessful ion-exchange purification of the protein obtained with the untagged pET23a (used for the C. virescens enzyme). Since pET28 has a His tag, the purification was carried out by immobilized metal ion affinity chromatography with a HiTrap IMAC FF 5-ml column, using the same FPLC system. First, the column was loaded with 0.1 M NiSO 4 to provide the metal ions where the tail histidines will be bound after sample injection in 10 mM Tris (pH 7.4), 0.3 M NaCl, and 20 mM imidazole. The retained enzyme was eluted with 0.5 M imidazole buffer, using a 0 to 50% gradient at a flow rate of 0.5 ml · min Ϫ1 . The UPO-containing fractions were concentrated in Amicon ultracentrifuge filters and applied to an SEC column (Superdex 75 10/300 GL) using an isocratic flow of 10 mM Tris (pH 7.4) and 0.15 M NaCl at 0.4 ml · min Ϫ1 . As in the previous case, protein elution was followed at 420 and 280 nm. The eventual effect of the poly-His tail on the enzyme activity was evaluated after its removal with thrombin (49).
Purification yield was followed by quantifying proteins at 280 nm (in an Thermo Scientific NanoDrop 2000) and activity by measuring 1 mM ABTS oxidation in 0.1 M tartrate (pH 3) using 1 mM H 2 O 2 (5 mM veratryl alcohol oxidation in 0.1 M Tris [pH 7.4] was also monitored) along the different purification steps. Enzyme purification was followed by 12% PAGE in the presence of 0.1% SDS and ϳ1% mercaptoethanol (present at a 5% concentration in the loading buffer), reducing disulfide bridges, and unfolding proteins (50). The pyridine hemochrome assay was used to assess the concentration of heme in the purified enzymes by adding pyridine to the enzyme reduced with sodium dithionite (51). Molar extinction coefficients were then calculated according to Beer's law using the above concentration of heme and the enzyme absorbance at 420 nm. The resulting 420 values were 114.2 mM Ϫ1 · cm Ϫ1 and 98.3 mM Ϫ1 · cm Ϫ1 for the recombinant C. virescens and D. caldariorum UPOs, respectively.
In addition to the above-described SDS-PAGE analyses, the native molecular masses of nondenatured purified UPOs were estimated by SEC in a Sepharose 12 10/300 column (GE Healthcare, USA) using 10 mM Tris (pH 7.4) with 0.15 M NaCl and ferritin (460 kDa), aldolase (170 kDa), bovine serum albumin (73.5 kDa), ovalbumin (44.0 kDa), carbonic anhydrase (29.3 kDa), and aprotinin (6.5 kDa) as molecular mass standards. Proper folding and binding of the cofactor were evaluated by inspecting the UV-vis spectrum of the resting state of the enzymes in 10 mM phosphate (pH 7.4) using a Cary 60 spectrophotometer. Formation of a complex between the chemically reduced enzyme (ferrous form) and CO, characteristic of active heme-thiolate enzymes as mentioned above, was assessed in 0.2 M phosphate (pH 8) after addition of Na 2 S 2 O 4 and CO flushing.
Peptide sequencing. Identification of the purified recombinant UPOs of C. virescens and D. caldariorum was confirmed by nLC-MS/MS analysis of the trypsin-digested peptides (using the Easy-nLC-LTQ Orbitrap Velos equipment) and comparison with the peptides from theoretical trypsin digestion of the C. virescens and D. caldariorum UPOs. N-terminal sequencing was carried out by automatic Edman degradation (52) (using the Procise 494 protein sequencer from Applied Biosystems). These assays were performed at the CIB (CSIC) Proteomics and Genomics Facility (member of the ProteoRed-ISCIII network) and the Protein Chemistry Facility, respectively. Molecular modeling and directed mutagenesis. The molecular structures of the C. virescens and D. caldariorum UPOs were modeled at the Swiss-Model server (https://swissmodel.expasy.org) (53) with the A. aegerita (PDB no. 2YP1) and M. rotula (PDB no. 5FUJ) UPO structures as templates. The F88L and T158F mutations were introduced in the C. virescens UPO gene using the Expand Long Template PCR kit from Roche (Basel, Switzerland). PCRs were run using ad hoc oligonucleotides (along with their reverse complementary counterparts) (mutated triplets in bold and substituted nucleotides underlined): (i) F88L mutation, 5=-AAC CGC CAT AAC CTG TTG GAA CAT GAT GCA TCT C-3=; and (ii) T158F mutation, 5=-ACT TAC ACC GTT CAG CAG CGT ATC TTT AGT TAC GGT GAA ACG-3=.
The PCRs (50-l volume) were carried out in an Eppendorf (Hamburg, Germany) Mastercycler pro-S using 30 ng of template DNA, 0.5 mM each deoxynucleoside triphosphate (dNTP), 125 ng of direct and reverse primers, 5 units of Expand Long Template PCR system polymerase mix (Roche), and the manufacturer's buffer. Reaction conditions were as follows: (i) initial denaturation step of 1 min at 95°C; (ii) 22 cycles of 30 s at 95°C, 30 s at 60°C, and 7 min at 68°C; and (iii) final elongation step of 7 min at 68°C. The mutated variants were produced in E. coli and purified as described above for the wild-type enzyme.
Kinetic curves were obtained at the optimal pH for each enzyme and substrate by varying the concentrations of veratryl alcohol (from 0.07 M to 4.9 mM, in 0.1 M acetate), benzyl alcohol (from 20 M to 0.1 M, in 0.1 M Tris), naphthalene (from 15.1 M to 0.5 mM, in 0.1 M tartrate with 5% ethanol), and ABTS (from 0.07 M to 4.9 mM, in 0.1 M acetate) in the spectrophotometer cuvette (maintained at 25°C) and adding a final concentration of 1 mM H 2 O 2 and 20 to 700 nM enzyme. The reactions were triggered by the H 2 O 2 addition, the time was adjusted to ensure that the reaction rate was in linear phase, and values were calculated as the change in absorbance over time. All reactions were carried out in triplicate.
The plotting and analysis of the above-described curves were done with SigmaPlot 11.0 by nonlinear least-squares fitting the k obs values to equation 1 (Michaelis-Menten model), except for (i) veratryl and benzyl alcohol oxidation by D. caldariorum UPO and ABTS oxidation by the C. virescens UPO and the F88L variant of C. virescens UPO, which were adjusted to equation 2 (inhibition model) and (ii) naphthalene oxidation by both UPOs and benzyl alcohol oxidation by the C. virescens UPO, which were adjusted to equation 3 (Hill model): where k cat is the turnover number (corresponding to the number of substrate molecules that are oxidized by one cofactor molecule per second), S the substrate concentration, K m the Michaelis constant, K i the inhibition constant, y 0 the y axis intercept, n H the Hill constant, and K 0.5 the substrate concentration for half saturation. Fatty acid oxyfunctionalization. Oleic acid was used to investigate the oxyfunctionalization abilities of the two new UPOs. For this purpose, 30-min reactions were performed with 0.1 mM substrate, 0.8 to 1.4 M enzyme, and 1.25 to 2.50 mM H 2 O 2 in 50 mM phosphate (pH 7) containing 20% acetone. Chromatographic analyses were performed with QP2010 Ultra GC-MS equipment, using a fused-silica DB-5HT capillary column from J&W Scientific (USA). The oven was heated from 120°C (1 min) to 300°C (15 min) at 5°C · min Ϫ1 . The injection was performed at 300°C, and the transfer line was kept at 300°C. Compounds were identified by mass fragmentography and by comparison of their mass spectra with those of authentic standards. Quantifications were obtained from total-ion peak areas, using external standard curves and molar response factors of the same or similar compounds.
Data availability.
All data are open-access available in the main article and its supplemental material.

SUPPLEMENTAL MATERIAL
Supplemental material is available online only. SUPPLEMENTAL FILE 1, PDF file, 0.6 MB.