Impact Factor 4.151

Frontiers journals are at the top of citation and impact metrics

Original Research ARTICLE

Front. Genet., 13 February 2018 | https://doi.org/10.3389/fgene.2018.00044

The Druggable Pocketome of Corynebacterium diphtheriae: A New Approach for in silico Putative Druggable Targets

Syed S. Hassan1†, Syed B. Jamal2†, Leandro G. Radusky3†, Sandeep Tiwari2†, Asad Ullah1, Javed Ali4, Behramand1, Paulo V. S. D. de Carvalho2, Rida Shams1, Sabir Khan5, Henrique C. P. Figueiredo6, Debmalya Barh2,7, Preetam Ghosh8, Artur Silva9, Jan Baumbach10, Richard Röttger10, Adrián G. Turjanski3,11* and Vasco A. C. Azevedo2*
  • 1Department of Chemistry, Islamia College University Peshawar, Peshawar, Pakistan
  • 2PG Program in Bioinformatics, Institute of Biological Sciences, Federal University of Minas Gerais, Belo Horizonte, Brazil
  • 3Departamento de Química Biológica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
  • 4Department of Chemistry, Kohat University of Science and Technology, Kohat, Pakistan
  • 5Department of Analytical Chemistry, Institute of Chemistry, São Paulo State University, São Paulo, Brazil
  • 6AQUACEN, National Reference Laboratory for Aquatic Animal Diseases, Ministry of Fisheries and Aquaculture, Federal University of Minas Gerais, Belo Horizonte, Brazil
  • 7Centre for Genomics and Applied Gene Technology, Institute of Integrative Omics and Applied Biotechnology, Purba Medinipur, India
  • 8Department of Computer Science, Virginia Commonwealth University, Richmond, VA, United States
  • 9Institute of Biological Sciences, Federal University of Pará, Belém, Brazil
  • 10Department of Mathematics and Computer Science, University of Southern Denmark, Odense, Denmark
  • 11INQUIMAE/UBA-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina

Diphtheria is an acute and highly infectious disease, previously regarded as endemic in nature but vaccine-preventable, is caused by Corynebacterium diphtheriae (Cd). In this work, we used an in silico approach along the 13 complete genome sequences of C. diphtheriae followed by a computational assessment of structural information of the binding sites to characterize the “pocketome druggability.” To this end, we first computed the “modelome” (3D structures of a complete genome) of a randomly selected reference strain Cd NCTC13129; that had 13,763 open reading frames (ORFs) and resulted in 1,253 (∼9%) structure models. The amino acid sequences of these modeled structures were compared with the remaining 12 genomes and consequently, 438 conserved protein sequences were obtained. The RCSB-PDB database was consulted to check the template structures for these conserved proteins and as a result, 401 adequate 3D models were obtained. We subsequently predicted the protein pockets for the obtained set of models and kept only the conserved pockets that had highly druggable (HD) values (137 across all strains). Later, an off-target host homology analyses was performed considering the human proteome using NCBI database. Furthermore, the gene essentiality analysis was carried out that gave a final set of 10-conserved targets possessing highly druggable protein pockets. To check the target identification robustness of the pipeline used in this work, we crosschecked the final target list with another in-house target identification approach for C. diphtheriae thereby obtaining three common targets, these were; hisE-phosphoribosyl-ATP pyrophosphatase, glpX-fructose 1,6-bisphosphatase II, and rpsH-30S ribosomal protein S8. Our predicted results suggest that the in silico approach used could potentially aid in experimental polypharmacological target determination in C. diphtheriae and other pathogens, thereby, might complement the existing and new drug-discovery pipelines.

Introduction

Corynebacterium diphtheriae belong to the class Actinomycetales and is a Gram-positive bacterium, a non-spore forming, non-motile and facultative anaerobe with pleomorphic cell shape and high GC content (∼53%) relative to the Firmicutes (Cerdeno-Tarraga et al., 2003; Trost et al., 2012). This bacterium is the causative agent of diphtheria, a severe human respiratory manifestation characterized by bacterial adhesion to host pharyngeal cell (pharyngitis and pseudomembranous inflammation). The pathogen target one or both tonsils that further disseminate at a later stage thereby resulting in complete airway obstruction and death (Hadfield et al., 2000). The cornerstone in diphtheria therapy involves the hyper immune antiserum-antitoxin produced in equines that neutralizes the C. diphtheriae toxin while among antibiotics are the broad-spectrum penicillin or erythromycin. However, recent emergence of numerous drug-resistant strains subsequently decreases the efficacy of current therapeutics (Barraud et al., 2011). Importantly, the World Health Organization recommends widespread DTPs immunization against toxigenic diphtheria strains as the only effective approach to counteract the infection. Although, a slight change in mortality has been observed since the availability and administration of antibiotics, specifically penicillin or erythromycin (Adler et al., 2013). C. diphtheriae has four biotypes: gravis, mitis, intermedius and belfanti that are non-sporulating, unencapsulated, non-motile and pleomorphic bacilli. They are subsequently classified on the basis of differences in colonial morphology, haemolytic potential, fermentation reactions and severity of the resulting disease (Gerald et al., 2009; Brooks et al., 2010). The infant mortality rate in an epidemic that resulted in thousands of casualties dropped gradually in countries where living standards were improved and immunization programs were introduced. Nevertheless, it still remains a significant pathogen around the globe (Hodes, 1979). The ‘strangling angel’ effects on children that scaled from wing-shaped disarticulation and pseudo-membranes formation in the oropharynx, triggered acute obstruction of airways and resulted in mortalities (Hodes, 1979; Hart et al., 1996; Jamal et al., 2017b). A plethora of cases were recently reported and still expected from both non-lethal and lethal diphtheria in different parts of the world due to significant population displacements via immigration. An adequate supervision necessitates quick measures to discover additional diphtheria antitoxin, antibiotic and therapeutic treatment (Pizza et al., 2000).

