Promiscuous and adaptable enzymes fill "holes" in the tetrahydrofolate pathway in Chlamydia species

The Lreu_1276 protein from Limosilactobacillus reuteri represents a third family of dihydroneopterin triphosphate pyrophosphohydrolases in bacteria

Tetrahydrofolate is a cofactor involved in C1 metabolism including biosynthesis pathways for adenine and serine. In the classical tetrahydrofolate biosynthesis pathway, the steps removing three phosphate groups from the precursor 7,8-dihydroneopterin triphosphate (DHNTP) remain unclear in many bacteria. DHNTP pyrophosphohydrolase hydrolyzes pyrophosphate from DHNTP and produces 7,8-dihydroneopterin monophosphate. Although two structurally distinct DHNTP pyrophosphohydrolases have been identified in the intestinal bacteria Lactococcus lactis and Escherichia coli, the distribution of their homologs is limited. Here, we aimed to identify a third DHNTP pyrophosphohydrolase gene in the intestinal lactic acid bacterium Limosilactobacillus reuteri. In a gene operon including genes involved in dihydrofolate biosynthesis, we focused on the lreu_1276 gene, annotated as Ham1 family protein or XTP/dITP diphosphohydrolase, as a candidate encoding DHNTP pyrophosphohydrolase. The Lreu_1276 recombinant protein was prepared using E. coli and purified. Biochemical analyses of the reaction product revealed that the Lreu_1276 protein displays significant pyrophosphohydrolase activity toward DHNTP. The optimal reaction temperature and pH were 35°C and around 7, respectively. Substrate specificity was relatively strict among 17 tested compounds. Although previously characterized DHNTP pyrophosphohydrolases prefer Mg2+, the Lreu_1276 protein exhibited maximum activity in the presence of Mn2+, with a specific activity of 28.2 ± 2.0 µmol min-1 mg-1 in the presence of 1 mM Mn2+. The three DHNTP pyrophosphohydrolases do not share structural similarity to one another, and the distribution of their homologs does not overlap, implying that the Lreu_1276 protein represents a third structurally novel DHNTP pyrophosphohydrolase in bacteria.

IMPORTANCE

The identification of a structurally novel DHNTP pyrophosphohydrolase in L. reuteri provides valuable information in understanding tetrahydrofolate biosynthesis in bacteria that possess lreu_1276 homologs. Interestingly, however, even with the identification of a third family of DHNTP pyrophosphohydrolases, there are still a number of bacteria that do not harbor homologs for any of the three genes while possessing other genes involved in the biosynthesis of the pterin ring structure. This suggests the presence of an unrecognized DHNTP pyrophosphohydrolase gene in bacteria. As humans do not harbor DHNTP pyrophosphohydrolase, the high structural diversity of enzymes responsible for a reaction in tetrahydrofolate biosynthesis may provide an advantage in designing inhibitors targeting a specific group of bacteria in the intestinal microbiota.

Metabolism and physiology of pathogenic bacterial obligate intracellular parasites

Bacterial obligate intracellular parasites (BOIPs) represent an exclusive group of bacterial pathogens that all depend on invasion of a eukaryotic host cell to reproduce. BOIPs are characterized by extensive adaptation to their respective replication niches, regardless of whether they replicate within the host cell cytoplasm or within specialized replication vacuoles. Genome reduction is also a hallmark of BOIPs that likely reflects streamlining of metabolic processes to reduce the need for de novo biosynthesis of energetically costly metabolic intermediates. Despite shared characteristics in lifestyle, BOIPs show considerable diversity in nutrient requirements, metabolic capabilities, and general physiology. In this review, we compare metabolic and physiological processes of prominent pathogenic BOIPs with special emphasis on carbon, energy, and amino acid metabolism. Recent advances are discussed in the context of historical views and opportunities for discovery.

Temporal and spatial differences in the vaginal microbiome of Chinese healthy women

Background

Up the reproductive tract, there are large differences in the composition of vaginal microbes. Throughout the menstrual cycle, the structure of the vaginal microbiome shifts. Few studies have examined both in combination. Our study was designed to explore trends in the microbiome of different parts of the vagina in healthy women over the menstrual cycle.

Methods

We performed metagenomic sequencing to characterize the microbiome differences between the cervical orifice and mid-vagina throughout the menstrual cycle.

Results

Our results showed the vaginal microbiome of healthy women in the cervical orifice and the mid-vagina was similar during the periovulatory and luteal phases, with Lactobacillus being the dominant bacteria. In the follicular phase, Acinetobacter was detected in the cervical orifice. From the follicular phase to the luteal phase, the community state types (all five community status types were defined as CSTs) in samples No. 10 and No. 11 changed from CST III to CST I. In addition, the composition of the vaginal microbiome in healthy women from different regions of China was significantly different. We also detected viruses including Human alphaherpesvirus 1 (HSV-1) during periovulatory phase.

Conclusion

This study is valuable for understanding whether the microbial composition of the vagina is consistent in different parts of the menstrual cycle.

Excision of a Protein-Derived Amine for p-Aminobenzoate Assembly by the Self-Sacrificial Heterobimetallic Protein CADD

Chlamydia protein associating with death domains (CADD), the founding member of a recently discovered class of nonheme dimetal enzymes termed hemeoxygenase-like dimetaloxidases (HDOs), plays an indispensable role in pathogen survival. CADD orchestrates the biosynthesis of p-aminobenzoic acid (pABA) for integration into folate via the self-sacrificial excision of a protein-derived tyrosine (Tyr27) and several additional processing steps, the nature and timing of which have yet to be fully clarified. Nuclear magnetic resonance (NMR) and proteomics approaches reveal the source and probable timing of amine installation by a neighboring lysine (Lys152). Turnover studies using limiting O2 have identified a para-aminobenzaldehyde (pABCHO) metabolic intermediate that is formed on the path to pABA formation. The use of pABCHO and other probe substrates shows that the heterobimetallic Fe/Mn form of the enzyme is capable of oxygen insertion to generate the pABA-carboxylate.

Streptomyces lividans 66 produces a protease inhibitor via a tRNA-utilizing enzyme interacting with a C-minus NRPS

Small peptide aldehydes (SPAs) with protease inhibitory activity are naturally occurring compounds shown to be synthesized by non-ribosomal peptide synthetases (NRPS). SPAs are widely used in biotechnology and have been utilized as therapeutic agents. They are also physiologically relevant and have been postulated to regulate the development of their producing microorganisms. Previously, we identified an NRPS-like biosynthetic gene cluster (BGC) in Streptomyces lividans 66 that lacked a condensation (C) domain but included a tRNA-utilizing enzyme (tRUE) belonging to the leucyl/phenylalanyl (L/F) transferase family. This system was predicted to direct the synthesis of a novel SPA, which we named livipeptin. Using evolutionary genome mining approaches, here, we confirm the presence of L/F transferase tRUEs within the genomes of diverse Streptomyces and related organisms, including fusions with the anticipated C-minus NRPS-like protein. We then demonstrate genetic functional cooperation between the identified L/F-transferase divergent tRUE homolog with the C-minus NRPS, leading to the synthesis of a metabolic fraction with protease inhibitory activity. Semisynthetic assays in the presence of RNAse revealed that the productive interaction between the tRUE and the C-minus NRPS enzymes is indeed tRNA dependent. We expect our findings to boost the discovery of SPAs, as well as the development of protease-mediated biotechnologies, by exploiting the uncovered genetic basis for synthesizing putative acetyl-leu/phe-arginine protease inhibitors. Furthermore, these results will facilitate the purification and structural elucidation of livipeptin, which has proven difficult to chemically characterize.

