A subsystems-based approach to the identification of drug targets in bacterial pathogens

Fluorometric Analysis of Carrier-Protein-Dependent Biosynthesis through a Conformationally Sensitive Solvatochromic Pantetheinamide Probe

Carrier proteins (CPs) play a fundamental role in the biosynthesis of fatty acids, polyketides, and non-ribosomal peptides, encompassing many medicinally and pharmacologically relevant compounds. Current approaches to analyze novel carrier-protein-dependent synthetic pathways are hampered by a lack of activity-based assays for natural product biosynthesis. To fill this gap, we turned to 3-methoxychromones, highly solvatochromic fluorescent molecules whose emission intensity and wavelength are heavily dependent on their immediate molecular environment. We have developed a solvatochromic carrier-protein-targeting probe which is able to selectively fluoresce when bound to a target carrier protein. Additionally, the probe displays distinct responses upon CP binding in carrier-protein-dependent synthases. This discerning approach demonstrates the design of solvatochromic fluorophores with the ability to identify biosynthetically active CP-enzyme interactions.

Novel Antimycobacterial Compounds Suppress NAD Biogenesis by Targeting a Unique Pocket of NaMN Adenylyltransferase

Conventional treatments to combat the tuberculosis (TB) epidemic are falling short, thus encouraging the search for novel antitubercular drugs acting on unexplored molecular targets. Several whole-cell phenotypic screenings have delivered bioactive compounds with potent antitubercular activity. However, their cellular target and mechanism of action remain largely unknown. Further evaluation of these compounds may include their screening in search for known antitubercular drug targets hits. Here, a collection of nearly 1400 mycobactericidal compounds was screened against Mycobacterium tuberculosis NaMN adenylyltransferase ( MtNadD), a key enzyme in the biogenesis of NAD cofactor that was recently validated as a new drug target for dormant and active tuberculosis. We found three chemotypes that efficiently inhibit MtNadD in the low micromolar range in vitro. SAR and cheminformatics studies of commercially available analogues point to a series of benzimidazolium derivatives, here named N2, with bactericidal activity on different mycobacteria, including M. abscessus, multidrug-resistant M. tuberculosis, and dormant M. smegmatis. The on-target activity was supported by the increased resistance of an M. smegmatis strain overexpressing the target and by a rapid decline in NAD(H) levels. A cocrystal structure of MtNadD with N2-8 inhibitor reveals that the binding of the inhibitor induced the formation of a new quaternary structure, a dimer-of-dimers where two copies of the inhibitor occupy symmetrical positions in the dimer interface, thus paving the way for the development of a new generation of selective MtNadD bioactive inhibitors. All these results strongly suggest that pharmacological inhibition of MtNadD is an effective strategy to combat dormant and resistant Mtb strains.

Umbrella Sampling and X-ray Crystallographic Analysis Unveil an Arg-Asp Gate Facilitating Inhibitor Binding Inside Phosphopantetheine Adenylyltransferase Allosteric Cleft

Phosphopantetheine adenylyltransferase (PPAT) is a rate-limiting enzyme essential for biosynthesis of coenzyme A (CoA), which in turn is responsible to regulate the secretion of exotoxins via type III secretion system in Pseudomonas aeruginosa, causing severe health concerns ranging from nosocomial infections to respiratory failure. Acetyl coenzyme A (AcCoA) is a newly reported inhibitor of PPAT, believed to regulate the cellular levels of CoA and thereby the pathogenesis. Very little is known so far regarding the mechanistic details of AcCoA binding inside PPAT-binding cleft. Herein, we have used extensive umbrella sampling simulations to decipher mechanistic insight into the inhibitor accommodation inside the binding cavity. We found that R90 and D94 residues act like a gate near the binding cavity to accommodate and stabilize the incoming ligand. Mutational models concerning these residues also show considerable difference in AcCoA-binding thermodynamics. To substantiate our findings, we have solved the first crystal structure of apo-PPAT from P. aeruginosa, which also found to agree with the simulation results. Collectively, these results describe the mechanistic details of accommodation of inhibitor molecule inside PPAT-binding cavity and also offer valuable insight into regulating cellular levels of CoA/AcCoA and thus controlling the pathogenicity.

Analysis of drug binding pockets and repurposing opportunities for twelve essential enzymes of ESKAPE pathogens

Background

