Fatty acyl-AMP ligases (FAALs) channelize fatty acids towards biosynthesis of virulent lipids in mycobacteria and other pharmaceutically or ecologically important polyketides and lipopeptides in other microbes. They do so by bypassing the ubiquitous coenzyme A-dependent activation and rely on the acyl carrier protein-tethered 4'-phosphopantetheine (holo-ACP). The molecular basis of how FAALs strictly reject chemically identical and abundant acceptors like coenzyme A (CoA) and accept holo-ACP unlike other members of the ANL superfamily remains elusive. We show that FAALs have plugged the promiscuous canonical CoA-binding pockets and utilize highly selective alternative binding sites. These alternative pockets can distinguish adenosine 3',5'-bisphosphate-containing CoA from holo-ACP and thus FAALs can distinguish between CoA and holo-ACP. These exclusive features helped identify the omnipresence of FAAL-like proteins and their emergence in plants, fungi, and animals with unconventional domain organizations. The universal distribution of FAALs suggests that they are parallelly evolved with FACLs for ensuring a CoA-independent activation and redirection of fatty acids towards lipidic metabolites.

In Lysobacter enzymogenes OH11, RpfB1 and RpfB2 were predicted to encode acyl coenzyme A (CoA) ligases. RpfB1 is located in the Rpf gene cluster. Interestingly, we found an RpfB1 homolog (RpfB2) outside this canonical gene cluster, and nothing is known about its functionality or mechanism. Here, we report that rpfB1 and rpfB2 can functionally replace EcFadD in the Escherichia colifadD mutant JW1794. RpfB activates long-chain fatty acids (n-C16:0 and n-C18:0) for the corresponding fatty acyl-CoA ligase (FCL) activity in vitro, and Glu-361 plays critical roles in the catalytic mechanism of RpfB1 and RpfB2. Deletion of rpfB1 and rpfB2 resulted in significantly increased heat-stable antifungal factor (HSAF) production, and overexpression of rpfB1 or rpfB2 completely suppressed HSAF production. Deletion of rpfB1 and rpfB2 resulted in increased L. enzymogenes diffusible signaling factor 3 (LeDSF3) synthesis in L. enzymogenes Overall, our results showed that changes in intracellular free fatty acid levels significantly altered HSAF production. Our report shows that intracellular free fatty acids are required for HSAF production and that RpfB affects HSAF production via FCL activity. The global transcriptional regulator Clp directly regulated the expression of rpfB1 and rpfB2 In conclusion, these findings reveal new roles of RpfB in antibiotic biosynthesis in L. enzymogenesIMPORTANCE Understanding the biosynthetic and regulatory mechanisms of heat-stable antifungal factor (HSAF) could improve the yield in Lysobacter enzymogenes Here, we report that RpfB1 and RpfB2 encode acyl coenzyme A (CoA) ligases. Our research shows that RpfB1 and RpfB2 affect free fatty acid metabolism via fatty acyl-CoA ligase (FCL) activity to reduce the substrate for HSAF synthesis and, thereby, block HSAF production in L. enzymogenes Furthermore, these findings reveal new roles for the fatty acyl-CoA ligases RpfB1 and RpfB2 in antibiotic biosynthesis in L. enzymogenes Importantly, the novelty of this work is the finding that RpfB2 lies outside the Rpf gene cluster and plays a key role in HSAF production, which has not been reported in other diffusible signaling factor (DSF)/Rpf-producing bacteria.

The polytopic endoplasmic reticulum (ER)-localized enzyme 3-hydroxy-3-methylglutaryl CoA reductase catalyzes a rate-limiting step in the synthesis of cholesterol and nonsterol isoprenoids. Excess sterols cause the reductase to bind to ER membrane proteins called Insig-1 and Insig-2, which are carriers for the ubiquitin ligases gp78 and Trc8. The resulting gp78/Trc8-mediated ubiquitination of reductase marks it for recognition by VCP/p97, an ATPase that mediates subsequent dislocation of reductase from ER membranes into the cytosol for proteasomal degradation. Here we report that in vitro additions of the oxysterol 25-hydroxycholesterol (25-HC), exogenous cytosol, and ATP trigger dislocation of ubiquitinated and full-length forms of reductase from membranes of permeabilized cells. In addition, the sterol-regulated reaction requires the action of Insigs, is stimulated by reagents that replace 25-HC in accelerating reductase degradation in intact cells, and is augmented by the nonsterol isoprenoid geranylgeraniol. Finally, pharmacologic inhibition of deubiquitinating enzymes markedly enhances sterol-dependent ubiquitination of reductase in membranes of permeabilized cells, leading to enhanced dislocation of the enzyme. Considered together, these results establish permeabilized cells as a viable system in which to elucidate mechanisms for postubiquitination steps in sterol-accelerated degradation of reductase.

