Geranylgeranyl reductase (GGR) is a flavin-containing redox enzyme that hydrogenates a variety of unactivated polyprenyl substrates, which are further processed mostly for lipid biosynthesis in archaea or chlorophyll biosynthesis in plants. To date, only a few GGR genes have been confirmed to reduce polyprenyl substrates in vitro or in vivo.
In this work, we aimed to expand the confirmed GGR activity space by searching for novel genes that function under amenable conditions for microbial mesophilic growth in conventional hosts such as Escherichia coli or Saccharomyces cerevisiae. 31 putative GGRs were selected to test for potential reductase activity in vitro on farnesyl pyrophosphate, geranylgeranyl pyrophosphate, farnesol (FOH), and geranylgeraniol (GGOH). We report the discovery of several novel GGRs exhibiting significant activity toward various polyprenyl substrates under mild conditions (i.e., pH 7.4, T = 37 °C), including the discovery of a novel bacterial GGR isolated from Streptomyces coelicolor. In addition, we uncover new mechanistic insights within several GGR variants, including GGR-mediated phosphatase activity toward polyprenyl pyrophosphates and the first demonstration of completely hydrogenated GGOH and FOH substrates.
These collective results enhance the potential for metabolic engineers to manufacture a variety of isoprenoid-based biofuels, polymers, and chemical feedstocks in common microbial hosts such as E. coli or S. cerevisiae.
Manufacturing of terpenoid-based compounds has been studied extensively in synthetic biology. The two biosynthetic pathways for terpene monomer biosynthesis are the mevalonate and 1-deoxy-d-xylulose 5-phosphate pathways, where pyruvate is ultimately converted into either of the C5 terpene building blocks, isopentenyl pyrophosphate or dimethylallyl pyrophosphate [1, 2]. These monomer units are subsequently fused by various prenyl transferases to make geranyl pyrophosphate (GPP, C10), farnesyl pyrophosphate (FPP, C15), and geranylgeranyl pyrophosphate (GGPP, C20) . The structural diversity of terpenes allows for a broad range of uses in areas including dietary supplements, polymer feedstocks, pharmaceuticals and cosmetics, household cleaners, and fuels [4,5,6,7,8]. Much of this structural diversity is achieved via downstream cyclization and redox steps on GPP, FPP, and GGPP using a plethora of terpene synthases [9,10,11]. Combinations of these core isoprenoid pyrophosphate intermediates serve as starting points for cholesterol biosynthesis, antibiotic biosynthesis, cofactor biosynthesis, and protein prenylation [12,13,14,15,16].
While microbes including E. coli and S. cerevisiae have emerged as robust hosts in the production of terpenoids, producing specially tailored natural products will require the use of novel chemistries and biosynthetic pathways. For example, isoprenoids have been considered as a promising precursor of alternative fuels, but reduction of isoprenoid double bonds is required to decrease the reactivity and sensitivity to oxidation and make them better fuels. Enzymatic alkene hydrogenation, however, is typically assisted by adjacent electron-withdrawing groups as observed in examples including old yellow enzyme, fatty acid enoyl reductases, and enone reductases [17,18,19,20].
Reduction of unactivated substrates like prenyl pyrophosphates typically involves oxidoreductases from the geranylgeranyl reductase (GGR) family. GGR generates fully saturated isoprenoid intermediates in archaeal membrane biosynthesis [21, 22]. In archaea, GGR’s native activity is believed to fully reduce all prenyl groups within the C20 isoprenoid chain of 2,3-di-O-geranylgeranylglyceryl phosphate (DGGGP) before carbon–carbon bond formation of reduced C20 isoprenoid chains form fully reduced C40 precursors needed for membrane synthesis [23, 24]. Moreover, in various organisms such as eukaryotes, bacteria, and archaea, GGRs also have been demonstrated to reduce a variety of prenylated substrates, including chlorophyll, tocopherol, dolichol, and menaquinone [25,26,27,28]. However, very few GGRs have been confirmed as oxidoreductases, and most enzymes having prenyl reductase activity were derived from species that thrive under extremophilic conditions or utilize photosynthesis for energy transduction [25,26,27,28,29,30,31,32]. To date, only two crystal structures have been solved for GGRs from archaeal organisms. Reducing equivalents are thought to be derived from a NAD(P)H/ferredoxin reductase, in which electron transfer is conducted throughout the protein and modulated by a conserved active site cysteine within the cofactor binding domain, located directly behind the FAD isoalloxazine ring .
Biomanufacturing of reduced isoprenoid compounds requires a reductase activity under biologically relevant conditions required by bacterial and yeast strains (i.e., at 30–37 °C, at pH 7). In this study, we sought to increase the diversity space of GGRs by testing several dozen putative GGR sequences across a broad phylogeny, and we proceeded to test their associated substrate promiscuities under conditions ideal for microbial manufacturing (Scheme 1). Herein, we present significant insights on GGR activities that encompass newly confirmed GGR enzymes, novel substrate activities, and promiscuous catalysis.
Results and discussion
Selection and expression of potential GGR candidates
The Interpro database (EMBL-EBI) predicts over 8000 proteins present within the GGR family (InterPro code: IPR011777), with many homologous genes containing sequence identities as low as 20–30%. After sequence alignment, a phylogeny tree includes 1787 sequences of predicted GGR from the InterPro database. A few GGRs within this database have been confirmed by other groups to reduce a wide variety of large prenylated substrates, including GGPP, DGGGP, geranylgeranylchlorophyll, menaquinone, and dolichol [25,26,27,28,29,30,31,32]. To investigate the in vitro prenyl reductase potential of other genes within the GGR family, we selected some with conserved sequence homologies to known GGRs and other more distant sequences. As observed in Fig. 1, it was possible to observe some subgroups with conserved sequences (e.g., Mc, Sa, Pf or Hl, Hv1, Hs). Most of the sequences in the predicted GGR family, however, are very divergent. Our selection was then based on kingdom and diversity of species (e.g., archaea, algae, plant, cyanobacteria, and bacteria), on environmental diversity (e.g., temperatures, pH, aerobic or anaerobic), as well as particular characteristics of some strains (e.g., Corynebacterium terpenotabidum or Gordonia polyisoprenivorans are actinomycetes capable of degrading squalene and rubber reciprocally). A few GGRs were also selected more randomly for their atypical sequences (Fig. 1).
The 31 selected genes were codon optimized for E. coli expression and were all successfully transformed into E. coli. Initial expression attempts were not successful for many proteins using E. coli BL21 (DE3). However, by E. coli BL21 (DE3) strain harboring the commercially available pG-KJE8 plasmid overexpressing several E. coli chaperones, 24 of 31 strains overexpressed soluble proteins at the target masses for each protein, with each protein’s presence in cell lysates confirmed by western blot containing the anti-His tag antibody (Fig. 2). Out of them, only 12 proteins (Af, Cn, Ma, Mb, Mi, Mr2, Pf, Sa, Sc, Syn, Ta, and Tn) were obtained in sufficiently large quantities needed for activity assays after a standard purification and concentration process conducted at pH 7.4.
In vitro activity with isoprenoid alcohols
The 12 soluble proteins successfully isolated were tested for reductase activity on GGOH and FOH, and products obtained after enzymatic incubation were analyzed by GC–MS.
Out of the 12 purified GGRs, five were discovered to enzymatically reduce geranylgeraniol (GGOH). Neat GGOH substrate eluted at a retention time (RT) of 8.4 ± 0.1 min (Fig. 3), with a directly proportional TIC response ranging from 0 to 200 µM (Additional file 1: Figure S1). Upon incubation with any of five putative GGRs isolated from Archaeoglobus fulgidus (AfGGR), Methanocaldococcus infernus (MiGGR), Pyrolobus fumarii (PfGGR), Thermococcus nautili (TnGGR), or Sulfolobus acidocaldarius (SaGGR), several peaks eluting earlier than 8.4 min were observed (Fig. 3). These peaks were assigned to structures of H2-GGOH (RT = 8.1 ± 0.1 min), H4-GGOH (RT = 7.9 ± 0.1 min), and H6-GGOH (RT = 7.7 ± 0.1 min). Moreover, as protein concentration was increased, substrate consumption accelerated (Additional file 1: Figure S2) with a concomitant increase in the formation of the various product peaks (data not showed), confirming enhanced isoprenoid reduction in the presence of higher concentration of enzymes. Out of all enzymes tested, SaGGR was the most active toward GGOH, with a specific activity of at least 50 ± 10 nmol terpenoids reduced per milligram of enzyme per hour (Fig. 4 and Table 2). Typically, 70% of the initial GGOH would be recovered regardless of the varying amounts of reduced product formed. Hence, we assumed that all unrecovered substrate was unreduced, and the turnover numbers presented herein most likely represent a lower bound for reductase activity.
The H2-GGOH and H4-GGOH peaks have respective prevalent ion abundances at 261 and 263 m/z, which can be achieved by loss of a 31 Da [M–CH2OH] fragment during ionization and subsequent formation of a resonance-stabilized singly or doubly reduced geranylgeranyl fragment. Such fragments most likely originate from the prenyl units distal from the alcohol group being reduced first, in accordance with previous mechanistic proposals performed using various substrates on a variety of GGRs [29, 32, 33]. Moreover, the H6-GGOH peak matches with a phytol peak from the NIST database with > 90% probability, further reinforcing a mechanism of serial reduction of substrate beginning with the δ-prenyl group. Interestingly, several GGRs exhibit unknown side-products, with the most prevalent behavior observed between the H2-GGOH and H4-GGOH peaks in Pyrolobus fumarii GGR (RT = 8.0 min) (Fig. 3). This peak contains aberrant patterns for prenyl units within m/z window of 50–100, and we suggest these are H4-GGOH regioisomers in which one or both internal prenyl units are reduced first, which was suggested from the NIST database with > 80% probability (Additional file 1: Figure S3).
Of most noteworthy interest is the product eluted at 7.5 ± 0.1 min RT from assays containing GGRs from Sulfolobus acidocaldarius (Fig. 3). The mass spectra are matched against the 3,7,11,15-teramethylhexadecan-1-ol compound, a complete hydrogenation product of GGOH, in the NIST database with > 88% probability. SaGGR, among others, has been demonstrated to reduce 3 out of 4 prenyl units of GGPP at best as observed in this work and others [29, 33]. Because a complete reduction is not observed in isoprenoid pyrophosphate substrates (Fig. 5) but is observed in the isoprenoid alcohol (Fig. 3), it seems that the absence of phosphate groups might facilitate enhanced diffusion of the α-prenyl group to the flavin reducing site in the alcohol substrates, leading to a fully reduced product. To our knowledge, this is the first evidence of any nonnative isoprenoid substrate undergoing full reduction by any known or putative GGR enzyme.
Similarities in reducing activity were also prevalent using farnesol as a substrate. The unreduced FOH substrate eluted with a RT of 8.0 ± 0.1 min, with the putative singly (H2-FOH) and doubly reduced (H4-FOH) farnesol eluting at 7.6 ± 0.1 min and 7.4 ± 0.1 min, respectively (Fig. 6). Farnesol ionization was also directly proportional to concentration ranging from 0 to 200 µM (Additional file 1: Figure S1). The accompanying mass spectrum for H2-FOH reveals a similar ionization pattern to that observed in H2-GGOH via the prevalence of a strong 193 m/z peak. This parallels the H2-GGOH peak pattern containing one less prenyl group (m/z = 70 Da). This suggests that the terminal isoprenoid unit is also reduced first in farnesol, conserving the enzymatic reduction mechanism regardless of substrate. The H4-FOH peak at 7.4 min more closely resembles the H6-GGOH peak, with identical peak groupings near the 71, 81, and 123 m/z parent fragments.
Unlike GGOH, all GGRs appeared to have similar levels of FOH products under standard assay conditions, exhibiting an average specific activity of 7 ± 2 nmol terpenoid groups reduced mg−1 enzyme h−1 (Fig. 4 and Table 1). Notably, reduction patterns in TnGGR on FOH differ slightly from the other GGRs under standard assay conditions, as its major product is H2-FOH instead of H4-FOH (Fig. 4). No fully reduced farnesol peaks were observed under standard assay conditions nor at enzyme concentrations as high as 150 µM at pH 7.4. However, SaGGR generated a modest amount of H6-FOH when incubated under the enzyme’s optimal conditions at 50 °C and pH 5.5  (Additional file 1: Figure S4).
Compared to GGOH, emergent side products are less prevalent in the farnesol TICs. Whereas multiple peaks were observed between the singly and doubly reduced GGOH (Fig. 3), a single prevalent peak at 7.8 min elutes between FOH and H2-FOH, mainly observed when incubated in the presence of Pf and MiGGR (Fig. 6). The associated mass spectrum is tentatively assigned to a regioisomer of H2-FOH where the middle prenyl unit is reduced first (Additional file 1: Figure S5). The spectrum matches the NIST database for (E)-3,7,11-trimethyldodeca-2,10-dien-1-ol with a probability of 85% (Additional file 1: Figure S5). Many of the aberrant mass groupings between 50 and 100 m/z (Additional file 1: Figure S5) parallel those observed in the GGOH reaction incubated with PfGGR (Additional file 1: Figure S3). However, the 7.8-min peak does not contain the m/z 193 ion. This ion could be formed by cyclization of a [M–CH2OH] fragment containing a reduced terminal prenyl group and is absent in products where the middle group is reduced first due structural rigidity associated with the remnant α- and γ-prenyl groups. This observation, coupled with the aberrant TIC product profile observed with GGOH, suggests a promiscuous mechanism in which PfGGR has been observed to reduce prenyl monomers out of order with respect to their polymeric structural order.
In vitro activity with isoprenoid pyrophosphates
The 12 soluble GGRs successfully purified were tested for reductase activity on FPP and GGPP, and products were detected by LC–MS–TOF. Both farnesyl pyrophosphate (FPP, m/z = 381.123 ± 0.001 Da) and geranylgeranyl pyrophosphate (GGPP, m/z = 449.183 ± 0.002 Da) standards eluted with a retention time of 1.70 ± 0.05 min (Fig. 5); both substrates produced linear standard curves over a concentration range of 0–120 µM (Additional file 1: Figure S6). When incubated with GGR under standard assay conditions, reduced isoprenoid products were observed to co-elute with fully oxidized substrate under isocratic LC conditions. Therefore, only normalized LC–MS–TOF spectra were utilized to distinguish the relative levels reduced and oxidized compounds that co-elute after incubating with GGRs isolated from various species (Additional file 1: Table S1).
Interestingly, all proteins in this study discovered to enzymatically reduce prenyl pyrophosphates revealed co-eluting side products indicative of substrate or product hydrolysis of one phosphate moiety (Fig. 5). Increased abundances of farnesyl monophosphate (FP, m/z = 301.177 ± 0.001 Da) or geranylgeranyl monophosphate (GGP, m/z = 369.213 ± 0.003 Da) only emerged when incubated with enzyme; minimal hydrolysis was observed in GGPP or FPP standards and relative GGPP/GGP and FPP/FP ratios remained constant as a function of time in negative controls ran without enzyme (Additional file 1: Figure S7). Structural studies of SaGGR crystallized with GGPP revealed three distinct substrate binding modes with varying degrees of phosphorylation within each binding position . Within the catalytically relevant binding mode, both phosphate moieties are resolved. In the other two binding modes, however, either one or zero phosphate group was structurally resolved. This was attributed to dephosphorylation during the crystallization process . Herein, we observed a time-dependent emergence of hydrolyzed monophosphate products via LC–MS–TOF; yet it still requires further characterization how the enzyme facilitates this phenomenon while conducting substrate reduction.
Reductase activity on FPP and GGPP varied from what was observed on alcohol substrates (Fig. 7). Indeed, none of the GGRs tested could significantly reduce all vinyl groups within FPP or GGPP even when GGRs were incubated under the optimal condition for enzyme activity (at 50 °C and pH 5.5) (data not shown). Out of the five GGRs found to reduce FOH, only PfGGR, TnGGR, and SaGGR could reduce FPP. On the other hand, three GGRs isolated from Streptomyces coelicolor (ScGGR), Methanosarcina acetivorans (MaGGR), and Thermoplasma acidophilum (TaGGR) were found to reduce GGPP along with the five GGRs demonstrating reductase activity toward GGOH. Most GGRs that have been isolated thus far were from archaea; to our knowledge, ScGGR is the first bacterial GGR demonstrated to reduce GGPP. Due to unexpected hydrolysis of one phosphate moiety under standard assay conditions, specific activities for reduction were not quantified for any GGR. However, relative reductase activities can be gleaned by quantifying the proportion of reduced and unreduced compounds present within intact or hydrolyzed mass groupings (Fig. 7). For example, the relative ion intensities of each singly reduced product (H2-FPP or H2-FP) present is normalized to the sum of FPP, H2-FPP, H4-FPP, H6-FPP, FP, H2-FP, H4-FP, and H6-FP extracted ion intensities.
Interestingly, all GGRs revealed a larger proportion of reduced products present as hydrolyzed moieties than non-hydrolyzed moieties (Fig. 7). To discern a correlation between enzymatic phosphate hydrolysis and enzymatic reduction of prenyl monophosphates, we assayed PfGGR and SaGGR as a function of time, as these enzymes are representative of low and high amounts of monophosphate found after standard assay incubation, respectively (Additional file 1: Figure S7). Indeed, substrate hydrolysis appears to react as a first-order exponential process which occurs more quickly in SaGGR than for PfGGR (Additional file 1: Figure S8). Moreover, the relative abundances of reduced monophosphate products increase over time in both assays, indicating that most GGRs can still reduce FP or GGP as a substrate during substrate hydrolysis (Additional file 1: Figure S9). Such an inference is reasonable considering many GGRs assayed enzymatically reduce terpenoid alcohols and pyrophosphates.
Promiscuous hydrolysis complicates any interpretations regarding which enzymes are most active toward a given substrate due to the inability to quantify the MS response of terpenoid phosphates. However, it can be inferred that all GGRs can reduce between 5 and 10 nmol prenyl groups of FPP or GGPP mg−1 enzyme h−1. The turnover number would be modestly elevated for GGPP reduction, as all C20 species are extracted as some partially reduced product within error after 1 h. Such turnover numbers are in line with other reports on GGRs with a variety of substrates [32, 33].
Structural insights and mechanistic implications
Several synthetic approaches are currently being explored to perform selective hydrogenation on a few substrates [34,35,36]. Biological systems such as enoyl-CoA reductase and old yellow enzyme exhibit a similar oxidoreductase activity to GGR, yet benefit from active sites that enhance the electron-withdrawing nature of α,β-unsaturated carbonyl substrates [17,18,19,20]. Patented ene-reductases utilizing old yellow enzyme as a scaffold enhance reductase activity on a variety of substrates by evolving active sites complementary to a variety of electron withdrawing groups among a diverse variety of α,β-unsaturated substrates . However, an evolved GGR active site designed for isoprenoid reduction would probably require significant divergence from these scaffolds since they do not utilize electron-withdrawing activation for alkene reduction .
Of the eight proteins that were identified as GGRs active toward terpenoid alcohols and/or terpenoid pyrophosphates, five (Sa-, Pf-, Af-, Mi-, and TnGGRs) were isolated from archaeal organisms that optimally thrive under hyperthermophilic conditions (i.e., T ≥ 80 °C). SaGGR, TaGGR, and AfGGR have been identified to reduce various large intermediates (i.e., larger than 20 carbons) associated with archaeal lipid biosynthesis, with GGPP or GGOH serving as the smallest substrates known to undergo prenyl reduction [27, 29, 32]. In this study, we have significantly expanded the known GGR substrate activity profiles, demonstrating multiple prenyl group reduction in GGOH and FOH within all five hyperthermophilic GGRs.
In addition to the five GGRs active on alcohols, TaGGR, MaGGR, and ScGGR also sufficiently reduced GGPP or GGP (Figs. 5 and 7). However, only PfGGR, SaGGR, and TnGGR were found to reduce the smaller FPP or FP substrates. Because the relative amount of H2-, H4-, and H6-GGOH increase in relative abundances within the monophosphate mass groupings relative to the pyrophosphate mass groupings, it can be inferred that prenyl monophosphates are also substrates reduced by several GGRs (Additional file 1: Figure S9). This seems suitable given the ability of several GGRs to reduce prenyl alcohols.
A structural alignment of all eight active GGRs reveals very little commonalities among all protein sequences with known crystal structures: SaGGR and TaGGR, with PfGGR ca. 46% identical to SaGGR and MaGGR, MiGGR, and AfGGR ca. 40–46% identical to TaGGR (Fig. 8). SaGGR and TaGGR contain three domains: an FAD binding domain, a catalytic domain, and a C-terminal domain [29, 30]. While sequence identities remain low among all demonstrably active GGRs, certain key structural motifs remain conserved within their predicted FAD binding domains and catalytic domains. Of the two known crystal structures of active GGRs, both contain an active site cysteine (Cys47 in SaGGR; Cys45 in TaGGR) thought to serve as a critical redox modulator within the active site during reduction. All GGRs shown to reduce either isoprenoid alcohols or pyrophosphates contain this critical cysteine within their cofactor binding domains, suggestive of a conserved electron transfer mechanism. In addition, all sequences predicted catalytic domains contain the YXWXFP (SaGGR residues 215–220) and GGG motifs (SaGGR residues 298–300) believed to modulate substrate interactions and assist in substrate diffusion through the reduction center.
Protein structures of aligned sequences were predicted using either SaGGR or TaGGR as a template. While there is a fair amount of expected structural divergence among the structures’ surfaces, a comparison of the active sites reveals a fair degree of similarity in topology (Fig. 9). However, some of the structural motifs strictly conserved among all archaeal GGRs exhibit significant divergence within ScGGR, the only known GGR to be isolated from a bacterial organism. While all archaeal GGRs studied in this work possess a YXWXFPX7-8GXG motif, the terminal glycine is mutated to isoleucine in ScGGR. Even more interestingly, the GGG motif has significantly diverged to REG in ScGGR. In several GGRs from photosynthetic organisms with demonstrated capability to reduce prenylated chlorophyll, Rhodobacter sphaeroides, Synechocystis sp. PCC 6803, and Arabidopsis thaliana, this motif was found to be GEG [26, 38, 39]. It seems that non-archaeal GGRs utilize preferentially charged residues within this critical catalytic region to either enhance polar interactions on prenylated substrates containing polar groups or to introduce critical hydrogen-bonding interactions that help maintain the integrity of the substrate tunnel during reduction (Fig. 9).
Mechanistic interpretations from other groups propose that the prenyl group closest to the pyrophosphate moiety (α-prenyl group) remains oxidized in GGPP and FPP. This observation additionally applies to their monophosphate counterparts in this work, FP and GGP. All enzymes tested to date seem to conserve this characteristic of avoiding reduction at the α-position on phosphate intermediates, aligning with current paradigms that auxiliary prenyl reductases are responsible for reducing this group in archaea and eukaryotes .
To our knowledge, full isoprenoid reduction by GGR has only been observed with its natural C40 isoprenoid substrate DGGGP. In this work, we observed full reduction for the first time on smaller (i.e., C20 or C15) isoprenoid alcohol substrates, namely GGOH and FOH with SaGGR (Figs. 3, 4, 6, Additional file 1: Figure S4). Interestingly, the absence of phosphate groups appears to assist in full substrate reduction. Analysis of the catalytically relevant binding mode of GGPP in SaGGR reveals that binding site residues His55 and Asn90 could provide hydrogen bonding interactions with phosphate moieties that could prevent the α-prenyl group from being reduced . Alcohol substrates may not interact as strongly with these residues, facilitating a degree of full reduction unobserved in pyrophosphate substrates. Why some enzymes reduce isoprenoid alcohol and pyrophosphate substrates, while others only reduce isoprenoid pyrophosphates requires further structural characterization.
In this study, we have significantly expanded the possible activities among proteins demonstrated to enzymatically reduce prenyl pyrophosphates or prenyl alcohols. We have demonstrated (1) the discovery of four novel protein sequences (PfGGR, MiGGR, ScGGR, and TnGGR) that have confirmed GGR activity in vitro in addition to expanded observed activities among previously characterized GGRs; (2) that several GGRs can reduce C15 terpenoid substrates, substrates smaller than reported substrates for GGR activity; (3) the complete reduction of double bonds on any C20 or C15 isoprenoid using SaGGR; (4) reductase activity on terpenoid monophosphates formed from hydrolysis of pyrophosphate substrates under reducing conditions in vitro; (5) the quantification of reductase specific activity on terpenoid alcohols; and (6) the confirmed isoprenoid reductase activity of the second known non-archaeal enzyme, as observed in the GGR isolated from Streptomyces coelicolor.
This demonstration of protein expression and reductase activity at neutral pH and low temperature highlights their potential suitability for integration into S. cerevisiae or E. coli. Moreover, the confirmation of reduction on C15 isoprenoids instantly expands the metabolic engineering potential for organisms producing sterol and squalene-derived isoprenoids. There are still unresolved issues to address for a direct application of these newly discovered GGRs to manufacture reduced isoprenoids. For example, more engineering will be needed on these enzymes to avoid enzymatic hydrolysis of isoprenoid pyrophosphates and to improve their activities especially at mesophilic condition. Nonetheless, this study demonstrated significant substrate promiscuity among these GGRs and could potentially open new pathways for isoprenoid-based polymers, chemicals, or biofuels by allowing for upstream reduction of various intermediates within the heavily utilized MEV or DXP terpene biosynthesis pathways.
All chemicals and reagent were purchased from Sigma-Aldrich (St. Louis, MO), unless otherwise indicated. (E,E)-farnesol was purchased from Alfa Aesar (Haverhill, MA) and glycerol from VWR (Westchester, PA). Solvents for high performance liquid chromatography (HPLC) were purchased from HoneyWell Burdick and Jackson (Morristown, NJ) and were of HPLC grade or higher. Ammonium carbonate (30–33% NH3 basis) was purchased from Fluka Analytical Sigma-Aldrich (St. Louis, MO). Restriction enzymes and polymerases were purchased from New England Biolabs (Ipswich, MA).
Sequence analysis and GGR homology
Multiple sequence alignments for potential GGR hits were generated using MUSCLE v. 3.8.31 and visualized using Geneious 7.0.6 [41, 42]. Sequences were curated manually, and phylogeny trees were computed using the maximum likelihood tree within the RAxML Software package, v. 8.1.24 under the LG plus gamma model of evolution (PROTGAMMALG in the RAxML model section) . The MRE-based bootstrapping criterion was automatically determined for phylogeny tree construction. Annotation of the tree was performed in Itol . After verification of GGR activity, the active enzymes underwent a second multiple sequence alignment and modeled for their predicted protein structures via SWISS-MODEL-PDB using either SaGGR or TaGGR as templates . Active site geometries and local structures for all proteins were visualized using Chimera .
Plasmid synthesis and transformation
The gene encoding SaGGR was amplified by PCR from the pSKB3-SaGGR plasmid using the forward (5′-GATATACATATGAAGGAACTTAAATATGACGTTCTG-3′) and reverse (5′-GTCGACGGAGCTCGAACTTAAACTTTTGTTAAACTCTGTTAGAAC-3′) primers synthesized by Integrated DNA Technologies . The PCR fragment was digested at the NdeI and SacI restriction sites and cloned into the pET-24a vector using the rapid DNA ligation kit (Roche). All other putative GGR genes were synthesized by GeneWiz (NJ, USA) and similarly cloned into the pET-24a vector at the same restriction sites. All gene constructs are available through the JBEI registry at http://public-registry.jbei.org (Table 1 and Additional file 1: Table S2).
Ten nanogram of each plasmid was transformed by heat shock at 42 °C for 1 min into chemically competent E. coli BL21 cells harboring the pG-KJE8 plasmid encoding DnaK, DnaJ, GrpE, GroES, and GroEL protein chaperones (Takara Bio Inc., Shiga, Japan). Transformed cells were recovered in 1 mL of Lysogeny Broth (LB) medium (VWR) and incubated for 1 h at 37 °C with shaking at 200 rpm. Following recovery, cells were plated on LB Agar containing 50 mg/L of kanamycin (VWR) and 30 mg/L of chloramphenicol (VWR), and incubated overnight at 37 °C. Select colonies were grown overnight in LB medium containing 50 mg/L of kanamycin and 30 mg/L of chloramphenicol and stored in 20% glycerol (VWR) at − 80 °C for future use.
Cell culture, protein expression, and protein purification
Overnight seed cultures of 1 mL each were inoculated into 400 mL of Terrific Broth (TB) medium supplemented with 50 mg/L kanamycin and 30 mg/L chloramphenicol and incubated at 37 °C and 200 rpm. At an OD600 of 0.2–0.3, chaperone overexpression was induced with 5 ng/mL tetracycline (VWR) and 2.5 mM arabinose (Sigma-Aldrich). After the OD600 reached ≥ 1.0, GGR expression was induced with 0.1 mM IPTG (VWR) and incubated at 18 °C overnight. Cells were pelleted at 6000×g for 10 min and immediately lysed using 20 mM phosphate buffer, pH 8.0 containing 1 mg/mL lysozyme, 20 mM imidazole, 200 mM NaCl, and 0.1 mM PMSF protease inhibitor (Sigma-Aldrich). After sonication for 10 min, the remaining cell debris was pelleted at 15,000×g for 45 min.
Protein expression was tested for each construct using SDS-PAGE and Western blot. For SDS-PAGE analysis, protein samples were normalized for concentration using absorbance at 280 nm. Lysates were diluted with 2× SDS loading dye buffer (Life Technologies, CA, USA) containing 10 mM DTT (Sigma-Aldrich) and incubated at 98 °C for 20 min. 10 µL of denatured lysate samples was loaded onto an 8–16% Tris–Glycine–SDS gradient gel (Bio-Rad), and separated using a voltage of 180 V in Tris–Glycine–SDS running buffer (Bio-Rad). Gels were either directly stained using GelCode Blue Safe Protein Stain (Thermo-Fisher) or transferred to a nitrocellulose membrane using the trans-Blot Turbo system (Life Technologies, CA, USA) for analysis by Western blot. Membranes were washed in TBS buffer (50 mM Tris, 150 mM NaCl, pH 7.4) and blocked overnight at 4 °C with 25 mL of 3% BSA in TBS-Tween20 (Sigma-Aldrich). The monoclonal mouse anti-His primary antibody (Sigma-Aldrich) was diluted 5000-fold, and an alkaline phosphatase-conjugated goat anti-mouse secondary antibody was diluted 10,000-fold in TBS-Tween20 containing 1% BSA. Membranes were incubated with antibodies for 1 h each at room temperature and washed three times in TBS-Tween20 after each antibody incubation. The membrane was then incubated in 10 mL of SigmaFast BCIP/NBT Alkaline Phosphatase detection solution (Sigma-Aldrich) for 10 min.
To further characterize those putative GGRs that showed significant protein expression, the cells harboring them were cultured in 400 mL of TB-Kan/Cm media and lysed as previously described. Their respective crude lysates were loaded directly onto a 1 mL HisTrap FastFlow column (GE Healthcare), washed with ten column volumes of 20 mM phosphate buffer containing 20 mM imidazole and 200 mM NaCl at pH 7.4, then eluted with the same buffer containing 240 mM imidazole. For enzyme kinetics, purified enzymes were buffer exchanged using 20 mM phosphate buffer at pH 7.4 and concentrated to 200–800 μM using 30 KDa molecular weight cutoff spin concentrators (EMD Millipore). Purified proteins were stored in 10% (v:v) glycerol and snap frozen in liquid nitrogen. Protein purity and sizes were verified by SDS-PAGE and protein concentrations were quantified by absorbance at 280 nm using each protein’s calculated extinction coefficient via the ExPASY ProtParam tool.
In vitro enzyme kinetics assays
Validation of enzymatic substrate reduction was determined by incubating all assays in triplicate for each respective substrate and putative GGR for 1 h at 37 °C. All assays were performed at pH 7.4 in 100 mM sodium phosphate buffer containing 30–150 µM enzyme, 200 µM FAD (Sigma-Aldrich), and 65 mM sodium dithionite (Sigma-Aldrich). Standard assays for alcohol reduction were incubated with 100 µM enzyme and 500 µM (E,E)-farnesol (Alfa-Aesar) or (E,E,E)-geranylgeraniol (Sigma-Aldrich); pyrophosphate assays were performed at 100 µM FPP or GGPP (Sigma-Aldrich). Alcohol-based assays were quenched by liquid extraction using a 3:1 (v:v) LC-grade ethyl acetate solution containing 100 µM dodecanol as a GC internal standard (Sigma-Aldrich). The organic layer was extracted and stored at − 20 °C until analysis by GC–MS. Pyrophosphate assays were similarly quenched using LC-grade n-butanol (Sigma-Aldrich) 1:1 (v:v) and centrifuged at 15,000×g for 2 min. The n-butanol layer was dried for 45 min at ambient temperature using a Labconco speedvac, reconstituted in 25 µL of a 62:38 (v:v) acetonitrile/50 mM ammonium carbonate solution, and stored at − 20 °C until further analysis by LC–MS–TOF . Characterization of enzymatic hydrolysis of isoprenoid pyrophosphate substrates by SaGGR and PfGGR was performed by quenching the enzyme reactions at 0, 2, 5, 10, 20, 40, and 60 min of incubation.
Analysis of alcohol reduction by GC–MS
Product identification and quantification of farnesol and hydrofarnesol derivatives were modified from previous detection methods . All GC–MS analyses were determined using an Agilent 6890 gas chromatography instrument coupled to an Agilent 5973 mass selective detector. 1 µL of extracted samples was injected in splitless mode onto an Agilent CycloSil-B column, with helium used as a carrier gas flowing at 1.0 mL/min. Following injection, the oven was held at 50 °C for 30 s, then increased to 175 °C at 35 °C/min. Farnesol and hydrofarnesols were resolved by increasing the temperature 4 °C/min up to 200 °C, then increased to 300 °C at a rate of 35 °C/min where it was held for 1.5 min. Geranylgeraniol and its hydrogenated derivatives were analyzed using the same injection method. After injection, the oven was held at 50 °C for 30 s, then increased to 235 °C at 35 °C/min. Hydrogeranylgeraniols were separated by increasing the oven temperature 4 °C/min to 250 °C, then ramped to 300 °C at a rate of 35 °C/min where it was held for 1.5 min.
The EI–MS detection was initiated after a solvent delay of 5.0 min. Detection and classification of hydrofarnesols were performed in scan mode at 9.8 scans/s ranging from 50 to 250 m/z in positive ion mode. For geranylgeraniol, the same scan parameters were implemented except for the mass range, which was expanded to 50–300 m/z in positive ion mode. The electron multiplier voltage was set to a gain factor of 1, with the MS ion source and quadrupole set to 230 °C and 150 °C, respectively.
Total ion chromatograms (TIC) were integrated using Agilent Technologies Masshunter software, version 6. Product formation was determined from the TIC area for C15 or C20 alcohol products eluting at each respective retention time. Absolute product concentrations were determined from standard curves (0–200 µM) of either farnesol or geranylgeraniol assuming the TIC area of each reduced product ionizes with an equivalent efficiency to that of the unreduced substrate (Additional file 1: Figure S1). Subsequently, enzyme turnover numbers for isoprenoid reduction were calculated as the total number of nanomoles of prenyl units reduced per milligram of enzyme in 1 h.
Analysis of pyrophosphate reduction by LC–MS–TOF
The separation of FPP, GGPP, and their reduced forms was conducted on a ZIC-pHILIC column (150 mm length, 2.1 mm internal diameter, and 5 µm particle size, Merck) using an Agilent Technologies 1200 Series Rapid Resolution high-performance liquid chromatography (HPLC) system. Solvents for HPLC were purchased from HoneyWell and were of HPLC grade or higher. The mobile phases used for this analysis were (A) 50 mM ammonium carbonate (Fluka, 30–33% NH3 basis) in water and (B) acetonitrile. Analytes were eluted isocratically with a mobile phase composition of 62% B at a flow rate of 0.2 mL/min. The total run time of the method was 6.5 min. The temperature of the sample tray was maintained at 6 °C using an Agilent FC/ALS Thermostat. The column compartment was set to 40 °C. A sample injection volume of 2 µL was used throughout .
The HPLC system was coupled to an Agilent Technologies 6210 time-of-flight mass spectrometer (LC–TOF–MS) by a 1/3 post-column split. Contact between both instrument set-ups was established using a LAN card to trigger the MS into operation upon the initiation of a run cycle from the MassHunter workstation (Agilent Technologies). Electrospray ionization (ESI) was conducted in the negative ion mode and a capillary voltage of − 3500 V was utilized. MS experiments were carried out in full scan mode, at 0.86 spectra/second for the detection of [M−H]− ions. The instrument was tuned for a range of 50–1700 m/z. Prior to LC–TOF–MS analysis, the TOF–MS was calibrated via an ESI-L low concentration tuning mix (Agilent Technologies).
Data acquisition and processing were performed by the Agilent Technologies MassHunter software package. Product formation was determined using extracted ion chromatogram abundances (± 0.02 Da) for each molecule’s [M−H]− mass (Additional file 1: Table S1). Substrate and product hydrolysis of SaGGR and PfGGR was characterized as a function of time by measuring the relative ratios of prenyl pyrophosphates (FPP/GGPP and reduced products) and monophosphates (FP/GGP and reduced products) at quenched fractions collected at 0, 2, 5, 10, 20, 40, and 60 min. Relative reductase reactivity among GGRs was determined by measuring the fractional abundance of singly, doubly, or triply reduced products to the total ion abundance present for intact and hydrolyzed moieties . Integrated areas for hydrolyzed monophosphate products were assumed to have the same ionization intensities as their pyrophosphate counterparts, as determined by their standard curves measured from 0 to 120 µM (Additional file 1: Figure S6).
liquid chromatography–mass spectrometry
time of flight
gas chromatography–mass spectrometry
total ion chromatogram
Rohmer M. The discovery of a mevalonate-independent pathway for isoprenoid biosynthesis in bacteria, algae and higher plants. Nat Prod Rep. 1999;16(5):565–74.
Ajikumar PK, Xiao WH, Tyo KEJ, Wang Y, Simeon F, Leonard E, Mucha O, Phon TH, Pfeifer B, Stephanopoulos G. Isoprenoid pathway optimization for taxol precursor overproduction in Escherichia coli. Science. 2010;330(6000):70–4.
Murakami M, Shibuya K, Nakayama T, Nishino T, Yoshimura T, Hemmi H. Geranylgeranyl reductase involved in the biosynthesis of archaeal membrane lipids in the hyperthermophilic archaeon Archaeoglobus fulgidus. FEBS J. 2007;274(3):805–14.
Sato S, Murakami M, Yoshimura T, Hemmi H. Specific partial reduction of geranylgeranyl diphosphate by an enzyme from the thermoacidophilic archaeon Sulfolobus acidocaldarius yields a reactive prenyl donor, not a dead-end product. J Bacteriol. 2008;190(11):3923–9.
Mizoguchi T, Isaji M, Yamano N, Harada J, Fujii R, Tamiaki H. Molecular structures and functions of chlorophylls-a esterified with geranylgeranyl, dihydrogeranylgeranyl, and tetrahydrogeranylgeranyl groups at the 17-propionate residue in a diatom, Chaetoceros calcitrans. Biochemistry. 2017;56(28):3682–8.
Xu QP, Eguchi T, Mathews II, Rife CL, Chiu HJ, Farr CL, Feuerhelm J, Jaroszewski L, Klock HE, Knuth MW, et al. Insights into substrate specificity of geranylgeranyl reductases revealed by the structure of digeranylgeranylglycerophospholipid reductase, an essential enzyme in the biosynthesis of archaeal membrane lipids. J Mol Biol. 2010;404(3):403–17.
Isobe K, Ogawa T, Hirose K, Yokoi T, Yoshimura T, Hemmi H. Geranylgeranyl reductase and ferredoxin from Methanosarcina acetivorans are required for the synthesis of fully reduced archaeal membrane lipid in Escherichia coli cells. J Bacteriol. 2014;196(2):417–23.
Ogawa T, Isobe K, Mori T, Asakawa S, Yoshimura T, Hemmi H. A novel geranylgeranyl reductase from the methanogenic archaeon Methanosarcina acetivorans displays unique regiospecificity. FEBS J. 2014;281(14):3165–76.
Keller Y, Bouvier F, D’Harlingue A, Camara B. Metabolic compartmentation of plastid prenyllipid biosynthesis—evidence for the involvement of a multifunctional geranylgeranyl reductase. Eur J Biochem. 1998;251(1–2):413–7.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.
Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, Kiefer F, Cassarino TG, Bertoni M, Bordoli L, et al. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42(W1):W252–8.
Rodriguez S, Kirby J, Denby CM, Keasling JD. Production and quantification of sesquiterpenes in Saccharomyces cerevisiae, including extraction, detection and quantification of terpene products and key related metabolites. Nat Protoc. 2014;9(8):1980–96.
CWM, FM, AC and TSL conceived the project. CWM, FM, BMG, TSL drafted the manuscript, reviewed and edited the final manuscript. CWM, FM, BMG, EEKB, VTB, AVR, and RA performed the experiments and analyzed data. All authors read and approved the final manuscript.
We thank Dr. Jay Keasling for support during the early stage of this project and for his editing assistance.
TSL has financial interest in Maple Bio.
Availability of data and materials
Consent for publication
Ethics approval and consent to participate
This work was part of the DOE Joint BioEnergy Institute (https://www.jbei.org) supported by the US Department of Energy, Office of Science, Office of Biological and Environmental Research, through contract DE-AC02-05CH11231 between Lawrence Berkeley National Laboratory and the US Department of Energy. This work, funded by Total, was performed as part of a collaborative program between JBEI and Total.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Corey W. Meadows and Florence Mingardon contributed equally to this work
Authors and Affiliations
Joint BioEnergy Institute, 5885 Hollis Street, 4th floor, Emeryville, CA, 94608, USA
Corey W. Meadows, Brett M. Garabedian, Edward E. K. Baidoo, Veronica T. Benites, Andria V. Rodrigues & Taek Soon Lee
Biological Systems & Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
Corey W. Meadows, Brett M. Garabedian, Edward E. K. Baidoo, Veronica T. Benites, Andria V. Rodrigues & Taek Soon Lee
Total Raffinage Chimie, 2 Pl. Jean Millier, 92400, Courbevoie, France
Florence Mingardon, Raya Abourjeily & Angelique Chanal
TIC for neat GGOH (RT = 8.4 min, top) and FOH (RT = 8.0 min, middle) substrates. The standard curve for quantifying farnesol (circles) and geranylgeraniol (squares) by GC–MS (bottom) exhibited a linear response for both substrates between 0 and 200 µM. Figure S2. Verification of accelerated substrate reduction as a function of enzyme concentration for GGOH (left) and FOH (right) for the Af (circles), Mi (squares), Tn (filled triangles), Sa (filled upside down triangles), and Pf (unfilled triangles) GGR enzymes. Specific activities are quoted in Table 1. Figure S3. Comparison of mass spectra between a side product containing one internal prenyl group reduced within H4-GGOH with an 8.0 min retention time (Top, black) and the assigned product with the terminal prenyl group reduced in H4-GGOH eluting at 7.7 min (Bottom, green). Figure S4. (Top) normalized TIC of farnesol activity assay incubated for 2 h with SaGGR at 50 °C, pH 5.5 showing a modest abundance of fully reduced farnesol (RT = 7.0 ± 0.1 min). For reference, FOH and H4-FOH elute at retention times of 7.9 and 7.3 min, respectively. H2-FOH (RT = 7.5 min) was not observed in any quantifiable abundance. All substrate and cofactor concentrations were held constant. Figure S5. Comparison of mass spectra between the middle prenyl group reduced within the putative H2-FOH side product eluting at 7.8 min retention time (Top, black) and the assigned product with the terminal prenyl group reduced in H2-FOH eluting at 7.6 min (Bottom, green). Figure S6. Standard curve for quantifying FPP (circles) and GGPP (squares) by LC–MS–TOF. Figure S7. MS-TOF Spectrum of 100 µM FPP and GGPP standards (Top). (Bottom) relative abundances of GGPP and GGP (left) or FPP and FP (middle) after incubation under standard assay conditions; negative controls containing all assay components without enzyme (right) rule out the possibility of spontaneous hydrolysis of substrate, as the ratio of pyrophosphate (dark gray) to monophosphate (light gray) products remain constant as a function of time. Reduced products within each mass grouping are included in the total abundance. Figure S8. Demonstration that first-order substrate hydrolysis catalyzed by either SaGGR (squares) or PfGGR (triangles) in either FPP (left) or GGPP (right). The no enzyme control (circles) contained all assay components except enzyme. Figure S9. Timecourse comparison of the standard assay for SaGGR (top row) and PfGGR (bottom row) on either GGPP or FPP substrates. Pyrophosphate abundances are shown in the left column and monophosphate abundances are shown in the right column under each substrate. Products are expressed as having zero reductions (blue), one reduction (green), two reductions (red), or three reductions (orange). Table S1. Masses used to analyze various products formed from GGR standard assays incubated with prenyl pyrophosphates. aMasses reported are for the deprotonated [M−H]− parent ion in negative mode detection. Table S2. Table of plasmids used in the present study. The strains harboring individual plasmid are available at the public registry of the Joint BioEnergy Institute (https://public-registry.jbei.org/) under the ID’s listed in the righthand column.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Meadows, C.W., Mingardon, F., Garabedian, B.M. et al. Discovery of novel geranylgeranyl reductases and characterization of their substrate promiscuity.
Biotechnol Biofuels11, 340 (2018). https://doi.org/10.1186/s13068-018-1342-2