The rapid increase in antibiotic resistance by various bacterial pathogens underlies the significance of developing new therapies and exploring different drug targets. A fraction of bacterial pathogens abbreviated as ESKAPE by the European Center for Disease Prevention and Control have been considered a major threat due to the rise in nosocomial infections. Here, we compared putative drug binding pockets of twelve essential and mostly conserved metabolic enzymes in numerous bacterial pathogens including those of the ESKAPE group and Mycobacterium tuberculosis. The comparative analysis will provide guidelines for the likelihood of transferability of the inhibitors from one species to another.

Methods

Nine bacterial species including six ESKAPE pathogens, Mycobacterium tuberculosis along with Mycobacterium smegmatis and Eschershia coli, two non-pathogenic bacteria, have been selected for drug binding pocket analysis of twelve essential enzymes. The amino acid sequences were obtained from Uniprot, aligned using ICM v3.8-4a and matched against the Pocketome encyclopedia. We used known co-crystal structures of selected target enzyme orthologs to evaluate the location of their active sites and binding pockets and to calculate a matrix of pairwise sequence identities across each target enzyme across the different species. This was used to generate sequence maps.

Results

High sequence identity of enzyme binding pockets, derived from experimentally determined co-crystallized structures, was observed among various species. Comparison at both full sequence level and for drug binding pockets of key metabolic enzymes showed that binding pockets are highly conserved (sequence similarity up to 100%) among various ESKAPE pathogens as well as Mycobacterium tuberculosis. Enzymes orthologs having conserved binding sites may have potential to interact with inhibitors in similar way and might be helpful for design of similar class of inhibitors for a particular species. The derived pocket alignments and distance-based maps provide guidelines for drug discovery and repurposing. In addition they also provide recommendations for the relevant model bacteria that may be used for initial drug testing.

Discussion

Comparing ligand binding sites through sequence identity calculation could be an effective approach to identify conserved orthologs as drug binding pockets have shown higher level of conservation among various species. By using this approach we could avoid the problems associated with full sequence comparison. We identified essential metabolic enzymes among ESKAPE pathogens that share high sequence identity in their putative drug binding pockets (up to 100%), of which known inhibitors can potentially antagonize these identical pockets in the various species in a similar manner.

Biogenesis and Homeostasis of Nicotinamide Adenine Dinucleotide Cofactor

Universal and ubiquitous redox cofactors, nicotinamide adenine dinucleotide (NAD) and its phosphorylated analog (NADP), collectively contribute to approximately 12% of all biochemical reactions included in the metabolic model of Escherichia coli K-12. A homeostasis of the NAD pool faithfully maintained by the cells results from a dynamic balance in a network of NAD biosynthesis, utilization, decomposition, and recycling pathways that is subject to tight regulation at various levels. A brief overview of NAD utilization processes is provided in this review, including some examples of nonredox utilization. The review focuses mostly on those aspects of NAD biogenesis and utilization in E. coli and Salmonella that emerged within the past 12 years. The first pyridine nucleotide cycle (PNC) originally identified in mammalian systems and termed the Preiss-Handler pathway includes a single-step conversion of niacin (Na) to NaMN by nicotinic acid phosphoribosyltransferase (PncB). In E. coli and many other prokaryotes, this enzyme, together with nicotinamide deamidase (PncA), compose the major pathway for utilization of the pyridine ring in the form of amidated (Nm) or deamidated (Na) precursors. The existence of various regulatory mechanisms and checkpoints that control the NAD biosynthetic machinery reflects the importance of maintaining NAD homeostasis in a variety of growth conditions. Among the most important regulatory mechanisms at the level of individual enzymes are a classic feedback inhibition of NadB, the first enzyme of NAD de novo biosynthesis, by NAD and a metabolic regulation of NadK by reduced cofactors.

Cloning, expression, purification, crystallization and preliminary X-ray diffraction studies of NAD synthetase from methicillin-resistant Staphylococcus aureus

Staphylococcus aureus is an important human and animal pathogen that causes a wide range of infections. The prevalence of multidrug-resistant S. aureus strains in both hospital and community settings makes it imperative to characterize new drug targets to combat S. aureus infections. In this context, enzymes involved in NAD metabolism and synthesis are significant drug targets as NAD is a central player in several cellular processes. NAD synthetase catalyzes the last step in the biosynthesis of nicotinamide adenine dinucleotide, making it a crucial intermediate enzyme linked to the biosynthesis of several amino acids, purine and pyrimidine nucleotides, coenzymes and antibiotics.

Pyrazinamide: the importance of uncovering the mechanisms of action in mycobacteria

