[Analysis of Conformational Features of Watson-Crick Duplex Fragments by Molecular Mechanics and Quantum Mechanics Methods].

It is generally accepted that the important characteristic features of the Watson-Crick duplex originate from the molecular structure of its subunits. However, it still remains to elucidate what properties of each subunit are responsible for the significant characteristic features of the DNA structure. The computations of desoxydinucleoside monophosphates complexes with Na-ions using density functional theory revealed a pivotal role of DNA conformational properties of single-chain minimal fragments in the development of unique features of the Watson-Crick duplex. We found that directionality of the sugar-phosphate backbone and the preferable ranges of its torsion angles, combined with the difference between purines and pyrimidines. in ring bases, define the dependence of three-dimensional structure of the Watson-Crick duplex on nucleotide base sequence. In this work, we extended these density functional theory computations to the minimal' fragments of DNA duplex, complementary desoxydinucleoside monophosphates complexes with Na-ions. Using several computational methods and various functionals, we performed a search for energy minima of BI-conformation for complementary desoxydinucleoside monophosphates complexes with different nucleoside sequences. Two sequences are optimized using ab initio method at the MP2/6-31++G** level of theory. The analysis of torsion angles, sugar ring puckering and mutual base positions of optimized structures demonstrates that the conformational characteristic features of complementary desoxydinucleoside monophosphates complexes with Na-ions remain within BI ranges and become closer to the corresponding characteristic features of the Watson-Crick duplex crystals. Qualitatively, the main characteristic features of each studied complementary desoxydinucleoside monophosphates complex remain invariant when different computational methods are used, although the quantitative values of some conformational parameters could vary lying within the limits typical for the corresponding family. We observe that popular functionals in density functional theory calculations lead to the overestimated distances between base pairs, while MP2 computations and the newer complex functionals produce the structures that have too close atom-atom contacts. A detailed study of some complementary desoxydinucleoside monophosphate complexes with Na-ions highlights the existence of several energy minima corresponding to BI-conformations, in other words, the complexity of the relief pattern of the potential energy surface of complementary desoxydinucleoside monophosphate complexes. This accounts for variability of conformational parameters of duplex fragments with the same base sequence. Popular molecular mechanics force fields AMBER and CHARMM reproduce most of the conformational characteristics of desoxydinucleoside monophosphates and their complementary complexes with Na-ions but fail to reproduce some details of the dependence of the Watson-Crick duplex conformation on the nucleotide sequence.

[Interactions of DNA bases with individual water molecules. Molecular mechanics and quantum mechanics computation results vs. experimental data].

To elucidate details of the DNA-water interactions we performed the calculations and systemaitic search for minima of interaction energy of the systems consisting of one of DNA bases and one or two water molecules. The results of calculations using two force fields of molecular mechanics (MM) and correlated ab initio method MP2/6-31G(d, p) of quantum mechanics (QM) have been compared with one another and with experimental data. The calculations demonstrated a qualitative agreement between geometry characteristics of the most of local energy minima obtained via different methods. The deepest minima revealed by MM and QM methods correspond to water molecule position between two neighbor hydrophilic centers of the base and to the formation by water molecule of hydrogen bonds with them. Nevertheless, the relative depth of some minima and peculiarities of mutual water-base positions in' these minima depend on the method used. The analysis revealed insignificance of some differences in the results of calculations performed via different methods and the importance of other ones for the description of DNA hydration. The calculations via MM methods enable us to reproduce quantitatively all the experimental data on the enthalpies of complex formation of single water molecule with the set of mono-, di-, and trimethylated bases, as well as on water molecule locations near base hydrophilic atoms in the crystals of DNA duplex fragments, while some of these data cannot be rationalized by QM calculations.

Electron correlated ab initio study of amino group flexibility for improvement of molecular mechanics simulations on nucleic acid conformations and interactions.

High level ab initio studies demonstrate substantial conformational flexibility of amino groups of nucleic acid bases. This flexibility is important for biological functions of DNA. Existing force field models of molecular mechanics do not describe this phenomenon due to a lack of quantitative experimental data necessary for an adjustment of empirical parameters. We have performed extensive calculations of nucleic acid bases at the MP2/6-31G(d,p) level of ab initio theory for broad set of amino group configurations. Two-dimensional maps of energy and geometrical characteristics as functions of two amino hydrogen torsions have been constructed. We approximate the maps by polynomial expressions, which can be used in molecular mechanics calculations. Detailed considerations of these maps enable us to propose a method for determination of numerical coefficients in the developed formulae using restricted sets of points obtained via higher-level calculations.

[Possible configuration of dimers of Gua-Cyt bases. Computation by molecular mechanic models and density functional theory]. / Vozmozhnye konfiguratsii dimerov osnovanii DNK Gua-Cyt. Raschety metodami mekhaniki i teorii funktsional plotnosti.