The emerging concepts of polypharmacology, differential genome analyses, and reverse vaccinology, comparative and subtractive microbial genomics have largely contributed by establishing complementary traditions for fast identification of novel targets in post-genomic era (Perumal et al., 2007; Barh et al., 2013). Comparative homology modeling (Baumbach, 2010; Rottger et al., 2013) has widely been used in expanding the structural space of pathogens (Chong et al., 2006; Asif et al., 2009).

These practices are being used for the identification of conserved targets in a several human and animal pathogens like C. tuberculosis (Hassan et al., 2014; Radusky et al., 2015), Mycobacterium. tuberculosis (Asif et al., 2009), Burkholderia pseudomallei (Chong et al., 2006), Neisseria gonorrhoea (Barh and Kumar, 2009), Helicobacter pylori (Dutta et al., 2006), Pseudomonas aeruginosa (Sakharkar et al., 2004; Perumal et al., 2007), and Salmonella typhi (Rathi et al., 2009).

In this work, a combination of in silico tools was primarily used to predict the core proteome of C. diphtheriae species to associate genomic information based on the 3D structures. The predicted proteomes were modeled (pan-modelome) using a methodology adapted by Hassan et al. (2014). From a structural point of view, druggability is the probability of small drug-like molecules binding to a given target protein with high affinity (<1 μM). We report for the first time the structural druggability assessment for multi-strain C. diphtheriae proteomes using a pan-druggability prediction pipeline based on the open source pocket detection code “fpocket”. The method integrates several physicochemical descriptors to estimate the pocket druggability on a genomic scale with suitable features that enable binding of a drug-like compound (Kinnings et al., 2010).

Materials and Methods

Initial Dataset Construction

All ORFs (Open Reading Frames) of the 13 completely sequenced genomes of C. diphtheriae were obtained from the NCBI database1. Table 1 shows the statistical data of all strains used in this study where the strain NCTC13129 has 2,272 reported ORFs and was randomly selected as a reference genome for modelome prediction and further analyses.

TABLE 1
www.frontiersin.org

TABLE 1. Summary of Corynebacterium diphtheriae strains used in this study and their respective modeling statistics for druggability analyses.

General Concept: Modelome Prediction

The binding affinity of small drug-like molecules to the active site of putative biological targets (druggable protein cavities) formulated a basis for this work, a slightly modified protocol of Radusky et al. (2015) (Figure 1). All genome ORF sequences of the 13 C. diphtheriae strains were subjected to the MHOLline workflow2 in.faa file format for 3D structure prediction. MHOLline utilizes multi fasta files of amino acids as an input data and then uses HMMTOP, BLAST, BATS, Modeller and Procheck programs for the detailed analyses. The program HMMTOP detects transmembrane regions. The BLAST algorithm is used to identify template structure by performing a random search against the Protein Data Bank. BATS (Blast Automatic Targeting for Structures) carry out the refinement in the template search; it is a key step for the model construction. BATS refinement identifies sequences that make the modeling possible by selecting a template from BLAST output file using their BATS scores, expectation values, identity and sequence similarity as criteria as well as considering the number of gaps and the alignment coverage. BATS select the best template for 3D model generation and perform automated alignment used by the MODELLER program. Furthermore, it gathers all the BLAST output files into four distinctive groups, i.e., G0, G1, G2, and G3, according to the following criteria; G0 = Not aligned sequence, G1 = E-value > 10e-5 or Identity < 15%, G2 = E-value ≤ 10e-5 and Identity ≥ 25% AND LVI ≤ 0.7, G3 = E-value ≤ 10e-5 and Identity ≤ 15% and <25% OR LVI > 0.7, Where LVI is the Length Variation Index, a MHOLline concept of coverage (LVI ≤ 0.1 is equivalent to a coverage ≥ 90%). Once the template is selected based on BATS results, MODELLER program is used for the generation of 3D protein model. There is no significant correlation statistically, between the number of templates used during model building and the overall quality of a model. In the next step, another MHOLline tool called FILTERS, categorizes the BATS selected sequences (G2) into distinct quality model subgroups, based on identity and LVI value. The subgroups ranges Very High to Very Low. To evaluate the overall quality and accuracy of the model, Ramachandran plot is obtained which explained the stereochemical quality of the model. Precisely, the MHOLine generates an aggregate structural information for all the submitted sequences in the fasta format, Ramachandran plot and other properties like structural quality and enzymatic functions are also determined. Further details can be obtained by visiting MHOLline homepage (Hassan et al., 2014; Webb and Sali, 2016; Jamal et al., 2017a). For all modeled structures, structural properties were figured as: (i) the Druggability Score (DS) for each pocket and (ii) the active site residues (if available) according to the template structures available at the protein databank RCSB-PDB3 (Berman et al., 2003).

