Next Article in Journal
Tailoring the Electrical Energy Storage Capability of Dielectric Polymer Nanocomposites via Engineering of the Host–Guest Interface by Phosphonic Acids
Next Article in Special Issue
Experimental and Theoretical Study on the Homodimerization Mechanism of 3-Acetylcoumarin
Previous Article in Journal
The Design of 3D-Printed Polylactic Acid–Bioglass Composite Scaffold: A Potential Implant Material for Bone Tissue Engineering
Previous Article in Special Issue
Thiocoumarins: From the Synthesis to the Biological Applications
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 21
10.3390/molecules27217230
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Challenges in the Heterologous Production of Furanocoumarins in Escherichia coli

by
Joana L. Rodrigues
1,2,*,
Daniela Gomes
1,2 and
Lígia R. Rodrigues
1,2
1
CEB—Centre of Biological Engineering, University of Minho, 4710-057 Braga, Portugal
2
LABBELS—Associate Laboratory, 4710-057 Braga, Portugal
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(21), 7230; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27217230
Submission received: 6 September 2022 / Revised: 13 October 2022 / Accepted: 20 October 2022 / Published: 25 October 2022
(This article belongs to the Special Issue Coumarin and Its Derivatives II)
Browse Figures
Versions Notes

Abstract

:
Coumarins and furanocoumarins are plant secondary metabolites with known biological activities. As they are present in low amounts in plants, their heterologous production emerged as a more sustainable and efficient approach to plant extraction. Although coumarins biosynthesis has been positively established, furanocoumarin biosynthesis has been far more challenging. This study aims to evaluate if Escherichia coli could be a suitable host for furanocoumarin biosynthesis. The biosynthetic pathway for coumarins biosynthesis in E. coli was effectively constructed, leading to the production of umbelliferone, esculetin and scopoletin (128.7, 17.6, and 15.7 µM, respectively, from tyrosine). However, it was not possible to complete the pathway with the enzymes that ultimately lead to furanocoumarins production. Prenyltransferase, psoralen synthase, and marmesin synthase did not show any activity when expressed in E. coli. Several strategies were tested to improve the enzymes solubility and activity with no success, including removing potential N-terminal transit peptides and expression of cytochrome P450 reductases, chaperones and/or enzymes to increase dimethylallylpyrophosphate availability. Considering the results herein obtained, E. coli does not seem to be an appropriate host to express these enzymes. However, new alternative microbial enzymes may be a suitable option for reconstituting the furanocoumarins pathway in E. coli. Nevertheless, until further microbial enzymes are identified, Saccharomyces cerevisiae may be considered a preferred host as it has already been proven to successfully express some of these plant enzymes.

1. Introduction

Coumarins are phenolic compounds produced in plants as secondary metabolites. Most known coumarins include coumarin, umbelliferone, scopoletin and esculetin and have several interesting biological activities, including anti-inflammatory, anticancer, neuroprotective, antifungal, antibacterial, antidiabetic, antiepileptic, cardiovascular protective, antiulcerogenic, and antidiarrheal, among others [1,2,3,4,5,6,7,8,9]. For example, several coumarins have been studied in clinical trials and demonstrated potential in the treatment of leukemia, renal, and prostate cancers, among others [10,11,12]. In addition, the synthetic coumarin derivative warfarin is frequently prescribed as an anticoagulant for thrombosis treatment [13]. Umbelliferone derivatives such as furanocoumarins have also reported several therapeutic properties [14]. Furanocoumarins (frequently simply referred as psoralens), such as xanthotoxin and bergapten, have been used in clinics for several years to treat severe cases of psoriasis and vitiligo and as a first line of treatment of mycosis fungoides [15]. Currently, there are several pharmaceuticals on the market (e.g., Uvadex®, Oxsoralen®, 8-MOP®, and Melanocyl®) with these active principles approved by the EMA/FDA that are used in PUVA therapy consisting of the combination of psoralens with UVA radiation. Furanocoumarin’s potential is also being evaluated for treating cancer, graft-versus-host disease and in the prevention and treatment of solid organ transplantation rejection, among others [16,17,18]. In addition to their application as pharmaceuticals, several coumarin derivatives including furanocoumarins could be used as agrochemicals due to their phytotoxic, antibacterial, antifungal, insecticide, and herbicide properties [19,20,21]. Coumarins and their derivatives, as other plant secondary metabolites, are produced in very low amounts in plants. Therefore, the extraction process, besides being environmentally unfriendly, is also costly and ineffective [14]. Consequently, more efficient and sustainable production processes need to be developed. Heterologous production has been suggested as a potentially viable solution as it is not limited by plant availability, seasonality, or environmental factors and is a greener approach [22,23]. Additionally, metabolic engineering and synthetic biology tools allow the construction of more and more efficient strains capable of producing the heterologous value-added compounds with higher titers and yields [24,25]. This type of approach, although with its own challenges [14,26,27], has been successfully used to produce several compounds such as curcuminoids, flavonoids, and sesquiterpenoids produced in plants in low amounts [22,28,29]. Coumarins heterologous production has already been explored in Escherichia coli [30,31,32] and Saccharomyces cerevisiae [33]. The heterologous production of these compounds in microbes starts by the introduction of the phenylpropanoids pathway where tyrosine ammonia lyase (TAL) converts tyrosine to p-coumaric acid (Figure 1). Then, 4-coumarate-CoA ligase (4CL) converts p-coumaric acid to p-coumaroyl-CoA. This enzyme is also able to convert other hydroxycinnamic acids, such as ferulic and caffeic acid to the respective CoA esters. Finally, a 2-oxoglutarate-dependent dioxygenase, which can be a feruloyl-CoA 6′-hydroxylase (F6′H) or a p-coumaroyl-CoA 2′-hydroxylase (C2′H), depending on the substrate specificity, converts p-coumaric acid, caffeic acid and ferulic acid to umbelliferone, esculetin and scopoletin, respectively. So far, up to 356.59 mg/L of umbelliferone were obtained in a E. coli tyrosine overproducing strain [31]. The other coumarins have also been produced but in lower amounts compared to umbelliferone [30,32].
The heterologous production of coumarin derivatives, in particular of umbelliferone derivatives has been much less explored. As far as we know, furanocoumarins have never been produced in microbial hosts. Some of the steps have been individually tested but the pathway has never been fully assembled [14]. So far, there is only one preliminary study that shows that umbelliferone can be converted in E. coli to demethylsuberosin by a plant prenyltransferase (PT) and then to marmesin by a microbial marmesin synthase (MS) [34]. Previous reports demonstrated that PTs could not be expressed in E. coli [35], being usually evaluated in the plant Nicotiana benthamiana [35,36,37]. Additionally, MS and psoralen synthase (PS) characterization studies [38,39,40] have also suggested that E. coli is not a good host to express these enzymes. Therefore, these steps need to be investigated further to elucidate if the heterologous production of these compounds in E. coli can be considered a viable alternative.
Therefore, the main goal of this work was to evaluate the possibility of producing furanocoumarins in E. coli through the in-depth study of several steps of the pathway. In this work, the heterologous production of coumarins was established and used as a basis for the study of the next steps in the furanocoumarin’s pathway (PT, MS, and PS steps). The heterologous production of umbelliferone, esculetin and scopoletin was well succeeded, leading to 128.7, 17.6, and 15.7 µM, respectively, from tyrosine. However, it was not possible to complete the biosynthetic pathway with a PT enzyme to produce demethylsuberosin from tyrosine. The PT, MS, and PS enzymes from the furanocoumarin’s pathway tested did not demonstrate catalytic activity in vivo or in vitro. Strategies, such as the removal of potential N-terminal transit peptides and expression of a cytochrome P450 reductase, chaperones, and/or enzymes to increase the availability of dimethylallylpyrophosphate (DMAPP) have been attempted to improve the enzyme’s solubility and activity, although without success. Therefore, it was concluded that E. coli might not be a suitable host for the expression of this pathway using plant enzymes. Future research may focus in using microbial enzymes or other hosts such as S. cerevisiae possibly adopting co-culture/sequential culture approaches.

2. Results and Discussion

2.1. Heterologous Production of Umbelliferone using C2′H

To convert p-coumaroyl-CoA to umbelliferone, we selected the enzyme IbF6′H2-2-1 from Ipomoea batatas [41]. This isoenzyme previously demonstrated to have similar affinity to p-coumaroyl-CoA and feruloyl-CoA as substrates [41] when expressed in E. coli BL21 as it was verified for other isoenzymes. Since the F6′H nomenclature is normally used on enzymes with a higher F6′H activity than C2′H, we herein refer to this enzyme simply as C2′H.
C2′H enzymes are produced in the cytoplasm of the plant whereby they do not have predictable signal sequences in the N- or C- terminals [42] as can be easily confirmed using different bioinformatic tools (LocTree3, iPSORT, SignalP, TargetP). However, it has been reported that C2′H/F6′H have low solubility when expressed in E. coli. Therefore, we tested the production of umbelliferone with C2′H not only in the commonly used high-copy plasmid pRSFDuet-1, but also in pET28GST-LIC that contains a glutathione S-transferase (GST)-tag. This tag has been previously used to increase F6′H stability and solubility [32]. Herein, C2′H cloned in different plasmids was combined with 4CL1 from Arabidopsis thaliana (pAC-4CL1) to produce umbelliferone from supplemented p-coumaric acid (0.7 mM). It was not possible to observe any production when C2′H was expressed in pRSFDuet-1 without any fusion tag, and the production was very residual when His or S-tags were used (Figure 2a). When C2′H was expressed with GST-tag, the umbelliferone production was significant demonstrating that this tag had a positive effect in C2′H expression although no significant differences were observed in the SDS-PAGE (sodium dodecyl sulfate polyacrylamide gel electrophoresis) gel (Figure S1). C2′H expressed in both plasmids was mostly in insoluble fraction.
The highest production herein obtained, 254.9 µM (41.33 mg/L), was higher than the one obtained by Lin et al. [30] (4.3 mg/L) that used a 4CL2 from A. thaliana and a F6′H from I. batatas (IbF6′H2-1-1) that also has a high specificity for p-coumaroyl-CoA [41]. However, it was considerably lower compared to the umbelliferone production (82.9 mg/L) obtained from p-coumaric acid in the Yang et al. [32] study. In this study, Oryza sativa 4CL and IbF6′H2-1-8 from I. batatas were used. IbF6′H2-1-8, according to in vitro studies, presents a higher catalytic efficiency and a higher preference for p-coumaroyl-CoA than feruloyl-CoA [41].
Previous studies reported the coumarins production in different culture media. The protein production phase is always performed in LB. However, in some cases in the coumarins production phase this LB is replaced by fresh LB, M9 (2% glucose) or YM9 (M9 2% glucose + 0.2% yeast extract) [30,31,32]. Therefore, we decided to also evaluate if the production phase or the complete experiment in M9 could be more favorable. Our previous studies in hydroxycinnamic acids and curcuminoids production demonstrated that M9 [43,44,45,46,47] could be an interesting option to produce these compounds. Preliminary studies (Figure 2a) were performed using LB in the protein production phase and then, fresh LB in the coumarin production phase, as it was demonstrated by Yang et al. [32] that fresh LB obtained more scopoletin than M9 or YM9. Our study also confirmed that fresh LB allows to produce more umbelliferone than fresh M9 or than using the same LB or M9 from the beginning to the end of the experiment (Figure 2b). In these assays, 1 mM of p-coumaric acid was used as precursor. The improved results with fresh LB compared to when the same LB is used for both production phases (protein and coumarins) makes sense as new nutrients are available, and the toxicity is lower as the by-products produced by E. coli during the protein production phase are removed. As expected, the production in only M9 minimal medium was very low since the bacterial growth is very slow in this medium and the bacteria takes more time to adapt.

