Homodimers of Cytosine and 1-MethylCytosine. A DFT study of geometry, relative stability and H-NMR shifts in gas-phase and selected solvents

Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking

Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.

Using Density Functional Theory To Study Neutral and Ionized Stacked Thymine Dimers

Stacking interactions in thymine dimers are studied with density functional theory. According to our calculations, six dimers of comparable stability can be prepared at low temperature, but dimerization is impossible at room temperature due to the large entropy contribution that accompanies it. Analysis of vibrational anharmonic coupling terms shows that each of the dimers exhibits distinct vibrational dynamics. Properties of electron density in the intermolecular region are used to analyze neutral stacked species and their ionized forms. Bond paths and critical points in the intermolecular region are identified, but a simple relationship between binding energy and total electron density in the intermolecular critical points could not be found due to an uneven electron distribution in the binding region. The reduced density gradient was confirmed to be a useful tool for analysis of weak stacking interactions. Those interactions also affect vertical and adiabatic ionization energies, which are computed to be slightly lower for the dimers compared to the monomer.

Theoretical and experimental NMR studies on muscimol from fly agaric mushroom (Amanita muscaria)

In this article we report results of combined theoretical and experimental NMR studies on muscimol, the bioactive alkaloid from fly agaric mushroom (Amanita muscaria). The assignment of (1)H and (13)C NMR spectra of muscimol in DMSO-d6 was supported by additional two-dimensional heteronuclear correlated spectra (2D NMR) and gauge independent atomic orbital (GIAO) NMR calculations using density functional theory (DFT). The effect of solvent in theoretical calculations was included via polarized continuum model (PCM) and the hybrid three-parameter B3LYP density functional in combination with 6-311++G(3df,2pd) basis set enabled calculation of reliable structures of non-ionized (neutral) molecule and its NH and zwitterionic forms in the gas phase, chloroform, DMSO and water. GIAO NMR calculations, using equilibrium and rovibrationally averaged geometry, at B3LYP/6-31G* and B3LYP/aug-cc-pVTZ-J levels of theory provided muscimol nuclear magnetic shieldings. The theoretical proton and carbon chemical shifts were critically compared with experimental NMR spectra measured in DMSO. Our results provide useful information on its structure in solution. We believe that such data could improve the understanding of basic features of muscimol at atomistic level and provide another tool in studies related to GABA analogs.