Folates are tripartite molecules comprising pterin, para-aminobenzoate (PABA), and glutamate moieties, which are essential cofactors involved in DNA and amino acid synthesis. The obligately intracellular Chlamydia species have lost several biosynthetic pathways for essential nutrients which they can obtain from their host but have retained the capacity to synthesize folate. In most bacteria, synthesis of the pterin moiety of folate requires the FolEQBK enzymes, while synthesis of the PABA moiety is carried out by the PabABC enzymes. Bioinformatic analyses reveal that while members of Chlamydia are missing the genes for FolE (GTP cyclohydrolase) and FolQ, which catalyze the initial steps in de novo synthesis of the pterin moiety, they have genes for the rest of the pterin pathway. We screened a chlamydial genomic library in deletion mutants of Escherichia coli to identify the “missing genes” and identified a novel enzyme, TrpFCtL2, which has broad substrate specificity. TrpFCtL2, in combination with GTP cyclohydrolase II (RibA), the first enzyme of riboflavin synthesis, provides a bypass of the first two canonical steps in folate synthesis catalyzed by FolE and FolQ. Notably, TrpFCtL2 retains the phosphoribosyl anthranilate isomerase activity of the original annotation. Additionally, we independently confirmed the recent discovery of a novel enzyme, CT610, which uses an unknown precursor to synthesize PABA and complements E. coli mutants with deletions of pabA, pabB, or pabC. Thus, Chlamydia species have evolved a variant folate synthesis pathway that employs a patchwork of promiscuous and adaptable enzymes recruited from other biosynthetic pathways.

BackgroundCurrent sequence-based approaches to identify enzyme functional shifts, such as enzyme promiscuity, have proven to be highly dependent on a priori functional knowledge, hampering our ability to reconstruct evolutionary history behind these mechanisms. Hidden Markov Model (HMM) profiles, broadly used to classify enzyme families, can be useful to distinguish between closely related enzyme families with different specificities. The (βα)8-isomerase HisA/PriA enzyme family, involved in L-histidine (HisA, mono-substrate) biosynthesis in most bacteria and plants, but also in L-tryptophan (HisA/TrpF or PriA, dual-substrate) biosynthesis in most Actinobacteria, has been used as model system to explore evolutionary hypotheses and therefore has a considerable amount of evolutionary, functional and structural knowledge available. We searched for functional evolutionary intermediates between the HisA and PriA enzyme families in order to understand the functional divergence between these families.ResultsWe constructed a HMM profile that correctly classifies sequences of unknown function into the HisA and PriA enzyme sub-families. Using this HMM profile, we mined a large metagenome to identify plausible evolutionary intermediate sequences between HisA and PriA. These sequences were used to perform phylogenetic reconstructions and to identify functionally conserved amino acids. Biochemical characterization of one selected enzyme (CAM1) with a mutation within the functionally essential N-terminus phosphate-binding site, namely, an alanine instead of a glycine in HisA or a serine in PriA, showed that this evolutionary intermediate has dual-substrate specificity. Moreover, site-directed mutagenesis of this alanine residue, either backwards into a glycine or forward into a serine, revealed the robustness of this enzyme. None of these mutations, presumably upon functionally essential amino acids, significantly abolished its enzyme activities. A truncated version of this enzyme (CAM2) predicted to adopt a (βα)6-fold, and thus entirely lacking a C-terminus phosphate-binding site, was identified and shown to have HisA activity.ConclusionAs expected, reconstruction of the evolution of PriA from HisA with HMM profiles suggest that functional shifts involve mutations in evolutionarily intermediate enzymes of otherwise functionally essential residues or motifs. These results are in agreement with a link between promiscuous enzymes and intragenic epistasis. HMM provides a convenient approach for gaining insights into these evolutionary processes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0378-1) contains supplementary material, which is available to authorized users.

