An elastic rod model for semi-flexible polymers is presented. Theory for a continuum rod is reviewed, and it is shown that a popular discretised model used in numerical simulations gives the correct continuum limit. Correlation functions relating to both bending and twisting of the rod are derived for both continuous and discrete cases, and results are compared with numerical simulations. Finally, two possible implementations of the discretised model in the multi-purpose molecular dynamics software package LAMMPS are described.
REFERENCES
1.
J. F.
Marko
and E. D.
Siggia
, Macromolecules
28
, 8759
(1995
).2.
C.
Bustamante
, J. F.
Marko
, E. D.
Siggia
, and S.
Smith
, Science
265
, 1599
(1994
).3.
P. L.
Hansen
and R.
Podgornik
, J. Chem. Phys.
114
, 8637
(2001
).4.
A.
Rosa
, T. X.
Hoang
, D.
Marenduzzo
, and A.
Maritan
, Macromolecules
36
, 10095
(2003
).5.
N. M.
Toan
, D.
Marenduzzo
, and C.
Micheletti
, Biophys. J.
89
, 80
(2005
).6.
C. R.
Calladine
and H.
Drew
, Understanding DNA
(Academic Press
, San Diego
, 1992
).7.
J. F.
Marko
and E. D.
Siggia
, Science
265
, 506
(1994
).8.
J. F.
Marko
and E. D.
Siggia
, Macromolecules
27
, 981
(1994
).9.
N.
Gilbert
and J.
Allan
, Curr. Opin. Genet. Dev.
25
, 15
(2014
).10.
J. F.
Marko
and E. D.
Siggia
, Phys. Rev. E
52
, 2912
(1995
).11.
S.
Forth
, C.
Deufel
, M. Y.
Sheinin
, B.
Daniels
, J. P.
Sethna
, and M. D.
Wang
, Phys. Rev. Lett.
100
, 148301
(2008
).12.
L.
Postow
, C. D.
Hardy
, J.
Arsuaga
, and N. R.
Cozzarelli
, Genes Dev.
18
, 1766
(2004
).13.
T. B. K.
Le
, M. V.
Imakaev
, L. A.
Mirny
, and M. T.
Laub
, Science
342
, 731
(2013
).14.
L. F.
Liu
and J. C.
Wang
, Proc. Natl. Acad. Sci. U.S.A.
84
, 7024
(1987
).15.
Y. P.
Tsao
, H. Y.
Wu
, and L. F.
Liu
, Cell.
56
, 111
(1989
).16.
S.
Deng
, R. A.
Stein
, and N. P.
Higgins
, Mol. Microbiol.
57
, 1511
(2005
).17.
S. P.
Mielke
, W. H.
Fink
, V. V.
Krishnan
, N.
Groenbech-Jensen
, and C. J.
Benham
, J. Chem. Phys.
121
, 8104
(2004
).18.
C.
Naughton
et al., Nat. Struct. Mol. Biol.
20
, 387
(2013
).19.
20.
J.
Roca
, Chromosoma
120
, 323
(2011
).21.
N.
Mondal
and J. D.
Parvin
, Nature (London)
413
, 435
(2001
).22.
C.
Bouchiat
and M.
Mézard
, Eur. Phys. J. E
2
, 377
(2000
).23.
B. S.
Fujimoto
and J. M.
Schurr
, Nature (London)
344
, 175
(1990
).24.
J.
Lipfert
, J. W. J.
Kerssemakers
, T.
Jager
, and N. H.
Dekker
, Nat. Methods
7
, 977
(2010
).25.
T. R.
Powers
, Rev. Mod. Phys.
82
, 1607
(2010
).26.
W. K.
Olson
, Curr. Opin. Struct. Biol.
6
, 242
(1996
).27.
R. K. Z.
Tan
and S. C.
Harvey
, J. Mol. Biol.
205
, 573
(1989
).28.
G. C.
Rollins
, A. S.
Petrov
, and S. C.
Harvey
, Biophys. J.
94
, L38
(2008
).29.
T. A.
Knotts
, N.
Rathore
, D. C.
Schwartz
, and J. J.
de Pablo
, J. Chem. Phys.
126
, 084901
(2007
).30.
J.
Huang
, T.
Schlick
, and A.
Vologodskii
, Proc. Natl. Acad. Sci. U.S.A.
98
, 968
(2001
).31.
T.
Schlick
and W. K.
Olson
, Science
257
, 1110
(1992
).32.
T.
Schlick
, Curr. Opin. Struct. Biol.
5
, 245
(1995
).33.
S.
Rappaport
and Y.
Rabin
, Macromolecules
37
, 7847
(2004
).34.
G.
Chirico
and J.
Langowski
, Biopolymers
34
, 415
(1994
).35.
G.
Chirico
and J.
Langowski
, Biophys. J.
71
, 955
(1996
).36.
S.
Allison
, R.
Austin
, and M.
Hogan
, J. Chem. Phys.
90
, 3843
(1989
).37.
38.
S. M.
Rappaport
and Y.
Rabin
, J. Phys. A
40
, 4455
(2007
).39.
In the case of a linear rod there are N position vectors describing the vertices of N − 1 segments, with tangents
$\mathbf {\hat{u}}_i$
with i = 1, …, N − 1. For a circular rod there are N segments and $\mathbf {\hat{u}}_N=(\mathbf {r}_1-\mathbf {r}_N)/a$
.40.
D.
Marenduzzo
, C.
Micheletti
, H.
Seyed-allaei
, A.
Trovato
, and A.
Maritan
, J. Phys. A
38
, L277
(2005
).41.
F.
Benedetti
, J.
Dorier
, Y.
Burnier
, and A.
Stasiak
, Nucl. Aci. Res.
42
, 2848
(2013
).42.
C. A.
Brackley
, J.
Allan
, D.
Keszenman-Pereyra
, and D.
Marenduzzo
, “How to make chromatin: Topological enzymes aid chromatin self-assembly in silico
” (unpublished).43.
The two potentials are however equivalent in the limit where all of αi, βi, and γi are small, a regime which is forced by the Hamiltonian considering either of the dihedrals.
44.
The angles ζi are minimised by choosing a large alignment energy κa; in practice, this can be set as large as possible subject to numerical stability.
45.
Source code for the LAMMPS orientation angle style discussed in Sec. V B is available on request from the authors.
46.
J. D.
Moroz
and P.
Nelson
, Proc. Natl. Acad. Sci. U.S.A.
94
, 14418
(1997
).47.
S. A.
Koehler
and T. R.
Powers
, Phys. Rev. Lett.
85
, 4827
(2000
).48.
C. W.
Wolgemuth
, T. R.
Powers
, and R. E.
Goldstein
, Phys. Rev. Lett.
84
, 1623
(2000
).49.
R. D.
Kamien
, Rev. Mod. Phys.
74
, 953
(2002
).50.
G. H.
Fredrickson
, The Equilibrium Theory of Inhomogeneous Polymers
(Oxford University Press
, New York
, 2006
).51.
If R is a rotation matrix for rotating a vector
$\mathbf {\hat{u}}$
by an angle δϕ about some axis $\mathbf {\hat{p}}$
, then $\delta \mathbf {\hat{u}}=(R-I)\mathbf {\hat{u}}$
. In the limit δϕ → 0, $\mathbf {\hat{p}}\times \mathbf {\hat{u}}\delta \phi =(R-I)\mathbf {\hat{u}}$
.© 2014 AIP Publishing LLC.
2014
AIP Publishing LLC
You do not currently have access to this content.