← Back to branch summary

~ubuntu-branches/debian/sid/lammps/sid

« back to all changes in this revision

Viewing changes to tools/moltemplate/examples/coarse_grained_examples/chromosome_metaphase_Naumova2013/README_FIRST.txt

  • Committer: Package Import Robot
  • Author(s): Anton Gladky
  • Date: 2015-04-29 23:44:49 UTC
  • mfrom: (5.1.3 experimental)
  • Revision ID: package-import@ubuntu.com-20150429234449-mbhy9utku6hp6oq8
Tags: 0~20150313.gitfa668e1-1
Upload into unstable.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tools/moltemplate/examples/coarse_grained_examples/chromosome_metaphase_Naumova2013/README_FIRST.txt
 
1
This is an implementation of the "two-stage" model used by Maxim Imakaev
 
2
in the Naumova et Al 2013 Science paper on metaphase chromatin.
 
3
(Download the supplemental materials section and scroll down to the section:
 
4
 "Two-stage process: linear compaction - axial compression")
 
5
 
 
6
---- SMALL MODIFICATION ----
 
7
 
 
8
Unlike that study, I did not use "softened" Lennard-Jones potentials
 
9
(which allow the chains to pass through each other).
 
10
 
 
11
--- Why use moltemplate? ---
 
12
 
 
13
Honestly, you don't need to use moltemplate to build this polymer.
 
14
It is almost counter-productive to use moltemplate to build this kind of
 
15
polymer because it is so simple.  (The polymer has only 1 bead per atom.
 
16
It just makes it more complicated to introduce all these extra 
 
17
files including monomer.lt, condensin.lt and system.lt, especially considering
 
18
that system.lt is a complex file which is generated by a separate script.)
 
19
 
 
20
However building the sytem using moltemplate may pay off if you
 
21
replace each point-like monomer with a multi-atom molecule later on.
 
22
(Right now, using moltemplate to build this system is sort of overkill.
 
23
 I'll post an example of building more complex models of chromatin eventually.)
 
24
 
 
25
Anyway, the two-stage model at the end of Naumova et al Science 2013 uses the "30nm-fiber" model, whose details are (somewhat vaguely) described in the supplemental materials section.
 
26
 
 
27
---- 10-nm fiber model: ----
 
28
 
 
29
For the 10nm model, 
 
30
    n=128000, 
 
31
    L=200, 
 
32
    U(alpha)=5*(1 - cos(alpha))
 
33
    bond_length=1.0     (=10nm)
 
34
    sigma=1.0 (particle radius = 10nm)
 
35
 
 
36
---- 30-nm fiber model: ----
 
37
 
 
38
"The 30nm-like fiber was modeled by increasing the volume of each monomer and the amount of DNA represented by each monomer by a factor of 4.25, while keeping other parameters the same at the monomer level."
 
39
 
 
40
I interpret this to mean that, for the 30nm model, 
 
41
    n=128000/4.25~=30117   (however I rounded up to 32768=2^15)
 
42
    L=200/4.25~=47         (however I rounded up to 51)
 
43
    U(alpha)=1.17647*(1 - cos(alpha))    (5/4.25=1.17647)
 
44
 
 
45
To increase the volume by a factor o 4.25, I increase both the diameter of each
 
46
bead (the "sigma" parameter), and the bond-lengths connecting them from 
 
47
1.0 (corresponding to 10nm) to 4.25^(1/3)~=1.6198 (corresponding to 16.198nm).
 
48