The endoplasmic reticulum (ER)-associated degradation (ERAD) pathway in the yeast Saccharomyces cerevisiae is mediated by two membrane-bound ubiquitin ligases, Doa10 and Hrd1. These enzymes are found in distinct multiprotein complexes that allow them to recognize and target a variety of substrates for proteasomal degradation. Although multiprotein complexes containing mammalian ERAD ubiquitin ligases likely exist, they have yet to be identified and characterized in detail. Here, we identify two ER membrane proteins, SPFH2 and TMUB1, as associated proteins of mammalian gp78, a membrane-bound ubiquitin ligase that bears significant sequence homology with mammalian Hrd1 and mediates sterol-accelerated ERAD of the cholesterol biosynthetic enzyme HMG-CoA reductase. Co-immunoprecipitation studies indicate that TMUB1 bridges SPFH2 to gp78 in ER membranes. The functional significance of these interactions is revealed by the observation that RNA interference (RNAi)-mediated knockdown of SPFH2 and TMUB1 blunts both the sterol-induced ubiquitination and degradation of endogenous reductase in HEK-293 cells. These studies mark the initial steps in the characterization of the mammalian gp78 ubiquitin ligase complex, the further elucidation of which may yield important insights into mechanisms underlying gp78-mediated ERAD.

The physiological function of biotin requires biotin protein ligase activity in order to attach the coenzyme to its cognate proteins, which are enzymes involved in central metabolism. The model intracellular pathogen Francisella novicida is unusual in that it encodes two putative biotin protein ligases rather than the usual single enzyme. F. novicida BirA has a ligase domain as well as an N-terminal DNA-binding regulatory domain, similar to the prototypical BirA protein in E. coli. However, the second ligase, which we name BplA, lacks the N-terminal DNA binding motif. It has been unclear why a bacterium would encode these two disparate biotin protein ligases, since F. novicida contains only a single biotinylated protein. In vivo complementation and enzyme assays demonstrated that BirA and BplA are both functional biotin protein ligases, but BplA is a much more efficient enzyme. BirA, but not BplA, regulated transcription of the biotin synthetic operon. Expression of bplA (but not birA) increased significantly during F. novicida infection of macrophages. BplA (but not BirA) was required for bacterial replication within macrophages as well as in mice. These data demonstrate that F. novicida has evolved two distinct enzymes with specific roles; BplA possesses the major ligase activity, whereas BirA acts to regulate and thereby likely prevent wasteful synthesis of biotin. During infection BplA seems primarily employed to maximize the efficiency of biotin utilization without limiting the expression of biotin biosynthetic genes, representing a novel adaptation strategy that may also be used by other intracellular pathogens.

Phenylpropanoids play important roles in plant physiology and the enzyme 4-coumarate: coenzyme A ligase (4CL) catalyzes the formation of thioesters. Despite extensive characterization in various plants, the functions of 4CLs in the liverwort Marchantia paleacea remain unknown. Here, four 4CLs from M. paleacea were isolated and functionally analyzed. Heterologous expression in Escherichia coli indicated the presence of different enzymatic activities in the four enzymes. Mp4CL1 and Mp4CL2 were able to convert caffeic, p-coumaric, cinnamic, ferulic, dihydro-p-coumaric, and 5-hydroxyferulic acids to their corresponding CoA esters, while Mp4CL3 and Mp4CL4 catalyzed none. Mp4CL1 transcription was induced when M. paleacea thalli were treated with methyl jasmonate (MeJA). The overexpression of Mp4CL1 increased the levels of lignin in transgenic Arabidopsis. In addition, we reconstructed the flavanone biosynthetic pathway in E. coli. The pathway comprised Mp4CL1, co-expressed with chalcone synthase (CHS) from different plant species, and the efficiency of biosynthesis was optimal when both the 4CL and CHS were obtained from the same species M. paleacea.

