Theory of the influence of oligonucleotide chain conformation on double helix stability

Publisher Website
Google Scholar
Abstract
We consider theoretical aspects of reactions that form base pairs in a double helix. The equilibrium constant for such reactions depends on the probability of finding the two bases in the correct orientation for pairing. This probability can be expressed in terms of the spatial and angular distribution of one micleotide around the other. In this paper we use Monte‐Carlo techniques to calculate the distribution of distances between chosen phosphates in nonhclical oligonucleotide backbones, using crystallographic data for bond lengths and angles, and a screened Coulomb potential for phosphate–phosphate interactions. The model chosen is one that predicts correctly the observed dimensions of an unperturbed polynucleotide chain. Knowledge of distance distribution functions permits calculation of the dependence on loop size of the probability of closing a single backbone strand into a hairpin helix. Our results agree roughly, although not exactly, with the semiempirical ring‐weighting functions determined by Schefller. Elson, and Baldwin. Further results are a comparison of intramolecular and bimolecular helix nucleation equilibrium constants and a calculation of the stacking free energy in a double helix.