Noncatalytic chalcone isomerase-fold proteins in Humulus lupulus are auxiliary components in prenylated flavonoid biosynthesis.

Xanthohumol (XN) and demethylxanthohumol (DMX) are specialized prenylated chalconoids with multiple pharmaceutical applications that accumulate to high levels in the glandular trichomes of hops (Humulus lupulus L.). Although all structural enzymes in the XN pathway have been functionally identified, biochemical mechanisms underlying highly efficient production of XN have not been fully resolved. In this study, we characterized two noncatalytic chalcone isomerase (CHI)-like proteins (designated as HlCHIL1 and HlCHIL2) using engineered yeast harboring all genes required for DMX production. HlCHIL2 increased DMX production by 2.3-fold, whereas HlCHIL1 significantly decreased DMX production by 30%. We show that CHIL2 is part of an active DMX biosynthetic metabolon in hop glandular trichomes that encompasses a chalcone synthase (CHS) and a membrane-bound prenyltransferase, and that type IV CHI-fold proteins of representative land plants contain conserved function to bind with CHS and enhance its activity. Binding assays and structural docking uncover a function of HlCHIL1 to bind DMX and naringenin chalcone to stabilize the ring-open configuration of these chalconoids. This study reveals the role of two HlCHILs in DMX biosynthesis in hops, and provides insight into their evolutionary development from the ancestral fatty acid-binding CHI-fold proteins to specialized auxiliary proteins supporting flavonoid biosynthesis in plants.

Hops (n class="Species">Humulus lupulus L., Canpan>nabaceae) is a dioecious perennial vine, whose female cones are a key ingredient that provide unique flavor anpan>d aroma for brewing beer. pan> class="Chemical">Essential oils, bitter acids, and prenylchalcones account for the major three categories of specialized metabolites that are highly accumulated in the glandular trichomes (lupulin glands) of female cones, while different combinations of these compounds dictate the bittering and finishing of beer (1–3). Trace amounts of prenylated flavanones have also been detected in hops (4). Recent studies have demonstrated that hop terpenophenolics (a term for both bitter acids and prenylchalcones) exhibit diverse bioactivities with a high potential for pharmaceutical applications (5–8) (Fig. 1), with the prenylchalcones exhibiting higher bioactivity than the prenylflavanones, mainly due to the α,β- unsaturated ketone functional group in chalcones (9–11). Among these prenylchalcones, xanthohumol (XN, 3′-prenyl-6′-O- methylchalconaringenin) (Fig. 1) has received much attention due to its cancer-preventive, antiinflammatory, and antioxidant properties (3, 12–16). XN exhibits more powerful antioxidant activity than resveratrol, the well-known antioxidant found naturally in red wine (17).

Fig. 1.

The proposed biosynthetic pathway for terpenophenolics in hop glandular trichomes. The XN pathway, starting from phenylalanine, is boxed in blue. The possible R-groups in bitter acids are isobutyryl, isopropyl, and butan-2-yl groups. DD-acylphloroglucinols, di-dimethylallylated acylphloroglucinol; Iso-XN, isoxanthohumol; MEP pathway, plastid-localized methylerythritol phosphate pathway; OMT, O-methyltransferase; VPS, valerophenone synthase.

The proposed biosynthetic pathway for n class="Chemical">terpenophenolics in hop glanpan>dular trichomes. The pan> class="Chemical">XN pathway, starting from phenylalanine, is boxed in blue. The possible R-groups in bitter acids are isobutyryl, isopropyl, and butan-2-yl groups. DD-acylphloroglucinols, di-dimethylallylated acylphloroglucinol; Iso-XN, isoxanthohumol; MEP pathway, plastid-localized methylerythritol phosphate pathway; OMT, O-methyltransferase; VPS, valerophenone synthase.

To understand the molecular basis for the biosynthesis of n class="Chemical">terpenophenolics in hop trichomes, a total of more thanpan> 22,000 expressed sequepan> class="Chemical">nce tags (ESTs) from several hop trichome-specific cDNA libraries have been deposited in the TrichOME database [www.planttrichome.org (18)], and numerous large RNAseq datasets from different hop tissues or cultivars have also been made publically available. By mining the hops transcriptome data, we and others have functionally identified several key terpenophenolic biosynthetic enzymes from hop glandular trichomes (1, 18–23); these include carboxyl CoA ligase (CCL) genes and two aromatic prenyltransferase (PT) genes (HlPT1L and HlPT2) (22, 23). We have shown that HlPT2 physically interacts with HlPT1L to form an active metabolon that catalyzes the major prenylations in the β-bitter acid pathway with high efficiency: PT1L catalyzes the first prenylation step and PT2 catalyzes the subsequent two prenylation steps. We then successfully reconstructed the whole β-bitter acid pathway by coexpressing two CoA ligases (HlCCL2 and HlCCL4), the polyketide synthase valerophenone synthase (HlVPS), and the dimethylallyl diphosphate (DMAPP)-consuming PT complex in an optimized yeast system (DD104 strain, in which the endogenous farnesyl pyrophosphate synthase activity was down-regulated by site-mutation of K197G) (23).

In n class="Chemical">XN biosynthesis, pan> class="Chemical">p-coumaroyl-CoA is produced by the sequential actions of l-Phe ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), and p-coumaroyl-CoA ligase (HlCCL1 from hops) (Fig. 1). Chalcone synthase (CHS; EC 2.3.1.74) then catalyzes the condensation of p-coumaroyl-CoA with malonyl-CoA to form naringenin chalcone (NC). A trichome-specific CHS gene, CHS_H1, has been identified from hops (24). NC is then prenylated by HlPT1L, and further methylated by an O-methyltransferase (HlOMT1) to form XN (20, 23, 25). Paradoxically, previous transcriptome data had indicated that several chalcone isomerase (CHI; EC 5.5.1.6) genes represent the most abundant ESTs from hop glandular trichomes that accumulate massive amounts of chalcone (1, 21). While bona fide CHI enzymes were the first to be identified in the context of plant flavonoid metabolism, CHI-fold family proteins are more widespread in other domains of life, such as fungi and bacteria (26). The plant CHI family can be classified into four subfamilies (type I to type IV) according to their phylogenetic relationships. Types I and II CHIs are bona fide catalysts having CHI enzymatic activity. Type I CHIs widely exist in vascular plants and are responsible for the production of plant flavonoids (27–29). Type II CHI proteins appear to be legume-specific and are involved in isoflavonoid production (27, 30). Type III CHIs are widely present in land plants and green algae, while type IV CHIs are restricted to land plants. Structural analysis showed that all CHIs share a similar backbone conformation (25, 30). However, type III and type IV CHIs do not possess bona fide CHI activity, which led to the renaming of both types of CHIs as CHI-like proteins (CHIL). Recently, type III CHI folds from Arabidopsis were shown to bind fatty acids in vitro, and play a role in fatty acid metabolism in planta (25). Type III CHIs are thus divided into three fatty-acid binding protein subfamilies (FAP1, FAP2, and FAP3). A loss-of-function mutation in type IV CHI from Japanese morning glory (Ipomoea nil) led to low amounts of anthocyanin, although the underlying mechanism remains unknown (31). Previous phylogeny and sequence analyses suggest that bona fide CHIs are diverged from type IV CHIs, which evolved from the common ancestor FAP3 (type III CHIs). However, the role of type IV CHIs in flavonoid biosynthesis during land plant evolution remains to be determined (25, 32, 33).

