Using reliability information to annotate RNA secondary structures.

M Zuker; A B Jacobson

doi:10.1017/s1355838298980116

Home/Archive/PubMed

PubMed · 9622126

Using reliability information to annotate RNA secondary structures.

M Zuker·A B Jacobson

View original source ↗Export metadata ↓

Abstract

A number of heuristic descriptors have been developed previously in conjunction with the mfold package that describe the propensity of individual bases to participate in base pairs and whether or not a predicted helix is "well-determined." They were developed for the "energy dot plot" output of mfold. Two descriptors, P-num and H-num, are used to measure the level of promiscuity in the association of any given nucleotide or helix with alternative complementary pairs. The third descriptor, S-num, measures the propensity of bases to be single-stranded. In the current work, we describe a series of programs that were developed in order to annotate individual structures with "well-definedness" information. We use color annotation to present the information. The programs can annotate PostScript files that are created by the mfold package or the PostScript secondary structure plots produced by the Weiser and Noller program XRNA (Weiser B, Noller HF, 1995, XRNA: Auto-interactive program for modeling RNA, The Center for Molecular Biology of RNA, Santa Cruz, California: University of California; Internet: ftp://fangio.ucsc.edu/pub/XRNA). In addition, these programs can annotate ss files that serve as input to XRNA. The annotation package can also handle structure comparison with a reference structure. This feature can be used to compare predicted structure with a phylogenetically deduced model, to compare two different predicted foldings, and to identify conformational changes that are predicted between wild-type and mutant RNAs. We provide several examples of application. Predicted structures of two RNase P RNAs were colored with P-num information and further annotated with comparative information. The comparative model of a 16S rRNA was annotated with P-num information from mfold and with base pair probabilities obtained from the Vienna RNA folding package. Further annotation adds comparisons with the optimal foldings obtained from mfold and the Vienna package, respectively. The results of all of these analyses are discussed in the context of the reliability of structure prediction.

Explore related subjects

Computer Simulation Endoribonucleases Models, Molecular Nucleic Acid Conformation RNA RNA, Bacterial RNA, Catalytic RNA, Protozoan RNA, Ribosomal, 16S Reproducibility of Results Ribonuclease P Software Journal Article Research Support, U.S. Gov't, Non-P.H.S.Research Support, U.S. Gov't, P.H.S.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

M Zuker, A B Jacobson. 1998. Using reliability information to annotate RNA secondary structures.. https://doi.org/10.1017/s1355838298980116

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

KEEP EXPLORING

Related citations

PubMed2025

Generating correlated data for omics simulation.

Simulation of realistic omics data is a key input for benchmarking studies that help users obtain optimal computational pipelines. Omics data involves large numbers of measured features on each sample and these measures are generally correlated with each other. However, simulation too often ignores these correlations, perhaps due to computational and statistical hurdles of doing so. To alleviate this, we describe three approaches for generating omics-scale data with correlated measures which mimic real datasets. These approaches are all based on a Gaussian copula approach with a covariance matrix that decomposes into a diagonal part and a low-rank part. This decomposition allows for extremely efficient simulation, overcoming a hurdle for adoption of past methods. We use these approaches to demonstrate the importance of including correlation in two benchmarking applications. First, we show that variance of results from the popular DESeq2 method increases when dependence is included. Second, we demonstrate that CYCLOPS, a method for inferring circadian time of collection from transcriptomics, improves in performance when given gene-gene dependencies in some circumstances. We provide an R package, dependentsimr, that has efficient implementations of these methods and can generate dependent data with arbitrary marginal distributions, including discrete (binary, ordered categorical, Poisson, negative binomial), continuous (normal), or with an empirical distribution.

Computer Simulation↗

PubMed2017

Addressing current challenges in cancer immunotherapy with mathematical and computational modelling.

The goal of cancer immunotherapy is to boost a patient's immune response to a tumour. Yet, the design of an effective immunotherapy is complicated by various factors, including a potentially immunosuppressive tumour microenvironment, immune-modulating effects of conventional treatments and therapy-related toxicities. These complexities can be incorporated into mathematical and computational models of cancer immunotherapy that can then be used to aid in rational therapy design. In this review, we survey modelling approaches under the umbrella of the major challenges facing immunotherapy development, which encompass tumour classification, optimal treatment scheduling and combination therapy design. Although overlapping, each challenge has presented unique opportunities for modellers to make contributions using analytical and numerical analysis of model outcomes, as well as optimization algorithms. We discuss several examples of models that have grown in complexity as more biological information has become available, showcasing how model development is a dynamic process interlinked with the rapid advances in tumour-immune biology. We conclude the review with recommendations for modellers both with respect to methodology and biological direction that might help keep modellers at the forefront of cancer immunotherapy development.

Computer Simulation↗

PubMed2002

Structure of docosahexaenoic acid-containing phospholipid bilayers as studied by (2)H NMR and molecular dynamics simulations.

Polyunsaturated phospholipids are known to be important with regard to the biological functions of essential fatty acids, for example, involving neural tissues such as the brain and retina. Here we have employed two complementary structural methods for the study of polyunsaturated bilayer lipids, viz. deuterium ((2)H) NMR spectroscopy and molecular dynamics (MD) computer simulations. Our research constitutes one of the first applications of all-atom MD simulations to polyunsaturated lipids containing docosahexaenoic acid (DHA; 22:6 cis-Delta(4,7,10,13,16,19)). Structural features of the highly unsaturated, mixed-chain phospholipid, 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (PDPC), have been studied in the liquid-crystalline (L(alpha)) state and compared to the less unsaturated homolog, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). The (2)H NMR spectra of polyunsaturated bilayers are dramatically different from those of less unsaturated phospholipid bilayers. We show how use of MD simulations can aid in interpreting the complex (2)H NMR spectra of polyunsaturated bilayers, in conjunction with electron density profiles determined from small-angle X-ray diffraction studies. This work clearly demonstrates preferred helical and angle-iron conformations of the polyunsaturated chains in liquid-crystalline bilayers, which favor chain extension while maintaining bilayer flexibility. The presence of relatively long, extended fatty acyl chains may be important for solvating the hydrophobic surfaces of integral membrane proteins, such as rhodopsin. In addition, the polyallylic DHA chains have a tendency to adopt back-bended (hairpin-like) structures, which increase the interfacial area per lipid. Finally, the material properties have been analyzed in terms of the response of the bilayer to mechanical stress. Simulated bilayers of phospholipids containing docosahexaenoic acid were less sensitive to the applied surface tension than were saturated phospholipids, possibly implying a decrease in membrane elasticity (area elastic modulus, bending rigidity). The above features distinguish DHA-containing lipids from saturated or monounsaturated lipids and may be important for their biological modes of action.

Computer Simulation↗