FIGURE 1
www.frontiersin.org

FIGURE 1. Overview of pipeline components for data analysis.

3D Protein Models in Non-reference Strains: Orthologs Identification

We used an applied bioinformatics procedure to find the conserved putative druggable targets across all the 13 C. diphtheriae strains at genome-scale by first predicting their 3D models. The ensemble methodology essentially is a filter of thousands of candidate genes to yield high-confidence 3D structural models from orthologous proteins in C. diphtheriae species. As aforementioned, the MHOLline resulted in 1,253 predicted structures for the randomly selected reference strain NCTC13129 that were later used as template structures for modeling the 3D structures in non-reference strains as well. Further, the BLASTp program was installed in a local machine and used to check if the ORFs of the reference proteome have orthologs in the remaining 12 strains using the following parameters; identity ≥ 85%, coverage > 80%. The protein sequences showing high identity values (> = 85%) for each reference and non-reference strain were considered as conserved and the modeled structures of reference strain were used again as templates to predict the 3D models for the aforementioned 12 non-reference strains. The core modelome was compared and evaluated for the quality of the obtained 3D structures. A reliable model has a probability of correct fold larger than 95% and coverage of over 50% with the template structure. For each sequence in the reference and non-reference strain that gave an identity hit of >85%, a mutation methodology was applied on each amino acid substitution using the MODELLER program. These models were then used to compute the druggability variation for the 13 strains of the C. diphtheriae species.

Evaluation of Structural Druggability

The protein structural druggability of each predicted 3D model was evaluated by determining the ability of putative pockets to bind drug-like molecule/s, using the fpocket (Finn et al., 2016) and the recently developed DrugScore (DS) programs (Velec et al., 2005). The later methodology is based on the Voronoi tessellation algorithm that identify and characterize pockets and compute suitable physicochemical descriptors such as polar and apolar surface area, hydrophobic-hydrophobic density and polarity score. In conjunction they yield a druggability value that ranges between 0 (non-druggable, ND) and 1 (highly druggable, HD). We categorized the druggability scores for all predicted 3D structures into four sets: non-druggable (ND; DS ≤ 0.2), poorly druggable (PD; DS ≥ 0.2 and ≤0.5), druggable (D; DS ≥ 0.5 and ≤0.7), and highly druggable (HD; DS ≥ 0.7) protein pockets, respectively. This distribution is in accordance to our previous work where the druggability score was computed for all pockets present in all unique proteins in the Protein Data Bank that were experimentally crystallized in complex with a drug-like compound/s (Radusky et al., 2014).

Identification of Active Site Residues

In order to identify the amino acid residues in the active site of the predicted druggable protein pocket/s, information were retrieved from the CSA database (Catalytic Site Atlas) (Furnham et al., 2014) and Pfam position site (Finn et al., 2016), respectively. A list of PDB_IDs was generated linked to a number of residues constituting the corresponding protein active sites. To map the active site residues to as many C. diphtheriae proteins as possible, each PDB_ID was used as a template in CSA and assigned to the modeled ORFs.

Host Homology, Essentiality, and Core-Modelomics of the Selected Targets

For off target prediction, the pool of global druggable (GD) proteins was piped into NCBI-BLASTp using default parameters (identity = 0% and/or no hit) against the human proteome to identify non-host homologs. Moreover, from the filtered list of 10 highly druggable non-host homologous target proteins, an approach based on subtractive genomics was implemented and applied to the GD targets that were essential to bacteria (Barh et al., 2011). Briefly, the set of target proteins of C. diphtheriae was submitted to the Database of Essential Genes (DEG, which contains experimentally validated essential genes from bacteria, archaea and eukaryotes) for homology analyses (Zhang et al., 2004). Again, we used BLASTp with E-value cut-off of 1e-05 a bit score ≥ 100 and identity ≥ 35% (Barh et al., 2011). The final list of putative targets based on criteria described earlier, contained 10 essential and non-host homologous target proteins. The obtained list was further subjected to ProtParam4 for molecular weight determination, biochemical pathway analysis to KEGG (Kyoto Encyclopedia of Genes and Genomes) (Kanehisa and Goto, 2000) using network enrichment (Alcaraz et al., 2012), virulence using PAIDB (Pathogenicity Island Database) (Yoon et al., 2007), functionality using UniProt (Universal Protein Resource) (Magrane and UniProt, 2011), and cellular localization using CELLO (subCELlular LOcalization predictor) (Yu et al., 2004). In addition, we merged the final set of 10 selected non-host homologous, essential and global druggable proteins with results obtained through experiments locally performed in our laboratories (Jamal et al., unpublished data) resulting in three common targets, which we selected as candidates.

Protein–Protein Interaction Network