Here, we used DD104 n class="Species">yeast strain to characterize two hop CHIL genes (pan> class="CellLine">HlCHIL1 and HlCHIL2). We demonstrate that although lacking CHI activity, both CHILs play a critical role in demethylxanthohumol (DMX)/XN production. HlCHIL2 enhanced the activities of CHS_H1 and PT1L through direct protein–protein interactions, suggesting that a functional metabolon composed of CHS, CHIL2, and PT1L, directs efficient production of prenylchalcones in hop glandular trichomes. On the other hand, CHIL1 appears to bind and stabilize the ring-open configuration of chalcones during XN biosynthesis. Our study not only reveals the functions of noncatalytic CHIs in flavonoid biosynthesis and the evolutionary trajectory of noncatalytic CHIs in land plants, but also provides valuable insight into engineering heterologous systems to produce high-valued prenylchalcones at a large scale.

Results

Characterization of Two CHIL Genes from Hop Glandular Trichomes.

We have previously demonstrated that pan> class="Chemical">XN biosynthesis occurs predominantly in the glandular trichomes of female hop flowers (1, 21–23). Searching the hop trichome-specific EST library led us to find 260 ESTs representing CHI genes, which were among the most abundant unigenes in the library. Two ORFs of CHIL genes were isolated and designated as HlCHIL1 and HlCHIL2, which encode polypeptides of 214 and 209 amino acids, respectively. HlCHIL1 shares 31% identity to FAP1 (At3g63170) from Arabidopsis thaliana (25) and 69% identity to a homologous protein from Cannabis sativa (GenBank accession no. JN679226), which is functionally uncharacterized to date (34). HlCHIL2 shares 69.7% identity to AtCHIL (At5g05270), which functions with bona fide CHI (At3g55120) to promote flavonoid production in Arabidopsis (35). A maximum-likelihood phylogenetic tree of plant CHIs was constructed to analyze the evolutionary relationship among plant CHIs: CHIL1 belonged to the type III subfamily (FAP1 clade), whereas CHIL2 falls into the type IV subfamily (Fig. 2). Sequence alignments showed that the conserved active residues in bona fide CHI sequences are missing in HlCHIL1 and HlCHIL2, suggesting both are noncatalytic CHIs (). Biochemical assays, using NC and isoliquiritigenin as substrates, at pH 6.4 and 7.5, supported this hypothesis: both CHIL1 and CHIL2 recombinant proteins did not exhibit CHI activity in vitro (). This result is consistent with the fact that prenylchalcones, rather than prenylflavanones, are predominantly accumulated in hop trichomes (21, 22). Quantitative RT-PCR analysis showed that expression of both HlCHIL1 and HlCHIL2 was trichome-specific (Fig. 2). Subcellular localization experiments indicated that CHIL1 and CHIL2 were localized to the cytoplasm, where the CCL1, CHS_H1, and OMT1 enzymes of XN biosynthesis are also located (Fig. 2). These results suggested that CHIL1 and CHIL2 are likely involved in XN biosynthesis.

Fig. 2.

Characterization of HlCHIL1 and HlCHIL2. (A) Phylogenetic analysis of CHI/CHIL proteins from plants using the maximum-likelihood method. A total of 43 CHI/CHIL proteins were obtained from 13 species representative of plant evolutionary history. To simplify the classification of CHI/CHIL proteins, four clades (types I to IV) are shown here. Bootstrap values (based on 1,000 replicates) >70% are shown for corresponding nodes. The two hops CHILs are marked with red dots, and functionally identified CHI/CHILs are marked with blue asterisks. Species abbreviations: At, Arabidopsis thaliana; Csa, Cannabis sativa; Gma, Glycine max; Hl, Humulus lupulus; Ini, Ipomoea nil; Lja, Lotus japonicus; Osa, Oryza sativa; Phy, Petunia hybrida; Ppa, Physcomitrella patens; Smo, Selaginella moellendorffii; Thy, Torenia hybrida; Vvi, Vitis vinifera; Zma, Zea mays. Protein sequences used in this analysis are listed in Dataset S1. (B) Quantitative RT-PCR analysis of two CHIL genes in different tissues of hop plants. Transcript levels are expressed relative to GADPH transcripts (n = 3). Cone bract, the glandular trichomes were removed; WAF, weeks after flowering. (C) Subcellular localization of XN-related enzymes in Arabidopsis leaf mesophyll protoplasts as revealed by laser confocal microscopy. Chloroplasts are revealed by red chlorophyll autofluorescence. (Scale bars, 5 µm.)

Characterization of n class="CellLine">HlCHIL1 anpan>d HlCHIL2. (A) Phylogenetic anpan>alysis of pan> class="Disease">CHI/CHIL proteins from plants using the maximum-likelihood method. A total of 43 CHI/CHIL proteins were obtained from 13 species representative of plant evolutionary history. To simplify the classification of CHI/CHIL proteins, four clades (types I to IV) are shown here. Bootstrap values (based on 1,000 replicates) >70% are shown for corresponding nodes. The two hops CHILs are marked with red dots, and functionally identified CHI/CHILs are marked with blue asterisks. Species abbreviations: At, Arabidopsis thaliana; Csa, Cannabis sativa; Gma, Glycine max; Hl, Humulus lupulus; Ini, Ipomoea nil; Lja, Lotus japonicus; Osa, Oryza sativa; Phy, Petunia hybrida; Ppa, Physcomitrella patens; Smo, Selaginella moellendorffii; Thy, Torenia hybrida; Vvi, Vitis vinifera; Zma, Zea mays. Protein sequences used in this analysis are listed in Dataset S1. (B) Quantitative RT-PCR analysis of two CHIL genes in different tissues of hop plants. Transcript levels are expressed relative to GADPH transcripts (n = 3). Cone bract, the glandular trichomes were removed; WAF, weeks after flowering. (C) Subcellular localization of XN-related enzymes in Arabidopsis leaf mesophyll protoplasts as revealed by laser confocal microscopy. Chloroplasts are revealed by red chlorophyll autofluorescence. (Scale bars, 5 µm.)

CHIL1 and CHIL2 Exert Contrasting Effects on DMX Production.

