# 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

- Aarón Ayllón Benítez
^{1}, - José Ginés Hernández Cifre
^{1}Email author, - Francisco Guillermo Díaz Baños
^{1}and - José García de la Torre
^{1}

**8**:11

https://doi.org/10.1186/s13628-015-0025-7

© Benítez et al. 2015

**Received: **22 April 2015

**Accepted: **18 November 2015

**Published: **1 December 2015

## Abstract

### Background

The 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.

### Results

That 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.

### Conclusions

The 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.

## Keywords

## Background

Ribonucleic acid, or RNA, is a key molecule for the synthesis of proteins in living cells. Transfer RNA (tRNA) carries an amino acid to the protein synthetic machinery of a cell, the ribosome, which is the particle responsible for translation of the genetic code. Ribosomes are mainly formed by the so-called ribosomal RNA (rRNA) along with a variety of protein subunits. As many biological macromolecules, the RNA functionality is related with its structure which, in turn, determines the RNA solution properties. Therefore, the latter are essential sources of information on RNA structure and dynamics [1]. However, the structural determination of even small RNA molecules like tRNAs (70 to 80 nucleotides long) or the 5S rRNA of *E. Coli* (120 nucleotides long [2]) is not an easy task, being most RNAs much larger molecules. For instance, besides 5S rRNA, the main classes of rRNA are 16S and 23S in prokaryotic cells, and 17S, 18S, 26S, and 28S in eukaryotic cells [2] (S stands for Svedberg, the unit of sedimentation coefficient; it is noteworthy that the notation for ribosomal RNAs, as for some other biomacromolecules, is made in terms of a hydrodynamic property).

The tRNA is a particularly simple case as it presents four short stems as sketched in Fig. 1
a. In usual conditions, those stems appear grouped forming two arms joined by a disordered region so that the crystal structure has the typical “L-letter” shape [4]. Thus, tRNA has been presented as a typical case among segmentally flexible macromolecules, composed by rigid parts connected by flexible hinges. Considerable effort has been devoted to the characterization of the structure and degree of flexibility of tRNA in dilute solution [5], by means of either measurements of solution properties [6–10] or atomistic molecular dynamics simulation [11] (particularly noteworthy are the pioneering works of Harvey, McCammon and coworkers [12–15]). Recent works and reviews on tRNA dynamics are those of Alexander et al. [16] and Agirrezabala and Valle [17]. Ribosomal RNAs are structurally quite more complex than tRNAs. For example, the secondary structure of the majority of the 5S rRNAs has five double helices, two hairpin loops, two internal loops, and a hinge region organized into three big domains (*α*, *β*, and *γ*) [18–20] as sketched in Fig. 1
b. Combining results from dilute solution experiments with numerical calculations, some models for the tertiary structure of 5S rRNA have been proposed [19, 20].

The level of detail, or coarse-graining, of a macromolecular model depends on the type of behavior that is intended to capture. Thus, a large variety of coarse-grained models and computational methods are used for biological macromolecules [21]. Since long, our group has developed procedures and devised computational tools to build coarse-grained models and predict solution properties of both rigid [22–24] and flexible macromolecules [25–27]. The latter are represented by bead-and-spring chains where, in the classical concept, each bead corresponds to a large piece of the macromolecular chain. Nonetheless, finer coarse-grained (bead-and-spring) models can be used whose elements represent, for instance, the repeating units – amino acid or nucleotide residues – of the biomacromolecules [28, 29]. Indeed, one of the present authors proposed forty years ago [30] to obtain hydrodynamic properties for double-helical DNA using what is nowadays regarded as a coarse-grained model with one bead per nucleotide and the approximate Kirkwood formula for friction coefficients. The problem was revisited, in further work [28], using modern bead-model hydrodynamics [22] for rigid DNA models. Subsequently, bending and twisting internal flexibility was included in the model by means of stiff, elastic springs [31]. The latter simulation of flexible, double-helical DNA, which is precursor of the present work, is in the spirit of a popular family of coarse-grained models for biomacromolecules and supramacromolecular structures, the so-called “elastic network model” [32], which has been applied to RNA structures [33, 34].

