Hampapathalu A. Nagarajaram

IMEx: Imperfect Microsatellite Extractor

MOTIVATION Microsatellites, also known as simple sequence repeats, are the tandem repeats of nucleotide motifs of the size 1-6 bp found in every genome known so far. Their importance in genomes is well known. Microsatellites are associated with various disease genes, have been used as molecular markers in linkage analysis and DNA fingerprinting studies, and also seem to play an important role in the genome evolution. Therefore, it is of importance to study distribution, enrichment and polymorphism of microsatellites in the genomes of interest. For this, the prerequisite is the availability of a computational tool for extraction of microsatellites (perfect as well as imperfect) and their related information from whole genome sequences. Examination of available tools revealed certain lacunae in them and prompted us to develop a new tool. RESULTS In order to efficiently screen genome sequences for microsatellites (perfect as well as imperfect), we developed a new tool called IMEx (Imperfect Microsatellite Extractor). IMEx uses simple string-matching algorithm with sliding window approach to screen DNA sequences for microsatellites and reports the motif, copy number, genomic location, nearby genes, mutational events and many other features useful for in-depth studies. IMEx is more sensitive, efficient and useful than the available widely used tools. IMEx is available in the form of a stand-alone program as well as in the form of a web-server. AVAILABILITY A World Wide Web server and the stand-alone program are available for free access at http://203.197.254.154/IMEX/ or http://www.cdfd.org.in/imex.

Support Vector Machine-based classification of protein folds using the structural properties of amino acid residues and amino acid residue pairs

MOTIVATION Fold recognition is a key step in the protein structure discovery process, especially when traditional sequence comparison methods fail to yield convincing structural homologies. Although many methods have been developed for protein fold recognition, their accuracies remain low. This can be attributed to insufficient exploitation of fold discriminatory features. RESULTS We have developed a new method for protein fold recognition using structural information of amino acid residues and amino acid residue pairs. Since protein fold recognition can be treated as a protein fold classification problem, we have developed a Support Vector Machine (SVM) based classifier approach that uses secondary structural state and solvent accessibility state frequencies of amino acids and amino acid pairs as feature vectors. Among the individual properties examined secondary structural state frequencies of amino acids gave an overall accuracy of 65.2% for fold discrimination, which is better than the accuracy by any method reported so far in the literature. Combination of secondary structural state frequencies with solvent accessibility state frequencies of amino acids and amino acid pairs further improved the fold discrimination accuracy to more than 70%, which is approximately 8% higher than the best available method. In this study we have also tested, for the first time, an all-together multi-class method known as Crammer and Singer method for protein fold classification. Our studies reveal that the three multi-class classification methods, namely one versus all, one versus one and Crammer and Singer method, yield similar predictions. AVAILABILITY Dataset and stand-alone program are available upon request.

Hansa: an automated method for discriminating disease and neutral human nsSNPs.

Variations are mostly due to nonsynonymous single nucleotide polymorphisms (nsSNPs), some of which are associated with certain diseases. Phenotypic effects of a large number of nsSNPs have not been characterized. Although several methods have been developed to predict the effects of nsSNPs as “disease” or “neutral,” there is still a need for development of methods with improved prediction accuracies. We, therefore, developed a support vector machine (SVM) based method named Hansa which uses a novel set of discriminatory features to classify nsSNPs into disease (pathogenic) and benign (neutral) types. Validation studies on a benchmark dataset and further on an independent dataset of well‐characterized known disease and neutral mutations show that Hansa outperforms the other known methods. For example, fivefold cross‐validation studies using the benchmark HumVar dataset reveal that at the false positive rate (FPR) of 20% Hansa yields a true positive rate (TPR) of 82% that is about 10% higher than the best‐known method. Hansa is available in the form of a web server at http://hansa.cdfd.org.in:8080. Hum Mutat 33:332–337, 2012.

Azadirachtin interacts with the tumor necrosis factor (TNF) binding domain of its receptors and inhibits TNF-induced biological responses.

