P. luminescens is a Gammaproteobacterium mutualistically associated with nematodes, and the pair prey on insect larvae (Waterfield et al., 2009). P. luminescens produces an assortment of bioactive molecules and antibiotics to regulate its mutualistic and pathogenic interactions (Vizcaino et al., 2014b; Challinor and Bode, 2015). Consequently, this genus, which includes one human pathogen, rivals the Streptomyces genus in terms of the number of secondary metabolic pathways in a given genome (Duchaud et al., 2003; Tobias et al., 2016). Using the MicroScope bioinformatics platform, we identified a genomic island (plu2792-plu2799, Figure 1—figure supplement 1), harboring mixed NRPS-pteridine synthase machinery (Figure 1). Protein sequence homology analysis demonstrated predicted pteridine biosynthetic enzymes, such as GTP cyclohydrolase (GTPCH) Ι, 6-pyruvoyltetrahydropterin synthase, pteridine reductase, and pteridine pyrophosphokinase. Interestingly, a NRPS carrier protein (thiolation domain, T) and a condensation domain (C) were genetically fused to a pyruvate dehydrogenase E2-like subunit (plu2796). The pyruvate dehydrogenase complex (E1, E2, and E3 subunits) is a well-studied set of enzymes that converts pyruvate to acetyl-CoA (Patel and Roche, 1990), serving as a key metabolic bridge between glycolysis and the citric acid cycle. A complementary E1-like subunit is also encoded in the pathway. This pathway was not detected by early versions of the antiSMASH algorithm (Blin et al., 2013) used for identifying the biosynthetic pathways of secondary metabolites; however, recent integration of the ClusterFinder algorithm which aids in identifying divergent biosynthetic gene clusters of unknown classes (Cimermancic et al., 2014; Weber et al., 2015) allowed detection of part of the pathway (plu2796-plu2798) as a likely pathway of secondary metabolic enzymes.

Figure 1 with 1 supplement see all

Download asset Open asset

We first cloned the hybrid NRPS-pteridine genomic island without its clustered regulatory protein (Plu2792) and placed it under the control of a phage T7 promoter for heterologous expression in Escherichia coli BAP1 (Pfeifer et al., 2001). Comparison of the culture broths from E. coli harboring the pathway relative to that of an empty vector control (pET28a) revealed a pathway-dependent yellow phenotype (Figure 2B). Initial comparative single quadrupole liquid chromatography-mass spectrometry (LC/MS) analysis of their butanol extracts revealed two dramatically increased peaks in the broth harboring the pathway (Figure 2A). Through extensive 1D- and 2D-NMR, high-resolution-electrospray ionization-quadrupole-time-of-flight-mass spectrometry (HR-ESI-QTOF-MS), and LC/MS co-injections with authentic standards (Figure 2C–E and Figure 2—figure supplements 1 and 2), the chemical structures of the small molecules attributed to these peaks were assigned as known 2-amino-4-oxopteridine small molecules, 7,8-dihydroxanthopterin (3) and pterin (4).

Figure 2 with 2 supplements see all

Download asset Open asset

To identify novel lower abundance advanced metabolites dependent on both the NRPS and pteridine biosynthetic enzymes, we individually constructed nonpolar genetic deletions of every biosynthetic gene in the pathway, maintaining transcriptional and translational control elements for comparative metabolomics. Relative to the control vector, we identified 224 molecular features in butanol extracts that were dependent on the presence of the wild-type pathway (Supplementary file 1A). These include small molecules that are encoded by the pathway, shunt metabolites emerging from the pathway, and host metabolites that are enhanced by the pathway from an undetectable to a statistically significant level. Of these, 114 were dependent on the GTPCH Ι homolog, Plu2793, which is predicted to convert GTP to 7,8-dihydroneopterin triphosphate, initiating pteridine synthesis (Burg and Brown, 1968). The predicted functional redundancy with primary metabolic enzymes is expected to account for the larger number of metabolic perturbations observed for plu2793. Thirty-seven molecular features were dependent on the atypical dehydrogenase E2-NRPS fusion enzyme, Plu2796, and only 12 were dependent on both. The number of wild-type molecular features that were dependent on the remaining biosynthetic enzymes, Plu2794, Plu2795, Plu2797, Plu2798, and Plu2799, are listed in Table 1.

We then conducted tandem MS analysis on the pathway-dependent molecular features for pathway-targeted molecular networking (Vizcaino et al., 2014a). Molecular networking is a powerful approach to network molecules based on tandem MS fragmentation similarities (Watrous et al., 2012). Pathway-targeted molecular networking focuses on networking of metabolites dependent on the presence of a functional pathway (e.g., wild-type versus a secondary metabolic pathway mutant). By assessing and mapping the relative production levels of the metabolites onto the larger wild-type pepteridine network (Figure 3—figure supplements 1 and 2), we could visualize how individual genetic mutations affect overall metabolite distributions at a systems level (Figure 3). White nodes represent metabolites abolished in a given mutant strain, thereby dramatically focusing metabolite discovery efforts. Some deletions led to a decrease in, rather than complete abolishment of, select molecular feature production, indicating that substantial metabolic crosstalk occurs between this specific hybrid secondary metabolic pathway and primary metabolism. With the observed crosstalk, we focused our structural characterization efforts on the smaller number of molecular features that were dependent on the atypical E2-NRPS machinery. We also included pterin (4, green node in Figure 3) as a standard in our network generation to aid in defining metabolites possessing pteridine structural scaffolds. Through this analysis, two prominent metabolites were selected for NMR-based structural characterization (Figure 3).

Figure 3 with 10 supplements see all

Download asset Open asset

Compounds 1 and 2 with molecular ions at m/z 224.1141 ([M+H]⁺, calcd m/z 224.1147, C₉H₁₄N₅O₂) and m/z 210.0987 ([M+H]⁺, calcd m/z 210.0991, C₈H₁₂N₅O₂), respectively, appeared to be structurally related in the molecular network with a 14 Da mass difference likely attributable to CH₂. A bacterial culture in M9-minimal medium (6 l) supplemented with casamino acids (5 g/l) was initiated for structural characterization. LC/MS analysis of the accumulated butanol extracts led to detection of compounds 1 (t_R = 7.3 min) and 2 (t_R = 5.8 min) (Figure 3—figure supplement 3). However, initial isolation attempts proved challenging because of the low solubility and high polarity of these molecules, which co-eluted with the amino acid supplements. Consequently, isolation of 1 and 2 proceeded through larger scale cultivation in M9-minimal medium (32 l). Although a much lower yield was observed, the isolation was streamlined in the absence of amino acid supplements. The clarified medium was lyophilized, and the dried residue was extracted with 50% aqueous methanol (2 l). Flash column chromatography (C₁₈) followed by several rounds of reverse-phase liquid chromatography separation led to isolation of pure compounds 1 (0.8 mg) and 2 (1.2 mg), which we named pepteridines A and B, respectively.