To explore the roles of n class="Chemical">HlCHILs in pan> class="Chemical">DMX biosynthesis, we initially reconstructed the DMX pathway by coexpressing CCL1, CHS_H1, and PT1L (PT genes used in this study were Arabidopsis codon-optimized sequences) in the DD104 yeast strain (23). On feeding with 0.5 mM p-coumarate, the engineered strain harboring CCL1/CHS_H1/PT1L was able to produce DMX, whereas the strain harboring CCL1/CHS_H1/PT2 could not produce DMX (). Moreover, addition of PT2 to the CCL1/CHS_H1/PT1L yeast strain did not further improve DMX production (). These results indicated that HlPT1L alone is responsible for the prenylation step in the DMX pathway (). Because DMX is converted to 6PN (6-prenylnaringenin) and 8PN (8-prenylnaringenin) in this yeast system, we report total prenylated chalcone/flavanone (DMX, 6PN, and 8PN). When we further introduced CHIL2 in the DD104 yeast strain harboring CCL1/CHS_H1/PT1L, the CCL1/CHS_H1/PT1L/CHIL2 strain produced approximate 2.3-fold higher total DMX (4.59 ± 0.98 μmol/L/OD; n = 3) than CCL1/CHS_H1/PT1L (2.0 ± 0.31 μmol/L/OD; n = 3). Introduction of CHIL2 also increased naringenin (N)/NC production by 1.3-fold (43.97 ± 5.8 μmol/L/OD vs. 34 ± 3.06 μmol/L/OD; n = 3) (Fig. 3). Unexpectedly, introduction of CHIL1 to the CCL1/CHS_H1/PT1L strain significantly reduced total DMX production by 30% (1.38 ± 0.32 μmol/L/OD; n = 3) (Fig. 3). We also checked the maximal production of N/NC and total DMX by adding p-coumarate to the culture (to final concentration 500 µM) every 24 h (29). The results showed that the DD104 yeast strain harboring CCL1/CHS_H1/PT1L/CHIL2 produced the N/NC and DMX at a relative flat rate (). Additionally, most of the N/NC (72%) was secreted into the culture medium, whereas 90% of the DMX remained inside the yeast cells ().

Fig. 3.

Effects of HlCHIL1 and HlCHIL2 on the production of NC and total DMX in engineered yeasts. (A) Chromatogram of selected ions of m/z 341.1384 for DMX/8PN/6PN using LC-qTOF-MS. Yeast strains harboring different gene combinations were grown in induction medium for 96 h before chemical extraction and analysis. PT1L used in this experiment was the Arabidopsis codon-optimized sequence. (B) Production of N/NC by yeast strains harboring no CHIL gene (control), CHIL1, and CHIL2. Data are means ± SD for at least three independent clones (t test, **P < 0.01). (C) Production of DMX/8PN/6PN by yeast strains harboring no CHIL gene (control), CHIL1 or CHIL2. Data are means ± SD for at least three independent clones (t test, **P < 0.01).

Type IV CHI-Folds Physically Interact with Plant CHSs to Enhance CHS Activity.

The enhan class="Chemical">ncing effect of CHIL2 on pan> class="Chemical">NC and DMX production leads to an assumption that CHIL2 might be physically interacting with other enzymes in the XN pathway. To test this hypothesis, we used the yeast two-hybrid (Y2H) technique to probe for interactions between CHIL2 and other enzymes in the XN pathway. The results showed that CHIL2 interact with CHS_H1, but not with other proteins in the XN pathway (Fig. 4). We further confirmed the interaction between CHIL2 and CHS_H1 using luciferase complementation imaging (Fig. 4) and coimmunoprecipitation (Co-IP) (Fig. 4). As HlPT1L is a membrane-bound protein, we used a split-ubiquitin MbY2H to test for possible interaction between CHIL2 and PT1L, and the results revealed that CHIL2 also physically interacted with PT1L protein (Fig. 4). Furthermore, LUC activity detected in Nicotiana benthamiana leaves confirmed the interaction between PT1L and CHIL2 (Fig. 4). We did not find any interaction between CHIL1 and CHIL2, PT1L, CCL1, CHS_H1, or OMT1 (Fig. 4 and ). We further evaluated the influence of CHIL2 on the enzymatic efficiency of CHS_H1 and PT1L through in vitro biochemical analysis. The combination of CHS_H1/CHIL2 had 1.5-fold higher NC production than CHS_H1 alone (). Microsomes prepared from yeast harboring CHS_H1/PT1L/CHIL2 displayed 1.4-fold higher DMX production than from yeast harboring only CHS_H1/PT1L (). The conversion rates of purified recombinant CHS_H1 (V value) were increased by ∼5.6-fold and 17.5-fold with CHIL2 for p-coumaroyl-CoA and malonyl-CoA, respectively, although the K value also increased by fivefold (Table 1). PT1L/CHIL2 had both a higher conversion rate and a lower K value than PT1L alone for both NC and DMAPP (Table 2).

Fig. 4.

Direct interactions between type IV CHIs and CHS. (A) Protein–protein interactions between CHIL2 and CHS_H1, CHIL1, and CCL1 in a Y2H system. The selective medium [SD-Trp-Leu-His containing 20 mM 3-amino-1,2,4-triazole (3AT)] was used for selecting for the interacting proteins. The CHIL2 gene was inserted into pGBKT7 vector (BD vector), and CHS_H1, CHIL1, and CCL1 were inserted into pGADT7 vector (AD vector). P.C., positive control; the pGBKT7-53 and pGADT7-T constructs used as positive controls were provided by the manufacturer. (B) Direct interactions between HlCHIL2 and CHS_H1 in N. benthamiana leaves. Luciferase image of N. benthamiana leaves coinfiltrated with the agrobacteria containing CHS_H1-nLuc and CHIL2-cLuc or CHIL2-nLuc/CHS_H1-cLuc combinations. (C) Reciprocal co-IP of MYC-tagged CHS_H1 and HA-tagged CHIL2 in yeast using HA- and Myc- antibodies. Total protein extracts were prepared from transgenic yeast strain harboring CHS_H1 and CHIL2, alone or together. (D) Protein–protein interactions between CHIL2 and PT1L in a split-ubiquitin Y2H system. Selection medium consisting of SD-Trp-Leu-His containing 20 mM 3-amino-1,2,4-triazole (3AT) was used for selecting the interactions between the proteins. The PT1L and PT2 constructs were used as positive controls. All PT genes were Arabidopsis codon-optimized sequences. (E) Direct interactions between HlCHIL2 and HlPT1L in N. benthamiana leaves. Luciferase image of N. benthamiana leaves coinfiltrated with agrobacteria containing PT1L-nLuc and CHIL2-cLuc or CHIL2-nLuc/PT1L-cLuc combinations. (F) Type IV CHIL–CHS interactions in P. patens (there are two type IV CHIL genes in the P. patens genome), S. moellendorffii, O. sativa (there are two type IV CHIL genes in O.sativa genome), and A. thaliana. (G) Production of N/NC by yeast strains harboring CHS alone or CHS-CHIL combination. Yeast strains harboring different gene combinations were grown in induction medium for 48 h before chemical extraction and analysis. Data are means ± SD for at least three independent clones (t test, **P < 0.01). N.C.1, negative control 1 (nLuc + cLuc); N.C.2, negative control 2 (Gene1-nLuc + cLuc); N.C.3, negative control 3 (nLuc + Gene2-cLuc); N.C.4, negative control 4 (Gene2-nLuc + cLuc); N.C.5, negative control 5 (nLuc + Gene1-cLuc).

Table 1.

Kinetic parameters for CHS_H1 and CHS_H1/CHIL2 complex

Complex	Substrate	Km (μM)	Kcat (s−1 × 10−3)	Kcat/Km (M−1s−1)
CHS_H1	p-Coumaroyl-CoA*	4.82 ± 0.99†	0.21 ± 0.02	44
	Malonyl-CoA‡	10.06 ± 2.99	0.24 ± 0.03	24
CHS_H1/CHIL2	p-Coumaroyl-CoA*	26.57 ± 3.78	1.2 ± 0.12	45
	Malonyl-CoA†	59.74 ± 8.49	4.2 ± 0.3	70