The role of azadirachtin, an active component of a medicinal plant Neem (Azadirachta indica), on TNF-induced cell signaling in human cell lines was investigated. Azadirachtin blocks TNF-induced activation of nuclear factor κB (NF-κB) and also expression of NF-κB-dependent genes such as adhesion molecules and cyclooxygenase 2. Azadirachtin inhibits the inhibitory subunit of NF-κB (IκBα) phosphorylation and thereby its degradation and RelA (p65) nuclear translocation. It blocks IκBα kinase (IKK) activity ex vivo, but not in vitro. Surprisingly, azadirachtin blocks NF-κB DNA binding activity in transfected cells with TNF receptor-associated factor (TRAF)2, TNF receptor-associated death domain (TRADD), IKK, or p65, but not with TNFR, suggesting its effect is at the TNFR level. Azadirachtin blocks binding of TNF, but not IL-1, IL-4, IL-8, or TNF-related apoptosis-inducing ligand (TRAIL) with its respective receptors. Anti-TNFR antibody or TNF protects azadirachtin-mediated down-regulation of TNFRs. Further, in silico data suggest that azadirachtin strongly binds in the TNF binding site of TNFR. Overall, our data suggest that azadirachtin modulates cell surface TNFRs thereby decreasing TNF-induced biological responses. Thus, azadirachtin exerts an anti-inflammatory response by a novel pathway, which may be beneficial for anti-inflammatory therapy.

MICdb: database of prokaryotic microsatellites