2.2. C2′H substrate Specificity

In order to confirm C2′H specificity, three different precursors were fed to the culture medium (LB + fresh LB): p-coumaric acid, caffeic acid, and ferulic acid (at 1 mM final concentration) that may be converted to umbelliferone, esculetin, and scopoletin, respectively (Figure 1). E. coli BL21 (DE3) carrying pAC-4CL1 and pET28GST-C2′H was able to produce the three different coumarins in different amounts (Figure 3). As expected, the production of umbelliferone and scopoletin was obtained in similar amounts, with umbelliferone being ~1.25-fold higher than scopoletin. Esculetin was also produced but in lower amounts, which is not in accordance with the literature [41] as it has been previously reported that IbF6′H2-2-1 (C2′H from this study) did not show any activity when caffeoyl-CoA was used as substrate in the in vitro reactions. It is possible that the activity was very low to be detected in the in vitro reactions. This was also verified by Yang et al. [32], who used two F6′H enzymes that produced esculetin in vivo and could not produce it in vitro [41]. Although the coumarins productions obtained in this study are lower than those reported by Yang et al. [32], thus suggesting that the enzyme from this study has lower activity towards the substrates, the proportion of coumarins obtained is very similar showing that they have a similar substrate specificity.

2.3. Heterologous Production of Coumarins from Tyrosine

To produce coumarins from the amino acid tyrosine, TAL from Rhodotorula glutinis, 4-coumarate 3-hydroxylase (C3H) from Saccharothrix espanaensis and caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) from Medicago sativa were used (Figure 1). All these genes were previously selected to construct caffeic acid and curcuminoids pathways by our research group [43,45,46]. TAL converts tyrosine to p-coumaric acid that is then converted to caffeic acid by C3H. 4CL1 is used to convert the hydroxycinnamic acids into their CoA esters and lastly, CCoAOMT converts caffeoyl-CoA to feruloyl-CoA. The combination of these genes allows to obtain all the three CoA esters (p-coumaroyl-CoA, caffeoyl-CoA, and feruloyl-CoA) that are finally converted to the respective coumarins by C2′H. In these experiments, tyrosine (3 mM) was used as a precursor. Depending on the combination of plasmids used, one, two, or three coumarins were produced (Figure 4). p-Coumaric acid was highly accumulated in all the experiments. This is consistent with previous results using TAL [43,44,45,46] as this enzyme is highly efficient. The other enzymes of the pathway are not as efficient. Nevertheless, it was possible to observe a small accumulation of caffeic and ferulic acid that justifies the low production of the respective coumarins. Umbelliferone was produced in higher amounts since its precursor is highly available. The production of scopoletin was relatively high taking into account the available amounts of ferulic acid, again demonstrating the high affinity of C2′H to feruloyl-CoA [41]. Comparing these results to other published studies, it is possible to conclude that the production of coumarins is in some cases lower which was expected as the production of coumarins directly from the hydroxycinnamic acids was already lower [32]. Then, the increase in the number of genes involved in the pathway also increased the metabolic burden leading to lower productions [44]. Nevertheless, the amount of umbelliferone produced was sufficient to continue our study and test the next step of the pathway catalyzed by a PT enzyme (Section 2.4).

2.4. The Prenyltransferase (PT) Step

The conversion of umbelliferone to demethylsuberosin catalyzed by a 6-PT is the entry point of the furanocoumarins pathway [14]. Plant PTs are known to be membrane-bound enzymes and, consequently, they are difficult to express in microorganisms. Therefore, its expression is not usually efficient in E. coli, with N. benthamiana often being used [35,36,37]. Recently, Bu et al. [34] showed that codon-optimized PT from Petroselinum crispum (PcPT) and PT2 from Pastinaca sativa (PsPT2) could be used to convert umbelliferone to demethylsuberosin and osthenol (angular furanocoumarin pathway), respectively, in Streptomyces xiamenensis and in E. coli. This was the first time that a plant PT enzyme was reported to be functional in a bacterium. Considering these results, we selected PcPT (codon-optimized for E. coli) to convert umbelliferone to demethylsuberosin in vivo. PcPT was cloned in pETDuet-TAL multiple cloning site 2 (MCS2) or in pCDFDuet-1 to be combined with the other plasmids (pAC-4CL1 and pET28GST-C2′H) and tyrosine was used as substrate. However, no demethylsuberosin production was observed. Therefore, p-coumaric acid and umbelliferone were also tested as precursors in vivo with no success. To rule out the idea that PcPT was not being expressed in enough amounts due to its presence in a medium-copy plasmid (pETDuet-1/pCDFDuet-1), we attempted the PcPT cloning in the high copy number pRSFDuet-1 plasmid and also in pET28GST-LIC to evaluate if GST-tag could improve the activity. PcPT activity was evaluated in vitro using these plasmids. However, again demethylsuberosin was not detected. In addition to PcPT, codon-optimized PT1 from P. sativa (PsPT1) known to convert umbelliferone to demethylsuberosin [35] (although not as efficiently as PcPT [36]) was also evaluated in vitro with no success. SDS-PAGE gels also did not reveal the desired proteins (data not shown) demonstrating that both PT enzymes were not being successfully or efficiently expressed in E. coli.
Plant PT enzymes, as membrane-bound proteins, contain transmembrane helices and are known to have transit peptides that guide them to transmembrane regions [48]. Therefore, to characterize Citrus limon PT using S. cerevisiae as a host, Munakata et al. [49] deleted 36 amino acids from the N-terminal region related to a transit peptide that was targeting the enzyme to the plastid. This strategy allowed them to successfully characterize the PT enzyme in yeast. It was suggested that the removal of the transit peptides can help the PT, the substrates, and the prenyl donors to be closer, improving the prenylation efficiency in microorganisms [48]. Hence, we decided to evaluate if this strategy could be used to improve the expression also in bacteria. We started by analyzing the PT sequences using several bioinformatic tools. The results were different depending on the tool used (Figures S2–S13). ChloroP (currently discontinued) did not predict any chloroplast transit peptide in PcPT. This was not expected as Karamat et al. [36] mentioned that ChloroP predicted a transit peptide of 48 aa targeting chloroplast. However, LocTree3 and BaCelLo predicted that the protein subcellular localization was the chloroplast/chloroplast membrane. iPSORT predicted a mitochondrial targeting peptide at N-terminal while TargetP and SignalP did not predict a signal peptide (mitochondrial or chloroplast). Since these results were not coherent, it was decided to consider a transited peptide with 48 aa. Bioinformatic analysis of P. sativa PT1 also presented contradictory results. ChloroP predicted that PsPT1 had a N-terminal chloroplast transit peptide with a length of 48 aa. LocTree3 and BaCelLo also pointed to a chloroplast transit peptide. However, iPSORT predicted a mitochondrial targeting peptide at N-terminal. TargetP also predicted a transfer peptide with a cleavage site at position 32–33 aa. However, SignalP-5.0 did not predict any signal peptide. Again, we decided to consider a transit peptide with 48 aa. Therefore, we deleted the first 144 bp from both genes. A second codon was added that codified alanine (Met-Ala) as it is a preferred second codon in E. coli, present in several highly expressed genes (e.g., β-galactosidase gene) [50]. GST was removed during the cloning in pET28GST-LIC to properly test the N-terminal modification. Unfortunately, the expression of truncated PTs did not show any improvement in in vivo and in vitro experiments, and in SDS-PAGE (data not shown).
Although the prenyl donor DMAPP availability was previously identified as a limiting step in the production of several aromatic compounds in microbes [48,51,52], we did not expect that it was so extreme that no production at all was observed. Indeed, Bu et al. [34] did not report any supplementation of DMAPP in the in vivo experiments with E. coli. Nevertheless, we repeated the in vivo experiments by expressing pMVA plasmid, containing several genes from the mevalonate (MVA) pathway responsible for the synthesis of DMAPP. The plasmid contained genes from Enterococcus faecalis (mvaS, hydroxymethylglutaryl-CoA synthase; mvaE, bifunctional acetoacetyl-CoA thiolase, and 3-hydroxy-3-methylglutaryl-CoA reductase), E. coli BL21 (idi, isopentenyl pyrophosphate—DMAPP isomerase) and S. cerevisiae (erg12, MVA kinase; erg8, phospho-MVA kinase; erg19, MVA pyrophosphate decarboxylase). Again, no improvement was observed in the in vivo experiments.

2.5. The Psoralen Synthase (PS) Step

