Sharmila Anishetty

Brief communication: Potential drug targets in Mycobacterium tuberculosis through metabolic pathway analysis

The emergence of multidrug resistant varieties of Mycobacterium tuberculosis has led to a search for novel drug targets. We have performed an insilico comparative analysis of metabolic pathways of the host Homo sapiens and the pathogen M. tuberculosis. Enzymes from the biochemical pathways of M. tuberculosis from the KEGG metabolic pathway database were compared with proteins from the host H. sapiens, by performing a BLASTp search against the non-redundant database restricted to the H. sapiens subset. The e-value threshold cutoff was set to 0.005. Enzymes, which do not show similarity to any of the host proteins, below this threshold, were filtered out as potential drug targets. We have identified six pathways unique to the pathogen M. tuberculosis when compared to the host H. sapiens. Potential drug targets from these pathways could be useful for the discovery of broad spectrum drugs. Potential drug targets were also identified from pathways related to lipid metabolism, carbohydrate metabolism, amino acid metabolism, energy metabolism, vitamin and cofactor biosynthetic pathways and nucleotide metabolism. Of the 185 distinct targets identified from these pathways, many are in various stages of progress at the TB Structural Genomics Consortium. However, 67 of our targets are new and can be considered for rational drug design. As a case study, we have built a homology model of one of the potential drug targets MurD ligase using WHAT IF software. The model could be further explored for insilico docking studies with suitable inhibitors. The study was successful in listing out potential drug targets from the M. tuberculosis proteome involved in vital aspects of the pathogens metabolism, persistence, virulence and cell wall biosynthesis. This systematic evaluation of metabolic pathways of host and pathogen through reliable and conventional bioinformatic methods can be extended to other pathogens of clinical interest.

Thalidomide attenuates nitric oxide-driven angiogenesis by interacting with soluble guanylyl cyclase

Background and purpose:  Nitric oxide (NO) promotes angiogenesis by activating endothelial cells. Thalidomide arrests angiogenesis by interacting with the NO pathway, but its putative targets are not known. Here, we have attempted to identify these targets.

Tripeptide analysis of protein structures

BackgroundAn efficient building block for protein structure prediction can be tripeptides. 8000 different tripeptides from a dataset of 1220 high resolution (≤ 2.0°A) structures from the Protein Data Bank (PDB) have been looked at, to determine which are structurally rigid and non-rigid. This data has been statistically analyzed, discussed and summarized. The entire data can be utilized for the building of protein structures.ResultsTripeptides have been classified into three categories: rigid, non-rigid and intermediate, based on the relative structural rigidity between Cα and Cβ atoms in a tripeptide. We found that 18% of the tripeptides in the dataset can be classified as rigid, 4% as non-rigid and 78% as intermediate. Many rigid tripeptides are made of hydrophobic residues, however, there are tripeptides with polar side chains forming rigid structures. The bulk of the tripeptides fall in the intermediate class while very small numbers actually fall in the non-rigid class. Structurally all rigid tripeptides essentially form two structural classes while the intermediate and non-rigid tripeptides fall into one structural class. This notion of rigidity and non-rigidity is designed to capture side chain interactions but not secondary structures.ConclusionsRigid tripeptides have no correlation with the secondary structures in proteins and hence this work is complementary to such studies. Tripeptide data may be used to predict plausible structures for oligopeptides and for denovo protein design.

In silico identification of putative metal binding motifs

Metal ion binding domains are found in proteins that mediate transport, buffering or detoxification of metal ions. In this study, we have performed an in silico analysis of metal binding proteins and have identified putative metal binding motifs for the ions of cadmium, cobalt, zinc, arsenic, mercury, magnesium, manganese, molybdenum and nickel. A pattern search against the UniProtKB/Swiss-Prot and UniProtKB/TrEMBL databases yielded true positives in each case showing the high-specificity of the motifs. Motifs were also validated against PDB structures and site directed mutagenesis studies.

Molecular dynamics simulations of human and dog gastric lipases: Insights into domain movements

Mammalian gastric lipases are stable and active under acidic conditions and also in the duodenal lumen. There has been considerable interest in acid stable lipases owing to their potential application in the treatment of pancreatic exocrine insufficiency. In order to gain insights into the domain movements of these enzymes, molecular dynamics simulations of human gastric lipase was performed at an acidic pH and under neutral conditions. For comparative studies, simulation of dog gastric lipase was also performed at an acidic pH. Analyses show, that in addition to the lid region, there is another region of high mobility in these lipases. The potential role of this novel region is discussed.

Functionally specified protein signatures distinctive for each of the different blue copper proteins.

