RNA model: Difference between revisions

From OxDNA
Jump to navigation Jump to search
Line 45: Line 45:




===Running simulation with the oxRNA model===
===Running a simulation with the oxRNA model===


The oxRNA model is integrated into the oxDNA simulation code. In particular, it is possible to use Virtual Move Monte Carlo (VMMC), Monte Carlo (MC) and Molecular Dynamics (MD) simulation algorithms using the same format of input file as for the DNA model. The format of the configuration files is also the same as for the DNA model, described in [[Documentation]]. When running simulations of the oxRNA model, the following additional line must be included in the input file to specify the RNA model:
The oxRNA model is integrated into the oxDNA simulation code. In particular, it is possible to use the Virtual Move Monte Carlo (VMMC), Monte Carlo (MC) and Molecular Dynamics (MD) simulation algorithms using the same format of input file as for the DNA model. The format of the configuration files is also the same as for the DNA model, described in [[Documentation]]. When running simulations of the oxRNA model, the following additional line must be included in the input file to specify that the RNA model is to be used:
<pre>
<pre>
interaction_type = RNA
interaction_type = RNA
</pre>
</pre>


The RNA model comes with two parametrizations, the average-base and sequence-dependent one. In the average-base parametrization, the <math>V_{\rm  H.B.}</math> interaction strengths are the same for all Watson-Crick and wobble base pairs (AU, GC, GU) and 0 for all other types of base pairs and the interaction strengths have the same strength for all possible pairs of nucleotides interacting with stacking interaction <math>V_{\rm stack}</math>.  
The RNA model comes with two parametrizations, the average-base and sequence-dependent one. In the average-base parametrization, the <math>V_{\rm  H.B.}</math> interaction strengths are the same for all Watson-Crick and wobble base pairs (AU, GC, GU) and 0 for all other types of base pairs, and the interaction strengths have the same strength for all possible pairs of nucleotides interacting with the stacking interaction <math>V_{\rm stack}</math>.  
In the sequence-dependent version of the model, the interaction strengths of <math>V_{\rm stack}</math> and <math>V_{\rm H.B.}</math> depend on the type of interacting bases (interactions for <math>V_{\rm H.B.}</math> are still 0 for base pairs other than  AU, GC or GU).
In the sequence-dependent version of the model, the interaction strengths of <math>V_{\rm stack}</math> and <math>V_{\rm H.B.}</math> depend on the type of interacting bases (interactions for <math>V_{\rm H.B.}</math> are still 0 for base pairs other than  AU, GC or GU).


Line 61: Line 61:
</pre>
</pre>


Note that the file <tt>rna_sequence_dependent_parameters.txt</tt> needs to be located in the directory where you run the simulation, or full location of the file needs to be specified in <tt>seq_dep_file</tt> option.
Note that the file <tt>rna_sequence_dependent_parameters.txt</tt> needs to be located in the directory where you run the simulation, or the full location of the file needs to be specified in the <tt>seq_dep_file</tt> option.


Furthermore, the initial configuration files need to be generated so that the nucleotides are positioned in an arrangement that satisfies RNA potentials (for instance in the case of a duplex, they need to be initialized in an A-helical structure). For this purpose, a script  <tt>generate-RNA.py</tt> is provided in <tt>UTILS/</tt> subdirectory of the source code main directory.
Furthermore, the initial configuration files need to be generated so that the nucleotides are positioned in an arrangement that satisfies the RNA potentials (for instance in the case of a duplex, they need to be initialized in an A-helical structure). For this purpose, a script  <tt>generate-RNA.py</tt> is provided in the <tt>UTILS/</tt> subdirectory of the source code main directory.
For instance, if one wants to generate an initial configuration consisting of three strands, two of them complementary (with sequence 3'-GCAAGUCG-5' and its complementary) in a duplex configuration, and one single strand with sequence 3'-ACCCGU-5', one needs to create the following file text file, called for example <tt>sequences.txt</tt>:
For instance, if one wants to generate an initial configuration consisting of three strands, two of them complementary (with sequence 3'-GCAAGUCG-5' and its complementary) in a duplex configuration, and one single strand with sequence 3'-ACCCGU-5', one needs to create the following text file, called for example <tt>sequences.txt</tt>:
<pre>
<pre>
DOUBLE GCAAGUCG
DOUBLE GCAAGUCG
Line 75: Line 75:
</pre>
</pre>


which will create the configuration files. Those can then be used as initial configuration for the simulations. Other input command options that apply to oxDNA, such as use of external forces, apply with the same syntax to oxRNA as well and are described in more detail in the [[Documentation]].
which will create the configuration files. These can then be used as an initial configuration for a simulation. Other input file options that apply to oxDNA, such external_forces=1 (for the use of external forces), can be used with oxRNA with the same syntax (see [[Documentation]] for a full list and further details).