PS is responsible for the conversion of marmesin to psoralen. This enzyme was first characterized in 2007 by Larbat et al. [38]. This enzyme belongs to cytochrome P450 (CYP450) family and is known to be difficult to express in E. coli as this organism does not have in its genome any CYP450 or CYP450 reductases (CPR) (redox partners to enhance CYP450 activity). Actually, preliminary tests in a previous work [38] suggest that PS enzyme from Ammi majus (CYP71AJ1) was not successfully expressed in E. coli and in S. cerevisiae without any protein modification. Nevertheless, in this study we selected PS enzyme from A. majus as from the few enzymes already characterized [38,39], it is the one that showed more activity towards marmesin. Codon-optimized PS was cloned in pRSFDuet-1 and pET28GST-LIC. The expression was evaluated in SDS-PAGE (Figure S14). No protein of interest was observed when expressed in pRSFDuet-1. In pET28GST-PS, it was possible to observe a protein highly expressed compared to the time 0 control. The expected protein band in this plasmid was expected around 82 kDa. However, the gels show 2 bands very close to each other around 70–60 kDa. The bands very close suggest that PS protein migrated as a doublet protein on SDS-PAGE gels. A small band around 30 kDa is also observed that can be related to proteolysis even though a protease inhibitor was used. The differences in size may be related to gel-shifting that is common [53] and apparently even more common for membrane proteins [54]. PS, as a CYP450, is a membrane-bound protein. Bioinformatic analysis using SignalP-5.0 showed a possible signal sequence of 30 aa and iPSORT predicted a mitochondrial targeting peptide. TMHMM-2.0 also predicted a transmembrane helix at 12–29 aa, which is in good agreement with the SignalP-5.0 prediction. LocTree3 predicts an endoplasmic reticulum localization. ChloroP did not predict any signal peptide, however, BaCelLo predicted a chloroplast localization which is contradictory. Bioinformatic results are presented in supplementary material (Figures S15–S21). Although some of the bioinformatic tools suggest a 29–30 aa signal sequence, Larbat et al. [38] considered the complete sequence of 37 aa preceding the highly conserved proline-rich region (37-PPSPP-43) known to link the N-terminal transmembrane sequence to the CYP450 catalytic domain [55]. Larbat et al. [38] replaced the 37 aa N-terminal for the CYP73A1 N-terminal in order for PS being successfully expressed in S. cerevisiae. Similar strategies have been used with success to express other CYP450 in yeast [35]. Additionally, this type of strategy has also been successfully used to express this type of enzymes in E. coli [56]. The truncation or replacement of CYP450 and CPR N-terminal regions led to an increase in CYP450 and CPR solubility and expression [57,58]. N-terminal has been frequently replaced by Met-Ala (as previously discussed) [58,59] or by a ‘universal’ N-terminal sequence (e.g., 8RP, MALLLAVF; 2C3, MAKKTSSKGK; 2B1, MAKKTSSKGKLPPG(PS)) [58,59,60,61,62,63]. More complex N-terminal sequences containing 84–214 aa have also been used (e.g., 28-tag, Sumo, MBP) [59,61].
After the literature review, it was decided to test the replacement of 37 aa N-terminal for MA (Met-Ala), 8RP, 2C3, and 28-tag (84 aa). The new N-terminal sequences were codon optimized for E. coli. Modified PS genes were cloned in pRSFDuet-1 and pET28GST-LIC. The SDS-PAGE results can be observed in Figures S22–S25. No PS protein was observed when MA N-terminal was used. However, high ∆37PS expression in insoluble phase was observed when the other three N-terminals were used, especially when cloned in pET28GST-LIC plasmid. Expression in pRSFDuet-1 was generally lower with exception of when 2C3 was used. In addition, when pRSFDuet-1 was used, the bands of interest were observed at the expected place (53–57 kDa). When expressed in pET28GST-LIC, the results were similar to the ones previously obtained before the N-terminal replacement (Figure S14), with a very intense band near the correct size, suggesting 2 bands very close to each other in the 8RP and 28tag cases, and one smaller band at around 30 kDa. We did not find any explanation for the presence of this lower band since in these cases the other band(s) had the correct size making the previously suggested proteolysis less probable.
As mentioned, E. coli does not contain a CPR enzyme. Therefore, CPR2 from A. thaliana [59], successfully used in the production of baicalein and scutellarein in E. coli, was selected for this work. This electron transfer partner was expressed in the supplied plasmid pETDuet-CPR2 in combination with other plasmids carrying PS variants (pET28GST-PS, pRSFDuet-PS, pET28-MA-∆37PS, pRSFDuet-MA-∆37PS, pET28–8RP-∆37PS, pRSFDuet-8RP-∆37PS, pET28-2C3-∆37PS, pRSFDuet-2C3-∆37PS, pET28-28tag-∆37PS, and pRSFDuet-28tag-∆37PS). The protein extracts were used in in vitro reactions to evaluate PS and PS variants activity. Activity of insoluble fractions was also tested as aggregation and inclusion body formation may result in only a moderate loss of activity [64]. Unfortunately, it was not possible to detect the product psoralen in any of the combinations tested.
PS protein was visible in the protein gels in very high amounts but in the insoluble phase. It is frequent that the overexpression of heterologous proteins using strong promoters can result in misfolding and aggregation of the heterologous proteins as inclusion bodies [65,66]. All this induces E. coli heat-shock response. However, this response is usually not enough for the high-speed production of difficult to fold proteins [65,66]. Therefore, the overexpression of proteins that assist in protein folding such as chaperones (e.g., GroEL/GroES) in combination with the protein(s) of interest has been used to improve protein folding, stability, and activity [67]. Therefore, GroEL and GroES chaperones were cloned in pCDFDuet-1, which later was combined with pETDuet-CPR2 and pET28GST/pRSFDuet-PS (and PS variants). The chaperones expression was evaluated in an SDS-PAGE gel (Figure S26). It is possible to observe that GroEL is highly expressed in both fractions. Since GroEL and PS variants are almost the same size (55–57 kDa) it is not possible to observe if the expression of PS variants improved in the soluble phase or not (Figures S26–S27). Afterwards, the PS variants’ activity was evaluated in vitro by expressing the enzymes with GroESL and CPR2. However, no psoralen was detected either using soluble or insoluble fractions. In the future, other strategies may be considered such as the refolding of purified inclusion bodies or the use of weaker promoters [66,68,69].

2.6. The Marmesin Synthase (MS) Step

The MS step was first described in 1988 by Hamerski and Matern [70] but only in 2021 was an enzyme able to perform this step discovered [40]. Villard et al. [40] identified an enzyme capable of converting demethylsuberosin to marmesin with very high affinity for the first time. This enzyme, CYP76F112, was identified in Ficus carica and also belongs to CYP450 family. The authors were able to successfully express and characterize this enzyme in S. cerevisiae with no modifications.
We started by analyzing the CYP76F112 protein sequence using bioinformatics tools. In the case of this MS protein, the bioinformatic results were in general very consistent (Figures S28–S34). iPSORT, TargetP and SignalP predicted that MS has a secretory signal peptide at its N-terminal with a cleavage site between position 25 and 26. LocTree3 pointed to a localization in the endoplasmic reticulum membrane and TMHMM also predicts a transmembrane helix at position 7 to 26. Considering these results, we decided to assess native CYP76F112 expression and a version without the first 25 aa and with 8RP N-terminal. The protein expression was evaluated using SDS-PAGE; howver, no protein was observed with any version, even when MS proteins were co-expressed with GroESL chaperones and CPR2 (data not shown). The in vitro reactions also did not show any production of marmesin.

3. Conclusions

As demonstrated in this work, producing furanocoumarins in E. coli is still a challenge. While coumarins have been successfully produced in this host in the past and also in this work, there was not enough information regarding the following steps of the pathway (PT, MS, and PS). In our study, despite several attempts it was not possible to confirm the PT, PS, or MS activity when the enzymes were expressed in E. coli. Therefore, new approaches are required. In the future, PTs from microbial sources might be a better option for use in biotechnological applications [71]. Contrarily to plant PTs, these microbial enzymes are soluble and present a very broad substrate promiscuity in vitro for prenyl acceptors. So far, several microbial aromatic PTs that are involved in the production of secondary metabolites in microbes have been identified and successfully expressed in E. coli [51,72,73,74,75]. Although these enzymes are not involved in the production of furanocoumarins, at least one [51] was shown to catalyze the prenylation of umbelliferone to 7-dimethylallyl-umbelliferone. These enzymes, whose activity towards furanocoumarins needs to be analyzed, and new ones that will soon be identified, are expected to present a more suitable approach for the reconstruction of the furanocoumarin’s pathway in E. coli. Actually, the same reasoning can be used regarding the MS step. This plant enzyme was poorly explored since it has only been recently identified [40]. So far, everything indicates that this enzyme, as PS, is not a suitable option to express in E. coli. Therefore, the recently identified gene XimD from S. xiamenensis [34] that was able to produce marmesin in vivo in S. xiamenensis and in E. coli might be a better option.
Although E. coli has been considered an alternative chassis to produce several plant secondary metabolites, other hosts might be more suitable to engineer the furanocoumarins pathway. The S. cerevisiae genetic toolbox is also highly developed [76] and, as an eukaryote, it can perform post-translational modifications and contains intracellular compartments similar to plants [22,77]. This makes it a more suitable host for the expression of plant enzymes, including CYP450. This host has been used to express several PS [38,39] and, recently, also MS [40]. Although it was reported that it may not be a suitable host to express PTs from the furanocoumarins pathway [35,36], it has been used to express other membrane-bound PTs [52], including through PT N-terminal modification [36]. This strategy using this host may, in the future, achieve better results than the ones obtained in this study with E. coli. In addition, S. cerevisiae has been recently used to produce coumarins using a biosynthetic pathway [33]. However, the productions obtained are very low compared to the ones from E. coli. Actually, since the first part of the pathway involving coumarin biosynthesis has been successfully implemented in E. coli, and S. cerevisiae appears as a more promising host for the expression of the other enzymes of the pathway, a co-culture or sequential culture strategy using both hosts might be considered in the future. Co-culture studies using both hosts for polyphenolic compound production have been reported [78]. This interesting strategy might be viable as umbelliferone is transported to the extracellular medium after production and it has been reported that S. cerevisiae is able to uptake coumarins [79].
In conclusion, this work comprised an in-depth study of several steps from furanocoumarins biosynthetic pathway and several limitations have been pointed out. Hence, in the future, to achieve an efficient production of furanocoumarins in heterologous microbial hosts, all the hypotheses herein discussed should be considered.

4. Materials and Methods

4.1. Strains, Plasmids, and Chemicals

E. coli NZY5α (NZYTech—MB00401) was used for cloning and plasmid propagation. E. coli BL21 (DE3) (NZYTech—MB006) was used for the expression of the heterologous pathway. All plasmids used in this study are described in Table 1. pETDuet-TAL, pETDuet-TAL-C3H, and pCDFDuet-CCoAOMT constructions were previously described [43,45]. pAC-4CL1 (Addgene #35947), pET28GST-LIC (Addgene #26101), pETDuet-CPR2, and pYeDP60-CYP76F112 were kindly provided by Dr. Claudia Schmidt-Dannert [80], Dr. Cheryl Arrowsmith, Dr. Yong Wang [59], and Dr. Alain Hehn [40], respectively. pMVA (Addgene #121149) [75] and pETDuet-PcPT-XimD-XimE (Addgene #172654) [34] were a gift from Dr. Min-Juan Xu. C2′H, PsPT1 and PS were codon-optimized using ATGenium algorithm and synthesized by NZYTech (Lisbon, Portugal) and provided in pNZY29 and pHTP0 plasmids. The DNA sequences of the codon-optimized C2′H, PsPT1 and PS genes are presented in Supplementary Table S1.
Super optimal broth with catabolite repression (SOC), lysogeny broth (LB) Miller medium, isopropyl β-D-1-thiogalactopyranoside (IPTG), and 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-Gal) were purchased from NZYTech. p-Coumaric acid, caffeic acid and tyrosine were obtained from Sigma-Aldrich; ferulic acid, umbelliferone and scopoletin from Acros; esculetin and psoralen from Alfa Aesar; and demethylsuberosin and marmesin from Chengdu Biopurify Phytochemicals (Sichuan, China). Glucose (Acros), NH4Cl, NaCl, CaCO3 (Panreac), Na2HPO4 (Chem-Lab, Zedelgem, Belgium), MgSO4, CoCl2, FeCl3, CuCl2, ZnCl2, nicotinic acid, KH2PO4 (Riel-deHaën), thiamine (Thermo Fisher Scientific, Loughborough, United Kingdom), NaMoO4, pyridoxine, H2BO3, folic acid, biotin (Merck, Kenilworth, New Jersey, USA), riboflavin, and pantothenic acid (Sigma Aldrich) were used to prepare the M9 minimal medium. Ampicillin (VWR), chloramphenicol, kanamycin (NZYTech), spectinomycin (Alfa Aesar, Kandel, Germany), and sucrose (Labkem) were used for strain selection. Acetonitrile, trifluoroacetic acid and ethyl acetate were purchased from Fisher-Scientific. DMAPP (Merck), dithiothreitol (DTT), Tris (Fisher BioReagents, Fair Lawn, New Jersey, USA), NaH2PO4 (Scharlau), protease inhibitor (NZYTech), nicotinamide adenine dinucleotide phosphate (NADPH) (Panreac, Barcelona, Spain), and MgCl2 (VWR) were supplemented in enzymatic assays.

