Identification of potential drug targets in Salmonella Typhimurium using metabolic modelling and experimental validation

almonella enterica serovar Typhimurium is an established model organism for Gram-negative, intracellular pathogens. Owing to the rapid spread of resistance to antibiotics among this group of pathogens, new approaches to identify suitable target proteins are required. Based on the genome sequence of Salmonella Typhimurium and associated databases, a genome-scale metabolic model was constructed. Output was based on an experimental determination of the biomass of Salmonella when growing in glucose minimal medium. Linear programming was used to simulate variations in energy demand, while growing in glucose minimal medium. By grouping reactions with similar flux responses, a sub-network of 34 reactions responding to this variation was identified (the catabolic core). This network was used to identify sets of one and two reactions, that when removed from the genome-scale model interfered with energy and biomass generation. 11 such sets were found to be essential for the production of biomass precursors. Experimental investigation of 7 of these showed that knock-outs of the associated genes resulted in attenuated growth for 4 pairs of reactions, while 3 single reactions were shown to be essential for growth.

The pgip family in soybean and three other legume species: evidence for a birth-and-death model of evolution

Background Polygalacturonase-inhibiting proteins (PGIPs) are leucine-rich repeat (LRR) plant cell wall glycoproteins involved in plant immunity. They are typically encoded by gene families with a small number of gene copies whose evolutionary origin has been poorly investigated. Here we report the complete characterization of the full complement of the pgip family in soybean (Glycine max [L.] Merr.) and the characterization of the genomic region surrounding the pgip family in four legume species. Results BAC clone and genome sequence analyses showed that the soybean genome contains two pgip loci. Each locus is composed of three clustered genes that are induced following infection with the fungal pathogen Sclerotinia sclerotiorum (Lib.) de Bary, and remnant sequences of pgip genes. The analyzed homeologous soybean genomic regions (about 126 Kb) that include the pgip loci are strongly conserved and this conservation extends also to the genomes of the legume species Phaseolus vulgaris L., Medicago truncatula Gaertn. and Cicer arietinum L., each containing a single pgip locus. Maximum likelihood-based gene trees suggest that the genes within the pgip clusters have independently undergone tandem duplication in each species. Conclusions The paleopolyploid soybean genome contains two pgip loci comprised in large and highly conserved duplicated regions, which are also conserved in bean, M. truncatula and C. arietinum. The genomic features of these legume pgip families suggest that the forces driving the evolution of pgip genes follow the birth-and-death model, similar to that proposed for the evolution of resistance (R) genes of NBS-LRR-type.

Temperature and oxygen dependent metabolite utilization by Salmonella enterica Serovars Derby and Mbandaka

Salmonella enterica is a zoonotic pathogen of clinical and veterinary significance, with over 2500 serovars. In previous work we compared two serovars displaying host associations inferred from isolation statistics. Here, to validate genome sequence data and to expand on the role of environmental metabolite constitution in host range determination we use a phenotypic microarray approach to assess the ability of these serovars to metabolise ~500 substrates at 25°C with oxygen (aerobic conditions) to represent the ex vivo environment and at 37°C with and without oxygen (aerobic/anaerobic conditions) to represent the in vivo environment. A total of 26 substrates elicited a significant difference in the rate of metabolism of which only one, D-galactonic acid-g-lactone, could be explained by the presence (S. Mbandaka) or the absence (S. Derby) of metabolic genes. We find that S. Mbandaka respires more efficiently at ambient temperatures and under aerobic conditions on 18 substrates including: glucosominic acid, saccharic acid, trehalose, fumaric acid, maltotriose, N-acetyl-D-glucosamine, N-acetyl-beta-D-mannosamine, fucose, L-serine and dihydroxy-acetone; whereas S. Derby is more metabolically competent anaerobically at 37°C for dipeptides, glutamine-glutamine, alanine-lysine, asparagine-glutamine and nitrogen sources glycine and nitrite. We conclude that the specific phenotype cannot be reliably predicted from the presence of metabolic genes directly relating to the metabolic pathways under study.

Structural and phylogenetic relationship of ORF 31 from the Anticarsia gemmatalis MNPV to poly (ADP-ribose) polymerases (PARP)