Certainly, the preferred simulation method for short time scales has been, following the above mentioned pioneering works [11, 12, 15], the classical molecular dynamics technique. There are many commercial and public domain tools that implement that technique [21] which has been applied to a variety of RNA structures [35–37]. However, the large size of the macromolecules and the typical long time scales of their global dynamics are not usually accessible to conventional (atomic-level) molecular dynamics simulations. An alternative is to apply the Brownian dynamics simulation technique to some type of coarse-grained model. Thus, in this work, we have attempted to devise a scheme for predicting global solution properties, that would also mimic the internal dynamics, of RNAs using the simple coarse-grained model with one bead per nucleotide, with intra-molecular interactions represented by a minimalist elastic network, along with the Brownian dynamics technique.

## Methods

### Coarse-grained model

The RNA molecules (tRNA and rRNA) were modeled by bead-and-spring chains with one bead per nucleotide. In that model, nucleotides along the RNA strand form a string of touching beads all of them with the same hydrodynamic radius, *σ*=3.15 Å. That value is determined according to the features of the double-helical regions which are modeled following previous works for double-stranded DNA [28, 31], taking into account that RNA double helices appear usually in the A form whereas DNA is commonly found in B (or Z) form. The A form is a right handed double helix with the following features [38]: major groove width 3.8 Å, minor groove width 10.9 Å, helix radius 10 Å, rise 2.8 Å, inclination 16.1°, number of base pairs per helical pitch 11. On the other hand, the phase angle, *ϕ*, which define the gap between the complementary helices and therefore the shape and size of the grooves [28], was set to *ϕ*= 200°.

*i*(within every double-helical region) to:

- 1.its first neighbors along its helix piece (beads
*i*±1), which keeps connectivity and bond equilibrium length. - 2.
its second neighbors along its helix piece (beads

*i*±2), which accounts for bending interactions. - 3.
its third neighbors along its helix piece (beads

*i*±3), which accounts for torsional interactions. - 4.
its counterpart in the complementary helix piece (bead

*i*+*N*_{ bp }), which accounts for interactions between nucleotides forming the base pair. - 5.
the first neighbors of its counterpart in the complementary helix piece (beads

*i*+(*N*_{ bp }±1)), which is necessary in order to keep the double helical structure.

*N*
_{
bp
} is the number of base pairs within a given double-helical region. In the above description, we use the term “helix piece” instead of “strand” because every nucleotide in the RNA molecule actually belongs to a unique strand. Figure 2
a shows the model for a double helix (in the straight, equilibrium conformation) displaying all the connectors involving one of the innermost beads (the rest of connectors was removed for the sake of clarity). On the other hand, each bead in a loop is only connected to their two first neighbors (beads *i*±1).

where *l* is the instantaneous spring length, *l*
_{
e
} is the equilibrium spring length, and *H* is the spring constant. Following previous works [31], we set \(H=200k_{B}T/l_{e,1}^{2}\) for all connectors regardless the type of interaction represented. In that expression, *k*
_{
B
}
*T* is the Boltzmann factor and *l*
_{
e,1} is the equilibrium length of the connector binding bead *i* to its first neighbor (*l*
_{
e,1}=2*σ*=6.30 Å, in our touching-beads model). For the chosen value of *H*, the root-mean-square fluctuation in the spring length, \(\sqrt {\langle l^{2} \rangle - \langle {l_{e}^{2}} \rangle }\), is 5 % of *l*
_{
e
}.

where *r*
_{
ij
} is the instantaneous distance between the interacting beads, *ε*
_{
LJ
}=0.1*k*
_{
B
}
*T*, and *σ*
_{
LJ
}=0.8*l*
_{
e,1} [39].

### Simulation, hydrodynamics, and solution properties

