Computationally assisted vibrational spectroscopy of nucleic acid bases. 2. Thymine

As in the case of cytosine [Phys. Chem. Chem. Phys. 2023, 25, 24121-24128], Raman and infrared (IR) spectra of aqueous thymine and its N-deuterated derivative, thymine-d2 have been computationally reproduced and interpreted with the use of the recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions [J. Phys. Chem. B, 2020, 124, 6664-6670]. A cluster model of a solute surrounded by 30 water molecules is shown to be sufficient to reproduce experimental vibrational frequencies and relative Raman intensities with the use of B3LYP-D3/def2-TZVP or B3LYP-D3/aug-cc-pVDZ simulations. Analogous PBE-D3 computations provided a less good, but still reasonably accurate, modeling of Raman spectra. It is shown that strong changes of frequencies and relative intensities of the Raman bands of thymine, caused by its hydration, can be interpreted mainly as a result of hydrogen bonding with 6 nearest water molecules. Non-negligible improvement of the quality of simulations for larger clusters comprising water molecules that do not have direct contacts with the solute, suggests that spectroscopic effects of hydration should be ascribed to the joined action of solute-solvent and solvent-solvent interactions. Nevertheless, the moderate number of water molecules required for successful simulations of the Raman spectra of aqueous thymine, suggests that the vibrational modes and derivatives of the polarizability of the solute are mainly locally influenced, while the effect of bulk water is rather modest.

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.

Nucleic Acids under Stress: Understanding and Simulating Nucleobase Fragmentation Pathways

The effects of radiations on nucleic acids and their constituents is widely studied across several research fields using different experimental and theoretical protocols. While a large number of studies were performed in this context, many fundamental physical and chemical effects are still being investigated, particularly involving the effect of the biological environment. As an example, the interpretation of experimental nucleic acid bases mass spectra, and hence inferring their reactivity in complex environment still poses great challenge. This Minireview summarizes recent theoretical advancements aiming to predict and interpret the reactivity of nucleic acid bases. We focus not only on the understanding of the inherent fragmentation pathways of isolated nucleobases but also on the modeling of a realistic nano-environments highlighting the importance of molecular dynamics simulations and the non-innocent role of the environment and also the possibility to open novel fragmentation pathways.

Water-biomolecule clusters studied by photoemission spectroscopy and multilevel atomistic simulations: hydration or solvation?

The properties of mixed water-uracil nanoaggregates have been probed by core electron-photoemission measurements to investigate supramolecular assembly in the gas phase driven by weak interactions. The interpretation of the measurements has been assisted by multilevel atomistic simulations, based on semi-empirical tight-binding and DFT-based methods. Our protocol established a positive-feedback loop between experimental and computational techniques, which has enabled a sound and detailed atomistic description of such complex heterogeneous molecular aggregates. Among biomolecules, uracil offers interesting and generalized skeletal features; its structure encompasses an alternation of hydrophilic H-bond donor and acceptor sites and hydrophobic moieties, typical in biomolecular systems, that induces a supramolecular core-shell-like organization of the mixed clusters with a water core and an uracil shell. This structure is far from typical models of both solid-state hydration, with water molecules in defined positions, or liquid solvation, where disconnected uracil molecules are completely surrounded by water.

Cardiolipin Structure and Oxidation Are Affected by Ca2+ at the Interface of Lipid Bilayers

Ca²⁺-overload contributes to the oxidation of mitochondrial membrane lipids and associated events such as the permeability transition pore (MPTP) opening. Numerous experimental studies about the Ca²⁺/cardiolipin (CL) interaction are reported in the literature, but there are few studies in conjunction with theoretical approaches based on ab initio calculations. In the present study, the lipid fraction of the inner mitochondrial membrane was modeled as POPC/CL large unilamellar vesicles (LUVs). POPC/CL and, comparatively, POPC, and CL LUVs were challenged by singlet molecular oxygen using the anionic porphyrin TPPS4 as a photosensitizer and by free radicals produced by Fe²⁺-citrate. Calcium ion favored both types of lipid oxidation in a lipid composition-dependent manner. In membranes containing predominantly or exclusively POPC, Ca²⁺ increased the oxidation at later reaction times while the oxidation of CL membranes was exacerbated at the early times of reaction. Considering that Ca²⁺ interaction affects the lipid structure and packing, density functional theory (DFT) calculations were applied to the Ca²⁺ association with totally and partially protonated and deprotonated CL, in the presence of water. The interaction of totally and partially protonated CL head groups with Ca²⁺ decreased the intramolecular P-P distance and increased the hydrophobic volume of the acyl chains. Consistently with the theoretically predicted effect of Ca²⁺ on CL, in the absence of pro-oxidants, giant unilamellar vesicles (GUVs) challenged by Ca²⁺ formed buds and many internal vesicles. Therefore, Ca²⁺ induces changes in CL packing and increases the susceptibility of CL to the oxidation promoted by free radicals and excited species.

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Intermolecular interactions in crystals of benzene and its mono- and dinitro derivatives: study from the energetic viewpoint

