Improving performance of mammalian microRNA target prediction

Background  

MicroRNAs (miRNAs) are single-stranded non-coding RNAs known to regulate a wide range of cellular processes by silencing the gene expression at the protein and/or mRNA levels. Computational prediction of miRNA targets is essential for elucidating the detailed functions of miRNA. However, the prediction specificity and sensitivity of the existing algorithms are still poor to generate meaningful, workable hypotheses for subsequent experimental testing. Constructing a richer and more reliable training data set and developing an algorithm that properly exploits this data set would be the key to improve the performance current prediction algorithms.     

Analysis and prediction of carbohydrate binding sites

An analysis of the characteristic properties of sugar binding sites was performed on a set of 19 sugar binding proteins. For each site six parameters were evaluated: solvation potential, residue propensity, hydrophobicity, planarity, protrusion and relative accessible surface area. Three of the parameters were found to distinguish the observed sugar binding sites from the other surface patches. These parameters were then used to calculate the probability for a surface patch to be a carbohydrate binding site. The prediction was optimized on a set of 19 non-homologous carbohydrate binding structures and a test prediction was carried out on a set of 40 protein-carbohydrate complexes. The overall accuracy of prediction achieved was 65%. Results were in general better for carbohydrate-binding enzymes than for the lectins, with a rate of success of 87%.     

ProPred: prediction of HLA-DR binding sites.

ProPred is a graphical web tool for predicting MHC class II binding regions in antigenic protein sequences. The server implement matrix based prediction algorithm, employing amino-acid/position coefficient table deduced from literature. The predicted binders can be visualized either as peaks in graphical interface or as colored residues in HTML interface. This server might be a useful tool in locating the promiscuous binding regions that can bind to several HLA-DR alleles. AVAILABILITY: The server is available at http://www.imtech.res.in/raghava/propred/ CONTACT: raghava@imtech.res.in SUPPLEMENTARY INFORMATION: http://www.imtech.res.in/raghava/propred/page3.html     

PSSM-based prediction of DNA binding sites in proteins

Background  

Detection of DNA-binding sites in proteins is of enormous interest for technologies targeting gene regulation and manipulation. We have previously shown that a residue and its sequence neighbor information can be used to predict DNA-binding candidates in a protein sequence. This sequence-based prediction method is applicable even if no sequence homology with a previously known DNA-binding protein is observed. Here we implement a neural network based algorithm to utilize evolutionary information of amino acid sequences in terms of their position specific scoring matrices (PSSMs) for a better prediction of DNA-binding sites.     

miRTar Hunter: A prediction system for identifying human microRNA target sites

MicroRNAs (miRNAs) are important regulators of gene expression and play crucial roles in many biological processes including apoptosis, differentiation, development, and tumorigenesis. Recent estimates suggest that more than 50% of human protein coding genes may be regulated by miRNAs and that each miRNA may bind to 300–400 target genes. Approximately 1,000 human miRNAs have been identified so far with each having up to hundreds of unique target mRNAs. However, the targets for a majority of these miRNAs have not been identified due to the lack of large-scale experimental detection techniques. Experimental detection of miRNA target sites is a costly and time-consuming process, even though identification of miRNA targets is critical to unraveling their functions in various biological processes. To identify miRNA targets, we developed miRTar Hunter, a novel computational approach for predicting target sites regardless of the presence or absence of a seed match or evolutionary sequence conservation. Our approach is based on a dynamic programming algorithm that incorporates more sequence-specific features and reflects the properties of various types of target sites that determine diverse aspects of complementarities between miRNAs and their targets. We evaluated the performance of our algorithm on 532 known human miRNA:target pairs and 59 experimentally-verified negative miRNA:target pairs, and also compared our method with three popular programs for 481 miRNA:target pairs. miRTar Hunter outperformed three popular existing algorithms in terms of recall and precision, indicating that our unique scheme to quantify the determinants of complementary sites is effective at detecting miRNA targets. miRTar Hunter is now available at http://203.230.194.162/~kbkim.     

miRNA及其靶位点多态性的研究进展

微RNA(microRNA,miRNA)是一类进化上保守、长度为21～23nt的非编码单链小RNA,参与个体发育、器官形成、细胞增殖、分化和细胞凋亡等生物学过程,并在其中发挥重要的调节作用。近年来研究发现,miRNA及其靶位点的多态将引起不同类型的疾患。文章主要从miRNA及其靶位点的多态类型,以及由多态性引起的相关疾病等方面来阐述miRNA的最新进展。     

CaMELS: In silico prediction of calmodulin binding proteins and their binding sites

