BarMap models RNA folding on changing energy landscapes by linking macrostates between landscape snapshots, enabling efficient analysis of co-transcriptional and externally perturbed folding scenarios.
Published as BarMap: RNA Folding on Dynamic Energy Landscapes, RNA 2010.
BarMap computes RNA folding kinetics on energy landscapes that change over time, for example while the molecule is being transcribed, by mapping the barrier trees of successive landscapes onto one another and propagating the population from one to the next.
RNA folding rarely happens on a fixed energy landscape. During transcription, degradation, ligand binding, or temperature shifts, the available structure space changes over time, and with it the kinetic pathways. This paper addresses exactly that nonstationary setting by extending landscape-based RNA kinetics from a single static landscape to a sequence of related landscapes connected through time.
The central idea of BarMap is to treat each time point or perturbation step as a snapshot landscape with its own macrostates, and then define a mapping between neighboring snapshots. Population densities can then be transferred from one coarse-grained landscape to the next instead of starting the kinetics calculation from scratch each time. In effect, the expensive landscape analysis becomes a preprocessing step, while the temporal evolution is handled by moving populations across linked barrier-tree representations.
Co-transcriptional folding is the obvious case: the chain grows by one nucleotide at a time and each elongation step reshapes the landscape slightly. The same formalism applies to temperature changes, refolding after cleavage, and mechanically constrained scenarios. In this way the energy-landscape description of RNA folding kinetics is extended to dynamic landscapes, i.e. to the situations regulatory RNAs encounter in cells.
From a methodological perspective, this is a useful complement to trajectory-based simulators. It does not aim to reconstruct every folding path explicitly. Instead, it keeps the coarse-grained view at the level of basins, barriers, and saddle points, while allowing those basins themselves to change over time. That makes it possible to study larger and more complex nonstationary folding scenarios than a direct simulation strategy would usually permit.
BarMap is therefore one of the key papers in the co-transcriptional RNA folding, RNA folding kinetics, and energy landscapes cluster. It formalizes the idea that kinetics should often be thought of as motion across a sequence of changing landscapes, not just as diffusion on a single fixed one.
Abstract of the paper
Dynamical changes of RNA secondary structures play an important role in the function of many regulatory RNAs. Such kinetic effects, especially in time-variable and externally triggered systems, are usually investigated by means of extensive and expensive simulations of large sets of individual folding trajectories. Here we describe the theoretical foundations of a generic approach that not only allows the direct computation of approximate population densities but also reduces the efforts required to analyze the folding energy landscapes to a one-time preprocessing step. The basic idea is to consider the kinetics on individual landscapes and to model external triggers and environmental changes as small but discrete changes in the landscapes. A ‘‘barmap’’ links macrostates of temporally adjacent landscapes and defines the transfer of population densities from one ‘‘snapshot’’ to the next. Implemented in the BarMap software, this approach makes it feasible to study folding processes at the level of basins, saddle points, and barriers for many nonstationary scenarios, including temperature changes, cotranscriptional folding, refolding in consequence to degradation, and mechanically constrained kinetics, as in the case of the translocation of a polymer through a pore.