Structural elucidation of pepteridine A (1) and B (2) was achieved through interpretation of 1D- (¹H and ¹³C) and 2D-NMR (gCOSY, gHSQCAD, gHMBCAD, and NOESY) spectral data (Figure 3—figure supplements 4–9 and Table 2). Briefly, ¹H NMR spectral data coupled with gHSQCAD of 1 were used to show that three NH or NH₂ exchangeable protons, three methylene protons, and one methyl signal were present, suggesting that compound 1 is composed of one hydrogenated pteridine scaffold and one distinguished acyl group. Interpretation of gCOSY and gHMBCAD NMR spectra established the two partial structures to be 2-amino-5,6,7,8-tetrahydropteridin-4(3H)-one and a propionyl group. Key HMBC correlations from the propionyl protons and the methylene protons at C-6 in the pteridine scaffold demonstrated the sharing of an amide carbonyl. These correlations established the connectivity of the propionyl group at N-5, thereby characterizing 1 as a 2-amino-5-propionyl-5,6,7,8-tetrahydropteridin-4(3H)-one (Figure 3). Interpretation of NOESY NMR data and comparison of the proton chemical shifts at the pteridine C-6 to cis- and trans-amide analogs supported the cis-amide conformation of the acyl group at N-5 in solution (Figure 3—figure supplement 10) (Lanyon-Hogg et al., 2015). The critical difference between 1 and 2 was the absence of a CH₂ signal in the NMR spectral data, supporting the assumption that 2 contains an acetyl group in place of the propionyl group in 1, and further HMBC NMR analysis and HR-ESI-QTOF-MS data confirmed this.

Combining the bioinformatic, genetic, and comparative metabolomic analyses with the new structures suggests a hybrid biosynthetic route to the pepteridines, consisting of pteridine synthesis from GTP, acyl-synthesis via oxidative decarboxylation of α-keto acids, and NRPS-dependent condensation of these distinct substrates (Figure 4). The presence of predicted GTPCH Ι, 6-pyruvoyltetrahydropterin synthase, and pteridine reductase homologs in the pathway supports the formation and redox control of a tetrahydropterin substrate. The presence of the pyruvate-dehydrogenase-like E1 and E2-NRPS fusion enzymes supports derivation of the pepteridine acyl-appendages from α-keto acids through an analogous dehydrogenase mechanism. It is likely, however, that the pathway interacts with other pteridine and dehydrogenase biosynthetic enzymes encoded in the genome of the heterologous host, E. coli, such as the E3 subunit of the pyruvate dehydrogenase complex, which is required for lipoamide regeneration but absent from the pathway. Rather than producing an acyl-CoA, as in the pyruvate dehydrogenase reaction, the CoA-derived carrier protein arm (T-domain) could be directly primed in place of CoA (Figure 4). Analogous carrier protein priming mechanisms have been proposed for branched chain fatty acid substrate utilization in the formation of N-acyl-amides and the pristinamycin IIa streptogramin antibiotic (Craig and Brady, 2011; Brachmann et al., 2012); and for a glycolicacyl-NRPS extender unit in formation of the naphthyridinomycin antitumor antibiotic (Peng et al., 2012). In contrast to these pathways, we propose that the resulting loaded acyl-carrier protein in pepteridine biosynthesis would be unusually condensed with a free tetrahydropterin substrate by the atypical NRPS C domain to install the cis-amide acyl-linkage. It is currently unclear whether the C domain catalyzes direct cis-amide bond formation or whether a trans-amide is formed and then isomerized to the observed cis-conformation in the pepteridines.

Figure 4

Download asset Open asset

To gain additional support for our proposed biosynthesis, we analyzed pepteridine production in our full series of pathway mutant strains. Differential LC/HR-ESI-QTOF-MS analysis of butanol extracts demonstrated that pepteridine A (1) detection was completely abolished in the GTPCH Ι homolog mutant (Δplu2793), the 6-pyruvoyl-tetrahydropterin synthase homolog mutant (Δplu2794), and the atypical E2-NRPS mutant (Δplu2796), indicating a genetic requirement for both pteridine and NRPS biosynthetic machineries (Figure 5A). Production of pepteridine B (2) was substantially reduced in the Δplu2793 and Δplu2794 strains, and completely abolished in the Δplu2796 mutant strain (Figure 5C). The minor residual production of 2 in these specific pteridine knockout strains is consistent with our observation that pterin (4) substrates can be detected at low abundance in the control strain (pET28a), further highlighting the metabolic crosstalk between this pathway and primary metabolism. To further demonstrate the biosynthetic dependence on the NRPS machinery, a Plu2796 S434A point mutant was generated in the wild-type pathway construct to site-specifically inactivate the NRPS carrier protein domain. Assessment of pepteridine production from this construct in comparison with wild-type, Δplu2796, and control constructs confirmed dependency on the NRPS machinery (Figure 5B,D).

Figure 5

Download asset Open asset