Two micrograms of purified recombinant CHS_H1 were used in each CHS_H1 assay; and 2 μg of purified recombinant CHS_H1 and 1.1 μg of purified recombinant CHIL2 (CHS_H1:CHIL2 = 1:1) were used in each CHS_H1/Hl CHIL2 assay.

A fixed concentration of 150 μM malonyl-CoA was used as substrate.

The data are presented as means ± SD (n = 3).

A fixed concentration of 40 μM p-Coumaroyl-CoA was used as substrate.

Table 2.

Kinetic parameters for PT1L and PT1L/CHIL2 complex

Complex	Substrate	Km (μM)	Vmax (μmol/min/g)	Vmax/Km
PT1L	NC*	5.74 ± 0.68†	0.1 ± 0.01	0.0174
	DMAPP‡	75.12 ± 15.7	0.13 ± 0.02	0.0017
PT1L/CHIL2	NC*	5.01 ± 0.43	0.12 ± 0.01	0.024
	DMAPP‡	62.78 ± 6.76	0.23 ± 0.01	0.0037

Total membrane-bound proteins (10 μg in each assay) were prepared from the yeast harboring PT1L or PT1L/CHIL2.

A fixed concentration of 50 μM DMAPP was used as substrate.

The data are presented as means ± SD (n = 3).

A fixed concentration of 200 μM NC was used as substrate.

Direct interactions between n class="Disease">type IV CHIs anpan>d pan> class="Gene">CHS. (A) Protein–protein interactions between CHIL2 and CHS_H1, CHIL1, and CCL1 in a Y2H system. The selective medium [SD-Trp-Leu-His containing 20 mM 3-amino-1,2,4-triazole (3AT)] was used for selecting for the interacting proteins. The CHIL2 gene was inserted into pGBKT7 vector (BD vector), and CHS_H1, CHIL1, and CCL1 were inserted into pGADT7 vector (AD vector). P.C., positive control; the pGBKT7-53 and pGADT7-T constructs used as positive controls were provided by the manufacturer. (B) Direct interactions between HlCHIL2 and CHS_H1 in N. benthamiana leaves. Luciferase image of N. benthamiana leaves coinfiltrated with the agrobacteria containing CHS_H1-nLuc and CHIL2-cLuc or CHIL2-nLuc/CHS_H1-cLuc combinations. (C) Reciprocal co-IP of MYC-tagged CHS_H1 and HA-tagged CHIL2 in yeast using HA- and Myc- antibodies. Total protein extracts were prepared from transgenic yeast strain harboring CHS_H1 and CHIL2, alone or together. (D) Protein–protein interactions between CHIL2 and PT1L in a split-ubiquitin Y2H system. Selection medium consisting of SD-Trp-Leu-His containing 20 mM 3-amino-1,2,4-triazole (3AT) was used for selecting the interactions between the proteins. The PT1L and PT2 constructs were used as positive controls. All PT genes were Arabidopsis codon-optimized sequences. (E) Direct interactions between HlCHIL2 and HlPT1L in N. benthamiana leaves. Luciferase image of N. benthamiana leaves coinfiltrated with agrobacteria containing PT1L-nLuc and CHIL2-cLuc or CHIL2-nLuc/PT1L-cLuc combinations. (F) Type IV CHIL–CHS interactions in P. patens (there are two type IV CHIL genes in the P. patens genome), S. moellendorffii, O. sativa (there are two type IV CHIL genes in O.sativa genome), and A. thaliana. (G) Production of N/NC by yeast strains harboring CHS alone or CHS-CHIL combination. Yeast strains harboring different gene combinations were grown in induction medium for 48 h before chemical extraction and analysis. Data are means ± SD for at least three independent clones (t test, **P < 0.01). N.C.1, negative control 1 (nLuc + cLuc); N.C.2, negative control 2 (Gene1-nLuc + cLuc); N.C.3, negative control 3 (nLuc + Gene2-cLuc); N.C.4, negative control 4 (Gene2-nLuc + cLuc); N.C.5, negative control 5 (nLuc + Gene1-cLuc).

Kinetic parameters for n class="Disease">CHS_H1 and CHS_H1/CHIL2 complex

Total membrane-bound proteins (10 μg in each assay) were prepared from the n class="Species">yeast harboring PT1L or PT1L/CHIL2.

Given that stress-protective n class="Chemical">flavonoids, pan> class="Disease">type IV CHIs, and CHSs are widely distributed in land plants (25, 36, 37), we assume that type IV CHIs share the conserved function of increasing flavonoid production by binding and increasing the enzymatic efficiency of plant CHSs. To test this hypothesis, we probed the possible type IV CHI–CHS interactions in four plant species [Physcomitrella patens (Bryophytes), Selaginella moellendorffii (Lycophytes), Oryza sativa (monocots in Euphyllophytes), and Arabidopsis thaliana (dicots in Euphyllophytes)], which are located at the major nodes during plant evolution. The physical interactions were detected for all six tested type IV CHI–CHS pairs (there are two type IV CHIL genes in the P. patens and O. sativa genomes) (Fig. 4). All tested type IV CHI proteins enhanced the CHS activity by 1.5- to 3-fold, as indicated by N/NC production in the engineered yeast system (Fig. 4).

To test the function of CHIL2 in planpan>ta, we also genepan> class="Species">rated transgenic Arabidopsis plants overexpressing CHIL2, PT1L, and CHIL2-PT1L in the tt5-1 background [TT5 (At3g55120), lacking a functional CHI protein; Ler ecotype], which accumulates approximate 60-fold more NC compared with Ler wild-type (). LC-MS analysis showed that the NC content was increased in CHIL2 overexpression lines by about 1.6-fold (), indicating that CHIL2 increased CHS activity in planta. This result drove us to test whether HlCHIL2 physically interacted with the CHS from A. thaliana (At5g13930) and other above-mentioned plant CHSs. We found direct interactions between HlCHIL2 and all tested CHS proteins from P. patens, S. moellendorffii, O. sativa, and A. thaliana (). However, we were unable to detect DMX or prenylated flavonoids in PT1L or CHIL2-PT1L overexpressing Arabidopsis lines, despite the detection of PT1L mRNA in these transgenic plants ().

CHIL1 Binds to DMX and NC to Stabilize the Chalconoids’ Ring-Opened Structure.