Due to Ca²⁺‐dependent binding and the sequence diversity of Calmodulin (CaM) binding proteins, identifying CaM interactions and binding sites in the wet‐lab is tedious and costly. Therefore, computational methods for this purpose are crucial to the design of such wet‐lab experiments. We present an algorithm suite called CaMELS (CalModulin intEraction Learning System) for predicting proteins that interact with CaM as well as their binding sites using sequence information alone. CaMELS offers state of the art accuracy for both CaM interaction and binding site prediction and can aid biologists in studying CaM binding proteins. For CaM interaction prediction, CaMELS uses protein sequence features coupled with a large‐margin classifier. CaMELS models the binding site prediction problem using multiple instance machine learning with a custom optimization algorithm which allows more effective learning over imprecisely annotated CaM‐binding sites during training. CaMELS has been extensively benchmarked using a variety of data sets, mutagenic studies, proteome‐wide Gene Ontology enrichment analyses and protein structures. Our experiments indicate that CaMELS outperforms simple motif‐based search and other existing methods for interaction and binding site prediction. We have also found that the whole sequence of a protein, rather than just its binding site, is important for predicting its interaction with CaM. Using the machine learning model in CaMELS, we have identified important features of protein sequences for CaM interaction prediction as well as characteristic amino acid sub‐sequences and their relative position for identifying CaM binding sites. Python code for training and evaluating CaMELS together with a webserver implementation is available at the URL: http://faculty.pieas.edu.pk/fayyaz/software.html#camels .     

ProPred1: prediction of promiscuous MHC Class-I binding sites

SUMMARY: ProPred1 is an on-line web tool for the prediction of peptide binding to MHC class-I alleles. This is a matrix-based method that allows the prediction of MHC binding sites in an antigenic sequence for 47 MHC class-I alleles. The server represents MHC binding regions within an antigenic sequence in user-friendly formats. These formats assist user in the identification of promiscuous MHC binders in an antigen sequence that can bind to large number of alleles. ProPred1 also allows the prediction of the standard proteasome and immunoproteasome cleavage sites in an antigenic sequence. This server allows identification of MHC binders, who have the cleavage site at the C terminus. The simultaneous prediction of MHC binders and proteasome cleavage sites in an antigenic sequence leads to the identification of potential T-cell epitopes. AVAILABILITY: Server is available at http://www.imtech.res.in/raghava/propred1/. Mirror site of this server is available at http://bioinformatics.uams.edu/mirror/propred1/ Supplementary information: Matrices and document on server are available at http://www.imtech.res.in/raghava/propred1/page2.html     

Genome-wide de novo prediction of cis-regulatory binding sites in prokaryotes

Although cis-regulatory binding sites (CRBSs) are at least as important as the coding sequences in a genome, our general understanding of them in most sequenced genomes is very limited due to the lack of efficient and accurate experimental and computational methods for their characterization, which has largely hindered our understanding of many important biological processes. In this article, we describe a novel algorithm for genome-wide de novo prediction of CRBSs with high accuracy. We designed our algorithm to circumvent three identified difficulties for CRBS prediction using comparative genomics principles based on a new method for the selection of reference genomes, a new metric for measuring the similarity of CRBSs, and a new graph clustering procedure. When operon structures are correctly predicted, our algorithm can predict 81% of known individual binding sites belonging to 94% of known cis-regulatory motifs in the Escherichia coli K12 genome, while achieving high prediction specificity. Our algorithm has also achieved similar prediction accuracy in the Bacillus subtilis genome, suggesting that it is very robust, and thus can be applied to any other sequenced prokaryotic genome. When compared with the prior state-of-the-art algorithms, our algorithm outperforms them in both prediction sensitivity and specificity.     

Genome-wide identification of protein binding sites on RNAs in mammalian cells

RNA-binding proteins (RBPs) are proteins that bind to the RNA and participate in forming ribonucleoprotein complexes. They have crucial roles in various biological processes such as RNA splicing, editing, transport, maintenance, degradation, intracellular localization and translation. The RBPs bind RNA with different RNA-sequence specificities and affinities, thus, identification of protein binding sites on RNAs (R-PBSs) will deeper our understanding of RNA-protein interactions. Currently, high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP, also known as CLIP-Seq) is one of the most powerful methods to map RNA-protein binding sites or RNA modification sites. However, this method is only used for identification of single known RBPs and antibodies for RBPs are required. Here we developed a novel method, called capture of protein binding sites on RNAs (RPBS-Cap) to identify genome-wide protein binding sites on RNAs without using antibodies. Double click strategy is used for the RPBS-Cap assay. Proteins and RNAs are UV-crosslinked in vivo first, then the proteins are crosslinked to the magnetic beads. The RNA elements associated with proteins are captured, reverse transcribed and sequenced. Our approach has potential applications for studying genome-wide RNA-protein interactions.     

Functional linkage of adenine nucleotide binding sites in mammalian muscle 6-phosphofructokinase

6-Phosphofructokinases (Pfk) are homo- and heterooligomeric, allosteric enzymes that catalyze one of the rate-limiting steps of the glycolysis: the phosphorylation of fructose 6-phosphate at position 1. Pfk activity is modulated by a number of regulators including adenine nucleotides. Recent crystal structures from eukaryotic Pfk revealed several adenine nucleotide binding sites. Herein, we determined the functional relevance of two adenine nucleotide binding sites through site-directed mutagenesis and enzyme kinetic studies. Subsequent characterization of Pfk mutants allowed the identification of the activating (AMP, ADP) and inhibitory (ATP, ADP) allosteric binding sites. Mutation of one binding site reciprocally influenced the allosteric regulation through nucleotides interacting with the other binding site. Such reciprocal linkage between the activating and inhibitory binding sites is in agreement with current models of allosteric enzyme regulation. Because the allosteric nucleotide binding sites in eukaryotic Pfk did not evolve from prokaryotic ancestors, reciprocal linkage of functionally opposed allosteric binding sites must have developed independently in prokaryotic and eukaryotic Pfk (convergent evolution).