BackgroundProteins having similar functions from different sources can be identified by the occurrence in their sequences, a conserved cluster of amino acids referred to as pattern, motif, signature or fingerprint. The wide usage of protein sequence analysis in par with the growth of databases signifies the importance of using patterns or signatures to retrieve out related sequences. Blue copper proteins are found in the electron transport chain of prokaryotes and eukaryotes. The signatures already existing in the databases like the type 1 copper blue, multiple copper oxidase, cyt b/b6, photosystem 1 psaA&B, psaG&K, and reiske iron sulphur protein are not specified signatures for blue copper proteins as the name itself suggests. Most profile and motif databases strive to classify protein sequences into a broad spectrum of protein families. This work describes the signatures designed based on the copper metal binding motifs in blue copper proteins. The common feature in all blue copper proteins is a trigonal planar arrangement of two nitrogen ligands [each from histidine] and one sulphur containing thiolate ligand [from cysteine], with strong interactions between the copper center and these ligands.ResultsSequences that share such conserved motifs are crucial to the structure or function of the protein and this could provide a signature of family membership. The blue copper proteins chosen for the study were plantacyanin, plastocyanin, cucumber basic protein, stellacyanin, dicyanin, umecyanin, uclacyanin, cusacyanin, rusticyanin, sulfocyanin, halocyanin, azurin, pseudoazurin, amicyanin and nitrite reductase which were identified in both eukaryotes and prokaryotes. ClustalW analysis of the protein sequences of each of the blue copper proteins was the basis for designing protein signatures or peptides. The protein signatures and peptides identified in this study were designed involving the active site region involving the amino acids bound to the copper atom. It was highly specific for each kind of blue copper protein and the false picks were minimized. The set of signatures designed specifically for the BCPs was entirely different from the existing broad spectrum signatures as mentioned in the background section.ConclusionsThese signatures can be very useful for the annotation of uncharacterized proteins and highly specific to retrieve blue copper protein sequences of interest from the non redundant databases containing a large deposition of protein sequences.

In silico analysis of DosR regulon proteins of Mycobacterium tuberculosis

One of the challenges faced by Mycobacterium tuberculosis (M. tuberculosis) in dormancy is hypoxia. DosR/DevR of M. tuberculosis is a two component dormancy survival response regulator which induces the expression of 48 genes. In this study, we have used DosR regulon proteins of M. tuberculosis H37Rv as the query set and performed a comprehensive homology search against the non-redundant database. Homologs were found in environmental mycobacteria, environmental bacteria and archaebacteria. Analysis of genomic context of DosR regulon revealed that they are distributed as nine blocks in the genome of M. tuberculosis with many transposases and integrases in their vicinity. Further, we classified DosR regulon proteins into eight functional categories. One of the hypothetical proteins Rv1998c could probably be a methylisocitrate lyase or a phosphonomutase. Another hypothetical protein, Rv0572 was found only in mycobacteria. Insights gained in this study can potentially aid in the development of novel therapeutic interventions.

Comparative pan genome analysis of oral Prevotella species implicated in periodontitis

Prevotella is part of the oral bacterial community implicated in periodontitis. Pan genome analyses of eight oral Prevotella species, P. dentalis, P. enoeca, P. fusca, P. melaninogenica, P. denticola, P. intermedia 17, P. intermedia 17-2 and P. sp. oral taxon 299 are presented in this study. Analysis of the Prevotella pan genome revealed features such as secretion systems, resistance to oxidative stress and clustered regularly interspaced short palindromic repeat (CRISPR)-Cas systems that enable the bacteria to adapt to the oral environment. We identified the presence of type VI secretion system (T6SS) in P. fusca and P. intermedia strains. For some VgrG and Hcp proteins which were not part of the core T6SS loci, we used gene neighborhood analysis and identified putative effector proteins and putative polyimmunity loci in P. fusca and polymorphic toxin systems in P. intermedia strains. Earlier studies have identified the presence of Por secretion system (PorSS) in P. gingivalis, P. melaninogenica and P. intermedia. We noted the presence of their homologs in six other oral Prevotella studied here. We suggest that in Prevotella, PorSS is used to secrete cysteine proteases such as interpain and C-terminal domain containing proteins with a “Por_secre_tail” domain. We identified subtype I-B CRISPR-Cas system in P. enoeca. Putative CRISPR-Cas system subtypes for 37 oral Prevotella and 30 non-oral Prevotella species were also predicted. Further, we performed a BLASTp search of the Prevotella proteins which are also conserved in the red-complex pathogens, against the human proteome to identify potential broad-spectrum drug targets. In summary, the use of a pan genome approach enabled identification of secretion systems and defense mechanisms in Prevotella that confer adaptation to the oral cavity.

Molecular dynamics simulations of lectin domain of FimH and immunoinformatics for the design of potential vaccine candidates

Adhesion of uropathogenic E. coli (UPEC) to uroepithelial cell receptors is facilitated through the lectin domain of FimH adhesin. In the current study, Molecular Dynamics (MD) simulations were performed for the lectin domain of FimH from UPEC J96. The high affinity state lectin domain was found to be stable and rigid during the simulations. Further, based on conserved subsequences around one of the disulfide forming cysteines, two sequence motifs were designed. An immunoinformatics approach was utilized to identify linear and discontinuous epitopes for the lectin domain of FimH. We propose that the accessibility of predicted epitopes should also be assessed in a dynamic aqueous environment to evaluate the potential of vaccine candidates. Since MD simulation data enables assessing the accessibility in a dynamic environment, we evaluated the accessibility of the top ranked discontinuous and linear epitopes using structures obtained at every nanosecond (ns) in the 1-20 ns MD simulation timeframe. Knowledge gained in this study has a potential utility in the design of vaccine candidates for Urinary Tract Infection (UTI).

Understanding mutations and protein stability through tripeptides.

A novel methodology to predict the local conformational changes in a protein as a consequence of missense mutations is proposed. A pentapeptide at the locus of mutation plays the dominant role and it is analyzed in terms of tripeptides. A measure for spatial and temporal fluctuations in a pentapeptide is devised and validated. The method does not involve any prior knowledge of structural templates from sequence homology studies. Structural deformations can be predicted with about 70–80% reliability in any protein. Disease causing mutations and benign mutations have been addressed. In particular, p53, retinoblastoma protein and lipoprotein lipase are studied in detail.