Because CHIL1 significantly reduced n class="Chemical">DMX production in engineered pan> class="Species">yeast cells (Fig. 3), although no physical interaction between CHIL1 and other XN-related proteins could be detected (), it is possible that CHIL1 protein binds to DMX, similar to the way that AtFAPs bind with fatty acids in A. thaliana (25), and that the bound DMX cannot be released from this metabolite–protein complex by the extraction method used here. To test this hypothesis, we incubated purified His-tagged CHIL1 with the chemical extract (DMX/6PN/8PN) prepared from yeast harboring CCL1/CHS_H1/CHI2/PT1L. Most of the DMX was eventually converted to 6PN and 8PN after incubation with His-tagged proteins for 8-h at 4 °C (see for the experimental workflow). However, CHIL1 protein enriched DMX levels by 14.2-fold [the ratio of (DMX)/(total DMX) in CHIL1-bound chemical mixture was 19.7%; n = 3, whereas the ratio of (DMX)/(total DMX) in solution was only 1.39%] (Fig. 5). In contrast, no CHIL2-associated DMX or 8PN were detected under the same incubation conditions, and the MagneHis Ni-Beads showed weak nonspecific binding activity to 6PN only (Fig. 5 ). Differential scanning fluorimetry assays (38) using purified HlCHIL1 protein and standard chemicals further showed that CHIL1 protein binds DMX and NC with K (dissociation constant) values of 25.08 μM (n = 3) and 43.56 μM (n = 3), respectively (Fig. 5 ), whereas no such binding could be detected with N, 6PN, 8PN, or XN. We therefore conclude that the 6PN and 8PN detected in CHIL1-bound compounds are most likely converted from DMX during sample preparation given the known spontaneous isomerization of DMX in vitro. CHIL1 protein showed weaker binding activity to lauric acid (C12:0) and myristic acid (C14:0), than to longer-chain fatty acids (C16:0, C18:0, and C18:3) (). It should be noted that both DMX (ΔTmmax = 10.23 °C, n = 3) and NC (ΔTmmax = 6.40 °C, n = 3) show a stronger stabilizing effect on CHIL1 protein than do lauric acid (ΔTmmax = 5.0 °C, n = 2) or myristic acid (ΔTmmax = 3.55 °C, n = 2). However, CHIL1 preferentially binds lauric acid (K = 11.81 µM, n = 2) and myristic acid (K = 3.45 µM, n = 2) compared with DMX (K = 25.08 µM, n = 3) and NC (K = 43.56 µM, n = 3) (). These results demonstrate that CHIL1 protein not only binds to DMX and NC, but also stabilizes the ring-open configuration of these compounds.

Fig. 5.

Identification of DMX binding of CHIL1. (A) Analysis of ligands associated with purified His-tagged HlCHIL1 protein (linked to MagneHis Ni-Particles) separated and detected by LC-QQQ-MS/MS (MRM condition for DMX/6PN/8PN: 341.0 → 165.0). MagneHis Ni-Beads alone were used as negative controls in these assays. All samples, including the chemical extract only, were incubated at 4 °C for 8 h before LC-MS analysis (for details, see ). (B) Relative quantification of DMX in protein-bound and supernatant fractions of different CHIL–small-molecule interaction assays. N.D., not determined due to low content of DMX in the samples (t test, **P < 0.01). (C) Representative dose–response curves for HlCHI1 generated using thermal shift response to increasing DMX concentration. ΔTmmax and K values were calculated as 10.23 °C and 25.08 µM, respectively (n = 3). (D) Representative dose–response curve for HlCHIL1 generated using thermal shift response to increasing NC concentration. ΔTmmax and K values were calculated as 6.40 °C and 43.56 µM, respectively (n = 3). (E) Homology model of HlCHIL1 (marine) with docked DMX ligand (yellow sticks) overlaid with AtFAP1 (PDB ID code 4DOO; magenta) bound to lauric acid (green sticks). The homology model was built using AtFAP1 as a template. (F) Zoomed-in view of HlCHIL1 with docked DMX, showing relevant protein side chains as gray sticks. Selective hydrogen atoms are shown in purple, H-bonding interactions as dotted black lines, and potential π-stacking or cation–π stabilization interactions as solid black arrows. Val139 and Ala143 are highlighted in red due to their important role in DMX binding.

Identification of n class="Chemical">DMX binding of CHIL1. (A) Analysis of liganpan>ds associated with purified pan> class="Chemical">His-tagged HlCHIL1 protein (linked to MagneHis Ni-Particles) separated and detected by LC-QQQ-MS/MS (MRM condition for DMX/6PN/8PN: 341.0 → 165.0). MagneHis Ni-Beads alone were used as negative controls in these assays. All samples, including the chemical extract only, were incubated at 4 °C for 8 h before LC-MS analysis (for details, see ). (B) Relative quantification of DMX in protein-bound and supernatant fractions of different CHIL–small-molecule interaction assays. N.D., not determined due to low content of DMX in the samples (t test, **P < 0.01). (C) Representative dose–response curves for HlCHI1 generated using thermal shift response to increasing DMX concentration. ΔTmmax and K values were calculated as 10.23 °C and 25.08 µM, respectively (n = 3). (D) Representative dose–response curve for HlCHIL1 generated using thermal shift response to increasing NC concentration. ΔTmmax and K values were calculated as 6.40 °C and 43.56 µM, respectively (n = 3). (E) Homology model of HlCHIL1 (marine) with docked DMX ligand (yellow sticks) overlaid with AtFAP1 (PDB ID code 4DOO; magenta) bound to lauric acid (green sticks). The homology model was built using AtFAP1 as a template. (F) Zoomed-in view of HlCHIL1 with docked DMX, showing relevant protein side chains as gray sticks. Selective hydrogen atoms are shown in purple, H-bonding interactions as dotted black lines, and potential π-stacking or cation–π stabilization interactions as solid black arrows. Val139 and Ala143 are highlighted in red due to their important role in DMX binding.

To further investigate the structn class="Chemical">ural basis for pan> class="Chemical">DMX binding, a homology model was generated for HlCHIL1 using AtFAP1 (PDB ID code 4DOO) as a template (25). Subsequent docking simulations for DMX revealed a similar binding pocket compared with lauric acid (C12:0) in the AtFAP1 model (Fig. 5). The best-fitting DMX pose positioned the molecule such that the 4′- and 6′-hydroxyl groups of ring A formed hydrogen bonds with the backbone carbonyl of V139 and the side chain of Y186, respectively, whereas M155 and M196 could be found within reasonable hydrogen bonding distance of the 2′-hydroxyl of ring A and the 4′-hydroxyl of ring B, respectively (Fig. 5). Additional binding interactions are established via π-stacking of F115 with ring A and potential cation–π stabilization via R42 with ring B (Fig. 5). This binding orientation places the Cβ atom of DMX 5.0 Å away from the 2′-hydroxyl, rationalizing both the ability of HlCHIL1 to stabilize the ring-open conformation and the absence of any CHI activity.

CHIL1 and Its Homologs Are Functionally Diverged from FAP1 Proteins.

To trace the evolutionary n class="Chemical">history of pan> class="CellLine">HlCHIL1, we identified HlCHIL1 homologs in diverse plant species by BLAST (National Center for Biotechnology Information and Phytozome database, phytozome.jgi.doe.gov/) using 50% identity as cut-off. The comprehensive phylogenetic analysis of the FAP1 subfamily shows that 33 HlCHIL1 homologs form a distinct branch which separated from FAP1 proteins after the emergence of the angiosperms, as no close homologs of HlCHIL1 can be found in ancient species, such as gymnosperms, pteridophytes, and bryophytes (). Interestingly, most AtFAP1 homologs, ranging from bryophytes to flowering plants, show a clear N-terminal plastidial signal peptide () (predicted by TargetP software, www.cbs.dtu.dk/services/TargetP/), suggesting that AtFAP1 homologs shared a conserved biochemical function involving fatty-acid metabolism, as demonstrated previously (25). However, all of the closest homologs of HlCHIL1 do not have a clear N-terminal signal peptide. This, together with the biochemical characterization of HlCHIL1 (binding with DMX for structural stabilization), suggests a divergent function of HlCHIL1 homologs from the conserved fatty-acid binding activity of FAP1s.