We performed Brownian dynamics (BD) simulations of bead-and-spring models of tRNA and 5S rRNA in order to obtain their solution properties in a solvent with viscosity *η*
_{
s
}=0.01 poise at *T*=293 K. The deterministic forces in the BD algorithm are the derivatives of the potentials defined by Eqs. (1) and (2). Both simulation and analysis were carried out with the software package SIMUFLEX [27] that implements a procedure based on the Ermak-McCammon algorithm [40]. That BD algorithm takes into account the hydrodynamic interaction (HI) effect through the Rotne-Prager-Yamakawa tensor [41]. Although the rigorous inclusion of the HI effect is necessary in order to generate the right macromolecular dynamics, BD simulations without HI (computationally less expensive) can be used as a “smart Monte Carlo” procedure [42] because they sample correctly the equilibrium conformational space. In such an approach, an ensemble of equilibrium conformations is generated and properties are calculated by means of the so-called rigid body treatment [43]. In the rigid body treatment, conformations are considered as instantaneously rigid and the hydrodynamic theory for rigid bead models is applied [23]. According to that procedure, conformational properties (like the radius of gyration) can be rigorously calculated whereas some overall hydrodynamic entities (like the translational diffusion and the intrinsic viscosity) can be acceptably estimated [44].

The initial conformations of the RNA models (including helices and loops) were built using an in-house program. After an equilibration period, long BD simulations with a duration about 100 times the relaxation time of the molecules were carried out. The resulting trajectories were divided into five pieces in order to calculate averages and standard deviations. We computed three properties: the radius of gyration *R*
_{
g
}, the translational diffusion coefficient *D*
_{
t
}, and the harmonic relaxation time *τ*
_{
h
}.

The macroscopic translational diffusion of a Brownian particle (consequence of its Brownian motion) is characterized by the time correlation function of the center of mass or Einstein equation 〈*d*
^{2}〉=6*D*
_{
t
}
*Δ*
*t*, where 〈*d*
^{2}〉 is the mean squared displacement of the center of mass, *Δ*
*t* is the displacement duration, and *D*
_{
t
} is the translational diffusion coefficient of the particle. Long BD-HI simulations were carried out to obtain *D*
_{
t
} using that equation. On the other hand, BD-noHI simulations were used as a “smart Monte Carlo” procedure to generate equilibrium conformations and apply the rigid body treatment to obtain *τ*
_{
h
} and *R*
_{
g
}. The harmonic relaxation time *τ*
_{
h
} characterizes the overall tumbling of the molecule and can be compared to experimental rotational relaxation times. In case of a rigid body with arbitrary shape, the harmonic mean of its five characteristic rotational relaxation times is calculated as *τ*
_{
h
}=1/(6*D*
_{
r
}) [23], where *D*
_{
r
} is a rotational diffusion coefficient obtained from the eigenvalues of the rotational diffusion tensor. Although the rigid body hydrodynamics is not strictly applicable to the reorientational dynamics of flexible structures (like tRNA and rRNA), it allows for the calculation of an ensemble average 〈*D*
_{
r
}〉 from which to calculate a value *τ*
_{
h
}=1/(6〈*D*
_{
r
}〉). Finally, since *R*
_{
g
} is a conformational property, it is correctly obtained from a conformational average over a BD-noHI trajectory.

## Results and discussion

^{phe}(PDB code 1EHZ), was built according to the structure sketched in Fig. 1 a [45]. After a BD-noHI simulation for equilibrating the initial conformation, BD simulations with and without HI were performed. Figure 3 a shows two of the many conformations adopted by that molecule along the Brownian trajectory once the equilibrium has been reached (the “Additional file 1” in the supplementary material displays the movie of the trajectory). Flexibility is noticed in the variation of the angle between arms (inter-arm angle

*θ*) as better appreciated in Fig. 4 a. There, it is depicted the fluctuation along the Brownian trajectory of the angle

*θ*subtended by the two arm vectors, one defined along the “acceptor stem” and the other one defined along the “anticodon stem”. As observed, fluctuations are quick and the value of the angle varies in a relatively wide range, in agreement with some molecular dynamics simulations [11]. Figure 4 b is an histogram that shows the frequency of occurrence of the different

*θ*values as taken at intervals of 10°. Its broadness is indicative of the flexibility of the structure. Clearly, there is a most probable equilibrium value around

*θ*∼ 90° (the precise average value being

*θ*= 88°), which is coincident with that found in the crystal structure. Therefore, overall (time-averaged) hydrodynamic properties can be reasonably estimated using rigid bead models built upon the atomic coordinates coming from PDB files [46, 47].