Pyrazinamide (PZA) is still one of the key drugs used in current therapeutic regimens for tuberculosis (TB). Despite its importance for TB therapy, the mode of action of PZA remains unknown. PZA has to be converted to its active form pyrazinoic acid (POA) by the nicotinamidase PncA and is then excreted by an unknown efflux pump. At acidic conditions, POA is protonated to HPOA and is reabsorbed into the cell where it causes cellular damage. For a long time, it has been thought that PZA/POA has no defined target of action, but recent studies have shown that both PZA and POA have several different targets interfering with diverse biochemical pathways, especially in the NAD(+) and energy metabolism. PZA resistance seems to depend not only on a defective pyrazinamidase but is also rather a result of the interplay of many different enzyme targets and transport mechanisms.

Metabolic profiling of alternative NAD biosynthetic routes in mouse tissues

NAD plays essential redox and non-redox roles in cell biology. In mammals, its de novo and recycling biosynthetic pathways encompass two independent branches, the "amidated" and "deamidated" routes. Here we focused on the indispensable enzymes gating these two routes, i.e. nicotinamide mononucleotide adenylyltransferase (NMNAT), which in mammals comprises three distinct isozymes, and NAD synthetase (NADS). First, we measured the in vitro activity of the enzymes, and the levels of all their substrates and products in a number of tissues from the C57BL/6 mouse. Second, from these data, we derived in vivo estimates of enzymes'rates and quantitative contributions to NAD homeostasis. The NMNAT activity, mainly represented by nuclear NMNAT1, appears to be high and nonrate-limiting in all examined tissues, except in blood. The NADS activity, however, appears rate-limiting in lung and skeletal muscle, where its undetectable levels parallel a relative accumulation of the enzyme's substrate NaAD (nicotinic acid adenine dinucleotide). In all tissues, the amidated NAD route was predominant, displaying highest rates in liver and kidney, and lowest in blood. In contrast, the minor deamidated route showed higher relative proportions in blood and small intestine, and higher absolute values in liver and small intestine. Such results provide the first comprehensive picture of the balance of the two alternative NAD biosynthetic routes in different mammalian tissues under physiological conditions. This fills a gap in the current knowledge of NAD biosynthesis, and provides a crucial information for the study of NAD metabolism and its role in disease.

Metabolic and bactericidal effects of targeted suppression of NadD and NadE enzymes in mycobacteria

Mycobacterium tuberculosis remains a major cause of death due to the lack of treatment accessibility, HIV coinfection, and drug resistance. Development of new drugs targeting previously unexplored pathways is essential to shorten treatment time and eliminate persistent M. tuberculosis. A promising biochemical pathway which may be targeted to kill both replicating and nonreplicating M. tuberculosis is the biosynthesis of NAD(H), an essential cofactor in multiple reactions crucial for respiration, redox balance, and biosynthesis of major building blocks. NaMN adenylyltransferase (NadD) and NAD synthetase (NadE), the key enzymes of NAD biosynthesis, were selected as promising candidate drug targets for M. tuberculosis. Here we report for the first time kinetic characterization of the recombinant purified NadD enzyme, setting the stage for its structural analysis and inhibitor development. A protein knockdown approach was applied to validate bothNadD and NadE as target enzymes. Induced degradation of either target enzyme showed a strong bactericidal effect which coincided with anticipated changes in relative levels of NaMN and NaAD intermediates (substrates of NadD and NadE, respectively) and ultimate depletion of the NAD(H) pool. A metabolic catastrophe predicted as a likely result of NAD(H) deprivation of cellular metabolism was confirmed by (13)C biosynthetic labeling followed by gas chromatography-mass spectrometry (GC-MS) analysis. A sharp suppression of metabolic flux was observed in multiple NAD(P)(H)-dependent pathways, including synthesis of many amino acids (serine, proline, aromatic amino acids) and fatty acids. Overall, these results provide strong validation of the essential NAD biosynthetic enzymes, NadD and NadE, as antimycobacterial drug targets.

IMPORTANCE

To address the problems of M. tuberculosis drug resistance and persistence of tuberculosis, new classes of drug targets need to be explored. The biogenesis of NAD cofactors was selected for target validation because of their indispensable role in driving hundreds of biochemical transformations. We hypothesized that the disruption of NAD production in the cell via genetic suppression of the essential enzymes (NadD and NadE) involved in the last two steps of NAD biogenesis would lead to cell death, even under dormancy conditions. In this study, we confirmed the hypothesis using a protein knockdown approach in the model system of Mycobacterium smegmatis. We showed that induced proteolytic degradation of either target enzyme leads to depletion of the NAD cofactor pool, which suppresses metabolic flux through numerous NAD(P)-dependent pathways of central metabolism of carbon and energy production. Remarkably, bactericidal effect was observed even for nondividing bacteria cultivated under carbon starvation conditions.