When comparing the amino acids, which are critical for binding n class="Chemical">fatty acids in Atpan> class="Gene">FAP1 or DMX in HlCHIL1, most amino acids in the fatty acids binding cleft of FAP1 proteins, including the Arg-Tyr pair that tethers the carboxylate group, are conserved in CHIL1 homologs (Fig. 6). It is noteworthy that Phe203 and Phe207 in AtFAP1 are substituted with Val139 and Ala143 in HlCHIL1, which results in an enlarged binding pocket for accommodating the aromatic ring A of DMX. To test the importance of these two amino acids for DMX binding, we generated V139F and A143F CHIL1 single mutants and V139F/A143F double mutant. Binding assays clearly showed that single CHIL1 mutant (CHIL1V139F and CHIL1A143F) had similar binding properties toward DMX (ΔTmmax = 10.15 °C and 10.87 °C for CHIL1V139F and CHIL1A143F, respectively; n = 3) and NC (ΔTmmax and K are 5.50 °C and 7.34 °C for CHIL1V139F and CHIL1A143F, respectively; n = 2), compared with those of wild-type CHIL1 (ΔTmmax = 10.23 °C for DMX and 6.40 °C for NC; n = 3) (Fig. 6). However, both DMX (ΔTmmax = 5.26 °C, n = 3) and NC (ΔTmmax = 4.70 °C, n = 2) show much weaker stabilizing effect on CHIL1V139F/A143F compared with CHIL1 (Fig. 6), suggesting that V139 and A143 of CHIL1 have a joint effect on DMX/NC binding to CHIL1. It is noteworthy that CHIL1V139F shows higher affinity to DMX (K = 10.05 µM) and NC (K = 7.61 µM) than CHIL1 (Fig. 6). Altogether, we tentatively designate these HlCHIL1 homologs as polyketide (a term covering flavonoid) binding proteins (PBP).

Fig. 6.

Sequence analysis of HlCHIL1 homologs and chemical binding identification of CHIL1 mutants. (A) Thirty-three HlCHIL1-close proteins (>50% identity) were added for the phylogenetic analysis. The GenBank number for each HlCHIL1-close protein is provided in parenthesis, next to the Latin name of the plant species. Bootstrap values (based on 1,000 replicates) >70% are shown for corresponding nodes. HlCHIL1 and CsaCHIL are marked in red dots. Twenty residues lining the DMX binding site, based on the HlCHIL1 modeling information, are denoted by red asterisks. Twenty-six residues lining the fatty-acid binding site, based on the AtFAP1 structure information, are denoted by blue asterisks. The Arg-Tyr pair (R103 and Y116 in AtFAP1) that tethers the carboxylate group in FAP1s and the Val139 and Ala143 that are critical for accommodation of DMX in HlCHIL1 are boxed in yellow and red, respectively. For species abbreviations, see the legend to Fig. 1. Protein sequences used in this analysis are listed in Dataset S2. (B) Representative dose–response curve for CHIL1, CHIL1V139F, CHIL1A143F, and CHIL1V139F/A143F generated using thermal shift response to increasing DMX or NC concentration.

Sequen class="Chemical">nce anpan>alysis of pan> class="CellLine">HlCHIL1 homologs and chemical binding identification of CHIL1 mutants. (A) Thirty-three HlCHIL1-close proteins (>50% identity) were added for the phylogenetic analysis. The GenBank number for each HlCHIL1-close protein is provided in parenthesis, next to the Latin name of the plant species. Bootstrap values (based on 1,000 replicates) >70% are shown for corresponding nodes. HlCHIL1 and CsaCHIL are marked in red dots. Twenty residues lining the DMX binding site, based on the HlCHIL1 modeling information, are denoted by red asterisks. Twenty-six residues lining the fatty-acid binding site, based on the AtFAP1 structure information, are denoted by blue asterisks. The Arg-Tyr pair (R103 and Y116 in AtFAP1) that tethers the carboxylate group in FAP1s and the Val139 and Ala143 that are critical for accommodation of DMX in HlCHIL1 are boxed in yellow and red, respectively. For species abbreviations, see the legend to Fig. 1. Protein sequences used in this analysis are listed in Dataset S2. (B) Representative dose–response curve for CHIL1, CHIL1V139F, CHIL1A143F, and CHIL1V139F/A143F generated using thermal shift response to increasing DMX or NC concentration.

Discussion

It is well-known that n class="Chemical">flavonoids are ubiquitously present in pan> class="Species">land plants, where they may have played a key role in land colonization during plant evolution (39). CHS, the first enzyme of the flavonoid biosynthetic pathway, has been intensively studied at the biochemical and molecular levels during past three decades. Recently, a proteolytic regulator (At1g23390, encoding a Kelch domain-containing F-box protein) controlling CHS stability was functionally identified from A. thaliana (40). Here, we have identified another regulator of plant CHS, the type IV CHI-fold proteins, which physically interact with CHS to increase its activity, thereby enhancing flavonoid production. This characterization explains the previous observations that loss-of-function of type IV CHI from Japanese morning glory (I. nil) and Arabidopsis leads to decreased amounts of flavonoids (31, 35). Presumably through protein–protein interactions, the type IV CHI remodels the active-site cavity of CHS for synthesizing NC. The crystal structure of the HlCHIL2/HlCHS_H1 will definitely elucidate the underlying molecular mechanism.

The n class="Gene">CHS-enhapan> class="Chemical">ncing property of the type IV CHI-fold proteins is conserved from mosses to flowering plants (Fig. 4). Although phylogenetic analysis suggested that the type IV CHI-fold proteins were probably the ancestor of plant functional CHIs, the function of type IV CHI-fold proteins remained unknown (25, 32, 33). Based on our data, we propose that binding with CHS is a conserved function for all plant CHI-fold proteins, type I and type II CHI-fold proteins inherited this trait from type IV CHI-fold proteins besides the gained CHI activity during enzyme evolution. The physical interaction between bona fide CHI and CHS has been demonstrated in plants (41). Given the importance of flavonoids, enhancing CHS activity likely enhances plant fitness. Phylogenetic analysis also shows that type IV CHIs are likely derived from FAP3 proteins, which are involved in fatty acid metabolism in plants. Therefore, it will be worth testing whether plant FAP3 proteins physically interact with enzymes of fatty acid biosynthesis, especially the ketoacyl-[acyl-carrier-protein] synthase III, which is the progenitor of plant CHSs (42).

HlCHIL2 physically interacts with PT1L to form an active metabolon for n class="Chemical">DMX production in hop glanpan>dular trichomes. PT1L, as a membranpan>e-bound protein, therefore fupan> class="Chemical">nctions to anchor the complex with CHIL2, which also interacts with CHS_H1, to the membrane. In Arabidopsis, a type IV CHI-fold protein (AtCHIL, At5g05270) was demonstrated to bind with bona fide CHI protein (type I CHI, At3g55120) for flavonoid production (35). These results suggest that type IV CHI-fold proteins may bind with enzymes downstream of CHS in a species-dependent manner; this needs to be experimentally validated.

