We study DNA self-assembly and DNA computation using a coarse-grained DNA model within the directional dynamic bonding framework {[}C. Svaneborg, Comp. Phys. Comm. 183, 1793 (2012){]}. In our model, a single nucleotide or domain is represented by a single interaction site. Complementary sites can reversibly hybridize and dehybridize during a simulation. This bond dynamics induces a dynamics of the angular and dihedral bonds, that model the collective effects of chemical structure on the hybridization dynamics. We use the DNA model to perform simulations of the self-assembly kinetics of DNA tetrahedra, an icosahedron, as well as strand displacement operations used in DNA computation.
During the last decade coarse-grained nucleotide models have emerged that allow us to DNA and RNA on unprecedented time and length scales. Among them is oxDNA, a coarse-grained, sequence-specific model that captures the hybridisation transition of DNA and many structural properties of single- and double-stranded DNA. oxDNA was previously only available as standalone software, but has now been implemented into the popular LAMMPS molecular dynamics code. This article describes the new implementation and analyses its parallel performance. Practical applications are presented that focus on single-stranded DNA, an area of research which has been so far under-investigated. The LAMMPS implementation of oxDNA lowers the entry barrier for using the oxDNA model significantly, facilitates future code development and interfacing with existing LAMMPS functionality as well as other coarse-grained and atomistic DNA models.
The flexibility and stiffness of small DNA play a fundamental role ranging from several biophysical processes to nano-technological applications. Here, we estimate the mechanical properties of short double-stranded DNA (dsDNA) having length ranging from 12 base-pairs (bps) to 56 bps, paranemic crossover (PX) DNA, and hexagonal DNA nanotubes (DNTs) using two widely used coarse-grain models $-$ Martini and oxDNA. To calculate the persistence length ($L_p$) and the stretch modulus ($gamma$) of the dsDNA, we incorporate the worm-like chain and elastic rod model, while for DNT, we implement our previously developed theoretical framework. We compare and contrast all the results with previously reported all-atom molecular dynamics (MD) simulation and experimental results. The mechanical properties of dsDNA ($L_p$ $sim$ 50nm, $gamma sim$ 800-1500 pN), PX DNA ($gamma sim$ 1600-2000 pN) and DNTs ($L_p sim 1-10 mu$m, $gamma sim$ 6000-8000 pN) estimated using Martini soft elastic network and oxDNA are in very good agreement with the all-atom MD and experimental values, while the stiff elastic network Martini reproduces order of magnitude higher values of $L_p$ and $gamma$. The high flexibility of small dsDNA is also depicted in our calculations. However, Martini models proved inadequate to capture the salt concentration effects on the mechanical properties with increasing salt molarity. OxDNA captures the salt concentration effect on small dsDNA mechanics. But it is found to be ineffective to reproduce the salt-dependent mechanical properties of DNTs. Also, unlike Martini, the time evolved PX DNA and DNT structures from the oxDNA models are comparable to the all-atom MD simulated structures. Our findings provide a route to study the mechanical properties of DNA nanostructures with increased time and length scales and has a remarkable implication in the context of DNA nanotechnology.
Integral equation theory is applied to a coarse-grained model of water to study potential of mean force between hydrophobic solutes. Theory is shown to be in good agreement with the available simulation data for methane-methane and fullerene-fullerene potential of mean force in water; the potential of mean force is also decomposed into its entropic and enthalpic contributions. Mode coupling theory is employed to compute self-diffusion coefficient of water, as well as diffusion coefficient of a dilute hydrophobic solute; good agreement with molecular dynamics simulation results is found.
Here we describe the development of the MALTS software which is a generalised tool that simulates Lorentz Transmission Electron Microscopy (LTEM) contrast of thin magnetic nanostructures. Complex magnetic nanostructures typically have multiple stable domain structures. MALTS works in conjunction with the open access micromagnetic software Object Oriented Micromagnetic Framework or MuMax. Magnetically stable trial magnetisation states of the object of interest are input into MALTS and simulated LTEM images are output. MALTS computes the magnetic and electric phases accrued by the transmitted electrons via the Aharonov-Bohm expressions. Transfer and envelope functions are used to simulate the progression of the electron wave through the microscope lenses. The final contrast image due to these effects is determined by Fourier Optics. Similar approaches have been used previously for simulations of specific cases of LTEM contrast. The novelty here is the integration with micromagnetic codes via a simple user interface enabling the computation of the contrast from any structure. The output from MALTS is in good agreement with both experimental data and published LTEM simulations. A widely-available generalized code for the analysis of Lorentz contrast addresses is a much needed step towards the use of LTEM as a standardized laboratory technique.