Physical properties of over 8000 intramolecular hydrogen bonds (iHBs), including 2901 ones of the types OH···O, OH···N, NH···O and OH···C, in 4244 conformers of the DNA-related molecules (four canonical 2'-deoxyribonucleotides, 1,2-dideoxyribose-5-phosphate, and 2-deoxy-D-ribose in its furanose, pyranose and linear forms) have been investigated using quantum theory of atoms in molecules (QTAIM) and vibrational analysis. It has been found that for all iHBs with positive red-shift of the proton donating group stretching frequency the shift value correlates with ρ(cp)-the electron charge density at the (3,-1)-type bond critical point. Combining QTAIM and spectroscopic data new relationships for estimation of OH···O, OH···N, NH···O and OH···C iHB enthalpy of formation (kcal mol(-1)) with RMS error below 0.8 kcal mol(-1) have been established: E(OH···O) = -3.09 + 239·ρ(cp), E(OH···N) = 1.72 + 142·ρ(cp), E(NH···O) = -2.03 + 225·ρ(cp), E(OH···C) = -0.29 + 288·ρ(cp), where ρ(cp) is in e a(0)(-3) (a(0)- the Bohr radius). It has been shown that XHY iHBs with red-shift values over 40 cm(-1) are characterized by the following minimal values of the XHY angle, ρ(cp) and nubla(2)ρ(cp): 112°, 0.005 e a(0)(-3) and 0.016 e a(0)(-5), respectively. New relationships have been used to reveal the strongest iHBs in canonical 2'-deoxy- and ribonucleosides and the O(5')H···N(3) H-bond in ribonucleoside guanosine was found to have the maximum energy (8.1 kcal mol(-1)).

No abstract

We present the procedure for transforming delocalized molecular orbitals into the localized property-optimized orbitals (LPOs) designed for building the most accurate, in the Frobenius norm sense, approximation to the first-order reduced density matrix in form of the sum of localized monoatomic and diatomic terms. In this way, a decomposition of molecular properties into contributions associated with individual atoms and the pairs of atoms is obtained with the a priori known upper bound for the decomposition accuracy. Additional algorithm is proposed for obtaining the set of "the Chemist's LPOs" (CLPOs) containing a single localized orbital, with nearly double occupancy, per a pair of electrons. CLPOs form an idealized Lewis structure optimized for the closest possible reproduction of one-electron properties derived from the original many-electron wavefunction. The computational algorithms for constructing LPOs and CLPOs from a general wavefunction are presented and their implementation within the open-source freeware program JANPA (http://janpa.sourceforge.net/) is discussed. The performance of the proposed procedures is assessed using the test set of density matrices of 33 432 small molecules obtained at both Hartree-Fock and second-order Moller-Plesset theory levels and excellent agreement with the chemist's Lewis-structure picture is found. K E Y W O R D S algorithms and software, chemist's Lewis-structure picture, localized orbitals, properties decomposition, quantum-chemical methods | INTRODUCTIONLocalized molecular orbitals resulting from a unitary transformation of occupied canonical molecular orbitals (MOs) [1] play essential role in physical chemistry as the "building blocks" or "descriptors" through which the complicated electronic structure of atoms and molecules can be modeled or interpreted in a comprehensible way. In addition, the localized orbitals concentrated in a limited spatial region of a molecule proved useful in making the high-level correlated quantum-chemical methods more computationally tractable. [2][3][4][5][6][7][8] Although the concept of an orbital itself has been much methodologically debatable and is too "fuzzy" [34] to be defined more precisely than just as the function of coordinates of a single electron somehow related to the system's wavefunction or electron density, this concept still remains virtually the best one proposed so far for moving the ideas of valence electrons and electron pairs (including bonding and antibonding orbitals, lone pairs etc.), which are the key elements of the chemist's Lewis-structure picture, [35] from qualitative to a quantum-mechanical ground. From this perspective, the localized orbitals are used to decompose (typically in an approximate manner) true many-electron wavefunction of the molecule into the components allowing a chemically meaningful interpretation.To transform delocalized orbitals (either the canonical MOs obtained as a solutions of self-consistent field Hartree-Fock or Kohn-Sham equations, [36] or the Lowding natural orbitals [37] ...

Relaxed force constants (RFCs) and vibrational root-mean-square deviations have been evaluated by the original calculation method for conformational parameters of the DNA structural units and their constituents: nucleic acid bases (uracile, thymine, cytosine, adenine and guanine) and their 'building blocks' (benzene, pyrimidine, imidazole and purine molecules), as well as the DNA backbone structural units: tetrahydrofuran, 1,2-dideoxyribose, methanol and orthophosphoric acid. It has been found that the RFCs for nomenclature torsions beta, gamma, epsilon; and sugar pseudorotation angle P in 1,2-dideoxyribose are sensible to the molecule conformation and their values are in the range of 1-25 kcal/(mole·rad²) obeying the inequality K(γ)> K(ε) > K(ρ) > K(β). The RFCs values for endocyclic torsions of nucleic acid bases six-member rings lie within 15-45 kcal/(mole·rad²) in pyrimidines and within 20-60 kcal/(mole·rad²) in purines. It is shown that the quantum zero-point motion effectively neglects the amino group non-planarity in cytosine, adenine and partially in guanine.

No abstract