We perform a spatially resolved simulation study of an AND gate based on DNA strand displacement using several lengths of the toehold and the adjacent domains. DNA strands are modelled using a coarse-grained dynamic bonding model {[}C. Svaneborg, Com
p. Phys. Comm. 183, 1793 (2012){]}. We observe a complex transition path from the initial state to the final state of the AND gate. This path is strongly influenced by non-ideal effects due to transient bubbles revealing undesired toeholds and thermal melting of whole strands. We have also characterized the bound and unbound kinetics of single strands, and in particular the kinetics of the total AND operation and the three distinct distinct DNA transitions that it is based on. We observe a exponential kinetic dependence on the toehold length of the competitive displacement operation, but that the gate operation time is only weakly dependent on both the toehold and adjacent domain length. Our gate displays excellent logical fidelity in three input states, and quite poor fidelity in the fourth input state. This illustrates how non-ideality can have very selective effects on fidelity. Simulations and detailed analysis such as those presented here provide molecular insights into strand displacement computation, that can be also be expected in chemical implementations.
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 s
ingle 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.
We have extended the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) to support directional bonds and dynamic bonding. The framework supports stochastic formation of new bonds, breakage of existing bonds, and conversion between bon
d types. Bond formation can be controlled to limit the maximal functionality of a bead with respect to various bond types. Concomitant with the bond dynamics, angular and dihedral interactions are dynamically introduced between newly connected triplets and quartets of beads, where the interaction type is determined from the local pattern of bead and bond types. When breaking bonds, all angular and dihedral interactions involving broken bonds are removed. The framework allows chemical reactions to be modeled, and use it to simulate a simplistic, coarse-grained DNA model. The resulting DNA dynamics illustrate the power of the present framework.
Recently we developed a formalism for the scattering from linear and acyclic branched structures build of mutually non-interacting sub-units.{[}C. Svaneborg and J. S. Pedersen, J. Chem. Phys. 136, 104105 (2012){]} We assumed each sub-unit has referen
ce points associated with it. These are well defined positions where sub-units can be linked together. In the present paper, we generalize the formalism to the case where each reference point can represent a distribution of potential link positions. We also present a generalized diagrammatic representation of the formalism. Scattering expressions required to model rods, polymers, loops, flat circular disks, rigid spheres and cylinders are derived. and we use them to illustrate the formalism by deriving the generic scattering expression for micelles and bottle brush structures and show how the scattering is affected by different choices of potential link positions.
We present a formalism for the scattering of an arbitrary linear or acyclic branched structure build by joining mutually non-interacting arbitrary functional sub-units. The formalism consists of three equations expressing the structural scattering in
terms of three equations expressing the sub-unit scattering. The structural scattering expressions allows a composite structures to be used as sub-units within the formalism itself. This allows the scattering expressions for complex hierarchical structures to be derived with great ease. The formalism is furthermore generic in the sense that the scattering due to structural connectivity is completely decoupled from internal structure of the sub-units. This allows sub-units to be replaced by more complex structures. We illustrate the physical interpretation of the formalism diagrammatically. By applying a self-consistency requirement we derive the pair distributions of an ideal flexible polymer sub-unit. We illustrate the formalism by deriving generic scattering expressions for branched structures such as stars, pom-poms, bottle-brushes, and dendrimers build out of asymmetric two-functional sub-units.
We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, $G(t)$, into the plateau regime
for chains with $Z=40$ entanglements and into the terminal relaxation regime for $Z=10$. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter -free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.
