From 20.309 Course Miki
20.309: Biological Instrumentation and Measurement
DNA in solution
 Consider a solution containing equal quantities of complementary single stranded DNA (ssDNA) oligonucleotides and .

 Complementary ssDNA strands bond to form double stranded DNA (dsDNA). The reaction is governed by the equation

 The forward reaction in which two ssDNA oligos combine to form dsDNA is called annealing. The reverse process is called thermal denaturation or melting.

 At low temperatures, dsDNA is favored. As the temperature rises, dsDNA increasingly separates into its component ssDNA oligos. (Think about why with respect to enthalpic and entropic considerations.)

 The melting temperature, , is defined to be the point where half of the dsDNA is denatured.

 Short sequences of about 1040 base pairs (such as those used in the DNA Melting lab) tend to denature all at once. Longer sequences may melt in segments.

 Less energy is required to split the double hydrogen bond of AT pairs than the triple bond of GC pairs. Thus, AT rich sequences tend to melt at lower temperatures than GC rich ones.^{[1]}

Several web tools are available to predict the melting temprature. (See, for example, DINA Melt or Oligocalc.)
Fundamental equilibrium relationships
 The concentrations of the reaction products are related by the equilibrium constant:

 The value of is a function of temperature. We can equate the fundamental definition of the standard free energy change with its relationship to the equilibrium constant in solution:


 where
 is the standard change in free energy
 is the standard enthalpy change
 is the temperature
 is the standard entropy change
 is the gas constant
 Solving for :


Note that the above equation can be differentiated with respect to temperature to yield the (perhaps once!) familiar van't Hoff equation.
Theoretical relation of dsDNA fraction and thermodynamic parameters
 In the lab, the fraction of dsDNA will be measured with a fluorescent dye that preferentially binds to dsDNA. As such, it will be useful to derive an equation that relates the fraction of dsDNA to temperature and the thermodynamic parameters.

 Let represent the concentration of either single stranded oligonucleotide: .

 Similarly, let be the concentration of double stranded DNA:

 is the total concentration of DNA strands.

 Let be the fraction of total DNA that is double stranded


 Therefore,

 Now we can solve for in terms of and :


 At the melting point, by definition and thus .

 Substituting for K_{eq} from equation 1:


 Taking the log of both sides and solving for :


We now know temperature as a function of dsDNA fraction for a given total DNA strand concentration and DNA identity (therefore ΔH, etc.).
Simulating DNA melting for tractable calculation of dsDNA fraction
 For simulating DNA melting experiments, it will be convenient to have an expression for in terms of . Unfortunately, this gets pretty yucky. On the bright side, Matlab and Python are good at calculating yuck.

 Taking the log of both sides of equation 2 (after resubstiting in equation 1 for simplicity) and using the quadratic formula (eliminating the nonphysical root):


 Substituting from equation 1 gives the desired result.


See the pages DNA Melting Part 1: Simulating DNA Melting  Basics. And if you're interested in a Python implementation see Python:Simulating DNA Melting
References
 ↑ Breslauer et al., Predicting DNA duplex stability from the base sequence PNAS 83: 3746, 1986