In biological systems, proteins work in a homogenous environment rather than individual, hence it is important to study protein–protein interactions (PPIs) for C. diphtheriae metabolism. The identified drug targets were evaluated to study their potential biological, functional and metabolic roles for proteomic interactions. The selected drug targets were used to develop intra-species protein–protein interactome using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database (Szklarczyk et al., 2015). STRING is an online network analyses tool that provides essential information regarding interactions of the desired proteins.

Results and Discussion

Prediction of Structural Homology Based Models

The complete modelome of the reference strain NCTC13129 was computed; consisting of 13,763 ORF, with 1,253 (∼9%) resulting models. Taking the original models of reference strain NCTC13129 as templates, we then generated 438 conserved models in the 12 remaining strains using the MODELLER software (Sali and Blundell, 1993). Afterward, the target-template alignments have been computed using a BLAST E-value cut-off of 10-6 in order to build the model structures using the MODELLER software (Sali and Blundell, 1993; Webb and Sali, 2016). For each target-template alignment, ten different target models were built, and their quality measures have been assessed using GA341 (Melo and Feytmans, 1998; Melo and Sali, 2007) and QMEAN (Benkert et al., 2009), keeping models with GA341 reliability scores ≥ 0.7 (Melo et al., 2002), leading to a final set of 401 protein models. All these proteins are tabulated in Supplementary Table S1.

Pocketome Druggability and Active Site Residues of C. diphtheriae

The list of 401 targets protein drastically reduced to 137 after druggability analyses using the aforementioned fpocket and the recently developed DS programs. A summary of only highly druggable (HD) targets with drug score remained ≥ 0.7 were considered as global druggable. The calculated structural druggability scores are given in Supplementary Table S2.

In Figure 2, a comparison of calculated druggability score distribution across all structures of C. diphtheriae reference and other strains is shown. Although the distribution has a small shift to higher values, we used the same bounds to define the sets of druggable proteins (Figure 2). A protein target, which remained druggable in all strains, was classified as Globally Druggable (GD).

FIGURE 2
www.frontiersin.org

FIGURE 2. Histogram of druggability score along the Corynebacterium diphtheriae reference/other strains. Representation of all ligand-bound structures in the PDB (blue pointed line), all modeled structures of C. diphtheriae reference strain (red line, 1,253 models), and all modeled structures of the C. diphtheriae non-reference genome strains (green line, 401 core proteins in 12 non-reference strains).

Non-host Homology, Essentiality and Core-Modelomics Analyses

As aforementioned, the list of 137 global druggable proteins (Supplementary Table S2) was computed to the corresponding human host proteome that resulted in the identification of a very small set of only 10 non-host homologous proteins; remaining 127 as host-homologous. The non-host homologous targets were selected following a very stringent criterion, i.e., no identity, no hits. This list of final 10 essential and non-host homologous targets in C. diphtheriae is given in Table 2. We further report the involvement of these putative targets in metabolic pathways, biological processes, cellular localization, molecular weights and most importantly their potential role as virulence factors. Out of 10 targets, 7 targets were found as pathogen virulence factors using the PAIDB database based on homology method. Further, we predicted the subcellular localization of these final target proteins using support vector machines, a methodology that is based on n-peptide composition of proteins, adapted in a related work by Yu et al. (2004), we obtained three high-confidential candidates, which are; hisE, glpX, and rpsH. Interestingly, these three high-confidential candidates were identified as essential and non-host homologous targets in our previous work by Jamal et al. (2017a). These proteins were subjected to molecular docking analysis against four different ligand libraries and a set of some potent molecules were suggested for active inhibition of these proteins (Jamal et al., 2017a).

TABLE 2
www.frontiersin.org

TABLE 2. List of global druggable, non-host homologous and essential putative targets, their functional annotation together with other information.

hisE (Phosphoribosyl-ATP Pyrophosphatase)

hisE is the second enzyme in histidine-biosynthetic pathway hydrolysing irreversibly phosphoribosyl-ATP to phosphoribosyl-AMP and pyrophosphate. The protein is encoded by the hisE gene, fused to hisI in many bacteria, fungi and plants but is present as a separate gene in some bacteria and archaea. Since it is seen in in vitro experiments that hisE is essential for microorganism growth, we assume it a potential drug target in C. diphtheriae. It is also reported as a drug candidate for tuberculosis (Javid-Majd et al., 2008).

glpX (Fructose 1,6-Bisphosphatase II)

It is one of the main enzyme for gluconeogenesis that catalyses the hydrolysis of fructose 1,6-bisphosphate to form fructose 6-phosphate and orthophosphate. In glycolysis, phosphofructokinase catalysis the reverse reaction, and the product, fructose 6-phosphate, are important precursors in various biosynthetic pathways (Horecker et al., 1975). Gluconeogenesis is an important metabolic pathway in all organisms and plays a key role by allowing the cells to synthesize glucose from non-carbohydrate precursors, such as glycerol, organic acids and amino acids. FBPases are members of lithium sensitive phosphatases a large superfamily which includes three families of inositol phosphatases and FBPases (phosphoesterase clan CL0171, AA sequences 3167 from Pfam data base). They are already reported as target for the treatment of non-insulin dependent diabetes and development of new drugs (Wright et al., 2002; Sassetti and Rubin, 2003).

rpsH (30S Ribosomal Protein S8)

