Uracil-water interaction revisited - in search of single H-bonded secondary minima

Monohydrated uracil (UW) complexes are stabilized by both O⋯HO and NH⋯O hydrogen bonds (H-bonds), simultaneously participating in forming three stable cyclic structures. The role and contribution of these individual H-bonds (O⋯HO and NH⋯O) to the stability of the three UW complexes are still not understood, because of the technical problems in obtaining their optimized structures by standard geometry optimization. The present study explores a non-standard approach to identify three single H-bonded local minima structures without imaginary frequency using DFT (M06-2X, B3LYP and B3LYP-D3), MP2 and CCSD(T) theories and Dunning's correlation-consistent aug-cc-pVTZ basis set, in both vacuum and aqueous media (CPCM method). The results reveal that these new structures are very shallow secondary minima between two deep wells or next to a deep well of primary minima (double H-bonded structures) in the potential energy surface. The H-bond energy of these single H-bonded complexes is found to be less sensitive to a wide range (about 15-20 degrees) of O⋯HO and NH⋯O angles, and the linearity is preferred in the stable three single H-bonded structures. The technical method used to locate such a shallow minimum is described in detail and may be useful for identifying local minima in other cases where consecutive multiple H-bonded structures are global minima. Energy decomposition (using symmetry adapted perturbation theory (SAPT)) of interaction energy, electron redistribution, and relevant vibrational modes are discussed.

Influence of a Methyl Group on the Unidirectional Flow of Vibrational Energy in an Adenine-Thymine Base Pair

The role of a methyl group in intramolecular vibrational energy redistribution (IVR) of the hydrogen-bonded adenine-thymine base pair has been studied using classical dynamics procedures. Energy transferred to the doorway bond thymine-NH from the vibrationally excited H₂O(v) efficiently redistributes among various bonds of the base pair through vibration-to-vibration coupling, depositing a large fraction of the available energy in the terminal bond adenine-NH. On the other hand, the extent of energy flow in the reverse direction from the excited adenine-NH to thymine-NH is insignificant, indicating IVR in adenine-thymine resulting from the intermolecular interaction with a vibrationally excited H₂O molecule, is direction-specific. The unidirectional flow is due to the coupling of stretch-torsion vibrations of a methyl group with conjugated bonds on the thymine ring, when the methyl rotor is present and is adjacent to the vibrationally excited thymine-NH. The insignificance of energy flow from the terminal-to-terminal bond in the reverse direction is attributed to the absence of a methyl group on the adenine moiety, even though the molecule has many CC and CN bonds coupled to their neighbors.

Rotational spectrum and internal dynamics of the hydrogen-bonded pyrrole-pyridine aromatic pair

Non-covalent interactions determine the three-dimensional structure and activity of biological molecules. In this work, the pyrrole-pyridine complex considered as a model of the NH⋯N hydrogen-bonded Watson-Crick base pairs has been generated in a supersonic expansion and characterized by chirped pulse Fourier transform microwave spectroscopy. The analysis of the unconventional spectral pattern of the 1:1 pyrrole-pyridine adduct and its ¹³C and ¹⁵N isotopologues reveal a non-planar complex, with a bent NH⋯N hydrogen bond and large amplitude motion of the pyrrole subunit. The bent structure is likely to arise from the stablishment of the secondary CH⋯N interaction between pyridine and pyrrole moieties.

Vibrational Energy Flow in the Uracil-H2O Complexes

In the uracil-H₂O complex, the vibrational energy initially stored in the OH(v = 1) stretch efficiently transfers to the first overtone-bending mode under a near-resonant condition. The relaxation of the overtone vibration redistributes its energy to uracil and the two hydrogen bonds in the intermolecular zone, which consists of the OH bond and the bonds between nearby C, N, O, and H atoms of uracil. The uracil NH bond and the hydrogen bond it formed with the H₂O molecule, N-H···O, store the major portion of the energy released by the relaxing bending mode, thus forming a localized hot band in the intermolecular zone. Energy transfer to the bonds beyond the zone is found to be not significant. The excited uracil NH is found to transfer its energy to the bending mode, thus indicating that the hydrogen bond of N-H···O is the principal energy pathway in both directions. In the presence of efficient near-resonant energy transfer pathways, the time evolution of the centers of mass distance shows the phenomenon of beats. One global and two different local minima energy structures are considered. The results of energy transfer do not differ significantly, suggesting that the two hydrogen bonds in all three structures have similar contributions to the energy transfer.

Exploration of the theobromine-water dimer: comparison with DNA microhydration

Understanding the molecular basis of the appearance of life on Earth is an exciting research field. Many factors may have influenced the election of the molecules used by living beings and evolution may have modified those original compounds. In an attempt to understand the role played by intermolecular interactions in the election of CGAT as the alphabet of life, we present here a thorough experimental and computational study on the interaction of theobromine with water. Theobromine is a xanthine derivative, structurally related to the nucleobases, and also present in many living beings. The experimental results demonstrate that the most stable isomer of theobromine-water was formed and detected in supersonic expansions. This isomer very well resembles the structure of the dimers between nucleobases and water, offering similar values of binding energy. A comparison between the results obtained for theobromine-water with those reported in the literature for monohydrates of nucleobases is also offered.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

Where's water? The many binding sites of hydantoin

Investigation of an astrobiologically relevant prebiotic molecule, hydantoin, by broadband rotational spectroscopy and its plethora of water complexes.

Solvation of Uracil and Its Derivatives by DMSO: A DFT-Supported 1H NMR and 13C NMR Study