4.2. Construction of Plasmids

C2′H was amplified from pNZY29-C2′H and cloned in pRSFDuet-1 and pET28GST plasmids. PcPT and ∆48PcPT were amplified from pETDuet-PcPT-XimD-XimE. PsPT1 and ∆48PsPT1 were amplified from pHTP0-PsPT1. PcPT was cloned in pETDuet-TAL, pCDFDuet-1, pRSFDuet-1 and pET28GST-LIC. PS and PS variants were amplified from pHTP0-PS. PsPT1, ∆48PcPT, ∆48PsPT1, PS, MA-∆37PS, 8RP-∆37PS, 2C3-∆37PS, and 28tag-∆37PS were cloned in pRSFDuet-1 and pET28GST-LIC. The GroESL operon was amplified from an E. coli BL21 colony using colony-PCR. Afterwards, the chaperones were cloned as an operon in pCDFDuet-1 MCS1. CYP76F112 and 8RP-∆25CYP76F112 were amplified from pYeDP60-CYP76F112 and cloned in pRSFDuet-1. All the primers (Metabion/Eurofins) used in this study are presented in Supplementary Table S2.
Plasmid DNA was isolated using NucleoSpin Plasmid Miniprep Kit (Macherey-Nagel, Düren, Germany). The genes were amplified using Phusion High Fidelity DNA Polymerase (Thermo Fisher Scientific) and purified from agarose using NucleoSpin Gel and PCR Clean-up Kit (Macherey-Nagel). DNA was quantified using NanoDrop One (Thermo) and digested with the appropriate restriction enzymes (Supplementary Table S2). After digestion, DNA was purified using NucleoSpin Gel and PCR Clean-up Kit. T4 DNA ligase (Thermo) was used for ligation. Chemical transformation was performed using E. coli NZY5α competent cells. The constructed plasmids were verified by colony-PCR (Speedy supreme NZYTaq 2x Green Master Mix, NZYTech), digestion and sequencing (Eurofins). After confirmation, the plasmids were transformed in E. coli BL21 (DE3) competent cells. All enzymes and kits were used according to the manufacturer’s instructions.

4.3. Coumarins/Furanocoumarins Production

LB culture medium was used for plasmid propagation in E. coli BL21 (DE3) and for inoculums (200 rpm, 37 °C). For coumarin production, cultures were grown at 37 °C in 50 mL LB (250 mL flasks) from an optical density at 600 nm (OD600) of 0.1 to 0.6. Afterwards, the heterologous protein expression was induced with IPTG (0.1 mM final concentration), and the culture was incubated at 26 °C for 5 h. Then, the cells were collected by centrifugation and suspended in coumarin production medium (LB or M9 minimal medium). In two cases, cells were maintained in the same LB used for protein production or grew from the beginning in M9. IPTG (0.1 mM final concentration), precursors (p-coumaric acid, 0.7 mM–1 mM; caffeic acid, 1 mM; ferulic acid, 1 mM; or tyrosine, 3 mM) and antibiotics (ampicillin, 100 µg/mL; kanamycin, 50 µg/mL; spectinomycin, 100 µg/mL; and/or chloramphenicol, 25 µg/mL) were supplemented at time 0 depending on the specific experiment. M9 minimal medium contained (per L): 40 g glucose, 3 g KH2PO4, 6 g Na2HPO4, 0.5 g NaCl, 1 g NH4Cl, 110 mg MgSO4, 15 mg CaCl2, 340 mg thiamine and 5 g CaCO3. Vitamins (12.2 mg nicotinic acid, 10.8 mg pantothenic acid, 2.8 mg pyridoxine, 0.84 mg riboflavin, 0.12 mg biotin, and 0.084 mg folic acid) and trace elements (54 mg FeCl3, 4 mg NaMoO4, 4 mg ZnCl2, 4 mg CoCl2, 2 mg CuCl2, and 1 mg H2BO3) were also supplemented to M9.
Assays targeting demethylsuberosin production in vivo using PT enzymes were performed as coumarins production assays.

4.4. Coumarin/Furanocoumarin Extraction

The whole broth was extracted as coumarins can be detected in significant amounts in the supernatant and inside the cells. Samples (2 mL whole broth) were mixed with 2 mL of ethyl acetate, vortexed and centrifuged (15,000× g, 5 min). The supernatant was concentrated by solvent evaporation, resuspended in 200 µL of acetonitrile and analyzed by ultra-high-performance liquid chromatography (UHPLC). The same method was used to evaluate demethylsuberosin production in vivo, but a higher volume of whole broth samples was extracted.

4.5. UHPLC Analysis

Hydroxycinnamic acids (p-coumaric acid, caffeic acid, and ferulic acid), coumarins (umbelliferone, esculetin, and scopoletin) and furanocoumarins (demethylsuberosin, marmesin, and psoralen) were quantified by UHPLC. Supernatant was used to quantify the hydroxycinnamic acids while coumarins and furanocoumarins were extracted as described in 4.4. Samples from in vitro reactions were directly analyzed without any extraction.
Quantification was performed using the Shimadzu Nexera-X2 system (Shimadzu Corporation, Kyoto, Japan) (LC-30AD pump unit, DGU-20A 5R degasser unit, CTO-20AC column oven, CBM-20A system controller, SPD-M20A detector unit, SIL-30AC autosampler unit) and a Kinetex 2.6 µm Polar C18 100 Å LC column (150 mm × 4.6 mm) (Phenomenex, Alcobendas, Spain). Water with 0.1% (v/v) of formic acid and pure acetonitrile were used as mobile phase A and B, respectively. The following gradient was used at a flow rate of 1 mL/min for the determination of hydroxycinnamic acids, coumarins and furanocoumarins: 10–80% mobile phase B for 16 min, 80–10% for 1 min, and 10% mobile phase B for an additional 3 min. p-Coumaric acid, caffeic acid, and ferulic acid were detected at 310 nm at the retention time of 6.4, 5.1 and 7.7 min, respectively. Umbelliferone, esculetin and scopoletin were detected at 324, 340, and 342 nm and at a retention time of 6.6, 4.9, and 7.5 min, respectively. Finally, demethylsuberosin, marmesin and psoralen were detected at 334, 335, and 294 nm and at a retention time of 12.5, 8.3 and 10.0 min, respectively.

4.6. Protein Analysis

E. coli BL21 (DE3) cells carrying the plasmids of interest were grown in LB at 37 °C to an OD600 of 0.6. Then, IPTG was added at a final concentration of 1 mM to induce protein expression, the temperature was decreased to 26 °C and the culture was incubated for 6 h. Samples (10 mL) were taken at 0 and 6 h of induction and centrifuged. The pellets were resuspended in 1 mL of Tris-HCl buffer (10 mM, pH 7.8) supplemented with protease inhibitor (NZYTech) and lysed by sonication on ice (35% amplitude, 3 s ON plus 9 s OFF for a total of 5 min ON) using a microtip probe linked to a Vibra-cell processor (Sonics). After sonication, samples were centrifuged and resuspended, and the protein of the soluble and insoluble fractions was quantified using the Bradford reagent (Panreac). The expression levels were evaluated through SDS-PAGE (4% stacking gel and 10% running gel). Soluble and insoluble protein fractions were mixed with 2 x sample buffer (65.8 mM Tris–HCl pH 6.8, 26.3% glycerol, 0.01% bromophenol blue, 2.1% SDS, 5% β-mercaptoethanol) and denaturated at 95 °C in the heating block for 5 min. The protein ladders used were NZYColour Protein Marker II (NZYTech) and Blue Prestained Protein Standard—Broad Range (NEB). The gel was stained with Coomassie Blue R-250 for 15 min and de-stained using distilled water.

4.7. Enzymatic Assays

The protein extracts were obtained and quantified as described in Section 4.6. Afterwards, soluble and insoluble fractions were used to evaluate PT (Section 4.7.1), PS (Section 4.7.2) and MS (Section 4.7.3) activity in vitro. The reactions occurred at 27 or 30 °C with or without agitation for 4 h. Samples were taken at times 0, 2, and 4 h and analyzed as described in Section 4.5.

4.7.1. PT In Vitro Reactions

PT enzymatic assays were performed based on [35,36,49]. One milliliter reactions were prepared in 100 mM Tris-HCl pH 8 (supplemented with protease inhibitor) containing 200 µM DMAPP, 200 µM umbelliferone, 10 mM MgCl2, and 1, 10, or 100 µg protein extract (soluble or insoluble phase). A quantity of 1 mM DTT was added in some tests.

4.7.2. PS In Vitro Reactions

PS enzymatic assays were performed based on [38,39]. One milliliter reactions were prepared in 100 mM sodium phosphate pH 7 (supplemented with protease inhibitor) containing 1 mM NADPH, 0.1 mM marmesin and 1, 10, or 100 µg protein extract (soluble or insoluble phase).

4.7.3. MS In Vitro Reactions

MS enzymatic assays were performed based on [40,70]. One milliliter reactions were prepared in 100 mM sodium phosphate pH 7 (supplemented with protease inhibitor) containing 200 µM NADPH, 50 µM demethylsuberosin, and 1, 10, or 100 µg protein extract (soluble or insoluble phase).

4.8. Bioinformatics