A genetic strategy to identify targets for the development of drugs that prevent bacterial persistence

Antibacterial drug development suffers from a paucity of targets whose inhibition kills replicating and nonreplicating bacteria. The latter include phenotypically dormant cells, known as persisters, which are tolerant to many antibiotics and often contribute to failure in the treatment of chronic infections. This is nowhere more apparent than in tuberculosis caused by Mycobacterium tuberculosis, a pathogen that tolerates many antibiotics once it ceases to replicate. We developed a strategy to identify proteins that Mycobacterium tuberculosis requires to both grow and persist and whose inhibition has the potential to prevent drug tolerance and persister formation. This strategy is based on a tunable dual-control genetic switch that provides a regulatory range spanning three orders of magnitude, quickly depletes proteins in both replicating and nonreplicating mycobacteria, and exhibits increased robustness to phenotypic reversion. Using this switch, we demonstrated that depletion of the nicotinamide adenine dinucleotide synthetase (NadE) rapidly killed Mycobacterium tuberculosis under conditions of standard growth and nonreplicative persistence induced by oxygen and nutrient limitation as well as during the acute and chronic phases of infection in mice. These findings establish the dual-control switch as a robust tool with which to probe the essentiality of Mycobacterium tuberculosis proteins under different conditions, including those that induce antibiotic tolerance, and NadE as a target with the potential to shorten current tuberculosis chemotherapies.

Quinolinate salvage and insights for targeting NAD biosynthesis in group A streptococci

The essential coenzyme NAD plays important roles in metabolic reactions and cell regulation in all organisms. As such, NAD synthesis has been investigated as a source for novel antibacterial targets. Cross-species genomics-based reconstructions of NAD metabolism in group A streptococci (GAS), combined with focused experimental testing in Streptococcus pyogenes, led to a better understanding of NAD metabolism in the pathogen. The predicted niacin auxotrophy was experimentally verified, as well as the essential role of the nicotinamidase PncA in the utilization of nicotinamide (Nm). PncA is dispensable in the presence of nicotinate (Na), ruling it out as a viable antibacterial target. The function of the "orphan" NadC enzyme, which is uniquely present in all GAS species despite the absence of other genes of NAD de novo synthesis, was elucidated. Indeed, the quinolinate (Qa) phosphoribosyltransferase activity of NadC from S. pyogenes allows the organism to sustain growth when Qa is present as a sole pyridine precursor. Finally, the redundancy of functional upstream salvage pathways in GAS species narrows the choice of potential drug targets to the two indispensable downstream enzymes of NAD synthesis, nicotinate adenylyltransferase (NadD family) and NAD synthetase (NadE family). Biochemical characterization of NadD confirmed its functional role in S. pyogenes, and its potential as an antibacterial target was supported by inhibition studies with previously identified class I inhibitors of the NadD enzyme family. One of these inhibitors efficiently inhibited S. pyogenes NadD (sp.NadD) in vitro (50% inhibitory concentration [IC(50)], 15 μM), exhibiting a noncompetitive mechanism with a K(i) of 8 μM.

Glutamine versus ammonia utilization in the NAD synthetase family

NAD is a ubiquitous and essential metabolic redox cofactor which also functions as a substrate in certain regulatory pathways. The last step of NAD synthesis is the ATP-dependent amidation of deamido-NAD by NAD synthetase (NADS). Members of the NADS family are present in nearly all species across the three kingdoms of Life. In eukaryotic NADS, the core synthetase domain is fused with a nitrilase-like glutaminase domain supplying ammonia for the reaction. This two-domain NADS arrangement enabling the utilization of glutamine as nitrogen donor is also present in various bacterial lineages. However, many other bacterial members of NADS family do not contain a glutaminase domain, and they can utilize only ammonia (but not glutamine) in vitro. A single-domain NADS is also characteristic for nearly all Archaea, and its dependence on ammonia was demonstrated here for the representative enzyme from Methanocaldococcus jannaschi. However, a question about the actual in vivo nitrogen donor for single-domain members of the NADS family remained open: Is it glutamine hydrolyzed by a committed (but yet unknown) glutaminase subunit, as in most ATP-dependent amidotransferases, or free ammonia as in glutamine synthetase? Here we addressed this dilemma by combining evolutionary analysis of the NADS family with experimental characterization of two representative bacterial systems: a two-subunit NADS from Thermus thermophilus and a single-domain NADS from Salmonella typhimurium providing evidence that ammonia (and not glutamine) is the physiological substrate of a typical single-domain NADS. The latter represents the most likely ancestral form of NADS. The ability to utilize glutamine appears to have evolved via recruitment of a glutaminase subunit followed by domain fusion in an early branch of Bacteria. Further evolution of the NADS family included lineage-specific loss of one of the two alternative forms and horizontal gene transfer events. Lastly, we identified NADS structural elements associated with glutamine-utilizing capabilities.

