Publication details

Stacked and H-Bonded Cytosine Dimers. Analysis of the Intermolecular Interaction Energies by Parallel Quantum Chemistry and Polarizable Molecular Mechanics.

Investor logo
Authors

GRESH Nohad ŠPONEROVÁ Judit DEVEREUX Mike GKIONIS Konstantinos COURCY DE Benoit PIQUEMAL Jean-Philip ŠPONER Jiří

Year of publication 2015
Type Article in Periodical
Magazine / Source Journal of Physical Chemistry B
MU Faculty or unit

Central European Institute of Technology

Citation
Web http://pubs.acs.org/doi/pdf/10.1021/acs.jpcb.5b01695
Doi http://dx.doi.org/10.1021/acs.jpcb.5b01695
Field Physical chemistry and theoretical chemistry
Keywords DENSITY-FUNCTIONAL THEORY; DISTRIBUTED MULTIPOLE ANALYSIS; PERTURBATION-THEORY APPROACH; BASIS-SET CONVERGENCE; KOHN-SHAM ORBITALS; NUCLEIC-ACID BASES; AMBER FORCE-FIELD; DYNAMICS SIMULATIONS; CHEMICAL COMPUTATIONS; CHARGE-DISTRIBUTION
Description Until now, atomistic simulations of DNA and RNA and their complexes have been executed using well calibrated but conceptually simple pair-additive empirical potentials (force fields). Although such simulations Provided Many valuable results, it is well established that simple force fields also introduce errors into the description, underlying the need for development of alternative anisotropic, polarizable molecular mechanics (APMM) potentials. One of the most abundant forces in all kinds of nucleic acids topologies is base stacking. Intra- and interstrand stacking is assumed to be the most essential factor affecting local conformational variations of B-DNA. However, stacking also contributes to formation of all kinds of noncanonical nucleic acids structures, such as quadruplexes or folded RNAs. The present study focuses on 14 stacked cytosine (Cyt) dimers and the doubly H-bonded dimer. We evaluate the extent to which an APMM procedure, SIBFA, could account quantitatively for the results of high-level quantum chemistry (QC) on the total interaction energies, and the individual energy contributions and their nonisotropic behaviors. Good agreements are found at both uncorrelated HF and correlated DFT and CCSD(T) levels. Resorting in SIBFA to distributed QC multipoles and to an explicit representation of the lone pairs is essential to respectively account for the anisotropies of the Coulomb and of the exchange-repulsion QC contributions.
Related projects:

You are running an old browser version. We recommend updating your browser to its latest version.

More info