*R*

_{ g },

*D*

_{ t }and

*τ*

_{ h }, along with values calculated previously using rigid body hydrodynamics [47] and experimental results. The rigid body hydrodynamic calculations were carried out with the program HYDROPRO [24] which is able to build a bead model from the atomic (PDB) coordinates. There is a general good agreement among HYDROPRO calculations, BD simulations and experimental results. The advantage of our BD scheme is the possibility of evaluating, in addition to global solution properties, the internal dynamics of the molecule as illustrated by Fig. 4.

Experimental and calculated values (via rigid body hydrodynamics and BD simulations) of solution properties of yeast tRNA ^{phe}

On the other hand, a bead-and-spring model of 5S rRNA of *E. Coli* (PDB code IC2X:C) was built according to the structure sketched in Fig. 1
b [20]. Again, after equilibrating the initial conformation, BD simulations with and without HI were performed. Figure 3
b illustrates two of the many conformations adopted by our 5S rRNA model along the equilibrium Brownian trajectory (the “Additional file 2” in the supplementary material displays the movie of the trajectory). The flexibility is revealed by the variation of the angle between the domains (*α*, *β*, and *γ*) as well as the expansion of the loops whereas the helical regions keep a rigid structure.

*R*

_{ g },

*D*

_{ t }and

*τ*

_{ h }coming from our BD simulations along with some values selected from Skibinska et al. [20]. In that work, authors collect experimental data from several sources and present their own experimental data. Also, they make some rigid body hydrodynamic calculations (partly carried out with the program HYDROPRO [24]) and propose a likely tertiary structure for the 5S rRNA.

Experimental and calculated values (via rigid body hydrodynamics and BD simulations) of solution properties of 5S rRNA of *E. Coli*

As observed, the results coming from our BD simulations agree relatively well with both experimental data and rigid body hydrodynamic predictions. In particular, our predictions for *D*
_{
t
} and *R*
_{
g
} are coincident with some of the experimental measurements collected in Table 2 [48, 49]. If we use our predicted *D*
_{
t
} value (6.2×10^{−7} cm ^{2}/s) and the Svedberg equation (\(s/D_{t}=M(1-\bar {v}\rho)/RT\)), we obtain a value for the sedimentation coefficient *s*≃5 S what further supports the validity of the scheme employed in this work. In order to obtain that value of *s*, we assume for the 5S rRNA a molecular weight *M*=40000 Da [50] and a partial specific volume \(\bar {v}=0.53\) cm ^{3}/g [51] as well as a solvent (water) density *ρ*≃1 g/cm ^{3}.

On the other hand our calculated *τ*
_{
h
} is smaller than the experimental value obtained by Skibinska et al. [20], although close to the value predicted by their rigid body hydrodynamic calculations. As already mentioned, the rigid body hydrodynamic treatment does not capture rigorously the reorientational dynamics of flexible molecules and it is not the most adequate procedure to calculate their rotational relaxation times. For the case of 5S rRNA, the rigid body treatment seems to underestimate the rotational relaxation time. In any case, theoretical predictions are reasonably good if one considers the existing discrepancy among the experimental data themselves. In the “Additional file 1” (yeast tRNA ^{phe}) and the “Additional file 2” (5S rRNA) of the supplementary material, one can appreciate the internal and rotational dynamics of the RNA molecules modeled in this work.

## Conclusions

We have devised a computational scheme to model RNA molecules, which present a partially disordered structure, and calculate their solution properties. According to the model, each nucleotide is replaced by a bead so that the secondary structure is represented in detail. The beads are connected by Hooke springs that form a type of elastic network. The model dynamics can be simulated by any Brownian dynamics algorithm. The algorithm and the analysis procedures employed in this work are implemented in our public-domain software SIMUFLEX available at http://leonardo.inf.um.es/macromol. The previously mentioned program HYDROPRO is also available at that web site.

The good agreement of our simulation results with some experimental data as well as with other theoretical results obtained by alternative numerical approaches ensures the validity of our scheme for the prediction of solution properties of RNAs. It is also indicative that the model is able to capture the tertiary folds of RNA molecules in solution. The advantage of the Brownian dynamics simulation technique (including hydrodynamic interaction) respect to other procedures is the possibility of simulating the long time internal dynamics of big macromolecules in real time.