Combining cheminformatics methods and pathway analysis to identify molecules with whole-cell activity against Mycobacterium tuberculosis

PURPOSE

New strategies for developing inhibitors of Mycobacterium tuberculosis (Mtb) are required in order to identify the next generation of tuberculosis (TB) drugs. Our approach leverages the integration of intensive data mining and curation and computational approaches, including cheminformatics combined with bioinformatics, to suggest biological targets and their small molecule modulators.

METHODS

We now describe an approach that uses the TBCyc pathway and genome database, the Collaborative Drug Discovery database of molecules with activity against Mtb and their associated targets, a 3D pharmacophore approach and Bayesian models of TB activity in order to select pathways and metabolites and ultimately prioritize molecules that may be acting as substrate mimics and exhibit activity against TB.

RESULTS

In this study we combined the TB cheminformatics and pathways databases that enabled us to computationally search >80,000 vendor available molecules and ultimately test 23 compounds in vitro that resulted in two compounds (N-(2-furylmethyl)-N'-[(5-nitro-3-thienyl)carbonyl]thiourea and N-[(5-nitro-3-thienyl)carbonyl]-N'-(2-thienylmethyl)thiourea) proposed as mimics of D-fructose 1,6 bisphosphate, (MIC of 20 and 40 μg/ml, respectively).

CONCLUSION

This is a simple yet novel approach that has the potential to identify inhibitors of bacterial growth as illustrated by compounds identified in this study that have activity against Mtb.

Annotator: postprocessing software for generating function-based signatures from quantitative mass spectrometry

Mass spectrometry is used to investigate global changes in protein abundance in cell lysates. Increasingly powerful methods of data collection have emerged over the past decade, but this has left researchers with the task of sifting through mountains of data for biologically significant results. Often, the end result is a list of proteins with no obvious quantitative relationships to define the larger context of changes in cell behavior. Researchers are often forced to perform a manual analysis from this list or to fall back on a range of disparate tools, which can hinder the communication of results and their reproducibility. To address these methodological problems, we developed Annotator, an application that filters validated mass spectrometry data and applies a battery of standardized heuristic and statistical tests to determine significance. To address systems-level interpretations, we incorporated UniProt and Gene Ontology keywords as statistical units of analysis, yielding quantitative information about changes in abundance for an entire functional category. This provides a consistent and quantitative method for formulating conclusions about cellular behavior, independent of network models or standard enrichment analyses. Annotator allows for "bottom-up" annotations that are based on experimental data and not inferred by comparison to external or hypothetical models. Annotator was developed as an independent postprocessing platform that runs on all common operating systems, thereby providing a useful tool for establishing the inherently dynamic nature of functional annotations, which depend on results from ongoing proteomic experiments. Annotator is available for download at http://people.cs.uchicago.edu/∼tyler/annotator/annotator_desktop_0.1.tar.gz .

Discovery of a gene involved in a third bacterial protoporphyrinogen oxidase activity through comparative genomic analysis and functional complementation

Tetrapyrroles are ubiquitous molecules in nearly all living organisms. Heme, an iron-containing tetrapyrrole, is widely distributed in nature, including most characterized aerobic and facultative bacteria. A large majority of bacteria that contain heme possess the ability to synthesize it. Despite this capability and the fact that the biosynthetic pathway has been well studied, enzymes catalyzing at least three steps have remained "missing" in many bacteria. In the current work, we have employed comparative genomics via the SEED genomic platform, coupled with experimental verification utilizing Acinetobacter baylyi ADP1, to identify one of the missing enzymes, a new protoporphyrinogen oxidase, the penultimate enzyme in heme biosynthesis. COG1981 was identified by genomic analysis as a candidate protein family for the missing enzyme in bacteria that lacked HemG or HemY, two known protoporphyrinogen oxidases. The predicted amino acid sequence of COG1981 is unlike those of the known enzymes HemG and HemY, but in some genomes, the gene encoding it is found neighboring other heme biosynthetic genes. When the COG1981 gene was deleted from the genome of A. baylyi, a bacterium that lacks both hemG and hemY, the organism became auxotrophic for heme. Cultures accumulated porphyrin intermediates, and crude cell extracts lacked protoporphyrinogen oxidase activity. The heme auxotrophy was rescued by the presence of a plasmid-borne protoporphyrinogen oxidase gene from a number of different organisms, such as hemG from Escherichia coli, hemY from Myxococcus xanthus, or the human gene for protoporphyrinogen oxidase.