SIGNIFICANCE

The discovery of natural products biosynthetic genes marks a significant advancement in our understanding of these metabolites, for example of their evolution, activity, and biosynthesis, but also opens biotechnological opportunities and knowledge to advance genome mining approaches. We made this possible by uncovering a new biosynthetic pathway in Streptomyces lividans 66 shown to direct the synthesis of a strong protease inhibitor, termed livipeptin, following unprecedented biosynthetic rules and genes. Thus, by shedding light on the genetic mechanisms predicted to govern the production of acetyl-leu/phe-arginine protease inhibitors, including the elusive livipeptin, this study enables novel protease-mediated biotechnologies as well as approaches for discovering protease inhibitors from genome data.

The Chlamydia trachomatis p-aminobenzoate synthase CADD is a manganese-dependent oxygenase that uses its own amino acid residues as substrates

CADD (chlamydia protein associating with death domains) is a p-aminobenzoate (pAB) synthase involved in a noncanonical route for tetrahydrofolate biosynthesis in Chlamydia trachomatis. Although previously implicated to employ a diiron cofactor, here, we show that pAB synthesis by CADD requires manganese and the physiological cofactor is most likely a heterodinuclear Mn/Fe cluster. Isotope-labeling experiments revealed that the two oxygen atoms in the carboxylic acid portion of pAB are derived from molecular oxygen. Further, mass spectrometry-based proteomic analyses of CADD-derived peptides demonstrated a glycine substitution at Tyr27, providing strong evidence that this residue is sacrificed for pAB synthesis. Additionally, Lys152 was deaminated and oxidized to aminoadipic acid, supporting its proposed role as a sacrificial amino group donor.

Host Cell Amplification of Nutritional Stress Contributes To Persistence in Chlamydia trachomatis

Persistence, a viable but non-replicating growth state, has been implicated in diseases caused by Chlamydia trachomatis. Starvation of distinct nutrients produces a superficially similar persistent state, implying convergence on a common intracellular environment. We employed host-pathogen dual RNA-sequencing under both iron- and tryptophan-starved conditions to systematically characterize the persistent chlamydial transcriptome and to define common contributions of the host cell transcriptional stress response in shaping the intracellular environment. The transcriptome of the infected host cells was highly specific to each nutritional stress, despite comparable effects on chlamydial growth and development in each condition. In contrast, the chlamydial transcriptomes between nutritional conditions were highly similar, suggesting some overlap in host cell responses to iron limitation and tryptophan starvation that contribute to a common persistent phenotype. We demonstrate that a commonality in the host cell responses is the suppression of GTP biosynthesis, a nucleotide for which Chlamydia are auxotrophic. Pharmacological inhibition of host IMP dehydrogenase (IMPDH1), which catalyzes the rate-limiting step in de novo guanine nucleotide synthesis, resulted in comparable GTP depletion to both iron and tryptophan starvation and induced chlamydial persistence. Moreover, IMPDH1 inhibition and iron starvation acted synergistically to control chlamydial growth. Thus, host cell reduction in GTP levels amplifies the nutritional stress to intracellular chlamydiae in infection-relevant models of persistence, illustrating the determinative role the infected host cell plays in bacterial stress responses. IMPORTANCE Bacteria respond to nutritional stress through universal and unique mechanisms. Genome reduction in the Chlamydiaceae, a consequence of coevolution with their obligate eukaryotic hosts, has reduced their repertoire of stress response mechanisms. Here, we demonstrate that the infected host cell may provide the context within which universal stress responses emerge for Chlamydia trachomatis. We report that during starvation of the essential nutrients iron or tryptophan, a common response of the infected epithelial cell is the suppression of GTP biosynthesis, which induces a persistent developmental state in the pathogen. Thus, chlamydial persistence results from the combined effects of primary stresses on the pathogen and the host, with the latter eliciting a secondary host cell response that intensifies the inhospitable intracellular environment.

In-Silico Functional Metabolic Pathways Associated to Chlamydia trachomatis Genital Infection

The advent of high-throughput technologies, such as 16s rDNA sequencing, has significantly contributed to expanding our knowledge of the microbiota composition of the genital tract during infections such as Chlamydia trachomatis. The growing body of metagenomic data can be further exploited to provide a functional characterization of microbial communities via several powerful computational approaches. Therefore, in this study, we investigated the predicted metabolic pathways of the cervicovaginal microbiota associated with C. trachomatis genital infection in relation to the different Community State Types (CSTs), via PICRUSt2 analysis. Our results showed a more rich and diverse mix of predicted metabolic pathways in women with a CST-IV microbiota as compared to all the other CSTs, independently from infection status. C. trachomatis genital infection further modified the metabolic profiles in women with a CST-IV microbiota and was characterized by increased prevalence of the pathways for the biosynthesis of precursor metabolites and energy, biogenic amino-acids, nucleotides, and tetrahydrofolate. Overall, predicted metabolic pathways might represent the starting point for more precisely designed future metabolomic studies, aiming to investigate the actual metabolic pathways characterizing C. trachomatis genital infection in the cervicovaginal microenvironment.

Self-sacrificial tyrosine cleavage by an Fe:Mn oxygenase for the biosynthesis of para-aminobenzoate in Chlamydia trachomatis

Chlamydia protein associating with death domains (CADD) is involved in the biosynthesis of para-aminobenzoate (pABA), an essential component of the folate cofactor that is required for the survival and proliferation of the human pathogen Chlamydia trachomatis. The pathway used by Chlamydiae for pABA synthesis differs from the canonical multi-enzyme pathway used by most bacteria that relies on chorismate as a metabolic precursor. Rather, recent work showed pABA formation by CADD derives from l-tyrosine. As a member of the emerging superfamily of heme oxygenase-like diiron oxidases (HDOs), CADD was proposed to use a diiron cofactor for catalysis. However, we report maximal pABA formation by CADD occurs upon the addition of both iron and manganese, which implicates a heterobimetallic Fe:Mn cluster is the catalytically active form. Isotopic labeling experiments and proteomics studies show that CADD generates pABA from a protein-derived tyrosine (Tyr27), a residue that is ∼14 Å from the dimetal site. We propose that this self-sacrificial reaction occurs through O₂ activation by a probable Fe:Mn cluster through a radical relay mechanism that connects to the "substrate" Tyr, followed by amination and direct oxygen insertion. These results provide the molecular basis for pABA formation in C. trachomatis, which will inform the design of novel therapeutics.

Dioxygen Activation and Nδ,Nε-Dihydroxylation Mechanism Involved in the Formation of N-Nitrosourea Pharmacophore in Streptozotocin Catalyzed by Nonheme Diiron Enzyme SznF

