5,257 research outputs found
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
RNAiFold2T: Constraint Programming design of thermo-IRES switches
Motivation: RNA thermometers (RNATs) are cis-regulatory ele- ments that
change secondary structure upon temperature shift. Often involved in the
regulation of heat shock, cold shock and virulence genes, RNATs constitute an
interesting potential resource in synthetic biology, where engineered RNATs
could prove to be useful tools in biosensors and conditional gene regulation.
Results: Solving the 2-temperature inverse folding problem is critical for RNAT
engineering. Here we introduce RNAiFold2T, the first Constraint Programming
(CP) and Large Neighborhood Search (LNS) algorithms to solve this problem.
Benchmarking tests of RNAiFold2T against existent programs (adaptive walk and
genetic algorithm) inverse folding show that our software generates two orders
of magnitude more solutions, thus allow- ing ample exploration of the space of
solutions. Subsequently, solutions can be prioritized by computing various
measures, including probability of target structure in the ensemble, melting
temperature, etc. Using this strategy, we rationally designed two thermosensor
internal ribosome entry site (thermo-IRES) elements, whose normalized
cap-independent transla- tion efficiency is approximately 50% greater at 42?C
than 30?C, when tested in reticulocyte lysates. Translation efficiency is lower
than that of the wild-type IRES element, which on the other hand is fully
resistant to temperature shift-up. This appears to be the first purely
computational design of functional RNA thermoswitches, and certainly the first
purely computational design of functional thermo-IRES elements. Availability:
RNAiFold2T is publicly available as as part of the new re- lease RNAiFold3.0 at
https://github.com/clotelab/RNAiFold and http:
//bioinformatics.bc.edu/clotelab/RNAiFold, which latter has a web server as
well. The software is written in C++ and uses OR-Tools CP search engine.Comment: 24 pages, 5 figures, Intelligent Systems for Molecular Biology (ISMB
2016), to appear in journal Bioinformatics 201
Limiting the valence: advancements and new perspectives on patchy colloids, soft functionalized nanoparticles and biomolecules
Limited bonding valence, usually accompanied by well-defined directional
interactions and selective bonding mechanisms, is nowadays considered among the
key ingredients to create complex structures with tailored properties: even
though isotropically interacting units already guarantee access to a vast range
of functional materials, anisotropic interactions can provide extra
instructions to steer the assembly of specific architectures. The anisotropy of
effective interactions gives rise to a wealth of self-assembled structures both
in the realm of suitably synthesized nano- and micro-sized building blocks and
in nature, where the isotropy of interactions is often a zero-th order
description of the complicated reality. In this review, we span a vast range of
systems characterized by limited bonding valence, from patchy colloids of new
generation to polymer-based functionalized nanoparticles, DNA-based systems and
proteins, and describe how the interaction patterns of the single building
blocks can be designed to tailor the properties of the target final structures
Explicit factorization of external coordinates in constrained Statistical Mechanics models
If a macromolecule is described by curvilinear coordinates or rigid
constraints are imposed, the equilibrium probability density that must be
sampled in Monte Carlo simulations includes the determinants of different
mass-metric tensors. In this work, we explicitly write the determinant of the
mass-metric tensor G and of the reduced mass-metric tensor g, for any molecule,
general internal coordinates and arbitrary constraints, as a product of two
functions; one depending only on the external coordinates that describe the
overall translation and rotation of the system, and the other only on the
internal coordinates. This work extends previous results in the literature,
proving with full generality that one may integrate out the external
coordinates and perform Monte Carlo simulations in the internal conformational
space of macromolecules. In addition, we give a general mathematical argument
showing that the factorization is a consequence of the symmetries of the metric
tensors involved. Finally, the determinant of the mass-metric tensor G is
computed explicitly in a set of curvilinear coordinates specially well-suited
for general branched molecules.Comment: 22 pages, 2 figures, LaTeX, AMSTeX. v2: Introduccion slightly
extended. Version in arXiv is slightly larger than the published on
RNA Folding with Soft Constraints: Reconciliation of Probing Data and Thermodynamic Secondary Structure Prediction
Thermodynamic folding algorithms and structure probing experiments are commonly used to determine the secondary structure of RNAs. Here we propose a formal framework to reconcile information from both prediction algorithms and probing experiments. The thermodynamic energy parameters are adjusted using ‘pseudo-energies’ to minimize the discrepancy between prediction and experiment. Our framework differs from related approaches that used pseudo-energies in several key aspects. (i) The energy model is only changed when necessary and no adjustments are made if prediction and experiment are consistent. (ii) Pseudo-energies remain biophysically interpretable and hold positional information where experiment and model disagree. (iii) The whole thermodynamic ensemble of structures is considered thus allowing to reconstruct mixtures of suboptimal structures from seemingly contradicting data. (iv) The noise of the energy model and the experimental data is explicitly modeled leading to an intuitive weighting factor through which the problem can be seen as folding with ‘soft’ constraints of different strength. We present an efficient algorithm to iteratively calculate pseudo-energies within this framework and demonstrate how this approach can be used in combination with SHAPE chemical probing data to improve secondary structure prediction. We further demonstrate that the pseudo-energies correlate with biophysical effects that are known to affect RNA folding such as chemical nucleotide modifications and protein binding.Austrian Science Fund. Erwin Schrodinger Fellowship (J2966-B12
Bayesian estimates of free energies from nonequilibrium work data in the presence of instrument noise
The Jarzynski equality and the fluctuation theorem relate equilibrium free
energy differences to non-equilibrium measurements of the work. These relations
extend to single-molecule experiments that have probed the finite-time
thermodynamics of proteins and nucleic acids. The effects of experimental error
and instrument noise have not previously been considered. Here, we present a
Bayesian formalism for estimating free-energy changes from non-equilibrium work
measurements that compensates for instrument noise and combines data from
multiple driving protocols. We reanalyze a recent set of experiments in which a
single RNA hairpin is unfolded and refolded using optical tweezers at three
different rates. Interestingly, the fastest and farthest-from-equilibrium
measurements contain the least instrumental noise, and therefore provide a more
accurate estimate of the free energies than a few slow, more noisy,
near-equilibrium measurements. The methods we propose here will extend the
scope of single-molecule experiments; they can be used in the analysis of data
from measurements with AFM, optical, and magnetic tweezers.Comment: 8 page
Recommended from our members
PATTERNA: transcriptome-wide search for functional RNA elements via structural data signatures.
Establishing a link between RNA structure and function remains a great challenge in RNA biology. The emergence of high-throughput structure profiling experiments is revolutionizing our ability to decipher structure, yet principled approaches for extracting information on structural elements directly from these data sets are lacking. We present PATTERNA, an unsupervised pattern recognition algorithm that rapidly mines RNA structure motifs from profiling data. We demonstrate that PATTERNA detects motifs with an accuracy comparable to commonly used thermodynamic models and highlight its utility in automating data-directed structure modeling from large data sets. PATTERNA is versatile and compatible with diverse profiling techniques and experimental conditions
- …