Simulation of conformational possibilities of DNA via calculation of nonbonded interactions of complementary dinucleoside phosphate complexes

Using classical potential functions, we carried out potential-energy calculations on the complementary deoxydinucleoside phosphate complexes dApdA:dUpdU, dUpdA:dUpdA, and dApdU:dApdU. All dihedral and bond angles, except those of the nitrogen bases, were varied. The resulting minimum-energy conformations of the complexes are close to DNA A- and B-family conformations, with a typical arrangement of the nitrogen bases. The dihedral and bond angles of one of the molecules forming the complex can thereby differ by several degrees from those of the other molecule. For different base sequences, some dihedral and bond angles may vary over a range of several degrees without appreciably changing the total energy of the complex. Some low-energy conformations of the complexes corresponding to other regions of the conformational space are also found. The biological consequences of possible changes in dihedral and bond angles, occurring on interaction with other molecules, are discussed.