Skip to main content
Article thumbnail
Location of Repository

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

By Boian S. Alexandrov, Vladimir Gelev, Yevgeniya Monisova, Ludmil B. Alexandrov, Alan R. Bishop, Kim Ø. Rasmussen and Anny Usheva


No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats

Topics: Computational Biology
Publisher: Oxford University Press
OAI identifier:
Provided by: PubMed Central

Suggested articles


  1. (1989). A computer program for choosing optimal oligonudeotides for filter hybridization, sequencing and in vitro amplification of DNA.
  2. (1972). A quantum mechanical muscle model.
  3. (1995). A second generation force field for the simulation of proteins, nucleic acids, and organic molecules.
  4. (1998). A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics.
  5. (2000). An initiation element in the yeast CUP1 promoter is recognized by RNA polymerase II in the absence of TATA box-binding protein if the DNA is negatively supercoiled.
  6. (2006). Base-stacking and base-pairing contributions into thermal stability of the DNA double helix.
  7. (2004). Bubble nucleation and cooperativity in DNA melting.
  8. (2006). Bubble statistics and dynamics in doublestranded DNA.
  9. (2004). DNA dynamically directs its own transcription initiation.
  10. (1993). Entropy-driven DNA denaturation.
  11. (1998). Experimental tests of the Peyrard-Bishop model applied to the melting of very short DNA chains.
  12. (1995). Extracting thermodynamic data from equilibrium melting curves for oligonucleotide order-disorder transitions.
  13. (1992). How many numbers are required to specify sequence-dependent properties of polynucleotides?
  14. (1996). Improved nearest-neighbor parameters for predicting DNA duplex stability.
  15. (1989). Molecular Cloning: A Laboratory Manual,
  16. (2008). Nature and magnitude of aromatic stacking of nucleic acid bases.
  17. (2006). Nature of base stacking: reference quantumchemical stacking energies in ten unique B-DNA base-pair steps.
  18. (1980). Nature of the open state in long polynucleotide double helices: possibility of soliton excitations.
  19. (2007). OligoCalc: an online oligonucleotide properties calculator.
  20. (1966). Phase transitions in onedimension and the helix-coil transition in polyamino acids.
  21. (1970). Physicochemical studies on polydeoxyribonucleotides containing defined repeating nucleotide sequences.
  22. (2009). Pre-melting dynamics of DNA and its relation to specific functions.
  23. (1986). Predicting DNA duplex stability from the base sequence.
  24. (2006). Probing the mechanical unzipping of DNA.
  25. (2008). Profiling the thermodynamic softness of adenoviral promoters.
  26. (2004). Sequence-specific thermal fluctuations identify start sites for DNA transcription.
  27. (1999). Statistical mechanical simulation of polymeric DNA melting with MELTSIM.
  28. (1989). Statistical mechanics of a nonlinear model for DNA denaturation.
  29. (2004). Temperature-dependent signatures of coherent vibrational openings in DNA.
  30. (2001). The sequence of the human genome.
  31. (1973). The theory of contraction of proteins under their excitation.
  32. (1996). YY1 transcriptional initiator: protein interactions and association with a DNA site containing unpaired strands.

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.