We then supplemented bacterial cultures harboring the wild-type pathway with varying concentrations of free α-ketobutyrate and pyruvate and, as expected, observed enhanced production of the respective pepteridines in a dose-dependent manner (Figure 6A,B and Figure 6—figure supplement 1). Additionally, universally ¹³C-labeled α-ketobutyrate (¹³C₄) supplementation led to ¹³C₃-labeling of pepteridine A (1), as determined by HR-ESI-QTOF-MS (m/z 227.1244 ([M+H]⁺, calcd m/z 227.1248, ¹³C₃¹²C₆H₁₄N₅O₂), further supporting the proposed biosynthesis (Figure 6C and Figure 6—figure supplement 2) and providing a forward path for protein biochemistry studies on this unprecedented enzymatic system. Importantly, these studies link pepteridine production to α-ketoacid substrate availability (i.e., pyruvate and α-ketobutyrate) rather than free acyl-CoA substrates (i.e., acetyl- and propionyl-CoA).

Figure 6 with 2 supplements see all

Download asset Open asset

Photorhabdus bacteria are both pathogens to insects and mutualists to a specific nematode host (Clarke, 2008; Waterfield et al., 2009; Clarke, 2014). Photorhabdus asymbiotica can also cause infections in humans (Gerrard et al., 2004). To enhance their fitness for these variable host-bacteria objectives, Photorhabdus bacteria undergo phenotypic variation, which is controlled by a stochastic invertible promoter switch (Somvanshi et al., 2012). The orientation of the promoter regulates the formation of P- and M-form phenotypic variants. The P-form, typically the dominant variant in wild-type cultures, is pathogenic to insects. The P-form switches to the M-form, a small colony variant, which adheres to specific cells in the nematode intestine. It is thought that the M-form phenotype participates in colonization of its mutualistic nematode host. To determine whether the pepteridines are produced in P. luminescens and in which variant, we analyzed butanol extracts of P. luminescens genetically ‘locked’ in the M- and P-forms using LC/HR-ESI-QTOF-MS. In the M9-base medium, both pepteridines A and B could be detected in the pathogenic P-form phenotypic variant (Figure 7). However, under identical conditions, no production was observed in the M-form. These studies link pepteridine structure to phenotypic variant status and suggest that the pepteridines may participate in P-form biological activities.

Figure 7 with 1 supplement see all

Download asset Open asset

To initiate functional cellular studies, we deleted the pepteridine genetic locus in P. luminescens TT01 on a wild-type background (Δlocus) and a ΔhexA background (ΔhexA/locus) using allelic-exchange mutagenesis. HexA is a LysR-type transcriptional repressor that regulates stilbene production and participates in the insect pathogen, nematode mutualist transition (Joyce and Clarke, 2003; Kontnik et al., 2010), and global transcriptomic analysis indicates that pepteridine genes are upregulated in ΔhexA/hfq strains (Tobias et al., 2016). We cultivated the wild-type and mutant strains in a culture medium based on the high concentrations of free proteinogenic amino acids found in insect hemolymph (Crawford et al., 2010). Twenty-four hour cultures were centrifuged, the cells were rapidly lysed, and the protein fractions were trypsinized for quantitative proteomic analysis. Relative fold changes (wild-type vs Δlocus and ΔhexA vs ΔhexA/locus) were calculated based on quantitative LC/MS-MS analyses (Figure 8). On a wild-type background, deletion of the locus had little effect on the proteome (Figure 8A). However, on a ΔhexA background and consistent with the transcriptomic data, more dramatic proteomic effects were observed, supporting regulation of the pepteridine genomic locus by HexA (Figure 8B). Proteins that participate in pyrone quorum sensing and secondary metabolism were downregulated in the ΔhexA/locus mutant relative to the ΔhexA control. For example, deletion of the pepteridine genetic locus led to a 5-fold decrease of the enzyme Plu4844 in the ΔhexA background (Figure 8B). Plu4844 participates in synthesis of pyrone autoinducers (Brachmann et al., 2013). We similarly observed a 3–4-fold decrease of the enzymes Plu2817, Plu2204, and Plu4187. These enzymes are involved in biosynthesis of the potent phenoloxidase inhibitor rhabduscin (Plu2817) (Crawford et al., 2010, 2012); cinnamic acid (Plu2204, 2207, 2208), a key substrate of the multipotent stilbenes (Joyce et al., 2008); and the polyketide anthraquinone pigments (Plu4186, 4187, 4188, 4192) (Brachmann et al., 2007). Collectively, our proteomic data support that the pepteridine genomic locus is regulated by HexA and positively affects pyrone quorum sensing and select secondary metabolic pathways.

Figure 8

Download asset Open asset

Figure 8—source data 1

Proteins increased in WT vs. WTΔlocus strains by LC-MS/MS.

Intensities from label-free quantification calculated using MaxQuant, averaged for biological triplicate samples, then log₂ transformed. Proteins presented exhibited ≥2 fold increased average intensity in WT compared with WTΔlocus. Two-tailed t-test performed using a cutoff of FDR = 0.01 and 2-fold signal difference, or p<0.05.

https://doi.org/10.7554/eLife.25229.029

Download elife-25229-fig8-data1-v3.xlsx

Figure 8—source data 2

Proteins increased in WTΔlocus vs. WT strains by LC-MS/MS.

Intensities from label-free quantification calculated using MaxQuant, averaged for biological triplicate samples, then log₂ transformed. Proteins presented exhibited ≥2 fold increased average intensity in WTΔlocus compared with WT. Two-tailed t-test performed using a cutoff of FDR = 0.01 and 2-fold signal difference, or p<0.05.

https://doi.org/10.7554/eLife.25229.030

Download elife-25229-fig8-data2-v3.xlsx

Figure 8—source data 3

Proteins increased in ΔhexA vs. ΔhexAΔlocus strains by LC-MS/MS.

Intensities from label-free quantification calculated using MaxQuant, averaged for biological triplicate samples, then log₂ transformed. Proteins presented exhibited ≥2 fold increased average intensity in ΔhexA compared with ΔhexAΔlocus. Two-tailed t-test performed using a cutoff of FDR = 0.01 and 2-fold signal difference, or p<0.05 (+ indicates statistical significance).

https://doi.org/10.7554/eLife.25229.031

Download elife-25229-fig8-data3-v3.xlsx

Figure 8—source data 4

Proteins increased in ΔhexAΔlocus vs. ΔhexA strains by LC-MS/MS.

Intensities from label-free quantification calculated using MaxQuant, averaged for biological triplicate samples, then log₂ transformed. Proteins presented exhibited ≥2 fold increased average intensity in ΔhexAΔlocus compared with ΔhexA. Two-tailed t-test performed using a cutoff of FDR = 0.01 and 2-fold intensity difference, or p<0.05 (+ indicates statistical significance).

https://doi.org/10.7554/eLife.25229.032

Download elife-25229-fig8-data4-v3.xlsx

Figure 8—source data 5

All proteins observed in WT and WTΔlocus strains by LC-MS/MS.

Intensities from label-free quantification calculated using MaxQuant and log₂ transformed. Intensities for biological triplicate samples are presented. Average fold change shown for WT compared with WTΔlocus. Two-tailed t-test performed using a cutoff of FDR = 0.01 and 2-fold signal difference, or p<0.05 (+ indicates statistical significance).

https://doi.org/10.7554/eLife.25229.033

Download elife-25229-fig8-data5-v3.xlsx

Figure 8—source data 6

All proteins observed in ΔhexA and ΔhexAΔlocus strains by LC-MS/MS.

Intensities from label-free quantification calculated using MaxQuant and log₂ transformed. Intensities for biological triplicate samples are presented. Average fold change shown for ΔhexA compared with ΔhexAΔlocus. Two-tailed t-test performed using a cutoff of FDR = 0.01 and 2-fold signal difference, or p<0.05 (+ indicates statistical significance).

https://doi.org/10.7554/eLife.25229.034

Download elife-25229-fig8-data6-v3.xlsx

Bacterial natural products represent a rich source of lead structures in small molecule drug discovery efforts and serve as excellent molecular probes in biology. Novel classes of putative biosynthetic enzymes, including enzymes belonging to new types of “hybrid" metabolic pathways, can be identified in expanding genome sequence databases through genome mining approaches. Integration of two or more distinct types of metabolic pathways provides combinatorial biosynthetic routes to expand the structural and functional properties of natural products. In this study, using genome synteny analysis, we identified an unprecedented genomic island in P. luminescens that harbored both NRPS-like and pteridine biosynthetic machineries. This hybrid biosynthetic pathway is responsible for synthesis of new pteridine metabolites, the pepteridines. Pepteridine A (1) and B (2), which respectively harbor atypical cis-amide functionalized C₃ and C₂ acyl-chains off a free pterin substrate, a particularly unusual feature for a NRPS, were extensively characterized by NMR. The genetic determinants of pepteridine production were assessed in a heterologous host and networked at the systems level, providing global pathway-dependent maps to visualize relative metabolite distributions among wild-type and mutant strains. The pepteridine metabolites were also identified in the wild-type host, P. luminescens, specifically associated with its phenotypic variant linked to pathogenesis (the P-form). Quantitative proteomic analysis further demonstrated that when the LysR-type transcriptional repressor is derepressed (ΔhexA), the pepteridine genetic locus positively affects pyrone qurourm sensing and select secondary metabolic pathways (i.e., these pathways were downregulated in the ΔhexA/locus strain relative to the ΔhexA strain). The pepteridines likely participate in chemical signaling. Indeed, in the Gram-negative plant pathogen Agrobacterium tumefaciens, pterin signaling was recently implicated in biofilm regulation (Feirer et al., 2015). The pepteridines discovered here by targeting atypical genomic sequence space, to the best of our knowledge, represent the first metabolites to be characterized from a hybrid NRPS-pteridine biosynthetic gene cluster. Characterization of these types of atypical hybrid pathways expands our view on the combinatorial biosynthetic potential available in nature’s metabolic toolbox.

The genomic island spanning genes plu2793 to plu2799 from Photorhabdus luminescens TT01 were amplified via PCR using Phusion High-Fidelity DNA Polymerase (New England Biolabs (NEB), USA) and primers 2796-cluster-5 and 2796-cluster-3 (Supplementary file 1B). Reactions (50 µl) were assembled according to the manufacturer’s protocol with inclusion of 5% (v/v) dimethyl sulfoxide and the use of 1 µl of confluent P. luminescens TT01 as template. Thermal cycling was carried out on a C1000 Touch thermal cycler equipped with a Dual 48/48 Fast Reaction Module (Bio-Rad, USA). Amplification was assessed using a 0.75% agarose gel in 1 × TAE (40 mM Tris pH 7.6, 20 mM Acetic Acid, and 1 mM EDTA) stained with GelGreen (Biotium, USA). Amplification products were purified using the QIAquick PCR Purification Kit (Qiagen, USA) according to the manufacturer’s protocol. The purified product and pET24b (EMDMillipore - Novagen, USA) were individually digested with NdeI and XhoI (NEB), and purified using the QIAquick PCR Purification Kit. Ligation of the products was carried out using T4 DNA ligase (NEB). Ligation mixtures were directly transformed into 50 µl of MAX Efficiency DH5α chemically competent cells (Invitrogen, USA) and recovered in 200 µl super optimal broth with catabolite repression (SOC, 2% (w/v) tryptone, 0.5% (w/v) yeast extract, 10 mM NaCl, 2.5 mM KCl, 10 mM MgCl₂, 10 mM MgSO₄, and 20 mM glucose) following the manufacturer’s protocol. Successful transformants were selected by plating 150 µl of the transformation outgrowth onto lysogeny broth (LB) agar (BD, USA; 1% (w/v) tryptone, 0.5% (w/v) yeast extract, 1% (w/v) NaCl, and 1.5% (w/v) agar) plates supplemented with 25 µg/ml kanamycin (American Bioanalytical, USA) and overnight growth at 37°C. Single, well-defined colonies were picked and grown overnight as suspension cultures at 37°C and 250 rpm in 5 ml LB (BD; 1% (w/v) tryptone, 0.5% (w/v) yeast extract, and 1% (w/v) NaCl) supplemented with 25 µg/ml kanamycin. Plasmids were harvested using the QIAprep Spin Miniprep Kit (Qiagen) according to the manufacturer’s protocol. Incorporation of the biosynthetic pathway encompassing genes plu2793-plu2799 was assessed by end sequencing using the T7 promoter and T7 terminator primers (Genewiz, USA; or Keck Foundation Biotechnology Resource Laboratory at Yale University, USA). Primer walking fully validated the inserted genetic sequence along with the plasmid’s transcriptional and translational control elements proximal to the region of insertion. The construct was named pEplu2796. Sequencing primers are listed in Supplementary file 1C.

Escherichia coli BAP1 (Pfeifer et al., 2001) chemically competent cells were prepared according to standard molecular biology protocols and transformed using 1 µl of the desired construct. Transformants were selected via overnight growth at 37°C on LB agar plates supplemented with 25 µg/ml kanamycin. Resulting single colonies were picked for overnight liquid culture at 37°C and 250 rpm in LB media supplemented with 25 µg/ml kanamycin. Glycerol stocks were subsequently prepared for long-term storage.

The E. coli BAP1 strain carrying the selected pathway was spread onto an LB agar plate supplemented with 25 µg/ml kanamycin and incubated at 37°C overnight. A single colony was inoculated into 5 ml LB liquid media containing 25 µg/ml kanamycin and grown overnight at 37°C and 250 rpm. 5 ml of M9 minimal medium (Amresco, USA; 0.6% (w/v) Na₂HPO₄, 0.3% (w/v) KH₂PO₄, 0.05% (w/v) NaCl, 0.1% (w/v) NH₄Cl, 0.2% (w/v) glucose, 2 mM MgSO₄, and 0.1 mM CaCl₂) supplemented with 5 g/l casamino acids (Amresco) and 25 µg/ml kanamycin was inoculated at 1:200 with the overnight growth, and grown at 37°C and 250 rpm until the optical density (OD₆₀₀) reached 0.5–0.6. Growth was slowed by placing the culture on ice. Production was induced via the addition of isopropyl-1-thio-β-D-galactopyranoside (IPTG, American Bioanalytical) at a final concentration of 0.1 mM and growth proceeded at 25°C and 250 rpm. After 72 hr, the culture was centrifuged (2000 × g, 20 min, 4°C) and the supernatant was then extracted 1:1 with water-saturated butanol (1 × 5 ml). The organic fraction was dried under reduced pressure on a Genevac (USA) HT-4X evaporation system for 2 hr. The crude material, which was re-suspended in 100 µl of 50% aqueous methanol, was analyzed via single quadrapole LC/MS (Column: Phenomenex Kinetex C₁₈ (100 Å) 5 µm (4.6 × 250 mm) column; Flow rate: 0.7 ml/min; Mobile phase composition: water:acetonitrile (ACN) gradient solvent system containing 0.1% formic acid: 0–30 min, 5–100% ACN; hold for 5 min, 100% ACN; 0.1 min, 100–5% ACN; hold for 1.9 min, 5% ACN; 6.1 min re-equilibration post-time, 5% ACN).

All of the genes in the pepteridine gene cluster were individually deleted using PCR to yield a complete series of clean deletion constructs designed to minimize polar effects on the operon and preserve necessary internal transcriptional and translational elements. PCR reactions (50 µl) were prepared using the pEplu2796 construct DNA as a template (~120 ng/µl). The primer pairs used for these reactions are listed in Supplementary file 1B (e.g., Delplu2793F and Delplu2793R - primer pair for the creation of a Δplu2793 construct). Primers were designed according to established protocols (Liu and Naismith, 2008). Briefly, each primer comprised two regions: one 3’ region, which binds immediately after the genetic region to be deleted (bold-faced), and one 5’ region, which binds immediately before the genetic region to be deleted. Standard molecular biology protocols as described above were employed. Deletion of the desired gene was validated by sequencing using the sequencing primer (Supplementary file 1C) immediately upstream of the deleted region (Δplu2793, pET28aUpstream; Δplu2794, Seq T7P*1; Δplu2795, Seq 2; Δplu2796, Seq 6; Δplu2797, Seq 11; Δplu2798, Seq 14; Δplu2799, Seq 15). Constructs possessing the desired deletion sequences were then fully sequence validated.

Comparative metabolomic analyses and pathway-targeted molecular networking were performed following established workflows with minimal modifications, detailed subsequently (Figure 3, Figure 3—figure supplements 1 and 2, Table 1 and Supplementary file 1A) (Vizcaino et al., 2014a; Vizcaino and Crawford, 2015). Biosynthetic pathway expression was carried out in M9 minimal media supplemented with 5 g/l casamino acids, 0.5 g/l L-phenylalanine (Sigma-Aldrich, USA), and 25 µg/ml kanamycin. Expression cultures were inoculated at 1:1000 from stationary phase cultures grown overnight. Seven expression cultures were prepared. Six of these were biological replicates and one was a technical replicate that was solely used to monitor growth. As each sample set reached the desired OD₆₀₀ of 0.5–0.6, it was briefly placed on ice to cool. Once all sample sets had reached density and had sufficiently cooled, expression was then induced with 0.1 mM IPTG, and the culture sets were allowed to grow for 72 hr at 25°C and 250 rpm. The cell mass was then pelleted via centrifugation at 2000 × g for 20 min at 4°C. The supernatant from each sample was collected (~5 ml) and subsequently extracted with 6 ml of water-saturated butanol. The organic layer was collected and dried under reduced pressure on a Genevac HT-4X evaporation system. Dried extracts were stored under nitrogen at −80°C until use. Samples were prepared for HR-ESI-QTOF-MS analysis by resuspension in 200 µl methanol (LC-MS grade, Fisher, USA). Insoluble debris was removed by centrifugation at 20,000 × g for 5 min and 50 µl of the supernatant was placed in an HPLC vial for analysis. The remaining sample was dried and stored as previously described for future analysis. The Agilent iFunnel 6550 QTOF system was used for sample analysis. 2 µl of sample was injected and analyzed at 25°C and 0.7 ml/min on a Phenomenex Kinetex C₁₈ (100 Å) 5 µm (4.6 × 250 mm) column with a water:ACN gradient solvent system containing 0.1% formic acid: 0–30 min, 5–100% ACN; hold for 5 min, 100% ACN; 0.1 min, 100–5% ACN; hold for 1.9 min, 5% ACN; 6.1 min re-equilibration post-time, 5% ACN. Mass spectra were acquired in the range of 25–1700 m/z at a scan rate of 1 spectra/s using Dual Agilent Jet Stream (AJS) ESI in positive mode. Source parameters were set as follows: drying gas temp, 225°C; drying gas flow, 12 l/min; nebulizer pressure, 35 psig; sheath gas temperature, 275°C; sheath gas flow, 12 l/min; fragmentor voltage, 125 V; skimmer voltage, 65 V; OCT 1 RF Vpp, 750 V; capillary voltage, 3500 V; nozzle voltage, 1000 V. Data were acquired and analyzed using MassHunter Workstation Data Acquisition (Version B.05.01, Build 5.01.5125.1, Agilent Technologies) and MassHunter Qualitative Analysis (Version B.06.00, Agilent), respectively. Comparison of the six biological replicates for each sample set allowed for selection of the five with the highest overall chromatographic similarity for further processing. Molecular feature extraction was performed with the following parameter alterations: peak spacing tolerance, 0.0025 m/z plus 7 ppm; limit assigned charge states to a maximum of, 1; restrict charge states to, 1. In importing data into Mass Profiler Professional (MPP, Agilent), ‘minimum number of ions’ was set to 1. In performing the initial abundance analysis on the samples, a minimal normalized abundance of 18 was set. For controls, a molecular feature found in any one of the five samples was tabulated. Conservatively for wild-type, all five biological replicates had to possess the molecular feature for it to be tabulated. To ensure that all molecular features present in mutant strains throughout the five biological replicates were tabulated, a feature was tabulated if found in any one of the five biological replicates. Molecular features found in control samples were removed from both the wild-type and individual mutant feature lists. Wild-type pathway molecular features were subsequently compared withmutant pathway molecular features. Features present in the wild-type and absent in a given mutant were deemed to be dependent on that given mutant. An inclusion list was generated around these features for subsequent pathway-targeted tandem MS (MS²) analysis. MS² analysis was run following the previously described acquisition method with some changes: nebulizer pressure, 50 psig; fragmentor voltage, 200 V. Auto MS² data collection was used with spectra being acquired at 1 spectra/s in the mass range of 25–1700 m/z. Fixed collision energies were set at 0, 25, 40, 45, 50, and 100 V. Auto MS² precursors were limited to those determined to be fully pathway-dependent or dependent on a given enzyme in the pathway. A maximum of 20 precursor ions were analyzed per cycle and a minimum precursor threshold was set at 10,000 counts (absolute) or 0.01% (relative). Isotope models were inactivated and scan speed was varied based on precursor abundance with a target of 25,000 counts/spectrum. The MS² accumulation time limit was employed. All target masses were allowed within a 10 ppm range and a 0.5 min retention time margin of error. Additionally, Auto MS² was used without a preferred ion list on the wild-type sample to establish an untargeted molecular network. The collision energy dependence of the fragmentation patterns of select molecular features was assessed, and we determined that 40 V was optimal for fragmentation of the pepteridines. After data conversion to the mzXML file format, each file was edited to remove non-optimal collision energies while maintaining the preferred 40 V collision energy. This was done as the presence of multiple fragmentation patterns for a single molecular feature can result in the presence of said feature multiple times in the molecular network. In establishing molecular networks through the Global Natural Products Social Molecular Networking Platform (GnPS, http://gnps.ucsd.edu), the following parameters were varied from their default values: parent mass tolerance, 0.001 Da; ion tolerance, 0.5 Da; min cos, 0.5; all filtering was disabled and filter below STD DEV and min peak intensity were both set to 0. Cytoscape (Cytoscape Consortium, USA) was used to visualize and edit the final networks. Parent ion masses networked in GnPS that did not match those present on our initial inclusion list based on precise retention time and high-resolution mass data, were removed from the network. Additionally, all features present in the wild-type inclusion list were validated against the control. Any features that were discovered to be false positives were subsequently removed from the network. When single nodes lacking connectivity were generated during this process, these nodes were removed from the network. The charge state for each feature was evaluated manually against the computationally generated charge state and any inaccuracies were corrected. For all networks, the thickness of the edges connecting nodes is representative of the strength of interaction between the nodes, with thicker lines denoting stronger interactions. The minimum thickness was set at 2 pts representing a cosine score of 0.5 while the thickest possible edge possessing a cosine score of 1 was set to 15 pts.

To decouple the NRPS machinery from that of the fused dehydrogenase E2 subunit in Plu2796, a point mutation in the wild-type pathway was generated (pEplu2796-S434A). In NRPS logic, the biosynthetic pathway intermediates are ferried through the various rounds of enzymatic modification via attachment to a phosphopantetheine arm. This phosphopantetheine tether is post-translationally coupled to a conserved Ser residue. Removal of this Ser destroys carrier protein activity. Bioinformatic assessment of Plu2796 demonstrated a phosphopantetheinyl binding site in the T domain. This conserved Ser was identified as S434 by pattern matching with known amino acid motifs indicative of a phosphopantetheinyl attachment site. Using the mutagenesis strategy outlined in ‘Construction of scar-less deletion constructs’ with the modifications detailed subsequently, S434 was mutated to A434. This mutation prevents the attachment of phosphopantetheine thereby inactivating the NRPS functionality of the Plu2796 enzyme. The primers PluS434AF and PluS434AR were used in the PCR reactions. The wild-type pEplu2796 vector was used as a template. Successful point mutagenesis was assessed via sequencing using the Seq 8 primer (Supplementary file 1C). The pEplu2796-S434A construct was then fully sequence validated over the region harboring the biosynthetic pathway. This construct was transformed into E. coli BAP1 (see ‘Preparation of constructs for heterologous expression studies’). Subsequently, LC-HR-ESI-QTOF-MS analysis was performed on this construct in comparison with the wild-type (pEplu2796), negative control (pET28a), and Δ2796 constructs according to the instrument parameters and chromatography method outlined in ‘Comparative metabolomic profiling and molecular networking.’ Comparison of extracted ion chromatograms demonstrated the dependence of pepteridine A and B on the NRPS thiolation domain present in Plu2796. Extracted ion chromatograms of the S434A mutant match with those for the full deletion of Plu2796. All experiments were conducted in triplicate.

A 5 ml LB liquid culture supplemented with 25 µg/ml kanamycin was initiated by inoculation of a single colony of the E. coli BAP1 strain carrying the wild-type pepteridine construct, pEplu2796. After overnight growth at 37°C and 250 rpm, the culture was used to seed additional 32 × 5 ml LB cultures, which were further incubated at 37°C and 250 rpm for 18 hr. Each of the 32, 5 ml cultures was used to inoculate one of 32, 1 l cultures containing M9 minimal medium supplemented with 25 µg/ml kanamycin. These cultures were incubated at 37°C and 250 rpm until the OD₆₀₀ reached 0.5–0.6. Pathway expression was initiated via the addition of 0.1 mM IPTG, and the cultures were further incubated at 25°C and 250 rpm for 72 hr. The whole 32 l culture volume was centrifuged at 14,000 g and 4°C for 30 min, and the supernatant was lyophilized (~7 days). The dried sample was extracted with a total of 2 l of 50% aqueous methanol, filtered, and evaporated under reduced pressure to yield the crude material (approximately 4.0 g).

The crude organic extract (4.0 g) was subjected to flash C₁₈ column chromatography (300 g) with a step gradient elution (0%, 10%, 20%, 40%, 60%, and 100% methanol in water). The 20% methanol fraction was further fractionated using an Agilent Prepstar HPLC system (Agilent Polaris C₁₈-A 5 µm (21.2 × 250 mm) column) with a linear gradient elution (5–50% methanol in water over 60 min, 8 ml/min, 1 min fraction collection window). A combined fraction (27 + 28) was subsequently isolated by reversed-phase HPLC (Phenomenex Luna C₁₈ C₈(2) (100 Å) 10 µm (10.0 × 250 mm) column) with a linear gradient elution (5–100% ACN in water over 60 min, 4 ml/min) to yield impure compound 1. Compound 1 (t_R = 19.85 min, 0.8 mg) was then purified by reversed phase C₈-HPLC (Phenomenex Luna C₈(2) (100 Å) 10 µm (10.0 × 250 mm) column) with isocratic separation (10% methanol in water, 2 ml/min), followed by final purification over an Agilent Phenyl-Hexyl 5 µm (9.4 × 250 mm) column using isocratic separation (5% ACN in water, 2 ml/min) to yield 0.8 mg of pure 1. For compound 2, the 10% methanol fraction from the flash C₁₈ column chromatography step was separated on a semi-preparative reversed-phase HPLC system (Agilent Polaris C₁₈-A 5 µm (21.2 × 250 mm) column) using isocratic separation (5% ACN in water with 0.01% TFA, 8 ml/min, 1 min fraction collection window). A pooled fraction (13 + 14) was then isolated by reversed-phase HPLC (Phenomenex Luna C₁₈(2) (100 Å) 5 µm (4.6 × 150 mm) column) using an isocratic solvent system (5% ACN in water over 30 min, 2 ml/min) to afford impure compound 2 (t_R = 15.86 min). Final purification of compound 2 (1.2 mg, t_R = 6.51 min) was carried out over an Agilent Phenyl-Hexyl 5 µm (9.4 × 250 mm) column using isocratic separation (5% ACN in water, 2 ml/min). Additionally, a combined fraction (10 + 11) containing compound 3 was further isolated by reversed-phase HPLC (Phenomenex Luna C₈(2) (100 Å) 10 µm (10.0 × 250 mm) column) using an isocratic solvent system (10% ACN in water over 30 min, 2 ml/min). Compound 3 (t_R = 9.5 min) was then purified via an Agilent Phenyl-Hexyl 5 µm (9.4 × 250 mm) column eluting with an isocratic solvent system (3% methanol in water, 2 ml/min) to yield 3 mg of pure material.

Pepteridine A (1): white powder; UV (CH₃OH) λ_max (log ε) 282 (3.68), 222 (3.85) nm; ¹H and ¹³C NMR spectra, see Table 2; HR-ESI-QTOF-MS [M+H]⁺m/z 224.1141 (calcd for C₉H₁₄N₅O₂, 224.1147).

Pepteridine B (2): white powder; UV (CH₃OH) λ_max (log ε) 282 (3.48), 222 (3.75) nm; ¹H and ¹³C NMR spectra, see Table 2; HR-ESI-QTOF-MS [M+H]⁺m/z 210.0985 (calcd for C₈H₁₂N₅O₂, 210.0991).

7,8-Dihydroxanthopterin (3): white powder; UV (CH₃OH) λ_max (log ε) 310 (2.99), 272, (3.25), 222 (3.54) nm; ¹H NMR (DMSO-d₆, 600 MHz) δ 10.32 (1H, br s, 3-NH), 8.80 (1H, br s, 5-NH), 6.60 (1H, br s, 8-NH), 6.14 (2H, br s, 2-NH₂), 3.73 (2H, s, H-7); ¹³C NMR (DMSO-d₆, 100 MHz) δ 161.9 (C-6), 154.4 (C-4), 152.0 (C-2), 150.7 (C-8a), 93.4 (C-4a), 46.0 (C-7); HR-ESI-QTOF-MS [M+H]⁺m/z 182.0631 (calcd for C₆H₈N₅O₂, 182.0678).

Pterin (4): pale yellow powder; UV (CH₃OH) λ_max (log ε) 330 (2.48), 272 (3.26), 226 (3.46) nm; HR-ESI-QTOF-MS [M+H]⁺m/z 164.0574 (calcd for C₆H₆N₅O, 164.0572).

The chemical structures of metabolites 1–4 were identified by analyses of NMR and HR-ESI-QTOF-MS data; and spectral comparison with validated standards. The structure of compound 3, which was elucidated by 1D- and 2D- (gCOSY and gHMBCAD) NMR experiments, was identified as the 7,8-dihydroxanthopterin (3) and supported by comparing its NMR spectral data with a commercial standard of 3. The structure of pterin (4) was unambiguously characterized by HR-ESI-QTOF-MS data, HPLC co-injection with standard pterin, and by comparison of the UV absorption spectral data with that of a pterin commercial standard. Compound 1 was isolated as a white powder. The molecular formula was assigned as C₉H₁₃N₅O₂ ([M+H]⁺ at m/z 224.1141) based on HR-ESI-QTOF-MS spectroscopic data. The ¹H NMR spectral data recorded in DMSO-d₆ displayed two NH protons [δ_H 10.07 (1H, br s), 6.96 (1H, br s)], an NH₂ proton [δ_H6.25 (2H, br s)], three methylene groups [(δ_H 4.54, 2.33), (δ_H 3.30, 2.79), (δ_H 2.57, 2.15)], and one methyl triplet (δ_H 0.88). Interpretation of the HSQC data coupled with ¹³C NMR spectral data showed a total of nine signals, which allowed us to assign all protons to the four directly bonded carbons (δ_C 41.9, 38.5, 26.5, and 9.7), together with the resonances of five quaternary carbons (δ_C 174.0, 157.3, 154.6, 153.1 and 93.1), suggesting the presence of a 2-amino-4-oxo-tetrahydropteridine moiety and a propionyl group. The sequential COSY correlations from a NH proton (δ_H 6.96) to a methylene group (δ_H 4.54, 2.33) established a dimethylene diamine-type partial structure and the COSY cross-peaks between a triplet methyl signal (δ_H 0.88) and a methylene group (δ_H 2.57, 2.15) also supported the presence of the propionyl group. The HMBC correlations from a NH proton (δ_H 6.96) and a methylene group (δ_H 2.57, 2.15) to a quaternary carbon (δ_C 93.1) allowed us to construct a 1,2,3,4-tetrahydropyrazine ring system and the three bond HMBC correlation sharing a carbonyl amide (δ_C 174.0) from a methyl signal (δ_H 0.88) and a methylene group (δ_H 4.54, 2.33) led to attachment of the propionyl group to the tetrahydropyrazine ring via N-acylation at the 5-position of the pteridine ring. NOESY interpretation and the presence of the shifted methylene protons (δ_H 4.54, 2.33), which are adjacent to the amide group, supported formation of the N-acyl group corresponding to the cis-amide conformation in the rotamer system. The full structure of compound 1 was unambiguously characterized by the HR-ESI-QTOF-MS data analysis. Compound 2 was also isolated as a white powder. HR-ESI-QTOF-MS analysis showed that the molecular formula of 2 was C₈H₁₁N₅O₂ ([M+H]⁺m/z 210.0985) possessing a 14 Da mass difference, which could arise from loss of CH₂ from 1. The 1D- (¹H and ¹³C) NMR spectral data of 2 were almost identical to that of 1 except for the absence of the CH₂ signal, indicating the presence of an acetyl group instead of a propionyl group characterized from the NMR interpretation of compound 1, which is supported by HMBC correlation from a singlet methyl proton (δ_H 1.97) and a methylene group (δ_H 4.52, 2.32) to a carbonyl amide (δ_C 170.6).

α-Ketobutyrate and pyruvate were purchased from Sigma-Aldrich. E. coli BAP1 transformed with the pepteridine pathway (pEplu2796) was plated onto LB agar containing 25 µg/ml kanamycin and grown overnight at 37°C. A single, well-defined colony was selected and cultured in 5 ml LB medium supplemented with 25 µg/ml kanamycin overnight at 37°C and 250 rpm. Five milliliters of M9 minimal media supplemented with 25 µg/ml kanamycin and either filter sterilized α-ketobutyrate or pyruvate at multiple concentrations were inoculated with 25 µl of overnight culture (1:200) and incubated at 37°C until the OD₆₀₀ was between 0.5 and 0.6. The M9 cultures were then induced by addition of 0.1 mM IPTG and incubated at 25°C and 250 rpm for 72 hr. The cultures were centrifuged (2000 g, 20 min, 4°C), and the supernatants were extracted with water-saturated butanol (1 × 5 ml). The organic layers were dried under reduced pressure using a Genevac HT-4X evaporation system, and the crude materials were resuspended in 100 µl of 50% aqueous methanol. The resuspensions were then analyzed by single quadrapole LC/MS (Column: Phenomenex Kinetex C₁₈ (100 Å) 5 µm (4.6 × 250 mm) column; flow rate: 0.7 ml/min; mobile phase composition: water:ACN gradient solvent system containing 0.1% formic acid: 0–30 min, 5–100% ACN; hold for 5 min, 100% ACN; 0.1 min, 100–5% ACN; hold for 1.9 min, 5% ACN; 6.1 min re-equilibration post-time, 5% ACN). The feeding experiments were conducted in triplicate.

¹³C₄-α-ketobutyrate was purchased from Cambridge Isotope Laboratories (USA). A single colony of E. coli BAP1 carrying the pepteridine pathway (pEplu2796) was inoculated into 5 ml LB medium supplemented with 25 µg/ml kanamycin. Three biological replicates were prepared, and the cultures were incubated at 37°C and 250 rpm overnight. Five milliliters of M9 minimal media supplemented with 25 µg/ml kanamycin and filter sterilized ¹³C₄-α-ketobutyrate (1.2 g/l) was inoculated with 25 µl of overnight culture and incubated at 37°C until the OD₆₀₀ was between 0.5 and 0.6. The M9 cultures were then induced by addition of 0.1 mM IPTG and incubated at 25°C and 250 rpm for 72 hr. The cultures were centrifuged (2000 × g, 20 min, 4°C), and the supernatants were extracted with water-saturated butanol (1 × 5 ml). The organic layers were dried under reduced pressure using a Genevac HT-4X evaporation system, and the crude materials were resuspended in 100 µl of 50% aqueous methanol. The resuspensions were then analyzed on the Agilent iFunnel 6650 QTOF system (Column: Phenomenex Kinetex C₁₈ (100 Å) 5 µm (4.6 × 250 mm) column; flow rate: 0.7 ml/min; mobile phase composition: water:ACN gradient solvent system containing 0.1% formic acid: 0–30 min, 5–100% ACN; hold for 5 min, 100% ACN; 0.1 min, 100–5% ACN; hold for 1.9 min, 5% ACN; 6.1 min re-equilibration post-time, 5% ACN).

Genetically locked P. luminescens in the M- and P- forms (Somvanshi et al., 2012) were grown on Luria-Bertani agar plates at 30°C. Single colonies of M- and P-forms were inoculated into individual 5 ml LB liquid medium and cultivated in a shaking incubator for 48 hr (30°C, 250 rpm). Five milliliters of each culture broth were then transferred to 1 l scale of M9 minimal medium (0.2% (w/v) glucose, 2 mM MgSO₄, and 0.1 mM CaCl₂) supplemented with 5 g/l casamino acids and cultivated at 30°C and 250 rpm. After 72 hr, 1 l each of M- and P-form culture broths were centrifuged at 14,000 × g for 20 min, and the supernatant was then extracted with butanol (2 × 1 l). The butanol-soluble fractions were separately dried under reduced pressure. The samples were subsequently resuspended in 5 ml methanol and 1–2 µl of each sample was injected for HR-ESI-QTOF-MS analysis (Column: Phenomenex Kinetex C₁₈ (100 Å) 5 µm (4.6 × 250 mm) column; flow rate: 0.7 ml/min; mobile phase composition: water:acetonitrile (ACN) gradient solvent system containing 0.1% formic acid: 0–30 min, 5–100% ACN; hold for 5 min, 100% ACN; 0.1 min, 100–5% ACN; hold for 1.9 min, 5% ACN; 6.1 min re-equilibration post-time, 5% ACN). Extracted ion count chromatograms were obtained by extracting with m/z 224.1147 and 210.0991 corresponding to pepteridines A and B, respectively, with a 10 ppm mass window. The butanol-soluble fraction sample from genetically locked P-form was used in co-injection experiments with standards for confirmation by HR-ESI-QTOF-MS using the same analytical methods.

Allelic-exchange mutagenesis was used to excise the pepteridine genetic locus in P. luminescens TT01 (pluCDS3313775R-plu2799)(Crawford et al., 2010; Kontnik et al., 2010). Approximately 1.5 kB of both upstream and downstream genomic sequences were amplified using the primer pairs Up-F and Up-R; and Dwn-F and Dwn-R, respectively. PCR reactions (50 µl) were prepared using Q5 High-Fidelity DNA Polymerase according to the manufacturer’s protocols with the inclusion of 5% (v/v) DMSO and 1 μl of confluent, wild-type P. luminescens TT01 as template. Products were purified using the QIAquick PCR Purification Kit following the manufacturer’s protocol. Overlap extension PCR was then used to fuse the fragments (Ho et al., 1989). The homologous recombination cassette and pDS132 were both digested with SacI-HF (NEB) (Philippe et al., 2004). rSAP (NEB) was included in the pDS132 digest to dephosphorylate the vector. The linearized vector and cassette were ligated using T4 DNA ligase. Ligation products were transformed into chemically competent E. coli DH5α λpir using a standard heat-shock transformation procedure. Positive transformats were selected by plating 150 μl of the outgrowth onto LB agar plates supplemented with 25 μg/ml chloramphenicol (American Bioanalytical) and grown overnight at 37°C. Colony PCR (cPCR) was used to screen for positive constructs. cPCR reactions were carried out as described above using the primer pair pDS-F and pDS-R. A positive construct was subsequently sequence validated using primers pDS-F, pDS-R, CDSSeq-1, CDSSeq-2, 99Seq-1, and 99Seq-2. This construct was transformed by heat-shock into the diaminopimelic acid (DAP) auxotrophic donor strain E. coli WM6026 λpir (Blodgett et al., 2007) as described above with additional inclusion of 0.3 mM DAP in the LB agar plates. This strain along with wild-type P. luminescens TT01 and P. luminescens Rif^rΔhexA (generously provided by Professor David Clarke) were grown to confluency in LB media at 37°C and 30°C, respectively, and 250 rpm. All overnight growths were subcultured 1:1000 into LB and grown to OD₆₀₀ ≈ 0.6; mixed at 1:1 and 1:4 (donor:recipient) ratios; and filtered through a 0.4 µm sterile filter. Filter mating was carried out overnight at 30°C on LB agar supplemented with 0.3 mM DAP. The outgrowth was then resuspended in LB and streaked onto LB agar containing chloramphenicol (25 μg/ml). Colored colonies indicative of Photorhabdus were selected and streaked onto LB agar supplemented with 5% sucrose for SacB counterselection. Positive colonies were re-streaked three times and single colonies were validated by cPCR. Positive cPCR products were purified using the QIAquick PCR Purification Kit and sequence validated using the primers LocUp, LocDown, CDSSeq-1, CDSSeq-2, 99Seq-1, and 99Seq-2.

P. luminescens wild-type and ΔhexA strains encoding or lacking the pepteridine synthesis pathway were grown in 5 ml hemolymph-mimetic medium (5 g yeast extract, 10 g NaCl, and proteinogenic amino acids based on hemolymph concentrations) for twenty-four hr at 30°C and 250 rpm. 200 μl pellets of strains were prepared in biological triplicate and frozen at −80°C. Upon thawing, sample processing (protein extraction, alkylation, and trypsin digest) was performed as in Lajoie et al. (2013). Peptides were desalted using a C₁₈ MacroSpin column (The Nest Group). Samples were then dried in a centrifugal vacuum concentrator and dissolved in 22 μl 70% formic acid:0.1% trifluoroacetic acid mixed 3:8 volumetrically. In the same buffer, samples were further diluted to 0.5 μg/μl based on A₂₈₀ measurements, and 2.5 μg of peptides (5 μl) were analyzed using an ACQUITY UPLC M-Class (Waters) paired with a Q Exactive Plus (Thermo) mass spectrometer. Column parameters, gradient profiles, and settings for mass spectrometry are described in Ferdaus et al. (2016).

The P. luminescens annotated proteome was downloaded from Uniprot (proteome ID UP000002514) and used for mass spectra searches using MaxQuant v.1.5.1.2 (Cox & Mann) with Cys carbamidomethyl fixed modification and Asn/Gln deamidation, Met oxidation, Ser/Thr/Tyr phosphorylation, and N-terminal acetylation variable modifications. Search parameters considered peptides resulting from two or fewer missed tryptic cleavages and peptides at least five residues long with a 1% false discovery rate. Using Perseus software v.1.5.0.15 (Tyanova et al., 2016), averaged label-free quantification intensities for sample replicates were compared between two samples by two-tailed T test, with significance cutoffs for protein intensity differences established by either permutation-based false discovery rate (FDR = 0.01 with a 2-fold minimum intensity difference) or p<0.05, both reported in Figure 8—source data 1–6.

UV/Vis spectra were obtained on an Agilent (Agilent Technologies, USA) Cary 300 UV-visible spectrophotometer with a path length of 10 mm. ¹H and 2D- (gCOSY, gHSQCAD, gHMBCAD and NOESY) NMR spectral data were recorded on an Agilent 600 MHz NMR spectrometer equipped with a cold probe. ¹³C NMR spectral data were recorded at 100 MHz on an Agilent NMR spectrometer. Flash column chromatography was carried out on Lichroprep RP-18 (40–63 µm, Merck, USA). Routine HPLC analysis was performed on an Agilent 1260 Infinity system with a Phenomenex (USA) Luna C₁₈(2) (100 Å) 5 µm (4.6 × 150 mm) column and a Photo Diode Array (PDA) detector. The separation and purification of metabolites were performed using an Agilent Prepstar HPLC system with an Agilent Polaris C₁₈-A 5 µm (21.2 × 250 mm) column, a Phenomenex Luna C₁₈(2) or C₈(2) (100 Å) 10 µm (10.0 × 250 mm) column, and an Agilent Phenyl-Hexyl 5 µm (9.4 × 250 mm) column. Low-resolution electrospray ionization mass spectrometry (ESI-MS) data were measured on an Agilent 6120 Quadrupole LC/MS system with a Phenomenex Kinetex C₁₈ (100 Å) 5 µm (4.6 × 250 mm) column. High-resolution ESI-MS data were obtained using an Agilent iFunnel 6550 QTOF (quadrupole time-of-flight) MS instrument fitted with an electrospray ionization (ESI) source coupled to an Agilent 1290 Infinity HPLC system.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Our work on the discovery of bioactive metabolites in host-bacteria interactions was supported by the National Institutes of Health (National Cancer Institute grant 1DP2-CA186575 and National Institute of General Medical Sciences grant R00-GM097096), the Searle Scholars Program (grant 13-SSP-210), and the Damon Runyon Cancer Research Foundation (grant DRR 39–16). CEP was supported in part by a National Institutes of Health Chemistry Biology Interface Training Grant (5T32GM067543-12). JR was supported by the National Institute of General Medical Sciences grant R01GM117230 and KWB was supported by National Science Foundation grant DGE-1122492. We thank Yue-Wei Ge for preliminary analysis of pterin and 7,8-dihydroxanthopterin.

1. Received: January 17, 2017
2. Accepted: February 28, 2017
3. Accepted Manuscript published: March 15, 2017 (version 1)
4. Accepted Manuscript updated: March 15, 2017 (version 2)
5. Version of Record published: April 7, 2017 (version 3)

© 2017, Park et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.