Topological distributions and the torsional rigidity of DNA. A Monte Carlo study of DNA circles.

Distributions of the linking number of circular DNA molecules, defined as the sum of twist and the writhing number, are obtained by Monte Carlo simulations of small, randomly closed DNA circles. We estimate the relative contributions of fluctuations in twist and writhe to the linking number distribution, as functions of DNA size. Published experimental data on topoisomer distributions in circular DNA molecules are interpreted to estimate the torsional rigidity of DNA in solution. We show that ignoring the writhe component of the linking number distribution, even for DNA circles as small as 250 base-pairs, leads to an underestimate for the torsional stiffness of the double helix. The value of the torsional modulus obtained from this analysis, C = 3.4 X 10(-19) erg cm, is from 10 to 40% larger than that estimated by others and more than twice as large as the values obtained from fluorescence depolarization or other time-resolved spectroscopic measurements. We also develop further the theoretical treatment of ring closure probabilities for DNA described in the previous article. It is shown that the torsional part of the ring closure probability, phi 0,1 (tau 0) is a periodic function of DNA length that contributes strongly to the ring closure probability for short chains but makes negligible contributions for chains over 1000 base-pairs in length.