Towards genome-scale signalling network reconstructions

Biological signalling networks allow living organisms to issue an integrated response to current conditions and make limited predictions about future environmental changes. Small-scale dynamic models of signalling cascades, including mitogen-activated protein kinase cascades, have been developed to generate hypotheses about signal transduction. Owing to technical limitations, these models and the hypotheses they generate have focused on a limited subset of signalling molecules. Now that we can simultaneously measure a substantial portion of the molecular components of a cell, we can begin to develop and test systems-level models of cellular signalling and regulatory processes, therefore gaining insights into the 'thought' processes of a cell.

Genomics-driven reconstruction of acinetobacter NAD metabolism: insights for antibacterial target selection

Enzymes involved in the last steps of NAD biogenesis, nicotinate mononucleotide adenylyltransferase (NadD) and NAD synthetase (NadE), are conserved and essential in most bacterial species and are established targets for antibacterial drug development. Our genomics-based reconstruction of NAD metabolism in the emerging pathogen Acinetobacter baumannii revealed unique features suggesting an alternative targeting strategy. Indeed, genomes of all analyzed Acinetobacter species do not encode NadD, which is functionally replaced by its distant homolog NadM. We combined bioinformatics with genetic and biochemical techniques to elucidate this and other important features of Acinetobacter NAD metabolism using a model (nonpathogenic) strain Acinetobacter baylyi sp. ADP1. Thus, a comparative kinetic characterization of PncA, PncB, and NadV enzymes allowed us to suggest distinct physiological roles for the two alternative, deamidating and nondeamidating, routes of nicotinamide salvage/recycling. The role of the NiaP transporter in both nicotinate and nicotinamide salvage was confirmed. The nondeamidating route was shown to be transcriptionally regulated by an ADP-ribose-responsive repressor NrtR. The NadM enzyme was shown to possess dual substrate specificity toward both nicotinate and nicotinamide mononucleotide substrates, which is consistent with its essential role in all three routes of NAD biogenesis, de novo synthesis as well as the two salvage pathways. The experimentally confirmed unconditional essentiality of nadM provided support for the choice of the respective enzyme as a drug target. In contrast, nadE, encoding a glutamine-dependent NAD synthetase, proved to be dispensable when the nondeamidating salvage pathway functioned as the only route of NAD biogenesis.

Complexes of bacterial nicotinate mononucleotide adenylyltransferase with inhibitors: implication for structure-based drug design and improvement

Bacterial nicotinate mononucleotide adenylyltransferase encoded by the essential gene nadD plays a central role in the synthesis of the redox cofactor NAD(+). The NadD enzyme is conserved in the majority of bacterial species and has been recognized as a novel target for developing new and potentially broad-spectrum antibacterial therapeutics. Here we report the crystal structures of Bacillus anthracis NadD in complex with three NadD inhibitors, including two analogues synthesized in the present study. These structures revealed a common binding site shared by different classes of NadD inhibitors and explored the chemical environment surrounding this site. The structural data obtained here also showed that the subtle changes in ligand structure can lead to significant changes in the binding mode, information that will be useful for future structure-based optimization and design of high affinity inhibitors.

'Unknown' proteins and 'orphan' enzymes: the missing half of the engineering parts list--and how to find it

