Triterpenoids are high-value plant metabolites with numerous applications in medicine, agriculture, food, and home and personal care products. However, plants produce triterpenoids in low abundance, and their complex structures make their chemical synthesis prohibitively expensive and often impossible. As such, the yeast Saccharomyces cerevisiae has been explored as an alternative means of production. An important triterpenoid is oleanolic acid because it is the precursor to many bioactive triterpenoids of commercial interest, such as QS-21 which is being evaluated as a vaccine adjuvant in clinical trials against HIV and malaria. Oleanolic acid is derived from 2,3-oxidosqualene (natively produced by yeast) via a cyclisation and a multi-step oxidation reaction, catalysed by a Î²-amyrin synthase and a cytochrome P450 of the CYP716A subfamily, respectively. Although many homologues have been characterised, previous studies have used arbitrarily chosen Î²-amyrin synthases and CYP716As to produce oleanolic acid and its derivatives in yeast. This study presents the first comprehensive comparison of Î²-amyrin synthase and CYP716A enzyme activities in yeast. Strains expressing different homologues are compared for production, revealing 6.3- and 4.5-fold differences in Î²-amyrin and oleanolic acid productivities and varying CYP716A product profiles, which are important to consider when engineering strains for the production of bioactive oleanolic acid derivatives.
Citation: Dale MP, Moses T, Johnston EJ, Rosser SJ (2020) A systematic comparison of triterpenoid biosynthetic enzymes for the production of oleanolic acid in Saccharomyces cerevisiae. PLoS ONE 15(5):
Editor: Alvaro Galli, CNR, ITALY
Received: November 7, 2019; Accepted: April 4, 2020; Published: May 1, 2020
Copyright: Â© 2020 Dale et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was funded by grants awarded by the Industrial Biotechnology Innovation Centre (Project-2016-150; www.ibioic.com; author TM) and the Biotechnology and Biological Sciences Research Council (BBSRC; IBCatalyst: Project No. 102297; author EJJ). MD was funded by the East of Scotland Bioscience Doctoral Training Partnership (EASTBIO DTP) with BBSRC for a doctoral fellowship. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Oleanolic acid is a plant specialized metabolite of considerable commercial interest that belongs to the oleanane class of triterpenoids. The oleanane scaffold is composed of 30 carbons arranged in a pentacyclic structure that is decorated with various functional groups to obtain vast structural and functional diversity. Several plant species produce oleanolic acid , which can accumulate as itself or as oxidized and glycosylated derivatives called saponins. Oleanolic acid is a particularly important triterpenoid, being the precursor to a range of bioactive triterpenoids such as the vaccine adjuvant QS-21, which is used in a vaccine against feline leukaemia virus and is being investigated in human vaccines against HIV and malaria . Semi-synthetic derivatives of oleanolic acid have been investigated in the treatment of various diseases, including leukaemia, chronic kidney disease, type II diabetes, osteoporosis and Friedreichâ€™s ataxia [3,4], while glycosylated derivatives of oleanolic acid have been shown to act as feeding deterrents to common crop pests .
As such, there is considerable interest in the production of oleanolic acid derivatives in industrially relevant quantities. Plants typically produce triterpenoids in low amounts in response to stress, and the complex stereochemistry of triterpenoids makes their chemical synthesis economically prohibitive and often impossible. Therefore, use of the yeast Saccharomyces cerevisiae, which is genetically tractable and widely used in industrial fermentations, has been explored for triterpenoid production [6â€“9]. Oleanolic acid biosynthesis occurs through the cyclisation of 2,3-oxidosqualene (natively produced by yeast) to produce Î²-amyrin, followed by the oxidation of carbon-28 (C-28) methyl group of Î²-amyrin to a carboxyl group. This oxidation reaction is predominantly catalysed by cytochrome P450 monooxygenases (P450s) of the CYP716A subfamily in conjunction with a cytochrome P450 reductase (CPR) partner, and proceeds through three steps, producing the intermediates erythrodiol and oleanolic aldehyde, which have a hydroxyl and aldehyde group at C-28, respectively (Fig 1A). Oleanolic acid can thus be produced in yeast by the co-expression of a Î²-amyrin synthase (BAS), a CYP716A and a CPR, with the expression of additional P450s and UDP-dependent glycosyltransferases allowing the production of bioactive oleanolic acid derivatives.
(a) Oleanolic acid biosynthesis pathway from 2,3-oxidosqualene. (b-j) GC/QTOF-MS confirms Î²-amyrin production by GhBAS, and oleanolic acid production by CYP716A147 and CYP716A48. (b, d) Total ion chromatograms of strains expressing GhBAS, CYP716A147 and CYP716A48, with control strains expressing GgBAS or empty vector, and a mixture of ergosterol, Î²-amyrin, erythrodiol and oleanolic acid standards at 50 Î¼g/mL. (c) Mass spectra of the GhBAS Î²-amyrin peak matched to the standard. (e-h) Mass spectra of erythrodiol and oleanolic acid peaks from the CYP716A strains matched to standards. (i-j) Mass spectra of predicted oleanolic aldehyde peaks matched those previously reported.
The choice of enzyme homologues can significantly impact productivity, as exemplified by Î²-amyrin titres of 19 and 36 mg/L when BASs from Medicago truncatula (MtBAS) and Glycyrrhiza glabra (GgBAS) were expressed in yeast . In addition, titres of 107 and 1.9 mg/L Î²-amyrin have been reported for GgBAS and PgBAS1 (from Panax ginseng) . Despite these broad-ranging yields, no systematic comparison of BAS and CYP716A productivities in yeast have been reported. The advantage of such comparisons has previously been demonstrated for the CYP93E subfamily of P450s, which catalyse the three-step oxidation of the Î²-amyrin C-24 methyl group to a carboxyl group. A 52-fold difference in accumulation of 24-hydroxy-Î²-amyrin was reported between strains expressing eight different CYP93E variants . As observed for the CYP93E subfamily, a comparative analysis of CYP716A enzymes would be particularly relevant due to the multi-step nature of the reaction catalysed, because the ratios of accumulated intermediates and product could differ between homologues. We therefore decided to conduct a comprehensive comparison of the productivities of BAS and CYP716A homologues in yeast, to identify the most productive variants for oleanolic acid production. Such variants could then be used to engineer high yields of oleanolic acid and its derivatives from yeast.
2. Materials and methods
Coprostan-3-ol, ethyl acetate and pyridine were purchased from Sigma-Aldrich. n-Hexane was obtained from Merck, N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) from Thermo Scientific, and potassium hydroxide pellets from Scientific Laboratory Supplies. Î²-Amyrin, erythrodiol and oleanolic acid standards were purchased from Extrasynthese.
The synthetic defined medium was composed of 6.9 g/L yeast nitrogen base without amino acids (Sigma) and 0.79 g/L complete supplement mixture or 0.77 g/L single drop-out of uracil (Formedium) with carbon source (either 2% á´…-glucose [Fisher Chemical], or 2% á´…-(+)-galactose [Acros Organics]).
2.3 Plasmid construction
All genes except AsBAS were codon-optimised for S. cerevisiae, synthesised by GeneArt and assembled into low-copy number plasmids expressing the URA3 auxotrophic marker . AsBAS could not be synthesised and the native Avena strigosa sequence was a gift from the Osbourn lab (John Innes Centre, UK).
Plasmid backbones were initially constructed using parts provided in the MoClo Yeast Toolkit (MoClo-YTK) . Plasmids pCPS1ULA, pCP12ULA and pCP2EULA each consisted of a CEN6/ARS4 low-copy origin, URA3 selection marker, ColE1 origin and Amicillin resistance gene for maintanence in E. coli, as well as a cloning site for either a BAS, CYP716A or ATR2 transcriptional unit (consisting of a promoter, ORF and terminator), respectively. Plasmid pMGSEULS instead had an E. coli Spectinomycin resistance gene and a site into which three transcriptional units (BAS, CYP716A and ATR2) can be cloned to form an oleanolic acid biosynthetic pathway.
BAS, CYP716A and ATR2 transcriptional units, and oleanolic pathways, were subsequently assembled into these vectors using the MoClo-YTK system . BAS genes were cloned into pCPS1ULA with the PGAL1 promoter and TDH1 terminator, CYP716A genes into pCP12ULA with PGAL10 and the ENO2 terminator, and ATR2 into pCP2EULA with PGAL3 and the PGK1 terminator. The resulting transcriptional units were then cloned into pMGSEULS to form oleanolic acid pathways, in which BAS, CYP716A and ATR2 are co-expressed.
2.4 Generation and culturing of strains
Yeast strains were generated by transforming plasmids into S. cerevisiae BY4741 (Î²-amyrin production strains by transforming pCPS1ULA-BAS plasmids, and oleanolic acid strains by transforming pMGSEULS-OA plasmids). Empty pCPS1ULA vector was used as a control. All strains were grown at 30Â°C in synthetic defined medium lacking uracil and containing either 2% á´…-glucose (for selection of transformants, maintenance of strains, or precultures prior to inoculation in culturing medium) or 2% á´…-(+)-galactose (for quantifying triterpenoid production). For production of triterpenoids, strains were precultured for 24 h, inoculated to a starting OD600 of 0.2 in galactose-containing medium, and grown for 96 h in glass vials.
2.5 Triterpenoid extraction
Cells from 1 mL culture were separated from media, resuspended in 500 Î¼L lysis buffer (40% potassium hydroxide, 50% ethanol, 80 nM coprostan-3-ol) and boiled for 10 min. Lysates were cooled to ambient temperature and extracted twice with 500 Î¼L n-hexane, and then twice with 500 Î¼L ethyl acetate. Organic phases were pooled, 600 Î¼L was evaporated to dryness in a GeneVac EZ-2 Elite evaporator, and trimethylsilylated with 50 Î¼L pyridine: N-methyl-N-trimethylsilyltrifluoroacetamide (1:4) for gas chromatography quadrupole time-of-flight mass spectrometry (GC/QTOF-MS) analysis.
2.6 Gas chromatographyâ€“Mass spectrometry analysis
The triterpene profiles were analysed on an Agilent 7890B gas chromatogram (GC) coupled to an Agilent 7200B quadrupole time-of-flight mass spectrometer (QTOF-MS) with GERSTEL multipurpose sampler (MPS) robotics (Anatune). Trimethysilylated samples (1 Î¼L) were injected at a split ratio of 10:1, with a split flow of 10 mL/min into a DB-5ms 40 m Ã— 250 Î¼m Ã— 0.25 Î¼m GC column (Agilent Technologies). Helium was used as the carrier gas at a flow rate of 1 mL/min, the inlet was set to 250Â°C and the GC oven was programmed to 100Â°C for 2 min, followed by ramping at 25Â°C/min to 300Â°C, where it was held for 25 min. The ion source were set to 230Â°C, 35 Î¼A filament current, 70 eV electron energy and the mass range of 60â€“900 m/z was scanned at an acquisition rate of 4 spectra/s with a solvent delay of 13 min. Total ion chromatograms and mass spectra were analysed using the Agilent MassHunter Qualitative Analysis B.07.00 software, and peak areas were calculated using the Agile 2 integrator method. GC/Q-TOF MS peak areas were used to obtain standard curves for Î²-amyrin, erythrodiol and oleanolic acid. The concentration range of 0 to 120 Î¼g/mL, with increments of 10 Î¼g/mL, was used and the linear region of the curves allowed quantifying concentration of the triterpenes.
3.1 Characterization of new BAS and CYP716A enzymes
We first identified new BAS and CYP716A homologues. BLASTp searches on NCBI using GgBAS from G. glabra  and CYP716A12 from M. truncatula  as queries identified hits with > 70% amino acid sequence identity. GhBAS from Gossypium hirsutum (cotton), CYP716A48 from Olea europaea (olive) and CYP716A147 from Theobroma cacao (cocoa tree) were selected for further analysis (S1 Table). While G. hirsutum and T. cacao produce triterpenoids [14,15], genes encoding oleanane biosynthetic enzymes have not been previously characterised. Similarly, O. europaea accumulates oleanolic acid , but a P450 involved in its biosynthesis has not yet been described.
In silico analysis of the GhBAS sequence (S1 Fig) identified the highly conserved DCTAE motif involved in initiation of the cyclization reaction [16,17], the MWCYCR motif known to confer product specificity , and the QW motifs that are speculated to be involved in BAS protein stability [19â€“21]. Sequence alignment of CYP716A147 and CYP716A48 with functionally characterised CYP716A proteins (S2 Fig) identified the proton transfer groove, ExxR motif, PERF motif and haem-binding loop, along with the N-terminal transmembrane anchor.
The selected genes were functionally characterised in yeast. GhBAS was expressed alone, while CYP716A48 and CYP716A147 were co-expressed with GgBAS and ATR2, a CPR from Arabidopsis thaliana. While yeast encodes its own native CPR (NCP1), it is often unable to support the catalytic activity of heterologously expressed plant P450s . As such, a plant CPR gene is usually co-expressed with plant P450s, typically one of the A. thaliana CPRs, ATR1 or ATR2 [7,8,23], of which ATR2 has been reported to have a higher catalytic activity [24,25].
Compared to a strain carrying an empty vector (MD-N1), the strain expressing GhBAS produced Î²-amyrin, with an identical retention time and fragmentation pattern as an authentic standard (Fig 1B and 1C). Similarly, peaks with retention times and fragmentation patterns matching standards of erythrodiol and oleanolic acid were observed from both CYP716A-expressing strains (Fig 1Dâ€“1H). An additional peak not present in MD-N1 was observed in these strains (Fig 1D), with mass spectra matching that reported for an oleanolic aldehyde standard (Fig 1I and 1J) . Thus, GhBAS is a monofunctional BAS that cyclises 2,3-oxidosqualene to Î²-amyrin, while CYP716A147 and CYP716A48 are C-28 oxidases that oxidise Î²-amyrin to oleanolic acid via erythrodiol and oleanolic aldehyde intermediates.
3.2 Strains expressing AaBAS or CqBAS1 produce the highest titre of Î²-amyrin
We screened the amount of Î²-amyrin produced in yeast by 12 BAS proteins from different plant species (S1 Table). These included the BAS variants most commonly used in engineering studies (GgBAS, MtBAS and PgBAS1), as well as AsBAS, BvBAS, AaBAS, CqBAS1, EtBAS, LjBAS, PtBAS, SlBAS and the newly characterised GhBAS.
All strains except that expressing AsBAS displayed reduced growth compared to the control strain MD-N1 (Fig 2A). Furthermore, all strains except that expressing PgBAS produced quantifiable titres of Î²-amyrin (> 1.7 mg/L), with 3.7- and 6.3-fold differences in titre (mg/L) and productivity (mg/L/OD600), respectively, between the most and least productive strains in our study (Fig 2B and 2C). Overall, the highest titre of 10.8 Â± 1.0 mg/L Î²-amyrin was produced by the strains expressing AaBAS or CqBAS1. The PtBAS, LjBAS and EtBAS strains closely followed, with Î²-amyrin titres of 9.0 Â± 0.7, 8.2 Â± 1.0 and 8.0 Â± 0.2 mg/L, respectively. The SlBAS strain produced the lowest titre of 2.9 Â± 0.3 mg/L Î²-amyrin.
Comparison of BAS (a-c) and CYP716A (d-g) proteins for Î²-amyrin and oleanolic acid production. (a, d) OD600 at end of culture (grey). (b, e) Titres (mg/L) of Î²-amyrin (red), erythrodiol (yellow) and oleanolic acid (cyan). (f) Oleanolic aldehyde peak areas obtained from GC/QTOF-MS analysis. (c, g) Î²-Amyrin (dark red) and oleanolic acid (dark cyan) productivities (mg/L/OD600). None, empty vector strain MD-N1 in a-c and strain expressing AaBAS only in d-g. Mean Â± SD for n = 3 is presented.
Interestingly, in addition to Î²-amyrin, a second peak was observed in the total ion chromatogram of the BvBAS strain (S3A Fig). This peak appeared to be a combination of two peaks, suggesting the presence of two additional compounds, and was unique to this strain, not appearing in the other BAS strains or the MD-N1 control strain. Comparison of its mass spectrum with that of an authentic lupeol standard indicated that one of the additional compounds contained within the peak most likely corresponds to this triterpene (S3B Fig).
3.3 Strain expressing CYP716AL1 produces the highest titre of oleanolic acid
We next compared the activities of 16 CYP716A family proteins from a range of plant species (S1 Table) for their ability to oxidise Î²-amyrin to oleanolic acid. These included CYP716A12, which is commonly used in engineering studies, the newly identified CYP716A147 and CYP716A48, as well as CYP716A15, CYP716AL1, CYP716A52v2, CYP716A44, CYP716A75, CYP716A79, CYP716A80, CYP716A83, CYP716A110, CYP716A140, CYP716A179, CYP716A244 and CYP716A1.
All strains showed comparable growth to the control strain expressing AaBAS alone (Fig 2D) but displayed marked differences in triterpene production (Fig 2E), with oleanolic acid titre (mg/L) and productivity (mg/L/OD600) differing 4.8- and 4.5-fold, respectively, between strains (Fig 2E). Oxidation of Î²-amyrin was observed in all strains except that expressing CYP716A140, with oleanolic acid produced in 11 of the 16 strains (cyan bars, Fig 2E). The highest and lowest titres of 14.3 Â± 1.6 and 3.0 Â± 0.0 mg/L oleanolic acid were produced by strains expressing CYP716AL1 and CYP716A52v2, respectively. CYP716AL1 also showed the highest oleanolic acid productivity (2.6 mg/L/OD600) and CYP716A52v2 the least (0.6 mg/L/OD600) (Fig 2G). Notably, all strains (except that expressing CYP716A140) accumulated varying amounts of erythrodiol (yellow bars, Fig 2E) and oleanolic aldehyde (green bars, Fig 2F; standard unavailable for quantification). Interestingly, CYP716A140 has previously been shown to oxidise Î²-amyrin to oleanolic acid, in a study using the M. truncatula CPR MTR1 and a non-codon optimised CYP716A140 sequence . These variances in the expression assay may explain the differences in activity observed for CYP716A140 in our study. The CYP716A75 strain produced the highest erythrodiol titre (24.6 Â± 3.9 mg/L), with trace amounts of oleanolic aldehyde and no oleanolic acid. Strains expressing CYP716A79, CYP716A110 and CYP716A1 accumulated 1.2 Â± 0.1, 1.2 Â± 0.1 and 4.5 Â± 0.7 mg/L erythrodiol, respectively, with no oleanolic aldehyde or acid.
3.4 Combining optimal enzymes to boost oleanolic acid production
To assess the impact of enzyme variant selection on oleanolic acid titres in yeast, we generated strains expressing SlBAS (the lowest Î²-amyrin titre) with CYP716AL1 (the highest oleanolic acid titre) or CYP716A52v2 (the lowest oleanolic acid titre) genes. Growth and oleanolic acid production was compared to strains expressing AaBAS (the highest Î²-amyrin titre) with CYP716AL1 or CYP716A52v2. All strains showed comparable growth after 96 h of culturing (Fig 3A). The expression of SlBAS with CYP716AL1 or CYP716A52v2 resulted in 3.9 Â± 0.2 and 2.8 Â± 0.0 mg/L oleanolic acid, respectively (Fig 3B). The combination of AaBAS with CYP716AL1 showed the highest titre (8.5 Â± 0.2 mg/L) and productivity (1.40 Â± 0.0 mg/L/OD600) of oleanolic acid (Fig 3B and 3C).
In this study, a comprehensive comparison of 12 BAS and 16 CYP716A variants for the production of oleanolic acid in yeast was performed. These included enzymes commonly used in engineering studies (GgBAS, MtBAS, PgBAS1 and CYP716A12), as well as one BAS and two CYP716As newly characterised in this study.
These enzymes displayed a range of activities and product specificities, and impacted growth rates to different extents, when expressed in yeast. Notably, AaBAS and CqBAS1 produced higher titres of Î²-amyrin (10.8 mg/L), and CYP716AL1 produced higher titres of oleanolic acid (14.3 mg/L), than the BAS and CYP716A enzymes commonly used in engineering studies (GgBAS, MtBAS, PgBAS1 and CYP716A12). In addition, CYP716A75 predominantly produced the triterpenoid erythrodiol, suggesting an inherent ability of CYP716A75 to mainly catalyse the hydroxylation of Î²-amyrin to erythrodiol, as previously reported . By contrast, CYP716AL1 and CYP716A12 produced substantially less erythrodiol, favouring the production of oleanolic acid. The use of P450s with different product specificities, such as those characterised here, will enable the production of a range of triterpenoids with specific modifications to the triterpene scaffold.
Interestingly, the strain expressing AsBAS was the only strain to show similar growth to the control strain MD-N1 carrying an empty vector. All other BAS strains displayed reduced growth compared to MD-N1. While experiments were not undertaken to determine the cause of this, it is notable that AsBAS was the only BAS gene that wasnâ€™t codon-optimised for expression in yeast, which could have led to lower production of the AsBAS enzyme compared to the other BASs, possibly resulting in reduced burden on the cells. Protein sequence alignment also showed that AsBAS differed significantly in sequence to the other BAS proteins, having 47â€“50% sequence identity to them; by contrast, the other BAS proteins had 69â€“93% sequence identity to one another.
While the expression of BAS genes reduced growth by differing amounts, the additional expression of CYP716A and ATR2 genes appeared to have much less impact on cell growth, with the growth of oleanolic acid-producing strains varying far less than the BAS strains. The observed reductions in growth may result from (i) amino acid burden caused by expression of the recombinant proteins, (ii) aberrant activities of the transgenic enzymes, (iii) accumulation of triterpenoid product and/or (iv) depletion of native metabolites (e.g. redirection of metabolic flux from sterol to triterpenoid biosynthesis) [3,29]. Most likely, a combination of these factors and others could explain the differences in growth observed between the strains.
The strains used in this study were derived from a BY4741 parent strain transformed with plasmids expressing the triterpenoid biosynthetic gene(s). These therefore represent unoptimised strains and numerous strain engineering approaches could be taken to improve titres. In an early study, overexpression of a truncated version of the native tHMG1 (a rate-limiting enzyme of the mevalonate pathway required for the production of 2,3-oxidosqualene), and downregulation of the native ERG7 (lanosterol synthase which uses 2,3-oxidosqualene as substrate), increased Î²-amyrin titres in yeast expressing AaBAS by 1.5-fold (to 6 mg/L), compared with expression of AaBAS alone . Higher titres of 50 mg/L Î²-amyrin were later achieved by expressing GgBAS in a strain overexpressing tHMG1 and with repression of the ERG7 gene . The strain was cultured in the presence of methyl-Î²-cyclodextrin, which improved titres, possibly by sequestering Î²-amyrin in the culture media and preventing toxic build-up within the cells. Meanwhile, another study reported oleanolic acid titres of 71 mg/L . Two copies of GgBAS were co-expressed with CYP716A12 and AtATR1, together with overexpression of tHMG1 and the native genes ERG9 and ERG1, which produce squalene from farnesyl pyrophosphate (FPP, a product of the mevalonate pathway) and convert it to 2,3-oxidosqualene, respectively. In addition, the strains were grown in rich medium (YPD), which may have facilitated growth and the production of triterpenoids. An organelle engineering approach increased titres of Î²-amyrin and its derivative 28-O-glucose-medicagenic acid by 8- and 16-fold compared to an unengineered strain . Knock-out of the PAH1 gene, encoding phosphatidic acid phosphatase, resulted in proliferation of the ER and a corresponding increase in the intracellular levels of triterpenoid biosynthetic enzymes, boosting strain productivity. The use of fermenter apparatus has also resulted in increased triterpenoid production. For example, oleanolic acid titres of 607 mg/L were achieved in a fermenter, compared with 186 mg/L in shake flasks using the same strain .
This study compares BAS and CYP716A enzmye productivities in yeast, reports a range of product specificities between CYP716A enzymes, and identifies that expression of AsBAS has minimal impact on yeast growth while producing relatively high titres of Î²-amyrin. This work should therefore serve as an important resource for the production of industrially relevant triterpenoids in yeast.
S1 Fig. Multiple sequence alignment of GhBAS with functionally characterized BASs.
The QW (orange), MWCYCR (blue) and DCTAE (red) motifs are highlighted. The consensus is shown under the alignment as identical residues (*), conserved substitutions (.) and semi-conserved substitutions (:). Protein sequences were aligned using Clustal Omega (www.ebi.ac.uk/Tools/msa/clustalo/) using default parameters. The BAS protein sequences were retrieved from GenBank: AaBAS, ACA13386.1; AsBAS, CAC84558.1; BvBAS, AFF27505.1; CqBAS1, ANY30852.1; EtBAS, BAE43642.1; GgBAS, BAA89815.1; GhBAS, XP_016749748.1; LjBAS, BAE53429.1; MtBAS, CAD23247.1, PgBAS, BAA33461.1, PtBAS, ABL07607.1; SlBAS, ADU52574.1; TcBAS, XP_017979420.
S2 Fig. Multiple sequence alignment of CYP716A proteins.
Protein sequences were aligned using Clustal Omega (www.ebi.ac.uk/Tools/msa/clustalo/) with default parameters. Conserved sequences are highlighted as transmembrane domain (green) predicted by the TMHMM server (http://www.cbs.dtu.dk/services/TMHMM); proline cluster that forms a hinge between the transmembrane and cytoplasmic domains (P450 consensus PPxP, purple); A/G-G-X-D/E-T-T/S motif (proton transfer groove, light blue); ExxR motif (blue); PERF motif (red); FxxGxRxCxG motif (heme binding loop, orange). The sequences contained an additional â€˜GSâ€™ at the C-terminus, but were otherwise identical to the GenBank sequences. The protein sequences were retrieved from GenBank: CYP716A1, NP_198460.1; CYP716A12, CBN88269.1; CYP716A15, BAJ84106; CYP716A44, AK329870.1; CYP716A48, BAP59949.1; CYP716A52v2, AFO63032; CYP716A75, AHF22088.1; CYP716A79, ANY30854.1; CYP716A80, ALR73782.1; CYP716A83, AOG74832.1; CYP716A110, AOG74847.1; CYP716A140, AOG74836.1; CYP716A147, XP_007023618; CYP716A179, BAW34647.1; CYP716A244, APZ88353.1; CYP716AL1, AEX07773.
S3 Fig. BvBAS is a multifunctional OSC producing multiple triterpenes.
(A) Total ion chromatogram (TIC) of strains expressing AaBAS, BvBAS and the MD-N1 control strain (expressing no BAS). While both AaBAS and BvBAS have a peak corresponding to Î²-amyrin, the BvBAS strain has a second unique peak (*) that is potentially a combination of two peaks and which is not present in either the AaBAS or control strains. (B) Mass spectrum at retention time marked by asterix (*) (top). Comparison with a mass spectrum of an authentic lupeol standard (bottom) suggests one of the compounds contained within the peak corresponds to lupeol.
Pollier J, Goossens A. Oleanolic acid. Phytochemistry. 2012;77:10â€“5. pmid:22377690
Sun H, Xie Y, Ye Y. Advances in saponin-based adjuvants. Vaccine. 2009;27(12):1787â€“96. pmid:19208455
Moses T, Pollier J, Thevelein JM, Goossens A. Bioengineering of plant (tri)terpenoids: from metabolic engineering of plants to synthetic biology in vivo and in vitro. New Phytol. 2013;200(1):27â€“43. pmid:23668256
Strawser C, Schadt K, Hauser L, McCormick A, Wells M, Larkindale J, et al. Pharmacological therapeutics in Friedreich ataxia: the present state. Expert Rev Neurother. 2017;17(9):895â€“907. pmid:28724340
Augustin JM, Drok S, Shinoda T, Sanmiya K, Nielsen JK, Khakimov B, et al. UDP-Glycosyltransferases from the UGT73C Subfamily in Barbarea vulgaris Catalyze Sapogenin 3-O-Glucosylation in Saponin-Mediated Insect Resistance. Plant Physiol. 2012;160(4):1881â€“95. pmid:23027665
Moses T, Pollier J, Almagro L, Buyst D, Van Montagu M, PedreÃ±o MA, et al. Combinatorial biosynthesis of sapogenins and saponins in Saccharomyces cerevisiae using a C-16Î± hydroxylase from Bupleurum falcatum. Proc Natl Acad Sci U S A. 2014;111(4):1634â€“9. pmid:24434554
Dai Z, Wang B, Liu Y, Shi M, Wang D, Zhang X, et al. Producing aglycons of ginsenosides in bakersâ€™ yeast. Sci Rep. 2014;4:3698. pmid:24424342
Arendt P, Miettinen K, Pollier J, Rycke R De, Callewaert N, Goossens A. An endoplasmic reticulum-engineered yeast platform for overproduction of triterpenoids. Metab Eng. 2017;40:165â€“75. pmid:28216107
Fukushima EO, Seki H, Sawai S, Suzuki M, Ohyama K, Saito K, et al. Combinatorial Biosynthesis of Legume Natural and Rare Triterpenoids in Engineered Yeast. Plant Cell Physiol. 2013;54(5):740â€“9. pmid:23378447
Moses T, Thevelein JM, Goossens A, Pollier J. Comparative analysis of CYP93E proteins for improved microbial synthesis of plant triterpenoids. Phytochemistry. 2014;108:47â€“56. pmid:25453910
Lee ME, DeLoache WC, Cervantes B, Dueber JE. A Highly Characterized Yeast Toolkit for Modular, Multipart Assembly. ACS Synth Biol. 2015;4(9):975â€“86. pmid:25871405
Hayashi H, Huang P, Kirakosyan A, Inoue K, Hiraoka N, Ikeshiro Y, et al. Cloning and Characterization of a cDNA Encoding Î²-Amyrin Synthase Involved in Glycyrrhizin and Soyasaponin Biosyntheses in Licorice. Biol Pharm Bull. 2001;24(8):912â€“6. pmid:11510484
Carelli M, Biazzi E, Panara F, Tava A, Scaramelli L, Porceddu A, et al. Medicago truncatula CYP716A12 Is a Multifunctional Oxidase Involved in the Biosynthesis of Hemolytic Saponins. Plant Cell. 2011;23:3070â€“81. pmid:21821776
Shakhidoyatov KM, Rashkes AM, Khidyrova NK. Components of cottonplant leaves, their functional role and biological activity. Chem Nat Compd. 1997;33(6):605â€“16.
Zainal B, Abdah M, Taufiq-Yap Y, Roslida A, Rosmin K. Anticancer Agents from Non-Edible Parts of Theobroma cacao. Nat Prod Chem Res. 2014;2(4).
Hoshino T. Î²-Amyrin biosynthesis: catalytic mechanism and substrate recognition. Org Biomol Chem. 2017;15:2869â€“91. pmid:28294269
Ito R, Masukawa Y, Hoshino T. Purification, kinetics, inhibitors and CD for recombinant Î²-amyrin synthase from Euphorbia tirucalli L and functional analysis of the DCTA motif, which is highly conserved among oxidosqualene cyclases. FEBS J. 2013;280:1267â€“80. pmid:23294602
Kushiro T, Shibuya M, Masuda K, Ebizuka Y. Mutational studies on triterpene synthases: Engineering lupeol synthase into Î²-amyrin synthase. J Am Chem Soc. 2000;122(29):6816â€“24.
Poralla K. The possible role of a repetitive amino acid motif in evolution of triterpenoid cyclases. Bioorg Med Chem Lett. 1994;4(2):285â€“90.
Salmon M, Thimmappa RB, Minto RE, Melton RE, Hughes RK, Oâ€™Maille PE, et al. A conserved amino acid residue critical for product and substrate specificity in plant triterpene synthases. Proc Natl Acad Sci. 2016;113(30):E4407â€“14. pmid:27412861
Wendt KU, Poralla K, Schulz GE. Structure and Function of a Squalene Cyclase. Science (80-). 1997;277:1811â€“6.
DeJong JHM, Liu Y, Bollon AP, Long RM, Jennewein S, Williams D, et al. Genetic engineering of taxol biosynthetic genes in Saccharomyces cerevisiae. Biotechnol Bioeng. 2006;93(2):212â€“24. pmid:16161138
Urban P, Mignotte C, Kazmaier M, Delorme F, Pompon D. Cloning, yeast expression, and characterization of the coupling of two distantly related Arabidopsis thaliana NADPH-cytochrome P450 reductases with P450 CYP73A5. J Biol Chem. 1997;272(31):19176â€“86. pmid:9235908
SchÃ¼ckel J, Rylott EL, Grogan G, Bruce NC. A Gene-Fusion Approach to Enabling Plant Cytochromes P450 for Biocatalysis. ChemBioChem. 2012;13(18):2758â€“63. pmid:23129550
Zhu M, Wang C, Sun W, Zhou A, Wang Y, Zhang G. Boosting 11-oxo-Î²-amyrin and glycyrrhetinic acid synthesis in Saccharomyces cerevisiae via pairing novel oxidation and reduction system from legume plants. Metab Eng [Internet]. 2018;45:43â€“50. Available from: pmid:29196123
Khakimov B, Kuzina V, Erthmann PÃ˜, Fukushima EO, Augustin J, Olsen CE, et al. Identification and genome organization of saponin pathway genes from a wild crucifer, and their use for transient production of saponins in Nicotiana benthamiana. Plant J. 2015;84(3):478â€“90. pmid:26333142
Miettinen K, Pollier J, Buyst D, Arendt P, Csuk R, Sommerwerk S, et al. The ancient CYP716 family is a major contributor to the diversification of eudicot triterpenoid biosynthesis. Nat Commun. 2017;8:14153. pmid:28165039
Moses T, Pollier J, Faizal A, Apers S, Pieters L, Thevelein JM, et al. Unraveling the Triterpenoid Saponin Biosynthesis of the African Shrub Maesa lanceolata. Mol Plant. 2015;8:122â€“35. pmid:25578277
Glick BR. Metabolic load and heterologous gene expression. Biotechnol Adv. 1995;13(2):247â€“61. pmid:14537822
Kirby J, Romanini DW, Paradise EM, Keasling JD. Engineering triterpene production in Saccharomyces cerevisiaeâ€“Î²-amyrin synthase from Artemisia annua. FEBS J. 2008;275:1852â€“9. pmid:18336574
Zhao Y, Fan J, Wang C, Feng X, Li C. Enhancing oleanolic acid production in engineered Saccharomyces cerevisiae. Bioresour Technol [Internet]. 2018;257:339â€“43. Available from: pmid:29526355