CALCULATION OF TEMPERATURE-DEPENDENCE OF INTERBASE BREATHING MOTION OF A GUANINE-CYTOSINE DNA DOUBLE HELIX WITH ADENINE-THYMINE INSERT

The modified self-consistent phonon approximation is generalized to calculate the critical temperature of a DNA double helix with block inserts of different base pair sequences. An iterative method based entirely on the Green's function is developed to compute efficiently the self-consistent anharmonic force constants and thermal mean-squared vibrational amplitudes in hydrogen bonds. The calculation is carried out for a long guanine-cytosine (G-C) type helix with an insert of d(A)4-d(T)4. The melting-associated behavior is predicted to initiate in major groove bonds in the inserted adenine-thymine (A-T) base pairs at 349 K. This temperature is above the mean-field melting temperature of poly(dA)-poly(dT) due to stabilization of the A-T bonds by the G-C helix and below the mean-field melting temperature of poly(dG)-poly(dC), where dA refers to repeating adenine bases on one strand and dT the repeating thymine bases on the other strand.