Acyl-coenzyme A (CoA) ligases catalyze the activation of carboxylic acids via a two-step reaction of adenylation followed by thioesterification. Here, we report the discovery of a non-adenylating acyl-CoA ligase PtmA2 and the functional separation of an acyl-CoA ligase reaction. Both PtmA1 and PtmA2, two acyl-CoA ligases from the biosynthetic pathway of platensimycin and platencin, are necessary for the two steps of CoA activation. Gene inactivation of ptmA1 and ptmA2 resulted in the accumulation of free acid and adenylate intermediates, respectively. Enzymatic and structural characterization of PtmA2 confirmed its ability to only catalyze thioesterification. Structural characterization of PtmA2 revealed it binds both free acid and adenylate substrates and undergoes the established mechanism of domain alternation. Finally, site-directed mutagenesis restored both the adenylation and complete CoA activation reactions. This study challenges the currently accepted paradigm of adenylating enzymes and inspires future investigations on functionally separated acyl-CoA ligases and their ramifications in biology.

Pseudomonas species strain SBV1 can rapidly grow on medium containing β-valine as a sole nitrogen source. The tertiary amine feature of β-valine prevents direct deamination reactions catalyzed by aminotransferases, amino acid dehydrogenases, and amino acid oxidases. However, lyase- or aminomutase-mediated conversions would be possible. To identify enzymes involved in the degradation of β-valine, a PsSBV1 gene library was prepared and used to complement the β-valine growth deficiency of a closely related Pseudomonas strain. This resulted in the identification of a gene encoding β-valinyl-coenzyme A ligase (BvaA) and two genes encoding β-valinyl-CoA ammonia lyases (BvaB1 and BvaB2). The BvaA protein demonstrated high sequence identity to several known phenylacetate CoA ligases. Purified BvaA enzyme did not convert phenyl acetic acid but was able to activate β-valine in an adenosine triphosphate (ATP)- and CoA-dependent manner. The substrate range of the enzyme appears to be narrow, converting only β-valine and to a lesser extent, 3-aminobutyrate and β-alanine. Characterization of BvaB1 and BvaB2 revealed that both enzymes were able to deaminate β-valinyl-CoA to produce 3-methylcrotonyl-CoA, a common intermediate in the leucine degradation pathway. Interestingly, BvaB1 and BvaB2 demonstrated no significant sequence identity to known CoA-dependent ammonia lyases, suggesting they belong to a new family of enzymes. BLAST searches revealed that BvaB1 and BvaB2 show high sequence identity to each other and to several enoyl-CoA hydratases, a class of enzymes that catalyze a similar reaction with water instead of amine as the leaving group.

The protein degradation machinery of Mycobacterium tuberculosis includes a proteasome and a ubiquitin-like protein (Pup). Proteasome accessory factor A (PafA) attaches Pup to proteins to target them for degradation by the proteasome. Free Pup is unstable and never observed in extracts of M. tuberculosis, an observation that led us to hypothesize that PafA may need alternative sources of Pup. Here, we show that PafA can move Pup from one proteasome substrate, inositol 1-phosphate synthetase (Ino1), to two different proteins, malonyl coenzyme A (CoA)-acyl carrier protein transacylase (FabD) and lonely guy (Log). This apparent "transpupylation" reaction required a previously unrecognized depupylase activity in PafA, and, surprisingly, this depupylase activity was much more efficient than the activity of the dedicated depupylase Dop (deamidase of Pup). Thus, PafA can potentially use both newly synthesized Pup and recycled Pup to doom proteins for degradation.IMPORTANCE Unlike eukaryotes, which contain hundreds of ubiquitin ligases, Pup-containing bacteria appear to have a single ligase to pupylate dozens if not hundreds of different proteins. The observation that PafA can depupylate and transpupylate in vitro offers new insight into how protein stability is regulated in proteasome-bearing bacteria. Importantly, PafA and the dedicated depupylase Dop are each required for the full virulence of Mycobacterium tuberculosis Thus, inhibition of both enzymes may be extremely attractive for the development of therapeutics against tuberculosis.

