106,694 research outputs found
A new procedure to analyze RNA non-branching structures
RNA structure prediction and structural motifs analysis are challenging tasks in the investigation of RNA function. We propose a novel procedure to detect structural motifs shared between two RNAs (a reference and a target). In particular, we developed two core modules: (i) nbRSSP_extractor, to assign a unique structure to the reference RNA encoded by a set of non-branching structures; (ii) SSD_finder, to detect structural motifs that the target RNA shares with the reference, by means of a new score function that rewards the relative distance of the target non-branching structures compared to the reference ones. We integrated these algorithms with already existing software to reach a coherent pipeline able to perform the following two main tasks: prediction of RNA structures (integration of RNALfold and nbRSSP_extractor) and search for chains of matches (integration of Structator and SSD_finder)
Ab initio RNA folding
RNA molecules are essential cellular machines performing a wide variety of
functions for which a specific three-dimensional structure is required. Over
the last several years, experimental determination of RNA structures through
X-ray crystallography and NMR seems to have reached a plateau in the number of
structures resolved each year, but as more and more RNA sequences are being
discovered, need for structure prediction tools to complement experimental data
is strong. Theoretical approaches to RNA folding have been developed since the
late nineties when the first algorithms for secondary structure prediction
appeared. Over the last 10 years a number of prediction methods for 3D
structures have been developed, first based on bioinformatics and data-mining,
and more recently based on a coarse-grained physical representation of the
systems. In this review we are going to present the challenges of RNA structure
prediction and the main ideas behind bioinformatic approaches and physics-based
approaches. We will focus on the description of the more recent physics-based
phenomenological models and on how they are built to include the specificity of
the interactions of RNA bases, whose role is critical in folding. Through
examples from different models, we will point out the strengths of
physics-based approaches, which are able not only to predict equilibrium
structures, but also to investigate dynamical and thermodynamical behavior, and
the open challenges to include more key interactions ruling RNA folding.Comment: 28 pages, 18 figure
Evolutionary Dynamics and Optimization: Neutral Networks as Model-Landscapes for RNA Secondary-Structure Folding-Landscapes
We view the folding of RNA-sequences as a map that assigns a pattern of base
pairings to each sequence, known as secondary structure. These preimages can be
constructed as random graphs (i.e. the neutral networks associated to the
structure ). By interpreting the secondary structure as biological
information we can formulate the so called Error Threshold of Shapes as an
extension of Eigen's et al. concept of an error threshold in the single peak
landscape. Analogue to the approach of Derrida & Peliti for a of the population
on the neutral network. On the one hand this model of a single shape landscape
allows the derivation of analytical results, on the other hand the concept
gives rise to study various scenarios by means of simulations, e.g. the
interaction of two different networks. It turns out that the intersection of
two sets of compatible sequences (with respect to the pair of secondary
structures) plays a key role in the search for ''fitter'' secondary structures.Comment: 20 pages, uuencoded compressed postscript-file, Proc. of ECAL '95
conference, to appear., email: chris @ imb-jena.d
Chance and Necessity in Evolution: Lessons from RNA
The relationship between sequences and secondary structures or shapes in RNA
exhibits robust statistical properties summarized by three notions: (1) the
notion of a typical shape (that among all sequences of fixed length certain
shapes are realized much more frequently than others), (2) the notion of shape
space covering (that all typical shapes are realized in a small neighborhood of
any random sequence), and (3) the notion of a neutral network (that sequences
folding into the same typical shape form networks that percolate through
sequence space). Neutral networks loosen the requirements on the mutation rate
for selection to remain effective. The original (genotypic) error threshold has
to be reformulated in terms of a phenotypic error threshold. With regard to
adaptation, neutrality has two seemingly contradictory effects: It acts as a
buffer against mutations ensuring that a phenotype is preserved. Yet it is
deeply enabling, because it permits evolutionary change to occur by allowing
the sequence context to vary silently until a single point mutation can become
phenotypically consequential. Neutrality also influences predictability of
adaptive trajectories in seemingly contradictory ways. On the one hand it
increases the uncertainty of their genotypic trace. At the same time neutrality
structures the access from one shape to another, thereby inducing a topology
among RNA shapes which permits a distinction between continuous and
discontinuous shape transformations. To the extent that adaptive trajectories
must undergo such transformations, their phenotypic trace becomes more
predictable.Comment: 37 pages, 14 figures; 1998 CNLS conference; high quality figures at
http://www.santafe.edu/~walte
A Statistical Analysis of RNA Folding Algorithms Through Thermodynamic Parameter Perturbation
Computational RNA secondary structure prediction is rather well established.
However, such prediction algorithms always depend on a large number of
experimentally measured parameters. Here, we study how sensitive structure
prediction algorithms are to changes in these parameters. We find that already
for changes corresponding to the actual experimental error to which these
parameters have been determined 30% of the structure are falsly predicted and
the ground state structure is preserved under parameter perturbation in only 5%
of all cases. We establish that base pairing probabilities calculated in a
thermal ensemble are a viable though not perfect measure for the reliability of
the prediction of individual structure elements. A new measure of stability
using parameter perturbation is proposed, and its limitations discussed.Comment: 6 pages, 3 figures, 1 table submitted to Nucleic Acids Researc
Probing complex RNA structures by mechanical force
RNA secondary structures of increasing complexity are probed combining single
molecule stretching experiments and stochastic unfolding/refolding simulations.
We find that force-induced unfolding pathways cannot usually be interpretated
by solely invoking successive openings of native helices. Indeed, typical
force-extension responses of complex RNA molecules are largely shaped by
stretching-induced, long-lived intermediates including non-native helices. This
is first shown for a set of generic structural motifs found in larger RNA
structures, and then for Escherichia coli's 1540-base long 16S ribosomal RNA,
which exhibits a surprisingly well-structured and reproducible unfolding
pathway under mechanical stretching. Using out-of-equilibrium stochastic
simulations, we demonstrate that these experimental results reflect the slow
relaxation of RNA structural rearrangements. Hence, micromanipulations of
single RNA molecules probe both their native structures and long-lived
intermediates, so-called "kinetic traps", thereby capturing -at the single
molecular level- the hallmark of RNA folding/unfolding dynamics.Comment: 9 pages, 9 figure
Reconstructing phylogeny from RNA secondary structure via simulated evolution
DNA sequences of genes encoding functional RNA molecules (e.g., ribosomal RNAs) are commonly used in phylogenetics (i.e. to infer evolutionary history). Trees derived from ribosomal RNA (rRNA) sequences, however, are inconsistent with other molecular data in investigations of deep branches in the tree of life. Since much of te functional constraints on the gene products (i.e. RNA molecules) relate to three-dimensional structure, rather than their actual sequences, accumulated mutations in the gene sequences may obscure phylogenetic signal over very large evolutionary time-scales. Variation in structure, however, may be suitable for phylogenetic inference even under extreme sequence divergence. To evaluate qualitatively the manner in which structural evolution relates to sequence change, we simulated the evolution of RNA sequences under various constraints on structural change
Genomic and structural investigation on dolphin morbillivirus (DMV) in Mediterranean fin whales (Balaenoptera physalus).
Dolphin morbillivirus (DMV) has been deemed as one of the most relevant threats for fin whales (Balaenoptera physalus) being responsible for a mortality outbreak in the Mediterranean Sea in the last years. Knowledge of the complete viral genome is essential to understand any structural changes that could modify virus pathogenesis and viral tissue tropism. We report the complete DMV sequence of N, P/V/C, M, F and H genes identified from a fin whale and the comparison of primary to quaternary structure of proteins between this fin whale strain and some of those isolated during the 1990-'92 and the 2006-'08 epidemics. Some relevant substitutions were detected, particularly Asn52Ser located on F protein and Ile21Thr on N protein. Comparing mutations found in the fin whale DMV with those occurring in viral strains of other cetacean species, some of them were proven to be the result of diversifying selection, thus allowing to speculate on their role in host adaptation and on the way they could affect the interaction between the viral attachment and fusion with the target host cells
Pathways and kinetic barriers in mechanical unfolding and refolding of RNA and proteins
Using self-organized polymer models, we predict mechanical unfolding and
refolding pathways of ribo-zymes, and the green fluorescent protein. In
agreement with experiments, there are between six and eight unfolding
transitions in the Tetrahymena ribozyme. Depending on the loading rate, the
number of rips in the force-ramp unfolding of the Azoarcus ribozymes is between
two and four. Force-quench refolding of the P4-P6 subdomain of the Tetrahymena
ribozyme occurs through a compact intermediate. Subsequent formation of
tertiary contacts between helices P5b-P6a and P5a/P5c-P4 leads to the native
state. The force-quench refolding pathways agree with ensemble experiments. In
the dominant unfolding route, the N-terminal a helix of GFP unravels first,
followed by disruption of the N terminus b strand. There is a third
intermediate that involves disruption of three other strands. In accord with
experiments, the force-quench refolding pathway of GFP is hierarchic, with the
rate-limiting step being the closure of the barrel.Comment: 33 pages 7 figure
- …