The protein rpsH is one of the key RNA-binding protein having a central position within the small ribosomal subunit. It interacts widely with 16S rRNA and is fundamental for the correct folding of the central domain of the ribosomal rRNA. Furthermore, this protein regulates the synthesis of various other ribosomal proteins by binding to mRNA. In the two RNA molecules, it binds exactly to very similar sites. rpsH has a medium size and recently it has been discovered that rpsH play vital role as a significant primary RNA-binding protein in the 30S subunit. Mutations in S8 within the protein are shown to result in defective ribosomal assembly. The S8-binding site within 16S rRNA in Escherichia coli, has been investigated independently by a number of techniques including protein crosslinking, nuclease protection, hydroxyl-radical foot printing, RNA–RNA modification and chemical probing. The 30S ribosomal protein S8 is also one of the principal regulatory elements that control ribosomal protein synthesis by the translational feedback inhibition mechanism discovered by Yates et al. (1980). It regulates the expression of spc operon that encodes the 10 ribosomal proteins L5, L6, L14, L15, L18, L24, L30, S5, S8, and S14, respectively (Davies et al., 1996).

Protein–Protein Interaction Network for Proposed Targets

Protein–protein interaction of target proteins with each other have been constructed showing two proteins, rplJ (VN94_02905/50S ribosomal protein L10) and rpsH (30S ribosomal protein S8) to be interacting directly based on highest confidence score 0.9. The confidence score is the approximate probability that a predicted link exists between two enzymes in the same metabolic map in the KEGG database. The evidences for this interaction are gene fusion, co-occurrence, co-expression, experimental and databases (Supplementary Figure S1).

Conclusion

We performed a comprehensive in silico study of the druggability scores on all sequenced genomes of C. diphtheria resulting in a list of intra-strain highly druggable pockets of 10 ORFs non-homologous in human hosts. Previously, we have implemented a similar approach using other bioinformatics tools for the identification of putative therapeutic targets in C. diphtheriae that relied primarily on the modelome construction followed by filtering the obtained data for conserved targets (Jamal et al., 2017a). In that work, a final set of eight essential and non-host homologs targets were subjected to virtual screening using different compound libraries but lacked a detailed overview of the druggable protein pockets of the selected targets. Here, we further extrapolated our work to the druggable pocketome at species level and then at the end compared our final data set obtained in this work with the aforementioned published data. The comparison showed that any of the two approaches for putative targets identification in pathogenic microorganisms might provide an easy-to-handle protocol in future drug discovery projects. Our pipeline is expandable and can be applied to other bacterial species as well. In the future, we will work on Cytoscape plugins to allow for mapping essential druggable non-homologous genes to biological networks interactively for follow-up systems biology investigation (Baumbach and Apeltsin, 2008). We believe that our approach has the potential to aid in designing drugs and/or vaccines, and in developing protein inhibitors as well as discovering new lead compounds.

Author Contributions

SH, SJ, LR, ST planned the whole work. SH, SJ, LR, ST, and PdC analyzed the data. SH, SJ, LR, ST drafted the manuscript. SH, HF, VA, AS, B, DB, PG, JB, RR, and AT reviewed and analyzed the manuscript. RS, PdC, SJ, and ST performed the literature review and formatting the tables/figures. AU, SK, JB, RR, and JA provided useful comments/suggestions for the improvement of the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work is a part of a mutual collaboration among the two leading research/academic institutions of Brazil and Argentine, 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, and the Structural Bioinformatics Group, Institute of Physical Chemistry of Materials, Environment and Energy, University of Buenos Aires, Argentine, respectively. We acknowledge the support of all co-authors who have worked and collaborated in accomplishing this project. We are grateful to our Argentinian collaborators, especially Prof. Dr. Adrian Gustavo Turjanski and Leandro G. Radusky for their collaboration, timely assistance in theoretical/biological discussion about the project and manuscript preparation. We acknowledge the Brazilian funding agencies CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brasil). Syed Babar Jamal acknowledges the “TWAS-CNPq Postgraduate Fellowship Program” for granting a fellowship for doctoral studies.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00044/full#supplementary-material

FIGURE S1 | Protein–protein interactome for the identified common conserved targets. Ribosomal pathway protein (VN94_02915 rplJ, rpsH) was interacting with each other either directly or indirectly. The lines color indicates evidences of interactions that were predicted or experimentally validated.

TABLE S1 | Tabulated are the 401 conserved target proteins with appropraite 3D structures that were used for further druggability etc. analyse.

TABLE S2 | Summary of core Global Druggable (GD) targets along the 13 strains of C. diphtheriae with predicted druggability scores (137 GD targets).

Footnotes

  1. ^ https://www.ncbi.nlm.nih.gov/genome/genomes/1025?
  2. ^ http://www.mholline.lncc.br
  3. ^ www.wwpdb.org
  4. ^ http://web.expasy.org/protparam/

References

Adler, N. R., Mahony, A., and Friedman, N. D. (2013). Diphtheria: forgotten, but not gone. Intern. Med. J. 43, 206–210. doi: 10.1111/imj.12049

PubMed Abstract | CrossRef Full Text | Google Scholar

