Quantum chemical MP2 results on some hydrates of cytosine: binding sites, energies and the first hydration shell

Computational analysis of the vibrational spectra and structure of aqueous cytosine

The recently developed efficient protocol for the explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is used to describe aqueous solutions of cytosine. The same cluster model of a solute surrounded by the first solvation shell of solvent molecules was shown to be sufficient to reproduce experimental vibrational frequencies and relative IR and Raman intensities. An equally good quality of Raman spectra was provided by B3LYP-D3/def2-TZVP and B3LYP-D3/aug-cc-pVDZ simulations. Computations using the PBE functional were sufficient for modeling of the IR spectra but failed in the simulations of Raman scattering. It is shown that strong changes of frequencies and relative intensities of Raman and IR bands of cytosine, caused by its hydration, cannot be completely assigned to the influence of hydrogen bonds (HBs) with 7 or 8 closest water molecules. They are rather ascribed to the combined effect of solute-solute and solute-solvent HBs with the participation of at least 30 water molecules separating cytosine from the bulk solvent. This suggests that the vibrational modes and derivatives of the polarizability and dipole moment of the solute are mainly locally influenced by its first hydration shell, while the influence of bulk water is rather modest.

Infrared spectra of PEHA molecule and its resistance to oxidation in water and methanol media at 298.15 K: solvent cluster size dependency

The present work investigates the infrared spectra and solvation free energies (SE) of PEHA ((E)-2-(Pyridin-2-yl) ethyl 3-(3,4-dihydroxyphenyl) acrylate) and their impact on the oxidation. The latter was examined through the ionization potential parameter (IP). These investigations were carried out by the DFT method at B3LYP/6-31G(d) for optimization and frequency calculations and corrected for BSSE. X3LYP/6-311++G(2d,2p) was employed for single-point energy calculations. Water and methanol cluster sizes were used for solvation through the explicit solvent model. Thus, the infrared spectra show that the overview frequencies of PEHA compare well with the experimental results. The intense infrared absorptions of complexes are due to the stretching of O-H bonds of solvent clusters in the range 2600-3850 cm^-1. The binding energy per solvent molecule of complexes was calculated and shows that water and methanol clusters mimic the liquid state as from 5 to 10 solvent molecules. The SE of PEHA increases with the increase of the cluster size of water and methanol in the direction of the limit. The latter was censured by the solvation done using the combined explicit-implicit solvent model. As for IP parameter, the results are largely above the IP limit and lower than the IP from gas phase. Thus, water and methanol media have an effect of lowering the IP of PEHA compound. Consequently, both media favour the oxidation of PEHA molecule, which facilitates its metabolism in human organism.

Ultrafast Excited-State Dynamics of Hydrogen-Bonded Cytosine Microsolvated Clusters with Protic and Aprotic Polar Solvents

Microsolvation effects on the ultrafast excited-state deactivation dynamics of cytosine (Cy) were studied in hydrogen-bonded Cy clusters with protic and aprotic solvents using mass-resolved femtosecond pump-probe ionization spectroscopy. Two protic solvents, water (H₂O) and methanol (MeOH), and one aprotic solvent, tetrahydrofuran (THF), were investigated, and transients of Cy·(H₂O)_1-6, Cy·(MeOH)_1-3, and Cy·THF microsolvated clusters produced in supersonic expansions were measured. With the aid of electronic structure calculations, we assigned the observed dynamics to the low-energy isomers of various Cy clusters and discussed the microsolvation effect on the excited-state deactivation dynamics. With the protic solvents only the microsolvated clusters of Cy keto tautomer were observed. The observed decay time constants of Cy·(H₂O) _n are 0.5 ps for n = 1 and ∼0.2-0.25 ps for n = 2-6. For Cy·(MeOH) _n clusters, the decay time constant for n = 1 cluster is similar to that of the Cy monohydrate, but for n = 2 and 3 the decays are about a factor of 2 slower than the corresponding microhydrates. With the aprotic solvent, THF, hydrogen-bonded complexes of both keto and enol tautomers are present in the beam. The keto-Cy·THF shows a decay similar to that of the keto-Cy monomer, whereas the enol-Cy·THF exhibits a 2-fold slower decay than the enol-Cy monomer, suggesting an increase in the barrier to excited-state deactivation upon binding of one THF molecule to the enol form of Cy.

Efficient structural elucidation of microhydrated biomolecules through the interrogation of hydrogen bond networks

Pattern analysis of H-bond networks through a graph-theoretic method is very effective in determining the global minima of microhydrated biomolecules.