Structure and activity analyses of Escherichia coli K-12 NagD provide insight into the evolution of biochemical function in the haloalkanoic acid dehalogenase superfamily

Systematic Metabolic Engineering for Enhanced Cytidine 5'-Monophosphate Production in Escherichia coli

Cytidine 5'-monophosphate (5'-CMP) is a key intermediate in various nucleotide derivatives and is widely used in the food and pharmaceutical industries. In this study, the titer of 5'-CMP in engineered Escherichia coli CM006, which blocks the degradation pathway of 5'-CMP, increased 134.7-fold compared to E. coli MG1655. Integrated expression of the 5'-CMP diphosphate hydrolase gene nudG significantly enhanced the 5'-CMP titer, reduced orotic acid accumulation, and alleviated feedback inhibition in the 5'-CMP synthesis pathway. By blocking cytidine degradation and increasing carbon flux from cytidine to 5'-CMP, the 5'-CMP titer in E. coli CM011 was increased to 338.2 mg/L. Furthermore, by enhancing the supply of the precursor 5-phospho-α-D-ribose 1-diphosphate (PRPP), the 5'-CMP titer in E. coli CM012 reached 417.9 mg/L, the highest value reported to date. The results and methods presented in this study are of great significance for the sustainable industrial production of 5'-CMP.

The activity of the ribonucleotide monophosphatase UmpH is controlled by interaction with the GlnK signaling protein in Escherichia coli

The PII signaling proteins are ubiquitous in prokaryotes serving as crucial metabolic hubs in different metabolic pathways because of their ability to sense and integrate signals of the cellular nitrogen, carbon, and energy levels. In this study, we used ligand fishing assays to identify the ribonucleotide monophosphatase UmpH enzyme as a novel target of the PII signaling protein GlnK in Escherichia coli. In vitro analyses showed that UmpH interacts specifically with the PII protein GlnK but not with its paralog protein GlnB. The UmpH-GlnK complex is modulated by the GlnK uridylylation status and by the levels of the GlnK allosteric effectors ATP, ADP, and 2-oxoglutarate. Upon engaging interaction with GlnK, UmpH becomes less active toward its substrate uridine 5'-monophosphate. We suggest a model where GlnK will physically interact to reduce the UmpH activity during the transition from N-starvation to N-sufficient conditions. Such a mechanism may help the cells to reprogram the fate of uridine 5'-monophosphate from catabolism to anabolism avoiding futile cycling of key nutrients.

MbovP0725, a secreted serine/threonine phosphatase, inhibits the host inflammatory response and affects metabolism in Mycoplasma bovis

Mycoplasma species are able to produce and release secreted proteins, such as toxins, adhesins, and virulence-related enzymes, involved in bacteria adhesion, invasion, and immune evasion between the pathogen and host. Here, we investigated a novel secreted protein, MbovP0725, from Mycoplasma bovis encoding a putative haloacid dehalogenase (HAD) hydrolase function of a key serine/threonine phosphatase depending on Mg2+ for the dephosphorylation of its substrate pNPP, and it was most active at pH 8 to 9 and temperatures around 40°C. A transposon insertion mutant strain of M. bovis HB0801 that lacked the protein MbovP0725 induced a stronger inflammatory response but with a partial reduction of adhesion ability. Using transcriptome sequencing and quantitative reverse transcription polymerase chain reaction analysis, we found that the mutant was upregulated by the mRNA expression of genes from the glycolysis pathway, while downregulated by the genes enriched in ABC transporters and acetate kinase-phosphate acetyltransferase pathway. Untargeted metabolomics showed that the disruption of the Mbov_0725 gene caused the accumulation of 9-hydroxyoctadecadienoic acids and the consumption of cytidine 5'-monophosphate, uridine monophosphate, and adenosine monophosphate. Both the exogenous and endogenous MbvoP0725 protein created by purification and transfection inhibited lipopolysaccharide (LPS)-induced IL-1β, IL-6, and TNF-α mRNA production and could also attenuate the activation of MAPK-associated pathways after LPS treatment. A pull-down assay identified MAPK p38 and ERK as potential substrates for MbovP0725. These findings define metabolism- and virulence-related roles for a HAD family phosphatase and reveal its ability to inhibit the host pro-inflammatory response.

IMPORTANCE

Mycoplasma bovis (M. bovis) infection is characterized by chronic pneumonia, otitis, arthritis, and mastitis, among others, and tends to involve the suppression of the immune response via multiple strategies to avoid host cell immune clearance. This study found that MbovP0725, a haloacid dehalogenase (HAD) family phosphatase secreted by M. bovis, had the ability to inhibit the host pro-inflammatory response induced by lipopolysaccharide. Transcriptomic and metabolomic analyses were used to identify MbovP0725 as an important phosphatase involved in glycolysis and nucleotide metabolism. The M. bovis transposon mutant strain T8.66 lacking MbovP0725 induced a higher inflammatory response and exhibited weaker adhesion to host cells. Additionally, T8.66 attenuated the phosphorylation of MAPK P38 and ERK and interacted with the two targets. These results suggested that MbovP0725 had the virulence- and metabolism-related role of a HAD family phosphatase, performing an anti-inflammatory response during M. bovis infection.

Characterization of haloacid dehalogenase superfamily acid phosphatase from Staphylococcus lugdunensis

The haloacid dehalogenase superfamily implicated in bacterial pathogenesis comprises different enzymes having roles in many metabolic pathways. Staphylococcus lugdunensis, a Gram-positive bacterium, is an opportunistic human pathogen causing infections in the central nervous system, urinary tract, bones, peritoneum, systemic conditions and cutaneous infection. The haloacid dehalogenase superfamily proteins play a significant role in the pathogenicity of certain bacteria, facilitating invasion, survival, and proliferation within host cells. The genome of S. lugdunensis encodes more than ten proteins belonging to this superfamily. However, none of them have been characterized. The present work reports the characterization of one of the haloacid dehalogenase superfamily proteins (SLHAD1) from Staphylococcus lugdunensis. The functional analysis revealed that SLHAD1 is a metal-dependent acid phosphatase, which catalyzes the dephosphorylation of phosphorylated metabolites of cellular pathways, including glycolysis, gluconeogenesis, nucleotides, and thiamine metabolism. Based on the substrate specificity and genomic analysis, the physiological function of SLHAD1 in thiamine metabolism has been tentatively assigned. The crystal structure of SLHAD1, lacking 49 residues at the C-terminal, was determined at 1.7 Å resolution with a homodimer in the asymmetric unit. It was observed that SLHAD1 exhibited time-dependent cleavage at a specific point, occurring through a self-initiated process. A combination of bioinformatics, biochemical, biophysical, and structural studies explored unique features of SLHAD1. Overall, the study revealed a detailed characterization of a critical enzyme of the human pathogen Staphylococcus lugdunensis, associated with several life-threatening infections.

Systematic Engineering of Escherichia coli for Efficient Production of Cytidine 5'-Monophosphate

Cytidine 5'-monophosphate (CMP) was widely applied in the food and pharmaceutical industries. Currently, CMP is mainly produced by enzyme catalysis. However, the starting materials for enzyme catalysis were relatively expensive. Therefore, seeking a low-cost production process for CMP was attractive. In this study, Escherichia coli (E. coli) was systematically modified to produce CMP. First, a the cytidine-producing strain was constructed by deleting cdd, rihA, rihB, and rihC. Second, the genes involved in the pyrimidine precursor competing pathway and negative regulation were deleted to increase cyti dine biosynthesis. Third, the deletion of the genes that caused the loss of CMP phosphatase activity led to the accumulation of CMP, and the overexpression of the rate-limiting step genes and feedback inhibition resistance genes greatly increased the yield of CMP. The yield of CMP was further increased to 1013.6 mg/L by blocking CMP phosphorylation. Ultimately, the yield of CMP reached 15.3 g/L in a 50 L bioreactor. Overall, the engineered E. coli with a high yield of CMP was successfully constructed and showed the potential for industrial production.

Glucose consumption rate-dependent transcriptome profiling of Escherichia coli provides insight on performance as microbial factories

Background

The modification of glucose import capacity is an engineering strategy that has been shown to improve the characteristics of Escherichia coli as a microbial factory. A reduction in glucose import capacity can have a positive effect on production strain performance, however, this is not always the case. In this study, E. coli W3110 and a group of four isogenic derivative strains, harboring single or multiple deletions of genes encoding phosphoenolpyruvate:sugar phosphotransferase system (PTS)-dependent transporters as well as non-PTS transporters were characterized by determining their transcriptomic response to reduced glucose import capacity.

Results