The first step of anaerobic benzoate degradation is the formation of benzoyl-coenzyme A by benzoate-coenzyme A ligase (BCL). The anaerobic route is steered by benzoyl-CoA reductase, which promotes benzoyl-CoA breakdown, which is subsequently oxidized. In certain bacteria at low oxygen conditions, the aerobic metabolism of monoaromatic hydrocarbons occurs through the degradation Box pathway. These pathways have undergone experimental scrutiny in Alphaproteobacteria and Betaproteobacteria and have also been explored bioinformatically in representative Betaproteobacteria. However, there is a gap in our knowledge regarding the distribution of the benzoyl-CoA pathway and the evolutionary forces propelling its adaptation beyond that of representative bacteria. To address these questions, we used bioinformatic procedures to identify the BCLs and the lower pathways that transform benzoyl-CoA. These procedures included the identification of conserved motifs. As a result, we identified two motifs exclusive to BCLs, describing some of the catalytic properties of this enzyme. These motifs helped to discern BCLs from other aryl-CoA ligases effectively. The predicted BCLs and the enzymes of lower pathways were used as genomic markers for identifying aerobic, anaerobic, or hybrid catabolism, which we found widely distributed in Betaproteobacteria. Despite these enhancements, our approach failed to distinguish orthologs from a small cluster of paralogs exhibiting all the specified features to predict an ortholog. Nonetheless, the conducted phylogenetic analysis and the properties identified in the genomic context aided in formulating hypotheses about how this redundancy contributes to refining the catabolic strategy employed by these bacteria to degrade the substrates.

In plants, 4-coumarate:coenzyme A ligases (4CLs), comprising some of the adenylate-forming enzymes, are key enzymes involved in regulating lignin metabolism and the biosynthesis of flavonoids and other secondary metabolites. Although several 4CL-related proteins were shown to play roles in secondary metabolism, no comprehensive study on 4CL-related genes in the pear and other Rosaceae species has been reported. In this study, we identified 4CL-related genes in the apple, peach, yangmei, and pear genomes using DNATOOLS software and inferred their evolutionary relationships using phylogenetic analysis, collinearity analysis, conserved motif analysis, and structure analysis. A total of 149 4CL-related genes in four Rosaceous species (pear, apple, peach, and yangmei) were identified, with 30 members in the pear. We explored the functions of several 4CL and acyl-coenzyme A synthetase (ACS) genes during the development of pear fruit by quantitative real-time PCR (qRT-PCR). We found that duplication events had occurred in the 30 4CL-related genes in the pear. These duplicated 4CL-related genes are distributed unevenly across all pear chromosomes except chromosomes 4, 8, 11, and 12. The results of this study provide a basis for further investigation of both the functions and evolutionary history of 4CL-related genes.

4-Coumarate coenzyme A ligase and 4-hydroxybenzoate coenzyme A ligase from the hornwort Anthoceros agrestis expressed in E. coli were characterized on biochemical and molecular levels and showed interesting substrate specificities. Acyl-activating enzymes are associated with the biosynthesis or degradation of various metabolic products such as lipids, amino acids, sugars, and natural compounds. In this work, cDNA sequences encoding 4-coumarate coenzyme A ligase (4CL) and 4-hydroxybenzoate coenzyme A ligase (4HBCL) were amplified from the hornwort Anthoceros agrestis. The coding sequences were expressed in E. coli and purified by Ni-chelate chromatography. The CoA ligases exhibited different substrate specificities. 4CL catalyzed the activation of 4-coumaric acid, 3-coumaric acid, 2-coumaric acid, caffeic acid, isoferulic acid, ferulic acid, and cinnamic acid but lacked activities towards sinapic acid and benzoic acids. In contrast, 4HBCL preferred 4-hydroxybenzoic acid and benzoic acid, but also accepted other benzoic acid derivatives except salicylic acid and 3-aminosalicylic acid. Furthermore, 4HBCL also activated isoferulic acid, cinnamic acid, 2-coumaric acid, 3-coumaric acid, 4-coumaric acid and caffeic acid, but lacked affinity for ferulic acid and sinapic acid. These substrate specificities could be related to the phenolic compounds identified in Anthoceros agrestis.