¹H NMR and ¹³C NMR spectra of uracil, thymine, 5-hydroxymethyluracil, 5,6-dihydrouracil, and 5,6-dihydrothymine in DMSO-d₆ solutions have been measured. Additionally, molecular structures as well as NMR parameters of these compounds and their various solvates have been calculated using DFT B3LYP/6-311++G(2d,p) PCM(DMSO) method. The analysis of the chemical shift data for these compounds has shown that, indeed, in DMSO solutions they occur as equilibrium mixtures of free molecules and solvates in which solute and solvent molecules are joined by NH···O or OH···O hydrogen bonds. The populations of particular species present in the solutions have been estimated. Moreover, it has been found that 5,6-dihydrothymine exists in DMSO solution preferentially in conformation with the methyl group occupying the pseudoequatorial position. This finding is based on the molecular energy calculations and remains in full agreement with the interpretation of NMR data and theoretical calculations of NMR parameters.

Clustering of Uracil Molecules on Ice Nanoparticles

We generate a molecular beam of ice nanoparticles (H₂O)_N, N̅ ≈ 130-220, which picks up several individual gas phase uracil (U) or 5-bromouracil (BrU) molecules. The mass spectra of the doped nanoparticles prove that the uracil and bromouracil molecules coagulate to clusters on the ice nanoparticles. Calculations of U and BrU monomers and dimers on the ice nanoparticles provide theoretical support for the cluster formation. The (U)_mH⁺ and (BrU)_mH⁺ intensity dependencies on m extracted from the mass spectra suggest a smaller tendency of BrU to coagulate compared to U, which is substantiated by a lower mobility of bromouracil on the ice surface. The hydrated U_m·(H₂O)_nH⁺ series are also reported and discussed. On the basis of comparison with the previous experiments, we suggest that the observed propensity for aggregation on ice nanoparticles is a more general trend for biomolecules forming strong hydrogen bonds. This, together with their mobility, leads to their coagulation on ice nanoparticles which is an important aspect for astrochemistry.

Gas-Phase IR Spectroscopy of Nucleobases

IR spectroscopy of nucleobases in the gas phase reflects simultaneous advances in both experimental and computational techniques. Important properties, such as excited state dynamics, depend in subtle ways on structure variations, which can be followed by their infrared signatures. Isomer specific spectroscopy is a particularly powerful tool for studying the effects of nucleobase tautomeric form and base pair hydrogen-bonding patterns.

Characterizing the dark state in thymine and uracil by double resonant spectroscopy and quantum computation

Double resonant spectroscopy characterizes both grounds state and dark excited state of uracil and thymine.

Microsolvation of 2-azetidinone: a model for the peptide group–water interactions

The geometries of the 2-azetidinone–(H₂O)_n clusters, determined by rotational spectroscopy, show the preference of water to interact with the CO group and the effects of cooperative hydrogen bonding.

Microsolvation of 2-thiouracil: molecular structure and spectroscopic parameters of the thiouracil-water complex

State-of-the-art quantum-chemical computations have been employed to accurately determine the equilibrium structure and interaction energy of the 2-thiouracil-water complex, thus extending available reference data for biomolecule solvation patterns. The coupled-cluster level of theory in conjunction with a triple-ζ basis set has been considered together with extrapolation to the basis set limit, performed by employing second-order Møller-Plesset perturbation theory, and inclusion of core-correlation and diffuse-function corrections. On the basis of the comparison of experiment and theory for 2-thiouracil [ Puzzarini et al. Phys. Chem. Chem. Phys. 2013 , 15 , 16965 - 16975 ], structural changes due to water complexation have been pointed out. Molecular and spectroscopic properties of the 2-thiouracil-water complex have then been studied by means of the composite computational approach introduced for the molecular structure evaluation. Among the results achieved, we mention the accurate determination of the molecular dipole moment and of the spectroscopic parameters required for predicting the rotational spectrum.

Rotational spectroscopy meets theory

Rotational spectroscopy is known to be a technique that is widely used to infer information on molecular structure and dynamics. In the last few decades, its role in the field of atmospheric and astrophysical investigations has rapidly grown. However, several are the challenging aspects in rotational spectroscopy, since the detection and analysis of spectra as well as interpretation of obtained results are not at all straightforward. Quantum chemistry has reached such an accuracy that can be used to disentangle these challenging situations by guiding the experimental investigation, assisting in the determination of the spectroscopic parameters, and extracting information of chemical interest. This perspective provides an overview of the theoretical background and computational requirements needed for the accurate evaluation of the spectroscopic parameters of relevance to rotational spectroscopy. The role of theory in guiding and supporting experiment is detailed through a few examples, and the interplay of experiment and theory is discussed in terms of the information of physical and chemical interest that can be derived.

Defining the hydrogen bond: An account (IUPAC Technical Report)

The term “hydrogen bond” has been used in the literature for nearly a century now. While its importance has been realized by physicists, chemists, biologists, and material scientists, there has been a continual debate about what this term means. This debate has intensified following some important experimental results, especially in the last decade, which questioned the basis of the traditional view on hydrogen bonding. Most important among them are the direct experimental evidence for a partial covalent nature and the observation of a blue-shift in stretching frequency following X–H···Y hydrogen bond formation (XH being the hydrogen bond donor and Y being the hydrogen bond acceptor). Considering the recent experimental and theoretical advances, we have proposed a new definition of the hydrogen bond, which emphasizes the need for evidence. A list of criteria has been provided, and these can be used as evidence for the hydrogen bond formation. This list is followed by some characteristics that are observed in typical hydrogen-bonding environments.