These strains were grown in bioreactors with M9 mineral salts medium containing 20 g/L of glucose, where they displayed specific growth rates ranging from 0.67 to 0.27 h⁻¹, and specific glucose consumption rates (qs) ranging from 1.78 to 0.37 g/g h. RNA-seq analysis revealed a transcriptional response consistent with carbon source limitation among all the mutant strains, involving functions related to transport and metabolism of alternate carbon sources and characterized by a decrease in genes encoding glycolytic enzymes and an increase in gluconeogenic functions. A total of 107 and 185 genes displayed positive and negative correlations with qs, respectively. Functions displaying positive correlation included energy generation, amino acid biosynthesis, and sugar import.

Conclusion

Changes in gene expression of E. coli strains with impaired glucose import capacity could be correlated with qs values and this allowed an inference of the physiological state of each mutant. In strains with lower qs values, a gene expression pattern is consistent with energy limitation and entry into the stationary phase. This physiological state could explain why these strains display a lower capacity to produce recombinant protein, even when they show very low rates of acetate production. The comparison of the transcriptomes of the engineered strains employed as microbial factories is an effective approach for identifying favorable phenotypes with the potential to improve the synthesis of biotechnological products.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01909-y.

The Catabolic System of Acetovanillone and Acetosyringone in Sphingobium sp. Strain SYK-6 Useful for Upgrading Aromatic Compounds Obtained through Chemical Lignin Depolymerization