Activation of fatty acids as acyl-adenylates by fatty acyl-AMP ligases (FAALs) in Mycobacterium tuberculosis is a variant of a classical theme that involves formation of acyl-CoA (coenzyme A) by fatty acyl-CoA ligases (FACLs). Here, we show that FAALs and FACLs possess similar structural fold and substrate specificity determinants, and the key difference is the absence of a unique insertion sequence in FACL13 structure. A systematic analysis shows a conserved hydrophobic anchorage of the insertion motif across several FAALs. Strikingly, mutagenesis of two phenylalanine residues, which are part of the anchorage, to alanine converts FAAL32 to FACL32. This insertion-based in silico analysis suggests the presence of FAAL homologues in several other non-mycobacterial genomes including eukaryotes. The work presented here establishes an elegant mechanism wherein an insertion sequence drives the functional divergence of FAALs from canonical FACLs.

Lipoic acid (LA) is a dithiol-containing cofactor that is essential for the function of alpha-keto acid dehydrogenase complexes. LA acts as a reversible acyl group acceptor and 'swinging arm' during acyl-coenzyme A formation. The cofactor is post-translationally attached to the acyl-transferase subunits of the multienzyme complexes through the action of octanoyl (lipoyl): N-octanoyl (lipoyl) transferase (LipB) or lipoic acid protein ligases (LplA). Remarkably, apicomplexan parasites possess LA biosynthesis as well as scavenging pathways and the two pathways are distributed between mitochondrion and a vestigial organelle, the apicoplast. The apicoplast-specific LipB is dispensable for parasite growth due to functional redundancy of the parasite's lipoic acid/octanoic acid ligases/transferases. In this study, we show that LplA1 plays a pivotal role during the development of the erythrocytic stages of the malaria parasite. Gene disruptions in the human malaria parasite P. falciparum consistently were unsuccessful while in the rodent malaria model parasite P. berghei the LplA1 gene locus was targeted by knock-in and knockout constructs. However, the LplA1((-)) mutant could not be cloned suggesting a critical role of LplA1 for asexual parasite growth in vitro and in vivo. These experimental genetics data suggest that lipoylation during expansion in red blood cells largely occurs through salvage from the host erythrocytes and subsequent ligation of LA to the target proteins of the malaria parasite.

Fatty acid uptake and metabolism are often dysregulated in cancer cells. Fatty acid activation is a critical step that allows these biomolecules to enter cellular metabolic pathways such as mitochondrial β-oxidation for ATP generation or the lipogenic routes that generate bioactive lipids such as the inositol phospholipids. Fatty acid activation by the addition of coenzyme A is catalysed by a family of enzymes called the acyl CoA synthetase ligases (ACSL). Furthermore, enhanced expression of particular ACSL isoforms, such as ACSL4, is a feature of some more aggressive cancers and may contribute to the oncogenic phenotype. This study focuses on ACSL3 and ACSL4, closely related structural homologues that preferentially activate palmitate and arachidonate fatty acids, respectively. In this study, immunohistochemical screening of multiple soft tissue tumour arrays revealed that ACSL3 and ACSL4 were highly, but differentially, expressed in a subset of leiomyosarcomas, fibrosarcomas and rhabdomyosarcomas, with consistent cytoplasmic and granular stainings of tumour cells. The intracellular localisations of endogenously expressed ACSL3 and ACSL4 were further investigated by detailed subcellular fractionation analyses of HT1080 fibrosarcoma and MCF-7 breast cancer cells. ACSL3 distribution closely overlapped with proteins involved in trafficking from the trans-Golgi network and endosomes. In contrast, the ACSL4 localisation pattern more closely followed that of calnexin which is an  endoplasmic reticulum resident chaperone. Confocal immunofluorescence imaging of MCF-7 cells confirmed the intracellular localisations of both enzymes. These observations reveal new information regarding the compartmentation of fatty acid metabolism in cancer cells.

Caloric restriction (CR) is an intervention that can increase maximal lifespan in organisms, but its application to humans remains challenging. A more feasible approach to achieve lifespan extension is to develop CR mimetics that target biochemical pathways affected by CR. Recent studies in the engineering and structural characterization of polyketide synthases (PKSs) have facilitated their use as biocatalysts to produce novel polyketides. Here, we show that by establishing a combinatorial biosynthetic route in Escherichia coli and exploring the substrate promiscuity of a mutant PKS from alfalfa, 413 potential anti-ageing polyketides were biosynthesized. In this approach, novel acyl-coenzyme A (CoA) precursors generated by promiscuous acid-CoA ligases were utilized by PKS to generate polyketides which were then fed to Caenorhabditis elegans to study their potential efficacy in lifespan extension. It was found that CR mimetics like resveratrol can counter the age-associated decline in mitochondrial function and increase the lifespan of C. elegans. Using the mitochondrial respiration profile of C. elegans supplemented for 8 days with 50 μM resveratrol as a blueprint, we can screen our novel polyketides for potential CR mimetics with improved potency. This study highlights the utility of synthetic enzymology in the development of novel anti-ageing therapeutics.

