Sven Jager1, Benjamin Schiller2, Philipp Babel1, Malte Blumenroth1, Thorsten Strufe2, Kay Hamacher3. 1. Department of Biology, TU Darmstadt, Schnittspahnstr. 2, 64283 Darmstadt, Germany. 2. Department of Computer Science, TU Dresden, Nöthnitzer Str. 46, 01187 Dresden, Germany. 3. Department of Biology, Department of Computer Science, Department of Physics, TU Darmstadt, Schnittspahnstr. 2, 64283 Darmstadt, Germany.
Abstract
BACKGROUND: In this work, we present a new coarse grained representation of RNA dynamics. It is based on adjacency matrices and their interactions patterns obtained from molecular dynamics simulations. RNA molecules are well-suited for this representation due to their composition which is mainly modular and assessable by the secondary structure alone. These interactions can be represented as adjacency matrices of k nucleotides. Based on those, we define transitions between states as changes in the adjacency matrices which form Markovian dynamics. The intense computational demand for deriving the transition probability matrices prompted us to develop StreAM-[Formula: see text], a stream-based algorithm for generating such Markov models of k-vertex adjacency matrices representing the RNA. RESULTS: We benchmark StreAM-[Formula: see text] (a) for random and RNA unit sphere dynamic graphs (b) for the robustness of our method against different parameters. Moreover, we address a riboswitch design problem by applying StreAM-[Formula: see text] on six long term molecular dynamics simulation of a synthetic tetracycline dependent riboswitch (500 ns) in combination with five different antibiotics. CONCLUSIONS: The proposed algorithm performs well on large simulated as well as real world dynamic graphs. Additionally, StreAM-[Formula: see text] provides insights into nucleotide based RNA dynamics in comparison to conventional metrics like the root-mean square fluctuation. In the light of experimental data our results show important design opportunities for the riboswitch.
BACKGROUND: In this work, we present a new coarse grained representation of RNA dynamics. It is based on adjacency matrices and their interactions patterns obtained from molecular dynamics simulations. RNA molecules are well-suited for this representation due to their composition which is mainly modular and assessable by the secondary structure alone. These interactions can be represented as adjacency matrices of k nucleotides. Based on those, we define transitions between states as changes in the adjacency matrices which form Markovian dynamics. The intense computational demand for deriving the transition probability matrices prompted us to develop StreAM-[Formula: see text], a stream-based algorithm for generating such Markov models of k-vertex adjacency matrices representing the RNA. RESULTS: We benchmark StreAM-[Formula: see text] (a) for random and RNA unit sphere dynamic graphs (b) for the robustness of our method against different parameters. Moreover, we address a riboswitch design problem by applying StreAM-[Formula: see text] on six long term molecular dynamics simulation of a synthetictetracycline dependent riboswitch (500 ns) in combination with five different antibiotics. CONCLUSIONS: The proposed algorithm performs well on large simulated as well as real world dynamic graphs. Additionally, StreAM-[Formula: see text] provides insights into nucleotide based RNA dynamics in comparison to conventional metrics like the root-mean square fluctuation. In the light of experimental data our results show important design opportunities for the riboswitch.