The weak intermolecular interactions and their role in mono- and dinitrobenzene crystal structure formation have been studied using quantum-chemical calculations.

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

A detailed quantum chemical investigation was undertaken to obtain the structure and energetics of cytosine hydrates Cyt·nH2O, with n = 1 to 7. The MP2(fc)/aug-cc-pVDZ level was used as the standard, with some DFT (B3LYP) and coupled cluster calculations, as well as calculations with the aug-cc-pVTZ basis set added for comparison. In a systematic search for microhydrated forms of cytosine, we have found that several structures have not yet been reported in the literature. The energies of different isomers, as well as binding energies are compared. When predicting the stability of a complex, we suggest using a scheme where the water molecules are extracted from a finite model of bulk water. Finally, based on energetic data, we suggest a rational definition of the first hydration shell; with this definition, it contains just six water molecules.

DNA spontaneous mutation and its role in the evolution of GC-content: assessing the impact of the genetic sequence

We use theoretical tools to investigate the possible role played by a DNA sequence in the base pair tautomerization phenomena.

Hydration properties of natural and synthetic DNA sequences with methylated adenine or cytosine bases in the R.DpnI target and BDNF promoter studied by molecular dynamics simulations

Adenine and cytosine methylation are two important epigenetic modifications of DNA sequences at the levels of the genome and transcriptome. To characterize the differential roles of methylating adenine or cytosine with respect to their hydration properties, we performed conventional MD simulations and free energy perturbation calculations for two particular DNA sequences, namely the brain-derived neurotrophic factor (BDNF) promoter and the R.DpnI-bound DNA that are known to undergo methylation of C5-methyl cytosine and N6-methyl adenine, respectively. We found that a single methylated cytosine has a clearly favorable hydration free energy over cytosine since the attached methyl group has a slightly polar character. In contrast, capping the strongly polar N6 of adenine with a methyl group gives a slightly unfavorable contribution to its free energy of solvation. Performing the same demethylation in the context of a DNA double-strand gave quite similar results for the more solvent-accessible cytosine but much more unfavorable results for the rather buried adenine. Interestingly, the same demethylation reactions are far more unfavorable when performed in the context of the opposite (BDNF or R.DpnI target) sequence. This suggests a natural preference for methylation in a specific sequence context. In addition, free energy calculations for demethylating adenine or cytosine in the context of B-DNA vs. Z-DNA suggest that the conformational B-Z transition of DNA transition is rather a property of cytosine methylated sequences but is not preferable for the adenine-methylated sequences investigated here.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

Guanine bases in DNA G-quadruplex adopt nonplanar geometries owing to solvation and base pairing

The effect of base pairing and solvation on pyramidalization of the glycosidic nitrogen found in the residues of parallel G-quadruplex with NDB ID UDF062 is analyzed and explained with theoretical calculations. The extent of the pyramidalization depends on the local geometry of the 2'-deoxyguanosine residues, namely on their glycosidic torsion and sugar pucker, which are predetermined by the 3D-architecture of G-quadruplex. Pyramidal inversion of the glycosidic nitrogen found in 2'-deoxyguanosines of G-quadruplex is induced owing to site-specifically coordinated solvent. Different adiabatic structural constraints used for fixing the base-to-sugar orientation of 2'-deoxyguanosine in geometry optimizations result in different extents of pyramidalization and induce pyramidal inversion of the glycosidic nitrogen. These model geometry constraints helped us analyze the effect of real constraints represented by explicit molecular environment of selected residues of the G-quadruplex. The maximal extent of the glycosidic nitrogen pyramidalization found in the high-resolution crystal structure corresponds to the calculation to deformation energy of only 1 kcal mol(-1). The out-of-plane deformations of nucleobases thus provide a way for compensating the site-specific external environmental stress on the G-quadruplex.