The energies of interactions between guanine and cytosine in various mutual positions were calculated by the methods of molecular mechanics with refined atom-atom potential functions and the quantum mechanics theory of density functional. Both methods indicate three types of mutual positions of bases in local energy minima. These types correspond to (1) nearly coplanar base positions with intermolecular hydrogen bond formation (base pairing); (2) arrangements of two bases in nearly parallel planes one above another (base stacking); and (3) nearly perpendicular positions of base planes. According to the calculations, the global energy minimum corresponds to the Watson-Crick base pair with three hydrogen bonds. A specific feature of the pair is a transition from many positions of type (2) to positions of type (1) without any energy barrier. This feature is revealed by both methods. Another special feature of this pair is a deviation, for most of mutual base positions, of the amine group atoms from the ring plane, the deviation being more pronounced for Gua. These features are important for understanding the conformational behavior of DNA fragments and the RNA structure.

[Interaction of caffeine with basic nucleic acids my a molecular mechanic method]. / Issledovanie vzaimodeistviia kofeina s osnovaniiami nukleinovykh kislot metodom molekuliarnoi mekhaniki.

To understand some aspects of the biological action of caffeine (CAF), the interaction energies for various mutual positions of CAF and DNA bases or basepairs were calculated. Three types of mutual CAF-base (CAF-basepair) arrangements corresponding to the minima of interaction energy were revealed. One type of minima correspond to the stacking arrangement of molecules. This type is important for interactions of CAF with DNA monomers and single-stranded fragments. The other two types of minima correspond to the formation of intermolecular hydrogen bonds. Some of these minima may occur during the interaction of CAF with the double helix. One of these types corresponds to the nearly in-plane position of molecules. The other type of minima correspond to the nearly perpendicular arrangement of molecule planes. The minima of the last type are supposed to be the most important for the interaction of CAF with the DNA duplex, and interaction energies for this type of minima have the most negative values.

[Interactions between nucleic acid bases. New parameters of potential functions and energy minima]. / Vzaimodeistviia mezhdu osnovaniiami nukleinovykh kislot. Novye parametry potentsialov i novye minimumy énergii.

The parameters of atom-atom potential functions suggested by one of the authors in 1979-1986 were slightly changed. The changes were performed to achieve a better agreement with experimental data of interaction energy values in global minima and hydrogen bond lengths. These changes resulted in better accord with experimental values of distances between the layers in DNA monomer crystals and between the base pairs in oligonucleotide duplexes. The refined potential functions were used to calculate the energy of interactions between nucleic acid bases in various mutual positions. The calculations revealed a few types of mutual base arrangements in minima of interaction energy for each pairwise base combination. A new type of minima was found, which correspond to a nearly perpendicular arrangement of base rings and the formation of the intermolecular hydrogen bond.

A study of the hydration of deoxydinucleoside monophosphates containing thymine, uracil and its 5-halogen derivatives: Monte Carlo simulation.

An extensive Monte Carlo simulation of hydration of various conformations of the dinucleoside monophosphates (DNP), containing thymine, uracil and its 5-halogen derivatives has been performed. An anti-anti conformation is the most energetically stable one for each of the DNPs. In the majority of cases the energy preference is determined by water-water interaction. For other dimers conformational energy is the most important factor, or both the factors are of nearly equal importance. The introduction of the methyl group into the 5-position of uracil ring most noticeably influences the conformational energy and leads to the decrease of its stabilizing contribution to the total interaction energy. The introduction of halogen atoms increases the relative content of anti-syn and syn-anti conformations of DNPs as compared to the parent ones due to the formation of an energetically more favorable water structure around these conformations. A correlation is observed between the Monte Carlo results for the halogenated DNPs and their experimental photoproduct distribution. The data obtained demonstrates a sequence dependence in the photochemistry of the halogenated dinucleoside monophosphates.

Monte Carlo simulation of DNA fragment hydration in the presence of alkaline cations using novel atom-atom potential functions.

The set of atom-atom potential functions specially adjusted to simulation of nucleic acid fragment hydration (Poltev, Grokhlina & Malenkov, J. Biomol. Struct. Dyn. 2, 413, 1984) is extended by including alkaline cation interactions. The choice of new potential functions was realized using experimental data on crystal hydrates of nucleotides and related compounds as well as thermodynamic data on ion solutions. The extended set of potential functions allows to reproduce many features of interactions between alkaline cations and nucleic acid fragments in water solutions. The sites of preferential cation localization near bases and phosphate groups were obtained and examined. The potential functions reproduce the dissociation tendency of cation-phosphate group and cation-base complexes in aqueous medium. Pathways of cation dissociations from nucleic acid components have been studied, and metastable water-bridged positions of cations near bases and phosphate group have been revealed.

