RNA secondary structural alignment with conditional random fields.

Kengo Sato; Yasubumi Sakakibara

doi:10.1093/bioinformatics/bti1139

Home/Archive/PubMed

PubMed · 16204111

RNA secondary structural alignment with conditional random fields.

Kengo Sato·Yasubumi Sakakibara

View original source ↗Export metadata ↓

Abstract

MOTIVATION: The computational identification of non-coding RNA regions on the genome is currently receiving much attention. However, it is essentially harder than gene-finding problems for protein-coding regions because non-coding RNA sequences do not have strong statistical signals. Since comparative sequence analysis is effective for non-coding RNA detection, efficient computational methods are expected for structural alignment of RNA sequences. Several methods have been proposed to accomplish the structural alignment tasks for RNA sequences, and we found that one of the most important points is to estimate an accurate score matrix for calculating structural alignments. RESULTS: We propose a novel approach for RNA structural alignment based on conditional random fields (CRFs). Our approach has some specific features compared with previous methods in the sense that the parameters for structural alignment are estimated such that the model can most probably discriminate between correct alignments and incorrect alignments, and has the generalization ability so that a satisfiable score matrix can be obtained even with a small number of sample data without overfitting. Experimental results clearly show that the parameter estimation with CRFs can outperform all the other existing methods for structural alignments of RNA sequences. Furthermore, structural alignment search based on CRFs is more accurate for predicting non-coding RNA regions than the other scoring methods. These experimental results strongly support our discriminative method employing CRFs to estimate the score matrix parameters. AVAILABILITY: The program which is implemented in C++ is available at http://phmmts.dna.bio.keio.ac.jp/ under the GNU public license.

Explore related subjects

Base Sequence Chromosome Mapping Computer Simulation Effect Modifier, Epidemiologic Models, Genetic Models, Statistical Molecular Sequence Data Nucleic Acid Conformation RNA, Untranslated Sequence Alignment Sequence Analysis, RNA Sequence Homology, Nucleic Acid Journal Article Research Support, Non-U.S. Gov't

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Kengo Sato, Yasubumi Sakakibara. 2005-09-01. RNA secondary structural alignment with conditional random fields.. https://doi.org/10.1093/bioinformatics%2Fbti1139

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2016

miR-503-3p promotes epithelial-mesenchymal transition in breast cancer by directly targeting SMAD2 and E-cadherin.

Although progress in clinical and basic research has significantly increased our understanding of breast cancer, little is known about the molecular mechanism underlying breast cancer metastasis. Identification of effective therapeutic targets to prevent breast cancer metastasis is urgently needed. The function of miR-503-3p has been investigated in other cancers, but its role in breast cancer remains undefined. Here, we found that miR-503-3p was overexpressed in breast cancer tissue and plasma compared with adjacent normal breast tissue and with plasma from healthy individuals. Moreover, we identified miR-503-3p to be an oncogene of breast cancer cell proliferation, migration and invasion. Upregulation of miR-503-3p in breast cancer cells inhibited expression of epithelial-mesenchymal transition (EMT)-related protein SMAD2 and the epithelial marker protein E-cadherin by directly binding to their mRNA 3' untranslated region, whereas increased expression of mesenchymal marker proteins, including vimentin and N-cadherin. Taken together, our findings support a critical role for miR-503-3p in induction of breast cancer EMT and suggest that plasma miR-503-3p may be a useful diagnostic biomarker for breast cancer.

Base Sequence↗

PubMed2007

Highly sensitive detection of cytotoxicity using a modified HSP70B' promoter.

We have previously found that the DNA fragment from nucleotides (nts) -287 to +110 in the HSP70B' gene is a functional promoter responding to Cadmium Chloride-induced cytotoxicity (Wada et al., Biotechnol Bioeng, 92, 410-415, 2005). In order to increase the cytotoxic response of this promoter, we first determined the location of the cytotoxic responding element (CRE) and then constructed tandem repeats of the CRE in front of the HSP70B' promoter. 5'- and 3'-deletion analysis revealed that the DNA fragment from nts -192 to -56 in the HSP70B' gene induces a significant response to cytotoxicity. When the AP-1 binding site in this region was mutated, the basal activity of HSP70B' gene promoter decreased but the cytotoxic response was unchanged. Thus, the CRE is located in nts -192 to -56 in the HSP70B' promoter, and the AP-1 binding site is not essential for the cytotoxic response. In addition, cells transfected with a luciferase construct carrying three tandem repeats of the CRE upstream of the HSP70B' promoter and containing AP-1 binding site mutation, showed a 2.28-fold higher response than that of no repeats. Moreover, the detection limit of Cadmium Chloride in the cells was 382 pmol/mL. Thus, highly sensitive sensor cells for Cadmium Chloride can be constructed using a HSP70B' promoter construct containing upstream tandem repeats of the CRE and mutation of the AP-1 binding site.

Base Sequence↗

PubMed2007

In vitro selection and characterization of RNA aptamers binding thyroxine hormone.

RNA possesses the ability to bind a wide repertoire of small molecules. Some of these binding interactions have been shown to be of primary importance in molecular biology. For example, several classes of mRNA domains, collectively referred to as riboswitches, have been shown to serve as RNA genetic control elements that sense the concentrations of specific metabolites (i.e. acting as direct sensors of chemical compounds). However, to date no RNA species binding a hormone has been reported. Here, we report that the use of an appropriate SELEX (systematic evolution of ligands by exponential enrichment) strategy results in the isolation of thyroxine-specific aptamers. Further biochemical characterization of these aptamers, including mutational studies, the use of transcripts with site-specific modified nucleotides, nuclease and chemical probing, binding-shift assays and CD, demonstrated that these RNA structures included a G-rich motif, reminiscent of a guanine quadruplex structure, adjacent to a helical region. The presence of the thyroxine appeared to be essential for the formation of the structural motif's scaffold. Moreover, the binding is shown to be specific to thyroxine (T4) and tri-iodothyronine (T3), the active forms of the hormone, whereas other inactive derivatives, including thyronine (T0), do not support complex formation. These results suggest that this aptamer specifically binds to the iodine moieties of the thyroxine, a previously unreported ability for an RNA molecule.

Base Sequence↗