Enzymes transit peptides and subcellular localization were predicted using ChloroP (https://services.healthtech.dtu.dk/service.php?ChloroP-1.1) (discontinued, accessed on 1 May 2020), iPSORT (http://ipsort.hgc.jp/) [81] (accessed on 1 May 2020), SignalP-5.0 (https://services.healthtech.dtu.dk/service.php?SignalP-5.0) [82] (accessed on 1 May 2020), TargetP-2.0 (https://services.healthtech.dtu.dk/service.php?TargetP-2.0) [83] (accessed on 1 May 2020), TMHMM-2.0 (https://services.healthtech.dtu.dk/service.php?TMHMM-2.0) [84] (accessed on 1 May 2020), Loctree3 (https://rostlab.org/services/loctree3/) [85] (accessed on 1 May 2020), and BaCelLo (http://gpcr.biocomp.unibo.it/bacello/) [86] (accessed from 1 May 2020 to 1 June 2022).

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/molecules27217230/s1, Table S1: Gene sequence of p-coumaroyl-CoA 2′-hydroxylase (C2′H), prenyltransferase 1 (PT1) and psoralen synthase (PS) with codon optimization; Table S2: Primers for PCR amplification of the genes of the biosynthetic pathway (forward (FW) and reverse (REV) primers) and sequencing; Figure S1: Protein SDS gels showing p-coumaroyl-CoA 2′-hydroxylase (C2′H) expression using different plasmids (pRSFDuet-C2′H and pET28GST-C2′H) at time zero (t0) of induction and after 6 h (t6) of induction. C2′H is expected around 40.41 kDa and GST + C2′H at 67.65 kDa. Arrows indicate where it is possible to observe the bands of interest. M: marker (NZYColour Protein Marker II—NZYTech); Figure S2: ChloroP results for PT from Petroselinum crispum; Figure S3: LocTree3 results for PT from Petroselinum crispum; Figure S4: BaCelLo results for PT from Petroselinum crispum; Figure S5: iPSORT results for PT from Petroselinum crispum; Figure S6: TargetP-2.0 results for PT from Petroselinum crispum; Figure S7: SignalP-5.0 results for PT from Petroselinum crispum; Figure S8: ChloroP results for PT1 from Pastinaca sativa; Figure S9: LocTree3 results for PT1 from Pastinaca sativa; Figure S10: BaCelLo results for PT1 from Pastinaca sativa; Figure S11: iPSORT results for PT1 from Pastinaca sativa; Figure S12: TargetP results for PT1 from Pastinaca sativa; Figure S13: SignalP-5.0 results for PT1 from Pastinaca sativa; Figure S14: Protein SDS gels showing psoralen synthase (PS) expression using different plasmids (pRSFDuet-PS and pET28GST-PS) at time zero (t0) of induction and after 6 h (t6) of induction. PS is expected around 55.9 kDa and GST+PS at 82.08 kDa. Arrows indicate where it is possible to observe the bands of interest. M: marker (NZYColour Protein Marker II—NZYTech); Figure S15: ChloroP results for PS from Ammi majus; Figure S16: iPSORT results for PS from Ammi majus; Figure S17: TMHMM-2.0 results for PS from Ammi majus; Figure S18: LocTree3 results for PS from Ammi majus; Figure S19: BaCelLo results for PS from Ammi majus; Figure S20: TargetP-2.0 results for PS from Ammi majus; Figure S21: SignalP-5.0 results for PS from Ammi majus; Figure S22: Protein SDS gels showing ∆37 aa psoralen synthase (PS) expression with MA N-terminal using different plasmids (pRSFDuet-MA-∆37PS and pET28-MA-∆37PS) at time zero (t0) of induction and after 6 h (t6) of induction. PS is expected around 52–54 kDa. M: marker (Blue Prestained Protein Standard, Broad Range—NEB); Figure S23: Protein SDS gels showing ∆37 aa psoralen synthase (PS) expression with 8RP N-terminal using different plasmids (pRSFDuet-8RP-∆37PS and pET28-8RP-∆37PS) at time zero (t0) of induction and after 6 h (t6) of induction. PS is expected around 53–55 kDa. Arrows indicate where it is possible to observe the bands of interest. M: marker (Blue Prestained Protein Standard, Broad Range—NEB); Figure S24: Protein SDS gels showing ∆37 aa psoralen synthase (PS) expression with 2C3 N-terminal using different plasmids (pRSFDuet-2C3-∆37PS and pET28-2C3-∆37PS) at time zero (t0) of induction and after 6 h (t6) of induction. PS is expected around 53–55 kDa. Arrows indicate where it is possible to observe the bands of interest. M: marker (Blue Prestained Protein Standard, Broad Range—NEB); Figure S25: Protein SDS gels showing ∆37 aa psoralen synthase (PS) expression with 28tag N-terminal using different plasmids (pRSFDuet-28tag-∆37PS and pET28-28tag-∆37PS) at time zero (t0) of induction and after 6 h (t6) of induction. PS is expected around 55–57 kDa. Arrows indicate where it is possible to observe the bands of interest. M: marker (NZYColour Protein Marker II—NZYTech); Figure S26: Protein SDS gels showing GroESL expression alone and when combined with cytochrome P450 reductase (CPR2) and 8RP-∆37PS at time zero (t0) of induction and after 6 h (t6) of induction. GroEL, GroES, CPR2 and 8RP-∆37PS are expected at round 57, 10, 79, and 55 kDa, respectively. Arrows indicate where it is possible to observe the bands of interest. M: marker (NZYColour Protein Marker II—NZYTech); Figure S27: Protein SDS gels showing GroESL expression combined with cytochrome P450 reductase (CPR2) and 2C3-∆37PS or 28tag-∆37PS at time zero (t0) of induction and after 6 h (t6) of induction. GroEL, GroES, CPR2, 2C3-∆37PS, and 28tag-∆37PS are expected at round 57, 10, 79, 55, and 56 kDa, respectively. Arrows indicate where it is possible to observe the bands of interest. M: marker (NZYColour Protein Marker II—NZYTech); Figure S28: ChloroP results for MS from Ficus carica; Figure S29: LocTree3 results for MS from Ficus carica; Figure S30: BaCelLo results for MS from Ficus carica; Figure S31: iPSORT results for MS from Ficus carica; Figure S32: TargetP-2.0 results for MS from Ficus carica; Figure S33: SignalP-5.0 results for MS from Ficus carica; Figure S34: TMHMM-2.0 results for MS from Ficus carica.

Author Contributions

Conceptualization, J.L.R. and L.R.R.; methodology, J.L.R. and D.G.; validation, J.L.R. and D.G.; formal analysis, J.L.R.; investigation, J.L.R.; writing—original draft preparation, J.L.R.; writing—review and editing, D.G. and L.R.R.; funding acquisition, L.R.R. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Portuguese Foundation for Science and Technology (FCT) under the scope of the strategic funding of UIDB/04469/2020 unit, and by LABBELS—Associate Laboratory in Biotechnology, Bioengineering and Microelectromechanical Systems, LA/P/0029/2020. D.G. acknowledges her grant SFRH/BD/04433/2020.

Acknowledgments

We acknowledge the support and expertise of Alain Hehn (University of Lorraine) regarding the furanocoumarin pathway.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, H.-J.; Jin, B.-R.; An, H.-J. Umbelliferone ameliorates benign prostatic hyperplasia by inhibiting cell proliferation and G1/S phase cell cycle progression through regulation of STAT3/E2F1 axis. Int. J. Mol. Sci. 2021, 22, 9019. [Google Scholar] [CrossRef] [PubMed]
  2. Cruz, L.F.; de Figueiredo, G.F.; Pedro, L.P.; Amorin, Y.M.; Andrade, J.T.; Passos, T.F.; Rodrigues, F.F.; Souza, I.L.A.; Gonçalves, T.P.R.; Lima, L.A.R.D.S.; et al. Umbelliferone (7-hydroxycoumarin): A non-toxic antidiarrheal and antiulcerogenic coumarin. Biomed. Pharmacother. 2020, 129, 110432. [Google Scholar] [CrossRef] [PubMed]
  3. Mahmoud, A.M.; Hozayen, W.G.; Hasan, I.H.; Shaban, E.; Bin-Jumah, M. Umbelliferone ameliorates CCl4-induced liver fibrosis in rats by upregulating PPARγ and attenuating oxidative stress, inflammation, and TGF-β1/Smad3 signaling. Inflammation 2019, 42, 1103–1116. [Google Scholar] [CrossRef]
  4. Zagaja, M.; Zagaja, A.; Szala-Rycaj, J.; Szewczyk, A.; Lemieszek, M.K.; Raszewski, G.; Andres-Mach, M. Influence of umbelliferone on the anticonvulsant and neuroprotective activity of selected antiepileptic drugs: An in vivo and in vitro study. Int. J. Mol. Sci. 2022, 23, 3492. [Google Scholar] [CrossRef] [PubMed]
  5. Zhang, L.; Xie, Q.; Li, X. Esculetin: A review of its pharmacology and pharmacokinetics. Phyther. Res. 2022, 36, 279–298. [Google Scholar] [CrossRef]
  6. Wang, Q.-H.; Qin, S.-W.; Jiang, J.-G. Improvement effects of esculetin on the formation and development of atherosclerosis. Biomed. Pharmacother. 2022, 150, 113001. [Google Scholar] [CrossRef]
  7. Lemos, A.S.O.; Florêncio, J.R.; Pinto, N.C.C.; Campos, L.M.; Silva, T.P.; Grazul, R.M.; Pinto, P.F.; Tavares, G.D.; Scio, E.; Apolônio, A.C.M.; et al. Antifungal activity of the natural coumarin scopoletin against planktonic cells and biofilms from a multidrug-resistant Candida tropicalis strain. Front. Microbiol. 2020, 11, 1525. [Google Scholar] [CrossRef] [PubMed]
  8. Napiroon, T.; Bacher, M.; Balslev, H.; Tawaitakham, K.; Santimaleeworagun, W.; Vajrodaya, S. Scopoletin from Lasianthus lucidus Blume (Rubiaceae): A potential antimicrobial against multidrug-resistant Pseudomonas aeruginosa. J. Appl. Pharm. Sci. 2018, 8, 1–6. [Google Scholar] [CrossRef] [Green Version]
  9. Parama, D.; Girisa, S.; Khatoon, E.; Kumar, A.; Alqahtani, M.S.; Abbas, M.; Sethi, G.; Kunnumakkara, A.B. An overview of the pharmacological activities of scopoletin against different chronic diseases. Pharmacol. Res. 2022, 179, 106202. [Google Scholar] [CrossRef]
  10. Kokron, O.; Maca, S.; Gasser, G.; Schmidt, P.R. Cimetidine and coumarin therapy of renal cell carcinoma. Oncology 1991, 48, 102–106. [Google Scholar] [CrossRef]
  11. Wang, C.-J.; Hsieh, Y.-J.; Chu, C.-Y.; Lin, Y.-L.; Tseng, T.-H. Inhibition of cell cycle progression in human leukemia HL-60 cells by esculetin. Cancer Lett. 2002, 183, 163–168. [Google Scholar] [CrossRef]
  12. Mohler, J.L.; Gomella, L.G.; Crawford, E.D.; Glode, L.M.; Zippe, C.D.; Fair, W.R.; Marshall, M.E. Phase II evaluation of coumarin (1,2-benzopyrone) in metastatic prostatic carcinoma. Prostate 1992, 20, 123–131. [Google Scholar] [CrossRef] [PubMed]
  13. Lim, G.B. Warfarin: From rat poison to clinical use. Nat. Rev. Cardiol. 2017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Rodrigues, J.L.; Rodrigues, L.R. Biosynthesis and heterologous production of furanocoumarins: Perspectives and current challenges. Nat. Prod. Rep. 2021, 38, 869–879. [Google Scholar] [CrossRef]
  15. Ling, T.C.; Clayton, T.H.; Crawley, J.; Exton, L.S.; Goulden, V.; Ibbotson, S.; McKenna, K.; Mohd Mustapa, M.F.; Rhodes, L.E.; Sarkany, R. British Association of Dermatologists and British Photodermatology Group guidelines for the safe and effective use of psoralen–ultraviolet A therapy 2015. Br. J. Dermatol. 2016, 174, 24–55. [Google Scholar] [CrossRef] [Green Version]
  16. Jagasia, M.; Scheid, C.; Socié, G.; Ayuk, F.A.; Tischer, J.; Donato, M.L.; Bátai, Á.; Chen, H.; Chen, S.-C.; Chin, T.; et al. Randomized controlled study of ECP with methoxsalen as first-line treatment of patients with moderate to severe cGVHD. Blood Adv. 2019, 3, 2218–2229. [Google Scholar] [CrossRef] [Green Version]
  17. Issa, M.Y.; Elshal, M.F.; Fathallah, N.; Abdelkawy, M.A.; Bishr, M.; Salama, O.; Abulfadl, Y.S. Potential anticancer activity of the furanocoumarin derivative xanthotoxin isolated from Ammi majus L. fruits: In vitro and in silico studies. Molecules 2022, 27, 943. [Google Scholar] [CrossRef]
  18. Steinack, C.; Robinson, C.A.; Nägeli, M.; Inci, I.; Benden, C. ECP as additional immunomodulation in idiopathic hyperammonemia and recurrent hypercapnic respiratory failure early post lung transplantation. J. Clin. Apher. 2021, 36, 186–188. [Google Scholar] [CrossRef]
  19. Meepagala, K.M.; Bracken, A.K.; Fronczek, F.R.; Johnson, R.D.; Wedge, D.E.; Duke, S.O. Furanocoumarin with phytotoxic activity from the leaves of Amyris elemifera (Rutaceae). ACS Omega 2021, 6, 401–407. [Google Scholar] [CrossRef]
  20. Britto, I.O.; Araújo, S.H.C.; Toledo, P.F.S.; Lima, G.D.A.; Salustiano, I.V.; Alves, J.R.; Mantilla-Afanador, J.G.; Kohlhoff, M.; Oliveira, E.E.; Leite, J.P. V Potential of Ficus carica extracts against Euschistus heros: Toxicity of major active compounds and selectivity against beneficial insects. Pest Manag. Sci. 2021, 77, 4638–4647. [Google Scholar] [CrossRef]
  21. Nebo, L.; Varela, R.M.; Molinillo, J.M.G.; Sampaio, O.M.; Severino, V.G.P.; Cazal, C.M.; Fernandes, M.F.D.G.; Fernandes, J.B.; Macías, F.A. Phytotoxicity of alkaloids, coumarins and flavonoids isolated from 11 species belonging to the Rutaceae and Meliaceae families. Phytochem. Lett. 2014, 8, 226–232. [Google Scholar] [CrossRef]
  22. Rainha, J.; Gomes, D.; Rodrigues, L.R.; Rodrigues, J.L. Synthetic biology approaches to engineer Saccharomyces cerevisiae towards the industrial production of valuable polyphenolic compounds. Life 2020, 10, 56. [Google Scholar] [CrossRef] [PubMed]
  23. Rodrigues, J.L.; Prather, K.L.J.; Kluskens, L.D.; Rodrigues, L.R. Heterologous production of curcuminoids. Microbiol. Mol. Biol. Rev. 2015, 79, 39–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rodrigues, J.L.; Ferreira, D.; Rodrigues, L.R. Synthetic biology strategies towards the development of new bioinspired technologies for medical applications. In Bioinspired Materials for Medical Applications; Mohapatra, S., Ranjan, S., Dasgupta, N., Mishra, R., Thomas, S., Eds.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 451–497. ISBN 9780081007464. [Google Scholar]
  25. Rodrigues, J.L.; Rodrigues, L.R. Synthetic biology: Perspectives in industrial biotechnology. In Current Developments in Biotechnology and Bioengineering: Foundations of Biotechnology and Bioengineering; Pandey, A., Teixeira, J., Eds.; Elsevier: Amsterdam, The Netherlands, 2017; pp. 239–269. ISBN 9780444636799. [Google Scholar]
  26. Sun, J.; Sun, W.; Zhang, G.; Lv, B.; Li, C. High efficient production of plant flavonoids by microbial cell factories: Challenges and opportunities. Metab. Eng. 2022, 70, 143–154. [Google Scholar] [CrossRef] [PubMed]
  27. Li, H.; Lyv, Y.; Zhou, S.; Yu, S.; Zhou, J. Microbial cell factories for the production of flavonoids–barriers and opportunities. Bioresour. Technol. 2022, 360, 127538. [Google Scholar] [CrossRef]
  28. Paddon, C.J.; Keasling, J.D. Semi-synthetic artemisinin: A model for the use of synthetic biology in pharmaceutical development. Nat. Rev. Microbiol. 2014, 12, 355–367. [Google Scholar] [CrossRef]
  29. Rainha, J.; Rodrigues, J.L.; Faria, C.; Rodrigues, L.R. Curcumin biosynthesis from ferulic acid by engineered Saccharomyces cerevisiae. Biotechnol. J. 2021, 17, 2100400. [Google Scholar] [CrossRef]
  30. Lin, Y.; Sun, X.; Yuan, Q.; Yan, Y. Combinatorial biosynthesis of plant-specific coumarins in bacteria. Metab. Eng. 2013, 18, 69–77. [Google Scholar] [CrossRef]
  31. Zhao, Y.; Jian, X.; Wu, J.; Huang, W.; Huang, C.; Luo, J.; Kong, L. Elucidation of the biosynthesis pathway and heterologous construction of a sustainable route for producing umbelliferone. J. Biol. Eng. 2019, 13, 44. [Google Scholar] [CrossRef] [Green Version]
  32. Yang, S.-M.; Shim, G.Y.; Kim, B.-G.; Ahn, J.-H. Biological synthesis of coumarins in Escherichia coli. Microb. Cell Fact. 2015, 14, 65. [Google Scholar] [CrossRef]
  33. Zhao, C.-H.; Zhang, R.-K.; Qiao, B.; Li, B.-Z.; Yuan, Y.-J. Engineering budding yeast for the production of coumarins from lignin. Biochem. Eng. J. 2020, 160, 107634. [Google Scholar] [CrossRef]
  34. Bu, X.; He, B.; Weng, J.; Jiang, C.; Zhao, Y.-L.; Li, S.-M.; Xu, J.; Xu, M. Constructing microbial hosts for the production of benzoheterocyclic derivatives. ACS Synth. Biol. 2020, 9, 2282–2290. [Google Scholar] [CrossRef] [PubMed]
  35. Munakata, R.; Olry, A.; Karamat, F.; Courdavault, V.; Sugiyama, A.; Krieger, C.; Silie, P.; Foureau, E.; Papon, N.; Grosjean, J. Molecular evolution of parsnip (Pastinaca sativa) membrane-bound prenyltransferases for linear and/or angular furanocoumarin biosynthesis. New Phytol. 2016, 211, 332–344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Karamat, F.; Olry, A.; Munakata, R.; Koeduka, T.; Sugiyama, A.; Paris, C.; Hehn, A.; Bourgaud, F.; Yazaki, K. A coumarin-specific prenyltransferase catalyzes the crucial biosynthetic reaction for furanocoumarin formation in parsley. Plant J. 2014, 77, 627–638. [Google Scholar] [CrossRef]
  37. Munakata, R.; Kitajima, S.; Nuttens, A.; Tatsumi, K.; Takemura, T.; Ichino, T.; Galati, G.; Vautrin, S.; Bergès, H.; Grosjean, J. Convergent evolution of the UbiA prenyltransferase family underlies the independent acquisition of furanocoumarins in plants. New Phytol. 2020, 225, 2166–2182. [Google Scholar] [CrossRef] [PubMed]
  38. Larbat, R.; Kellner, S.; Specker, S.; Hehn, A.; Gontier, E.; Hans, J.; Bourgaud, F.; Matern, U. Molecular cloning and functional characterization of psoralen synthase, the first committed monooxygenase of furanocoumarin biosynthesis. J. Biol. Chem. 2007, 282, 542–554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Larbat, R.; Hehn, A.; Hans, J.; Schneider, S.; Jugdé, H.; Schneider, B.; Matern, U.; Bourgaud, F. Isolation and functional characterization of CYP71AJ4 encoding for the first P450 monooxygenase of angular furanocoumarin biosynthesis. J. Biol. Chem. 2009, 284, 4776–4785. [Google Scholar] [CrossRef] [Green Version]
  40. Villard, C.; Munakata, R.; Kitajima, S.; Velzen, R.; Schranz, M.E.; Larbat, R.; Hehn, A. A new P450 involved in the furanocoumarin pathway underlies a recent case of convergent evolution. New Phytol. 2021, 231, 1923–1939. [Google Scholar] [CrossRef]
  41. Matsumoto, S.; Mizutani, M.; Sakata, K.; Shimizu, B.-I. Molecular cloning and functional analysis of the ortho-hydroxylases of p-coumaroyl coenzyme A/feruloyl coenzyme A involved in formation of umbelliferone and scopoletin in sweet potato, Ipomoea batatas (L.) Lam. Phytochemistry 2012, 74, 49–57. [Google Scholar] [CrossRef]
  42. Kai, K.; Mizutani, M.; Kawamura, N.; Yamamoto, R.; Tamai, M.; Yamaguchi, H.; Sakata, K.; Shimizu, B. Scopoletin is biosynthesized via ortho-hydroxylation of feruloyl CoA by a 2-oxoglutarate-dependent dioxygenase in Arabidopsis thaliana. Plant J. 2008, 55, 989–999. [Google Scholar] [CrossRef]
  43. Rodrigues, J.L.; Araújo, R.G.; Prather, K.L.J.; Kluskens, L.D.; Rodrigues, L.R. Production of curcuminoids from tyrosine by a metabolically engineered Escherichia coli using caffeic acid as an intermediate. Biotechnol. J. 2015, 10, 599–609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Rodrigues, J.L.; Gomes, D.; Rodrigues, L.R. A combinatorial approach to optimize the production of curcuminoids from tyrosine in Escherichia coli. Front. Bioeng. Biotechnol. 2020, 8, 59. [Google Scholar] [CrossRef] [PubMed]
  45. Rodrigues, J.L.; Araújo, R.G.; Prather, K.L.J.; Kluskens, L.D.; Rodrigues, L.R. Heterologous production of caffeic acid from tyrosine in Escherichia coli. Enzyme Microb. Technol. 2015, 71, 36–44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Rodrigues, J.L.; Couto, M.R.; Araújo, R.G.; Prather, K.L.J.; Kluskens, L.; Rodrigues, L.R. Hydroxycinnamic acids and curcumin production in engineered Escherichia coli using heat shock promoters. Biochem. Eng. J. 2017, 125, 41–49. [Google Scholar] [CrossRef] [Green Version]
  47. Couto, M.R.; Rodrigues, J.L.; Rodrigues, L.R. Optimization of fermentation conditions for the production of curcumin by engineered Escherichia coli. J. R. Soc. Interface 2017, 14, 20170470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. de Bruijn, W.J.C.; Levisson, M.; Beekwilder, J.; van Berkel, W.J.H.; Vincken, J.-P. Plant aromatic prenyltransferases: Tools for microbial cell factories. Trends Biotechnol. 2020, 38, 917–934. [Google Scholar] [CrossRef]
  49. Munakata, R.; Inoue, T.; Koeduka, T.; Karamat, F.; Olry, A.; Sugiyama, A.; Takanashi, K.; Dugrand, A.; Froelicher, Y.; Tanaka, R. Molecular cloning and characterization of a geranyl diphosphate-specific aromatic prenyltransferase from lemon. Plant Physiol. 2014, 166, 80–90. [Google Scholar] [CrossRef] [Green Version]
  50. Looman, A.C.; Bodlaender, J.; Comstock, L.J.; Eaton, D.; Jhurani, P.; de Boer, H.A.; van Knippenberg, P.H. Influence of the codon following the AUG initiation codon on the expression of a modified lacZ gene in Escherichia coli. EMBO J. 1987, 6, 2489–2492. [Google Scholar] [CrossRef]
  51. He, J.; Hu, Z.; Dong, Z.; Li, B.; Chen, K.; Shang, Z.; Zhang, M.; Qiao, X.; Ye, M. Enzymatic O-prenylation of diverse phenolic compounds by a permissive O-prenyltransferase from the medicinal mushroom Antrodia camphorata. Adv. Synth. Catal. 2020, 362, 528–532. [Google Scholar] [CrossRef]
  52. Levisson, M.; Araya-Cloutier, C.; De Bruijn, W.J.C.; Van Der Heide, M.; López, J.M.S.; Daran, J.-M.; Vincken, J.-P.; Beekwilder, J. Toward developing a yeast cell factory for the production of prenylated flavonoids. J. Agric. Food Chem. 2019, 67, 13478–13486. [Google Scholar] [CrossRef]
  53. Shirai, A.; Matsuyama, A.; Yashiroda, Y.; Hashimoto, A.; Kawamura, Y.; Arai, R.; Komatsu, Y.; Horinouchi, S.; Yoshida, M. Global analysis of gel mobility of proteins and its use in target identification. J. Biol. Chem. 2008, 283, 10745–10752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Rath, A.; Glibowicka, M.; Nadeau, V.G.; Chen, G.; Deber, C.M. Detergent binding explains anomalous SDS-PAGE migration of membrane proteins. Proc. Natl. Acad. Sci. USA 2009, 106, 1760–1765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Chen, C.-D.; Kemper, B. Different structural requirements at specific proline residue positions in the conserved proline-rich region of cytochrome P450 2C2. J. Biol. Chem. 1996, 271, 28607–28611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Barnes, H.J.; Arlotto, M.P.; Waterman, M.R. Expression and enzymatic activity of recombinant cytochrome P450 17 alpha-hydroxylase in Escherichia coli. Proc. Natl. Acad. Sci. USA 1991, 88, 5597–5601. [Google Scholar] [CrossRef] [Green Version]
  57. Stahlhut, S.G.; Siedler, S.; Malla, S.; Harrison, S.J.; Maury, J.; Neves, A.R.; Forster, J. Assembly of a novel biosynthetic pathway for production of the plant flavonoid fisetin in Escherichia coli. Metab. Eng. 2015, 31, 84–93. [Google Scholar] [CrossRef] [Green Version]
  58. Tamaki, S.; Yagi, M.; Nishihata, Y.; Yamaji, H.; Shigeri, Y.; Uno, T.; Imaishi, H. Modification of N-terminal amino acids of fungal benzoate hydroxylase (CYP53A15) for the production of p-hydroxybenzoate and optimization of bioproduction conditions in Escherichia coli. J. Microbiol. Biotechnol. 2018, 28, 439–447. [Google Scholar] [CrossRef] [Green Version]
  59. Li, J.; Tian, C.; Xia, Y.; Mutanda, I.; Wang, K.; Wang, Y. Production of plant-specific flavones baicalein and scutellarein in an engineered E. coli from available phenylalanine and tyrosine. Metab. Eng. 2019, 52, 124–133. [Google Scholar] [CrossRef]
  60. Dimaano, N.G.; Yamaguchi, T.; Fukunishi, K.; Tominaga, T.; Iwakami, S. Functional characterization of cytochrome P450 CYP81A subfamily to disclose the pattern of cross-resistance in Echinochloa phyllopogon. Plant Mol. Biol. 2020, 102, 403–416. [Google Scholar] [CrossRef]
  61. Christensen, U.; Vazquez-Albacete, D.; Søgaard, K.M.; Hobel, T.; Nielsen, M.T.; Harrison, S.J.; Hansen, A.H.; Møller, B.L.; Seppälä, S.; Nørholm, M.H.H. De-bugging and maximizing plant cytochrome P450 production in Escherichia coli with C-terminal GFP fusions. Appl. Microbiol. Biotechnol. 2017, 101, 4103–4113. [Google Scholar] [CrossRef]
  62. Rouck, J.E.; Biggs, B.W.; Kambalyal, A.; Arnold, W.R.; De Mey, M.; Ajikumar, P.K.; Das, A. Heterologous expression and characterization of plant Taxadiene-5α-Hydroxylase (CYP725A4) in Escherichia coli. Protein Expr. Purif. 2017, 132, 60–67. [Google Scholar] [CrossRef]
  63. Biggs, B.W.; Lim, C.G.; Sagliani, K.; Shankar, S.; Stephanopoulos, G.; De Mey, M.; Ajikumar, P.K. Overcoming heterologous protein interdependency to optimize P450-mediated Taxol precursor synthesis in Escherichia coli. Proc. Natl. Acad. Sci. USA 2016, 113, 3209–3214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. García-Fruitós, E.; González-Montalbán, N.; Morell, M.; Vera, A.; Ferraz, R.M.; Arís, A.; Ventura, S.; Villaverde, A. Aggregation as bacterial inclusion bodies does not imply inactivation of enzymes and fluorescent proteins. Microb. Cell Fact. 2005, 4, 27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Li, Z.; Nimtz, M.; Rinas, U. Global proteome response of Escherichia coli BL21 to production of human basic fibroblast growth factor in complex and defined medium. Eng. Life Sci. 2017, 17, 881–891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Rodrigues, J.L.; Rodrigues, L.R. Potential applications of the Escherichia coli heat shock response in synthetic biology. Trends Biotechnol. 2018, 36, 186–198. [Google Scholar] [CrossRef] [Green Version]
  67. Lee, A.R.; Li, L.; Shin, S.-Y.; Moon, J.S.; Eom, H.-J.; Han, N.S. Soluble expression of the fucosyltransferase gene from Helicobacter pylori in Escherichia coli by co-expression of molecular chaperones. Microbiol. Biotechnol. Lett. 2015, 43, 212–218. [Google Scholar] [CrossRef]
  68. Vallejo, L.F.; Rinas, U. Strategies for the recovery of active proteins through refolding of bacterial inclusion body proteins. Microb. Cell Fact. 2004, 3, 11. [Google Scholar] [CrossRef] [Green Version]
  69. Rodrigues, J.L.; Sousa, M.; Prather, K.L.J.; Kluskens, L.D.; Rodrigues, L.R. Selection of Escherichia coli heat shock promoters toward their application as stress probes. J. Biotechnol. 2014, 188, 61–71. [Google Scholar] [CrossRef] [Green Version]
  70. Hamerski, D.; Matern, U. Elicitor-induced biosynthesis of psoralens in Ammi majus L. suspension cultures. Eur. J. Biochem. 1988, 171, 369–375. [Google Scholar] [CrossRef]
  71. Gomes, D.; Rodrigues, L.R.; Rodrigues, J.L. Perspectives on the design of microbial cell factories to produce prenylflavonoids. Int. J. Food Microbiol. 2022, 367, 109588. [Google Scholar] [CrossRef]
  72. Chen, R.; Gao, B.; Liu, X.; Ruan, F.; Zhang, Y.; Lou, J.; Feng, K.; Wunsch, C.; Li, S.-M.; Dai, J.; et al. Molecular insights into the enzyme promiscuity of an aromatic prenyltransferase. Nat. Chem. Biol. 2017, 13, 226–234. [Google Scholar] [CrossRef]
  73. Ni, W.; Zheng, Z.; Liu, H.; Wang, P.; Wang, H.; Sun, X.; Yang, Q.; Fang, Z.; Tang, H.; Zhao, G. Combining mutagenesis on Glu281 of prenyltransferase NovQ and metabolic engineering strategies for the increased prenylated activity towards menadione. Appl. Microbiol. Biotechnol. 2020, 104, 4371–4382. [Google Scholar] [CrossRef] [PubMed]
  74. Yang, X.; Yang, J.; Jiang, Y.; Yang, H.; Yun, Z.; Rong, W.; Yang, B. Regiospecific synthesis of prenylated flavonoids by a prenyltransferase cloned from Fusarium oxysporum. Sci. Rep. 2016, 6, 24819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. He, B.-B.; Bu, X.-L.; Zhou, T.; Li, S.-M.; Xu, M.-J.; Xu, J. Combinatory biosynthesis of prenylated 4-hydroxybenzoate derivatives by overexpression of the substrate-promiscuous prenyltransferase XimB in engineered E. coli. ACS Synth. Biol. 2018, 7, 2094–2104. [Google Scholar] [CrossRef]
  76. Rainha, J.; Rodrigues, J.L.; Rodrigues, L.R. CRISPR-Cas9: A powerful tool to efficiently engineer Saccharomyces cerevisiae. Life 2021, 11, 13. [Google Scholar] [CrossRef] [PubMed]
  77. Gomes, D.; Rainha, J.; Rodrigues, L.R.; Rodrigues, J.L. Yeast synthetic biology approaches for the production of valuable polyphenolic compounds. In Synthetic Biology of Yeasts; Harzevili, F., Ed.; Springer International Publishing: Berlin/Heidelberg, Germany, 2022; pp. 119–156. [Google Scholar]
  78. Yuan, S.-F.; Yi, X.; Johnston, T.G.; Alper, H.S. De novo resveratrol production through modular engineering of an Escherichia coli–Saccharomyces cerevisiae co-culture. Microb. Cell Fact. 2020, 19, 143. [Google Scholar] [CrossRef]
  79. Häser, K.; Wenk, H.H.; Schwab, W. Biocatalytic production of dihydrocoumarin from coumarin by Saccharomyces cerevisiae. J. Agric. Food Chem. 2006, 54, 6236–6240. [Google Scholar] [CrossRef]
  80. Watts, K.T.; Lee, P.C.; Schmidt-Dannert, C. Biosynthesis of plant-specific stilbene polyketides in metabolically engineered Escherichia coli. BMC Biotechnol. 2006, 6, 22. [Google Scholar] [CrossRef] [Green Version]
  81. Bannai, H.; Tamada, Y.; Maruyama, O.; Nakai, K.; Miyano, S. Extensive feature detection of N-terminal protein sorting signals. Bioinformatics 2002, 18, 298–305. [Google Scholar] [CrossRef] [Green Version]
  82. Almagro Armenteros, J.J.; Tsirigos, K.D.; Sønderby, C.K.; Petersen, T.N.; Winther, O.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat. Biotechnol. 2019, 37, 420–423. [Google Scholar] [CrossRef] [Green Version]
  83. Armenteros, J.J.A.; Salvatore, M.; Emanuelsson, O.; Winther, O.; von Heijne, G.; Elofsson, A.; Nielsen, H. Detecting sequence signals in targeting peptides using deep learning. Life Sci. Alliance 2019, 2, e201900429. [Google Scholar] [CrossRef]
  84. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Goldberg, T.; Hecht, M.; Hamp, T.; Karl, T.; Yachdav, G.; Ahmed, N.; Altermann, U.; Angerer, P.; Ansorge, S.; Balasz, K.; et al. LocTree3 prediction of localization. Nucleic Acids Res. 2014, 42, W350–W355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Pierleoni, A.; Martelli, P.L.; Fariselli, P.; Casadio, R. BaCelLo: A balanced subcellular localization predictor. Bioinformatics 2006, 22, e408–e416. [Google Scholar] [CrossRef] [PubMed]
Molecules 27 07230 g001 550
Figure 1. Coumarin’s biosynthetic pathway in Escherichia coli (a). Linear furanocoumarin’s pathway in plants (b). 4CL1: 4-coumarate-CoA ligase 1; C2′H: p-coumaroyl-CoA 2′-hydroxylase; C3H: 4-coumarate 3-hydroxylase; CCoAOMT: caffeoyl-CoA 3-O-methyltransferase; DMAPP: dimethylallylpyrophosphate; F6′H: feruloyl-CoA 6′-hydroxylase; MS: marmesin synthase; NADPH: nicotinamide adenine dinucleotide phosphate; PS: psoralen synthase; PT: prenyltransferase; TAL: tyrosine ammonia lyase.
Figure 1. Coumarin’s biosynthetic pathway in Escherichia coli (a). Linear furanocoumarin’s pathway in plants (b). 4CL1: 4-coumarate-CoA ligase 1; C2′H: p-coumaroyl-CoA 2′-hydroxylase; C3H: 4-coumarate 3-hydroxylase; CCoAOMT: caffeoyl-CoA 3-O-methyltransferase; DMAPP: dimethylallylpyrophosphate; F6′H: feruloyl-CoA 6′-hydroxylase; MS: marmesin synthase; NADPH: nicotinamide adenine dinucleotide phosphate; PS: psoralen synthase; PT: prenyltransferase; TAL: tyrosine ammonia lyase.
Molecules 27 07230 g001
Molecules 27 07230 g002 550
Figure 2. Umbelliferone production from p-coumaric acid using p-coumaroyl-CoA 2′-hydroxylase (C2′H) expressed with different tags (His-tag, S-tag, and GST-tag) (a) and using different culture media (LB + LB, LB, LB + M9, and M9) (b). p-Coumaric acid was added at a concentration of 0.7 mM (a) and 1 mM (b).
Figure 2. Umbelliferone production from p-coumaric acid using p-coumaroyl-CoA 2′-hydroxylase (C2′H) expressed with different tags (His-tag, S-tag, and GST-tag) (a) and using different culture media (LB + LB, LB, LB + M9, and M9) (b). p-Coumaric acid was added at a concentration of 0.7 mM (a) and 1 mM (b).
Molecules 27 07230 g002
Molecules 27 07230 g003 550
Figure 3. Coumarin production in Escherichia coli using different hydroxycinnamic acids as precursors. Umbelliferone, esculetin, and scopoletin are catalyzed from p-coumaric acid, caffeic acid and ferulic acid when pAC-4CL1 and pET28GST-C2′H are expressed. 4CL1: 4-coumarate-CoA ligase 1; C2′H: p-coumaroyl-CoA 2′-hydroxylase.
Figure 3. Coumarin production in Escherichia coli using different hydroxycinnamic acids as precursors. Umbelliferone, esculetin, and scopoletin are catalyzed from p-coumaric acid, caffeic acid and ferulic acid when pAC-4CL1 and pET28GST-C2′H are expressed. 4CL1: 4-coumarate-CoA ligase 1; C2′H: p-coumaroyl-CoA 2′-hydroxylase.
Molecules 27 07230 g003
Molecules 27 07230 g004 550
Figure 4. Coumarin production in Escherichia coli using tyrosine as precursor and different combinations of plasmids. Coumarins concentrations (a), hydroxycinammic acids concentrations (b). 4CL1: 4-coumarate-CoA ligase 1; C2′H: p-coumaroyl-CoA 2′-hydroxylase; C3H: 4-coumarate 3-hydroxylase; CCoAOMT: caffeoyl-CoA 3-O-methyltransferase; TAL: tyrosine ammonia lyase.
Figure 4. Coumarin production in Escherichia coli using tyrosine as precursor and different combinations of plasmids. Coumarins concentrations (a), hydroxycinammic acids concentrations (b). 4CL1: 4-coumarate-CoA ligase 1; C2′H: p-coumaroyl-CoA 2′-hydroxylase; C3H: 4-coumarate 3-hydroxylase; CCoAOMT: caffeoyl-CoA 3-O-methyltransferase; TAL: tyrosine ammonia lyase.
Molecules 27 07230 g004
Table
Table 1. Plasmids used in this study.
Table 1. Plasmids used in this study.
PlasmidsConstructSource
pRSFDuet-1RSF1030 ori, lacI, double PT7lac, KanRNovagen
pCDFDuet-1CloDF13 ori, lacI, double PT7lac, SpecRNovagen
pNZY29-C2′HColE1(pBR322) ori, AmpR; pNZY29 carrying codon-optimized p-coumaroyl-CoA 2′-hydroxylase (C2′H) from Ipomoea batatasNZYTech
pHTP0-PsPT1ColE1(pBR322) ori, AmpR; pHTP0 carrying codon-optimized prenyltransferase 1 from Pastinaca sativa (PsPT1)NZYTech
pHTP0-PSColE1(pBR322) ori, AmpR; pHTP0 carrying codon-optimized psoralen synthase (PS) from Ammi majusNZYTech
pET28GST-LICColE1(pBR322) ori, lacI, PT7lac, KanR; carrying levansucrase (SacB) from Bacillus subtilis and glutathione S-transferase (GST)Addgene (#26101)
pAC-4CL1p15A ori, Plac, CmR, pACYC184-derived plasmid carrying 4-coumarate-CoA ligase 1 (4CL1) from Arabidopsis thalianaAddgene (#35947)
pETDuet-PcPT-XimD-XimEpBR322 ori, lacI, double PT7lac, AmpR; pETDuet-1 carrying codon-optimized prenyltransferase from Petroselinum crispum (PcPT) and marmesin synthase (XimD) and snoaL-like cyclase (XimE) from Streptomyces xiamenensis Addgene (#172654)
pMVApCDFDuet-1 carrying hydroxymethylglutaryl-CoA synthase (MvaS) and bifunctional acetoacetyl-CoA thiolase and 3-hydroxy-3-methylglutaryl-CoA reductase (MvaE) from Enterococcus faecalis, isopentenyl pyrophosphate—dimethylallylpyrophosphate isomerase (IDI) from E. coli BL21 and mevalonate (MVA) kinase (ERG12), phospho-MVA kinase (ERG8) and MVA pyrophosphate decarboxylase (ERG19) from Saccharomyces cerevisiaeAddgene (#121149)
pETDuet-CPR2pBR322 ori, lacI, double PT7lac, AmpR; pETDuet-1 carrying cytochrome P450 reductase 2 (CPR2) from A. thaliana[59]
pYeDP60-CYP76F112pUC ori, 2 µ ori, URA3, TetR, AmpR; pYeDP60 carrying marmesin synthase from Ficus carica (CYP76F112)[40]
pETDuet-TALpBR322 ori, lacI, double PT7lac, AmpR; pETDuet-1 carrying tyrosine ammonia lyase (TAL) from Rhodotorula glutinis[45]
pETDuet-TAL-C3HpETDuet-TAL carrying 4-coumarate 3-hydroxylase (C3H) from Saccharothrix espanaensis[45]
pCDFDuet-CCoAOMTpCDFDuet-1 carrying caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) from Medicago sativa[43]
pET28GST-C2′HpET28GST-LIC without SacB and carrying C2′H in frame with GST and His-tagThis study
pRSFDuet-C2′HpRSFduet-1 carrying C2′H without His-tag or S-tag in frameThis study
pRSFDuet-His-tag-C2′HpRSFDuet-1 carrying C2′H with His-tag in frameThis study
pRSFDuet-C2′H-S-tagpRSFDuet-1 carrying C2′H with S-tag in frameThis study
pETDuet-TAL-PcPTpETDuet-TAL carrying PcPTThis study
pCDFDuet-PcPTpCDFDuet-1 carrying PcPTThis study
pRSFDuet-PcPTpRSFDuet-1 carrying PcPTThis study
pET28GST-PcPTpET28GST-LIC without SacB and carrying PcPTThis study
pRSFDuet-∆48PcPTpRSFDuet-1 carrying truncated PcPTThis study
pET28-∆48PcPTpET28GST-LIC without SacB and carrying truncated PcPTThis study
pRSFDuet-PsPT1pRSFDuet-1 carrying PsPT1This study
pET28GST-PsPT1pET28GST-LIC without SacB and carrying PsPT1This study
pRSFDuet-∆48PsPT1pRSFDuet-1 carrying truncated PsPT1This study
pET28-∆48PsPT1pET28GST-LIC without SacB and carrying truncated PsPT1This study
pRSFDuet-PSpRSFDuet-1 carrying PSThis study
pET28GST-PSpET28GST-LIC without SacB and carrying PSThis study
pRSFDuet-MA-∆37PSpRSFDuet-1 carrying MA-∆37PSThis study
pET28-MA-∆37PSpET28GST-LIC without SacB and GST and carrying MA-∆37PSThis study
pRSFDuet-8RP-∆37PSpRSFDuet-1 carrying 8RP-∆37PSThis study
pET28-8RP-∆37PSpET28GST-LIC without SacB and GST and carrying 8RP-∆37PSThis study
pRSFDuet-2C3-∆37PSpRSFDuet-1 carrying 2C3-∆37PSThis study
pET28-2C3-∆37PSpET28GST-LIC without SacB and GST and carrying 2C3-∆37PSThis study
pRSFDuet-28tag-∆37PSpRSFDuet-1 carrying 28tag-∆37PSThis study
pET28-28tag-∆37PSpET28GST-LIC without SacB and GST and carrying 28tag-∆37PSThis study
pCDFDuet-GroESLpCDFDuet-1 carrying chaperones GroES and GroEL from E. coli BL21 in an operonThis study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodrigues, J.L.; Gomes, D.; Rodrigues, L.R. Challenges in the Heterologous Production of Furanocoumarins in Escherichia coli. Molecules 2022, 27, 7230. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27217230

AMA Style

Rodrigues JL, Gomes D, Rodrigues LR. Challenges in the Heterologous Production of Furanocoumarins in Escherichia coli. Molecules. 2022; 27(21):7230. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27217230

Chicago/Turabian Style

Rodrigues, Joana L., Daniela Gomes, and Lígia R. Rodrigues. 2022. "Challenges in the Heterologous Production of Furanocoumarins in Escherichia coli" Molecules 27, no. 21: 7230. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules27217230

Article Metrics

Back to TopTop