No Arabic abstract
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, Comp. 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.
Cytosine methylation has been found to play a crucial role in various biological processes, including a number of human diseases. The detection of this small modification remains challenging. In this work, we computationally explore the possibility of detecting methylated DNA strands through direct electrical conductance measurements. Using density functional theory and the Landauer-Buttiker method, we study the electronic properties and charge transport through an eight base-pair methylated DNA strand and its native counterpart. We first analyze the effect of cytosine methylation on the tight-binding parameters of two DNA strands and then model the transmission of the electrons and conductance through the strands both with and without decoherence. We find that the main difference of the tight-binding parameters between the native DNA and the methylated DNA lies in the on-site energies of (methylated) cytosine bases. The intra- and inter- strand hopping integrals between two nearest neighboring guanine base and (methylated) cytosine base also change with the addition of the methyl groups. Our calculations show that in the phase-coherent limit, the transmission of the methylated strand is close to the native strand when the energy is nearby the highest occupied molecular orbital level and larger than the native strand by 5 times in the bandgap. The trend in transmission also holds in the presence of the decoherence with the same rate. The lower conductance for the methylated strand in the experiment is suggested to be caused by the more stable structure due to the introduction of the methyl groups. We also study the role of the exchangecorrelation functional and the effect of contact coupling by choosing coupling strengths ranging from weak to strong coupling limit.
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.
When DNA molecules are heated they denature. This occurs locally so that loops of molten single DNA strands form, connected by intact double-stranded DNA pieces. The properties of this melting transition have been intensively investigated. Recently there has been a surge of interest in this question, caused by experiments determining the properties of partially bound DNA confined to nanochannels. But how does such confinement affect the melting transition? To answer this question we introduce, and solve a model predicting how confinement affects the melting transition for a simple model system by first disregarding the effect of self-avoidance. We find that the transition is smoother for narrower channels. By means of Monte-Carlo simulations we then show that a model incorporating self-avoidance shows qualitatively the same behaviour and that the effect of confinement is stronger than in the ideal case.
In this work, we model the zero-bias conductance for the four different DNA strands that were used in conductance measurement experiment [A. K. Mahapatro, K. J. Jeong, G. U. Lee, and D. B. Janes, Nanotechnology 18, 195202 (2007)]. Our approach consists of three elements: (i) ab initio calculations of DNA, (ii) Greens function approach for transport calculations, and (iii) the use of two parameters to determine the decoherence rates. We first study the role of the backbone. We find that the backbone can alter the coherent transmission significantly at some energy points by interacting with the bases, though the overall shape of the transmission stays similar for the two cases. More importantly, we find that the coherent electrical conductance is tremendously smaller than what the experiments measure. We consider DNA strands under a variety of different experimental conditions and show that even in the most ideal cases, the calculated coherent conductance is much smaller than the experimental conductance. To understand the reasons for this, we carefully look at the effect of decoherence. By including decoherence, we show that our model can rationalize the measured conductance of the four strands, both qualitatively and quantitatively. We find that the effect of decoherence on G:C base pairs is crucial in getting agreement with the experiments. However, the decoherence on G:C base pairs alone does not explain the experimental conductance in strands containing a number of A:T base pairs. Including decoherence on A:T base pairs is also essential. By fitting the experimental trends and magnitudes in the conductance of the four different DNA molecules, we estimate for the first time that the deocherence rate is 6 meV for G:C and 1.5 meV for A:T base pairs.
The threading of a polymer chain through a small pore is a classic problem in polymer dynamics and underlies nanopore sensing technology. However important experimental aspects of the polymer motion in a solid-state nanopore, such as an accurate measurement of the velocity variation during translocation, have remained elusive. In this work we analysed the translocation through conical quartz nanopores of a 7 kbp DNA double-strand labelled with six markers equally spaced along its contour. These markers, constructed from DNA hairpins, give direct experimental access to the translocation dynamics. On average we measure a 5% reduction in velocity during the translocation. We also find a striking correlation in velocity fluctuations with a decay constant of 100s of {mu}s. These results shed light on hitherto unresolved problems in the dynamics of DNA translocation and provide guidance for experiments seeking to determine positional information along a DNA strand.