Revision as of 17:10, 11 March 2014

Description of the oxRNA model

The RNA model, oxRNA, treats each RNA nucleotide as a single rigid body with multiple interaction sites, following the coarse-graining approach adopted for the DNA model.

The nucleotides interact with the following pairwise interaction potentials:

  1. Backbone connectivity Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm backbone~}} ,
  2. Excluded volume Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm exc~}} ,
  3. Hydrogen bonding Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm H.B.~}} ,
  4. Nearest-neighbor stacking Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm stack~}} ,
  5. Cross-stacking in a duplex Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm cross~st.}} ,
  6. Coaxial stacking Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm cx.~stack}} .

which are schematically illustrated in the picture:

Image duplex combined annotated.png


Simulation units

The code uses dimensionless energy, mass, length and timescales for convenience. The relationship between simulation units (SU) and SI units is given below.

Simulation unit Physical unit
1 unit of length 8.4x10Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ^{-10}} m
1 unit of energy 4.142x10Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ^{-20}} J
1 unit of temperature 3000 K
1 unit of force 49.3x10Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ^{-11}} N
1 unit of mass 5.34x10Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ^{-25}} kg
1 unit of time 3.06x10Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle ^{-12}} s


Running a simulation with the oxRNA model

The oxRNA model is integrated into the oxDNA simulation code. In particular, it is possible to use the Virtual Move Monte Carlo (VMMC), Monte Carlo (MC) and Molecular Dynamics (MD) simulation algorithms using the same format of input file as for the DNA model. The format of the configuration files is also the same as for the DNA model, described in Documentation. When running simulations of the oxRNA model, the following additional line must be included in the input file to specify that the RNA model is to be used:

interaction_type = RNA

The RNA model comes with two parametrizations, the average-base and sequence-dependent one. In the average-base parametrization, the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm H.B.}} interaction strengths are the same for all Watson-Crick and wobble base pairs (AU, GC, GU) and 0 for all other types of base pairs, and the interaction strengths have the same strength for all possible pairs of nucleotides interacting with the stacking interaction Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm stack}} . In the sequence-dependent version of the model, the interaction strengths of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm stack}} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm H.B.}} depend on the type of interacting bases (interactions for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{\rm H.B.}} are still 0 for base pairs other than AU, GC or GU).

The average-base parametrization is used by default. In order to use the sequence-dependent version of the model, the following options need to be added into the input file:

use_average_seq = 0
seq_dep_file = rna_sequence_dependent_parameters.txt

Note that the file rna_sequence_dependent_parameters.txt needs to be located in the directory where you run the simulation, or the full location of the file needs to be specified in the seq_dep_file option.

Furthermore, the initial configuration files need to be generated so that the nucleotides are positioned in an arrangement that satisfies the RNA potentials (for instance in the case of a duplex, they need to be initialized in an A-helical structure). For this purpose, a script generate-RNA.py is provided in the UTILS/ subdirectory of the source code main directory. For instance, if one wants to generate an initial configuration consisting of three strands, two of them complementary (with sequence 3'-GCAAGUCG-5' and its complementary) in a duplex configuration, and one single strand with sequence 3'-ACCCGU-5', one needs to create the following text file, called for example sequences.txt:

DOUBLE GCAAGUCG
ACCCGU

Note that the sequences are always specified in 3'-5' order. In order to create the initial configuration files generated.top and generated.conf with the duplex and single strand randomly placed in a simulation cube with side of length 20 in simulation units, run the script

generate-RNA.py sequences.txt generated 20.0

which will create the configuration files. These can then be used as an initial configuration for a simulation. Other input file options that apply to oxDNA, such external_forces=1 (for the use of external forces), can be used with oxRNA with the same syntax (see Documentation for a full list and further details).


For an example on how to use VMMC simulations to determine the melting temperature of an RNA duplex, please see the RNA duplex melting tutorial.

Visualization of RNA configurations

In order to visualize the configurations of the oxRNA model, one can use the traj2chimera.py script, as described for the oxDNA model. It is however necessary to first set environment variable OXRNA to 1 in order for the script to properly generate visual representation of oxRNA:

export OXRNA=1

The visualization of configuration specified in, for example, generated.top and generated.conf can be then obtained by running

traj2chimera.py generated.conf generated.top 

in the UTILS/ directory which creates files generated.conf.pdb and chimera.com which can then be visualized with Chimera software by running the following command:

chimera generated.conf.pdb chimera.com

or alternatively, you can load generated.conf.pdb in the Chimera software and then click on Tools->General Controls->Command line and specify

read chimera.com

in the command line, where chimera.com needs to be present in the directory where you started Chimera.


References

The model and its performance is discussed in detail in the following reference:

P. Šulc, F. Romano, T. E. Ouldridge, J. P. K. Doye, A. A. Louis: A nucleotide-level coarse-grained model of RNA, submitted