Prediction of solution properties and dynamics of RNAs by means of Brownian dynamics simulation of coarse-grained models: Ribosomal 5S RNA and phenylalanine transfer RNA

Abstract

BackgroundThe possibility of validating biological macromolecules with locally disordered domains like RNA against solution properties is helpful to understand their function. In this work, we present a computational scheme for predicting global properties and mimicking the internal dynamics of RNA molecules in solution. A simple coarse-grained model with one bead per nucleotide and two types of intra-molecular interactions (elastic interactions and excluded volume interactions) is used to represent the RNA chain. The elastic interactions are modeled by a set of Hooke springs that form a minimalist elastic network. The Brownian dynamics technique is employed to simulate the time evolution of the RNA conformations.ResultsThat scheme is applied to the 5S ribosomal RNA of E. Coli and the yeast phenylalanine transfer RNA. From the Brownian trajectory, several solution properties (radius of gyration, translational diffusion coefficient, and a rotational relaxation time) are calculated. For the case of yeast phenylalanine transfer RNA, the time evolution and the probability distribution of the inter-arm angle is also computed.ConclusionsThe general good agreement between our results and some experimental data indicates that the model is able to capture the tertiary structure of RNA in solution. Our simulation results also compare quite well with other numerical data. An advantage of the scheme described here is the possibility of visualizing the real time macromolecular dynamics.Electronic supplementary materialThe online version of this article (doi:10.1186/s13628-015-0025-7) contains supplementary material, which is available to authorized users.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

“…Finally, we note the work of Benı́tez et al, who used an ENM parameterization to describe the diffusion properties of RNA structures (tRNA and rRNA). In this case, the elastic networks were simulated using a Brownian dynamics , approach that introduces hydrodynamic interactions and thus realistic solvent friction.…”

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties

“…For DNA structures, four different lengths of duplexes (8-24 basepair), for which hydrodynamic properties have been published (33)(34)(35)(36) were used for code validation. RNA molecules in isolation are flexible (37,38), so we used two different protein-RNA complexes-histidinyl-tRNA synthetase (39,40) and bacterial 50s ribosomal subunit-to test the method on known structures containing RNA that are expected to be less flexible. These procedures resulted in three final data sets of folded structures to validate the convex hull method: one data set for which translational hydrodynamic properties of proteins have been experimentally determined (Table 1), one data set for which rotational properties of proteins have been experimentally determined (Table 2), and one data set of DNA duplexes and RNA-protein complexes for which several types of experimental data are available (Table 3).…”

Dynamic light scattering: a practical guide and applications in biomedical sciences

“…To validate SAXS models without any bias, we calculate hydrodynamic parameters including R h and sedimentation coefficient from the models using program HYDROPRO (Ortega et al 2011b), and compare them to experimentally determined properties by DLS and complementary hydrodynamic techniques. It is also possible to study DNA and polysaccharides using packages in the HYDRO family Benitez et al 2015). We have followed such practice for proteins, RNA, DNA, and their complexes (Ariyo et al 2015;Deo et al 2014Deo et al , 2015Dzananovic et al 2013Dzananovic et al , 2014Meier et al 2013;Patel et al 2010Patel et al , 2011Patel et al , 2012Patel et al , 2014Vadlamani et al 2015).…”