Acetovanillone is a major aromatic monomer produced in oxidative/base-catalyzed lignin depolymerization. However, the production of chemical products from acetovanillone has not been explored due to the lack of information on the microbial acetovanillone catabolic system. Here, the acvABCDEF genes were identified as specifically induced genes during the growth of Sphingobium sp. strain SYK-6 cells with acetovanillone and these genes were essential for SYK-6 growth on acetovanillone and acetosyringone (a syringyl-type acetophenone derivative). AcvAB and AcvF produced in Escherichia coli phosphorylated acetovanillone/acetosyringone and dephosphorylated the phosphorylated acetovanillone/acetosyringone, respectively. AcvCDE produced in Sphingobium japonicum UT26S carboxylated the reaction products generated from acetovanillone/acetosyringone by AcvAB and AcvF into vanilloyl acetic acid/3-(4-hydroxy-3,5-dimethoxyphenyl)-3-oxopropanoic acid. To demonstrate the feasibility of producing cis,cis-muconic acid from acetovanillone, a metabolic modification on a mutant of Pseudomonas sp. strain NGC7 that accumulates cis,cis-muconic acid from catechol was performed. The resulting strain expressing vceA and vceB required for converting vanilloyl acetic acid to vanillic acid and aroY encoding protocatechuic acid decarboxylase in addition to acvABCDEF successfully converted 1.2 mM acetovanillone to approximately equimolar cis,cis-muconic acid. Our results are expected to help improve the yield and purity of value-added chemical production from lignin through biological funneling. IMPORTANCE In the alkaline oxidation of lignin, aromatic aldehydes (vanillin, syringaldehyde, and p-hydroxybenzaldehyde), aromatic acids (vanillic acid, syringic acid, and p-hydroxybenzoic acid), and acetophenone-related compounds (acetovanillone, acetosyringone, and 4'-hydroxyacetophenone) are produced as major aromatic monomers. Also, base-catalyzed depolymerization of guaiacyl lignin resulted in vanillin, vanillic acid, guaiacol, and acetovanillone as primary aromatic monomers. To date, microbial catabolic systems of vanillin, vanillic acid, and guaiacol have been well characterized, and the production of value-added chemicals from them has also been explored. However, due to the lack of information on the microbial acetovanillone and acetosyringone catabolic system, chemical production from acetovanillone and acetosyringone has not been achieved. This study elucidated the acetovanillone/acetosyringone catabolic system and demonstrates the potential of using these genes for the production of value-added chemicals from these compounds.

Srg1, a putative protein phosphatase from the HAD-family, is involved in stress adaptation in Candida albicans

BACKGROUND

The cell stress response plays an important role in the survival of organisms. Studies have revealed that the pathogenic yeast Candida albicans that constantly encounters various environmental insults inside the host has emerged as an ideal system to understand the molecular mechanism (s) of stress response. In this study, we characterize a stress-inducible gene SRG1 which is a Halo Acid Dehalogenase (HAD) family member from C. albicans.

METHODS

We used confocal microscopy, site-directed mutagenesis, gene deletion techniques, and tandem-affinity purification and co-immunoprecipitation studies to functionally characterize SRG1.

RESULTS

The sub-cellular localization of Srg1 is predominantly cytoplasmic and includes punctate mitochondrial staining in the presence of salt. Protein purification studies coupled with LC-MS analysis showed that Srg1 is a phosphoprotein. The Srg1 mutant carrying S47A and S49A mutations failed to migrate to mitochondria in the presence of salt but retained its phosphatase activity. Srg1 migrates to the nucleus in ∆hog1 mutant cells indicating an unorthodox role for HAD family proteins in stress-mediated transcriptional response. Srg1 also interacts with Erg13, a component involved in the mitochondrial membrane lipid biosynthesis pathway.

CONCLUSIONS

A multistep relay mechanism that includes a positive modulation by the MAP kinase Hog1 and a negative modulation by the global repressor Tup1 controls SRG1 expression.

GENERAL SIGNIFICANCE

Taken together, our work contributes towards gaining a functional insight into a class of phosphatases that probably have evolved with novel specificities in the pathogenic yeast C. albicans to counteract stressful conditions.

Microbial 5'-nucleotidases: their characteristics, roles in cellular metabolism, and possible practical applications

5′-Nucleotidases (EC 3.1.3.5) are enzymes that catalyze the hydrolytic dephosphorylation of 5′-ribonucleotides and 5′-deoxyribonucleotides to their respective nucleosides and phosphate. Most 5′-nucleotidases have broad substrate specificity and are multifunctional enzymes capable of cleaving phosphorus from not only mononucleotide phosphate molecules but also a variety of other phosphorylated metabolites. 5′-Nucleotidases are widely distributed throughout all kingdoms of life and found in different cellular locations. The well-studied vertebrate 5′-nucleotidases play an important role in cellular metabolism. These enzymes are involved in purine and pyrimidine salvage pathways, nucleic acid repair, cell-to-cell communication, signal transduction, control of the ribo- and deoxyribonucleotide pools, etc. Although the first evidence of microbial 5′-nucleotidases was obtained almost 60 years ago, active studies of genetic control and the functions of microbial 5′-nucleotidases started relatively recently. The present review summarizes the current knowledge about microbial 5′-nucleotidases with a focus on their diversity, cellular localizations, molecular structures, mechanisms of catalysis, physiological roles, and activity regulation and approaches to identify new 5′-nucleotidases. The possible applications of these enzymes in biotechnology are also discussed.

Key points

• Microbial 5′-nucleotidases differ in molecular structure, hydrolytic mechanism, and cellular localization.

• 5′-Nucleotidases play important and multifaceted roles in microbial cells.

• Microbial 5′-nucleotidases have wide range of practical applications.

Enzymatic characterization of a thermostable phosphatase from Thermomicrobium roseum and its application for biosynthesis of fructose from maltodextrin

Phosphatases, which catalyze the dephosphorylation of compounds containing phosphate groups, are important members of the haloacid dehalogenase (HAD)-like superfamily. Herein, a thermostable phosphatase encoded by an open reading frame of Trd_1070 from Thermomicrobium roseum was enzymologically characterized. This phosphatase showed promiscuous activity against more than ten sugar phosphates, with high specific activity toward ribose 5-phosphate, followed by ribulose 5-phosphate and fructose 6-phosphate. The half-life of Trd_1070 at 70 °C and pH 7.0 was about 14.2 h. Given that the catalytic efficiency of Trd_1070 on fructose 6-phosphate was 49-fold higher than that on glucose 6-phosphate, an in vitro synthetic biosystem containing alpha-glucan phosphorylase, phosphoglucomutase, phosphoglucose isomerase, and Trd_1070 was constructed for the production of fructose from maltodextrin by whole-cell catalysis, resulting in 21.6 g/L fructose with a ratio of fructose to glucose of approximately 2:1 from 50 g/L maltodextrin. This in vitro biosystem provides an alternative method to produce fructose with higher fructose content compared with the traditional production method using glucose isomerization. Further discovery and enzymologic characterization of phosphatases may promote further production of alternative monosaccharides through in vitro synthetic biosystems.

A Proteomic View of Salmonella Typhimurium in Response to Phosphate Limitation

Salmonella enterica serovar Typhimurium (S. Typhimurium), an important foodborne pathogen, often encounters phosphate (P_i) shortage both in the environment and inside host cells. To gain a global view on its physiological responses to P_i starvation, we performed proteomic profiling of S. Typhimurium upon the shift from P_i-rich to P_i-low conditions. In addition to the Pho regulon, many metabolic processes were up-regulated, such as glycolysis, pentose phosphate pathway, pyrimidine degradation, glycogen, and trehalose metabolism, allowing us to chart an overview of S. Typhimurium carbon metabolism under P_i starvation. Furthermore, proteomic analysis of a mutant lacking phoB (that encodes a key regulator of P_i shortage response) suggested that only a small subset of the altered proteins upon P_i limitation was PhoB-dependent. Importantly, we present evidence that S. Typhimurium N-acetylglucosamine catabolism was induced under P_i-limiting conditions in a PhoB-dependent manner. Immunoblotting and β-galactosidase assays demonstrated that PhoB was required for the full activation of NagB, a key enzyme of this pathway, in response to low P_i. Thus, our study reveals that N-acetylglucosamine catabolism may represent an additional PhoB-regulated pathway to tackle bacterial P_i shortage.

Allosteric Control of Substrate Specificity of the Escherichia coli ADP-Glucose Pyrophosphorylase

The substrate specificity of enzymes is crucial to control the fate of metabolites to different pathways. However, there is growing evidence that many enzymes can catalyze alternative reactions. This promiscuous behavior has important implications in protein evolution and the acquisition of new functions. The question is how the undesirable outcomes of in vivo promiscuity can be prevented. ADP-glucose pyrophosphorylase from Escherichia coli is an example of an enzyme that needs to select the correct substrate from a broad spectrum of alternatives. This selection will guide the flow of carbohydrate metabolism toward the synthesis of reserve polysaccharides. Here, we show that the allosteric activator fructose-1,6-bisphosphate plays a role in such selection by increasing the catalytic efficiency of the enzyme toward the use of ATP rather than other nucleotides. In the presence of fructose-1,6-bisphosphate, the k_cat/S_0.5 for ATP was near ~600-fold higher that other nucleotides, whereas in the absence of activator was only ~3-fold higher. We propose that the allosteric regulation of certain enzymes is an evolutionary mechanism of adaptation for the selection of specific substrates.

Identification, Heterologous Expression, and Functional Characterization of Bacillus subtilis YutF, a HAD Superfamily 5'-Nucleotidase with Broad Substrate Specificity

5'-nucleotidases (EC 3.1.3.5) catalyze the hydrolytic dephosphorylation of 5'-ribonucleotides and 5'-deoxyribonucleotides as well as complex nucleotides, such as uridine 5'-diphosphoglucose (UDP-glucose), nicotinamide adenine dinucleotide and flavin adenine dinucleotide, to their corresponding nucleosides plus phosphate. These enzymes have been found in diverse species in intracellular and membrane-bound, surface-localized forms. Soluble forms of 5'-nucleotidases belong to the ubiquitous haloacid dehalogenase superfamily (HADSF) and have been shown to be involved in the regulation of nucleotide, nucleoside and nicotinamide adenine dinucleotide (NAD⁺) pools. Despite the important role of 5'-nucleotidases in cellular metabolism, only a few of these enzymes have been characterized in the Gram-positive bacterium Bacillus subtilis, the workhorse industrial microorganism included in the Food and Drug Administration’s GRAS (generally regarded as safe) list. In the present study, we report the identification of a novel 5'-nucleotidase gene from B. subtilis, yutF, which comprises 771 bp encoding a 256-amino-acid protein belonging to the IIA subfamily of the HADSF. The gene product is responsible for the major p-nitrophenyl phosphatase activity in B. subtilis. The yutF gene was overexpressed in Escherichia coli, and its product fused to a polyhistidine tag was purified and biochemically characterized as a soluble 5'-nucleotidase with broad substrate specificity. The recombinant YutF protein was found to hydrolyze various purine and pyrimidine 5'-nucleotides, showing preference for 5'-nucleoside monophosphates and, specifically, 5'-XMP. Recombinant YutF also exhibited phosphohydrolase activity toward nucleotide precursors, ribose-5-phosphate and 5-phosphoribosyl-1-pyrophosphate. Determination of the kinetic parameters of the enzyme revealed a low substrate specificity (K_m values in the mM concentration range) and modest catalytic efficiencies with respect to substrates. An initial study of the regulation of yutF expression showed that the yutF gene is a component of the yutDEF transcription unit and that YutF overproduction positively influences yutDEF expression.

Bacillus subtilis 5'-nucleotidases with various functions and substrate specificities

BACKGROUND

In Escherichia coli, nagD, yrfG, yjjG, yieH, yigL, surE, and yfbR encode 5'-nucleotidases that hydrolyze the phosphate group of 5'-nucleotides. In Bacillus subtilis, genes encoding 5'-nucleotidase have remained to be identified.

RESULTS

We found that B. subtilis ycsE, araL, yutF, ysaA, and yqeG show suggestive similarities to nagD. Here, we expressed them in E. coli to purify the respective His₆-tagged proteins. YcsE exhibited significant 5'-nucleotidase activity with a broader specificity, whereas the other four enzymes had rather weak but suggestive activities with various capacities and substrate specificities. In contrast, B. subtilis yktC shares high similarity with E. coli suhB encoding an inositol monophosphatase. YktC exhibited inositol monophosphatase activity as well as 5'-nucleotidase activity preferential for GMP and IMP. The ycsE, yktC, and yqeG genes are induced by oxidative stress and were dispensable, although yqeG was required to maintain normal growth on solid medium. In the presence of diamide, only mutants lacking yktC exhibited enhanced growth defects, whereas the other mutants without ycsE or yqeG did not.

CONCLUSIONS

Accordingly, in B. subtilis, at least YcsE and YktC acted as major 5'-nucleotidases and the four minor enzymes might function when the intracellular concentrations of substrates are sufficiently high. In addition, YktC is involved in resistance to oxidative stress caused by diamide, while YqeG is necessary for normal colony formation on solid medium.

Regulation of the Utilization of Amino Sugars by and : Same Genes, Different Control

Amino sugars are dual-purpose compounds in bacteria: they are essential components of the outer wall peptidoglycan (PG) and the outer membrane of Gram-negative bacteria and, in addition, when supplied exogenously their catabolism contributes valuable supplies of energy, carbon and nitrogen to the cell. The enzymes for both the synthesis and degradation of glucosamine (GlcN) and N-acetylglucosamine (GlcNAc) are highly conserved but during evolution have become subject to different regulatory regimes. <i>Escherichia coli</i> grows more rapidly using GlcNAc as a carbon source than with GlcN. On the other hand, <i>Bacillus subtilis,</i> but not other <i>Bacilli</i> tested, grows more efficiently on GlcN than GlcNAc. The more rapid growth on this sugar is associated with the presence of a second, GlcN-specific operon, which is unique to this species. A single locus is associated with the genes for catabolism of GlcNAc and GlcN in <i>E. coli,</i> although they enter the cell via different transporters. In <i>E. coli</i> the amino sugar transport and catabolic genes have also been requisitioned as part of the PG recycling process. Although PG recycling likely occurs in <i>B. subtilis,</i> it appears to have different characteristics.

Functional Diversity of Haloacid Dehalogenase Superfamily Phosphatases from Saccharomyces cerevisiae: BIOCHEMICAL, STRUCTURAL, AND EVOLUTIONARY INSIGHTS

The haloacid dehalogenase (HAD)-like enzymes comprise a large superfamily of phosphohydrolases present in all organisms. The Saccharomyces cerevisiae genome encodes at least 19 soluble HADs, including 10 uncharacterized proteins. Here, we biochemically characterized 13 yeast phosphatases from the HAD superfamily, which includes both specific and promiscuous enzymes active against various phosphorylated metabolites and peptides with several HADs implicated in detoxification of phosphorylated compounds and pseudouridine. The crystal structures of four yeast HADs provided insight into their active sites, whereas the structure of the YKR070W dimer in complex with substrate revealed a composite substrate-binding site. Although the S. cerevisiae and Escherichia coli HADs share low sequence similarities, the comparison of their substrate profiles revealed seven phosphatases with common preferred substrates. The cluster of secondary substrates supporting significant activity of both S. cerevisiae and E. coli HADs includes 28 common metabolites that appear to represent the pool of potential activities for the evolution of novel HAD phosphatases. Evolution of novel substrate specificities of HAD phosphatases shows no strict correlation with sequence divergence. Thus, evolution of the HAD superfamily combines the conservation of the overall substrate pool and the substrate profiles of some enzymes with remarkable biochemical and structural flexibility of other superfamily members.

Panoramic view of a superfamily of phosphatases through substrate profiling

Large-scale activity profiling of enzyme superfamilies provides information about cellular functions as well as the intrinsic binding capabilities of conserved folds. Herein, the functional space of the ubiquitous haloalkanoate dehalogenase superfamily (HADSF) was revealed by screening a customized substrate library against >200 enzymes from representative prokaryotic species, enabling inferred annotation of ∼35% of the HADSF. An extremely high level of substrate ambiguity was revealed, with the majority of HADSF enzymes using more than five substrates. Substrate profiling allowed assignment of function to previously unannotated enzymes with known structure, uncovered potential new pathways, and identified iso-functional orthologs from evolutionarily distant taxonomic groups. Intriguingly, the HADSF subfamily having the least structural elaboration of the Rossmann fold catalytic domain was the most specific, consistent with the concept that domain insertions drive the evolution of new functions and that the broad specificity observed in HADSF may be a relic of this process.

Deletion of PHO13, encoding haloacid dehalogenase type IIA phosphatase, results in upregulation of the pentose phosphate pathway in Saccharomyces cerevisiae

The haloacid dehalogenase (HAD) superfamily is one of the largest enzyme families, consisting mainly of phosphatases. Although intracellular phosphate plays important roles in many cellular activities, the biological functions of HAD enzymes are largely unknown. Pho13 is 1 of 16 putative HAD enzymes in Saccharomyces cerevisiae. Pho13 has not been studied extensively, but previous studies have identified PHO13 to be a deletion target for the generation of industrially attractive phenotypes, namely, efficient xylose fermentation and high tolerance to fermentation inhibitors. In order to understand the molecular mechanisms underlying the improved xylose-fermenting phenotype produced by deletion of PHO13 (pho13Δ), we investigated the response of S. cerevisiae to pho13Δ at the transcriptomic level when cells were grown on glucose or xylose. Transcriptome sequencing analysis revealed that pho13Δ resulted in upregulation of the pentose phosphate (PP) pathway and NADPH-producing enzymes when cells were grown on glucose or xylose. We also found that the transcriptional changes induced by pho13Δ required the transcription factor Stb5, which is activated specifically under NADPH-limiting conditions. Thus, pho13Δ resulted in the upregulation of the PP pathway and NADPH-producing enzymes as a part of an oxidative stress response mediated by activation of Stb5. Because the PP pathway is the primary pathway for xylose, its upregulation by pho13Δ might explain the improved xylose metabolism. These findings will be useful for understanding the biological function of S. cerevisiae Pho13 and the HAD superfamily enzymes and for developing S. cerevisiae strains with industrially attractive phenotypes.

Enzyme promiscuity: engine of evolutionary innovation

Catalytic promiscuity and substrate ambiguity are keys to evolvability, which in turn is pivotal to the successful acquisition of novel biological functions. Action on multiple substrates (substrate ambiguity) can be harnessed for performance of functions in the cell that supersede catalysis of a single metabolite. These functions include proofreading, scavenging of nutrients, removal of antimetabolites, balancing of metabolite pools, and establishing system redundancy. In this review, we present examples of enzymes that perform these cellular roles by leveraging substrate ambiguity and then present the structural features that support both specificity and ambiguity. We focus on the phosphatases of the haloalkanoate dehalogenase superfamily and the thioesterases of the hotdog fold superfamily.

Structure of a His170Tyr mutant of thermostable pNPPase from Geobacillus stearothermophilus

Using directed evolution based on random mutagenesis and heat-treated selection, a thermostable His170Tyr mutant of Geobacillus stearothermophilus thermostable p-nitrophenylphosphatase (TpNPPase) was obtained. The temperature at which the His170Tyr mutant lost 50% of its activity (T1/2) was found to be 4.40 K higher than that of wild-type TpNPPase, and the melting temperature of the His170Tyr mutant increased by 2.39 K. The crystal structure of the His170Tyr mutant was then determined at 2.0 Å resolution in the presence of a sodium ion and a sulfate ion in the active site. The cap domain of chain B shows a half-closed conformation. The hydrophobic side chain of the mutated residue, the hydroxyphenyl group, forms a hydrophobic contact with the methyl group of Ala166. This hydrophobic interaction was found using the Protein Interactions Calculator (PIC) web server with an interaction distance of 4.6 Å, and might be a key factor in the thermostabilization of the His170Tyr mutant. This study potentially offers a molecular basis for both investigation of the catalytic mechanism and thermostable protein engineering.

Structure-guided approach for detecting large domain inserts in protein sequences as illustrated using the haloacid dehalogenase superfamily

In multi-domain proteins, the domains typically run end-to-end, that is, one domain follows the C-terminus of another domain. However, approximately 10% of multi-domain proteins are formed by insertion of one domain sequence into that of another domain. Detecting such insertions within protein sequences is a fundamental challenge in structural biology. The haloacid dehalogenase superfamily (HADSF) serves as a challenging model system wherein a variable cap domain (∼5-200 residues in length) accessorizes the ubiquitous Rossmann-fold core domain, with variations in insertion site and topology corresponding to different classes of cap types. Herein, we describe a comprehensive computational strategy, CapPredictor, for determining large, variable domain insertions in protein sequences. Using a novel sequence-alignment algorithm in conjunction with a structure-guided sequence profile from 154 core-domain-only structures, more than 40,000 HADSF member sequences were assigned cap types. The resulting data set afforded insight into HADSF evolution. Notably, a similar distribution of cap-type classes across different phyla was observed, indicating that all cap types existed in the last universal common ancestor. In addition, comparative analyses of the predicted cap-type and functional assignments showed that different cap types carry out similar chemistries. Thus, while cap domains play a role in substrate recognition and chemical reactivity, cap-type does not strictly define functional class. Through this example, we have shown that CapPredictor is an effective new tool for the study of form and function in protein families where domain insertion occurs.

Chronophin dimerization is required for proper positioning of its substrate specificity loop

Mammalian phosphatases of the haloacid dehalogenase (HAD) superfamily have emerged as important regulators of physiology and disease. Many of these enzymes are stable homodimers; however, the role of their dimerization is largely unknown. Here, we explore the function of the obligatory homodimerization of chronophin, a mammalian HAD phosphatase known to dephosphorylate pyridoxal 5'-phosphate (PLP) and serine/threonine-phosphorylated proteins. The exchange of two residues in the murine chronophin homodimerization interface (chronophin(A194K,A195K)) yields a constitutive monomer both in vitro and in cells. The catalytic activity of monomeric chronophin toward PLP is strongly impaired. X-ray crystallographic studies of chronophin(A194K,A195K) revealed that dimer formation is essential for an intermolecular arginine-arginine-tryptophan stacking interaction that positions a critical histidine residue in the substrate specificity loop of chronophin for PLP coordination. Analysis of all available crystal structures of HAD hydrolases that are grouped together with chronophin in the C2a-type structural subfamily uncovered a highly conserved mode of dimerization that results in intermolecular contacts involving the substrate specificity loop. Our results explain how the dimerization of HAD hydrolases contributes to their catalytic efficiency and substrate specificity.

Evolutionary and structural analyses of mammalian haloacid dehalogenase-type phosphatases AUM and chronophin provide insight into the basis of their different substrate specificities

Mammalian haloacid dehalogenase (HAD)-type phosphatases are an emerging family of phosphatases with important functions in physiology and disease, yet little is known about the basis of their substrate specificity. Here, we characterize a previously unexplored HAD family member (gene annotation, phosphoglycolate phosphatase), which we termed AUM, for aspartate-based, ubiquitous, Mg(2+)-dependent phosphatase. AUM is a tyrosine-specific paralog of the serine/threonine-specific protein and pyridoxal 5'-phosphate-directed HAD phosphatase chronophin. Comparative evolutionary and biochemical analyses reveal that a single, differently conserved residue in the cap domain of either AUM or chronophin is crucial for phosphatase specificity. We have solved the x-ray crystal structure of the AUM cap fused to the catalytic core of chronophin to 2.65 Å resolution and present a detailed view of the catalytic clefts of AUM and chronophin that explains their substrate preferences. Our findings identify a small number of cap domain residues that encode the different substrate specificities of AUM and chronophin.

Consequences of domain insertion on sequence-structure divergence in a superfold

Although the universe of protein structures is vast, these innumerable structures can be categorized into a finite number of folds. New functions commonly evolve by elaboration of existing scaffolds, for example, via domain insertions. Thus, understanding structural diversity of a protein fold evolving via domain insertions is a fundamental challenge. The haloalkanoic dehalogenase superfamily serves as an excellent model system wherein a variable cap domain accessorizes the ubiquitous Rossmann-fold core domain. Here, we determine the impact of the cap-domain insertion on the sequence and structure divergence of the core domain. Through quantitative analysis on a unique dataset of 154 core-domain-only and cap-domain-only structures, basic principles of their evolution have been uncovered. The relationship between sequence and structure divergence of the core domain is shown to be monotonic and independent of the corresponding type of domain insert, reflecting the robustness of the Rossmann fold to mutation. However, core domains with the same cap type share greater similarity at the sequence and structure levels, suggesting interplay between the cap and core domains. Notably, results reveal that the variance in structure maps to α-helices flanking the central β-sheet and not to the domain-domain interface. Collectively, these results hint at intramolecular coevolution where the fold diverges differentially in the context of an accessory domain, a feature that might also apply to other multidomain superfamilies.

Pyrimidine homeostasis is accomplished by directed overflow metabolism

Cellular metabolism converts available nutrients into usable energy and biomass precursors. The process is regulated to facilitate efficient nutrient use and metabolic homeostasis. Feedback inhibition of the first committed step of a pathway by its final product is a classical means of controlling biosynthesis^1–4. In a canonical example, the first committed enzyme in the pyrimidine pathway in Escherichia coli is allosterically inhibited by cytidine triphosphate^1,4,5. The physiological consequences of disrupting this regulation, however, have not been previously explored. Here we identify an alternative regulatory strategy that enables precise control of pyrimidine pathway end-product levels, even in the presence of dysregulated biosynthetic flux. The mechanism involves cooperative feedback regulation of the near-terminal pathway enzyme uridine monophosphate kinase⁶. Such feedback leads to build-up of the pathway intermediate uridine monophosphate, which is in turn degraded by a conserved phosphatase, here termed UmpH, with previously unknown physiological function^7,8. Such directed overflow metabolism allows homeostasis of uridine triphosphate and cytidine triphosphate levels at the expense of uracil excretion and slower growth during energy limitation. Disruption of the directed overflow regulatory mechanism impairs growth in pyrimidine-rich environments. Thus, pyrimidine homeostasis involves dual regulatory strategies, with classical feedback inhibition enhancing metabolic efficiency and directed overflow metabolism ensuring end-product homeostasis.

Discovery of a glycerol 3-phosphate phosphatase reveals glycerophospholipid polar head recycling in Mycobacterium tuberculosis

Functional assignment of enzymes encoded by the Mycobacterium tuberculosis genome is largely incomplete despite recent advances in genomics and bioinformatics. Here, we applied an activity-based metabolomic profiling method to assign function to a unique phosphatase, Rv1692. In contrast to its annotation as a nucleotide phosphatase, metabolomic profiling and kinetic characterization indicate that Rv1692 is a D,L-glycerol 3-phosphate phosphatase. Crystal structures of Rv1692 reveal a unique architecture, a fusion of a predicted haloacid dehalogenase fold with a previously unidentified GCN5-related N-acetyltransferase region. Although not directly involved in acetyl transfer, or regulation of enzymatic activity in vitro, this GCN5-related N-acetyltransferase region is critical for the solubility of the phosphatase. Structural and biochemical analysis shows that the active site features are adapted for recognition of small polyol phosphates, and not nucleotide substrates. Functional assignment and metabolomic studies of M. tuberculosis lacking rv1692 demonstrate that Rv1692 is the final enzyme involved in glycerophospholipid recycling/catabolism, a pathway not previously described in M. tuberculosis.

The use of amino sugars by Bacillus subtilis: presence of a unique operon for the catabolism of glucosamine

B. subtilis grows more rapidly using the amino sugar glucosamine as carbon source, than with N-acetylglucosamine. Genes for the transport and metabolism of N-acetylglucosamine (nagP and nagAB) are found in all the sequenced Bacilli (except Anoxybacillus flavithermus). In B. subtilis there is an additional operon (gamAP) encoding second copies of genes for the transport and catabolism of glucosamine. We have developed a method to make multiple deletion mutations in B. subtilis employing an excisable spectinomycin resistance cassette. Using this method we have analysed the contribution of the different genes of the nag and gam operons for their role in utilization of glucosamine and N-acetylglucosamine. Faster growth on glucosamine is due to the presence of the gamAP operon, which is strongly induced by glucosamine. Although the gamA and nagB genes encode isozymes of GlcN6P deaminase, catabolism of N-acetylglucosamine relies mostly upon the gamA gene product. The genes for use of N-acetylglucosamine, nagAB and nagP, are repressed by YvoA (NagR), a GntR family regulator, whose gene is part of the nagAB yvoA(nagR) operon. The gamAP operon is repressed by YbgA, another GntR family repressor, whose gene is expressed divergently from gamAP. The nagAB yvoA synton is found throughout the Bacilli and most firmicutes. On the other hand the ybgA-gamAP synton, which includes the ybgB gene for a small protein of unknown provenance, is only found in B. subtilis (and a few very close relatives). The origin of ybgBA-gamAP grouping is unknown but synteny analysis suggests lateral transfer from an unidentified donor. The presence of gamAP has enabled B. subtilis to efficiently use glucosamine as carbon source.

HAD hydrolase function unveiled by substrate screening: enzymatic characterization of Arabidopsis thaliana subclass I phosphosugar phosphatase AtSgpp

This work presents the isolation and the biochemical characterization of the Arabidopsis thaliana gene AtSgpp. This gene shows homology with the Arabidopsis low molecular weight phosphatases AtGpp1 and AtGpp2 and the yeast counterpart GPP1 and GPP2, which have a high specificity for DL-glycerol-3-phosphate. In addition, it exhibits homology with DOG1 and DOG2 that dephosphorylate 2-deoxy-D-glucose-6-phosphate. Using a comparative genomic approach, we identified the AtSgpp gene as a conceptual translated haloacid dehalogenase-like hydrolase HAD protein. AtSgpp (locus tag At2g38740), encodes a protein with a predicted Mw of 26.7 kDa and a pI of 4.6. Its sequence motifs and expected structure revealed that AtSgpp belongs to the HAD hydrolases subfamily I, with the C1-type cap domain. In the presence of Mg(2+) ions, the enzyme has a phosphatase activity over a wide range of phosphosugars substrates (pH optima at 7.0 and K m in the range of 3.6-7.7 mM). AtSgpp promiscuity is preferentially detectable on D-ribose-5-phosphate, 2-deoxy-D-ribose-5-phosphate, 2-deoxy-D-glucose-6-phosphate, D-mannose-6-phosphate, D-fructose-1-phosphate, D-glucose-6-phosphate, DL-glycerol-3-phosphate, and D-fructose-6-phosphate, as substrates. AtSgpp is ubiquitously expressed throughout development in most plant organs, mainly in sepal and guard cell. Interestingly, expression is affected by abiotic and biotic stresses, being the greatest under Pi starvation and cyclopentenone oxylipins induction. Based on both, substrate lax specificity and gene expression, the physiological function of AtSgpp in housekeeping detoxification, modulation of sugar-phosphate balance and Pi homeostasis, is provisionally assigned.

Substrate ambiguity among the nudix hydrolases: biologically significant, evolutionary remnant, or both?

Many members of the nudix hydrolase family exhibit considerable substrate multispecificity and ambiguity, which raises significant issues when assessing their functions in vivo and gives rise to errors in database annotation. Several display low antimutator activity when expressed in bacterial tester strains as well as some degree of activity in vitro towards mutagenic, oxidized nucleotides such as 8-oxo-dGTP. However, many of these show greater activity towards other nucleotides such as ADP-ribose or diadenosine tetraphosphate (Ap(4)A). The antimutator activities have tended to gain prominence in the literature, whereas they may in fact represent the residual activity of an ancestral antimutator enzyme that has become secondary to the more recently evolved major activity after gene duplication. Whether any meaningful antimutagenic function has also been retained in vivo requires very careful assessment. Then again, other examples of substrate ambiguity may indicate as yet unexplored regulatory systems. For example, bacterial Ap(4)A hydrolases also efficiently remove pyrophosphate from the 5' termini of mRNAs, suggesting a potential role for Ap(4)A in the control of bacterial mRNA turnover, while the ability of some eukaryotic mRNA decapping enzymes to degrade IDP and dIDP or diphosphoinositol polyphosphates (DIPs) may also be indicative of new regulatory networks in RNA metabolism. DIP phosphohydrolases also degrade diadenosine polyphosphates and inorganic polyphosphates, suggesting further avenues for investigation. This article uses these and other examples to highlight the need for a greater awareness of the possible significance of substrate ambiguity among the nudix hydrolases as well as the need to exert caution when interpreting incomplete analyses.

Identification of a HAD superfamily phosphatase, HdpA, involved in 1,3-dihydroxyacetone production during sugar catabolism in Corynebacterium glutamicum

Corynebacterium glutamicum produces 1,3-dihydroxyacetone (DHA) as metabolite of sugar catabolism but the responsible enzyme is yet to be identified. Using a transposon mutant library, the gene hdpA (cgR_2128) was shown to encode a haloacid dehalogenase superfamily member that catalyzes dephosphorylation of dihydroxyacetone phosphate to produce DHA. Inactivation of hdpA led to a drastic decrease in DHA production from each of glucose, fructose, and sucrose, indicating that HdpA is the main enzyme responsible for DHA production from sugars in C. glutamicum. Confirmation of DHA production via dihydroxyacetone phosphatase finally confirms a long-speculated route through which bacteria produce DHA.

Multi-copy genes that enhance the yield of mammalian G protein-coupled receptors in Escherichia coli

Low yields of recombinant expression represent a major barrier to the physical characterization of membrane proteins. Here, we have identified genes that globally enhance the production of properly folded G protein-coupled receptors (GPCRs) in Escherichia coli. Libraries of bacterial chromosomal fragments were screened using two separate systems that monitor: (i) elevated fluorescence conferred by enhanced expression of GPCR-GFP fusions and (ii) increased binding of fluorescent ligand in cells producing more active receptor. Three multi-copy hits were isolated by both methods: nagD, encoding the ribonucleotide phosphatase NagD; a fragment of nlpD, encoding a truncation of the predicted lipoprotein NlpD, and the three-gene cluster ptsN-yhbJ-npr, encoding three proteins of the nitrogen phosphotransferase system. Expression of these genes resulted in a 3- to 10-fold increase in the yields of different mammalian GPCRs. Our data is consistent with the hypothesis that the expression of these genes may serve to maintain the integrity of the bacterial periplasm and to provide a favorable environment for proper membrane protein folding, possibly by inducing a fine-tuned stress response and/or via modifying the composition of the bacterial cell envelope.

Divergence of structure and function in the haloacid dehalogenase enzyme superfamily: Bacteroides thetaiotaomicron BT2127 is an inorganic pyrophosphatase

The explosion of protein sequence information requires that current strategies for function assignment evolve to complement experimental approaches with computationally based function prediction. This necessitates the development of strategies based on the identification of sequence markers in the form of specificity determinants and a more informed definition of orthologues. Herein, we have undertaken the function assignment of the unknown haloalkanoate dehalogenase superfamily member BT2127 (Uniprot accession code Q8A5 V9) from Bacteroides thetaiotaomicron using an integrated bioinformatics-structure-mechanism approach. The substrate specificity profile and steady-state rate constants of BT2127 (with a k(cat)/K(m) value for pyrophosphate of ~1 × 10(5) M(-1) s(-1)), together with the gene context, support the assigned in vivo function as an inorganic pyrophosphatase. The X-ray structural analysis of wild-type BT2127 and several variants generated by site-directed mutagenesis shows that substrate discrimination is based, in part, on active site space restrictions imposed by the cap domain (specifically by residues Tyr76 and Glu47). Structure-guided site-directed mutagenesis coupled with kinetic analysis of the mutant enzymes identified the residues required for catalysis, substrate binding, and domain-domain association. On the basis of this structure-function analysis, the catalytic residues Asp11, Asp13, Thr113, and Lys147 as well the metal binding residues Asp171, Asn172, and Glu47 were used as markers to confirm BT2127 orthologues identified via sequence searches. This bioinformatic analysis demonstrated that the biological range of BT2127 orthologue is restricted to the phylum Bacteroidetes/Chlorobi. The key structural determinants in the divergence of BT2127 and its closest homologue, β-phosphoglucomutase, control the leaving group size (phosphate vs glucose phosphate) and the position of the Asp acid/base in the open versus closed conformations. HADSF pyrophosphatases represent a third mechanistic and fold type for bacterial pyrophosphatases.

The X-ray crystallographic structure and specificity profile of HAD superfamily phosphohydrolase BT1666: comparison of paralogous functions in B. thetaiotaomicron

Analysis of the haloalkanoate dehalogenase superfamily (HADSF) has uncovered homologues occurring within the same organism that are found to possess broad, overlapping substrate specificities, and low catalytic efficiencies. Here we compare the HADSF phosphatase BT1666 from Bacteroides thetaiotaomicron VPI-5482 to a homologue with high sequence identity (40%) from the same organism BT4131, a known hexose-phosphate phosphatase. The goal is to find whether these enzymes represent duplicated versus paralogous activities. The X-ray crystal structure of BT1666 was determined to 1.82 Å resolution. Superposition of the BT1666 and BT4131 structures revealed a conserved fold and identical active sites suggestive of a common physiological substrate. The steady-state kinetic constants for BT1666 were determined for a diverse panel of phosphorylated metabolites to define its substrate specificity profile and overall level of catalytic efficiency. Whereas BT1666 and BT4131 are both promiscuous, their substrate specificity profiles are distinct. The catalytic efficiency of BT1666 (k(cat) /K(m) = 4.4 × 10(2) M(-1) s(-1) for the best substrate fructose 1,6-(bis)phosphate) is an order of magnitude less than that of BT4131 (k(cat) /K(m) = 6.7 × 10(3) M(-1) s(-1) for 2-deoxyglucose 6-phosphate). The seemingly identical active-site structures point to sequence variation outside the active site causing differences in conformational dynamics or subtle catalytic positioning effects that drive the divergence in catalytic efficiency and selectivity. The overlapping substrate profiles may be understood in terms of differential regulation of expression of the two enzymes or a conferred advantage in metabolic housekeeping functions by having a larger range of possible metabolites as substrates.

Characterization and regulation of a bacterial sugar phosphatase of the haloalkanoate dehalogenase superfamily, AraL, from Bacillus subtilis

AraL from Bacillus subtilis is a member of the ubiquitous haloalkanoate dehalogenase superfamily. The araL gene has been cloned, over-expressed in Escherichia coli and its product purified to homogeneity. The enzyme displays phosphatase activity, which is optimal at neutral pH (7.0) and 65 °C. Substrate screening and kinetic analysis showed AraL to have low specificity and catalytic activity towards several sugar phosphates, which are metabolic intermediates of the glycolytic and pentose phosphate pathways. On the basis of substrate specificity and gene context within the arabinose metabolic operon, a putative physiological role of AraL in the detoxification of accidental accumulation of phosphorylated metabolites has been proposed. The ability of AraL to catabolize several related secondary metabolites requires regulation at the genetic level. In the present study, using site-directed mutagenesis, we show that the production of AraL is regulated by a structure in the translation initiation region of the mRNA, which most probably blocks access to the ribosome-binding site, preventing protein synthesis. Members of haloalkanoate dehalogenase subfamily IIA and IIB are characterized by a broad-range and overlapping specificity anticipating the need for regulation at the genetic level. We provide evidence for the existence of a genetic regulatory mechanism controlling the production of AraL.

Enzyme promiscuity: a mechanistic and evolutionary perspective

Many, if not most, enzymes can promiscuously catalyze reactions, or act on substrates, other than those for which they evolved. Here, we discuss the structural, mechanistic, and evolutionary implications of this manifestation of infidelity of molecular recognition. We define promiscuity and related phenomena and also address their generality and physiological implications. We discuss the mechanistic enzymology of promiscuity--how enzymes, which generally exert exquisite specificity, catalyze other, and sometimes barely related, reactions. Finally, we address the hypothesis that promiscuous enzymatic activities serve as evolutionary starting points and highlight the unique evolutionary features of promiscuous enzyme functions.

Divergence of biochemical function in the HAD superfamily: D-glycero-D-manno-heptose-1,7-bisphosphate phosphatase (GmhB)

D-Glycero-d-manno-heptose-1,7-bisphosphate phosphatase (GmhB) is a member of the histidinol-phosphate phosphatase (HisB) subfamily of the haloalkanoic acid dehalogenase (HAD) enzyme superfamily. GmhB supports two divergent biochemical pathways in bacteria: the d-glycero-d-manno-heptose-1alpha-GDP pathway (in S-layer glycoprotein biosynthesis) and the l-glycero-d-manno-heptose-1beta-ADP pathway (in lipid A biosynthesis). Herein, we report the comparative analysis of substrate recognition in selected GmhB orthologs. The substrate specificity of the l-glycero-d-manno-heptose-1beta-ADP pathway GmhB from Escherichia coli K-12 was evaluated using hexose and heptose bisphosphates, histidinol phosphate, and common organophosphate metabolites. Only d-glycero-d-manno-heptose 1beta,7-bisphosphate (k(cat)/K(m) = 7 x 10(6) M(-1) s(-1)) and d-glycero-d-manno-heptose 1alpha,7-bisphosphate (k(cat)/K(m) = 7 x 10(4) M(-1) s(-1)) displayed physiologically significant substrate activity. (31)P NMR analysis demonstrated that E. coli GmhB selectively removes the C(7) phosphate. Steady-state kinetic inhibition studies showed that d-glycero-d-manno-heptose 1beta-phosphate (K(is) = 60 microM, and K(ii) = 150 microM) and histidinol phosphate (K(is) = 1 mM, and K(ii) = 6 mM), while not hydrolyzed, do in fact bind to E. coli GmhB, which leads to the conclusion that nonproductive binding contributes to substrate discrimination. High catalytic efficiency and a narrow substrate range are characteristic of a well-evolved metabolic enzyme, and as such, E. coli GmhB is set apart from most HAD phosphatases (which are typically inefficient and promiscuous). The specialization of the biochemical function of GmhB was examined by measuring the kinetic constants for hydrolysis of the alpha- and beta-anomers of d-glycero-d-manno-heptose 1beta,7-bisphosphate catalyzed by the GmhB orthologs of the l-glycero-d-manno-heptose 1beta-ADP pathways operative in Bordetella bronchiseptica and Mesorhizobium loti and by the GmhB of the d-glycero-d-manno-heptose 1alpha-GDP pathway operative in Bacteroides thetaiotaomicron. The results show that although each of these representatives possesses physiologically significant catalytic activity toward both anomers, each displays substantial anomeric specificity. Like E. coli GmhB, B. bronchiseptica GmhB and M. loti GmhB prefer the beta-anomer, whereas B. thetaiotaomicron GmhB is selective for the alpha-anomer. By determining the anomeric configuration of the physiological substrate (d-glycero-d-manno-heptose 1,7-bisphosphate) for each of the four GmhB orthologs, we discovered that the anomeric specificity of GmhB correlates with that of the pathway kinase. The conclusion drawn from this finding is that the evolution of the ancestor to GmhB in the HisB subfamily provided for specialization toward two distinct biochemical functions.

Biosynthesis of ethers: unusual or common natural events?

Ether bonds are found in a wide variety of natural products--mainly secondary metabolites--including lipids, oxiranes, terpenoids, flavonoids, polyketides, and carbohydrate derivatives, to name some representative examples. To furnish such a biodiversity of structures, a large number of different enzymes are involved in several different biosynthetic pathways. Depending on the compound and on the (micro) environment in which the reaction is performed, ethers are produced by very different (enzymatic) reactions, thus providing an impressive display of how Nature has combined evolution and thermodynamics to be able to produce a vast number of compounds. In addition, many of these compounds possess different biological activities of pharmacological interest. Moreover, some of these ethers (i.e., epoxides) have high chemical reactivity, and can be useful starting materials for further synthetic processes. This review aims to provide an overview of the different strategies that are found in Nature for the formation of these "bioethers". Both fundamental and practical insights of the biosynthetic processes will be discussed.

Prediction of function in protein superfamilies

Assignment of function for enzymes encoded in sequenced genomes is a challenging task. Predictions of enzyme function can be made using clues from superfamily assignment, structure, genome context, phylogenetic conservation, and virtual screening to identify potential ligands. Ultimately, confident assignment of function requires experimental verification as well as an understanding of the physiological role of an enzyme in the context of the metabolic network.

Markers of fitness in a successful enzyme superfamily

Haloalkanoic acid dehalogenase (HAD) superfamily members serve as the predominant catalysts of metabolic phosphate ester hydrolysis in all three superkingdoms of life. Collectively, the known structural, bioinformatic, and mechanistic data offer a glimpse of the variety of HAD enzymes that have evolved in the service of metabolic expansion. Factors that have contributed to superfamily dominance include a chemically versatile nucleophile, stability of the core superfold, structural modularity of the chemistry and specificity domains, conformational coupling conferred by the topology of the inserted specificity elements, and retention of a conserved mold for stabilization of the trigonal bipyramidal transition state.

Allosteric regulation of glucosamine-6-phosphate deaminase (NagB) and growth of Escherichia coli on glucosamine

Growth on N-acetylglucosamine (GlcNAc) produces intracellular N-acetylglucosamine-6-phosphate (GlcNAc6P), which affects the regulation of the catabolism of amino sugars in Escherichia coli in two ways. First, GlcNAc6P is the inducing signal for the NagC repressor, and thus it increases the expression of the enzymes of the nagE-nagBACD operon. Second, it is the allosteric activator of glucosamine-6P (GlcN6P) deaminase, NagB, and thus increases the catalytic capacity of this key enzyme in the metabolism of amino sugars. We showed previously that both the level of expression of the nagB gene and the transport of glucosamine were limiting the growth rate on GlcN (L. I. Alvarez-Añorve et al., J. Bacteriol. 187:2974-2982, 2005). We were unable to conclude if the lack of allosteric activation of wild-type NagB was also contributing to the slower growth rate on GlcN. Using a single-copy plasmid, with a constitutive promoter, we have separated the effects of GlcNAc6P on the NagB protein level and on deaminase activity. We show that over a range of intracellular NagB concentrations it is the quantity of the substrate, GlcN6P, which is limiting growth rather than the concentration of the allosteric activator, GlcNAc6P. On the other hand, the F174A mutant of NagB, which requires higher concentrations of GlcNAc6P for activity in vitro, grew better on GlcN in the presence of GlcNAc6P. However, wild-type NagB behaves as if it is already fully allosterically activated during growth on GlcN, and we present evidence suggesting that sufficient GlcNAc6P for allosteric activation is derived from the recycling of peptidoglycan.

Discovery of a dipeptide epimerase enzymatic function guided by homology modeling and virtual screening

We have developed a computational approach to aid the assignment of enzymatic function for uncharacterized proteins that uses homology modeling to predict the structure of the binding site and in silico docking to identify potential substrates. We apply this method to proteins in the functionally diverse enolase superfamily that are homologous to the characterized L-Ala-D/L-Glu epimerase from Bacillus subtilis. In particular, a protein from Thermotoga martima was predicted to have different substrate specificity, which suggests that it has a different, but as yet unknown, biological function. This prediction was experimentally confirmed, resulting in the assignment of epimerase activity for L-Ala-D/L-Phe, L-Ala-D/L-Tyr, and L-Ala-D/L-His, whereas the enzyme is annotated incorrectly in GenBank as muconate cycloisomerase. Subsequently, crystal structures of the enzyme were determined in complex with three substrates, showing close agreement with the computational models and revealing the structural basis for the observed substrate selectivity.

Human symbiont Bacteroides thetaiotaomicron synthesizes 2-keto-3-deoxy-D-glycero-D- galacto-nononic acid (KDN)

The proper functioning of the human intestine is dependent on its bacterial symbionts, the most predominant of which belong to the Phylum Bacteroidetes. These bacteria are known to use variable displays of multiple capsular polysaccharides (CPs) to aid in their survival and foraging within the intestine. Bacteroides thetaiotaomicron is a prominent human gut symbiont and a remarkably versatile glycophile. The structure determination of the CPs, encoded by the eight CP loci, is the key to understanding the mechanism of this organism's adaptation on a molecular level. Herein, we report the bioinformatics-based discovery and chemical demonstration of a biosynthetic pathway that forms and cytidylates 2-keto-3-deoxy-D-glycero-D-galacto-nononic acid (KDN), most likely for inclusion in the CP encoded by B. thetaiotaomicron CP locus 7.

Structure-function analysis of 2-keto-3-deoxy-D-glycero-D-galactonononate-9-phosphate phosphatase defines specificity elements in type C0 haloalkanoate dehalogenase family members

The phosphotransferases of the haloalkanoate dehalogenase superfamily (HADSF) act upon a wide range of metabolites in all eukaryotes and prokaryotes and thus constitute a significant force in cell function. The challenge posed for biochemical function assignment of HADSF members is the identification of the structural determinants that target a specific metabolite. The "8KDOP" subfamily of the HADSF is defined by the known structure and catalytic activity of 2-keto-3-deoxy-8-phospho-d-manno-octulosonic acid (KDO-8-P) phosphatase. Homologues of this enzyme have been uniformly annotated as KDO-8-P phosphatase. One such gene, BT1713, from the Bacteroides thetaiotaomicron genome was recently found to encode the enzyme 2-keto-3-deoxy-d-glycero-d-galacto-9-phosphonononic acid (KDN-9-P) phosphatase in the biosynthetic pathway of the 9-carbon alpha-keto acid, 2-keto-3-deoxy-d-glycero-d-galactonononic acid (KDN). To find the structural elements that provide substrate-specific interactions and to allow identification of genomic sequence markers, the x-ray crystal structures of BT1713 liganded to the cofactor Mg(2+)and complexed with tungstate or VO(3)(-)/Neu5Ac were determined to 1.1, 1.85, and 1.63 A resolution, respectively. The structures define the active site to be at the subunit interface and, as confirmed by steady-state kinetics and site-directed mutagenesis, reveal Arg-64(*), Lys-67(*), and Glu-56 to be the key residues involved in sugar binding that are essential for BT1713 catalytic function. Bioinformatic analyses of the differentially conserved residues between BT1713 and KDO-8-P phosphatase homologues guided by the knowledge of the structure-based specificity determinants define Glu-56 and Lys-67(*) to be the key residues that can be used in future annotations.

The catalytic scaffold of the haloalkanoic acid dehalogenase enzyme superfamily acts as a mold for the trigonal bipyramidal transition state

The evolution of new catalytic activities and specificities within an enzyme superfamily requires the exploration of sequence space for adaptation to a new substrate with retention of those elements required to stabilize key intermediates/transition states. Here, we propose that core residues in the large enzyme family, the haloalkanoic acid dehalogenase enzyme superfamily (HADSF) form a "mold" in which the trigonal bipyramidal transition states formed during phosphoryl transfer are stabilized by electrostatic forces. The vanadate complex of the hexose phosphate phosphatase BT4131 from Bacteroides thetaiotaomicron VPI-5482 (HPP) determined at 1.00 A resolution via X-ray crystallography assumes a trigonal bipyramidal coordination geometry with the nucleophilic Asp-8 and one oxygen ligand at the apical position. Remarkably, the tungstate in the complex determined to 1.03 A resolution assumes the same coordination geometry. The contribution of the general acid/base residue Asp-10 in the stabilization of the trigonal bipyramidal species via hydrogen-bond formation with the apical oxygen atom is evidenced by the 1.52 A structure of the D10A mutant bound to vanadate. This structure shows a collapse of the trigonal bipyramidal geometry with displacement of the water molecule formerly occupying the apical position. Furthermore, the 1.07 A resolution structure of the D10A mutant complexed with tungstate shows the tungstate to be in a typical "phosphate-like" tetrahedral configuration. The analysis of 12 liganded HADSF structures deposited in the protein data bank (PDB) identified stringently conserved elements that stabilize the trigonal bipyramidal transition states by engaging in favorable electrostatic interactions with the axial and equatorial atoms of the transferring phosphoryl group.

The X-ray crystallographic structure and activity analysis of a Pseudomonas-specific subfamily of the HAD enzyme superfamily evidences a novel biochemical function

The haloacid dehalogenase (HAD) superfamily is a large family of proteins dominated by phosphotransferases. Thirty-three sequence families within the HAD superfamily (HADSF) have been identified to assist in function assignment. One such family includes the enzyme phosphoacetaldehyde hydrolase (phosphonatase). Phosphonatase possesses the conserved Rossmanniod core domain and a C1-type cap domain. Other members of this family do not possess a cap domain and because the cap domain of phosphonatase plays an important role in active site desolvation and catalysis, the function of the capless family members must be unique. A representative of the capless subfamily, PSPTO_2114, from the plant pathogen Pseudomonas syringae, was targeted for catalytic activity and structure analyses. The X-ray structure of PSPTO_2114 reveals a capless homodimer that conserves some but not all of the intersubunit contacts contributed by the core domains of the phosphonatase homodimer. The region of the PSPTO_2114 that corresponds to the catalytic scaffold of phosphonatase (and other HAD phosphotransfereases) positions amino acid residues that are ill suited for Mg+2 cofactor binding and mediation of phosphoryl group transfer between donor and acceptor substrates. The absence of phosphotransferase activity in PSPTO_2114 was confirmed by kinetic assays. To explore PSPTO_2114 function, the conservation of sequence motifs extending outside of the HADSF catalytic scaffold was examined. The stringently conserved residues among PSPTO_2114 homologs were mapped onto the PSPTO_2114 three-dimensional structure to identify a surface region unique to the family members that do not possess a cap domain. The hypothesis that this region is used in protein-protein recognition is explored to define, for the first time, HADSF proteins which have acquired a function other than that of a catalyst.

Structural and mutational analyses reveal the functional role of active-site Lys-154 and Asp-173 of Salmonella typhimurium AphA protein

The Salmonella typhimurium class B nonspecific acid phosphatase (AphA protein) belongs to the L2-haloacid dehalogenase superfamily. The conserved Lys-154 interacts with substrate phosphate, nucleophile Asp-46, and Asp-173 in the wild-type AphA protein. Asp-173 also interacts with Mg(II) water ligand and with main-chain amide of loop-4. We report here the mutational analysis of Lys-154 and Asp-173, the crystal structures of the K154N and K154R mutants, and the results of electrostatic potential calculations. The K154N, K154R and D173N mutants display significant reduction in the phosphatase activity. Lys-154 may not be responsible for a juxtaposition of the substrate phosphate and the aspartyl nucleophile, but has an hitherto unknown functional role of rendering the substrate phosphorous atom electron deficient. Nearly 10,000-fold increase in the K(d) value for dissociation of the cofactor Mg(II) observed for the D173N mutant correlates well with theoretically estimated change in the binding free energy of Mg(II).

Diversification of function in the haloacid dehalogenase enzyme superfamily: The role of the cap domain in hydrolytic phosphoruscarbon bond cleavage

Phosphonatase functions in the 2-aminoethylphosphonate (AEP) degradation pathway of bacteria, catalyzing the hydrolysis of the C-P bond in phosphonoacetaldehyde (Pald) via formation of a bi-covalent Lys53ethylenamine/Asp12 aspartylphosphate intermediate. Because phosphonatase is a member of the haloacid dehalogenase superfamily, a family predominantly comprised of phosphatases, the question arises as to how this new catalytic activity evolved. The source of general acid-base catalysis for Schiff-base formation and aspartylphosphate hydrolysis was probed using pH-rate profile analysis of active-site mutants and X-ray crystallographic analysis of modified forms of the enzyme. The 2.9 A X-ray crystal structure of the mutant Lys53Arg complexed with Mg2+ and phosphate shows that the equilibrium between the open and the closed conformation is disrupted, favoring the open conformation. Thus, proton dissociation from the cap domain Lys53 is required for cap domain-core domain closure. The likely recipient of the Lys53 proton is a water-His56 pair that serves to relay the proton to the carbonyl oxygen of the phosphonoacetaldehyde (Pald) substrate upon addition of the Lys53. The pH-rate profile analysis of active-site mutants was carried out to test this proposal. The proximal core domain residues Cys22 and Tyr128 were ruled out, and the role of cap domain His56 was supported by the results. The X-ray crystallographic structure of wild-type phosphonatase reduced with NaBH4 in the presence of Pald was determined at 2.4A resolution to reveal N epsilon-ethyl-Lys53 juxtaposed with a sulfate ligand bound in the phosphate site. The position of the C2 of the N-ethyl group in this structure is consistent with the hypothesis that the cap domain N epsilon-ethylenamine-Lys53 functions as a general base in the hydrolysis of the aspartylphosphate bi-covalent enzyme intermediate. Because the enzyme residues proposed to play a key role in P-C bond cleavage are localized on the cap domain, this domain appears to have evolved to support the diversification of the HAD phosphatase core domain for catalysis of hydrolytic P-C bond cleavage.

Structure of Leishmania mexicana Phosphomannomutase Highlights Similarities with Human Isoforms

Phosphomannomutase (PMM) catalyses the conversion of mannose-6-phosphate to mannose-1-phosphate, an essential step in mannose activation and the biosynthesis of glycoconjugates in all eukaryotes. Deletion of PMM from Leishmania mexicana results in loss of virulence, suggesting that PMM is a promising drug target for the development of anti-leishmanial inhibitors. We report the crystallization and structure determination to 2.1 A of L. mexicana PMM alone and in complex with glucose-1,6-bisphosphate to 2.9 A. PMM is a member of the haloacid dehalogenase (HAD) family, but has a novel dimeric structure and a distinct cap domain of unique topology. Although the structure is novel within the HAD family, the leishmanial enzyme shows a high degree of similarity with its human isoforms. We have generated L. major PMM knockouts, which are avirulent. We expressed the human pmm2 gene in the Leishmania PMM knockout, but despite the similarity between Leishmania and human PMM, expression of the human gene did not restore virulence. Similarities in the structure of the parasite enzyme and its human isoforms suggest that the development of parasite-selective inhibitors will not be an easy task.