Meanwhile, CHIL1 binds and stabilizes the ring-open conformation of n class="Chemical">NC anpan>d pan> class="Chemical">DMX, the intermediates of XN biosynthesis. The DMX-binding ability of CHIL1 and the DMX metabolon presented here are consistent with the chemical phenomenon observed in hop glandular trichomes: high accumulation of DMX and XN, with almost no N/NC being detected (21, 22). The structural information and sequence analysis has led to the discovery of a protein clade (PBP) with unprecedented functions: HlCHIL1 homologs probably bind with the products of other type III polyketide synthases, at least those that contain an aromatic ring derived from the condensation of three molecules of malonyl CoA (Fig. 5). CsaCHIL (from C. sativa, Cannabaceae) is the closest homolog of HlCHIL1 (69% identity), and is specifically expressed in hemp glandular trichomes (34). Two polyketides, Δ9-tetrahydrocannabinolic acid (THCA) and cannabidiolic acid (CBDA), are synthesized and stored at high levels in hemp glandular trichomes. THCA and CBDA, like DMX, are unstable and will be nonenzymatically converted to the decarboxylated forms, Δ9-tetrahydrocannabinol and cannabidiol, respectively. We therefore hypothesize that CsaCHIL probably binds with THCA and/or CBDA to stabilize these compounds in hemp glandular trichomes. Furthermore, the endogenous content of free polyketides will be accordingly decreased due to increased protein-bound forms, which can alleviate the feedback inhibition of polyketide biosynthesis to result in high production of polyketides in a specific plant cell/tissue/organ.

In con class="Chemical">nclusion, we have fupan> class="Chemical">nctionally identified two CHIL genes, HlCHIL1 and HlCHIL2, which play a significant role in DMX biosynthesis in hop glandular trichomes: CHIL2 enhances the catalytic efficiency of CHS_H1 and PT1L through protein–protein interactions, whereas CHIL1 stabilizes the ring-open conformation of DMX. Moreover, the engineered yeast generated in this study produced up to 4.59 μmol/L/OD of DMX under flask-shake conditions, which will be a starting point for producing these valuable prenylchalcones at large scale using microbial cell factories.

Materials and Methods

Plant Materials, RNA Analysis, and Chemicals.

The growth of n class="Species">H. lupulus cv. Nugget, EST sequepan> class="Chemical">nce analysis, RNA isolation, and cDNA preparation from hop tissues were performed as described previously (1). A. thaliana (Col-0 ecotype) and N. benthamiana plants for transient transformation were maintained in a greenhouse under 16-h light/8-h dark, 22 °C conditions. All available commercial chemicals used in this study were purchased from Sigma-Aldrich, except for NC, which was purchased from ChromaDex. PIVP (phlorisovalerophenone) was synthesized as described previously (43). The purity and concentration of these chemicals were determined using LC-MS.

Hop CHI Gene Isolation.

To obtain the full-length sequences of HlCHIL1 and HlCHIL2 from hop trichomes, the ORFs of HlCHIL1 and HlCHIL2 obtained by RT-PCR were cloned into pEASY-Blunt vector (TransGen Biotech) and verified by sequencing of at least five independent clones (see for primer information).

Quantitative RT-PCR Analysis.

Real-time PCR analyses were performed using Ultra SYBR Mixture (CWBio) on a CFX96 Real-Time PCR Detection System (Bio-Rad) following the manufacturer’s instructions. Ct values were calculated using the Bio-Rad real-time analysis software. Compan class="Species">rative Ct method was used for relative gene expression anpan>alysis by normalizing to the referepan> class="Chemical">nce gene [glyceraldehyde-3-P dehydrogenase (GAPDH) from hop and Actin2 (At3g18780) from Arabidopsis]. Every PCR was repeated with three independent biological replicates. Primers are listed in . All raw real-time PCR analysis data are showed in Dataset S3.

Subcellular Localization.

The ORF of the hGFP gene was fused to the C-terminal of the CCL1, CHS_H1, CHIL2, CHIL1, and OMT1 ORFs, under control of the CaMV 35S promoter (pJIT163-hGFP vector). Arabidopsis leaf protoplast preparation, transformation, and image assay using laser scanning confocal microscopy were performed as described previously (44). Briefly, mesophyll protoplasts freshly isolated from rosette leaves of 4-wk-old Arabidopsis (Col-0 ecotype). The fresh prepared protoplasts are transfected with 10 μg plasmid using a PEG-calcium–mediated transfection method. Living cellular image of GFP fusion proteins are observed under Axio Imager Z2 fluorescence microscopy (Zeiss). Localization was determined by surveying more than 50 protoplasts. Primers are listed in .

Reconstitution of DMX Pathway in Yeast and Product Analysis.

The following vectors were used: pESC-n class="Chemical">Leu for pan> class="Gene">CCL1 and CHIL2, pESC-His for CHS_H1 and CHIL1, and pESC-Ura for PT1L and PT2. The different combinations of three constructs (either empty vector or the construct with inserts) were cotransformed into yeast strain DD104 using a high-efficiency lithium acetate transformation protocol (45). The resulting positive yeast clones were cultured in SD dropout medium (-Leu, -His, -Ura), and then harvested and induced with 2% galactose dropout medium (-Leu, -His, -Ura) for 4 d, then fed with 500 μM p-coumarate. Chemical extraction and LC-MS analysis were performed as previously described (22, 23). PIVP (1.5 μM) was added as internal standard. The contents of N/NC and DMX/6PN/8PN in the samples were quantified based on the concentration of the internal standard after normalization.

In Vitro Enzyme Assays.

n class="Species">Yeast microsome preparations were made as described previously (46). The soluble and microsomal proteins were used directly for enzymatic assays after quantification of protein by Bradford assay. The chalcone synthase assay reaction contained 50 mM Tris⋅HCl (pH 7), 0.1 mg/mL BSA, 30 μM malonyl-CoA, 10 μM p-coumaroyl-CoA, and 40 μg total soluble proteins. The reaction mixtures were incubated in a total volume of 200 μL at 30 °C for 60 min. The PT assays were performed as previously described (23). Control reactions were performed with boiled protein. For the determination of apparent K values for different substrates, 40 μM p-coumaroyl-CoA was added in the standard CHS_H1 reactions (200 μL volume) containing 2 μg purified recombinant CHS_H1, with or without 1.1 μg of purified recombinant CHIL2 (CHS_H1:CHIL2 = 1:1), and malonyl-CoA at a series of concentrations, or 150 μM malonyl-CoA was added in the standard CHS_H1 reactions containing 2 μg purified recombinant CHS_H1, with or without 1.1 μg of purified recombinant CHIL2 (CHS_H1:CHIL2 = 1:1), and p-coumaroyl-CoA at a series of concentrations. DMAPP (50 μM) was added in the standard PT reactions (200 μL) containing 10-μg membrane proteins and NC at a series of concentrations, or 200 μM NC was added in the standard PT reactions containing 10-μg membrane proteins and DMAPP at a series of concentrations. The enzymatic products were extracted and quantified as previously described (23). The apparent K data were calculated using Hanes plots (Hyper32, v1.0.0).

