- Open Access
An integrated structural proteomics approach along the druggable genome of Corynebacterium pseudotuberculosis species for putative druggable targets
- Leandro G Radusky†1,
- Syed Shah Hassan†3,
- Esteban Lanzarotti1,
- Sandeep Tiwari3,
- Syed Babar Jamal3,
- Javed Ali4,
- Amjad Ali5,
- Rafaela Salgado Ferreira6,
- Debmalya Barh7,
- Artur Silva8,
- Adrián G Turjanski1, 2Email author and
- Vasco AC Azevedo3Email author
© Radusky et al.; licensee BioMed Central Ltd. 2015
- Published: 26 May 2015
The bacterium Corynebacterium pseudotuberculosis (Cp) causes caseous lymphadenitis (CLA), mastitis, ulcerative lymphangitis, and oedema in a number of hosts, comprising ruminants, thereby intimidating economic and dairy industries worldwide. So far there is no effective drug or vaccine available against Cp. Previously, a pan-genomic analysis was performed for both biovar equi and biovar ovis and a Pathogenicity Islands (PAIS) analysis within the strains highlighted a large set of proteins that could be relevant therapeutic targets for controlling the onset of CLA. In the present work, a structural druggability analysis pipeline was accomplished along 15 previously sequenced Cp strains from both biovar equi and biovar ovis.
Methods and results
We computed the whole modelome of a reference strain Cp1002 (NCBI Accession: NC_017300.1) and then the homology models of proteins, of 14 different Cp strains, with high identity (≥ 85%) to the reference strain were also done. Druggability score of all proteins pockets was calculated and only those targets that have a highly druggable (HD) pocket in all strains were kept, a set of 58 proteins. Finally, this information was merged with the previous PAIS analysis giving two possible highly relevant targets to conduct drug discovery projects. Also, off-targeting information against host organisms, including Homo sapiens and a further analysis for protein essentiality provided a final set of 31 druggable, essential and non-host homologous targets, tabulated in table S4, additional file 1. Out of 31 globally druggable targets, 9 targets have already been reported in other pathogenic microorganisms, 3 of them (3-isopropylmalate dehydratase small subunit, 50S ribosomal protein L30, Chromosomal replication initiator protein DnaA) in C. pseudotuberculosis.
Overall we provide valuable information of possible targets against C. pseudotuberculosis where some of these targets have already been reported in other microorganisms for drug discovery projects, also discarding targets that might be physiologically relevant but are not amenable for drug binding. We propose that the constructed in silico dataset might serve as a guidance for the scientific community to have a better understanding while selecting putative therapeutic protein candidates as druggable ones as effective measures against C. pseudotuberculosis.
- Corynebacterium pseudotuberculosis (Cp)
- Druggable genome
- Structural proteomics approach
- Putative globally/conserved druggable/bindable targets
- Caseous lymphadenitis
Efforts to find new bacterial drug and/or vaccine targets are becoming indispensable due to the antimicrobial resistance, rapid loss of effectiveness in antibiotic treatment and the quantitative emergence of multi-resistant microbial strains that pose a global challenge and threat. Corynebacterium pseudotuberculosis (Cp) is a pathogen of great veterinary and economic importance, since it affects a broad spectrum of animal livestock worldwide, mainly sheep and goats, as well as mammals in numerous Asiatic, Arabic and African countries, North and Latin America and Australia . C. pseudotuberculosis is a Gram-positive, facultative intracellular and pleomorphic organism; it possesses fimbriae but is non-motile in nature . The rpoB gene analysis for CMNR group of bacteria (Corynebacterium, Mycobacterium, Nocardia and Rhodococcus), which has a great medical, veterinary and biotechnological importance, has shown a close phylogenetic relationship . A number of pathogenic strains from a wide range of hosts have already been sequenced, demonstrating the importance of this microorganism . The pathogen infects goat and sheep populations (biovar ovis), causing caseous lymphadenitis (CLA), a chronic contagious disease with abscess formation in superficial lymph nodes and subcutaneous tissues. Biovar equi, on the other hand, infects lung, kidney, liver and spleen in higher mammals like cow, camel, buffalo etc. thereby, threatening the life of herd animals [2, 4]. There are few reports in humans of symptoms similar to lymphadenitis abscesses, caused by an occupational exposure to the infected animals [5–7]. Bearing in mind the medical importance of C. pseudotuberculosis due to a lack of efficient medicines, here, we have made an effort and applied a computational strategy to search for new therapeutic molecular targets from this bacterium.
Homology modelling is a widely used technique that has proved good results to expand structural space of pathogens [8, 9]. We have designed and implemented a protein structure prediction pipeline using homology modelling based on Martí-Renom methodologies . The pipeline was applied to the randomly selected Cp1002 genome as a reference strain; the genomes of other strains were modelled using a mutation procedure, the sequences that present homology in both, i.e., the core genome, were already modelled in the reference strain Cp1002.
The main purpose of this study is to offer information based on recently reported structure-based prediction of protein druggability that might be valuable for target selection in drug design projects. Druggability is a concept used to describe the ability of a given protein to bind a drug-like molecule, which in turn modulates its function in a desired way. Purely, from a structural point of view, it is the likelihood that small drug-like molecules bind a given target protein with high affinity (< 1µM), a concept also referred to "bind ability". Early attempts to determine the druggable genome of an organism were based on counting the number of targets belonging to known druggable domains that have yielded values in range of 10-14% for the human genome. According to our knowledge, a structural druggability assessment for the C. pseudotuberculosis multi-strain proteome was never performed before.
Recently, we developed a fast method for druggability prediction based on the open source pocket detection code "fpocket", which combines several physicochemical descriptors to estimate the pocket druggability and that can be used on a genomic scale . Druggable proteins should have a pocket with suitable features that enable binding of a drug-like compound. After computing proteins that remains druggable along the 15 Cp genomes, 58 target candidates were selected.
The Cp genome has been reported to include seven putative pathogenicity islands (PAIS) , which contain several classical virulence factors, including genes for fimbrial subunits, adhesion factors, iron uptake and secreted toxins. Additionally, all of the virulence factors in the islands have characteristics that indicate the phenomenon of horizontal gene transfer. The importance of our dataset is enhanced with the emerged information from the literature regarding the PAIS, pan-genome and also with the pan-modelome strategies for target selection .
Initial dataset construction
All ORFs or possible proteins for all the strains of Cp were obtained by downloading the information available at the NCBI database (ftp://ftp.ncbi.nih.gov/genomes/Bacteria). The randomly selected strain used as a reference genome for further calculations is Cp1002, according to recent work by Hassan et al., 2014 . The Cp1002 genome has 2097 reported ORFs.
Pfam domain assignment
All the ORFs in the reference proteome were analysed with the HMMer program and were later assigned the Pfam families or domains, leading to a total of 2455 domains assignments from Pfam-A entries and 509 ORFs with no domain assigned. However, as expected, more than one ORF can be assigned to the same domain. Thus, considering this information, a total of 1327 unique (i.e. different) domains were assigned to a whole Cp reference genome. On average, the Cp reference genome has 1.87 domains per ORF and 1.58 unique domains per ORF.
Generation of structural homology-based models
The strain Cp1002 was used a basis for the structural study. For each sequence in this strain, several models were built using the following procedures: First, a PSI-Blast  search was performed against UniRef50, using 3 iterations and an e.value threshold of 0.0001, in order to compute a checkpoint that will be used as a profile for the target sequence. Second, PSI-Blast search is restarted using the aforementioned checkpoint against a template library. The template library consisted of all sequences from every individual protein chain in the Protein Data Bank (PDB), grouped at 95% sequence identity threshold using CD-HIT . Afterward, the target-template alignments were computed in the last step, with an e-value of ≤ 0.00001 to build the model structures using the MODELLER software [17, 18]. For each target-template alignment, ten different target models were built, and their quality measures were assessed using the GA341 [19, 20] and QMEAN methods . Only the models with GA341 score ≥ 0.7 (score for the reliability of a model derived from statistical potentials) , were retained. A reliable model has a probability of correct fold larger than 95% and coverage of over 50%. Every sequence belonging to the other 14 strains were compared to strain Cp1002 using BLASTp. For each sequence that gives a hit with sequence identity above 85%, a mutation methodology was applied on each amino acid substitution using MODELLER between this sequence and the sequence belonging to Cp1002.
Structural assessment of druggability
Extending this information along the different strains, those structures labelled as HD in the reference strain were analysed in the other ones. A protein target, which remained druggable in all the strains, was classified as Globally Druggable (GD), also, a protein druggable in all the virulent strains is classified as Virulent Druggable (VD). The GD, obviously, remained a subset of VD.
Active site identification
In order to identify the active site pocket, our pipeline implements two different analyses that rely on; i) the information from the CSA, and ii) a Pfam position site, as the importance criteria.
The data from CSA was downloaded from http://www.ebi.ac.uk/thornton-srv/databases/CSA/ that consisted of a list of PDB_IDs linked to a number of residues, which comprised the corresponding protein active site. To map the active site residues to as many Cp protein domains as possible, each PDB_ID used as a template for the homology modelling in CSA was assigned to the modelled ORFs.
As an alternative approach to determine the relevance of a given pocket (or residues), we looked for residues of a given Pfam family/domain that were located in an important position and are well conserved. Important positions were defined as those positions in the corresponding HMMer model whose information content was larger than a defined importance cut-off value (icov). The nature of the conserved amino acids in the corresponding position was determined by comparing each residue type emission probability (ep) with icov. If the ratio between ep and icov was larger than a conserved type cut-off value (ctcov), the corresponding residue type was assumed as conserved. Optimal values of icov and ctcov  were 0.27 and 0.24, respectively. These values were calculated as the ones that allowed labelling all the important residues of the CSA database in their respective protein's domains. Briefly, this approach is a strategy to extend the definition of catalytic site capturing all the sites described in CSA. CSA is a curated database with limited amount of data. Our strategy gives us a clue about new candidates of being catalytic sites, with some false positives.
By using these analyses, for each Pfam domain, the pipeline provided a list of position-residue type relevant residues, which could thus be mapped to all Cp ORFs with assigned Pfam domain.
Pathogenicity islands and pan-modelomics information
Information collected from other software pipelines is used in the present work to enhance the target selection process.
All the ORFs that belonged to a pathogenicity island (PAIS) were properly labelled with this indicator. The PAIS were previously computed and identified using the PIPS software  that has predicted 16 pathogenicity islands in C. pseudotuberculosis. PIPS predicts pathogenicity islands by taking into account some important features, i.e., flanking tRNA, codon usage bias, GC content deviation; transposases, virulence factors and their absence in non-pathogenic organism of the same genus or related species. Pathogenicity islands are large regions that were acquired through horizontal gene transfer that represent the genome plasticity of a species and possess a high concentration of virulence factors. Virulence factors are proteins whose function is related to bacterial virulence and pathogenicity. They help the pathogenic bacterium in adhesion mechanisms, invasion, and surviving through colonization, and replication inside the host as well as in immune system evasion .
A novel integrative approach has been adapted in a recent work for the identification of new therapeutic targets in C. pseudotuberculosis  where a final set of 10 proteins has been selected that was essential for the bacterial survival. Here, too, a focus was made only on the previously selected homologous target proteins in the reference Cp1002 genome, whose detailed analysis with our pipeline is reported.
Bacterial essential and non-host homologous proteins analysis
The pool of global highly druggable (GD) 58 target proteins was subjected to NCBI-BLASTp (e-value lesser than 1e-07, bit score ≥100, identity ≥ 50%) against human and ovine proteomes to identify non-host homologs targets, and coverage of all the sequences of the proteins greater than 80% . The exhaustive list of GD and Virulent highly druggable (VD) ORFs is available in supplementary material table S1, additional file 1.
Furthermore, from the filtered list of 41 highly druggable non-host homologous target proteins, an approach based on a subtractive genomics was made for the conserved GD targets that were essential to bacteria . Briefly, the set of proteins from C. pseudotuberculosis were submitted to the Database of Essential Genes (DEG, which contains experimentally validated essential genes from 20 bacteria) for homology analyses . The BLASTp cut-off values used were: e-value = 1e-05, bit score ≥100, identity ≥ 35% . The final list of targets based on previously described criteria contained 31 essential and non-host homologous target proteins. The list of putative targets was subjected to biochemical pathway analysis to KEGG (Kyoto Encyclopaedia of Genes and Genomes) , virulence using PIPS (Pathogenicity Island Prediction Software) , functionality using UniProt (Universal Protein Resource) , and cellular localization using CELLO (subCELlular LOcalization predictor) .
Models in the reference Cp1002 strain
The following subsections describe the database data summary after running the HMMer software and the modelling pipeline.
A total of 2598 models were built from the ORFs of the reference genome and then a part of the ORFs were assigned to the Pfam families. Here, 1206 unique ORFs were involved. The models from the Pfam family were 1546 in total, where only 879 unique Pfam families were used. The other 1051 models correspond to full ORF models (table 1).
Models in non-reference strains
Summary of modelling pipeline over non-reference strains.
Homolog ORFs to Cp1002
Druggability summary of 15 C. pseudotuberculosis strains
Summary of druggability for all C. pseudotuberculosis strains, Global highly druggable (GD) ORFs and Virulent highly druggable (VD) ORFs.
ND or PD
Global highly druggable (GD) ORFs
Virulent highly druggable (VD) ORFs
Bacterial essential and non-host homologous proteins analyses
An exhaustive literature review of candidate druggable targets
As aforementioned, the list of 58 most druggable proteins (table S1, additional file 1) was compared to the corresponding host proteomes, leading to the identification of 41 non-host homologous proteins (table S2, additional file 1) and 17 host homologous proteins (table S3, additional file 1). A final list of 31 essential and non-host homologous targets from C. pseudotuberculosis is given in table S4, additional file 1 after computing the 41 non-host homologous proteins to the database of essential genes (DEG), using the earlier mentioned default parameters.
Essential and non-host homologous bacterial protein already reported as drug target in other pathogenic microorganisms.
Reference Protein Loci
1,4-alpha-glucan branching enzyme GlgB
3-isopropylmalate dehydratase small subunit
50S ribosomal protein L30
Bifunctional protein folC
Riboflavin biosynthesis protein ribF
Chromosomal replication initiator protein DnaA
In the present work we have attempted to show a comprehensive study of the druggability scores along the known completely sequenced strains of C. pseudotuberculosis species to complement further the research work performed by our colleagues and collaborators. After our pipeline was executed, a list of cross-strain highly druggable 58 ORFs was obtained. These ORFs had the information about their membership in the set of virulent strains. We too have provided the information if the most druggable pockets of the ORF have highly conserved residues in the Pfam families they belong to, and if a catalytic site is reported in the template structure or structures that were used to build the homology model for each ORF. We expect that the constructed dataset might serve as a guide for the scientific community to have a better understanding while selecting protein candidates as therapeutic and druggable ones. All the data is available via web at http://www.inquimae.fcen.uba.ar/turjanski_adrian.
The information obtained here was also used to exhaustively analyse the target proteins actually considered as druggable targets in previous works. Our analyses have proposed a set of putative druggable targets in the veterinary pathogen C. pseudotuberculosis on one hand, while on the other hand it has also demonstrated the efficiency and the high-throughput nature of our pipeline used in this study. All this work is expandable and could be applied to the emergence of new strains of the same organism species as well as to the new organisms with the same characteristics.
This work resulted from a mutual collaboration among various, especially the two leading research and academic institutions of Argentine and Brazil, the Structural Bioinformatics Group, Institute of Physical Chemistry of Materials, Environment and Energy, University of Buenos Aires, Argentine and the Laboratory of Cellular and Molecular Genetics, PG Program in Bioinformatics, Institute of Biological Sciences, Federal University of Minas Gerais, Belo Horizonte, Minas Gerais, Brazil. We acknowledge an immense support from all the co-authors who have worked hard in accomplishing this task. We are grateful to our Brazilian collaborators, especially Prof. Dr. Vasco A.C Azevedo, Dr. Syed Shah Hassan (postdoctoral fellow) and other team members for their collaboration, timely consultation in theoretical and biological knowledge regarding the genetics background of the microorganism, C. pseudotuberculosis.
Publication costs for this article were funded by grant from the Argentinian National Agency PICT-2010-2805 and VA from the Programa Nacional de Pós-doutorado da CAPES-PNPD/2011: Concessão Institucional, Brazil.
This article has been published as part of BMC Genomics Volume 16 Supplement 5, 2015: Proceedings of the 10th International Conference of the Brazilian Association for Bioinformatics and Computational Biology (X-Meeting 2014). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S5.
- Hassan SS, Schneider MP, Ramos RT, Carneiro AR, Ranieri A, Guimaraes LC, et al: Whole-genome sequence of Corynebacterium pseudotuberculosis strain Cp162, isolated from camel. J Bacteriol. 2012, 194 (20): 5718-5719. 10.1128/JB.01373-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Dorella FA, Pacheco LG, Oliveira SC, Miyoshi A, Azevedo V: Corynebacterium pseudotuberculosis: microbiology, biochemical properties, pathogenesis and molecular studies of virulence. Vet Res. 2006, 37 (2): 201-218. 10.1051/vetres:2005056.View ArticlePubMedGoogle Scholar
- Soares SC, Trost E, Ramos RT, Carneiro AR, Santos AR, Pinto AC, et al: Genome sequence of Corynebacterium pseudotuberculosis biovar equi strain 258 and prediction of antigenic targets to improve biotechnological vaccine production. J Biotechnol. 2013, 167 (2): 135-141. 10.1016/j.jbiotec.2012.11.003.View ArticlePubMedGoogle Scholar
- Khamis A, Raoult D, La Scola B: Comparison between rpoB and 16S rRNA gene sequencing for molecular identification of 168 clinical isolates of Corynebacterium. J Clin Microbiol. 2005, 43 (4): 1934-1936. 10.1128/JCM.43.4.1934-1936.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Luis MA, Lunetta AC: [Alcohol and drugs: preliminary survey of Brazilian nursing research]. Rev Lat Am Enfermagem. 2005, 1219-1230. 13 Spec NoGoogle Scholar
- Peel MM, Palmer GG, Stacpoole AM, Kerr TG: Human lymphadenitis due to Corynebacterium pseudotuberculosis: report of ten cases from Australia and review. Clin Infect Dis. 1997, 24 (2): 185-191. 10.1093/clinids/24.2.185.View ArticlePubMedGoogle Scholar
- Williamson LH: Caseous lymphadenitis in small ruminants. Vet Clin North Am Food Anim Pract. 2001, 17 (2): 359-371. viiPubMedGoogle Scholar
- Kinnings SL, Xie L, Fung KH, Jackson RM, Xie L, Bourne PE: The Mycobacterium tuberculosis drugome and its polypharmacological implications. PLoS Comput Biol. 2010, 6 (11): e1000976-10.1371/journal.pcbi.1000976.PubMed CentralView ArticlePubMedGoogle Scholar
- Anand P, Sankaran S, Mukherjee S, Yeturu K, Laskowski R, Bhardwaj A, et al: Structural annotation of Mycobacterium tuberculosis proteome. PLoS One. 2011, 6 (10): e27044-10.1371/journal.pone.0027044.PubMed CentralView ArticlePubMedGoogle Scholar
- Marti-Renom MA, Stuart AC, Fiser A, Sanchez R, Melo F, Sali A: Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct. 2000, 29: 291-325. 10.1146/annurev.biophys.29.1.291.View ArticlePubMedGoogle Scholar
- Radusky L, Defelipe LA, Lanzarotti E, Luque J, Barril X, Marti MA, Turjanski AG: TuberQ: a Mycobacterium tuberculosis protein druggability database. Database (Oxford). 2014, 2014 (0): bau035-10.1093/database/bau035.View ArticleGoogle Scholar
- Ruiz JC, D'Afonseca V, Silva A, Ali A, Pinto AC, Santos AR, et al: Evidence for reductive genome evolution and lateral acquisition of virulence functions in two Corynebacterium pseudotuberculosis strains. PLoS One. 2011, 6 (4): e18551-10.1371/journal.pone.0018551.PubMed CentralView ArticlePubMedGoogle Scholar
- Bilofsky HS, Burks C: The GenBank genetic sequence data bank. Nucleic Acids Res. 1988, 16 (5): 1861-1863. 10.1093/nar/16.5.1861.PubMed CentralView ArticlePubMedGoogle Scholar
- Sayers EW, Barrett T, Benson DA, Bolton E, Bryant SH, Canese K, et al: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2011, 39 (Database issue): D38-D51.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Li W, Jaroszewski L, Godzik A: Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics. 2001, 17 (3): 282-283. 10.1093/bioinformatics/17.3.282.View ArticlePubMedGoogle Scholar
- Eswar N, Webb B, Marti-Renom MA, Madhusudhan MS, Eramian D, Shen MY, Pieper U, Sali A: Comparative protein structure modeling using MODELLER. Curr Protoc Protein Sci. Edited by: John E Coligan [et al]. 2007, Chapter 2 (Unit 2.9):Google Scholar
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234 (3): 779-815. 10.1006/jmbi.1993.1626.View ArticlePubMedGoogle Scholar
- Melo F, Feytmans E: Assessing protein structures with a non-local atomic interaction energy. J Mol Biol. 1998, 277 (5): 1141-1152. 10.1006/jmbi.1998.1665.View ArticlePubMedGoogle Scholar
- Melo F, Sali A: Fold assessment for comparative protein structure modeling. Protein Sci. 2007, 16 (11): 2412-2426. 10.1110/ps.072895107.PubMed CentralView ArticlePubMedGoogle Scholar
- Benkert P, Kunzli M, Schwede T: QMEAN server for protein model quality estimation. Nucleic Acids Res. 2009, 37 (Web Server issue): W510-W514.PubMed CentralView ArticlePubMedGoogle Scholar
- Melo F, Sanchez R, Sali A: Statistical potentials for fold assessment. Protein Sci. 2002, 11 (2): 430-448.PubMed CentralView ArticlePubMedGoogle Scholar
- Velec HF, Gohlke H, Klebe G: DrugScore(CSD)-knowledge-based scoring function derived from small molecule crystal data with superior recognition rate of near-native ligand poses and better affinity prediction. Journal of medicinal chemistry. 2005, 48 (20): 6296-6303. 10.1021/jm050436v.View ArticlePubMedGoogle Scholar
- Soares SC, Abreu VA, Ramos RT, Cerdeira L, Silva A, Baumbach J, et al: PIPS: pathogenicity island prediction software. PLoS One. 2012, 7 (2): e30848-10.1371/journal.pone.0030848.PubMed CentralView ArticlePubMedGoogle Scholar
- Barh D, Gupta K, Jain N, Khatri G, Leon-Sicairos N, Canizalez-Roman A, et al: Conserved host-pathogen PPIs. Globally conserved inter-species bacterial PPIs based conserved host-pathogen interactome derived novel target in C. pseudotuberculosis, C. diphtheriae, M. tuberculosis, C. ulcerans, Y. pestis, and E. coli targeted by Piper betel compounds. Integr Biol (Camb). 2013, 5 (3): 495-509. 10.1039/c2ib20206a.View ArticleGoogle Scholar
- Barh D, Jain N, Tiwari S, Parida BP, D'Afonseca V, Li L, et al: A novel comparative genomics analysis for common drug and vaccine targets in Corynebacterium pseudotuberculosis and other CMN group of human pathogens. Chem Biol Drug Des. 2011, 78 (1): 73-84. 10.1111/j.1747-0285.2011.01118.x.View ArticlePubMedGoogle Scholar
- Zhang R, Ou HY, Zhang CT: DEG: a database of essential genes. Nucleic Acids Res. 2004, 32 (Database issue): D271-D272.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoon SH, Park YK, Lee S, Choi D, Oh TK, Hur CG, Kim JF: Towards pathogenomics: a web-based resource for pathogenicity islands. Nucleic Acids Res. 2007, 35 (Database issue): D395-D400.PubMed CentralView ArticlePubMedGoogle Scholar
- Magrane M, Consortium U: UniProt Knowledgebase: a hub of integrated protein data. Database (Oxford). 2011, 2011: bar009-View ArticleGoogle Scholar
- Yu CS, Lin CJ, Hwang JK: Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 2004, 13 (5): 1402-1406. 10.1110/ps.03479604.PubMed CentralView ArticlePubMedGoogle Scholar
- Garg SK, Alam MS, Kishan KV, Agrawal P: Expression and characterization of alpha-(1,4)-glucan branching enzyme Rv1326c of Mycobacterium tuberculosis H37Rv. Protein Expr Purif. 2007, 51 (2): 198-208. 10.1016/j.pep.2006.08.005.View ArticlePubMedGoogle Scholar
- Sanchez-Martinez M, Marcos E, Tauler R, Field M, Crehuet R: Conformational compression and barrier height heterogeneity in the N-acetylglutamate kinase. J Phys Chem B. 2013, 117 (46): 14261-14272. 10.1021/jp407016v.View ArticlePubMedGoogle Scholar
- Anishetty S, Pulimi M, Pennathur G: Potential drug targets in Mycobacterium tuberculosis through metabolic pathway analysis. Comput Biol Chem. 2005, 29 (5): 368-378. 10.1016/j.compbiolchem.2005.07.001.View ArticlePubMedGoogle Scholar
- Mathieu M, Debousker G, Vincent S, Viviani F, Bamas-Jacques N, Mikol V: Escherichia coli FolC structure reveals an unexpected dihydrofolate binding site providing an attractive target for anti-microbial therapy. J Biol Chem. 2005, 280 (19): 18916-18922. 10.1074/jbc.M413799200.View ArticlePubMedGoogle Scholar
- Pillai B, Cherney MM, Diaper CM, Sutherland A, Blanchard JS, Vederas JC, James MNG: Structural insights into stereochemical inversion by diaminopimelate epimerase: An antibacterial drug target. Proc Natl Acad Sci U S A. 2006, 103 (23): 8668-8673. 10.1073/pnas.0602537103.PubMed CentralView ArticlePubMedGoogle Scholar
- Jia DF: [Novel targets for antibiotics discovery: riboswitches]. Yao Xue Xue Bao. 2013, 48 (9): 1361-1368.PubMedGoogle Scholar
- Hassan SS, Tiwari S, Guimarães LC, Jamal SB, Folador E, Sharma NB, et al: Proteome Scale Comparative Modeling for Conserved Drug and Vaccine Targets Identification in Corynebacterium pseudotuberculosis. BMC Genomics. 2014, 15 (Suppl 7): S3-PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.