Like other forms of engineering, metabolic engineering requires knowledge of the components (the 'parts list') of the target system. Lack of such knowledge impairs both rational engineering design and diagnosis of the reasons for failures; it also poses problems for the related field of metabolic reconstruction, which uses a cell's parts list to recreate its metabolic activities in silico. Despite spectacular progress in genome sequencing, the parts lists for most organisms that we seek to manipulate remain highly incomplete, due to the dual problem of 'unknown' proteins and 'orphan' enzymes. The former are all the proteins deduced from genome sequence that have no known function, and the latter are all the enzymes described in the literature (and often catalogued in the EC database) for which no corresponding gene has been reported. Unknown proteins constitute up to about half of the proteins in prokaryotic genomes, and much more than this in higher plants and animals. Orphan enzymes make up more than a third of the EC database. Attacking the 'missing parts list' problem is accordingly one of the great challenges for post-genomic biology, and a tremendous opportunity to discover new facets of life's machinery. Success will require a co-ordinated community-wide attack, sustained over years. In this attack, comparative genomics is probably the single most effective strategy, for it can reliably predict functions for unknown proteins and genes for orphan enzymes. Furthermore, it is cost-efficient and increasingly straightforward to deploy owing to a proliferation of databases and associated tools.

Targeting NAD biosynthesis in bacterial pathogens: Structure-based development of inhibitors of nicotinate mononucleotide adenylyltransferase NadD

The emergence of multidrug-resistant pathogens necessitates the search for new antibiotics acting on previously unexplored targets. Nicotinate mononucleotide adenylyltransferase of the NadD family, an essential enzyme of NAD biosynthesis in most bacteria, was selected as a target for structure-based inhibitor development. Using iterative in silico and in vitro screens, we identified small molecule compounds that efficiently inhibited target enzymes from Escherichia coli (ecNadD) and Bacillus anthracis (baNadD) but had no effect on functionally equivalent human enzymes. On-target antibacterial activity was demonstrated for some of the selected inhibitors. A 3D structure of baNadD was solved in complex with one of these inhibitors (3_02), providing mechanistic insights and guidelines for further improvement. Most importantly, the results of this study help validate NadD as a target for the development of antibacterial agents with potential broad-spectrum activity.

Controlled vocabularies for microbial virulence factors

Knowledge about pathogenesis is increasing dramatically, and most of this information is stored in the scientific literature or in sequence databases. This information can be made more accessible by the use of ontologies or controlled vocabularies. Recently, several ontologies, controlled vocabularies and databases have been developed or adapted for virulence factors and their roles in pathogenesis. Here, we discuss these systems, how they are being used in research and the challenges that remain for developing and applying ontologies for virulence factors.

Nicotinamide mononucleotide synthetase is the key enzyme for an alternative route of NAD biosynthesis in Francisella tularensis

Enzymes involved in the last 2 steps of nicotinamide adenine dinucleotide (NAD) cofactor biosynthesis, which catalyze the adenylylation of the nicotinic acid mononucleotide (NaMN) precursor to nicotinic acid dinucleotide (NaAD) followed by its amidation to NAD, constitute promising drug targets for the development of new antibiotics. These enzymes, NaMN adenylyltransferase (gene nadD) and NAD synthetase (gene nadE), respectively, are indispensable and conserved in nearly all bacterial pathogens. However, a comparative genome analysis of Francisella tularensis allowed us to predict the existence of an alternative route of NAD synthesis in this category A priority pathogen, the causative agent of tularaemia. In this route, the amidation of NaMN to nicotinamide mononucleotide (NMN) occurs before the adenylylation reaction, which converts this alternative intermediate to the NAD cofactor. The first step is catalyzed by NMN synthetase, which was identified and characterized in this study. A crystal structure of this enzyme, a divergent member of the NadE family, was solved at 1.9-A resolution in complex with reaction products, providing a rationale for its unusual substrate preference for NaMN over NaAD. The second step is performed by NMN adenylyltransferase of the NadM family. Here, we report validation of the predicted route (NaMN --> NMN --> NAD) in F. tularensis including mathematical modeling, in vitro reconstitution, and in vivo metabolite analysis in comparison with a canonical route (NaMN --> NaAD --> NAD) of NAD biosynthesis as represented by another deadly bacterial pathogen, Bacillus anthracis.

Methods of information geometry in computational system biology (consistency between chemical and biological evolution)

Interest in simulation of large-scale metabolic networks, species development, and genesis of various diseases requires new simulation techniques to accommodate the high complexity of realistic biological networks. Information geometry and topological formalisms are proposed to analyze information processes. We analyze the complexity of large-scale biological networks as well as transition of the system functionality due to modification in the system architecture, system environment, and system components. The dynamic core model is developed. The term dynamic core is used to define a set of causally related network functions. Delocalization of dynamic core model provides a mathematical formalism to analyze migration of specific functions in biosystems which undergo structure transition induced by the environment. The term delocalization is used to describe these processes of migration. We constructed a holographic model with self-poetic dynamic cores which preserves functional properties under those transitions. Topological constraints such as Ricci flow and Pfaff dimension were found for statistical manifolds which represent biological networks. These constraints can provide insight on processes of degeneration and recovery which take place in large-scale networks. We would like to suggest that therapies which are able to effectively implement estimated constraints, will successfully adjust biological systems and recover altered functionality. Also, we mathematically formulate the hypothesis that there is a direct consistency between biological and chemical evolution. Any set of causal relations within a biological network has its dual reimplementation in the chemistry of the system environment.