The MICdb (Microsatellites Database) (http://www.cdfd.org.in/micas) is a comprehensive relational database of non-redundant microsatellites extracted from fully sequenced prokaryotic genomes. The current version (1.0) of the database has been compiled from 83 genomes belonging to different phylogenetic groups. This database has been linked to MICAS, the web-based Microstatellite Analysis Server. MICAS provides a user-friendly front-end to systematically extract data on microsatellite tracts from genomes. The database contains the following information pertaining to the microsatellites: the regions (coding/non-coding, if coding, their GenBank annotations) containing microsatellite tracts; the frequencies of their occurrences, the size and the number of repeating motifs; and the sequences of the tracts. MICAS also provides an interface to Autoprimer, a primer design program to automatically design primers for selected microsatellite loci.

Microsatellite polymorphism across the M. tuberculosis and M. bovis genomes: implications on genome evolution and plasticity.

BackgroundMicrosatellites are the tandem repeats of nucleotide motifs of size 1–6 bp observed in all known genomes. These repeats show length polymorphism characterized by either insertion or deletion (indels) of the repeat units, which in and around the coding regions affect transcription and translation of genes.ResultsSystematic comparison of all the equivalent microsatellites in the coding regions of the three mycobacterial genomes, viz. Mycobacterium tuberculosis H37Rv, Mycobacterium tuberculosis CDC1551 and Mycobacterium bovis, revealed for the first time the presence of several polymorphic microsatellites. The coding regions affected by frame-shifts owing to microsatellite indels have undergone changes indicative of gene fission/fusion, premature termination and length variation. Interestingly, the genes affected by frame-shift mutations code for membrane proteins, transporters, PPE, PE_PGRS, cell-wall synthesis proteins and hypothetical proteins.ConclusionThis study has revealed the role of microsatellite indel mutations in imparting novel functions and a certain degree of plasticity to the mycobacterial genomes. There seems to be some correlation between microsatellite polymorphism and the variations in virulence, host-pathogen interactions mediated by surface antigen variations, and adaptation of the pathogens. Several of the polymorphic microsatellites reported in this study can be tested for their polymorphic nature by screening clinical isolates and various mycobacterial strains, for establishing correlations between microsatellite polymorphism and the phenotypic variations among these pathogens.

SilkSatDb: a microsatellite database of the silkworm, Bombyx mori

The SilkSatDb (silkmoth microsatellite database) (http://www.cdfd.org.in/silksatdb) is a relational database of microsatellites extracted from the available expressed sequence tags and whole genome shotgun sequences of the silkmoth, Bombyx mori. The database has been rendered with a simple and robust web-based search facility, developed using PHP. The SilkSatDb also stores information on primers developed and validated in the laboratory. Users can retrieve information on the microsatellite and the protocols used, along with informative figures and polymorphism status of those microsatellites. In addition, the interface is coupled with Autoprimer, a primer-designing program, using which users can design primers for the loci of interest.

Simple sequence repeats in mycobacterial genomes

Simple sequence repeats (SSRs) or microsatellites are the repetitive nucleotide sequences of motifs of length 1–6 bp. They are scattered throughout the genomes of all the known organisms ranging from viruses to eukaryotes. Microsatellites undergo mutations in the form of insertions and deletions (INDELS) of their repeat units with some bias towards insertions that lead to microsatellite tract expansion. Although prokaryotic genomes derive some plasticity due to microsatellite mutations they have in-built mechanisms to arrest undue expansions of microsatellites and one such mechanism is constituted by post-replicative DNA repair enzymes MutL, MutH and MutS. The mycobacterial genomes lack these enzymes and as a null hypothesis one could expect these genomes to harbour many long tracts. It is therefore interesting to analyse the mycobacterial genomes for distribution and abundance of microsatellites tracts and to look for potentially polymorphic microsatellites. Available mycobacterial genomes, Mycobacterium avium, M. leprae, M. bovis and the two strains of M. tuberculosis (CDC1551 and H37Rv) were analysed for frequencies and abundance of SSRs. Our analysis revealed that the SSRs are distributed throughout the mycobacterial genomes at an average of 220–230 SSR tracts per kb. All the mycobacterial genomes contain few regions that are conspicuously denser or poorer in microsatellites compared to their expected genome averages. The genomes distinctly show scarcity of long microsatellites despite the absence of a post-replicative DNA repair system. Such severe scarcity of long microsatellites could arise as a result of strong selection pressures operating against long and unstable sequences although influence of GC-content and role of point mutations in arresting microsatellite expansions can not be ruled out. Nonetheless, the long tracts occasionally found in coding as well as non-coding regions may account for limited genome plasticity in these genomes.

A low prevalence of MYH7/MYBPC3 mutations among Familial Hypertrophic Cardiomyopathy patients in India

Familial Hypertrophic Cardiomyopathy (FHC) is an autosomal dominant disorder affecting the cardiac muscle and exhibits varied clinical symptoms because of genetic heterogeneity. Several disease causing genes have been identified and most code for sarcomere proteins. In the current study, we have carried out clinical and molecular analysis of FHC patients from India. FHC was detected using echocardiography and by analysis of clinical symptoms and family history. Disease causing mutations in the β-cardiac myosin heavy chain (MYH7) and Myosin binding protein C3 (MYBPC3) genes were identified using Polymerase Chain Reaction-Deoxyribose Nucleic Acid (PCR-DNA) sequencing. Of the 55 patient samples screened, mutations were detected in only nineteen in the two genes; MYBPC3 mutations were identified in 12 patients while MYH7 mutations were identified in five, two patients exhibited double heterozygosity. All four MYH7 mutations were missense mutations, whereas only 3/9 MYPBC3 mutations were missense mutations. Four novel mutations in MYBPC3 viz. c.456delC, c.2128G>A (p.E710K), c.3641G>A (p.W1214X), and c.3656T>C (p.L1219P) and one in MYH7 viz. c.965C>T (p.S322F) were identified. A majority of missense mutations affected conserved amino acid residues and were predicted to alter the structure of the corresponding mutant proteins. The study has revealed a greater frequency of occurrence of MYBPC3 mutations when compared to MYH7 mutations.

Effects of disease causing mutations on the essential motions in proteins.

Abstract This study embodies a detailed comparative analysis of the essential motions of the Wild type and the eight different disease mutant forms of the Human CYPlbl. The mutations considered in this study have been implicated in Primary Congenital Glaucoma, an in-born, genetic disorder associated with eye-abnormality. The principal component analysis for Wild type and the Mutants was carried out using the stabilized molecular dynamics trajectories, which ranged from 35 to 45 nanoseconds. Investigations revealed the nature of the collective motions that characterize functionally relevant ‘essential motions’. The essential motions in Wild type are characterized by the collective motions of the Substrate Access Channel including the β-rich domain and the loops in the region of p450-reductase interaction. Comparative analysis of the essential motions of the Wild type and Mutants, especially those involving the functionally important regions indicated distinct differences in their magnitudes as well as the residue-wise distribution. The Mutants in general are associated with higher root mean square fluctuations, and involve some of the relatively intact core regions of the protein, in large collective motions. This study sheds light on the possible effects of disease causing mutations on the large functionally important collective motions in proteins.