Alcaraz, N., Friedrich, T., Kotzing, T., Krohmer, A., Muller, J., Pauling, J., et al. (2012). Efficient key pathway mining: combining networks and OMICS data. Integr. Biol. 4, 756–764. doi: 10.1039/c2ib00133k

PubMed Abstract | CrossRef Full Text | Google Scholar

Asif, S. M., Asad, A., Faizan, A., Anjali, M. S., Arvind, A., Neelesh, K., et al. (2009). Dataset of potential targets for Mycobacterium tuberculosis H37Rv through comparative genome analysis. Bioinformation 4, 245–248. doi: 10.6026/97320630004245

PubMed Abstract | CrossRef Full Text | Google Scholar

Barh, D., Gupta, K., Jain, N., Khatri, G., Leon-Sicairos, N., Canizalez-Roman, A., et al. (2013). 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. 5, 495–509. doi: 10.1039/c2ib20206a

PubMed Abstract | CrossRef Full Text | Google Scholar

Barh, D., Jain, N., Tiwari, S., Parida, B. P., D’Afonseca, V., Li, L., et al. (2011). 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. 78, 73–84. doi: 10.1111/j.1747-0285.2011.01118.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Barh, D., and Kumar, A. (2009). In silico identification of candidate drug and vaccine targets from various pathways in Neisseria gonorrhoeae. In Silico Biol. 9, 225–231.

PubMed Abstract | Google Scholar

Barraud, O., Badell, E., Denis, F., Guiso, N., and Ploy, M. C. (2011). Antimicrobial drug resistance in Corynebacterium diphtheriae mitis. Emerg. Infect. Dis. 17, 2078–2080. doi: 10.3201/eid1711.110282

PubMed Abstract | CrossRef Full Text | Google Scholar

Baumbach, J. (2010). On the power and limits of evolutionary conservation–unraveling bacterial gene regulatory networks. Nucleic Acids Res. 38, 7877–7884. doi: 10.1093/nar/gkq699

PubMed Abstract | CrossRef Full Text | Google Scholar

Baumbach, J., and Apeltsin, L. (2008). Linking cytoscape and the corynebacterial reference database CoryneRegNet. BMC Genomics 9:184. doi: 10.1186/1471-2164-9-184

PubMed Abstract | CrossRef Full Text | Google Scholar

Benkert, P., Kunzli, M., and Schwede, T. (2009). QMEAN server for protein model quality estimation. Nucleic Acids Res. 37, W510–W514. doi: 10.1093/nar/gkp322

PubMed Abstract | CrossRef Full Text | Google Scholar

Berman, H., Henrick, K., and Nakamura, H. (2003). Announcing the worldwide protein data bank. Nat. Struct. Biol. 10:980. doi: 10.1038/nsb1203-980

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. F., Carroll, K. C., Butel, J. S., Morse, S. A., and Mietzneron, T. A. (eds) (2010). “Aerobic non–spore-forming gram-positive bacilli: Corynebacterium, Listeria, Erysipelothrix, Actinomycetes, and related pathogens,” in Jawetz, Melnick, & Adelberg’s Medical Microbiology, 25 Edn, (New York, NY: McGraw-Hill).

Google Scholar

Cerdeno-Tarraga, A. M., Efstratiou, A., Dover, L. G., Holden, M. T., Pallen, M., Bentley, S. D., et al. (2003). The complete genome sequence and analysis of Corynebacterium diphtheriae NCTC13129. Nucleic Acids Res. 31, 6516–6523. doi: 10.1093/nar/gkg874

CrossRef Full Text | Google Scholar

Chong, C. E., Lim, B. S., Nathan, S., and Mohamed, R. (2006). In silico analysis of Burkholderia pseudomallei genome sequence for potential drug targets. In Silico Biol. 6, 341–346.

PubMed Abstract | Google Scholar

Davies, C., Ramakrishnan, V., and White, S. W. (1996). Structural evidence for specific S8-RNA and S8-protein interactions within the 30S ribosomal subunit: ribosomal protein S8 from Bacillus stearothermophilus at 1.9 A resolution. Structure 4, 1093–1104. doi: 10.1016/S0969-2126(96)00115-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Dutta, A., Singh, S. K., Ghosh, P., Mukherjee, R., Mitter, S., and Bandyopadhyay, D. (2006). In silico identification of potential therapeutic targets in the human pathogen Helicobacter pylori. In Silico Biol. 6, 43–47.

Google Scholar

Finn, R. D., Coggill, P., Eberhardt, R. Y., Eddy, S. R., Mistry, J., Mitchell, A. L., et al. (2016). The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 44, D279–D285. doi: 10.1093/nar/gkv1344

PubMed Abstract | CrossRef Full Text | Google Scholar

Furnham, N., Holliday, G. L., de Beer, T. A., Jacobsen, J. O., Pearson, W. R., and Thornton, J. M. (2014). The catalytic site atlas 2.0: cataloging catalytic sites and residues identified in enzymes. Nucleic Acids Res. 42, D485–D489. doi: 10.1093/nar/gkt1243

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerald, L., Mandell, J. E. B., and Raphael, D. (2009). Mandell, Douglas, and Bennett’s Principles and Practice of Infectious Diseases. Philadelphia, PA: ELSEVIER.

Google Scholar