Nicotinamide riboside (NR), vitamin B3, is a substrate for nicotinamide adenine dinucleotide (NAD+)-consuming enzymes and is a coenzyme for hydride-transfer enzymes, including adenosine diphosphate (ADP)-ribose transferases, poly (ADP-ribose) polymerases, cADP-ribose synthases, and sirtuins, which play a central role in the aging process, neurodegenerative processes, and myopathy. Since cancer cachexia is a disease condition presenting with weight loss, skeletal muscle atrophy, and loss of adipose tissue in patients with advanced cancer, we hypothesized that NR intake could ameliorate sarcopenia. In this study, we investigated whether preemptive administration of NR ameliorated C26 adenocarcinoma-induced cancer cachexia and explored anti-cachexic mechanisms focused on the changes in muscle atrophy, cachexic inflammation, and catabolic catastrophe. Dietary intake of the NR-containing pellet diet significantly attenuated cancer cachexia in a mouse model. Starting with significant inhibition of cachexic factors, tumor necrosis factor alpha, and interleukin-6, NR significantly inhibited muscle-specific ubiquitin-proteasome ligases, such as atrogin-1, muscle RING-finger protein-1 (MuRF-1), mitofusin-2, and peroxisome proliferator-activated receptor gamma coactivator-1-alpha (PCG-1α). Significant inhibition of epididymal fat lipolysis was noted with significant inhibition of adipose triglyceride lipase (ATGL) gene. Furthermore, NR administration significantly increased the levels of crucial enzymes involved in the biosynthesis of NAD+ and nicotinamide phosphoribosyl transferase and significantly inhibited the NAD+-sensitive deacetylase sirtuin 1 (SIRT1). Preemptive intake of NR in patients vulnerable to cachexia can be a preemptive option to ameliorate cancer cachexia.

Proper functioning of the protein-folding quality control network depends on the network's ability to discern diverse structural perturbations to the native states of its protein substrates. Despite the centrality of the detection of misfolded states to cell home-ostasis, very little is known about the exact sequence and structural features that mark a protein as being misfolded. To investigate these features, we studied the requirements for the degradation of the yeast kinetochore protein Ndc10p. Mutant Ndc10p is a substrate of a protein-folding quality control pathway mediated by the E3 ubiquitin (Ub) ligase Doa10p at the endoplasmic reticulum (ER)/nuclear envelope membrane. Analysis of Ndc10p mutant derivatives, employing a reverse genetics approach, identified an autonomous quality control-associated degradation motif near the C-terminus of the protein. This motif is composed of two indispensable hydrophobic elements: a hydrophobic surface of an amphipathic helix and a loosely structured hydrophobic C-terminal tail. Site-specific point mutations expose these elements, triggering ubiquitin-mediated and HSP70 chaperone-dependent degradation of Ndc10p. These findings substantiate the ability of the ER quality control system to recognize subtle perturbation(s) in the native structure of a nuclear protein.

Geranyl pyrophosphate (GPP) and p-hydroxybenzoate (PHB) are the basic precursors involved in shikonins biosynthesis. GPP is derived from mevalonate (MVA) and/or 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway(s), depending upon the metabolite and the plant system under consideration. PHB, however, is synthesized by only phenylpropanoid (PP) pathway. GPP and PHB are central moieties to yield shikonins through the synthesis of m-geranyl-p-hydroxybenzoate (GHB). Enzyme p-hydroxybenzoate-m-geranyltransferase (PGT) catalyses the coupling of GPP and PHB to yield GHB. The present research was carried out in shikonins yielding plant arnebia [Arnebia euchroma (Royle) Johnston], wherein no molecular work has been reported so far. The objective of the work was to identify the preferred GPP synthesizing pathway for shikonins biosynthesis, and to determine the regulatory genes involved in the biosynthesis of GPP, PHB and GHB.