During the past several years, the possibilities of selecting targets using computational approaches with integrated "omics" data such as genomics, proteomics, metabolomics etc have been increasing continuously. Amongst these, two in silico methods comparative genomics and subtractive genomics are being widely used for the prediction and identification of potential drug targets in numerous pathogenic bacteria and fungal species (Sakharkar et al. 2004; Perumal et al. 2007; Amineni et al. 2010; Abadio et al. 2011). In principle, these approaches rely on search for such genes which are absent in host but are present in the pathogen. Furthermore, these non-host homologues must be essential for the survival of the pathogen and a critical component in vital physicochemical and metabolic pathways, so that a designed drug or a lead compound specific to such target(s) will only impact on the pathogen's system, without hampering host physiology or any aspect of host biology. A great piece of aid has been provided to this area with the availability of complete genome sequences of several pathogenic microorganisms. Such steps aim to reduce the problem of searching for potential drug targets from a large list to selecting from a chosen few. Anti-bacterials are essentially inhibitors of certain bacterial enzymes, all enzymes specific to bacteria can be considered as potential drug targets (Galperin and Koonin 1999).
Get your grade
or your money back
using our Essay Writing Service!
The Bacillus cereus group comprises six species, Bacillus cereus, Bacillus thuringiensis, Bacillus anthracis, Bacillus weihenstephanensis, Bacillus mycoides and Bacillus pseudomycoides (Klee et al. 2010). B.anthracis is a Gram-positive, aerobic, spore-forming bacteria and is the causative agent behind a bacterial disease known as Anthrax. Although primarily a disease of animals, it can also infect human beings with fatal consequences (Baillie 2009; Keim et al. 2009). Isolates of B. anthracis forms a highly monophyletic clade, and are differentiated on the basis of single nucleotide polymorphisms (SNPs) and variable number of tandem repeats (VNTRs). The pathogen is able to cause edema and cell death by a tripartite toxin consisting of the protective antigen, the edema factor, and the lethal factor. The production of a poly-D-glutamic acid capsule allows it to escape the immune system (Keim et al. 2000; Mock and Mignot 2003; Van Ert et al. 2007).
The intense focus on the Ames strain during the anthrax letter attack investigation lead to a highly accurate sequence determination that may be completely free of errors. Therefore, it represents best publicly available genome sequence of B.anthracis. The DNA was generated from the frozen stock used for strain distribution to many different laboratories, worldwide. Hence, it is an important reference genome for research purposes. Overall, the ancestral Ames genome is 5,503,926 nucleotides in size with 5,775 protein coding genes identified. In addition, there are 33 ribosomal RNA genes (23S, 16S and 5S) arranged in 11 operons and, along with 95 tRNA genes, found exclusively on the chromosome. The chromosome itself represents about 95% of the genome with the two large plasmids containing the remaining coding capacity (pXO1=181,677 bp; pXO2=94,830 bp) (Read et al. 2003; Ravel et al. 2009).
It is anticipated that the identified common drug targets will expand our understanding of the molecular mechanisms of B.anthracis pathogenesis and facilitate the identification of novel drug candidates.
Results and discussion
Infectious diseases are the second leading cause of death worldwide . Though there is an increasing demand for new antimicrobial agents, their development lacks the much required progress due to requirement of huge investment, less market, short term usage in same patient, and high level of competition with newly developed agents (Spellberg et al. 2004). Developments in Bioinformatics have brought the algorithms, tools, and facilitated the automation of microbial genome sequencing, development of integrated databases over the internet, comparison of genomes, identification of gene product function, and paved the way for development of antimicrobial agents, vaccines, and rational drug design (Bansal 2005).
Here we report comparative genomics analysis of different metabolic pathways in B.anthracis. Total 90 genes were identified as non-homologous to human genome. DEG v6.0 is a database which contains experimentally determined essential genes from....... different bacteria. Furthermore, these 90 enzymes were subjected to BLAST against DEG v.6.0. onsidered bit score >100 and E.value < e^-10 to enhance the specificity of enzyme in B.anthracis. At present DEG do not contains information about essential genes of B.anthracis however, it contains essential genes data for B.subtilis which belongs to same genra... Following selection criteria Total 8 enzymes out of all were found to be essential for B.anthracis life cycle (Table 2). These targets were found to be potential targets and could be considered for rational drug design. Using metabolic pathway information as the starting point for the identification of potential targets has its advantages as each step in the pathway is validated as the essential function for the survival of the bacterium. The sub cellular localization analysis of all supposed essential and unique enzymes of B.anthracis were evaluated by cello server.
Always on Time
Marked to Standard
As it was suggested that membrane associated protein could be the better target for developing vaccines. After analysis 5 proteins were found to be located in the cytoplasm, 1 as trans membrane and 2 on plasma membrane. A result of comparative analysis of the metabolic pathways of the host and pathogen by using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database reveals some pathways that are unique to B.anthracis and are not present in the humans like D-Alanine metabolism, Peptidoglycan biosynthesis, Phosphotransferase system (PTS), bacterial secretion system and two-component system. In addition to these pathways other pathways like streptomycin biosynthesis, novobiocin biosynthesis,C5-Branched dibasic acid metabolism and many other pathways are not present in humans, and some enzymes of these pathways we found non homologous to human protein sequence, but they are not essential for survival of the B.anthracis. Non homologous sequences from these pathways are 93 and are present in table 1. Further we analyze our essential enzymes of DEG database result against the drug bank database, and we identified about 8 approved drug target and 24 small molecule. Note that we did not consider the hypothetical or putative protein sequences as drug target.
Unique metabolic pathways
All those pathways which are present in pathogen but absent in the host were considered as unique pathways. Based on our manual comparative analysis 28 pathways were identified as unique pathways. Major of them are the two component system (TCS), phosphotransferase system (PTS), bacterial secretion system, flageller assembly system. Since these pathways are not present in H.sapiens, it can be postulated that genes in these pathways are also unique to the pathogen and thus can serve as potential drug targets. However, not every enzyme or protein is always a favourable drug target if for example; it has low accessibility value and/or is not essential for the survival of the pathogen. Therefore, there is need to identify genes whose products regulate factors such as essential nutrients uptake, virulence and pathogenicity. For this purpose, we compared proteins from unique pathways against DEG. Although, presently DEG does not contain information about experimentally proven essential genes of B.anthracis however, information from other Gram-positive and Gram-negative bacteria and one member of genus i.e. B.subtilis is present. Comparative analysis against these bacteria led us to identify several essential genes of B.anthracis as key components of metabolic pathways. We also compared our results with experimental information taken from the literature wherever possible. General mechanism and potential drug targets identified during this study are discussed below
Phosphotransferase system (PTS)
Firstly discovered in E.coli (Kundig et al. 1964), the phosphoenolpyruvate (PEP) dependent phosphotransferase system (PTS) is now well established as a dynamic transport system among diverse range of bacterial species and is responsible for the uptake and phosphorylation of numerous carbohydrates. The basic composition of the PTS is similar in different bacterial species and consists of two general cytoplasmic components, enzyme 1 (E1) and histidine phosphorcarrier protein (HPr). Membrane bound sugar specific permeases (EII) form the third component of this system and show specificity towards carbohydrates. Therefore, different EIIs are present in bacteria. The PTS uses PEP as a source of energy and phosphoryl donor. Each EII complex consists of one or two integral membrane domains (domains C and D) which are hydrophobic in nature and two hydrophilic domains (domains A and B), which together are responsible for the transport of the carbohydrate across the bacterial membrane as well as its phosphorylation. In a sense, the EII complexes constitute parallel transport pathways connected to a common PEP/EI/HPr phosphoryl transfer pathway. E. coli contains at least 15 different EII complexes. A similar number of PTSs have also been reported in B. subtilis (Reizer et al. 1999; Lindner et al. 2002). A growing body of evidence suggests that the pathogens depend upon their hosts for the uptake of nutrients. Therefore, PTS regulation network not only controls carbohydrate uptake and metabolism but also interferes with the utilization of nitrogen and phosphorus and the virulence of certain pathogens (Galperin and Koonin 1999; Deutscher et al. 2006). These findings favor the targeting of PTS proteins as potential drug targets. During our in silico analysis, we identified 10 genes as non homologues to humans from the PTS of B.anthracis. These proteins were then checked against DEG database and two proteins, phosphoenolpyruvate-protein phosphotransferase (ptsI) and glucose-specific II ABC component (ptsG) (also present in amino sugar and nucleotide sugar metabolism) were found to be essential for B.anthracis.
This Essay is
a Student's Work
This essay has been submitted by a student. This is not an example of the work written by our professional essay writers.Examples of our work
Cell wall is one of the major components shared by all bacterial species. Presence of cell wall helps bacteria maintain their morphology as well as to withstand against unfavourable conditions such as in a nonisotonic environment. Typical bacterial cell wall is composed of several components such as peptidoglycan (PG), teichoic acids and proteins (Schleifer KH (1983). Peptidogylcan, which forms more than 70% of the weight of the cell wall, is a large molecule responsible for maintaining the morphology and balance in osmotic pressure. Biosynthesis of PG is a complex two steps process of assembly and polymerization. At the first step, assembly of the peptide moiety of the monomer unit takes place by the successive additions of L-alanine, D-glutamic acid, meso-diaminopimelic acid or L-lysine, and D-alanyl-D-alanine to UDP-N-acetylmuramic acid (UDP-MurNAc). These steps are catalyzed by specific peptide synthetases (ligases), which are designated as MurC, MurD, MurE, and MurF, respectively. The second stage of the PG biosynthesis takes place in the periplasmic space and is catalyzed by the penicillin-binding proteins. This involves transglycosylation and transpeptidation reactions of the disaccharide pentapeptide monomers.
Antibiotics such as Î²-lactams and glycopeptides either inhibit or are poor substrates for the final transpeptidation step of cell wall biosynthesis, resulting in a weakening of the cell wall, followed by rapid lysis and death. These drugs have been a popular mainstay in the treatment of bacterial infections at least in part due to this powerful bactericidal effect. However, there has been a significant increase in the incidence of resistance to these drug classes among important bacterial pathogens. Even so, cell wall biosynthesis remains a valid target for novel antibiotic development, especially for agents that can specifically inhibit any one of the series of essential enzymatic functions involved in assembly of the peptidoglycan.
MurA enzyme which catalyzes has been shown to be inhibited by fosfomycin. However, one potential drawback of MurA as a novel antibiotic target is the presence of two separate genes, murA1 and murA2, that encode for proteins with the same enzymatic activity in gram positive pathogens such as Staphylococcus aureus and Streptococcus pneumoniae. Mutagenic studies have showed that disruption of either murA1 or murA2 had no significant affects on cell growth but cells were unable to survive when both genes were removed (Du et al. 2000). Differences in the active site also make it difficult to develop MurA-specific antibiotic that could effectively inhibit both enzymes because the two homologs of murA share less than 60% identity among gram-positive species. Like other gram-positive bacteria, B.anthracis also contains two murA gene homologs and show high degree of intrinsic fosfomycin resistance. Therefore, there is a need to identify additional drug targets which could effectively block essential metabolic pathways and their respective genes. In our analysis, we compared 17 number of genes from PG biosynthesis pathway of B.anthrax against human genome. A total of 11 enzymes were identified as non homologues to humans. All these enzymes were then further classified following selection criteria (see materials and methods). Subcellular localization were determined for each non homologue and many of the identified proteins
It was identified that UDP-Mur- NAc:L-alanine ligase (MurC) is non homologous to human proteins. Furthermore, MurC is also present in D-glutamine and D-glutamate bacterial metabolic pathways. This makes The MurC, encoded by the murC gene, an interesting candidate for drug development against B.anthrax. Other ligases such as MurD, MurE and MurF from B.anthracis showed homology to.... . Although members of Mur family share less homology among them but they are more conserved when compared among different species of bacteria. We therefore performed sequence and structural comparison of MurC ligase with Tuber and Coryn.
Inhibition of cell wall synthesis serves as a primary antibiotic target for Gram-positive and Gram-negative bacteria and will continue to provide new targets for drug development. Therefore,
Two component systems
Since their appearance on earth, bacteria have evolved a variety of functions to respond to environmental changes. One such specific mechanism is the evolution of two-component signal transduction systems (TCS). A typical TCS is composed of a sensor kinase (histidine kinase, HK), which is capable of auto-phosphorylation in response to an environmental signal, and a response regulator (RR) that interacts with the phosphorylated HK (Gao and Stock 2009). A bacterium possesses multiple TCSs to respond appropriately to different environmental changes such as pH, nutrient level, redox state, osmotic pressure, quorum signals, and antibiotics. Some TCSs also control gene clusters that contribute to cell growth, virulence, biofilms, quorum sensing, etc. During the infectious cycle, pathogenic bacteria encounter different microenvironments and their ability to efficiently adapt towards their host organisms is frequently mediated by TCSs, which can, therefore, be considered as an essential prerequisite for their pathogenicity (Beier and Gross 2006). In addition, bacteria also exchange information between different TCSs to form a complex signal transduction network with greater sensitivity (Eguchi and Utsumi 2008; Mitrophanov and Groisman 2008). For example, P.aeruginosa, which inhabits diverse environments, is estimated to have 64 HKs and 72 RRs. Among the TCSs, 19 TCSs are involved in some way with virulence or antibiotic resistance. Even for Gram-positive S. pneumoniae, 10 of its 13 TCSs are involved in pathogenicity (Gooderham and Hancock 2009).
Since signal transduction in mammals occurs by a different mechanism, antimicrobial agents with potential to target enzymes of TCS have attracted much attention in the past years. However, evolution of multi drug resistant microbes is also on the rise therefore, there is need to identify additional drug targets in TCSs. Identification of such targets will be more useful if they found to be non homologue to the host proteins, essential for pathogen survival and common among diverse range of pathogenic bacteria. In case of B.anthracis, we identified 37 TCS proteins which showed no homology to humans. Subsequent sub cellular localization predictions and DEG analysis was performed for each non homologous protein (Table..). Three proteins were identified as essential for the pathogen namely, respiratory nitrate reductase alpha chain narG (also involved in nitrogen metabolism), respiratory nitrate reductase beta chain narH (also involved in nitrogen metabolism), and chromosomal replication initiation protein dnaA.
Bacterial secretion systems
Bacteria form a very wide diversity of biotic associations, ranging from biofilms to mutualistic or pathogenic associations with larger host organisms. Protein secretion pathways plays a key role in modulating all of these interactions including all the known folding and targeting routes of inner and outer membrane proteins as well as of proteins that are secreted by several specific export routes. In Gram-negative bacteria, six general classes of protein secretion systems have been identified. Each class shows considerable diversity and facilitate the process of entry of the secreted proteins inside host cells and modification of host physiology thus promoting colonization. On the other hand, in Gram-positive bacteria, secreted proteins are commonly translocated across the single membrane by the Sec pathway or the two-arginine (Tat) pathway. Although, Gram-positive bacteria share some of the same secretion systems as Gram-negative bacteria, others such as mycobacteria that have a hydrophobic, nearly impermeable cell wall, called the mycomembrane, a specialized type VII secretion system translocates proteins across both the membrane and the cell wall (Sandkvist 2001; Tseng et al. 2009). Therefore, the importance of secretion systems not only for bacterial viability but also for pathogenicity is well established and has a great potential as a target for development of new drugs (Briken 2008). KEGG database contains information about..... proteins of the B.anthracis secretion system. Each protein sequence was subjected to BLAST analysis and four were identified as non homologues to humans. Amongst them three were predicted as membrane localized and one into the cytoplasm. DEG based search for essential proteins suggested....... proteins as essential for B.anthracis namely, .........
In this study, we have performed a comparative metabolic pathway analysis of the host H. sapiens and the pathogen B. anthracis. The computational genomic approach has already facilitated the search for potential drug targets against many pathogens. Use of the DEG database is more efficient than conventional methods for identification of essential genes and facilitates the exploratory identification of the most relevant drug targets in the pathogen. In the list of DEG result enzymes are involved in peptidoglycan biosynthesis, phosphotransferase system(PTS), bacterial secretion, two component system and other many pathways The present study has thus led to the identification of several proteins that can be targeted for effective drug design and vaccine development against B.anthracis. The drugs developed against identified targets will be specific to the pathogen and may show lesser toxicity to the host. Although the number of essential genes in the metabolic pathways of B.anthracis, identified in the present study, is relatively small (only 8), these can be further characterized and their role in the survival of the bacteria can be verified. Since many of them have been reported to play a role in its virulence, a systematic approach to develop novel drugs against these targets can be adapted for treating B.anthracis infections. Further, homology modeling of these targets will help identify the best possible sites that can be targeted for drug design by simulation modeling.
Materials and methods
Identification of host and pathogen metabolic pathways
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used as a source (Kanehisa et al. 2006; Kanehisa et al. 2010) of metabolic pathway information. List of metabolic pathways and pathway identification numbers of the host H.sapiens and the pathogen B.anthracis was extracted from the KEGG database and saved locally. A manual comparative analysis was then performed and pathways which do not appear in the host but present in the pathogen according to KEGG database annotation were selected as pathways unique to B.anthracis. Genes in these unique pathways were identified and the respective protein sequences were obtained from Uniprot. Protein sequences were initially filtered on the basis of number of amino acids residues i.e. proteins of less than 100 amino acid residues were excluded.
Identification of potential drug targets
Two step comparison was performed between host and pathogen genomes for the identification of non homologues proteins. At first, only proteins from pathogen specific pathways were subjected to a BLASTP. Secondly, proteins from common pathways were compared by BLASTP. In each scenario, search was restricted to proteins from H.sapeins only, e-value inclusion threshold set to 0.005 and minimum bit score of 100. Proteins, which do not have hits below the e-value inclusion threshold of 0.005, were picked out as potential drug targets.
The selected proteins from the unique pathways were then BLASTP against DEG to identify essential genes. After pBLAST analysis the 235 non-homologous enzymes were further analyzed for essentiality to pathogen by DEG database (Zhang et al. 2004). Essential genes are those indispensable for the survival of an organism, and their functions are therefore, considered a foundation of life (Zhang et al. 2004).
59 enzymes out of 235 from different pathways were found to be essential. These imperative 59 enzymes were submitted in Drug Bank database against approve drug targets and small molecules legends.
Evaluating biological significance of the drug targets
Biological significance and subcellular locations of the non homologues genes were calculated by two methods; This is required to find out the surface membrane proteins which could be probable vaccine targets. The predicted membrane proteins were further analyzed in psortb v3.0  to confirm if the drug targets are identified as membrane proteins irrespective of the subcellular localization prediction methods
AMB designed the study and drafted the manuscript. ST performed computational analysis and generated results. MI and YT analyzed the data and finalized the manuscript.
We are thankful to the Director of the National Centre of Excellence in Molecular Biology (CEMB) for providing facilities.
Conflicts of Interest
Figure 1 - Sample figure title
Figure legend text
Figure 2 - Another sample figure title
Figure legend text.
Table 1 - Sample table title
Table legend text.
Table 2 - Another sample table title
Table legend text.