ORF 31 is a unique baculovirus gene in the genome of Anticarsia gemmatalis multiple nucleopolyhedrovirus isolate 2D (AgMNPV-2D). It encodes a putative polypeptide of 369 aa homologous to poly (ADP-ribose) polymerase (PARP) found in the genomes of several organisms. Moreover, we found a phylogenetic association with Group I PARP proteins and a 3D homology model of its conserved PARP C-terminal catalytic domain indicating that had almost an exact spatial superimposition of < 1 angstrom with other PARP available structures. The 5` end of ORF 31 mRNA was located at the first nucleotide of a CATT motif at position -27. Using real-time PCR we detected transcripts at 3 h post-infection (p.i.) increasing until 24 h p.i., which coincides with the onset of DNA replication, suggestive of a possible role in DNA metabolism.

Evolutionary placement of Xanthomonadales based on conserved protein signature sequences

Xanthomonadales comprises one of the largest phytopathogenic bacterial groups, and is currently classified within the gamma-proteobacteria. However, the phylogenetic placement of this group is not clearly resolved, and the results of different studies contradict one another. In this work, the evolutionary position of Xanthomonadales was determined by analyzing the presence of shared insertions and deletions (INDELs) in highly conserved proteins. Several distinctive insertions found in most of the members of the gamma-proteobacteria are absent in Xanthomonadales and groups such as Legionelalles, Chromatiales, Methylococcales, Thiotrichales and Cardiobacteriales. These INDELs were most likely introduced after the branching of Xanthomonadales from most of the gamma-proteobacteria and provide evidence for the phylogenetic placement of the early gamma-proteobacteria. Moreover, other proteins contain insertions exclusive to the Xanthomonadales order, confirming that this is a monophyletic group and provide important specific genetic markers. Thus, the data presented clearly support the Xanthomonadales group as an independent subdivision, and constitute one of the deepest branching lineage within the gamma-proteobacteria clade. (C) 2009 Elsevier Inc. All rights reserved.

Systematic structural studies of iron superoxide dismutases from human parasites and a statistical coupling analysis of metal binding specificity

Superoxide dismutases (SODs) are a crucial class of enzymes in the combat against intracellular free radical damage. They eliminate superoxide radicals by converting them into hydrogen peroxide and oxygen. In spite of their very different life cycles and infection strategies, the human parasites Plasmodium falciparum, Trypanosoma cruzi and Trypanosoma brucei are known to be sensitive to oxidative stress. Thus the parasite Fe-SODs have become attractive targets for novel drug development. Here we report the crystal structures of FeSODs from the trypanosomes T. brucei at 2.0 angstrom and T. cruzi at 1.9 angstrom resolution, and that from P. falciparum at a higher resolution (2.0 angstrom) to that previously reported. The homodimeric enzymes are compared to the related human MnSOD with particular attention to structural aspects which are relevant for drug design. Although the structures possess a very similar overall fold, differences between the enzymes at the entrance to the channel which leads to the active site could be identified. These lead to a slightly broader and more positively charged cavity in the parasite enzymes. Furthermore, a statistical coupling analysis (SCA) for the whole Fe/MnSOD family reveals different patterns of residue coupling for Mn and Fe SODs, as well as for the dimeric and tetrameric states. In both cases, the statistically coupled residues lie adjacent to the conserved core surrounding the metal center and may be expected to be responsible for its fine tuning, leading to metal ion specificity.

Sugarcane Phosphoribosyl Pyrophosphate Synthetase: Molecular Characterization of a Phosphate-independent PRS

Phosphoribosyl pyrophosphate synthetase (PRS-EC:2.7.6.1) is an important enzyme present in several metabolic pathways, thus forming a complex family of isoenzymes. However, plant PRS enzymes have not been extensively investigated. In this study, a sugarcane prs gene has been characterized from the Sugar Cane Expressed Sequence Tag Genome Project. This gene contains a 984-bp open reading frame encoding a 328-amino acid protein. The predicted amino acid sequence has 77% and 78% amino acid sequence identity to Arabidopsis thaliana and Spinacia oleracea PRS4, respectively. The assignment of sugarcane PRS as a phosphate-independent PRS isoenzyme (Class II PRS) is verified following enzyme assay and phylogenetic reconstruction of PRS homologues. To gain further insight into the structural framework of the phosphate independence of sugarcane PRS, a molecular model is described. This model reveals the formation of two conserved domains elucidating the structural features involved in sugarcane PRS phosphate independence. The recombinant PRS retains secondary structure elements and a quaternary arrangement consistent with known PRS homologues, based on circular dichroism measurements.

The Aedes aegypti larval transcriptome: a comparative perspective with emphasis on trypsins and the domain structure of peritrophins

The genome sequence of Aedes aegypti was recently reported. A significant amount of Expressed Sequence Tags (ESTs) were sequenced to aid in the gene prediction process. In the present work we describe an integrated analysis of the genomic and EST data, focusing on genes with preferential expression in larvae (LG), adults (AG) and in both stages (SG). A total of 913 genes (5.4% of the transcript complement) are LG, including ion transporters and cuticle proteins that are important for ion homeostasis and defense. From a starting set of 245 genes encoding the trypsin domain, we identified 66 putative LG, AG, and SG trypsins by manual curation. Phylogenetic analyses showed that AG trypsins are divergent from their larval counterparts (LG), grouping with blood-induced trypsins from Anopheles gambiae and Simulium vittatum. These results support the hypothesis that blood-feeding arose only once, in the ancestral Culicomorpha. Peritrophins are proteins that interlock chitin fibrils to form the peritrophic membrane (PM) that compartmentalizes the food in the midgut. These proteins are recognized by having chitin-binding domains with 6 conserved Cys and may also present mucin-like domains (regions expected to be highly O-glycosylated). PM may be formed by a ring of cells (type 2, seen in Ae. aegypti larvae and Drosophila melanogaster) or by most midgut cells (type 1, found in Ae. aegypti adult and Tribolium castaneum). LG and D. melanogaster peritrophins have more complex domain structures than AG and T. castaneum peritrophins. Furthermore, mucin-like domains of peritrophins from T. castaneum (feeding on rough food) are lengthier than those of adult Ae. aegypti (blood-feeding). This suggests, for the first time, that type 1 and type 2 PM may have variable molecular architectures determined by different peritrophins and/or ancillary proteins, which may be partly modulated by diet.

Highly polymorphic in silico-derived microsatellite loci in the potato-infecting fungal pathogen Rhizoctonia solani anastomosis group 3 from the Colombian Andes

Fourteen polymorphic microsatellite DNA markers derived from the draft genome sequence of Rhizoctonia solani anastomosis group 3 (AG-3), strain Rhs 1AP, were designed and characterized from the potato-infecting soil fungus R. solani AG-3. All loci were polymorphic in two field populations collected from Solanum tuberosum and S. phureja in the Colombian Andes. The total number of alleles per locus ranged from two to seven, while gene diversity (expected heterozygosity) varied from 0.11 to 0.81. Considering the variable levels of genetic diversity observed, these markers should be useful for population genetic analyses of this important dikaryotic fungal pathogen on a global scale.

Contributions and perspectives of chicken genomics in Brazil: from biological model to export commodity

Chicken is one of the most important sources of animal protein for human consumption, and breeding programmes have been responsible for constant improvements in production efficiency and product quality. Furthermore, chicken has largely contributed to fundamental discoveries in biology for the last 100 years. In this article we review recent developments in poultry genomics and their contribution to adding functional information to the already existing structural genomics, including the availability of the complete genome sequence, a comprehensive collection of mRNA sequences ( ESTs), microarray platforms, and their use to complement QTL mapping strategies in the identification of genes that underlie complex traits. Efforts of the Brazilian Poultry Genomics Programme in this area resulted in generation of a resource population, which was used for identification of Quantitative Trait Loci ( QTL) regions, generation of ESTs and candidate gene studies that contributed to furthering our understanding of the complex biological processes involved in growth and muscular development in chicken.

Structural bioinformatics study of EPSP synthase from Mycobacterium tuberculosis

The shikimate pathway is an attractive target for herbicides and antimicrobial agent development because it is essential in algae, higher plants, bacteria, and fungi, but absent from mammals. Homologues to enzymes in the shikimate pathway have been identified in the genome sequence of Mycobacterium tuberculosis. Among them, the EPSP synthase was proposed to be present by sequence homology. Accordingly, in order to pave the way for structural and functional efforts towards anti-mycobacterial agent development, here we describe the molecular modeling of 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase isolated from M. tuberculosis that should provide a structural framework on which the design of specific inhibitors may be based on. Significant differences in the relative orientation of the domains in the two models result in open and closed conformations. The possible relevance of this structural transition in the ligand biding is discussed. (C) 2003 Elsevier B.V. All rights reserved.

Functional shikimate dehydrogenase from Mycobacterium tuberculosis H37Rv: Purification and characterization

Tuberculosis (TB) poses a major worldwide public health problem. The increasing prevalence of TB, the emergence of multi-drug-resistant strains of Mycobacterium tuberculosis, the causative agent of TB, and the devastating effect of co-infection with HIV have highlighted the urgent need for the development of new antimycobacterial agents. Analysis of the complete genome sequence of M. tuberculosis shows the presence of genes involved in the aromatic amino acid biosynthetic pathway. Experimental evidence that this pathway is essential for M. tuberculosis has been reported. The genes and pathways that are essential for the growth of the microorganisms make them attractive drug targets since inhibiting their function may kill the bacilli. We have previously cloned and expressed in the soluble form the fourth shikimate pathway enzyme of the M. tuberculosis, the aroE-encoded shikimate dehydrogenase (mtSD). Here, we present the purification of active recombinant aroE-encoded M. tuberculosis shikimate dehydrogenase (mtSD) to homogeneity, N-terminal sequencing, mass spectrometry, assessment of the oligomeric state by gel filtration chromatography, determination of apparent steady-state kinetic parameters for both the forward and reverse directions, apparent equilibrium constant, thermal stability, and energy of activation for the enzyme-catalyzed chemical reaction. These results pave the way for structural and kinetic studies, which should aid in the rational design of mtSD inhibitors to be tested as antimycobacterial agents. (c) 2005 Elsevier B.V. All rights reserved.

Molecular models of NS3 protease variants of the Hepatitis C virus

Background: Hepatitis C virus (HCV) currently infects approximately three percent of the world population. In view of the lack of vaccines against HCV, there is an urgent need for an efficient treatment of the disease by an effective antiviral drug. Rational drug design has not been the primary way for discovering major therapeutics. Nevertheless, there are reports of success in the development of inhibitor using a structure-based approach. One of the possible targets for drug development against HCV is the NS3 protease variants. Based on the three-dimensional structure of these variants we expect to identify new NS3 protease inhibitors. In order to speed up the modeling process all NS3 protease variant models were generated in a Beowulf cluster. The potential of the structural bioinformatics for development of new antiviral drugs is discussed.Results: the atomic coordinates of crystallographic structure 1CU1 and 1DY9 were used as starting model for modeling of the NS3 protease variant structures. The NS3 protease variant structures are composed of six subdomains, which occur in sequence along the polypeptide chain. The protease domain exhibits the dual beta-barrel fold that is common among members of the chymotrypsin serine protease family. The helicase domain contains two structurally related beta-alpha-beta subdomains and a third subdomain of seven helices and three short beta strands. The latter domain is usually referred to as the helicase alpha-helical subdomain. The rmsd value of bond lengths and bond angles, the average G-factor and Verify 3D values are presented for NS3 protease variant structures.Conclusions: This project increases the certainty that homology modeling is an useful tool in structural biology and that it can be very valuable in annotating genome sequence information and contributing to structural and functional genomics from virus. The structural models will be used to guide future efforts in the structure-based drug design of a new generation of NS3 protease variants inhibitors. All models in the database are publicly accessible via our interactive website, providing us with large amount of structural models for use in protein-ligand docking analysis.

Preliminary radiation hybrid map for river buffalo chromosome 6 and comparison to bovine chromosome 3

We present the first radiation hybrid (RH) map of river buffalo (Bubalus bubalis) chromosome 6 (BBU6) developed with a recently constructed river buffalo whole-genome RH panel (BBURH5000). The preliminary map contains 33 cattle-derived markers, including 12 microsatellites, 19 coding genes and two ESTs, distributed across two linkage groups. Retention frequencies for markers ranged from 14.4% to 40.0%. Most of the marker orders within the linkage groups on BBU6 were consistent with the cattle genome sequence and RH maps. This preliminary RH map is the starting point for comparing gene order between river buffalo and cattle, presenting an opportunity for the examination of micro-rearrangements of these chromosomes. Also, resources for positional candidate cloning in river buffalo are enhanced.

The copper resistance operon copAB from Xanthomonas axonopodis pathovar citri: gene inactivation results in copper sensitivity

Xanthomonas axonopodis pv. citri (Xac) causes citrus canker and the completion of the Xac genome sequence has opened up the possibility of investigating basic cellular mechanisms at the genomic level. Copper compounds have been extensively used in agriculture to control plant diseases. The copA and copB genes, identified by annotation of the Xac genome, encode homologues of proteins involved in copper resistance. A gene expression assay by Northern blotting revealed that copA and copB are expressed as a unique transcript specifically induced by copper. Synthesis of the gene products was also induced by copper, reaching a maximum level at 4 h after addition of copper to the culture medium. CopA was a cytosolic protein and CopB was detected in the cytoplasmic membrane. The gene encoding CopA was disrupted by the insertion of a transposon, leading to mutant strains that were unable to grow in culture medium containing copper, even at the lowest CUSO4 concentration tested (0.25 mM), whereas the wild-type strain was able to grow in the presence of 1 mM copper. Cell suspensions of the wild-type and mutant strains in different copper concentrations were inoculated in lemon leaves to analyse their ability to induce citrus canker symptoms. Cells of mutant strains showed higher sensitivity than the wild-type strain in the presence of copper, i.e. they were not able to induce citrus canker symptoms at high copper concentrations and exhibited a more retarded growth in planta.