Unraveling Water Solvation Effects with Quantum Mechanics/Molecular Mechanics Semiclassical Vibrational Spectroscopy: The Case of Thymidine

We introduce a quantum mechanics/molecular mechanics semiclassical method for studying the solvation process of molecules in water at the nuclear quantum mechanical level with atomistic detail. We employ it in vibrational spectroscopy calculations because this is a tool that is very sensitive to the molecular environment. Specifically, we look at the vibrational spectroscopy of thymidine in liquid water. We find that the C═O frequency red shift and the C═C frequency blue shift, experienced by thymidyne upon solvation, are mainly due to reciprocal polarization effects, that the molecule and the water solvent exert on each other, and nuclear zero-point energy effects. In general, this work provides an accurate and practical tool to study quantum vibrational spectroscopy in solution and condensed phase, incorporating high-level and computationally affordable descriptions of both electronic and nuclear problems.

Mapping electronic decoherence pathways in molecules

Establishing the fundamental chemical principles that govern molecular electronic quantum decoherence has remained an outstanding challenge. Fundamental questions such as how solvent and intramolecular vibrations or chemical functionalization contribute to the decoherence remain unanswered and are beyond the reach of state-of-the-art theoretical and experimental approaches. Here we address this challenge by developing a strategy to isolate electronic decoherence pathways for molecular chromophores immersed in condensed phase environments that enables elucidating how electronic quantum coherence is lost. For this, we first identify resonance Raman spectroscopy as a general experimental method to reconstruct molecular spectral densities with full chemical complexity at room temperature, in solvent, and for fluorescent and non-fluorescent molecules. We then show how to quantitatively capture the decoherence dynamics from the spectral density and identify decoherence pathways by decomposing the overall coherence loss into contributions due to individual molecular vibrations and solvent modes. We illustrate the utility of the strategy by analyzing the electronic decoherence pathways of the DNA base thymine in water. Its electronic coherences decay in [Formula: see text]30 fs. The early-time decoherence is determined by intramolecular vibrations while the overall decay by solvent. Chemical substitution of thymine modulates the decoherence with hydrogen-bond interactions of the thymine ring with water leading to the fastest decoherence. Increasing temperature leads to faster decoherence as it enhances the importance of solvent contributions but leaves the early-time decoherence dynamics intact. The developed strategy opens key opportunities to establish the connection between molecular structure and quantum decoherence as needed to develop chemical strategies to rationally modulate it.

Quantum Vibrational Spectroscopy of Explicitly Solvated Thymidine in Semiclassical Approximation

In this paper, we demonstrate the possibility to perform spectroscopy simulations of solvated biological species taking into consideration quantum effects and explicit solvation. We achieve this goal by interfacing our recently developed divide-and-conquer approach for semiclassical initial value representation molecular dynamics with the polarizable AMOEBABIO18 force field. The method is applied to the study of solvation of the thymidine nucleoside in two different polar solvents, water and N,N-dimethylformamide. Such systems are made of up to 2476 atoms. Experimental evidence concerning the different behavior of thymidine in the two solvents is well reproduced by our study, even though quantitative estimates are hampered by the limited accuracy of the classical force field employed. Overall, this study shows that semiclassically approximate quantum dynamical studies of explicitly solvated biological systems are both computationally affordable and insightful.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Initial excited-state structural dynamics of thymine derivatives

Thymine is one of the pyrimidine nucleobases found in DNA. Upon absorption of UV light, thymine forms a number of photoproducts, including the cyclobutyl photodimer, the pyrimidine pyrimidinone [6-4] photoproduct and the photohydrate. Here, we use UV resonance Raman spectroscopy to measure the initial excited-state structural dynamics of the N(1)-substituted thymine derivatives N(1)-methylthymine, thymidine, and thymidine 5'-monophosphate in an effort to understand the role of the N1 substituent in determining the excited-state structural dynamics and the subsequent photochemistry. The UV resonance Raman spectrum of thymidine and thymidine 5'-monophosphate are similar to that of thymine, suggesting that large masses at N(1) effectively isolate the vibrations of the nucleobase. However, the UV resonance Raman spectrum of N(1)-methylthymine is significantly different, suggesting that the methyl group couples into the thymine ring vibrations. The resulting resonance Raman intensities and absorption spectra are self-consistently simulated with a time-dependent expression to quantitatively extract the initial excited-state slopes, homogeneous and inhomogeneous linewidths, and electronic parameters. These results are discussed in the context of the known photochemistry of thymine and its derivatives.

Initial Excited-State Structural Dynamics of Thymine Are Coincident with the Expected Photochemical Dynamics

To explore the excited-state structural dynamics of thymine, a DNA nucleobase, we measured the resonance Raman spectra of thymine in aqueous solution at wavelengths throughout the lowest-energy absorption band. Self-consistent analysis of the resulting resonance Raman excitation profiles and absorption spectrum using a time-dependent wave packet formalism yielded the excited-state structural dynamics. The photochemically relevant C=C stretching and C-H deformation vibrational modes were found to exhibit maximum resonance Raman intensity and structural change upon photoexcitation for thymine, suggesting that the initial dynamics of thymine lie along the photochemical reaction coordinate.

Combined use of Overhauser spectroscopy and NMR diffusion experiments for mapping the hydration structure of nucleosides: structure and dynamics of uridine in water

Complementary results from 13C intermolecular nuclear Overhauser effects (NOE), 1H-13C heteronuclear Overhauser spectroscopy (HOSEY) and 1H-NMR diffusion measurements were used for probing the structure of the first solvation shell of uridine in water. It is demonstrated that a cyclic dihydrate is formed. The two water molecules produce two hydrogen bonds with the two oxygen atoms from the pyrimidine ring and accept only one hydrogen bond from the amide proton. The dihydrate has only a short lifetime as compared with the rotational correlation time of the free nucleoside. The chemical exchange constant of the amide proton with water is then estimated by diffusion experiments. The results are consistent with previous data obtained for uracil in water and provide interesting information about water accessibility in nucleic acid bases.