Conserved Molecular Mechanism of TyrA Dehydrogenase Substrate Specificity Underlying Alternative Tyrosine Biosynthetic Pathways in Plants and Microbes

Introduction

L-Tyrosine (Tyr) is an aromatic amino acid required for protein synthesis in all organisms, but synthesized de novo in plants and microbes. Thus, in animals Tyr must be acquired through the diet or produced from L-phenylalanine (Phe) by Phe-hydroxylase (Fitzpatrick, 1999). In addition to protein synthesis, Tyr is used to synthesize animal neurotransmitters, such as dopamine and epinephrine (adrenaline) (Fernstrom and Fernstrom, 2007) and melanin skin pigments (Slominski et al., 2004). Tyr also serves as the precursor to numerous plant natural products with diverse functions such as electron carriers (e.g., plastoquinone and ubiquinone; Millner and Barber, 1984), defense (e.g., dhurrin and rosmarinic acid; Petersen, 2013; Gleadow and Møller, 2014), and pollinator attraction (e.g., betalain pigments; Gandía-Herrero and García-Carmona, 2013). Some of these natural products also serve medicinal and nutritional roles in humans such as antioxidants (vitamin E; Falk and Munné-Bosch, 2010), and analgesics (e.g., morphine; Sato et al., 2007).

Tyr is synthesized from prephenate, downstream of the shikimate pathway, by two alternative routes. In most microbes prephenate is first converted into 4-hydroxyphenylpyruvate (HPP) by a NAD⁺-dependent prephenate-specific TyrA dehydrogenase (PDH/TyrA_p), followed by transamination to form Tyr (Figure 1) (Jensen and Pierson, 1975; Bentley, 1990). In contrast, plants first transaminate prephenate to form arogenate, which is converted to Tyr by a NADP⁺-dependent arogenate-specific TyrA dehydrogenase (ADH/TyrA_a) (Gaines et al., 1982; Connelly and Conn, 1986; Rippert and Matringe, 2002). Plant TyrA_a and microbial TyrA_p catalyze the key regulatory step in Tyr biosynthesis, and their substrate specificity defines the Tyr biosynthetic routes via arogenate and prephenate intermediate, respectively (Figure 1).

FIGURE 1

Figure 1. L-Tyrosine biosynthesis via two pathways in plants and microbes. Two pathways for Tyr biosynthesis from prephenate downstream of the shikimate pathway are shown. The PDH/TyrA_p pathway is present in most microbes, whereas the ADH/TyrA_a pathway is ubiquitous in plants. Tyr serves as the precursor for protein synthesis and many downstream metabolites in plants and humans.

Notably, exceptions have been reported for the typical cofactor and substrate specificities of TyrA dehydrogenases in plants and microbes. Some microbes, for example, use NADP⁺ cofactor instead of NAD⁺ (Fazel et al., 1980; Subramaniam et al., 1994). Arogenate-specific TyrA_a enzymes have also been identified in some microbes, such as the α-proteobacteria Zymomonas mobilis and Phenylobacterium immobile (Mayer et al., 1985; Zhao et al., 1993). While all plants investigated have arogenate-specific TyrA_a, the legume family additionally possesses prephenate-specific TyrA_p enzymes (Gamborg and Keeley, 1966; Rubin and Jensen, 1979; Schenck et al., 2015). Using the unique presence of TyrA_p in legumes combined with structural analyses, a single acidic residue in the active site, Asp222, was recently shown to confer arogenate substrate specificity of plant TyrAs by directly interacting with the side chain amine of arogenate substrate (Schenck et al., 2017), which is absent in prephenate (Figure 1). Indeed, mutating Asp222 of diverse plant TyrA_a into the corresponding Asn residue in legume TyrA_p reduced their ADH activity and introduced PDH activity, suggesting that Asn222 played a key role in the recent evolution of legume-specific TyrA_p (Schenck et al., 2017). However, the early evolutionary origin and mechanism of microbial and plant TyrA_a enzymes are still unresolved.

Here, we used the Asp222 residue to trace the evolutionary history of TyrA_a enzymes in deep taxonomic lineages across Plantae and microbes. Structure-guided, phylogenetic analyses combined with biochemical characterization show that microbial TyrA orthologs closely-related to plant TyrA_a also have a corresponding Asp (or Glu) residue and prefer ADH activity. Contrarily, microbial TyrA, which are more distantly-related to plant TyrAs, contained a non-acidic Gln or Asn and preferred PDH activity. Furthermore, site-directed mutagenesis of an acidic Asp into a neutral Asn on a spirochaetes TyrA_a reduced ADH activity, while introducing PDH activity. The reciprocal mutation of Gln into an acidic Glu on an archaeon TyrA_p reduced PDH and introduced ADH activity. These data suggest that plants and some microbial TyrA orthologs share an evolutionarily conserved substrate specificity mechanism, and that acquisition of the key active site acidic residue was crucial in evolution of arogenate-specific TyrA_a enzymes in plants and closely-related microbes.

Materials and Methods

Identification of Microbial TyrA Orthologs

