Open Access Article
Victoria Sodréa,
Goran M. M. Rashida,
Boriana V. Yotzovaab and
Timothy D. H. Bugg
*a
aDepartment of Chemistry, University of Warwick, Coventry CV4 7AL, UK. E-mail: T.D.Bugg@warwick.ac.uk
bDickinson College, Carlisle, Pennsylvania, USA
First published on 9th June 2026
Pyridines are found in many pharmaceuticals and agrochemicals, but are synthesised from fossil fuel conversion. 2,4- and 2,5-pyridinedicarboxylic acids have been reported as products from bioconversion of renewable lignin feedstocks using engineered strains of Rhodococcus jostii RHA1 (Z. Mycroft et al., Green Chem., 2015, 17, 4974–4979), but previously it has been uncertain whether the formation of the pyridine ring was assisted by enzyme catalysis. The 4,5-extradiol ring fission product of protocatechuic acid, 4-carboxy-2-hydroxymuconate 6-semialdehyde (CHMS) shows structural similarity to α-ketoglutaric acid, the substrate for reductive amination by glutamate dehydrogenase (GDH). Testing of five glutamate dehydrogenase (GDH) isozymes from R. jostii RHA1 revealed that GDH5 catalyses NADH-dependent reductive amination of CHMS, and its cyclisation to form a dihydropyridine product. The dihydropyridine can be oxidised to 2,4-pyridinedicarboxylic acid using recombinant P. fluorescens dye-decolorizing peroxidase DyP1B, providing a route to substituted pyridines from a renewable feedstock.
Our group has previously published the conversion of polymeric lignin or lignocellulose by engineered Rhodococcus jostii RHA1 into pyridine-dicarboxylic acid (PDCA) bioproducts.9 Instead of the normal 3,4-oxidative cleavage of key intermediate protocatechuic acid, insertion of genes encoding either protocatechuate 4,5-dioxygenase (encoded by ligAB genes from Sphingobium sp. SYK-6) or protocatechuate 2,3-dioxygenase (encoded by praA from Paenibacillus sp. JJ-1b) generates new extradiol oxidative cleavage products 4-carboxy-2-hydroxymuconate 6-semialdehyde (CHMS) and 5-carboxy-2-hydroxymuconate 6-semialdehyde (5-CHMS), which are cyclised in the presence of ammonium chloride present in M9 minimal media, to generate 2,4-pyridine-dicarboxylic acid (2,4-PDCA) or 2,5-pyridine-dicarboxylic acid (2,5-PDCA) respectively, as shown in Fig. 1.9 These compounds are analogues of terephthalic acid found in plastics such as PET and PBAT, and they can be converted using biocatalysis into PDCA-containing bioplastics.10 Titres of 60–125 mg L−1 were initially obtained by insertion of plasmid-borne ligAB or praA genes into wild-type R. jostii RHA1,9 which have been improved to 330 mg L−1 by gene deletion of the competing beta-ketoadipate pathway, and gene insertion onto the R. jostii chromosome.11
The molecular mechanism responsible for cyclisation of the pyridine ring during PDCA formation in R. jostii RHA1 was not certain.9 Although the non-enzymatic reaction of extradiol ring cleavage products with ammonium salts to form picolinic acids is precedented,12,13 it is slow at pH 7, and the same cyclisation occurs much more slowly in growing Pseudomonas putida cells,14 suggesting that the pyridine cyclisation reaction in R. jostii RHA1 is probably assisted in some way by enzyme catalysis. However, there are few precedents for enzyme-catalysed nitrogen heterocycle formation from acyclic substrates. For 6-membered nitrogen heterocycles, there is precedent for the use of ω-transaminase to convert a 1,5-dicarbonyl substrate into a 1,5-aminoketone, which then cyclises to a cyclic imine, and can then be reduced to a piperidine product.15,16 Cyclisation of a 1,2-diamine with a dicarbonyl compound to form a pyrazine heterocycle is also precedented,17 but there are no known examples of enzyme-catalysed cyclisation to form a pyridine product. For 5-membered nitrogen heterocycles, pyrrole formation is catalysed by transaminase enzyme PigE in the biosynthesis of a prodigiosin in Serratia,18 and by 5-aminolevulinic acid dehydratase in porphyrin biosynthesis,19 and formation of L-proline from L-ornithine is catalysed by ornithine cyclodeaminase.20 Here we report the identification of an enzyme responsible for pyridine ring formation from CHMS and 5-CHMS in R. jostii RHA1.
| Enzyme name | Gene ID | Amino acids | Mr (kDa) | Expressed? | Active with CHMS? |
|---|---|---|---|---|---|
| GDH1A | ro04644 | 1572 | 171 | No | |
| GDH1B | ro03717 | 1542 | 168 | No | |
| GDH1C | ro01009 | 1527 | 165 | No | |
| GDH2 | ro05607 | 1060 | 114 | No | |
| GDH3 | ro03471 | 1130 | 123 | No | |
| GDH4 | ro00339 | 382 | 39.5 | Yes | No |
| GDH5 | ro00573 | 429 | 45.5 | Yes | Yes |
| GDH6 | ro03288 | 122 | 12.9 | Yes | Weakly |
| GDH7 | ro08786 | 423 | 44.7 | Yes | No |
| GDH8 | ro01405 | 447 | 48.1 | Yes | No |
For protein purification, cell pellets were resuspended in a final volume of 20 mL IMAC buffer (20 mM sodium phosphate buffer pH 7.4, 300 mM NaCl, 20 mM imidazole), supplemented with 1 mM phenylmethylsulfonyl fluoride (PMSF), 1 U mL−1 DNAse I, and 0.35 mg mL−1 lysozyme. For LigAB and PraA purification, the buffer was further supplemented with 1 mM 1,4-dithiothreitol (DTT). The resuspended cells were incubated for 30 min at room temperature, under agitation, and lysed using a E1061 cell disruptor (Constant Systems Ltd), at 20 kpsi. The lysed cells were centrifuged at 17
000 rpm, 4 °C, for 35 minutes. The supernatant containing soluble proteins was filtered through a 0.22 µm filter into a His GraviTrap™ column (Cytiva), pre-equilibrated with IMAC buffer. The column was washed with IMAC buffer containing increasing concentrations of imidazole: 20, 62.5, 125, 250, and 500 mM. The fractions from 250 and 500 mM washes were joined and concentrated using Amicon® Ultra-15 centrifugal filters, with 10 or 30 kDa cutoff, depending on the enzyme being purified. The concentrated fractions from IMAC purification were then buffer exchanged using PD-10 Sephadex™ G-25 columns (Cytiva) into 20 mM MOPS buffer pH 7.4 containing 150 mM NaCl and 15% glycerol, aliquoted, flash frozen and stored at −80 °C until use. The resulting purified proteins were analysed by SDS-PAGE and quantified by Bradford assay.
LigAB and PraA activities were assayed in 200 µL or 1 mL reactions containing 20 mM Tris-HCl buffer pH 8.0, 0.75–1 mM protocatechuic acid (PCA), and approximately 50 µg of the respective enzyme. The reactions were incubated at 30 °C for up to 30 min. The appearance of a bright yellow colour indicated dioxygenase activity, which was monitored by measuring the absorbance at 350 nm for 5-CHMS and 410 nm for CHMS.
Glutamate dehydrogenase activity was assayed using a Hidex Sense microtitre plate reader in 200 µL reactions containing 50 mM Tris-HCl pH 7.5, 1 mg ml−1 CHMS or 5-CHMS, 0.1 mM NADH (or NADPH), 20 mM NH4Cl, and 25–50 µg GDH5. Assays were incubated at 30 °C for up to 30 min and monitored by measuring the consumption of NADH at 340 nm. For determination of kinetic parameters, assays were carried out at a range of substrate concentrations: 0.625–10 mM for 5-CHMS and 0.03–1.0 mM for CHMS. The amount of NADH in µmoles was determined by fitting the data to the trendline equation derived from a NADH standard curve. Kinetic parameters (KM and kcat) were determined by nonlinear regression fitted to the Michaelis–Menten equation using GraphPad Prism 8.
GABA transaminase activity was assayed in 200 µL reactions containing 50 mM Tris-HCl pH 7.5, 1 mg ml−1 CHMS or 5-CHMS, 1 mM GABA or L-glutamic acid, 50 µM pyridoxal 5′-phosphate, and 50 µg GABA-T enzyme. Reactions were monitored at 380 nm (at this wavelength there was less background absorbance due to PLP) over 30 min for reduction in absorbance of CHMS or 5-CHMS, and compared with control assays lacking enzyme.
Pseudomonas fluorescens Pf-5 Dyp1B was expressed and purified as previously described.21 A total of 0.5 mg and 1 mM H2O2 were added to the GDH5 reaction products resuspended in buffer, to a final volume of 1 mL. The reactions were incubated at 30 °C for 1 hour, freeze-dried, and the resulting solids resuspended in HPLC-grade methanol for HPLC and LC-MS analysis.
:
20 (v/v) acetonitrile
:
ammonium acetate (5 mM) at a flow rate of 1 mL min−1, and a 5 µL sample injection volume. Detection was at 270 nm.
Reverse-phase high-performance liquid chromatography (RP-HPLC) was carried out using an Eclipse XDB-C18 column, 5 µm, 4.6 × 250 mm, (Agilent). The HPLC solvents were water + 0.1% formic acid (solvent A) and methanol + 0.1% formic acid (solvent B). The applied gradient was 5–100% B for 15 min, 100% B for 10 min, 100–5% B for 1 minute, and 5% B for 10 min, with a 1 mL min−1 flow rate and detection at 270 nm.
Liquid chromatography coupled to mass spectrometry (LC-MS) analyses were performed using a Bruker Amazon X ESI mass spectrometer, using the LC method given above. Dihydropyridines were identified by extracted ion analysis for m/z 170 (MH+), in positive ion mode. Pyridine-dicarboxylic acids were detected by extracted ion analysis for m/z 168 (MH+), in positive ion mode, and compared with authentic standards for 2,4- and 2,5-PDCA. NMR spectroscopy was performed on a Bruker Avance III HD 400 MHz instrument, in d4-methanol as solvent.
Analysis of the genome of R. jostii RHA1 revealed that this microbe contains 10 genes encoding putative glutamate dehydrogenase enzymes (see Table 1). Due to the high GC content of Rhodococcus jostii RHA1, protein expression of large proteins in Escherichia coli K12 is very challenging, so we designed constructs for the five GDH genes with enzyme Mr <100 kDa, for expression in Escherichia coli. Recombinant R. jostii RHA1 glutamate dehydrogenases 4–8 were successfully expressed as recombinant His6-tagged proteins, and purified by Ni2+ affinity chromatography (see SI Fig. S1).
The 2,3-extradiol cleavage product 5-CHMS was also generated by expression of recombinant PraA in E. coli, followed by Ni2+ affinity chromatography. Incubation of recombinant PraA with protocatechuic acid in 20 mM Tris pH 8.0 gave rise to the ring fission product 5-CHMS, absorbing at λmax 350 nm as reported previously,27 which was isolated as described above.
Recombinant R. jostii RHA1 glutamate dehydrogenases 4–8 were screened for activity with CHMS, by incubation of each enzyme with PCA, LigAB, ammonium chloride, and 0.1 mM NADH. In the case of GDHs 5 and 6, the formation of new products with similar retention times to 2,4-PDCA was observed by hydrophobic interaction liquid chromatography (HILIC), accompanied by consumption of NADH (see Fig. 3A).
Recombinant R. jostii RHA1 GDH5 expressed well in E. coli, whereas GDH6 expressed poorly (see SI Fig. S1), so GDH5 was selected for further characterisation. As shown in Fig. 3B, time-dependent consumption of NADH was observed at 340 nm upon incubation of RjGDH5 with CHMS, with low background rate. Assaying GDH5 with either NADH or NADPH revealed that both cofactors could be utilised, with a preference for NADH, with highest activity at pH 7.5–8.0 (see SI Fig. S4).
The maximal rate of NADH consumption in the presence of CHMS by GDH5 was >80% of that observed for the natural substrate α-ketoglutarate (see SI Fig. S4). GDH5 also showed activity with 5-CHMS as substrate in the presence of NADH (see SI Fig. S4). GDH5 showed saturation kinetic behaviour with CHMS and 5-CHMS (see SI Fig. S5), and the steady-state kinetic parameters are shown in Table 2. CHMS shows lower estimated KM than 5-CHMS, consistent with the closer structural resemblance to its natural substrate.
| Substrate | KM (µM) | kcat (min−1) |
|---|---|---|
| Kinetic plots shown in SI Fig. S5. | ||
| CHMS | 45 ± 4 | 0.54 ± 0.03 |
| 5-CHMS | 57 ± 5 | 0.38 ± 0.02 |
| α-Ketoglutarate | 8.2 ± 0.4 | 0.94 ± 0.04 |
We also considered the possibility that enzyme-catalysed transamination might be taking place at the aldehyde terminus of CHMS. There is a possible structural resemblance of the aldehyde portion of CHMS to succinate semialdehyde, which is transaminated by γ-aminobutyric acid (GABA) transaminase to form GABA (see Fig. 4). There are four GABA transaminase genes in R. jostii RHA1 (ro04469, ro04544, ro05598, ro08787, GABA-T 1–4 respectively), whereas most bacteria contain only one GABA transaminase. Three R. jostii RHA1 GABA transaminases (GABA-T 1–3 as above) were successfully expressed as recombinant proteins in E. coli. Enzyme assay in the presence of CHMS or 5-CHMS with addition of GABA or L-glutamic acid showed consumption of substrate at 380 nm (see SI Fig. S6), but with a specific activity of only ∼0.2 nmoles per min per mg protein, compared with 10 nmoles per min per mg protein for GDH5. Therefore, although there is some observable activity of CHMS with GABA transaminase enzymes, it is much less than that observed for GDH5.
C
N, but more than one set of peaks was observed, implying the existence of several isomers of the imine product (see SI Fig. S7).
Treatment of the dihydropyridine product with P. fluorescens peroxidase Dyp1B21 and hydrogen peroxide led to the formation of 2,4-pyridinedicarboxylic acid, observed at m/z 168 by LC-MS (see Fig. 5C). Co-injection of the GDH5/Dyp1B product with authentic 2,4-PDCA by HPLC verified their co-elution (see SI Fig. S8), confirming that the observed product is 2,4-PDCA. Analysis of the Dyp1B oxidation product by 1H NMR spectroscopy showed signals for the aromatic hydrogens of 2,4-PDCA at 8.61, 9.04 and 9.34 ppm (see SI Fig. S9).
We examined whether we could detect the non-enzymatic cyclisation of CHMS with ammonium chloride. No 2,4-PDCA was observed in the CHMS/GDH5 reaction after 30 min (see Fig. 5B, m/z 168), which suggested that the non-enzymatic cyclisation was slow. However, control reactions lacking GDH5 did gradually form a smaller peak for 2,4-PDCA (see Fig. 5D, m/z 168), to the extent of 25–30%, indicating the presence of a slower non-enzymatic reaction.
Reaction of GDH5 with 5-CHMS gave only weak signals at m/z 170 by LC-MS analysis (data not shown), consistent with 5-CHMS being a weaker substrate for GDH5, therefore there may be an alternative enzyme responsible for ammonia cyclisation of 5-CHMS to 2,5-PDCA in R. jostii RHA1.
Therefore, we conclude that GDH5 uses CHMS as substrate for reductive amination, resulting in several dihydropyridine isomers (m/z 170 MH+). These reaction products can be oxidised by Dyp1B, generating 2,4-PDCA (m/z 168 MH+). A summary of this proposed pathway is shown in Fig. 6.
We have shown that the dihydro-PDCA intermediate can be oxidised to the pyridine 2,4-PDCA by P. fluorescens Dyp1B.21 R. jostii RHA1 contains two dye-decolorising peroxidases, DypA and DypB, which could catalyse this oxidation in R. jostii RHA1,28 but this strain also contains multi-copper oxidases McoA and McoC whose overexpression in R. jostii RHA1 has been shown to enhance titre of 2,4-PDCA,29 which might also catalyse this oxidation. The cyclisation of CHMS appears to generate more than one cyclic imine product, which could all be oxidised to form 2,4-PDCA.
Nevertheless, RjGDH5 activity could benefit from optimisation via protein engineering, not only to enhance activity, but to tailor its catalytic site to accept a wider range of substrates.
While ammonia cyclisation can happen non-enzymatically, several observations indicate that the non-enzymatic reaction is slower and less efficient than the GDH5-catalysed reaction: low background reaction in the GDH5 enzyme assay (Fig. 3B), and linear dependence of CHMS disappearance at 450 nm vs. GDH5 concentration (Fig. S4, panel F); no apparent formation of 2,4-PDCA in the presence of GDH5 but absence of Dyp1B (Fig. 5, panel B). But 25–30% formation of 2,4-PDCA was observed in the presence of Dyp1B but absence of GDH5 (Fig. 5, panel D), indicating that the background non-enzymatic reaction can be detected. Our conclusion is that the proportion of 2,4-PDCA product formed by non-enzymatic reaction in a microbial conversion will depend upon the amount of GDH5 enzyme present in the cell, or in vitro will depend upon the amount of GDH5 enzyme added.
Extradiol ring fission products are known to be prone to keto–enol tautomerisation,30 which was observed for CHMS and 5-CHMS (data not shown), and can show toxicity towards a microbial host.14 Thus, enzymes that can catalyse ammonia cyclisation of aromatic ring fission products such as RjGDH5 offer a more specific and efficient route to pyridines than non-enzymatic cyclisation, especially under physiological conditions – which is important for biotechnological applications.
A recent study has reported high-yield production of pyrone and pyridine dicarboxylic acids from glucose as feedstock in engineered Corynebacterium glutamicum strains.31 A BLAST search of RjGDH5 amino acid sequence against the C. glutamicum genome identified 17 matches with 25–32% sequence identity. By comparison, only 2 matches were identified when the same BLAST search was applied against the genome of P. putida KT2440. Considering that intense glutamate metabolism is a feature of C. glutamicum, as well as the high number of GDH isomers encoded in its genome, it is possible that some of these enzymes may have contributed towards pyridine cyclisation of the final products.
The identification of an enzyme capable of catalysing cyclisation of a pyridine ring opens up the possibility of generating a wider range of pyridines, from renewable sources, and hence could provide a new biocatalytic route to substituted pyridine chemicals.
| This journal is © The Royal Society of Chemistry 2026 |