Difference between revisions of "Model introduction"

From OxDNA
Line 1: Line 1:
The model is discussed in detail in:
+
The model treats DNA as a string of rigid nucleotides, which interact through potentials which depend on the position and orientation of the nucleotides. The interactions are:
 +
#Sugar-phosphate backbone connectivity,
 +
#Excluded volume,
 +
#Hydrogen bonding,
 +
#Nearest-neighbour stacking,
 +
#Cross-stacking between base-pair steps in a duplex,
 +
#Coaxial stacking.
 +
 
 +
This interactions are illustrated in PICTURE. Orientational modulations of the stacking potential encourage the bases to form coplanar stacks, the twist arising from the different length scales of the backbone separation and the optimum stacking separation. The possibility of unstacking allows single strands to be very flexible. Hydrogen bonding  can occur between complementary bases when they are anti-aligned, leading to the formation of double helical structures.
 +
 
 +
In the original code, all complementary base pairs and stacking partners interact with the same strength (there is no attractive interaction between non-complementary bases). A preliminary sequence-dependent parameterisation of the hydrogen-bonding and stacking interactions is included as an option in the code release: a paper discussing this parameterisation and its effects is currently in preparation.
 +
 
 +
The model does not incorporate the differentiation between the major and minor grooves of DNA double helices, and incorporates electrostatics only through the short-ranged excluded volume. For this reason, it is only appropriate for the study of systems at high salt concentration, when electrostatic interactions are strongly screened.
 +
 
 +
+LENGTHSCALES IN THE MODEL
 +
 
 +
The model and its performance is discussed in detail in the following references (the thesis provides the most complete analysis):
  
 
T. E. Ouldridge, D.Phil. Thesis, University of Oxford, 2011.
 
T. E. Ouldridge, D.Phil. Thesis, University of Oxford, 2011.
 
[http://ora.ox.ac.uk/objects/uuid:b2415bb2-7975-4f59-b5e2-8c022b4a3719 Coarse-grained modelling of DNA and DNA self-assembly]
 
[http://ora.ox.ac.uk/objects/uuid:b2415bb2-7975-4f59-b5e2-8c022b4a3719 Coarse-grained modelling of DNA and DNA self-assembly]
  
T.E. Ouldridge, A.A. Louis and J.P.K. Doye, J. Chem. Phys, 134, 085101 (2011)
+
T. E. Ouldridge, A. A. Louis and J. P. K. Doye, J. Chem. Phys, 134, 085101 (2011)
 
[http://link.aip.org/link/?JCP/134/085101 Structural, mechanical and thermodynamic properties of a coarse-grained DNA model] ([http://arxiv.org/abs/arXiv:1009.4480 arXiv])
 
[http://link.aip.org/link/?JCP/134/085101 Structural, mechanical and thermodynamic properties of a coarse-grained DNA model] ([http://arxiv.org/abs/arXiv:1009.4480 arXiv])

Revision as of 15:57, 16 April 2012

The model treats DNA as a string of rigid nucleotides, which interact through potentials which depend on the position and orientation of the nucleotides. The interactions are:

  1. Sugar-phosphate backbone connectivity,
  2. Excluded volume,
  3. Hydrogen bonding,
  4. Nearest-neighbour stacking,
  5. Cross-stacking between base-pair steps in a duplex,
  6. Coaxial stacking.

This interactions are illustrated in PICTURE. Orientational modulations of the stacking potential encourage the bases to form coplanar stacks, the twist arising from the different length scales of the backbone separation and the optimum stacking separation. The possibility of unstacking allows single strands to be very flexible. Hydrogen bonding can occur between complementary bases when they are anti-aligned, leading to the formation of double helical structures.

In the original code, all complementary base pairs and stacking partners interact with the same strength (there is no attractive interaction between non-complementary bases). A preliminary sequence-dependent parameterisation of the hydrogen-bonding and stacking interactions is included as an option in the code release: a paper discussing this parameterisation and its effects is currently in preparation.

The model does not incorporate the differentiation between the major and minor grooves of DNA double helices, and incorporates electrostatics only through the short-ranged excluded volume. For this reason, it is only appropriate for the study of systems at high salt concentration, when electrostatic interactions are strongly screened.

+LENGTHSCALES IN THE MODEL

The model and its performance is discussed in detail in the following references (the thesis provides the most complete analysis):

T. E. Ouldridge, D.Phil. Thesis, University of Oxford, 2011. Coarse-grained modelling of DNA and DNA self-assembly

T. E. Ouldridge, A. A. Louis and J. P. K. Doye, J. Chem. Phys, 134, 085101 (2011) Structural, mechanical and thermodynamic properties of a coarse-grained DNA model (arXiv)