ASYMMETRIC PERIODIC BOUNDARY CONDITIONS FOR MOLECULAR DYNAMICS AND COARSE-GRAINED SIMULATIONS OF NUCLEIC ACIDS
Main Article Content
Abstract
Periodic boundary conditions are commonly applied in molecular dynamics simulations in the microcanonical (NVE), canonical (NVT) and isothermal-isobaric (NpT) ensembles. In their simplest application, a biological system of interest is placed in the middle of a solvation box, which is chosen ‘sufficiently large’ to minimize any numerical artefacts associated with the periodic boundary conditions. This practical approach brings limitations to the size of biological systems that can be simulated. Here, we study simulations of effectively infinitelylong nucleic acids, which are solvated in the directions perpendicular to the polymer chain, while periodic boundary conditions are also applied along the polymer chain. We study the effects of these asymmetric periodic boundary conditions (APBC) on the simulated results, including the mechanical properties of biopolymers and the properties of the surrounding solvent. To get some further insights into the advantages of using the APBC, a coarse-grained worm-like chain model is first studied, illustrating how the persistence length can be extracted from local properties of the polymer chain, which are less affected by the APBC than some global averages. This is followed by all-atom molecular dynamics simulations of DNA in ionic solutions, where we use the APBC to investigate sequence-dependent properties of DNA molecules and properties of the surrounding solvent.
Downloads
Metrics
Article Details
This work is licensed under a Creative Commons Attribution 4.0 International License.
You are free to:
- Share — copy and redistribute the material in any medium or format for any purpose, even commercially.
- Adapt — remix, transform, and build upon the material for any purpose, even commercially.
- The licensor cannot revoke these freedoms as long as you follow the license terms.
Under the following terms:
- Attribution — You must give appropriate credit , provide a link to the license, and indicate if changes were made . You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
- No additional restrictions — You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits.
Notices:
You do not have to comply with the license for elements of the material in the public domain or where your use is permitted by an applicable exception or limitation .
No warranties are given. The license may not give you all of the permissions necessary for your intended use. For example, other rights such as publicity, privacy, or moral rights may limit how you use the material.
References
(1) Vologodskii, A. Biophysics of DNA;
(2) Jacobson, D.; Saleh, O. Counting the
ions surrounding nucleic acids. Nucleic
Acids Research 2016, 45, 1596–1605.
(3) Mocci, F.; Laaksonen, A. Insight into
nucleic acid counterion interactions from
inside molecular dynamics simulations is
“worth its salt”. Soft Matter 2012, 8,
–9284.
(4) Várnai, P.; Zakrzewska, K. DNA and
its counterions: a molecular dynamics
study. Nucleic Acids Research 2004, 32,
–4280. (5) Zavadlav, J.; Podgornik,
R.; Praprotnik, M. Adaptive resolution
simulation of a DNA molecule in salt
solution. Journal of Chemical Theory and
Computation 2015, 11, 5035–5044. (6)
Zavadlav, J.; Sablic, J.; Podgornik, R.;
Praprotnik, M. Open-Boundary Molecular
Dynamics of a DNA Molecule in a Hybrid
Explicit/Implicit Salt Solution.
Biophysical Journal 2018, 114, 2352–
(7) Dans, P.; Walther, J.; Gómez, H.;
Orozco, M. Multiscale simulation of DNA. Current Opinion in Structural Biology
, 37, 29–45.
(8) Korolev, N.; Nordenskiöld, L.;
Lyubartsev, A. Multiscale coarse-grained
modelling of chromatin components: DNA
and the nucleosome. Advances in Colloid
and Interface Science 2016, 232, 36–48.
(9) Poppleton, R.; Matthies, M.; Mandal,
D.; Romano, F.; Šulc, P.; Rovigatti, L.
oxDNA: coarse-grained simulations of
nucleic acids made simple. The Journal of
Open Source Software 2023, 8, 4693.
(10) Sengar, A.; Ouldridge, T.; Henrich,
O.; Rovigatti, L.; Šulc, P. A Primer on the
oxDNA Model of DNA: When to Use it,
How to Simulate it and How to Interpret
the Results. Frontiers in Molecular
Biosciences 2021, 8 . (11) Kovaleva, N.;
Koroleva, I.; Mazo, M.; Zubova, E. The
“sugar” coarse-grained DNA model.
Journal of Molecular Modelling 2017, 23,
19
(12) Rolls, E.; Togashi, Y.; Erban, R.
Varying the resolution of the Rouse model
on temporal and spatial scales: application
to multiscale modelling of DNA dynamics.
Multiscale Modeling and Simulation 2017,
, 1672–1693.
(13) Maffeo, C.; Aksimentiev, A. MrDNA:
a multi-resolution model for predicting the
structure and dynamics of DNA systems.
Nucleic Acids Research 2020, 48, 5135–
(14) Minhas, V.; Sun, T.; Mirzoev,
A.; Korolev, N.; Lyubartsev, A.;
Nordenskiöld, L. Modeling DNA
Flexibility: Comparison of Force Fields
from Atomistic to Multiscale Levels.
Journal of Physical Chemistry 2020, 124,
–49.
(15) Lee, S.; Rasaiah, J. Molecular
dynamics simulation of ion mobility. 2.
alkali metal and halide ions using the
SPC/E model for water at 25◦ C. Journal of
Physical Chemistry 1996, 100, 1420–1425.