CHI assays were carried out as described previously (27, 47), using n class="Chemical">NC anpan>d pan> class="Chemical">isoliquiritigenin as substrates, at pH 6.4 and 7.5. Two microliters of each enzymatic reaction, after filtered through a 0.22-μm syringe filter, was loaded on a 1290 Infinity LC pump coupled to a 6495 triple quadrupole mass spectrometer equipped with a dual electrospray ion source operated in positive mode (LC-QQQ-MS/MS; Agilent). Chromatographic separation was performed on a ZORBAX Extend C18 column (50 mm × 2.1 mm ID, 1.8 µm; Agilent). Gradient condition of the mobile phase (solution A is water and solution B is methanol; flow rate is 0.35 mL/min) was set as follows: 0–6 min, a linear gradient of from 40% of B to 60% of B; 6.0–6.5 min, a linear gradient of from 60% of B to 98% of B; 6.5–7.5 min, 98% of B; then the system was equilibrated using the initial condition (40% of B) for 5 min before the next sample injection.

Yeast Two-Hybrid Assays.

n class="Species">Yeast two-hybrid assays were based on the Matchmaker pan> class="Gene">GAL4 two-hybrid system (Clontech). CHIL2 was inserted into pGBKT7 vector, and CCL1, CHS_H1, CHIL1 were inserted into pGADT7 vector. Constructs were cotransformed into the yeast strain AH109. The presence of the transgenes was confirmed by growth on SD/-Leu/-Trp plates. To assess protein interactions, the transformed yeast was suspended in liquid SD/-Leu/-Trp medium and cultured to OD600 = 1.0. Ten microliters of suspended yeast was spread on SD/-His/-Leu/-Trp medium with 20 mM 3-AT. Interactions were observed after 3 d of incubation at 30 °C. The split-ubiquitin MbYTH assays were performed with a DUAL membrane kit following the manufacturer’s instructions (Dualsystems Biotech). PT1L was cloned into the bait vector pTMBV4, and CHIL1, CHIL2, CHS_H1, and CCL1 were cloned into the pDL2-Nx vector. Constructs were cotransformed into the yeast strain DSY-1. The experiments were repeated for at least three times. The primers used in this study are listed in .

Luciferase Complementation Imaging Assay.

CHIL genes (AtCHIL, OsCHIL1, OsCHIL2, SmCHIL, PpCHILa, and PpCHILb), PT1L, and n class="Gene">CHS genes (pan> class="Disease">CHS_H1, AtCHS, OsCHS, SmCHS and PpCHS) were fused to the N- or C-terminal part of LUC to generate the corresponding nLUC and cLUC constructs. Agrobacterium tumefaciens GV3101 containing the corresponding n-LUC and c-LUC fusion constructs were simply mixed with P19, then infiltrated into leaves of N. benthamiana for transient expression. LUC activity was detected after 48–72 h. The luciferase complementation imaging assays were performed as previously described (48), using three biological replicates, by infiltration of leaves on separate plants. Primers used for the vector construction are shown in .

Co-IP Assays.

N-terminal HA-tagged CHIL2 and Myc-tagged n class="Gene">CHS_H1 were cloned into the pESC-pan> class="Chemical">Trp and pESC-HIS vectors respectively. HA-pESC-Trp and Myc-pESC-HIS plasmids which expressed HA tag or Myc tag alone were constructed as negative controls using a PCR-mediated method (for primer information, see ). The transformed yeast cells (DSY-1 strain) were harvested after galactose induction until the OD600 value reached 1.0, and disrupted with glass beads with the Mini-Bead-Beater (Biospec Products) in an extraction buffer (2 mL/g cells) containing 50 mM Tris⋅HCl (pH 7.5), 1 mM EDTA, 1%Triton X-100, 1 mM phenylmethylsulfonyl fluoride, 1 μg/mL leupeptin, and 1 μg/mL pepstatin. Glass beads and cell debris were removed by centrifugation at 12,000 × g for 10 min. The supernatant was incubated with the first antibody (anti-Myc) at 4 °C overnight with shaking, followed by incubation with Protein G Agarose beads for an additional 2 h. The beads were washed twice with extraction buffer supplemented with 300 mM NaCl and then the proteins were released by incubation for 10 min in SDS sample buffer at 98 °C and analyzed with immunoblotting using an anti-HA antibody.

Chemical Measurements in the Transgenic Arabidopsis Lines.

Empty pGreen binary vector (control) (49), CHIL2, PT1L, and CHIL2-PT1L were introduced into tt5-1 mutant n class="Species">Arabidopsis using the floral-dip method (50). Two-week-old seedlings, grown on 1/2 MS plates, were ground into a fine powder in liquid pan> class="Chemical">nitrogen, of which 200-mg fresh weight equivalents was extracted with 1.6 mL of ethyl acetate. After centrifugation for 10 min, the supernatant was evaporated to dryness and the residue dissolved in methanol and filtered through a 0.22-μm syringe filter. N/NC measurements were conducted using above-mentioned LC-QQQ-MS/MS program.

Protein–Metabolite Interaction Assay.

CHIL genes were subcloned into the pEasy-Blunt E2 Expression vector (Transgen). All constructs were transformed into n class="Species">Escherichia coli BL21(+) cells for prokaryotic expression, anpan>d the resulting pan> class="Chemical">His-tagged fusion proteins were purified using Ni-NTA affinity chromatography. Quantification and evaluation of the relative purity of the recombinant proteins was performed using SDS/PAGE with BSA as a standard. The in vitro reaction buffer contained 50 mM Tris⋅HCl, pH 7.5, 20% methanol, 8 μL yeast extracts (final concentration around 200 μM N/NC, 5 μM DMX), and 30 μg purified protein in a final volume of 500 μL. After incubation at 4 °C for 8 h, the protein in the buffer was extracted using MagneHis Protein Purification System (Promega). The compounds were extracted from the supernatant using ethyl acetate, while the MagneHis Ni-Beads were washed twice with 50 mM Tris⋅HCl (pH 7.5), and the bound chemicals were eluted with ethyl acetate. The chemicals obtained from the supernatant and Ni-Beads were analyzed by LC-QQQ-MS/MS, as described above.

Differential Scanning Fluorimetry Assays.

The parameters (K and ΔTmmax) of CHIL1 binding with DMX or NC were measured using a StepOne real-time PCR system (Applied Biosystems) as previously described (25, 38). The reaction assays contained 5.0 μL of Protein Thermal Shift Buffer, 2.5 μL of 8× Diluted Protein Thermal Shift Dye, 5 μM purified HlCHIL protein, in a final volume of 20 μL. The chemicals used in these assays were serially diluted from 250 μM to 0.12 μM (total 12 points), and K and ΔTmmax were calculated using GraphPad Prism 5 from at least three biological replicates.

Homology Modeling and Docking Simulations.

All structure modeling was performed using the Maestro software package from Schrödinger LLC. The homology model was genen class="Species">rated using Prime with Atpan> class="Gene">FAP1 (PDB ID code 4DOO) as the template model. Hydrogen atoms were added to the model and minimized before docking simulations. Potential binding pockets were identified with SiteMap and the DMX ligand was prepared using LigPrep. Individual DMX docking poses were obtained and scored using Glide, then assessed for selection of the best fit. PyMOL was used for the generation of the final model images.