The connection between gene loss and the functional adaptation of retained proteins is still poorly understood. We apply phylogenomics and metabolic modeling to detect bacterial species that are evolving by gene loss, with the finding that Actinomycetaceae genomes from human cavities are undergoing sizable reductions, including loss of L-histidine and L-tryptophan biosynthesis. We observe that the dual-substrate phosphoribosyl isomerase A or priA gene, at which these pathways converge, appears to coevolve with the occurrence of trp and his genes. Characterization of a dozen PriA homologs shows that these enzymes adapt from bifunctionality in the largest genomes, to a monofunctional, yet not necessarily specialized, inefficient form in genomes undergoing reduction. These functional changes are accomplished via mutations, which result from relaxation of purifying selection, in residues structurally mapped after sequence and X-ray structural analyses. Our results show how gene loss can drive the evolution of substrate specificity from retained enzymes.DOI: http://dx.doi.org/10.7554/eLife.22679.001

We investigate the evolution of co-occurring analogous enzymes involved in L-tryptophan and L-histidine biosynthesis in Actinobacteria Phylogenetic analysis of trpF homologues, a missing gene in certain clades of this lineage whose absence is complemented by a dual-substrate HisA homologue, termed PriA, found that they fall into three categories: (i) trpF-1, an L-tryptophan biosynthetic gene horizontally acquired by certain Corynebacterium species; (ii) trpF-2, a paralogue known to be involved in synthesizing a pyrrolopyrrole moiety and (iii) trpF-3, a variable non-conserved orthologue of trpF-1 We previously investigated the effect of trpF-1 upon the evolution of PriA substrate specificity, but nothing is known about the relationship between trpF-3 and priA After in vitro steady-state enzyme kinetics we found that trpF-3 encodes a phosphoribosyl anthranilate isomerase. However, mutation of this gene in Streptomyces sviceus did not lead to auxothrophy, as expected from the biosynthetic role of trpF-1 Biochemical characterization of a dozen co-occurring TrpF-2 or TrpF-3, with PriA homologues, explained the prototrophic phenotype, and unveiled an enzyme activity trade-off between TrpF and PriA. X-ray structural analysis suggests that the function of these PriA homologues is mediated by non-conserved mutations in the flexible L5 loop, which may be responsible for different substrate affinities. Thus, the PriA homologues that co-occur with TrpF-3 represent a novel enzyme family, termed PriB, which evolved in response to PRA isomerase activity. The characterization of co-occurring enzymes provides insights into the influence of functional redundancy on the evolution of enzyme function, which could be useful for enzyme functional annotation.

Neither the Pseudomonas aeruginosa aldehyde dehydrogenase encoded by the PA4189 gene nor its ortholog proteins have been biochemically or structurally characterized and their physiological function is unknown. We cloned the PA4189 gene, obtained the PA4189 recombinant protein, and studied its structure-function relationships. PA4189 is an NAD+-dependent aminoaldehyde dehydrogenase highly efficient with protonated aminoacetaldehyde and 3-aminopropionaldehyde, which are much more preferred to the non-protonated species as indicated by pH studies. Based on the higher activity with aminoacetaldehyde than with 3-aminopropionaldehyde, we propose that aminoacetaldehyde might be the PA4189 physiological substrate. Even though at the physiological pH of P. aeruginosa cells the non-protonated aminoacetaldehyde species will be predominant, and despite the competition of these species with the protonated ones, PA4189 would very efficiently oxidize ACTAL in vivo, producing glycine. To our knowledge, PA4189 is the first reported enzyme that might metabolize ACTAL, which is considered a dead-end metabolite because its consuming reactions are unknown. The PA4189 crystal structure reported here suggested that the charge and size of the active-site residue Glu457, which narrows the aldehyde-entrance tunnel, greatly define the specificity for small positively charged aldehydes, as confirmed by the kinetics of the E457G and E457Q variants. Glu457 and the residues that determine Glu457 conformation inside the active site are conserved in the PA4189 orthologs, which we only found in proteobacteria species. Also is conserved the PA4189 genomic neighborhood, which suggests that PA4189 participates in an uncharacterized metabolic pathway. Our results open the door to future efforts to characterize this pathway.