Hadfield, T. L., McEvoy, P., Polotsky, Y., Tzinserling, V. A., and Yakovlev, A. A. (2000). The pathology of diphtheria. J. Infect. Dis. 181(Suppl. 1), S116–S120. doi: 10.1086/315551

PubMed Abstract | CrossRef Full Text | Google Scholar

Hart, P. E., Lee, P. Y., Macallan, D. C., and Wansbrough-Jones, M. H. (1996). Cutaneous and pharyngeal diphtheria imported from the Indian subcontinent. Postgrad. Med. J. 72, 619–620. doi: 10.1136/pgmj.72.852.619

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, S. S., Tiwari, S., Guimaraes, L. C., Jamal, S. B., Folador, E., Sharma, N. B., et al. (2014). Proteome scale comparative modeling for conserved drug and vaccine targets identification in Corynebacterium pseudotuberculosis. BMC Genomics 15(Suppl. 7):S3. doi: 10.1186/1471-2164-15-S7-S3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodes, H. L. (1979). Diphtheria. Pediatr. Clin. North Am. 26, 445–459. doi: 10.1016/S0031-3955(16)33716-6

CrossRef Full Text

Horecker, B. L., Melloni, E., and Pontremoli, S. (1975). Fructose 1,6-bisphosphatase: properties of the neutral enzyme and its modification by proteolytic enzymes. Adv. Enzymol. Relat. Areas Mol. Biol. 42, 193–226.

PubMed Abstract | Google Scholar

Jamal, S. B., Hassan, S. S., Tiwari, S., Viana, M. V., Benevides, L. J., Ullah, A., et al. (2017a). An integrative in-silico approach for therapeutic target identification in the human pathogen Corynebacterium diphtheriae. PLOS ONE 12:e0186401. doi: 10.1371/journal.pone.0186401

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamal, S. B., Tiwari, S., Silva, A., and Azevedo, V. (2017b). Pathogenesis of Corynebacterium diphtheriae and available vaccines; an overview. Glob. J. Infect. Dis. Clin. Res. 3, 20–24.

Google Scholar

Javid-Majd, F., Yang, D., Ioerger, T. R., and Sacchettini, J. C. (2008). The 1.25 A resolution structure of phosphoribosyl-ATP pyrophosphohydrolase from Mycobacterium tuberculosis. Acta Crystallogr. D Biol. Crystallogr. 64(Pt 6), 627–635. doi: 10.1107/S0907444908007105

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27

CrossRef Full Text | Google Scholar

Kinnings, S. L., Xie, L., Fung, K. H., Jackson, R. M., Xie, L., and Bourne, P. E. (2010). The Mycobacterium tuberculosis drugome and its polypharmacological implications. PLOS Comput. Biol. 6:e1000976. doi: 10.1371/journal.pcbi.1000976

PubMed Abstract | CrossRef Full Text | Google Scholar

Magrane, M., and UniProt, C. (2011). UniProt Knowledgebase: a hub of integrated protein data. Database 2011:bar009. doi: 10.1093/database/bar009

PubMed Abstract | CrossRef Full Text | Google Scholar

Melo, F., and Feytmans, E. (1998). Assessing protein structures with a non-local atomic interaction energy. J. Mol. Biol. 277, 1141–1152. doi: 10.1006/jmbi.1998.1665

PubMed Abstract | CrossRef Full Text | Google Scholar

Melo, F., and Sali, A. (2007). Fold assessment for comparative protein structure modeling. Protein Sci. 16, 2412–2426. doi: 10.1110/ps.072895107

PubMed Abstract | CrossRef Full Text | Google Scholar

Melo, F., Sanchez, R., and Sali, A. (2002). Statistical potentials for fold assessment. Protein Sci. 11, 430–448. doi: 10.1002/pro.110430

PubMed Abstract | CrossRef Full Text | Google Scholar

Perumal, D., Lim, C. S., Sakharkar, K. R., and Sakharkar, M. K. (2007). Differential genome analyses of metabolic enzymes in Pseudomonas aeruginosa for drug target identification. In Silico Biol. 7, 453–465.

PubMed Abstract | Google Scholar

Pizza, M., Scarlato, V., Masignani, V., Giuliani, M. M., Arico, B., Comanducci, M., et al. (2000). Identification of vaccine candidates against serogroup B meningococcus by whole-genome sequencing. Science 287, 1816–1820. doi: 10.1126/science.287.5459.1816

CrossRef Full Text | Google Scholar

Radusky, L., Defelipe, L. A., Lanzarotti, E., Luque, J., Barril, X., Marti, M. A., et al. (2014). TuberQ: a Mycobacterium tuberculosis protein druggability database. Database 2014, bau035. doi: 10.1093/database/bau035

PubMed Abstract | CrossRef Full Text | Google Scholar

Radusky, L. G., Hassan, S., Lanzarotti, E., Tiwari, S., Jamal, S., Ali, J., et al. (2015). An integrated structural proteomics approach along the druggable genome of Corynebacterium pseudotuberculosis species for putative druggable targets. BMC Genomics 16(Suppl. 5):S9. doi: 10.1186/1471-2164-16-S5-S9

PubMed Abstract | CrossRef Full Text | Google Scholar

