A molecular orbital study of stability and the conformation of double-stranded DNA-like polymers

An optimized potential function for the calculation of nucleic acid interaction energies I. base stacking

An optimized potential function for base-stacking interactions is constructed. Stacking energies between the complementary pairs of a dimer are calculated as a function of the rotational angle and separation distance. Using several different sets of atomic charges, the electrostatic component in the monopole-monopole approximation (MMA) is compared to the more refined segmented multipole-multipole representation (SMMA); the general features of the stacking minima are found to be correctly reproduced with IEHT or CNDO atomic charges. The electrostatic component is observed to control the location of stacking minima.The MMA, in general, is not a reliable approximation of the SMMA in regions away from minima; however, the MMA is reliable in predicting the location and nature of stacking minima.The attractive part of the Lennard-Jones 6-12 potential is compared to and parameterized against the expressions for the second-order interaction terms composed of multipole-bond polarizability for the polarization energy and transition-dipole bond polariz abilities for approximation of the dispersion energy. The repulsive part of the Lennard-Jones potential is compared to a Kitaygorodski-type repulsive function; changing the exponent from its usual value of 12 to 11.7 gives significantly better agreement with the more refined repulsive function.Stacking minima calculated with the optimized potential method are compared with various perturbation-type treatments. The optimized potential method yields results that compare as well with melting data as do any of the more recent and expensive perturbation methods.

Efficient and accurate calculations on the electronic structure of B-type poly(dG).poly(dC) DNA by elongation method: first step toward the understanding of the biological properties of aperiodic DNA

Elongation method was applied to determine the electronic structures of B-type poly(dG).poly(dC) DNA at the ab initio molecular orbital level as a first step toward the calculation of aperiodic DNA. The discrepancy in total energy between the elongation method and a conventional calculation was negligibly small in the order of 10(-8) hartreeat. for 14 G-C base pair model. The local density of states for 10 G-C base pair model estimated by the elongation method well reproduced the results by the conventional calculation. It was found that the band gap of the whole system is mainly due to the energy difference between the valence band of guanine and the conduction band of cytosine. Moreover, the electron transfer path through stacking G-C base pairs rather than sugar-phosphate backbones has been confirmed by the authors' calculations.

On a simple model of low-frequency vibrations in DNA macromolecules

A simple quasi-continuity model of low-frequency dynamics of DNA macromolecules is presented. The model is based on the phenomenological theory by Volkov and Kosevich (J. Biomol. Struct. Dyn. 8, 1069 (1991)). We propose improvements to their theory by correcting the model energy and corresponding equations of motion, recalculating the model DNA parameters, and estimating the effects of hydration and counterion binding on the masses of DNA subunits. Comparing the calculated low-frequency vibration spectrum and the results of classical Raman scattering experiments from the literature, the values for the model force constants are determined. These values differ significantly from the estimations made by Volkov and Kosevich. A good quantitative agreement with experiment can be obtained both for the case when the "25 cm-1 mode" is an external one, and when it is considered as an intrahelical mode. Problems connected with the qualitative assignments of the calculated and experimentally observed modes are discussed.

The electrostatic contribution to DNA base-stacking interactions

Base-stacking and phosphate-phosphate interactions in B-DNA are studied using the finite difference Poisson-Boltzmann equation. Interaction energies and dielectric constants are calculated and compared to the predictions of simple dielectric models. No extant simple dielectric model adequately describes phosphate-phosphate interactions. Electrostatic effects contribute negligibly to the sequence and conformational dependence of base-stacking interactions. Electrostatic base-stacking interactions can be adequately modeled using the Hingerty screening function. The repulsive and dispersive Lennard-Jones interactions dominate the dependence of the stacking interactions on roll, tilt, twist, and propellor. The Lennard-Jones stacking energy in ideal B-DNA is found to be essentially independent of sequence.

Informational parameters of nucleic acid and molecular evolution

From the point of view of information theory, a statistical analysis of 2000 nucleic acid sequences (732 coding regions and 1177 non-coding regions) is given. The sequences are grouped into 20 categories. The probability-order-difference (POD) matrix is defined which is used to analyse the evolutionary distance of any two categories of sequences. The informational parameters D1, D2 and X = (1 + D1/D2)-1 and F are calculated for each sequence and averaged in each category. The statistical dependence of these parameters on molecular evolution is discussed. It is found that [X] is a good statistical quantity which describes the vocabulary compositions as well as the grammatical constructions of the genetic language. From the statistical analysis it is shown that [X] may play an important role in investigating the evolutionary level of nucleic acid molecules.

Origins of the specificity in the intercalation of ethidium into nucleic acids. A theoretical analysis

The sequence preferences observed in the intercalative binding of ethidium to dinucleoside phosphates have been theoretically examined. This specificity is for pyrimidine (3'--5') purine sequences as compared to their purine (3'--5') pyrimidine sequence isomers. It is shown that the stacking energies between the ethidium cation and the base pairs are fairly constant for all combinations of bases at the intercalation site. In contrast, the energy of unwinding the double helix to assume the geometry of the intercalation complex shows substantial sequence differences. Thus, the specificity observed is more readily explained in terms of these conformation energy changes than by preferential stacking interactions. These results imply that there may be a large class of intercalating molecules which exhibit similar pyrimidine (3'--5') purine sequence specificity.

Kinked DNA. Energetics and conditions favoring its formation

Theoretical calculations demonstrate that the formation of kinks in DNA in the manner proposed by Crick and Klug produces a conformation which corresponds to a local energy minimum when the phosphate backbone is charged. This conformation becomes the global minimum under certain conditions which neutralize the anionic backbone oxygens. A sequence preference for kink formation is predicted. The possible role of kinking in experimentally observed phenomena and experimental means of demonstrating the existence of kinks are discussed.