BlastP searches were performed using the amino acid sequences of previously characterized TyrA homologs from plants [soybean PDH; GmPDH1 (Schenck et al., 2015) and Arabidopsis ADH; AtADH2 (Rippert and Matringe, 2002)] and microbes [Synechocystis sp. PCC6803 ADH (Legrand et al., 2006), and E. coli PDH (Hudson et al., 1984)] as the query in the NCBI database. This yielded only closely-related plant and microbial TyrA orthologs (e.g., algae and, γ-proteobacteria), which were then used as the query to perform additional BlastP searches. Every 5th BlastP hit was selected to provide sequences from various microbial lineages and limit bias in sample selection. Data S1 contains all the sequence information for the TyrA orthologs used in Figure 2 and Figure S1. A structure-guided amino acid alignment was performed in PROMALS3D (Pei and Grishin, 2007) using the default parameters with structures of TyrA enzymes from plants and microbes with varying substrate specificities (G. max TyrA_p; GmPDH1; PDB # 5T8X, H. influenzae TyrA_p; HiPDH; 2PV7, and Synechocystis sp. PCC6803 TyrA_a; SynADH; PDB # 2F1K). The amino acid alignment from PROMALS3D was used to construct phylogenetic trees using MEGA7 (Kumar et al., 2016). The full amino acid alignment can be found in Data S2. The analyses involved 130 amino acid sequences and all sites with <75% coverage were eliminated from the analysis. A neighbor-joining method (Figure S1A; Saitou and Nei, 1987) was used to estimate evolutionary history using 1,000 bootstrap replicates (values shown at branches). The tree in Figure 2 is a representative tree using a subset of the sequences found in Figure S1. Additional phylogenetic analyses were performed using the Maximum Likelihood method based on the Jones-Taylor-Thornton (JTT) matrix-based model (Jones et al., 1992), which gave overall similar results (Figure S1B). All phylogenetic trees are drawn to scale, with branch lengths measured in the number of substitutions per site.

FIGURE 2

Figure 2. The conserved acidic residue at 222 among clade I TyrA orthologs from plants, algae, and closely-related bacteria. Structure-guided phylogenetic analysis of plant and microbial TyrAs. Three distinct clades are formed; clade I contains all plant TyrAs and closely-related microbes (blue), clade II contains bacteria, archaea, and fungi TyrAs (green), and clade III, which was used as an outgroup. Enzymes characterized in this study are marked by black arrows. Structures used to guide the alignment are labeled with their PDB IDs. Previously characterized TyrAs are labeled in red with their preferred PDH or ADH activity. Scale bar represents number of substitutions per branch length. A trimmed amino acid alignment of corresponding sequences shows a conserved acidic residue (Asp or Glu, highlighted in blue) among clade I, which is replaced with a non-acidic Asn or Gln residue (highlighted in green) in most clade II at the corresponding 222 position marked with a red star. Identical amino acids present in >50%, black shading; biochemically similar residues present in >50% of the sequences, gray shading.

Recombinant Protein Expression and Purification and Site Directed Mutagenesis

Full length coding sequences from Ochrobactrum intermedium LMG 3301 (EEQ93947.1; OiTyrA), Sediminispirochaeta smaragdinae DSM 11293 (ADK80640.1; SsTyrA), and Methanosaeta harundinacea (KUK94425.1; MhTyrA) were codon optimized for expression in E. coli, gene synthesized (Biomatik), and inserted into pET28a vector using EcoR1 and Nde1 sites in frame with an N-terminal 6x-His tag using a previously described cloning method (Schenck et al., 2015).

For site directed mutagenesis, plasmid template was diluted 100-fold, mixed with 0.04 U/μL Phusion DNA polymerase (Thermo), 0.2 mM dNTP's, 0.5 μM forward and reverse mutagenesis primers and 1x Phusion reaction buffer (Thermo), and then placed in a thermocycler for 98°C for 30 s followed by 20 cycles of 10 s at 98°C, 20 s at 70°C, 4.5 min at 72°C with a final extension at 72°C for 10 min. The sequence of primers used for mutagenesis were (5′-CATTCTGGCCGAAAGCCCGGAACTGTATAGTAGC-3′) and (5′-GTTCCGGGCTTTCGGCCAGAATGCGGCCCACAAAATC-3′) for MhTyrA and (5′-GTAACAATCCACTTCAGCTGTTTATAGACTTGCAAC-3′) and (5′-CTGAAGTGGATTGTTACACGTTTGTTCGCGCACCTG-3′) for SsTyrA (mutated codons are italicized). The PCR products were purified with QIAquick Gel Extraction Kit (Qiagen), treated with DpnI (Thermo) to digest methylated template DNA for 30 min at 37°C, and then transformed into E. coli XL1-Blue cells. Plasmids were sequenced to confirm that no errors were introduced during PCR and cloning.

For recombinant protein expression, E. coli Rosetta2 (DE3) cells (Novagen) transformed with the above plasmids were cultured as previously reported (Schenck et al., 2017). For protein purification, 20 mL of the E. coli supernatant expressing the appropriate plasmid was applied to a 1 mL HisTrap FF column for purification of the His-tagged recombinant protein using an ÄKTA FPLC system (GE Healthcare). After loading the supernatant, the column was washed with 20 column volumes of 90% buffer A (0.5 M NaCl, 0.2 M sodium phosphate, and 20 mM imidazole) and 10% buffer B (0.5 M NaCl, 0.2 M sodium phosphate, and 0.5 M imidazole) followed by elution with 100% buffer B. Fractions containing purified recombinant enzymes were pooled and desalted by Sephadex G50 column (GE Healthcare) size-exclusion chromatography into lysis buffer (Schenck et al., 2017). The purity of purified proteins were analyzed by SDS-PAGE using ImageJ software (Schneider et al., 2012). All protein purification steps were performed at 4°C unless stated otherwise.

ADH and PDH Assays

ADH and PDH assays were performed using purified recombinant enzymes for SsTyrA Wild-type (Wt) and D208N mutant, and MhTyrA Wt and Q227E mutant, while the E. coli cell lysate was used for OiTyrA as expression and purification of this enzyme was unsuccessful. Reactions contained 0.8 mM substrate (arogenate or prephenate) and 0.8 mM cofactor (NADP⁺ or NAD⁺) together with reaction buffer [25 mM HEPES pH 7.6, 50 mM KCl, 10% (v/v) ethylene glycol]. For OiTyrA assays containing cell lysates, reactions were incubated for 45 min and analyzed using HPLC as previously reported (Schenck et al., 2015). For pure enzymes, reactions were monitored every 10–15 s for reduced cofactor at A_340nm using a microplate reader (Tecan Genios) in a reaction volume of 30 μL. Kinetic parameters of purified recombinant enzymes were determined from assays containing varying concentrations of arogenate (39.1 μM−5 mM) or prephenate (46.9 μM−6 mM) substrate, 0.8 mM of the preferred cofactor and monitored every 10–15 s for production of reduced cofactor at A_340nm using a microplate reader (Tecan Genios). Kinetic parameters were determined by fitting initial velocity data to the Michaelis–Menten equation using Origin software (OriginLab) from technical replicate assays (n = 3). Arogenate substrate was prepared by enzymatic conversion of prephenate (Sigma-Aldrich) as previously reported (Maeda et al., 2010). Enzyme assays were quantified using the A_340nm of a standard curve of reduced cofactor (NADPH or NADH) and activity is expressed as Kat/mg (moles of product produced per second per mg protein). All enzyme assays were conducted at a reaction time and protein concentration that were in the linear range and proportional to reaction velocity.

Modeling Microbial TyrA Enzymes

Homology models were made using SWISS-MODEL (Biasini et al., 2014) with default parameters to predict the structures of divergent TyrA enzymes. Enzymes that are more closely-related to plants (e.g., SsTyrA and MhTyrA) were modeled using the GmPDH1 structure as the template, though this resulted in a poor model for BdTyrA, which falls within the outgroup. BdTyrA was additionally modeled using Synechocystis sp. PCC6803 ADH. Homology models were visualized using PyMOL.

Results

Phylogenetic Relationship of Plant and Microbial TyrAs

Previous studies suggested that plant TyrAs are not derived from an eukaryotic ancestor or through cyanobacterial endosymbiosis because they are most similar to other microbes including some proteobacteria (Bonner et al., 2008; Reyes-Prieto and Moustafa, 2012; Dornfeld et al., 2014; Schenck et al., 2017); however, their precise origin was unclear. To resolve the phylogenetic relationship of TyrA orthologs from divergent organisms including plants and microbes, here we performed structure-guided phylogenetic analyses using PROMALS3D to achieve alignment of TyrA orthologs with low sequence similarities (see section Materials and Methods; Pei and Grishin, 2007). Three distinct clades were identified that contain: plant TyrAs together with those from algae, spirochaetes, α- and δ-proteobacteria (clade I, shaded blue in Figure 2, Figure S1), TyrA orthologs from some archaea, fungi, γ-proteobacteria, and chloroflexi (clade II, shaded green), and TyrA orthologs from various microbes, which formed the outgroup and contains previously characterized microbial TyrA orthologs from Synechocystis sp. PCC 6803 and Aquifex aeolicus having very low sequence similarity (~30%) to plant TyrAs (clade III, Figure 2, Figure S1). Interestingly, TyrAs from some spirochaetes lineages (some of which are known to cause harmful human diseases like Lyme disease; Pritt et al., 2016) formed a subclade with plant and algae TyrAs within clade I using various phylogenetic methods (Figure 2, Figure S1). These data suggest that Plantae TyrA may have been acquired through horizontal gene transfer (HGT) from an ancestor of one of these closely-related microbes.

Microbial TyrA Orthologs Containing an Acidic 222 Residue Prefer ADH Over PDH Activity

The amino acid sequence alignment of TyrAs showed that the Asp222 residue, which is conserved across plant TyrA_a (Schenck et al., 2017) was also highly conserved in clade I (Figure 2). On the other hand, most sequences in clade II, including some archaea TyrA, have a non-acidic Gln residue at the corresponding 222 position (Figure 2), similar to legume TyrA_p enzymes (Schenck et al., 2017). Homology models of representative TyrA from clade I—Arabidopsis thaliana ADH (AtADH2, Plantea; Rippert and Matringe, 2002) and S. smaragdinae DSM 11293 (SsTyrA, spirochaetes)—and clade II—M. harundinacea (MhTyrA, archaea)—generated using GmPDH1 structure as the template indeed showed that their acidic and non-acidic residues, respectively, correspond to Asp222 in the active site of plant TyrA (Figure S2). These data together suggest that TyrAs from clade I are likely arogenate-specific TyrA_a enzymes, whereas more distantly-related microbial TyrAs from clade II are likely prephenate-specific TyrA_p enzymes.

To experimentally test if TyrAs from clade I have ADH activity, representative TyrA orthologs from two distinct subclades of clade I, spirochaetes (SsTyrA) and α-proteobacteria (O. intermedium; OiTyrA, Figure 2, Figure S1), were expressed in E. coli as recombinant enzymes and biochemically characterized. SsTyrA and OiTyrA were chosen as they are located at key phylogenetic boundaries within clade I and contain residues required for cofactor binding and catalysis (Figure S2, Data S2). Purified SsTyrA recombinant enzyme showed ADH activity with a slight preference for NAD⁺ over NADP⁺ cofactor; however, PDH activity was not detectable (Figure 3A). Similarly, the E. coli cell lysate expressing OiTyrA had ADH but not PDH activity and strongly preferred NAD⁺ over NADP⁺ cofactor (Figure 3B), although the purification of OiTyrA was not successful due to low expression. These results demonstrate that microbial TyrA orthologs from clade I, which contain an acidic residue at the corresponding 222 position (Figure 2, Data S2), are arogenate-specific TyrA_a enzymes.

FIGURE 3

Figure 3. Substrate and cofactor specificity of microbial TyrA orthologs. ADH and PDH assays were performed with 0.8 mM arogenate and prephenate, respectively, and 0.8 mM cofactor (NADP⁺, black; NAD⁺, gray). (A) Purified recombinant SsTyrA (spirochaetes) was used to test enzymatic activity, and shown as the average in nKat/mg protein ± SEM of n = 3. (B) α-proteobacteria TyrA (OiTyrA) cell lysate was used as purification of the recombinant enzyme was not successful. Average enzymatic activity is shown as pKat/mg protein ± SEM of n = 3 (C) Purified recombinant MhTyrA (archaea) was used to test enzymatic activity, and shown as the average in nKat/mg protein ± SEM of n = 3. N.D. no activity detected. Cofactor preference is indicated by the fold-change over the bars.

An Archaeon TyrA Containing a Non-acidic Residue Prefers PDH Over ADH Activity

To test if TyrA orthologs from clade II, which contain a non-acidic residue at the corresponding 222 position, are prephenate specific TyrA_p enzymes, a representative archaeon TyrA from M. harundinacea (MhTyrA) was biochemically characterized. MhTyrA was chosen as no TyrAs from its subclade of clade II have previously been characterized (Figure 2). Also, MhTyrA is a monofunctional enzyme, while some archaea, fungi, and γ-proteobacteria orthologs in clade II are bifunctional and have a chorismate mutase enzyme domain (Hudson et al., 1984; Shlaifer et al., 2017). MhTyrA was expressed in E. coli and the recombinant enzyme was purified to homogeneity using affinity-chromatography (Figure S3) and used for biochemical analyses. Unlike plant and microbial TyrA_a orthologs from clade I, MhTyrA showed strong PDH and very weak ADH activity (Figure 3C). Interestingly, MhTyrA strongly preferred NADP⁺ over NAD⁺ cofactor (Figure 3C), like plant TyrAs (Gaines et al., 1982; Connelly and Conn, 1986). These results suggest that TyrA orthologs from clade II that have a non-acidic residue at the corresponding 222 position are TyrA_p enzymes that strongly prefer prephenate over arogenate substrate.

A Single D208N Mutation Introduces PDH Activity in a Spirochaetes TyrA_a

To test if the acidic residue at the corresponding 222 position in spirochaetes SsTyrA_a is involved in its substrate specificity, the D208N mutation that converts the active site acidic residue into a non-acidic residue was introduced in SsTyrA. The purified recombinant SsTyrA_a D208N enzyme (Figure S3) drastically decreased its original ADH activity (Figure 4A), but now exhibited PDH activity (Figure 4B), with a slight alteration in cofactor specificity in the mutant compared to Wt (Figure S4).

FIGURE 4

Figure 4. Kinetic analysis of Wt and mutant enzymes of SsTyrA_a and MhTyrA_p. Kinetic analyses were performed with SsTyrA_a Wt (filled circle) and D208N mutant (open square) enzymes with arogenate (A) and SsTyrA_a D208N mutant (open square) with prephenate (B). Kinetic analyses were also performed with MhTyrA_p Wt (filled circle) and Q227E mutant (open square) enzymes using arogenate (C) and prephenate (D). Initial velocity values at each substrate concentration were fit to the Michaelis-Menten equation using Origin software. Kinetic analyses were conducted for MhTyrA_p Wt using 3.41 μg of purified recombinant enzyme, and 4.56 and 2.28 μg of purified recombinant Q227E using prephenate and arogenate, respectively. Kinetic analyses conducted for SsTyrA_a Wt used 0.30 μg of purified recombinant enzyme, and 0.39 μg of purified recombinant D208N. Data are means ± SEM (n = 3 independent experiments). Error bars smaller than symbols are not shown.

Kinetic analyses showed that SsTyrA_a Wt enzyme did not accept prephenate and had a K_m for arogenate of 901 μM (Table 1), which is substantially higher than previously reported plant and microbial TyrA enzymes (Rippert and Matringe, 2002; Bonvin et al., 2006; Ku et al., 2010; Schenck et al., 2015). The SsTyrA_a D208N mutant exhibited a seven-fold lower catalytic efficiency (k_cat/K_m) with arogenate than Wt, though the K_m was not altered (Table 1). SsTyrA_a D208N showed significant PDH activity that was absent in Wt, though its K_m was still much higher than previously reported TyrA enzymes (Rippert and Matringe, 2002; Bonvin et al., 2006; Ku et al., 2010; Schenck et al., 2015) and had poor catalytic efficiency (Table 1). These results suggest that a single mutation of the active site acidic residue to a non-acidic residue can alter the substrate specificity of spirochaetes TyrA, similar to plant TyrAs (Schenck et al., 2017).

TABLE 1

Table 1. Kinetic analysis of MhTyrA_p and SsTyrA_a Wt and mutant enzymes.

A Single Q227E Mutation Introduces ADH Activity in an Archaeon TyrA_p

To test if the non-acidic residue of MhTyrA_p at the corresponding 222 position (Gln227) is involved in prephenate substrate specificity, site-directed mutagenesis was performed on MhTyrA_p to replace Gln227 with acidic Glu and generate the MhTyrA_p Q227E mutant. The purified recombinant MhTyrA_p Q227E enzyme (Figure S3) showed decreased PDH activity (Figure 4D) with a substantial gain of ADH activity (Figure 4C, Table 1) without altering cofactor preference (Figure S4).

Further kinetic analyses showed that Wt MhTyrA_p had a K_m of 378 μM and turnover rate (k_cat) of 2.4 s⁻¹ using prephenate substrate and NADP⁺ cofactor (Figure 4, Table 1), which are comparable to previously characterized microbial TyrA_p enzymes (Bonvin et al., 2006; Ku et al., 2010). The very weak ADH activity of MhTyrA_p Wt (Figure 4, Table 1) precluded it from kinetic analysis using arogenate.

The Q227E mutant, on the other hand, exhibited almost 10-fold reduction in K_m for prephenate (2.7 mM), while the catalytic efficiency (k_cat/K_m) was reduced by 60-fold (0.1 vs. 6.4 mM⁻¹ s⁻¹, Figure 4D, Table 1). The Q227E mutant displayed substantial ADH activity compared to the Wt enzyme with a K_m for arogenate of 3.3 mM, similar to that of Q227E for prephenate (2.7 mM, Figure 4C, Table 1) though still 10-fold higher than that of the Wt enzyme for prephenate (Figure 4D, Table 1) and other previously characterized TyrA_a enzymes (Bonvin et al., 2006; Ku et al., 2010; Schenck et al., 2015, 2017). The Q227E mutant had roughly two-fold higher catalytic efficiency with arogenate than with prephenate (0.2 vs. 0.1 mM⁻¹ s⁻¹, Figure 3). These results demonstrate that the Q227E mutation can shift the substrate preference of MhTyrA_p from prephenate to arogenate, suggesting that the single residue is responsible for substrate specificity of archaea TyrA_p enzymes.

Discussion

Previous studies suggest that microbes predominantly use a PDH-mediated pathway to synthesize Tyr, whereas plants mainly use an ADH-mediated Tyr pathway (Jensen and Pierson, 1975; Bentley, 1990; Siehl, 1999; Rippert and Matringe, 2002; Maeda and Dudareva, 2012; Schenck et al., 2015, 2017). In this study, structure-guided phylogenetic analyses from diverse organisms identified ADH-like sequences in some bacteria, e.g., spirochaetes, α- and δ-proteobacteria, which form a monophyletic clade with plant TyrAs (Figure 2, Figure S1). Biochemical characterization further demonstrated that TyrAs from spirochaetes and α-proteobacteria indeed have ADH, but not PDH activity (Figures 3A,B). A native TyrA enzyme purified from the α-proteobacteria P. immobile, which belongs to the same α-proteobacteria genus found in clade I, was also previously shown to have ADH, but not PDH activity (Mayer et al., 1985). Therefore, our study revealed that arogenate-specific TyrA_a enzymes are more widely distributed in microbes than previously thought.

Previous evolutionary studies revealed that plant aromatic amino acid pathway enzymes are derived from a wide range of, and sometimes unexpected microbial origins (Richards et al., 2006; Reyes-Prieto and Moustafa, 2012; Dornfeld et al., 2014). For example, plant shikimate kinase is most likely derived from cyanobacteria endosymbiosis (Richards et al., 2006) whereas plant prephenate aminotransferase and arogenate dehydratase involved in Phe biosynthesis are sister to Chlorobi/Bacteroidetes orthologs (Dornfeld et al., 2014). However, the evolutionary origin of plant TyrAs is currently unknown. TyrAs from some spirochaetes were more closely-related to plant and algae TyrA_as than other microbial TyrAs from clade I (Figure 2, Figure S1) and, like Plantae TyrA_a enzymes, had a conserved acidic residue at the corresponding 222 position. BlastP searches across different spirochaetes genomes showed that plant-like TyrAs are restricted to the order Spirochaetales, and absent in Leptospirales, Brevinematales, and Brachyspirales (Figure S5; Gupta et al., 2013). Thus, the current result suggests that the common ancestor of algae and plants acquired a TyrA_a enzyme from a spirochaetes ancestor likely through a novel HGT event, rather than from an α-proteobacteria through mitochondria symbiosis (Gray et al., 1999).

The archaeon MhTyrA from clade II preferred PDH over ADH activity (Figure 3C) and had a non-acidic residue at the 222 position (Figure 2, Figure S1). This is consistent with previously-characterized clade II TyrA enzymes from γ-proteobacteria and fungi, which also preferred PDH over ADH activity (Mannhaupt et al., 1989; Christendat and Turnbull, 1999; Chiu et al., 2010) though they belonged to distinct subclades (Figure 2). As almost all TyrA sequences within clade II have a non-acidic residue (Gln or Asn) at the corresponding 222 position, except for Chloroflexi TyrAs (Figure 2, Figure S1), they are likely prephenate-specific TyrA_p enzymes. Previously characterized microbial TyrA_p enzymes had similar K_m with prephenate as MhTyrA_p (Figure 4, Table 1); however the catalytic efficiency of MhTyrA_p was lower than previously characterized TyrA_p enzymes (Bonvin et al., 2006; Ku et al., 2010). This implies that MhTyrA may also have alternative in vivo substrates and further genetic studies are needed to determine the in vivo function of MhThrA_p.

In plant TyrAs, an acidic residue at the corresponding 222 position confers ADH activity by directly interacting with the side chain amine of arogenate, and when mutated to a non-acidic Asn, switches to PDH activity (Schenck et al., 2017). Consistently, in OiTyrA_a and SsTyrA_a, which contain acidic residues at the corresponding 222 position prefer ADH activity (Figure 3). Furthermore, mutation of the corresponding acidic residue into a non-acidic residue on SsTyrA_a introduced novel PDH activity (Figure 4, Table 1). The reciprocal mutation (Gln to Glu) on MhTyrA_p reduced PDH activity while gaining ADH activity (Figure 4, Table 1), further supporting that the corresponding 222 position in microbial TyrA enzymes is also important for their substrate specificity. However, ADH activity of MhTyrA_p Q227E was only two-fold higher than its PDH activity and still 30-fold lower than PDH of MhTyrA_p Wt. Also, PDH activity of SsTyrA_a D208N was seven-fold lower than SsTyrA_a Wt with arogenate (Figure 4, Table 1). These results suggest that residues besides the corresponding 222 substrate specificity determining residue likely contribute to overall catalytic activity of microbial TyrA enzymes. These data together suggest that mutation of the non-acidic to an acidic residue at the corresponding 222 position played a key role in the evolution of arogenate-specific TyrA_a enzymes in microbes from clade I that gave rise to plant TyrAs.

The outgroup (clade III) appears to contain TyrA enzymes with both PDH and ADH activity (Xia and Jensen, 1990; Zhao et al., 1993; Bonvin et al., 2006; Legrand et al., 2006). Homology models of microbial TyrAs from the outgroup (e.g., Bifidobacterium dentium TyrA; BdTyrA) were compared to previously crystallized GmPDH1 (Schenck et al., 2017) and Synechocystis ADH (Legrand et al., 2006) to determine if the substrate specificity mechanism of TyrAs from clade I and II are also conserved in clade III TyrAs (Figure S6). The global conformations of these divergent TyrA proteins from clade I and III are similar in structure, though there are some differences, such as additional α-helices around the C-terminal dimerization domain (Figure S6). All structures have conserved catalytic Ser101 and His124 (Christendat and Turnbull, 1999; Sun et al., 2006) that directly interact with ring hydroxyl of arogenate and prephenate substrate (Schenck et al., 2017), suggesting that the key catalytic residues have been maintained across divergent TyrAs. However, the two loop regions surrounding and recognizing the substrate side chain by the 222 residue (Schenck et al., 2017), are not well conserved in clade III as compared to clade I TyrAs (Figure S6). This makes it difficult to confidently assign a corresponding residue in clade III TyrAs to the 222 position of clade I TyrAs (Figure 2, Data S2). Thus, clade III TyrAs likely use a different molecular mechanism(s) for their substrate specificity than plant and closely-related microbial TyrAs from clade I and II.

Alteration of substrate specificity of microbial TyrA enzymes provides another example of a single active site residue that modifies substrate specificity of plant and microbial enzymes (Louie et al., 2006; He et al., 2011; Fan et al., 2016). Substitutions between active site Val and Phe residues switch the acyl-CoA substrate specificity of Solanum acylsugar acyltransferases that belong to the BAHD acyltransferase family (Fan et al., 2016). A single mutation of the conserved active site Leu residue of isopropylmalate dehydrogenases involved in leucine biosynthesis is sufficient to convert their specificity to 3-(2′-methylthio)ethylmalate, an intermediate of the methionine chain-elongation pathway required for aliphatic glucosinolate biosynthesis (He et al., 2011). Mutating His89 of a microbial Tyr ammonia lyase (TAL) into Phe also switches its substrate specificity to prefer Phe instead of Tyr (Louie et al., 2006). Many of these examples are the result of gene duplications followed by neofunctionalization by a single key amino acid mutation, resulting in a novel substrate specificity and recruitment to a different metabolic pathway (Leong and Last, 2017). Despite the single amino acid mutation that switches substrate specificity of TyrA dehydrogenases, they are still involved in the same Tyr biosynthesis but alters the pathway architecture.

In conclusion, the current study revealed that arogenate-specific TyrA_a enzymes evolved in some bacterial lineages, through the acquisition of an acidic residue at the 222 position, which later gave rise to the TyrAs of algae and land plants likely through a novel HGT event. More recently, the same residue was mutated back to a non-acidic residue uniquely in legume plants, which resulted in prephenate-specific TyrA_p enzymes (Schenck et al., 2017). Thus, in the course of TyrA enzyme evolution, microbial TyrA_p were converted into microbial TyrA_a and then to legume-specific TyrA_p by altering the same active site residue from a non-acidic to an acidic, and then back to a non-acidic residue. Previous studies proposed that the ubiquitous presence of the ADH-mediated Tyr pathway among photosynthetic organisms is to avoid futile cycling of tocopherol and plastoquinone biosynthesis from HPP (Siehl, 1999; Graindorge et al., 2014). Identification of arogenate-specific TyrA among many non-photosynthetic microbes may require revisiting the biological significance of the ADH vs. PDH-mediated Tyr biosynthetic pathways in diverse organisms. Given that arogenate and prephenate substrate specificity of TyrAs can be readily converted by a single residue (Figure 4, Table 1; Schenck et al., 2017), there must be significant selection pressure to maintain the acidic 222 residue and thus ADH activity in many organisms. The molecular mechanism and the key amino acid residue regulating the biochemical properties of diverse TyrAs also enables the optimization of Tyr biosynthesis via two alternative Tyr biosynthetic pathways in both plants and microbes, for enhanced production of pharmaceutically important natural products derived from Tyr (e.g., morphine and vitamin E).

Author Contributions

CAS and YM performed experiments and analyzed data; CAS and HAM conceived the experiments, and CAS wrote the manuscript. All authors read and edited the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the National Science Foundation (IOS-1354971 to HAM).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2017.00073/full#supplementary-material

References

Bentley, R. (1990). The shikimate pathway - a metabolic tree with many branches. Crit. Rev. Biochem. Mol. Biol. 25, 307–384. doi: 10.3109/10409239009090615

PubMed Abstract | CrossRef Full Text | Google Scholar

Biasini, M., Bienert, S., Waterhouse, A., Arnold, K., Studer, G., Schmidt, T., et al. (2014). SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 42, W252–W258. doi: 10.1093/nar/gku340

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonner, C. A., Disz, T., Hwang, K., Song, J., Vonstein, V., Overbeek, R., et al. (2008). Cohesion group approach for evolutionary analysis of TyrA, a protein family with wide-ranging substrate specificities. Microbiol. Mol. Biol. Rev. 72, 13–53. doi: 10.1128/MMBR.00026-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonvin, J., Aponte, R. A., Marcantonio, M., Singh, S., Christendat, D., and Turnbull, J. L. (2006). Biochemical characterization of prephenate dehydrogenase from the hyperthermophilic bacterium Aquifex aeolicus. Protein Sci. 15, 1417–1432. doi: 10.1110/ps.051942206

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, H. J., Abdubek, P., Astakhova, T., Axelrod, H. L., Carlton, D., Clayton, T., et al. (2010). The structure of Haemophilus influenzae prephenate dehydrogenase suggests unique features of bifunctional TyrA enzymes. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 66, 1317–1325. doi: 10.1107/S1744309110021688

PubMed Abstract | CrossRef Full Text | Google Scholar

Christendat, D., and Turnbull, J. L. (1999). Identifying groups involved in the binding of prephenate to prephenate dehydrogenase from Escherichia coli. Biochemistry 38, 4782–4793. doi: 10.1021/bi982673o

PubMed Abstract | CrossRef Full Text | Google Scholar

Connelly, J. A., and Conn, E. E. (1986). Tyrosine biosynthesis in Sorghum bicolor: isolation and regulatory properties of arogenate dehydrogenase. Z. Naturforschung C 41, 69–78.

PubMed Abstract | Google Scholar

Dornfeld, C., Weisberg, A. J. K. C. R., Dudareva, N., Jelesko, J. G., and Maeda, H. A. (2014). Phylobiochemical characterization of class-Ib aspartate/prephenate aminotransferases reveals evolution of the plant arogenate phenylalanine pathway. Plant Cell 26, 3101–3114. doi: 10.1105/tpc.114.127407

PubMed Abstract | CrossRef Full Text | Google Scholar

Falk, J., and Munné-Bosch, S. (2010). Tocochromanol functions in plants: antioxidation and beyond. J. Exp. Bot. 61, 1549–1566. doi: 10.1093/jxb/erq030

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, P., Miller, A. M., Schilmiller, A. L., Liu, X., Ofner, I., Jones, A. D., et al. (2016). In vitro reconstruction and analysis of evolutionary variation of the tomato acylsucrose metabolic network. Proc. Natl. Acad. Sci. U.S.A. 113, E239–E248. doi: 10.1073/pnas.1517930113

PubMed Abstract | CrossRef Full Text | Google Scholar

Fazel, A. M., Bowen, J. R., and Jensen, R. A. (1980). Arogenate (pretyrosine) is an obligatory intermediate of L-tyrosine biosynthesis: confirmation in a microbial mutant. Proc. Natl. Acad. Sci. U.S.A. 77, 1270–1273. doi: 10.1073/pnas.77.3.1270

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernstrom, J. D., and Fernstrom, M. H. (2007). Tyrosine, phenylalanine, and catecholamine synthesis and function in the brain. J. Nutr. 137, 1539S–1547S. Available online at: http://jn.nutrition.org/content/137/6/1539S.long

PubMed Abstract | Google Scholar

Fitzpatrick, P. F. (1999). Tetrahydropterin-dependent amino acid hydroxylases. Annu. Rev. Biochem. 68, 355–381. doi: 10.1146/annurev.biochem.68.1.355

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaines, C. G., Byng, G. S., Whitaker, R. J., and Jensen, R. A. (1982). L-Tyrosine regulation and biosynthesis via arogenate dehydrogenase in suspension-cultured cells of Nicotiana silvestris Speg. et Comes. Planta 156, 233–240. doi: 10.1007/BF00393730

PubMed Abstract | CrossRef Full Text | Google Scholar

Gamborg, O. L., and Keeley, F. W. (1966). Aromatic metabolism in plants I. A study of the prephenate dehydrogenase from bean plants. Biochim. Biophys. Acta 115, 65–72. doi: 10.1016/0304-4165(66)90049-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Gandía-Herrero, F., and García-Carmona, F. (2013). Biosynthesis of betalains: yellow and violet plant pigments. Trends Plant Sci. 18, 334–343. doi: 10.1016/j.tplants.2013.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Gleadow, R. M., and Møller, B. L. (2014). Cyanogenic glycosides: synthesis, physiology, and phenotypic plasticity. Annu. Rev. Plant Biol. 65, 155–185. doi: 10.1146/annurev-arplant-050213-040027

PubMed Abstract | CrossRef Full Text | Google Scholar

Graindorge, M., Giustini, C., Kraut, A., Moyet, L., Curien, G., and Matringe, M. (2014). Three different classes of aminotransferases evolved prephenate aminotransferase functionality in arogenate-competent microorganisms. J. Biol. Chem. 289, 3198–3208. doi: 10.1074/jbc.M113.486480

PubMed Abstract | CrossRef Full Text | Google Scholar

Gray, M. W., Burger, G., and Lang, B. F. (1999). Mitochondrial evolution. Science 283, 1476–1481. doi: 10.1126/science.283.5407.1476

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, R. S., Mahmood, S., and Adeolu, M. (2013). A phylogenomic and molecular signature based approach for characterization of the phylum Spirochaetes and its major clades: proposal for a taxonomic revision of the phylum. Front. Microbiol. 4:322. doi: 10.3389/fmicb.2013.00322

CrossRef Full Text | Google Scholar

He, Y., Galant, A., Pang, Q., Strul, J. M., Balogun, S. F., Jez, J. M., et al. (2011). Structural and functional evolution of isopropylmalate dehydrogenases in the leucine and glucosinolate pathways of Arabidopsis thaliana. J. Biol. Chem. 286, 28794–28801. doi: 10.1074/jbc.M111.262519

PubMed Abstract | CrossRef Full Text | Google Scholar

Hudson, G. S., Wong, V., and Davidson, B. E. (1984). Chorismate mutase/prephenate dehydrogenase from Escherichia coli K12: purification, characterization, and identification of a reactive cysteine. Biochemistry 23, 6240–6249. doi: 10.1021/bi00320a054

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, R. A., and Pierson, D. L. (1975). Evolutionary implications of different types of microbial enzymology for L-tyrosine biosynthesis. Nature 254, 667–671. doi: 10.1038/254667a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, D. T., Taylor, W. R., and Thornton, J. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275–282. doi: 10.1093/bioinformatics/8.3.275

PubMed Abstract | CrossRef Full Text | Google Scholar

Ku, H., Park, S., Yang, I., and Kim, S. (2010). Expression and functional characterization of prephenate dehydrogenase from Streptococcus mutans. Process Biochem. 45, 607–612. doi: 10.1016/j.procbio.2009.12.006

CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Legrand, P., Dumas, R., Seux, M., Rippert, P., Ravelli, R., Ferrer, J. L., et al. (2006). Biochemical characterization and crystal structure of Synechocystis arogenate dehydrogenase provide insights into catalytic reaction. Structure 14, 767–776. doi: 10.1016/j.str.2006.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Leong, B. J., and Last, R. L. (2017). Promiscuity, impersonation and accommodation: evolution of plant specialized metabolism. Curr. Opin. Struct. Biol. 47, 105–112. doi: 10.1016/j.sbi.2017.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Louie, G. V., Bowman, M. E., Moffitt, M. C., Baiga, T. J., Moore, B. S., and Noel, J. P. (2006). Structural determinants and modulation of substrate specificity in phenylalanine-tyrosine ammonia-lyases. Chem. Biol. 13, 1327–1338. doi: 10.1016/j.chembiol.2006.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeda, H., and Dudareva, N. (2012). The shikimate pathway and aromatic amino acid biosynthesis in plants. Annu. Rev. Plant Biol. 63, 73–105. doi: 10.1146/annurev-arplant-042811-105439

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeda, H., Shasany, A. K., Schnepp, J., Orlova, I., Taguchi, G., Cooper, B. R., et al. (2010). RNAi suppression of Arogenate Dehydratase1 reveals that phenylalanine is synthesized predominantly via the arogenate pathway in petunia petals. Plant Cell 22, 832–849. doi: 10.1105/tpc.109.073247

PubMed Abstract | CrossRef Full Text | Google Scholar

Mannhaupt, G., Stucka, R., Pilz, U., Schwarzlose, C., and Feldmann, H. (1989). Characterization of the prephenate dehydrogenase-encoding gene, TYR1, from Saccharomyces cerevisae. Gene 85, 303–311. doi: 10.1016/0378-1119(89)90422-8

CrossRef Full Text | Google Scholar

Mayer, E., Waldner-Sander, S., Keller, B., Keller, E., and Lingens, F. (1985). Purification of arogenate dehydrogenase from Phenylobacterium immobile. FEBS Lett. 179, 208–212. doi: 10.1016/0014-5793(85)80519-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Millner, P. A., and Barber, J. (1984). Plastoquinone as a mobile redox carrier in the photosynthetic membrane. FEBS Lett. 169, 1–6. doi: 10.1016/0014-5793(84)80277-X

CrossRef Full Text | Google Scholar

Pei, J., and Grishin, N. V. (2007). PROMALS: towards accurate multiple sequence alignments of distantly related proteins. Bioinformatics 23, 802–808. doi: 10.1093/bioinformatics/btm017

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, M. (2013). Rosmarinic acid: new aspects. Phytochem. Rev. 12, 207–227. doi: 10.1007/s11101-013-9282-8

CrossRef Full Text | Google Scholar

Pritt, B. S., Mead, P. S., Johnson, D. K. H., Neitzel, D. F., Respicio-Kingry, L. B., Davis, J. P., et al. (2016). Identification of a novel pathogenic Borrelia species causing Lyme borreliosis with unusually high spirochaetaemia: a descriptive study. Lancet Infect. Dis. 16, 556–564. doi: 10.1016/S1473-3099(15)00464-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Reyes-Prieto, A., and Moustafa, A. (2012). Plastid-localized amino acid biosynthetic pathways of Plantae are predominantly composed of non-cyanobacterial enzymes. Sci. Rep. 2, 955–967. doi: 10.1038/srep00955

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, T. A., Dacks, J. B., Campbell, S. A., Blanchard, J. L., Foster, P. G., McLeod, R., et al. (2006). Evolutionary origins of the eukaryotic shikimate pathway: gene fusions, horizontal gene transfer, and endosymbiotic replacements. Eukaryot. Cell 5, 1517–1531. doi: 10.1128/EC.00106-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Rippert, P., and Matringe, M. (2002). Purification and kinetic analysis of the two recombinant arogenate dehydrogenase isoforms of Arabidopsis thaliana. Eur. J. Biochem. 269, 4753–4761. doi: 10.1046/j.1432-1033.2002.03172.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubin, J. L., and Jensen, R. A. (1979). Enzymology of L-tyrosine biosynthesis in mung bean (Vigna radiata [L.] Wilczek). Plant Physiol. 64, 727–734. doi: 10.1104/pp.64.5.727

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425.

PubMed Abstract | Google Scholar

Sato, F., Inui, T., and Takemura, T. (2007). Metabolic engineering in isoquinoline alkaloid biosynthesis. Curr. Pharm. Biotechnol. 8, 211–218. doi: 10.2174/138920107781387438

PubMed Abstract | CrossRef Full Text | Google Scholar

Schenck, C. A., Chen, S., Siehl, D. L., and Maeda, H. A. (2015). Non-plastidic, tyrosine-insensitive prephenate dehydrogenases from legumes. Nat. Chem. Biol. 11, 52–57. doi: 10.1038/nchembio.1693

PubMed Abstract | CrossRef Full Text | Google Scholar

Schenck, C. A., Holland, C. K., Schneider, M. R., Men, Y., Lee, S. G., Jez, J. M., et al. (2017). Molecular basis of the evolution of alternative tyrosine biosynthetic routes in plants. Nat. Chem. Biol. 13, 1029–1035. doi: 10.1038/nchembio.2414

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi: 10.1038/nmeth.2089

PubMed Abstract | CrossRef Full Text | Google Scholar

Shlaifer, I., Quashie, P. K., Kim, H. Y., and Turnbull, J. L. (2017). Biochemical characterization of TyrA enzymes from Ignicoccus hospitalis and Haemophilus influenzae: a comparative study of the bifunctional and monofunctional dehydrogenase forms. Biochim. Biophys. Acta 1865, 312–320. doi: 10.1016/j.bbapap.2016.12.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Siehl, D. (1999). “The biosynthesis of tryptophan, tyrosine, and phenylalanine from chorismate,” in Plant Amino Acids: Biochemistry and Biotechnology, ed B. Singh (New York, NY: CRC Press), 171–204.

Slominski, A., Tobin, D. J., Shibahara, S., and Wortsman, J. (2004). Melanin pigmentation in mammalian skin and its hormonal regulation. Physiol. Rev. 84, 1155–1228. doi: 10.1152/physrev.00044.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramaniam, P., Bhatnagar, R., Hooper, A., and Jensen, R. A. (1994). The dynamic progression of evolved character states for aromatic amino acid biosynthesis in gram-negative bacteria. Microbiology 140, 3431–3440. doi: 10.1099/13500872-140-12-3431

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, W., Singh, S., Zhang, R., Turnbull, J. L., and Christendat, D. (2006). Crystal structure of prephenate dehydrogenase from Aquifex aeolicus: insights into the catalytic mechanism. J. Biol. Chem. 281, 12919–12928. doi: 10.1074/jbc.M511986200

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, T. H., and Jensen, R. A. (1990). A single clycohexadienyl dehydrogenase specifies the prephenate dehydrogenase and arogenate dehydrogenase components of the dual pathways to L-tyrosine in Pseudomonas aeruginosa. J. Biol. Chem. 265, 20033–20036.

PubMed Abstract | Google Scholar

Zhao, G., Xia, T., Ingram, L. O., and Jensen, R. A. (1993). An allosterically insensitive class of cyclohexadienyl dehydrogenase from Zymomonas mobilis. Eur. J. Biochem. 212, 157–165. doi: 10.1111/j.1432-1033.1993.tb17646.x

PubMed Abstract | CrossRef Full Text | Google Scholar