In a previous work [52], we tested the adequacy of a computational scheme based on a residue level model and the Brownian dynamics simulation to predict the internal dynamics of partially disordered proteins. The present work on the more complex RNA molecule further confirms the success of that type of approach which is expected to be appropriate for the study of other disordered biological macromolecules.

## Availability of supporting data

All the supporting data are included as additional files.

## Declarations

### Acknowledgments

Financial support was provided by Ministerio de Economía y Competitividad within project CTQ2012-33717 including FEDER funds and Fundación Séneca de la Región de Murcia (grant 19353/PI/14).

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Bloomfield VA, Crothers DM, Tinoco I. Nucleic Acids. Structures, Properties and Functions. Sausalito California: University Science Books; 2000.Google Scholar
- Noller HF. Structure of ribosomal RNA. Ann Rev Biochem. 1984; 53:119–62.View ArticleGoogle Scholar
- Hyeon C, Dima RI, Thirumalai D. Size, shape and flexibility of RNA structures. J Chem Phys. 2006; 125:194905.View ArticleADSGoogle Scholar
- Kim SH, Suddath FL, Quigley GJ, McPherson A, Sussman JL, Wang AHJ, et al.Three-dimensional tertiary structure of yeast phenylalanine transfer RNA. Science. 1974; 185:435–40.View ArticleADSGoogle Scholar
- Hagerman P. Flexibility of RNA. Ann Rev Biophys Biomol Str. 1997; 26:139–56.View ArticleGoogle Scholar
- Wang CC, Ford Jr. N C, Fournier MJ. Laser light-scattering analysis of the dimerization of transfer ribonucleic acids with complementary anticodons. Biopolymers. 1981; 20:155–68.View ArticleMATHGoogle Scholar
- Patkowski A, Jen S, Chu B. Intensity-fluctuation spectroscopy and tRNA conformation. II. Changes of size and shape of tRNA in the melting process. Biopolymers. 1978; 17:2643–62.View ArticleGoogle Scholar
- Porschke D, Antosiewicz J. Permanent dipole moment of tRNA’s and variation of their structure in solution. Biophys J. 1990; 58:403–11.View ArticleGoogle Scholar
- Müller JJ, Damaschun G, Wilhelm P, Welfle H, Pilz I. Comparison of the structures of the native form of rat liver 5S rRNA and yeast tRNA: small-angle and wide-angle scattering. Int J Biol Macromol. 1982; 4:289–96.View ArticleGoogle Scholar
- Friederich MW, Vacano E, Hagerman PJ. Global flexibility of tertiary structure in RNA: yeast tRNA as a model system. Proc Natl Acad Sci U S A. 1998; 95:3572–7.View ArticleADSGoogle Scholar
- (McCammon JA, Harvey SC, editors.)1987. Dynamics of Proteins and Nucleic Acids. Cambridge: Cambridge University Press.Google Scholar
- Harvey SC, McCammon JA. Intramolecular flexibility in phenylalanine transfer RNA. Nature. 1981; 294:286–7.View ArticleADSGoogle Scholar
- Harvey SC, Prabhakaran M, Mao B, McCammon JA. Phenylalanine transfer RNA: molecular dynamics simulation. Science. 1984; 223:1189–91.View ArticleADSGoogle Scholar
- Tung C-S, Harvey SC, McCammon JA. Large-amplitude bending motions in phenylalanine transfer RNA. Biopolymers. 1984; 23:2173–93.View ArticleGoogle Scholar
- Harvey SC, Prabhakaran M, McCammon JA. Molecular-dynamics simulation of phenylalanine transfer RNA. I. Methods and general results. Biopolymers. 1985; 24:1169–88.View ArticleGoogle Scholar
- Alexander RW, Eargle J, Luthey-Schulten Z. Experimental and computational determination of tRNA dynamics. FEBS Lett. 2010; 584:376–86.View ArticleGoogle Scholar
- Agirrezabala X, Valle M. Structural insights into tRNA dynamics on the ribosome. Int J Mol Sci. 2015; 16:9866–95.View ArticleGoogle Scholar
- Fox JW, Wong K-P. The hydrodynamic shape, conformation, and molecular model of Escherichia coli ribosomal 5 S RNA. J Biol Chem. 1979; 254:10139–44.Google Scholar
- Barciszewska MZ, Erdmann VA, Barciszewski J. Ribosomal 5S RNA: tertiary structure and interactions with proteins. Biol Rev. 1996; 71:1–25.View ArticleGoogle Scholar
- Skibinska L, Banachowicz E, Gapiński J, Patkowski A, Barciszewski J. Structural similarity of E. coli 5S rRNA in solution and within the ribosome. Biopolymers. 2004; 73:316–25.View ArticleGoogle Scholar
- Takada S. Coarse-grained molecular simulations of large biomolecules. Current Opinion in Structural Biology. 2012; 22:130–7.View ArticleGoogle Scholar
- García de la Torre J, Navarro S, López Martínez MC, Díaz FG, López Cascales JJ. HYDRO: a computer software for the prediction of hydrodynamic properties of macromolecules. Biophys J. 1994; 67:530–31.View ArticleGoogle Scholar
- Carrasco B, García de la Torre J. Hydrodynamic properties of rigid particles. Comparison of different modelling and computational procedures. Biophys J. 1999; 76:3044–57.View ArticleGoogle Scholar
- García de la Torre J, Huertas ML, Carrasco B. Calculation of hydrodynamic properties of globular proteins from their atomic-level structures. Biophys J. 2000; 78:719–30.View ArticleGoogle Scholar
- García de la Torre J, Pérez Sánchez HE, Ortega A, Hernández Cifre JG, Fernandes MX, Díaz Baños FG, et al.Calculation of the solution properties of flexible macromolecules: methods and applications. Eur Biophys J. 2003; 32:477–86.View ArticleGoogle Scholar
- García de la Torre J, Ortega A, Pérez Sánchez HE, Hernández Cifre JG. MULTIHYDRO and MONTEHYDRO: Conformational search and Monte Carlo calculation of solution properties of rigid and flexible macromolecular models. Biophys Chem. 2005; 116:121–8.View ArticleGoogle Scholar
- García de la Torre J, Hernández Cifre JG, Ortega A, Rodrígez Schmidt R, Fernandes MX, Pérez Sánchez HE, et al.SIMUFLEX : Algorithms and tools for simulation of the conformation and dynamics of flexible molecules and nanoparticles in solution. J Chem Theor Comput. 2009; 5:2606–18.View ArticleGoogle Scholar
- García de la Torre J, Navarro S, López Martínez MC. Hydrodynamic properties of a double-helical model for DNA. Biophys J. 1994; 66:1573–9.View ArticleGoogle Scholar
- Tepper H, Voth GA. A coarse-grained model for double-helix molecules in solution: Spontaneous helix formation and equilibrium properties. J Chem Phys. 2005; 122:124906.View ArticleADSGoogle Scholar
- García de la Torre J, Horta A. Sedimentation coefficient and X-ray scattering of double helical model for DNA. J Phys Chem. 1976; 80:2028–35.View ArticleGoogle Scholar
- Huertas ML, Navarro S, López Martínez MC, García de la Torre J. Simulation of the conformation and dynamics of a double-helical model for DNA. Biophys J. 1997; 73:3142–53.View ArticleGoogle Scholar
- Tozzini V. Coarse-grained models for proteins. Curr Opin Struct Biol. 2005; 15:144–50.View ArticleGoogle Scholar
- Xia Z, Gradner DP, Gutell RR, Ren P. Coarse-grained model for simulation of RNA three-dimensional structures. J Phys Chem B. 2010; 114:13497–506.View ArticleGoogle Scholar
- Hyeon C, Denesyuk NA, Thirumali D. Development and applications of coarse grained models for RNA. Isr J Chem. 2014; 54:1358–73.View ArticleGoogle Scholar
- Réblová K, Lankas F, Rázga F, Krasovska MV, Koca J, Sponer J. Structure, dynamics, and elasticity of free 16S rRNA helix studied by molecular dynamics simulations. Biopolymers. 2006; 82:504–20.View ArticleGoogle Scholar
- McDowell SE, Spacková N, Sponer J, Sponer NG. Molecular dynamics simulations of RNA: An in silico single molecule approach. Biopolymers. 2006; 85:169–84.View ArticleGoogle Scholar
- Badu SR, Melnik R, Paliy M, Prabhakar S, Sebetci A, Shapiro BA. Modeling of RNA nanotubes using molecular dynamics simulation. Eur Biophys J. 2014; 43:555–64.View ArticleGoogle Scholar
- Tanaka Y, Fujii S, Hiroaki H, Sakata T, Tanaka T, Uesugi S, et al.A’-form RNA double helix in the single crystal structure of r(UGAGCUUCGGCUC). Nucleic Acids Res. 1999; 27:949–55.View ArticleGoogle Scholar
- Rey A, Freire JJ, García de la Torre J. Monte Carlo calculations for linear and star polymers. 3. Dimensions and hydrodynamic properties in good solvents. Macromolecules. 1987; 20:342–6.View ArticleADSGoogle Scholar
- Ermak DL, McCammon JA. Brownian dynamics with hydrodynamic interactions. J Chem Phys. 1978; 69:1352–60.View ArticleADSGoogle Scholar
- Rotne J, Prager S. Variational treatment of hydrodynamic interaction on polymers. J Chem Phys. 1969; 50:4831–7.View ArticleADSGoogle Scholar
- Rossky PJ, Doll JD, Friedman HL. Brownian dynamics as a smart Monte Carlo method. J Chem Phys. 1978; 69:4628–33.View ArticleADSGoogle Scholar
- Zimm BH. Chain molecule hydrodynamics by the Monte-Carlo method and the validity of the Kirkwood-Riseman approximation. Macromolecules. 1980; 13:592–602.View ArticleADSGoogle Scholar
- Fixman M. Variational bounds for polymer transport coefficients. J Chem Phys. 1983; 78:1588–93.View ArticleADSGoogle Scholar
- Shi H, Moore PB. The crystal structure of yeast phenylalanine tRNA at 1.93 Ȧ resolution: A classic structure revisited. RNA. 2000; 6:1091–105.View ArticleGoogle Scholar
- Antosiewicz J, Porschke D. An unusual electrooptical effect observed for DNA fragments and its apparent relation to a permanent electric moment associated with bent DNA. Biophys Chem. 1989; 33:19–30.View ArticleGoogle Scholar
- Fernandes MX, Ortega A, López Martínez MC, García de la Torre J. Calculation of hydrodynamic properties of small nucleic acids from their atomic structures. Nucleic Acids Res. 2002; 30:1782–8.View ArticleGoogle Scholar
- Müller JJ, Zalkova TN, Ziwer D, Misselwitz R, Gast K, Serdyuk IN, et al.Comparison of the structure of ribosomal 5S RNA from E. coli and from rat liver using x-ray scattering and dynamic light scattering. Eur Biophys J. 1986; 13:301–7.Google Scholar
- Brunel C, Romby P, Westhof E, Ehresmann C, Ehresmann J. Three-dimensional model of escherichia coli ribosomal 5 S RNA as deduced from structure probing in solution and computer modeling. J Mol Biol B. 1991; 221:293–308.View ArticleGoogle Scholar
- Erdmann VA, Barciszewski J, Szymansky M. Barciszewska MZ. 5S ribosomal RNA database. Nucleic Acids Res. 2002; 30:176–8.View ArticleGoogle Scholar
- Funari SS, Rapp G, Perbandt M, Dierks K, Vallazza M, Betzel C, et al.Structure of free thermus flavus 5S rRNA at 1.3 nm resolution from synchrotron solution scattering. J Biol Chem. 2000; 275:31283–8.View ArticleGoogle Scholar
- Amorós D, Ortega A, García de la Torre J. Prediction of hydrodynamic and other solution properties of partially disordered proteins with a simple, coarse-grained model. J Chem Theory Comput. 2013; 9:1678–85.View ArticleGoogle Scholar