SznF is a nonheme diiron-dependent enzyme that catalyzes the critical N-nitrosation involved in the formation of the N-nitrosourea moiety in the pancreatic cancer drug streptozotocin. The N-nitrosation contains two successive N-hydroxylation and N-nitrosation steps, which are carried out by two separate active sites, namely, the central domain and cupin domain. Recently, the crystal structure of SznF was obtained, and the central domain was proved to contain a diiron cofactor to catalyze the N-hydroxylation. In this work, to gain insights into the O₂ activation and the successive N-hydroxylation mechanism, on the basis of the high-resolution crystal structure, the enzyme-substrate complex models were constructed, and a series of combined QM/MM calculations were performed. Based on our calculations, the activation of O₂ starts from the diiron(II,III)-superoxo (S) to generate the diiron(IV)-oxo species (Q) via a diiron(III,III)-peroxo (P)-like transition state or unstable intermediate (P'), and species P' can be described as a hybridization of diiron(IV)-oxo species and diiron(III,III)-peroxo (P) owing to the long distances of Fe1-Fe2 (4.22 Å) and O1-O2 (1.89 Å), which is different from those of other nonheme diiron enzymes. In the following hydroxylation of N^δ and N^ε, the N^δ-hydroxylation was confirmed to occur first, agreeing with the experimental observations. Because the diiron(IV)-oxo species (Q) is responsible for hydroxylation, the reaction follows the H-abstraction/OH rebound mechanism, and the first abstraction occurs on the N^δ-H rather than N^ε-H, which may be attributed to the different orientation of Fe(IV)-oxo relative to N-H as well as the bond dissociation enthalpies of two N-H bonds. The hydroxylation of N-methyl-L-arginine does not employ the diiron(III,III)-hydroperoxo (P″) to trigger the electrophilic attack of the guanidine to directly form the N-O bond, as previously suggested. In addition, our calculations also revealed that the direct attack of the Fe(IV)═O unit to the N^δ of the substrate corresponds to a higher barrier than that in the H-abstraction/OH rebound mechanism. These results may provide useful information for understanding the formation of the di-hydroxylation intermediate involved in the biosynthesis of N-nitrosation.

BesC Initiates C-C Cleavage through a Substrate-Triggered and Reactive Diferric-Peroxo Intermediate

BesC catalyzes the iron- and O₂-dependent cleavage of 4-chloro-l-lysine to form 4-chloro-l-allylglycine, formaldehyde, and ammonia. This process is a critical step for a biosynthetic pathway that generates a terminal alkyne amino acid which can be leveraged as a useful bio-orthogonal handle for protein labeling. As a member of an emerging family of diiron enzymes that are typified by their heme oxygenase-like fold and a very similar set of coordinating ligands, recently termed HDOs, BesC performs an unusual type of carbon-carbon cleavage reaction that is a significant departure from reactions catalyzed by canonical dinuclear-iron enzymes. Here, we show that BesC activates O₂ in a substrate-gated manner to generate a diferric-peroxo intermediate. Examination of the reactivity of the peroxo intermediate with a series of lysine derivatives demonstrates that BesC initiates this unique reaction trajectory via cleavage of the C4-H bond; this process represents the rate-limiting step in a single turnover reaction. The observed reactivity of BesC represents the first example of a dinuclear-iron enzyme that utilizes a diferric-peroxo intermediate to capably cleave a C-H bond as part of its native function, thus circumventing the formation of a high-valent intermediate more commonly associated with substrate monooxygenations.

Horizontal-Acquisition of a Promiscuous Peptidoglycan-Recycling Enzyme Enables Aphids To Influence Symbiont Cell Wall Metabolism

During evolution, enzymes can undergo shifts in preferred substrates or in catalytic activities. An intriguing question is how enzyme function changes following horizontal gene transfer, especially for bacterial genes that have moved to animal genomes. Some insects have acquired genes that encode enzymes for the biosynthesis of bacterial cell wall components and that appear to function to support or control their obligate endosymbiotic bacteria. In aphids, the bacterial endosymbiont Buchnera aphidicola provides essential amino acids for aphid hosts but lacks most genes for remodeling of the bacterial cell wall. The aphid genome has acquired seven genes with putative functions in cell wall metabolism that are primarily expressed in the aphid cells harboring Buchnera. In analyses of aphid homogenates, we detected peptidoglycan (PGN) muropeptides indicative of the reactions of PGN hydrolases encoded by horizontally acquired aphid genes but not by Buchnera genes. We produced one such host enzyme, ApLdcA, and characterized its activity with both cell wall derived and synthetic PGN. Both ApLdcA and the homologous enzyme in Escherichia coli, which functions as an l,d-carboxypeptidase in the cytoplasmic PGN recycling pathway, exhibit turnover of PGN substrates containing stem pentapeptides and cross-linkages via l,d-endopeptidase activity, consistent with a potential role in cell wall remodeling. Our results suggest that ApLdcA derives its functions from the promiscuous activities of an ancestral LdcA enzyme, whose acquisition by the aphid genome may have enabled hosts to influence Buchnera cell wall metabolism as a means to control symbiont growth and division. IMPORTANCE Most enzymes are capable of performing biologically irrelevant side reactions. During evolution, promiscuous enzyme activities may acquire new biological roles, especially after horizontal gene transfer to new organisms. Pea aphids harbor obligate bacterial symbionts called Buchnera and encode horizontally acquired bacterial genes with putative roles in cell wall metabolism. Though Buchnera lacks cell wall endopeptidase genes, we found evidence of endopeptidase activity among peptidoglycan muropeptides purified from aphids. We characterized a multifunctional, aphid-encoded enzyme, ApLdcA, which displays l,d-endopeptidase activities considered promiscuous for the Escherichia coli homolog, for which these activities do not contribute to its native role in peptidoglycan recycling. These results exemplify the roles of enzyme promiscuity and horizontal gene transfer in enzyme evolution and demonstrate how aphids influence symbiont cell wall metabolism.

Comparative Genomic Analysis of the DUF34 Protein Family Suggests Role as a Metal Ion Chaperone or Insertase

Members of the DUF34 (domain of unknown function 34) family, also known as the NIF3 protein superfamily, are ubiquitous across superkingdoms. Proteins of this family have been widely annotated as "GTP cyclohydrolase I type 2" through electronic propagation based on one study. Here, the annotation status of this protein family was examined through a comprehensive literature review and integrative bioinformatic analyses that revealed varied pleiotropic associations and phenotypes. This analysis combined with functional complementation studies strongly challenges the current annotation and suggests that DUF34 family members may serve as metal ion insertases, chaperones, or metallocofactor maturases. This general molecular function could explain how DUF34 subgroups participate in highly diversified pathways such as cell differentiation, metal ion homeostasis, pathogen virulence, redox, and universal stress responses.

Structure and assembly of the diiron cofactor in the heme-oxygenase-like domain of the N-nitrosourea-producing enzyme SznF

The enzyme SznF assembles the N-nitrosourea pharmacophore of the drug streptozotocin. Its N-oxygenase domain resembles heme-oxygenase (HO) and belongs to an emerging superfamily of HO-like diiron enzymes (HDOs) with unstable metallocofactors that have resisted structural characterization. Insight into cofactor dynamics from our prior investigation of SznF’s N-oxygenation reactions suggested an approach that has yielded a structure of a functionally assigned HDO with its diiron cofactor intact. Conformational changes accompanying cofactor (dis)assembly explain its instability, and the observation of an unanticipated glutamate ligand that is conserved in only a subset of HDO sequences provides a potential basis for top-level assignment of enzymatic function. Our results thus provide a roadmap for structural and functional characterization of novel HDOs.

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates N^δ and N^ω′ of N^ω-methyl-l-arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme’s C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase–like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase–like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O₂-reactive Fe₂(II/II) complex and the structural basis for cofactor instability—a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon–hydrogen and carbon–carbon bonds to install olefins. Among ∼9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ∼25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

An Unusual Route for p-Aminobenzoate Biosynthesis in Chlamydia trachomatis Involves a Probable Self-Sacrificing Diiron Oxygenase

Chlamydia trachomatis lacks the canonical genes required for the biosynthesis of p-aminobenzoate (pABA), a component of essential folate cofactors. Previous studies revealed a single gene from C. trachomatis, the CT610 gene, that rescues Escherichia coli ΔpabA, ΔpabB, and ΔpabC mutants, which are otherwise auxotrophic for pABA. CT610 shares low sequence similarity to nonheme diiron oxygenases, and the previously solved crystal structure revealed a diiron active site. Genetic studies ruled out several potential substrates for CT610-dependent pABA biosynthesis, including chorismate and other shikimate pathway intermediates, leaving the actual precursor(s) unknown. Here, we supplied isotopically labeled potential precursors to E. coli ΔpabA cells expressing CT610 and found that the aromatic portion of tyrosine was highly incorporated into pABA, indicating that tyrosine is a precursor for CT610-dependent pABA biosynthesis. Additionally, in vitro enzymatic experiments revealed that purified CT610 exhibits low pABA synthesis activity under aerobic conditions in the absence of tyrosine or other potential substrates, where only the addition of a reducing agent such as dithiothreitol appears to stimulate pABA production. Furthermore, site-directed mutagenesis studies revealed that two conserved active site tyrosine residues are essential for the pABA synthesis reaction in vitro Thus, the current data are most consistent with CT610 being a unique self-sacrificing enzyme that utilizes its own active site tyrosine residue(s) for pABA biosynthesis in a reaction that requires O₂ and a reduced diiron cofactor.IMPORTANCE Chlamydia trachomatis is the most reported sexually transmitted infection in the United States and the leading cause of infectious blindness worldwide. Unlike many other intracellular pathogens that have undergone reductive evolution, C. trachomatis is capable of de novo biosynthesis of the essential cofactor tetrahydrofolate using a noncanonical pathway. Here, we identify the biosynthetic precursor to the p-aminobenzoate (pABA) portion of folate in a process that requires the CT610 enzyme from C. trachomatis We further provide evidence that CT610 is a self-sacrificing or "suicide" enzyme that uses its own amino acid residue(s) as the substrate for pABA synthesis. This work provides the foundation for future investigation of this chlamydial pABA synthase, which could lead to new therapeutic strategies for C. trachomatis infections.

A Peroxodiiron(III/III) Intermediate Mediating Both N-Hydroxylation Steps in Biosynthesis of the N-Nitrosourea Pharmacophore of Streptozotocin by the Multi-domain Metalloenzyme SznF

The alkylating warhead of the pancreatic cancer drug streptozotocin (SZN) contains an N-nitrosourea moiety constructed from N^ω-methyl-l-arginine (l-NMA) by the multi-domain metalloenzyme SznF. The enzyme's central heme-oxygenase-like (HO-like) domain sequentially hydroxylates N^δ and N^ω' of l-NMA. Its C-terminal cupin domain then rearranges the triply modified arginine to N^δ-hydroxy-N^ω-methyl-N^ω-nitroso-l-citrulline, the proposed donor of the functional pharmacophore. Here we show that the HO-like domain of SznF can bind Fe(II) and use it to capture O₂, forming a peroxo-Fe₂(III/III) intermediate. This intermediate has absorption- and Mössbauer-spectroscopic features similar to those of complexes previously trapped in ferritin-like diiron oxidases and oxygenases (FDOs) and, more recently, the HO-like fatty acid oxidase UndA. The SznF peroxo-Fe₂(III/III) complex is an intermediate in both hydroxylation steps, as shown by the concentration-dependent acceleration of its decay upon exposure to either l-NMA or N^δ-hydroxy-N^ω-methyl-l-Arg (l-HMA). The Fe₂(III/III) cluster produced upon decay of the intermediate has a small Mössbauer quadrupole splitting parameter, implying that, unlike the corresponding product states of many FDOs, it lacks an oxo-bridge. The subsequent decomposition of the product cluster to one or more paramagnetic Fe(III) species over several hours explains why SznF was previously purified and crystallographically characterized without its cofactor. Programmed instability of the oxidized form of the cofactor appears to be a unifying characteristic of the emerging superfamily of HO-like diiron oxidases and oxygenases (HDOs).

Chlamydia trachomatis Oligopeptide Transporter Performs Dual Functions of Oligopeptide Transport and Peptidoglycan Recycling

Peptidoglycan, the sugar-amino acid polymer that composes the bacterial cell wall, requires a significant expenditure of energy to synthesize and is highly immunogenic. To minimize the loss of an energetically expensive metabolite and avoid host detection, bacteria often recycle their peptidoglycan, transporting its components back into the cytoplasm, where they can be used for subsequent rounds of new synthesis. The peptidoglycan-recycling substrate binding protein (SBP) MppA, which is responsible for recycling peptidoglycan fragments in Escherichia coli, has not been annotated for most intracellular pathogens. One such pathogen, Chlamydia trachomatis, has a limited capacity to synthesize amino acids de novo and therefore must obtain oligopeptides from its host cell for growth. Bioinformatics analysis suggests that the putative C. trachomatis oligopeptide transporter OppABCDF (OppABCDF _Ct ) encodes multiple SBPs (OppA1 _Ct , OppA2 _Ct , and OppA3 _Ct ). Intracellular pathogens often encode multiple SBPs, while only one, OppA, is encoded in the E. coli opp operon. We hypothesized that the putative OppABCDF transporter of C. trachomatis functions in both oligopeptide transport and peptidoglycan recycling. We coexpressed the putative SBP genes (oppA1_Ct , oppA2_Ct , oppA3_Ct ) along with oppBCDF_Ct in an E. coli mutant lacking the Opp transporter and determined that all three chlamydial OppA subunits supported oligopeptide transport. We also demonstrated the in vivo functionality of the chlamydial Opp transporter in C. trachomatis Importantly, we found that one chlamydial SBP, OppA3 _Ct , possessed dual substrate recognition properties and is capable of transporting peptidoglycan fragments (tri-diaminopimelic acid) in E. coli and in C. trachomatis These findings suggest that Chlamydia evolved an oligopeptide transporter to facilitate the acquisition of oligopeptides for growth while simultaneously reducing the accumulation of immunostimulatory peptidoglycan fragments in the host cell cytosol. The latter property reflects bacterial pathoadaptation that dampens the host innate immune response to Chlamydia infection.

How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation

Promiscuity is the coincidental ability of an enzyme to catalyze its native reaction and additional reactions that are not biological functions in the same active site. Promiscuity plays a central role in enzyme evolution and is thus a useful property for protein and metabolic engineering. This review examines enzyme evolution holistically, beginning with evaluating biochemical support for four enzyme evolution models. As expected, there is strong biochemical support for the subfunctionalization and innovation-amplification-divergence models, in which promiscuity is a central feature. In many cases, however, enzyme evolution is more complex than the models indicate, suggesting much is yet to be learned about selective pressures on enzyme function. A complete understanding of enzyme evolution must also explain the ability of metabolic networks to integrate new enzyme activities. Hidden within metabolic networks are underground metabolic pathways constructed from promiscuous activities. We discuss efforts to determine the diversity and pervasiveness of underground metabolism. Remarkably, several studies have discovered that some metabolic defects can be repaired via multiple underground routes. In prokaryotes, metabolic innovation is driven by connecting enzymes acquired by horizontal gene transfer (HGT) into the metabolic network. Thus, we end the review by discussing how the combination of promiscuity and HGT contribute to evolution of metabolism in prokaryotes. Future studies investigating the contribution of promiscuity to enzyme and metabolic evolution will need to integrate deeper probes into the influence of evolution on protein biophysics, enzymology, and metabolism with more complex and realistic evolutionary models. ENZYMES: lactate dehydrogenase (EC 1.1.1.27), malate dehydrogenase (EC 1.1.1.37), OSBS (EC 4.2.1.113), HisA (EC 5.3.1.16), TrpF, PriA (EC 5.3.1.24), R-mandelonitrile lyase (EC 4.1.2.10), Maleylacetate reductase (EC 1.3.1.32).

Comprehensive Flux Modeling of Chlamydia trachomatis Proteome and qRT-PCR Data Indicate Biphasic Metabolic Differences Between Elementary Bodies and Reticulate Bodies During Infection

Metabolic adaptation to the host cell is important for obligate intracellular pathogens such as Chlamydia trachomatis (Ct). Here we infer the flux differences for Ct from proteome and qRT-PCR data by comprehensive pathway modeling. We compare the comparatively inert infectious elementary body (EB) and the active replicative reticulate body (RB) systematically using a genome-scale metabolic model with 321 metabolites and 277 reactions. This did yield 84 extreme pathways based on a published proteomics dataset at three different time points of infection. Validation of predictions was done by quantitative RT-PCR of enzyme mRNA expression at three time points. Ct’s major active pathways are glycolysis, gluconeogenesis, glycerol-phospholipid (GPL) biosynthesis (support from host acetyl-CoA) and pentose phosphate pathway (PPP), while its incomplete TCA and fatty acid biosynthesis are less active. The modeled metabolic pathways are much more active in RB than in EB. Our in silico model suggests that EB and RB utilize folate to generate NAD(P)H using independent pathways. The only low metabolic flux inferred for EB involves mainly carbohydrate metabolism. RB utilizes energy -rich compounds to generate ATP in nucleic acid metabolism. Validation data for the modeling include proteomics experiments (model basis) as well as qRT-PCR confirmation of selected metabolic enzyme mRNA expression differences. The metabolic modeling is made fully available here. Its detailed insights and models on Ct metabolic adaptations during infection are a useful modeling basis for future studies.

An N-nitrosating metalloenzyme constructs the pharmacophore of streptozotocin

Small molecules containing the N-nitroso group, such as the bacterial natural product streptozotocin, are prominent carcinogens1,2 and important cancer chemotherapeutics3,4. Despite the considerable importance of this functional group to human health, enzymes dedicated to the assembly of the N-nitroso unit have not been identified. Here we show that SznF, a metalloenzyme from the biosynthesis of streptozotocin, catalyses an oxidative rearrangement of the guanidine group of Nω-methyl-L-arginine to generate an N-nitrosourea product. Structural characterization and mutagenesis of SznF reveal two separate active sites that promote distinct steps in this transformation using different iron-containing metallocofactors. This biosynthetic reaction, which has little precedent in enzymology or organic synthesis, expands the catalytic capabilities of non-haem-iron-dependent enzymes to include N-N bond formation. We find that biosynthetic gene clusters that encode SznF homologues are widely distributed among bacteria-including environmental organisms, plant symbionts and human pathogens-which suggests an unexpectedly diverse and uncharacterized microbial reservoir of bioactive N-nitroso metabolites.

Folate biosynthesis pathway: mechanisms and insights into drug design for infectious diseases

Folate pathway is a key target for the development of new drugs against infectious diseases since the discovery of sulfa drugs and trimethoprim. The knowledge about this pathway has increased in the last years and the catalytic mechanism and structures of all enzymes of the pathway are fairly understood. In addition, differences among enzymes from prokaryotes and eukaryotes could be used for the design of specific inhibitors. In this review, we show a panorama of progress that has been achieved within the folate pathway obtained in the last years. We explored the structure and mechanism of enzymes, several genetic features, strategies, and approaches used in the design of new inhibitors that have been used as targets in pathogen chemotherapy.

Chlamydia trachomatis dapF Encodes a Bifunctional Enzyme Capable of Both d-Glutamate Racemase and Diaminopimelate Epimerase Activities

Peptidoglycan is a sugar/amino acid polymer unique to bacteria and essential for division and cell shape maintenance. The d-amino acids that make up its cross-linked stem peptides are not abundant in nature and must be synthesized by bacteria de novo d-Glutamate is present at the second position of the pentapeptide stem and is strictly conserved in all bacterial species. In Gram-negative bacteria, d-glutamate is generated via the racemization of l-glutamate by glutamate racemase (MurI). Chlamydia trachomatis is the leading cause of infectious blindness and sexually transmitted bacterial infections worldwide. While its genome encodes a majority of the enzymes involved in peptidoglycan synthesis, no murI homologue has ever been annotated. Recent studies have revealed the presence of peptidoglycan in C. trachomatis and confirmed that its pentapeptide includes d-glutamate. In this study, we show that C. trachomatis synthesizes d-glutamate by utilizing a novel, bifunctional homologue of diaminopimelate epimerase (DapF). DapF catalyzes the final step in the synthesis of meso-diaminopimelate, another amino acid unique to peptidoglycan. Genetic complementation of an Escherichia coli murI mutant demonstrated that Chlamydia DapF can generate d-glutamate. Biochemical analysis showed robust activity, but unlike canonical glutamate racemases, activity was dependent on the cofactor pyridoxal phosphate. Genetic complementation, enzymatic characterization, and bioinformatic analyses indicate that chlamydial DapF shares characteristics with other promiscuous/primordial enzymes, presenting a potential mechanism for d-glutamate synthesis not only in Chlamydia but also numerous other genera within the Planctomycetes-Verrucomicrobiae-Chlamydiae superphylum that lack recognized glutamate racemases.IMPORTANCE Here we describe one of the last remaining "missing" steps in peptidoglycan synthesis in pathogenic Chlamydia species, the synthesis of d-glutamate. We have determined that the diaminopimelate epimerase (DapF) encoded by Chlamydia trachomatis is capable of carrying out both the epimerization of DAP and the pyridoxal phosphate-dependent racemization of glutamate. Enzyme promiscuity is thought to be the hallmark of early microbial life on this planet, and there is currently an active debate as to whether "moonlighting enzymes" represent primordial evolutionary relics or are a product of more recent reductionist evolutionary pressures. Given the large number of Chlamydia species (as well as members of the Planctomycetes-Verrucomicrobiae-Chlamydiae superphylum) that possess DapF but lack homologues of MurI, it is likely that DapF is a primordial isomerase that functions as both racemase and epimerase in these organisms, suggesting that specialized d-glutamate racemase enzymes never evolved in these microbes.

Genomewide Transcriptional Responses of Iron-Starved Chlamydia trachomatis Reveal Prioritization of Metabolic Precursor Synthesis over Protein Translation

By utilizing an experimental approach that monitors the immediate global response of Chlamydia trachomatis to iron starvation, clues to long-standing issues in Chlamydia biology are revealed, including how Chlamydia adapts to this stress. We determined that this pathogen initiates a transcriptional program that prioritizes replenishment of nutrient stores over replication, possibly in preparation for rapid growth once optimal iron levels are restored. Transcription of genes for biosynthesis of metabolic precursors was generally upregulated, while those involved in multiple steps of translation were downregulated. We also observed an increase in transcription of genes involved in DNA repair and neutralizing oxidative stress, indicating that Chlamydia employs an “all-or-nothing” strategy. Its small genome limits its ability to tailor a specific response to a particular stress. Therefore, the “all-or-nothing” strategy may be the most efficient way of surviving within the host, where the pathogen likely encounters multiple simultaneous immunological and nutritional insults.

Iron is essential for growth and development of Chlamydia. Its long-term starvation in cultured mammalian cells leads to production of aberrant noninfectious chlamydial forms, also known as persistence. Immediate transcriptional responses to iron limitation have not been characterized, leaving a knowledge gap of how Chlamydia regulates its response to changes in iron availability. We used the fast-chelating agent 2,2′-bipyridyl (BPDL) to homogeneously starve Chlamydia trachomatis serovar L2 of iron, starting at 6 or 12 h postinfection. Immediate transcriptional responses were monitored after only 3 or 6 h of BPDL treatment, well before formation of aberrant Chlamydia. The first genomewide transcriptional response of C. trachomatis to iron starvation was subsequently determined utilizing RNA sequencing. Only 7% and 8% of the genome were differentially expressed in response to iron starvation at the early and middle stages of development, respectively. Biological pathway analysis revealed an overarching theme. Synthesis of macromolecular precursors (deoxynucleotides, amino acids, charged tRNAs, and acetyl coenzyme A [acetyl-CoA]) was upregulated, while energy-expensive processes (ABC transport and translation) were downregulated. A large fraction of differentially downregulated genes are involved in translation, including those encoding ribosome assembly and initiation and termination factors, which could be analogous to the translation downregulation triggered by stress in other prokaryotes during stringent responses. Additionally, transcriptional upregulation of DNA repair, oxidative stress, and tryptophan salvage genes reveals a possible coordination of responses to multiple antimicrobial and immunological insults. These responses of replicative-phase Chlamydia to iron starvation indicate a prioritization of survival over replication, enabling the pathogen to “stock the pantry” with ingredients needed for rapid growth once optimal iron levels are restored.

IMPORTANCE By utilizing an experimental approach that monitors the immediate global response of Chlamydia trachomatis to iron starvation, clues to long-standing issues in Chlamydia biology are revealed, including how Chlamydia adapts to this stress. We determined that this pathogen initiates a transcriptional program that prioritizes replenishment of nutrient stores over replication, possibly in preparation for rapid growth once optimal iron levels are restored. Transcription of genes for biosynthesis of metabolic precursors was generally upregulated, while those involved in multiple steps of translation were downregulated. We also observed an increase in transcription of genes involved in DNA repair and neutralizing oxidative stress, indicating that Chlamydia employs an “all-or-nothing” strategy. Its small genome limits its ability to tailor a specific response to a particular stress. Therefore, the “all-or-nothing” strategy may be the most efficient way of surviving within the host, where the pathogen likely encounters multiple simultaneous immunological and nutritional insults.

Mechanisms of Resistance to Folate Pathway Inhibitors in Burkholderia pseudomallei: Deviation from the Norm

The trimethoprim and sulfamethoxazole combination, co-trimoxazole, plays a vital role in the treatment of Burkholderia pseudomallei infections. Previous studies demonstrated that the B. pseudomallei BpeEF-OprC efflux pump confers widespread trimethoprim resistance in clinical and environmental isolates, but this is not accompanied by significant resistance to co-trimoxazole. Using the excluded select-agent strain B. pseudomallei Bp82, we now show that in vitro acquired trimethoprim versus co-trimoxazole resistance is mainly mediated by constitutive BpeEF-OprC expression due to bpeT mutations or by BpeEF-OprC overexpression due to bpeS mutations. Mutations in bpeT affect the carboxy-terminal effector-binding domain of the BpeT LysR-type activator protein. Trimethoprim resistance can also be mediated by dihydrofolate reductase (FolA) target mutations, but this occurs rarely unless BpeEF-OprC is absent. BpeS is a transcriptional regulator that is 62% identical to BpeT. Mutations affecting the BpeS DNA-binding or carboxy-terminal effector-binding domains result in constitutive BpeEF-OprC overexpression, leading to trimethoprim and sulfamethoxazole efflux and thus to co-trimoxazole resistance. The majority of laboratory-selected co-trimoxazole-resistant mutants often also contain mutations in folM, encoding a pterin reductase. Genetic analyses of these mutants established that both bpeS mutations and folM mutations contribute to co-trimoxazole resistance, although the exact role of folM remains to be determined. Mutations affecting bpeT, bpeS, and folM are common in co-trimoxazole-resistant clinical isolates, indicating that mutations affecting these genes are clinically significant. Co-trimoxazole resistance in B. pseudomallei is a complex phenomenon, which may explain why resistance to this drug is rare in this bacterium.

Burkholderia pseudomallei causes melioidosis, a tropical disease that is difficult to treat. The bacterium’s resistance to antibiotics limits therapeutic options. The paucity of orally available drugs further complicates therapy. The oral drug of choice is co-trimoxazole, a combination of trimethoprim and sulfamethoxazole. These antibiotics target two distinct enzymes, FolA (dihydrofolate reductase) and FolP (dihydropteroate synthase), in the bacterial tetrahydrofolate biosynthetic pathway. Although co-trimoxazole resistance is minimized due to two-target inhibition, bacterial resistance due to folA and folP mutations does occur. Co-trimoxazole resistance in B. pseudomallei is rare and has not yet been studied. Co-trimoxazole resistance in this bacterium employs a novel strategy involving differential regulation of BpeEF-OprC efflux pump expression that determines the drug resistance profile. Contributing are mutations affecting folA, but not folP, and folM, a folate pathway-associated gene whose function is not yet well understood and which has not been previously implicated in folate inhibitor resistance in clinical isolates.

Multifunctional enzymes from reduced genomes - model proteins for simple primordial metabolism?

Billions of years of evolution have yielded today's complex metabolic networks driven by efficient and highly specialized enzymes. In contrast, the metabolism of the earliest cellular life forms was likely much simpler with only a few enzymes of comparatively low activity. It has been speculated that these early enzymes had low specificities and in turn were able to perform multiple functions. In this issue of Molecular Microbiology, Ferla et al. describe examples of enzymes that catalyze chemically distinct reactions while using the same active site. Most importantly, the authors demonstrated that the comparatively weak activities of these multifunctional enzymes are each physiologically relevant. These findings contrast with simply promiscuous enzyme activities, which have been described numerous times but are not physiologically relevant. Ferla et al. elegantly combined initial bioinformatics searches for enzyme candidates with sound kinetic measurements, evolutionary considerations and even structural discussions. The phenomenon of multifunctionality appears to be a mechanism for bacteria with reduced genomes to compensate for their lack of certain enzymes. In the broader context of evolution, these organisms could be considered living model systems to study features of long-extinct early cellular life.

A glimpse into the genetic basis of symbiosis between Hydrogenophaga and their helper strains in the biodegradation of 4-aminobenzenesulfonate

We report the whole genome sequences of Hydrogenophaga intermedia S1 and Agrobacterium radiobacter S2, the first reported bacterial co-culture capable of degrading 4-aminobenzenesulfonate (4-ABS), a recalcitrant industrial waste product. To gain insights into the genetic basis for the syntrophic interaction between this symbiotic pair and also another recently reported Hydrogenophaga associated co-culture, Hydrogenophaga sp. PBC and Ralstonia sp. PBA, we performed detailed genetic analysis of these four strains focusing on the metabolic pathways associated with biotin, para-aminobenzoic acid (pABA), and protocatechuate metabolism. Both assembled Hydrogenophaga draft genomes are missing a majority of the genetic components associated in the biosynthetic pathway of pABA and biotin. Interestingly, a fused pABA synthase was found in R. sp PBA but not in A. radiobacter S2. Furthermore, using whole genome data, the taxonomic classification of R. sp. PBA and A. radiobacter S2 (both previously inferred from 16S rRNA gene) was re-investigated, providing new evidence to propose for their re-classification at the genus and species level, respectively

Evolution of substrate specificity in a retained enzyme driven by gene loss

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

Differential regulation of riboflavin supply genes in Vibrio cholerae

Background

Riboflavin is the precursor of important redox cofactors such as flavin mononucleotide (FMN) and flavin adenine dinucleotide, required for several biological processes. Vibrio cholerae, a pathogenic bacterium responsible for the cholera disease, possesses the ability to biosynthesize de novo as well as to uptake riboflavin through the riboflavin biosynthetic pathway (RBP) and the RibN importer, respectively. The intra-organism relationship between riboflavin biosynthesis and uptake functions has not been studied.

Results

This work determined the transcriptional organization of RBP genes and ribN in V. cholerae through reverse transcription polymerase chain reaction and analyzed their expression when growing with or without extracellular riboflavin using real time PCR. The RBP is organized in three transcriptional units, the major one containing ribD, ribE, ribA and ribH together with genes involved in functions not directly related to riboflavin biosynthesis such as nrdR and nusB. In addition, two independent monocistronic units contain ribA2 and ribB, the later conserving a putative FMN riboswitch. The ribN gene is encoded in operon with a gene coding for a predicted outer membrane protein and a gene encoding a protein with a glutaredoxin domain. Regulation analysis showed that among these transcriptional units, only ribB is negatively regulated by riboflavin and that its repression depends on the RibN riboflavin importer. Moreover, external riboflavin highly induced ribB transcription in a ΔribN strain. Also, a genomic database search found a negative correlation between the presence of nrdR and nusB and the FMN riboswitch in bacterial RBP operons.

Conclusions

Growing in the presence of riboflavin downregulates only a single element among the transcriptional units of riboflavin supply pathways. Thus, endogenous riboflavin biosynthesis seems to be negatively regulated by extracellular riboflavin through its specific effect on transcription of ribB in V. cholerae.

Electronic supplementary material

The online version of this article (doi:10.1186/s13099-017-0159-z) contains supplementary material, which is available to authorized users.

Underlying mechanisms for syntrophic metabolism of essential enzyme cofactors in microbial communities

Many microorganisms are unable to synthesize essential B vitamin-related enzyme cofactors de novo. The underlying mechanisms by which such microbes survive in multi-species communities are largely unknown. We previously reported the near-complete genome sequence of two ~18-member unicyanobacterial microbial consortia that maintain stable membership on defined medium lacking vitamins. Here we have used genome analysis and growth studies on isolates derived from the consortia to reconstruct pathways for biogenesis of eight essential cofactors and predict cofactor usage and precursor exchange in these communities. Our analyses revealed that all but the two Halomonas and cyanobacterial community members were auxotrophic for at least one cofactor. We also observed a mosaic distribution of salvage routes for a variety of cofactor precursors, including those produced by photolysis. Potentially bidirectional transporters were observed to be preferentially in prototrophs, suggesting a mechanism for controlled precursor release. Furthermore, we found that Halomonas sp. do not require cobalamin nor control its synthesis, supporting the hypothesis that they overproduce and export vitamins. Collectively, these observations suggest that the consortia rely on syntrophic metabolism of cofactors as a survival strategy for optimization of metabolic exchange within a shared pool of micronutrients.

Overlapping riboflavin supply pathways in bacteria

Riboflavin derivatives are essential cofactors for a myriad of flavoproteins. In bacteria, flavins importance extends beyond their role as intracellular protein cofactors, as secreted flavins are a key metabolite in a variety of physiological processes. Bacteria obtain riboflavin through the endogenous riboflavin biosynthetic pathway (RBP) or by the use of importer proteins. Bacteria frequently encode multiple paralogs of the RBP enzymes and as for other micronutrient supply pathways, biosynthesis and uptake functions largely coexist. It is proposed that bacteria shut down biosynthesis and would rather uptake riboflavin when the vitamin is environmentally available. Recently, the overlap of riboflavin provisioning elements has gained attention and the functions of duplicated paralogs of RBP enzymes started to be addressed. Results point towards the existence of a modular structure in the bacterial riboflavin supply pathways. Such structure uses subsets of RBP genes to supply riboflavin for specific functions. Given the importance of riboflavin in intra and extracellular bacterial physiology, this complex array of riboflavin provision pathways may have developed to contend with the various riboflavin requirements. In riboflavin-prototrophic bacteria, riboflavin transporters could represent a module for riboflavin provision for particular, yet unidentified processes, rather than substituting for the RBP as usually assumed.

Quantitative Protein Profiling of Chlamydia trachomatis Growth Forms Reveals Defense Strategies Against Tryptophan Starvation

Chlamydia trachomatis is one of the most common sexually transmitted bacterial pathogens in humans. The infection is often asymptomatic and can lead to chronic manifestations. The infectious elementary body and the replicating reticulate body are the two growth forms in the normal developmental cycle. Under the influence of interferon-γ, the normal cycle is disrupted because of tryptophan degradation, leading to a third persistent form, the aberrant reticulate body.

For the genital strain C. trachomatis D/UW-3/CX we established a quantitative, label-free proteomic approach, and identified in total 655 out of 903 (73%) predicted proteins, allowing the first quantitative comparison of all three growth forms. Inclusion membrane proteins and proteins involved in translation were more abundant in the reticulate body (RB)¹ and aberrant reticulate body (ARB) forms, whereas proteins of the type III Secretion System and the cell envelope were more abundant in the elementary body (EB) form, reflecting the need for these proteins to establish infection and for host interactions.

In the interferon-γ induced ARB proteome, the tryptophan synthase subunits were identified as biomarkers with a strong increase from less than 0.05% to 9% of the total protein content, reflecting an inherent defense strategy for the pathogen to escape interferon-γ mediated immune pressure. Furthermore, the total tryptophan content in the ARB form was 1.9-fold lower compared with the EB form, and we demonstrate that modulation of the protein repertoire toward lower abundance of proteins with high tryptophan content, is a mechanism which contributes to establish and maintain chlamydial persistence. Thus, quantitative proteomics provides insights in the Chlamydia defense mechanisms to escape interferon-γ mediated immune pressure.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

Mechanistic and Evolutionary Insights from the Reciprocal Promiscuity of Two Pyridoxal Phosphate-dependent Enzymes

Enzymes that utilize the cofactor pyridoxal 5′-phosphate play essential roles in amino acid metabolism in all organisms. The cofactor is used by proteins that adopt at least five different folds, which raises questions about the evolutionary processes that might explain the observed distribution of functions among folds. In this study, we show that a representative of fold type III, the Escherichia coli alanine racemase (ALR), is a promiscuous cystathionine β-lyase (CBL). Furthermore, E. coli CBL (fold type I) is a promiscuous alanine racemase. A single round of error-prone PCR and selection yielded variant ALR(Y274F), which catalyzes cystathionine β-elimination with a near-native Michaelis constant (K_m = 3.3 mm) but a poor turnover number (k_cat ≈10 h⁻¹). In contrast, directed evolution also yielded CBL(P113S), which catalyzes l-alanine racemization with a poor K_m (58 mm) but a high k_cat (22 s⁻¹). The structures of both variants were solved in the presence and absence of the l-alanine analogue, (R)-1-aminoethylphosphonic acid. As expected, the ALR active site was enlarged by the Y274F substitution, allowing better access for cystathionine. More surprisingly, the favorable kinetic parameters of CBL(P113S) appear to result from optimizing the pK_a of Tyr-111, which acts as the catalytic acid during l-alanine racemization. Our data emphasize the short mutational routes between the functions of pyridoxal 5′-phosphate-dependent enzymes, regardless of whether or not they share the same fold. Thus, they confound the prevailing model of enzyme evolution, which predicts that overlapping patterns of promiscuity result from sharing a common multifunctional ancestor.

The Rickettsia Endosymbiont of Ixodes pacificus Contains All the Genes of De Novo Folate Biosynthesis

Ticks and other arthropods often are hosts to nutrient providing bacterial endosymbionts, which contribute to their host’s fitness by supplying nutrients such as vitamins and amino acids. It has been detected, in our lab, that Ixodes pacificus is host to Rickettsia species phylotype G021. This endosymbiont is predominantly present, and 100% maternally transmitted in I. pacificus. To study roles of phylotype G021 in I. pacificus, bioinformatic and molecular approaches were carried out. MUMmer genome alignments of whole genome sequence of I. scapularis, a close relative to I. pacificus, against completely sequenced genomes of R. bellii OSU85-389, R. conorii, and R. felis, identified 8,190 unique sequences that are homologous to Rickettsia sequences in the NCBI Trace Archive. MetaCyc metabolic reconstructions revealed that all folate gene orthologues (folA, folC, folE, folKP, ptpS) required for de novo folate biosynthesis are present in the genome of Rickettsia buchneri in I. scapularis. To examine the metabolic capability of phylotype G021 in I. pacificus, genes of the folate biosynthesis pathway of the bacterium were PCR amplified using degenerate primers. BLAST searches identified that nucleotide sequences of the folA, folC, folE, folKP, and ptpS genes possess 98.6%, 98.8%, 98.9%, 98.5% and 99.0% identity respectively to the corresponding genes of Rickettsia buchneri. Phylogenetic tree constructions show that the folate genes of phylotype G021 and homologous genes from various Rickettsia species are monophyletic. This study has shown that all folate genes exist in the genome of Rickettsia species phylotype G021 and that this bacterium has the genetic capability for de novo folate synthesis.

Insights into the evolution of enzyme substrate promiscuity after the discovery of (βα)₈ isomerase evolutionary intermediates from a diverse metagenome

BACKGROUND

Current 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.

RESULTS

We 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.

CONCLUSION

As 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.

A substrate ambiguous enzyme facilitates genome reduction in an intracellular symbiont

BACKGROUND

Genome evolution in intracellular microbial symbionts is characterized by gene loss, generating some of the smallest and most gene-poor genomes known. As a result of gene loss these genomes commonly contain metabolic pathways that are fragmented relative to their free-living relatives. The evolutionary retention of fragmented metabolic pathways in the gene-poor genomes of endosymbionts suggests that they are functional. However, it is not always clear how they maintain functionality. To date, the fragmented metabolic pathways of endosymbionts have been shown to maintain functionality through complementation by host genes, complementation by genes of another endosymbiont and complementation by genes in host genomes that have been horizontally acquired from a microbial source that is not the endosymbiont. Here, we demonstrate a fourth mechanism.

RESULTS

We investigate the evolutionary retention of a fragmented pathway for the essential nutrient pantothenate (vitamin B5) in the pea aphid, Acyrthosiphon pisum endosymbiosis with Buchnera aphidicola. Using quantitative analysis of gene expression we present evidence for complementation of the Buchnera pantothenate biosynthesis pathway by host genes. Further, using complementation assays in an Escherichia coli mutant we demonstrate functional replacement of a pantothenate biosynthesis enzyme, 2-dehydropantoate 2-reductase (E.C. 1.1.1.169), by an endosymbiont gene, ilvC, encoding a substrate ambiguous enzyme.

CONCLUSIONS

Earlier studies have speculated that missing enzyme steps in fragmented endosymbiont metabolic pathways are completed by adaptable endosymbiont enzymes from other pathways. Here, we experimentally demonstrate completion of a fragmented endosymbiont vitamin biosynthesis pathway by recruitment of a substrate ambiguous enzyme from another pathway. In addition, this work extends host/symbiont metabolic collaboration in the aphid/Buchnera symbiosis from amino acid metabolism to include vitamin biosynthesis.

Variations in metabolic pathways create challenges for automated metabolic reconstructions: Examples from the tetrahydrofolate synthesis pathway

The availability of thousands of sequenced genomes has revealed the diversity of biochemical solutions to similar chemical problems. Even for molecules at the heart of metabolism, such as cofactors, the pathway enzymes first discovered in model organisms like Escherichia coli or Saccharomyces cerevisiae are often not universally conserved. Tetrahydrofolate (THF) (or its close relative tetrahydromethanopterin) is a universal and essential C1-carrier that most microbes and plants synthesize de novo. The THF biosynthesis pathway and enzymes are, however, not universal and alternate solutions are found for most steps, making this pathway a challenge to annotate automatically in many genomes. Comparing THF pathway reconstructions and functional annotations of a chosen set of folate synthesis genes in specific prokaryotes revealed the strengths and weaknesses of different microbial annotation platforms. This analysis revealed that most current platforms fail in metabolic reconstruction of variant pathways. However, all the pieces are in place to quickly correct these deficiencies if the different databases were built on each other's strengths.