Bifunctional NMN adenylyltransferase/ADP-ribose pyrophosphatase: structure and function in bacterial NAD metabolism

Bacterial NadM-Nudix is a bifunctional enzyme containing a nicotinamide mononucleotide (NMN) adenylyltransferase and an ADP-ribose (ADPR) pyrophosphatase domain. While most members of this enzyme family, such as that from a model cyanobacterium Synechocystis sp., are involved primarily in nicotinamide adenine dinucleotide (NAD) salvage/recycling pathways, its close homolog in a category-A biodefense pathogen, Francisella tularensis, likely plays a central role in a recently discovered novel pathway of NAD de novo synthesis. The crystal structures of NadM-Nudix from both species, including their complexes with various ligands and catalytic metal ions, revealed detailed configurations of the substrate binding and catalytic sites in both domains. The structure of the N-terminal NadM domain may be exploited for designing new antitularemia therapeutics. The ADPR binding site in the C-terminal Nudix domain is substantially different from that of Escherichia coli ADPR pyrophosphatase, and is more similar to human NUDT9. The latter observation provided new insights into the ligand binding mode of ADPR-gated Ca2+ channel TRPM2.

Transcriptional regulation of NAD metabolism in bacteria: genomic reconstruction of NiaR (YrxA) regulon

A comparative genomic approach was used to reconstruct transcriptional regulation of NAD biosynthesis in bacteria containing orthologs of Bacillus subtilis gene yrxA, a previously identified niacin-responsive repressor of NAD de novo synthesis. Members of YrxA family (re-named here NiaR) are broadly conserved in the Bacillus/Clostridium group and in the deeply branching Fusobacteria and Thermotogales lineages. We analyzed upstream regions of genes associated with NAD biosynthesis to identify candidate NiaR-binding DNA motifs and assess the NiaR regulon content in these species. Representatives of the two distinct types of candidate NiaR-binding sites, characteristic of the Firmicutes and Thermotogales, were verified by an electrophoretic mobility shift assay. In addition to transcriptional control of the nadABC genes, the NiaR regulon in some species extends to niacin salvage (the pncAB genes) and includes uncharacterized membrane proteins possibly involved in niacin transport. The involvement in niacin uptake proposed for one of these proteins (re-named NiaP), encoded by the B. subtilis gene yceI, was experimentally verified. In addition to bacteria, members of the NiaP family are conserved in multicellular eukaryotes, including human, pointing to possible NaiP involvement in niacin utilization in these organisms. Overall, the analysis of the NiaR and NrtR regulons (described in the accompanying paper) revealed mechanisms of transcriptional regulation of NAD metabolism in nearly a hundred diverse bacteria.

Transcriptional regulation of NAD metabolism in bacteria: NrtR family of Nudix-related regulators

A novel family of transcription factors responsible for regulation of various aspects of NAD synthesis in a broad range of bacteria was identified by comparative genomics approach. Regulators of this family (here termed NrtR for Nudix-related transcriptional regulators), currently annotated as ADP-ribose pyrophosphatases from the Nudix family, are composed of an N-terminal Nudix-like effector domain and a C-terminal DNA-binding HTH-like domain. NrtR regulons were reconstructed in diverse bacterial genomes by identification and comparative analysis of NrtR-binding sites upstream of genes involved in NAD biosynthetic pathways. The candidate NrtR-binding DNA motifs showed significant variability between microbial lineages, although the common consensus sequence could be traced for most of them. Bioinformatics predictions were experimentally validated by gel mobility shift assays for two NrtR family representatives. ADP-ribose, the product of glycohydrolytic cleavage of NAD, was found to suppress the in vitro binding of NrtR proteins to their DNA target sites. In addition to a major role in the direct regulation of NAD homeostasis, some members of NrtR family appear to have been recruited for the regulation of other metabolic pathways, including sugar pentoses utilization and biogenesis of phosphoribosyl pyrophosphate. This work and the accompanying study of NiaR regulon demonstrate significant variability of regulatory strategies for control of NAD metabolic pathway in bacteria.