Sequence-dependent binding of metal ion to DNA oligomeres. A comparison of molecular electrostatic potentials with NMR data.

Experimentally observed sequence-selective binding of metal ion to DNA oligonucleotides have been compared with variations of electrostatic potential (EP) along the helix. Calculations of EP have been performed for three atomic models of the oligonucleotide duplex [d(CGCGAATTCGCG)2] using several variants of EP calculations, including a solution of non-linear Poisson-Boltzmann equation (NPBE). N7 atom of guanine adjacent to adenine base was identified as a region with the most negative electrostatic potential in the major groove. The EP value for the Me ion binding site surpasses the value for N7 of other guanines by 10-26% depending on particular duplex conformation. Qualitatively, the sequence dependent variations of EP near guanine N7 atoms are in agreement with the sequence-selective behavior of Mn(II) and Zn(II) ions as revealed by NMR experiments. But the difference in EP between the two most negative regions near guanine N7 atoms does not exceed 1.25 kT/e. Simple model suggests that metal ions are capable to form ion-hydrate complexes with G-Pu steps of DNA duplex. These complexes are formed via one Me...G and five Me...water coordination bonds with water molecules hydrogen bonded to two adjacent purine bases in the same chain. We suppose that such a stereospecific structural possibility is the main factor which control the sequence-selectivity in the metal ion binding. A combination of both mechanisms allows to explain sequence specific Mn(II) and Zn(II) binding to a set of oligonucleotides.

Protein-nucleic acid recognition: simulation of base and "model" amino acids complexes in DMSO by the Monte Carlo method.

A computer simulation of guanine (G), cytosine (C), the G-C base pair, protonated C (CH+), acetic acid in neutral (AcOH) and deprotonated (AcO-) forms, G-AcO-, C-AcOH, and CH(+)-AcO- complexes, solvated in DMSO was carried out by the Monte Carlo method. It is shown that the G-C base pair formation in DMSO is energetically favorable. The G-AcO- complex formation is comparable with the formation of G-C base pair in energetically favorability. In this case the acetate anion can replace C in the G-C base pair. The formation of the C-AcOH complex is much less favorable than the formation of the G-C pair. However proton transfer from AcOH to C leads to the formation of the CH(+)-AcO- complex, which is the most favorable of all complexes studied. Here the acetic acid can replace G in a G-C base pair. The formation of G-AcO- and CH(+)-AcO- specific complexes detected in DMSO with the help of experiment and theory is a competitive process with respect to the formation of G-C base pairs, and can be considered the primary step in the real mechanism of protein-nucleic acid recognition.

The study of the stability of Watson-Crick nucleic acid base pairs in water and dimethyl sulfoxide: computer simulation by the Monte Carlo method.

An extensive computer simulation of nucleic acid bases and Watson-Crick base pairs in a water cluster and DMSO cluster is performed by the Monte Carlo method. It is demonstrated that the unfavorable energetics of pair formation in a water cluster is determined by the significant destabilizing contribution of solvent to the energy of complex formation. It is shown that the formation of coplanar base pairs in a DMSO cluster is favorable. The DMSO cluster stabilizes A-U and A-T base pairs and the insignificant destabilization of the G-C base pair by a DMSO cluster is much less than the stabilization which occurs due to the attraction between bases.

The study of possible A and B conformations of alternating DNA using a new program for conformational analysis of duplexes (CONAN).

A new program, CONAN has been designed for CONformational ANalysis of oligonucleotide duplexes with natural and modified bases. It allows to model both regular DNA fragments with different types of symmetry and irregular ones including bends, junctions, mismatched pairs and base lesions. Computations and minimization of the energy are performed in a space of internal structural variables chosen to build start structure easier and conveniently analyze the results obtained. These internal structural variables determine mutual base-base and base-sugar arrangement and sugar puckering. The analytical closure procedure is applied both to sugar rings and to backbone fragments between adjacent sugars. For more effective energy minimization, analytical gradient is calculated. The CONAN was applied to the search for low-energy conformations of poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC) duplexes. Extended regions of low-energy A and B conformations are revealed and characterized. These regions contain structures with different relative values of helical twist, tau, for pur-pyr and pyr-pur steps, namely, conformations with tau (pur-pyr) > tau (pyr-pur) and with tau (pur-pyr) < tau (pyr-pur). Two types of sugar puckering were found for B-form low-energy conformations, the first type with all C2'-endo sugar residues and the second one - with C2'-endo purines and O1'-endo pyrimidines. The calculated conformations are compared with X-ray diffraction data for crystals and fibers and NMR data for solution.

Modeling DNA hydration: comparison of calculated and experimental hydration properties of nuclic acid bases.

Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensen's TIP3P water model are widely used in the AMBER program, and (iii) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9-methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methylguanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration.