Asymmetric periodic boundary conditions for molecular dynamics and coarse-grained simulations of nucleic acids

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 infinitely-long 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.

[1]  M. Matthies,et al.  oxDNA: coarse-grained simulations of nucleic acids made simple , 2023, J. Open Source Softw..

[2]  G. Papoian,et al.  On Stretching, Bending, Shearing and Twisting of Actin Filaments II: Multi-Resolution Modelling , 2022, 2203.01284.

[3]  G. Papoian,et al.  On Stretching, Bending, Shearing, and Twisting of Actin Filaments I: Variational Models. , 2021, Journal of chemical theory and computation.

[4]  T. Ouldridge,et al.  A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results , 2021, Frontiers in Molecular Biosciences.

[5]  Yi Wang,et al.  Scalable molecular dynamics on CPU and GPU architectures with NAMD. , 2020, The Journal of chemical physics.

[6]  Miho M. Suzuki,et al.  Structural Dynamics of DNA Depending on Methylation Pattern , 2020, bioRxiv.

[7]  Christopher Maffeo,et al.  MrDNA: a multi-resolution model for predicting the structure and dynamics of DNA systems , 2020, Nucleic acids research.

[8]  Lars Nordenskiöld,et al.  Modeling DNA Flexibility: Comparison of Force Fields from Atomistic to Multiscale Levels. , 2019, The journal of physical chemistry. B.

[9]  R. Erban,et al.  Stochastic Modelling of Reaction–Diffusion Processes , 2019 .

[10]  R. Erban.  Coarse-graining molecular dynamics: stochastic models with non-Gaussian force distributions , 2019, Journal of Mathematical Biology.

[11]  Xiang-Jun Lu,et al.  Web 3DNA 2.0 for the analysis, visualization, and modeling of 3D nucleic acid structures , 2019, Nucleic Acids Res..

[12]  M. Praprotnik,et al.  Open-Boundary Molecular Dynamics of a DNA Molecule in a Hybrid Explicit/Implicit Salt Solution. , 2018, Biophysical journal.

[13]  R. Pérez,et al.  Understanding the mechanical response of double-stranded DNA and RNA under constant stretching forces using all-atom molecular dynamics , 2017, Proceedings of the National Academy of Sciences.

[14]  O. Saleh,et al.  Counting the ions surrounding nucleic acids , 2016, Nucleic acids research.

[15]  Radek Erban,et al.  Varying the Resolution of the Rouse Model on Temporal and Spatial Scales: Application to Multiscale Modeling of DNA Dynamics , 2016, Multiscale Model. Simul..

[16]  A. Lyubartsev,et al.  Multiscale coarse-grained modelling of chromatin components: DNA and the nucleosome. , 2016, Advances in colloid and interface science.

[17]  Modesto Orozco,et al.  Multiscale simulation of DNA. , 2016, Current opinion in structural biology.

[18]  Matej Praprotnik,et al.  Adaptive resolution simulation of a DNA molecule in salt solution. , 2015, Journal of chemical theory and computation.

[19]  Radek Erban,et al.  Coupling all-atom molecular dynamics simulations of ions in water with Brownian dynamics , 2015, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[20]  Alexander Vologodskii,et al.  Biophysics of DNA , 2015 .

[21]  Sten Rüdiger,et al.  Particle-Based Multiscale Modeling of Calcium Puff Dynamics , 2015, Multiscale Model. Simul..

[22]  Matej Praprotnik,et al.  Adaptive Resolution Simulation of MARTINI Solvents. , 2014, Journal of chemical theory and computation.

[23]  Matej Praprotnik,et al.  Adaptive resolution simulation of an atomistic protein in MARTINI water. , 2014, The Journal of chemical physics.

[24]  E. A. Zubova,et al.  The “sugar” coarse-grained DNA model , 2014, Journal of Molecular Modeling.

[25]  Aatto Laaksonen,et al.  Insight into nucleic acid counterion interactions from inside molecular dynamics simulations is “worth its salt” , 2012 .

[26]  Alexander Vologodskii,et al.  Sequence dependence of DNA bending rigidity , 2010, Proceedings of the National Academy of Sciences.

[27]  Xiang-Jun Lu,et al.  3DNA: a versatile, integrated software system for the analysis, rebuilding and visualization of three-dimensional nucleic-acid structures , 2008, Nature Protocols.

[28]  A. Lyubartsev,et al.  A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations. , 2003, Nucleic acids research.

[29]  A. Lyubartsev,et al.  Molecular dynamics simulations of DNA in solutions with different counter-ions. , 1998, Journal of biomolecular structure & dynamics.

[30]  J. Rasaiah,et al.  Molecular Dynamics Simulation of Ion Mobility. 2. Alkali Metal and Halide Ions Using the SPC/E Model for Water at 25 °C† , 1996 .

[31]  J. W. Halley,et al.  Molecular dynamics simulation of water beween two ideal classical metal walls , 1989 .

[32]  Aidan Twomey,et al.  On the Stochastic Modelling of Reaction-Diffusion Processes , 2007 .

[33]  Krystyna Zakrzewska,et al.  DNA and its counterions: a molecular dynamics study. , 2004, Nucleic acids research.