Generalized Poland-Scheraga model for DNA hybridization.
The Poland-Scheraga (PS) model for the helix-coil transition of DNA considers the statistical mechanics of the binding (or hybridization) of two complementary strands of DNA of equal length, with the restriction that only bases with the same index along the strands are allowed to bind. In this article, we extend this model by relaxing these constraints: We propose a generalization of the PS model that allows for the binding of two strands of unequal lengths N1 and N2 with unrelated sequences. We study in particular (i) the effect of mismatches on the hybridization of complementary strands, (ii) the hybridization of noncomplementary strands (as resulting from point mutations) of unequal lengths N1 and N2. The use of a Fixman-Freire scheme scales down the computational complexity of our algorithm from O(N1(2)N2(2) to O(N1N2). The simulation of complementary strands of a few kilo base pairs yields results almost identical to the PS model. For short strands of equal or unequal lengths, the binding displays a strong sensitivity to mutations. This model may be relevant to the experimental protocol in DNA microarrays, and more generally to the molecular recognition of DNA fragments. It also provides a physical implementation of sequence alignments.
Garel
Thomas
T
Service de Physique Théorique, CEA/DSM/SPhT, Unité de recherche associée au CNRS, 91191 Gif-sur-Yvette cedex, France. garel@spht.saclay.cea.fr
Orland
Henri
H
D000465:Algorithms; D000818:Animals; D020137:Base Pair Mismatch; D020029:Base Pairing; D001483:Base Sequence; D002875:Chromosomes / Q000737:chemistry; D003198:Computer Simulation; D004247:DNA / Q000737:chemistry* / Q000235:genetics / Q000378:metabolism; D004331:Drosophila melanogaster / Q000235:genetics; D008956:Models, Chemical; D009690:Nucleic Acid Conformation; D009693:Nucleic Acid Hybridization / Q000379:methods*; D017354:Point Mutation; D013056:Spectrophotometry, Ultraviolet; D013696:Temperature; D013816:Thermodynamics
D004247:DNA