Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userVarying the resolution of the Rouse model on temporal and spatial scales: application to multiscale modelling of DNA dynamics
70
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A multi-resolution bead-spring model for polymer dynamics is developed as a generalization of the Rouse model. A polymer chain is described using beads of variable sizes connected by springs with variable spring constants. A numerical scheme which can use different timesteps to advance the positions of different beads is presented and analyzed. The position of a particular bead is only updated at integer multiples of the timesteps associated with its connecting springs. This approach extends the Rouse model to a multiscale model on both spatial and temporal scales, allowing simulations of localized regions of a polymer chain with high spatial and temporal resolution, while using a coarser modelling approach to describe the rest of the polymer chain. A method for changing the model resolution on-the-fly is developed using the Metropolis-Hastings algorithm. It is shown that this approach maintains key statistics of the end-to-end distance and diffusion of the polymer filament and makes computational savings when applied to a model for the binding of a protein to the DNA filament.
rate research
Read More
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.
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.
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.
Conformational change of a DNA molecule is frequently observed in multiple biological processes and has been modelled using a chain of strongly coupled oscillators with a nonlinear bistable potential. While the mechanism and properties of conformational change in the model have been investigated and several reduced order models developed, the conformational dynamics as a function of the length of the oscillator chain is relatively less clear. To address this, we used a modified Lindstedt-Poincare method and numerical computations. We calculate a perturbation expansion of the frequency of the models nonzero modes, finding that approximating these modes with their unperturbed dynamics, as in a previous reduced order model, may not hold when the length of the DNA model increases. We investigate the conformational change to local perturbation in models of varying lengths, finding that for chosen input and parameters, there are two regions of DNA length in the model, first where the minimum energy required to undergo the conformational change increases with DNA length; and second, where it is almost independent of the length of the DNA model. We analyze the conformational change in these models by adding randomness to the local perturbation, finding that the tendency of the system to remain in a stable conformation against random perturbation decreases with an increase in the DNA length. These results should help to understand the role of the length of a DNA molecule in influencing its conformational dynamics.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