Rathi, B., Sarangi, A. N., and Trivedi, N. (2009). Genome subtraction for novel target definition in Salmonella typhi. Bioinformation 4, 143–150. doi: 10.6026/97320630004143

PubMed Abstract | CrossRef Full Text | Google Scholar

Rottger, R., Kalaghatgi, P., Sun, P., Soares Sde, C., Azevedo, V., Wittkop, T., et al. (2013). Density parameter estimation for finding clusters of homologous proteins–tracing actinobacterial pathogenicity lifestyles. Bioinformatics 29, 215–222. doi: 10.1093/bioinformatics/bts653

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakharkar, K. R., Sakharkar, M. K., and Chow, V. T. (2004). A novel genomics approach for the identification of drug targets in pathogens, with special reference to Pseudomonas aeruginosa. In Silico Biol. 4, 355–360.

PubMed Abstract | Google Scholar

Sali, A., and Blundell, T. L. (1993). Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 234, 779–815. doi: 10.1006/jmbi.1993.1626

PubMed Abstract | CrossRef Full Text | Google Scholar

Sassetti, C. M., and Rubin, E. J. (2003). Genetic requirements for mycobacterial survival during infection. Proc. Natl. Acad. Sci. U.S.A. 100, 12989–12994. doi: 10.1073/pnas.2134250100

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447–D452. doi: 10.1093/nar/gku1003

PubMed Abstract | CrossRef Full Text | Google Scholar

Trost, E., Blom, J., Soares Sde, C., Huang, I. H., Al-Dilaimi, A., Schroder, J., et al. (2012). Pangenomic study of Corynebacterium diphtheriae that provides insights into the genomic diversity of pathogenic isolates from cases of classical diphtheria, endocarditis, and pneumonia. J. Bacteriol. 194, 3199–3215. doi: 10.1128/JB.00183-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Velec, H. F., Gohlke, H., and Klebe, G. (2005). 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. J. Med. Chem. 48, 6296–6303. doi: 10.1021/jm050436v

PubMed Abstract | CrossRef Full Text | Google Scholar

Webb, B., and Sali, A. (2016). Comparative protein structure modeling using MODELLER. Curr. Protoc. Protein Sci. 86, 2.9.1–2.9.37. doi: 10.1002/cpps.20

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, S. W., Carlo, A. A., Carty, M. D., Danley, D. E., Hageman, D. L., Karam, G. A., et al. (2002). Anilinoquinazoline inhibitors of fructose 1,6-bisphosphatase bind at a novel allosteric site: synthesis, in vitro characterization, and X-ray crystallography. J. Med. Chem. 45, 3865–3877. doi: 10.1021/jm010496a

PubMed Abstract | CrossRef Full Text | Google Scholar

Yates, J. L., Arfsten, A. E., and Nomura, M. (1980). In vitro expression of Escherichia coli ribosomal protein genes: autogenous inhibition of translation. Proc. Natl. Acad. Sci. U.S.A. 77, 1837–1841. doi: 10.1073/pnas.77.4.1837

CrossRef Full Text | Google Scholar

Yoon, S. H., Park, Y. K., Lee, S., Choi, D., Oh, T. K., Hur, C. G., et al. (2007). Towards pathogenomics: a web-based resource for pathogenicity islands. Nucleic Acids Res. 35, D395–D400. doi: 10.1093/nar/gkl790

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, C. S., Lin, C. J., and Hwang, J. K. (2004). Predicting subcellular localization of proteins for Gram-negative bacteria by support vector machines based on n-peptide compositions. Protein Sci. 13, 1402–1406. doi: 10.1110/ps.03479604

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Ou, H. Y., and Zhang, C. T. (2004). DEG: a database of essential genes. Nucleic Acids Res. 32, D271–D272. doi: 10.1093/nar/gkh024

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Corynebacterium diphtheria, pocketome, druggable genome, structural proteomics, putative therapeutic targets, highly druggable (HD), global druggable (GD)

Citation: Hassan SS, Jamal SB, Radusky LG, Tiwari S, Ullah A, Ali J, Behramand, de Carvalho PVSD, Shams R, Khan S, Figueiredo HCP, Barh D, Ghosh P, Silva A, Baumbach J, Röttger R, Turjanski AG and Azevedo VAC (2018) The Druggable Pocketome of Corynebacterium diphtheriae: A New Approach for in silico Putative Druggable Targets. Front. Genet. 9:44. doi: 10.3389/fgene.2018.00044

Received: 23 November 2017; Accepted: 30 January 2018;
Published: 13 February 2018.

Edited by:

Thiago Motta Venancio, State University of Norte Fluminense, Brazil

Reviewed by:

A. Maxwell Burroughs, National Center for Biotechnology Information (NIH), United States
Tiago Sobreira, Purdue University, United States

Copyright © 2018 Hassan, Jamal, Radusky, Tiwari, Ullah, Ali, Behramand, de Carvalho, Shams, Khan, Figueiredo, Barh, Ghosh, Silva, Baumbach, Röttger, Turjanski and Azevedo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Vasco A. C. Azevedo, vasco@icb.ufmg.br Adrián G. Turjanski, adrian@qi.fcen.uba.ar

These authors have contributed equally to this work.