Four Tautomers of Isolated Guanine from Infrared Laser Spectroscopy in Helium Nanodroplets

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

Tautomeric equilibrium and spectroscopic properties of 8-azaguanine revealed by quantum chemistry methods

8-azaguanine is a triazolopyrimidine nucleobase analog possessing potent antibacterial and antitumor activities, and it has been implicated as a lead molecule in cancer and malaria therapy. Its intrinsic fluorescence properties can be utilized for monitoring its interactions with biological polymers like proteins or nucleic acids. In order to better understand these interactions, it is important to know the tautomeric equilibrium of this compound. In this work, the tautomeric equilibrium of all natural neutral and anionic compound forms (except highly improbable imino-enol tautomers) as well as their methyl derivatives and ribosides was revealed by quantum chemistry methods. It was shown that, as expected, tautomers protonated at positions 1 and 9 dominate neutral forms both in gas phase and in aqueous solution. 8-azaguanines methylated at any position of the triazole ring are protonated at position 1. The computed vertical absorption and emission energies are in very good agreement with the experimental data. They confirm the validity of the assumption that replacing the proton with the methyl group does not significantly change the positions of absorption and fluorescence peaks.

Accurate Structures and Spectroscopic Parameters of Guanine Tautomers in the Gas Phase by the Pisa Conventional and Explicitly Correlated Composite Schemes (PCS and PCS-F12)

A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase is proposed and validated for tautomeric equilibria. The main features of the new model are the inclusion of core-valence correlation in geometry optimizations by a double hybrid functional and the systematic use of wave-function composite methods in conjunction with cc-pVnZ-F12 basis sets with separate extrapolation of MP2 and post-MP2 contributions. The resulting Pisa composite scheme employing conventional (PCS) or explicitly correlated (PCS-F12) approaches is applied to the challenging problem of guanine tautomers in the gas phase. The results are in remarkable agreement with the experimental structures, relative stabilities, and spectroscopic signatures of different tautomers. The accuracy of the results obtained at reasonable cost by means of black-box parameter-free approaches paves the way toward systematic investigations of other molecular bricks of life also by non-specialists.

Vibrational spectroscopy simulation of solvation effects on a G-quadruplex

It is commonly believed that solvation effects on the vibrational properties of a solute are easily accounted for by simple rules of thumbs, that is, solvating a polar molecule in a polar medium has the only effect of red shifting all its spectroscopical features and, similarly, solvating a polar molecule in a nonpolar medium has the opposite effect. In this work, we use theoretical vibrational spectroscopy at quasi-classical and quantum approximate semiclassical level to gain atomistic insights about solvent-solute interactions for 2'-deoxyguanosine and the G-quadruplex. We employ the quasi-classical trajectory method to include full anharmonicity into our calculated spectra, and then introduce quantum nuclear effects by means of divide-and-conquer semiclassical spectroscopy calculations. Solvation is treated explicitly leading to a good reproducibility of the available experimental data and reliable predictions when an experimental reference is missing.Communicated by Ramaswamy H. Sarma.

Cryogenic IR and UV spectroscopy of isomer-selected cytosine radical cation

The UV photodissociation of cryogenic-cooled isomer-selected cytosine–silver complex leads to the production of cytosine radical cation without isomerization.

Structural Characterization of 6-Thioguanosine and Its Monohydrate in the Gas Phase

Detailed structural analysis of 6-thioguanosine (6TGs) in relation to its tautomerization and sugar conformation is performed in the gas phase using UV and IR spectroscopy combined with ab initio calculations. We have observed a thiol tautomer of 6TGs with its sugar moiety in the syn conformation that is stabilized by a strong intramolecular H-bonding between O5'H of the sugar and N3 atom of the guanine moiety. This observation is consistent with previous results for guanosine (Gs) in which the corresponding enol form is solely detected. We have also identified a monohydrate of 6TGs consisting of a thiol tautomer with the water linking guanine moiety and sugar OH group. It is demonstrated that hydration behavior of 6TGs is significantly different from that of Gs as a result of a weaker H-bonding ability of the thiol group.

Guanine Tautomerism in Ionic Complexes with Ag+ Investigated by IRMPD Spectroscopy and Mass Spectrometry

In this paper, we present the IRMPD spectra of three ionic complexes between guanine (G) and silver (Ag⁺): [GAg-H₂O]⁺, [GAgG]⁺ produced in the electrospray ionization source of the mass spectrometer, and [GAg]⁺ produced by collision induced dissociation of the [GAgG]⁺ complex. On the basis of the comparison of theoretically calculated IR spectra, we show that there are two isomers of each complex containing two different keto-amino (KA) tautomers of G (GKA(1,9) and GKA(1,7)). The observed isomers are the most stable structures in aqueous solution, and their experimentally estimated relative populations are in better agreement with the calculated relative populations in solution than in the gas phase, both at 298 K. We concluded that these observations suggest that GKA(1,9) and GKA(1,7) coexist in solution according to previous theoretical reports (Colominas, C.; et al. J. Am. Chem. Soc. 1996, 118, 6811). We were unable to find any evidence of the presence of the GEA(9), GKA(3,7), GKA(3,9), or GKA(7,9), whose relative stabilities in solution are strongly dependent on the theoretical method used to account for the solvent effect (Hanus, M.; et al. J. Am. Chem. Soc. 2003, 125, 7678).

Variable temperature FTIR spectra of polycrystalline purine nucleobases and estimating strengths of individual hydrogen bonds

In the first part of this work, we report the FTIR spectra of pure NH and isotopically substituted ND (10-15% D and 80-90% D) polycrystalline hypoxanthine, xanthine, adenine and guanine recorded in the 400-4000 cm^-1 range, as a function of temperature (10-300 K). We provide assignments of the stretching and out-of-plane bending amine (NH₂) and imine (NH) bands to the distinct H-bonds present in the crystal, based on the temperature sensitivity and isotopic exchange behavior. Empirical correlations between spectral and thermodynamic or structural parameters enabled us to estimate the energies and lengths of H-bonds in the studied nucleobase crystals and to correlate them with literature data. The empirical H-bonding energies are compared with H-bonding and stacking energies computed for hypoxanthine. In the second part, strategies for using the empirical correlations together with information extracted from quantum mechanical data (in particular from the Bader's quantum theory of atoms in molecules, QTAIM) for the evaluation of hydrogen bonding properties are discussed, and their advantages and drawbacks pointed out. The justification for a cooperative use of quantum-mechanical calculations with empirical spectra-energy correlations is discussed.

The thermochemical, structural, and spectroscopic analyses of the tautomers of sulfur and selenium modified emissive nucleobases

In this study, density functional theory (DFT) and time dependent density functional theory (TD-DFT) are used to investigate the stabilities and spectral properties [IR, UV–vis, and two-photon absorption (2PA)] of two sets of modified RNA nucleobase tautomers. The modifications introduce either a sulfur or selenium atom to form an isothiazolo[4,3-d]pyrimidine or isoselenazolo[4,3-d]pyrimidine heterocyclic core, respectively. The relative stabilities of both sets of modified tautomers determined with B3LYP/6-31++G(d,p) reveal that in water (with a polarizable continuum model), the 6-keto-2-amino tautomer of guanine and the rare 4-imino-2-keto tautomer of cytosine may be present at significant populations, whereas the 6-enol-2-amino tautomer of guanine is more common in the gas phase. The identification of these modified tautomers is possible due to the natural differences in their vibrational modes and hence IR spectra. Furthermore, the photophysical properties of both these sets of modified tautomers indicate that excitation and emission energies are shifted relative to their more abundant form in both one photon absorption and emission and 2PA spectra, as determined at the B3LYP/6-31++G(d,p) and CAM-B3LYP/aug-cc-pVDZ levels of theory, respectively. Even though the 2PA cross sections in water for all of the species are small (0.3–2.3 GM), the modified cytosine tautomer shows promise, as its cross section is larger than the more dominant form. The spectral separation between the dominant form and the tautomers of isoselenazole and isothiazole modified cytosine and guanine are relatively similar, suggesting both modifications could be useful in elucidating the tautomers from their more abundant counterparts.

Recent Progress in the Manipulation of Molecules with DC Electric Fields

The structure and reactivity of a molecule in the condensed phase are governed by its intermolecular interactions with the surrounding environment. The multipole expansion of each molecule in the condensed phase indicates that the intermolecular interactions are essentially electrostatic (e.g., ion-dipole, dipole-dipole, dipole-quadrupole, dipole-induced dipole). The electrostatic field is a fundamental language of intermolecular communications. Therefore, understanding the influence of the electrostatic field on a molecule, that is, the mechanisms by which an electrostatic field manipulates a molecule, from the perspective of molecular structure, energy states, and dynamics is indispensable for illustrating and, by extension, controlling the chemistry in molecular systems.In this Account, we describe the recent progress made in manipulation of molecular processes using an external DC electrostatic field. An electrostatic field with unprecedentedly high strength (≤4 × 10⁸ V/m) was applied in a controlled manner across a molecular film sample using the ice film nanocapacitor method. This field strength is comparable in magnitude to that of weak intermolecular interactions such as van der Waals interactions in the condensed phases. The samples were prepared using a thin film growing technique in vacuum to obtain the desired chemically tailored molecular systems. The examples of prepared systems included small molecules and molecular clusters isolated in cryogenic Ar matrices, frozen molecular films in amorphous or crystalline phase, and interfaces of multilayered molecular films. The response of the molecules to the external field was monitored by reflection-absorption infrared spectroscopy. This approach allowed us to investigate a variety of molecular systems with various intermolecular strength and environments under the influence of strong electrostatic fields. The range of observed molecular behaviors includes the manipulation of molecular orientation, intramolecular dynamics, and proton transfer reactions as an example of stereodynamic control of chemical reactivity. These observations improve our understanding of molecular behaviors in strong electric fields and broaden our perspective on electrostatic manipulation of molecules. This information is also relevant to a variety of research topics in physical and biological sciences where electric fields play a role in molecular and biological functions.

Representing molecular ground and excited vibrational eigenstates with nuclear densities obtained from semiclassical initial value representation molecular dynamics

We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in the protonated glycine molecule [Aieta et al., Nat Commun 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in the correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, also to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond lengths, angles, and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching, the typical harmonic nodal pattern is absent in the anharmonic distribution.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Semiclassical Vibrational Spectroscopy of Biological Molecules Using Force Fields

Semiclassical spectroscopy is a practical way to get an accurately approximate quantum description of spectral features starting from ab initio molecular dynamics simulations. The computational bottleneck for the method is represented by the cost of ab initio potential, gradient, and Hessian matrix estimates. This drawback is particularly severe for biological systems due to their unique complexity and large dimensionality. The main goal of this manuscript is to demonstrate that quantum dynamics and spectroscopy, at the level of semiclassical approximation, are doable even for sizable biological systems. To this end, we investigate the possibility of performing semiclassical spectroscopy simulations when ab initio calculations are replaced by computationally cheaper force field evaluations. Both polarizable (AMOEBABIO18) and nonpolarizable (AMBER14SB) force fields are tested. Calculations of some particular vibrational frequencies of four nucleosides, i.e., uridine, thymidine, deoxyguanosine, and adenosine, show that ab initio simulations are accurate and widely applicable. Conversely, simulations based on AMBER14SB are limited to harmonic approximations, but those relying on AMOEBABIO18 yield acceptable semiclassical values if the investigated conformation has been included in the force field parametrization. The main conclusion is that AMOEBABIO18 may provide a viable route to assist semiclassical spectroscopy in the study of large biological molecules for which an ab initio approach is not computationally affordable.

UV-promoted radical formation, and near-IR-induced and spontaneous conformational isomerization in monomeric 9-methylguanine isolated in low-temperature Ar matrices

Three low-energy isomers of 9-methylguanine, the amino-oxo (AO) form and two amino-hydroxy (AH1 and AH2) conformers, were trapped from the gas phase into low-temperature argon matrices. The AH1 and AH2 isomers, differing in the orientation of the OH group, were found to transform into each other upon excitation with near-IR light. The population of the AO form of the compound was not changed upon any near-IR irradiation of the matrix samples. Using monochromatic near-IR light, generated by a frequency-tunable laser source, it was possible to selectively induce the AH1 → AH2 or AH2 → AH1 conversion. Photoreversibility of this conformational transformation was then demonstrated. Exposure of matrix-isolated monomers of 9-methylguanine to broadband near-IR light also led to conformational conversions within the amino-hydroxy tautomeric form; the final stage of this process was always the same photostationary state independent of the initial ratio of AH1 and AH2 populations. Spontaneous conformational conversion, transforming the higher-energy AH2 form into the lower-energy AH1 isomer, was observed for matrix-isolated monomers of 9-methylguanine kept in the dark. The mechanism of this process must rely on quantum tunneling of the light hydrogen atom. Irradiation of matrix-isolated 9-methylguanine with UV laser light at λ = 288 or 285 nm led to a substantial consumption of the two AH forms, while the amount of AO isomer remained unchanged. On the other hand, a decrease in the population of the AO isomer occurred upon excitations at shorter wavelengths, λ = 280 or 275 nm. The spectral changes observed after UV-irradiation suggest the generation (and stabilization in the matrix) of a radical species, resulting from the photocleavage of the O-H or N1-H bonds, in the AH or AO isomer, respectively.

Vibrational investigation of nucleobases by means of divide and conquer semiclassical dynamics

In this work, we report a computational study of the vibrational features of four different nucleobases employing the divide-and-conquer semiclassical initial value representation molecular dynamics method. Calculations are performed on uracil, cytosine, thymine, and adenine. Results show that the overall accuracy with respect to experiments is within 20 wavenumbers, regardless of the dimensionality of the nucleobase. Vibrational estimates are accurate even in the complex case of cytosine, where two relevant conformers are taken into account. These results are promising in the perspective of future studies on more complex systems, such as nucleotides or nucleobase pairs.

Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods

We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).

Photodynamics of alternative DNA base isoguanine

Pump–probe experiments and quantum-chemical simulations of UV-excited isoguanine elucidate its tautomer dependent photochemical properties.

Spatially Separating the Conformers of a Dipeptide

Atomic-resolution-imaging approaches for single molecules, such as coherent X-ray diffraction at free-electron lasers, require the delivery of high-density beams of identical molecules. However, even very cold beams of biomolecules typically have multiple conformational states populated. We demonstrate the production of very cold (T_rot ≈2.3 K) molecular beams of intact dipeptide molecules, which were then spatially separated into the individual populated conformational states. This is achieved using the combination of supersonic expansion laser-desorption vaporization with electrostatic deflection in strong inhomogeneous fields. This represents the first demonstration of a conformer-separated and rotationally cold molecular beam of a peptide, which enables the investigation of conformer-specific chemistry using inherently non-conformer-specific techniques. It furthermore represents a milestone toward the direct structural imaging of individual biological molecules with atomic resolution by ultrafast diffractive-imaging methods.

Nickel Adatoms Induced Tautomeric Dehydrogenation of Thymine Molecules on Au(111)

Tautomerization of nucleobases may induce base mismatches resulting in the abnormal disturbance of gene replication and expression, which has therefore attracted widespread interests in many disciplines. Metal atoms participating in a variety of important biological processes are found to be able to affect the nucleobase tautomerization as evidenced by many theoretical and spectroscopic studies. To get the real-space evidence and to unravel the underlying mechanism for the metal-induced tautomerization, especially from the keto form to the enol one, the interplay of high-resolution scanning tunneling microscopy imaging/manipulation and density functional theory (DFT) calculations has been employed. We present a process showing the Ni adatom-induced keto-enol tautomeric dehydrogenation of thymine molecules on Au(111). The key to making such a process feasible is the Ni atoms which greatly lower the energy barrier for the tautomerization from keto to enol form, which is rationalized by extensive DFT-based transition-state search calculations.

An experimental and theoretical investigation of XPS and NEXAFS of 5-halouracils

The C, N and O 1s excitation and ionization processes of 5X-uracil (X = F, Cl, Br, and I) were investigated using near edge X-ray absorption fine structure (NEXAFS) and X-ray photoemission (XPS) spectroscopies. This study aims at the fine assessment of the effects of the functionalization of uracil molecules by halogen atoms having different electronegativity and bound to the same molecular site. Two DFT-based approaches, which rely on different paradigms, have been used to simulate the experimental spectra and to assign the corresponding features. The analysis of the screening of the core holes of the different atoms via electronic charge density plots has turned out to be a useful tool to illustrate the competition between the partially aromatic and partially conjugate properties of this class of molecules.

How nature covers its bases

The response of nucleobases to UV radiation depends on structure in subtle ways, as revealed by gas-phase experiments.

Excited-State Dynamics of Isocytosine: A Hybrid Case of Canonical Nucleobase Photodynamics

We present resonant two-photon ionization (R2PI) spectra of isocytosine (isoC) and pump-probe results on two of its tautomers. IsoC is one of a handful of alternative bases that have been proposed in scenarios of prebiotic chemistry. It is structurally similar to both cytosine (C) and guanine (G). We compare the excited-state dynamics with the Watson-Crick (WC) C and G tautomeric forms. These results suggest that the excited-state dynamics of WC form of G may primarily depend on the heterocyclic substructure of the pyrimidine moiety, which is chemically identical to isoC. For WC isoC we find a single excited-state decay with a rate of ∼10¹⁰ s^-1, while the enol form has multiple decay rates, the fastest of which is 7 times slower than for WC isoC. The excited-state dynamics of isoC exhibits striking similarities with that of G, more so than with the photodynamics of C.

Excited State Dynamics of 6-Thioguanine

Here we present the excited state dynamics of jet-cooled 6-thioguanine (6-TG), using resonance-enhanced multiphoton ionization (REMPI), IR-UV double resonance spectroscopy, and pump-probe spectroscopy in the nanosecond and picosecond time domains. We report data on two thiol tautomers, which appear to have different excited state dynamics. These decay to a dark state, possibly a triplet state, with rates depending on tautomer form and on excitation wavelength, with the fastest rate on the order of 10¹⁰ s^-1. We also compare 6-TG with 9-enolguanine, for which we observed decay to a dark state with a 2 orders of magnitude smaller rate. At increased excitation energy (∼+500 cm^-1) an additional pathway appears for the predominant thiol tautomer. Moreover, the excited state dynamics for 6-TG thiols is different from that recently predicted for thiones.

Conformational structures of jet-cooled acetaminophen-water clusters: a gas phase spectroscopic and computational study

Jet-cooled acetaminophen (AAP)-water clusters, AAP-(H₂O)₁, were investigated by mass-selected resonant two-photon ionization (R2PI), ultraviolet-ultraviolet hole-burning (UV-UV HB), infrared-dip (IR-dip), and infrared-ultraviolet hole-burning (IR-UV HB) spectroscopy. Each syn- and anti-AAP rotamer has three distinctive binding sites (-OH, >CO, and >NH) for a water molecule, thus 6 different AAP-(H₂O)₁ conformers are expected to exist in the molecular beam. The origin bands of the AAP(OH)-(H₂O)₁ and AAP(CO)-(H₂O)₁ conformers (including their syn- and anti-conformers) in the R2PI spectrum are shifted to red and blue compared to those of the AAP monomer, respectively. These frequency shifts upon complexation between a water molecule and a specific binding site of AAP are also predicted by theoretical calculations. The spectral assignments of the origin bands in the R2PI spectra and the IR vibrational bands in the IR-dip spectra of the four lowest-energy conformers of AAP-(H₂O)₁, [syn- and anti-AAP(OH)-(H₂O)₁ and syn- and anti-AAP(CO)-(H₂O)₁], are aided by ab initio and time-dependent density functional theory (TDDFT) calculations. Further investigation of the IR-dip spectra has revealed a hydrogen-bonded NH stretching mode, supporting the presence of the syn-AAP(NH)-(H₂O)₁ conformer. Moreover, by employing IR-UV HB spectroscopy, we have reconfirmed the existence of the syn-AAP(NH)-(H₂O)₁ conformer, which happened to be buried underneath the broad background contributed by the AAP(OH)-(H₂O)₁ conformers. These observations have led us to conclude that all of the possible conformers of AAP-(H₂O)₁ have been found in this study.

Mechanisms and energetics for N-glycosidic bond cleavage of protonated 2'-deoxyguanosine and guanosine

Experimental and theoretical investigations suggest that hydrolysis of N-glycosidic bonds generally involves a concerted SN2 or a stepwise SN1 mechanism. While theoretical investigations have provided estimates for the intrinsic activation energies associated with N-glycosidic bond cleavage reactions, experimental measurements to validate the theoretical studies remain elusive. Here we report experimental investigations for N-glycosidic bond cleavage of the protonated guanine nucleosides, [dGuo+H](+) and [Guo+H](+), using threshold collision-induced dissociation (TCID) techniques. Two major dissociation pathways involving N-glycosidic bond cleavage, resulting in production of protonated guanine or the elimination of neutral guanine are observed in competition for both [dGuo+H](+) and [Guo+H](+). The detailed mechanistic pathways for the N-glycosidic bond cleavage reactions observed are mapped via electronic structure calculations. Excellent agreement between the measured and B3LYP calculated activation energies and reaction enthalpies for N-glycosidic bond cleavage of [dGuo+H](+) and [Guo+H](+) in the gas phase is found indicating that these dissociation pathways involve stepwise E1 mechanisms in analogy to the SN1 mechanisms that occur in the condensed phase. In contrast, MP2 is found to significantly overestimate the activation energies and slightly overestimate the reaction enthalpies. The 2'-hydroxyl substituent is found to stabilize the N-glycosidic bond such that [Guo+H](+) requires ∼25 kJ mol(-1) more than [dGuo+H](+) to activate the glycosidic bond.

Theoretical Investigation on the Substituent Effect of Halogen Atoms at the C8 Position of Adenine: Relative Stability, Vibrational Frequencies, and Raman Spectra of Tautomers

We have theoretically investigated the substituent effect of adenine at the C8 position with a substituent X = H, F, Cl, and Br by using the density functional theory (DFT) at the B3LYP/6-311+G(d, p) level. The aim is to study the substituent effect of halogen atoms on the relative stability, vibrational frequencies, and solvation effect of tautomers. Our calculated results show that for substituted adenine molecules the N9H8X tautomer to be the most stable structure in gas phase at the present theoretical level. Here N9H8X denotes the hydrogen atom binds to the N9 position of imidazole ring and X denotes H, F, Cl, and Br atoms. The influence of the induced attraction of the fluorine substituent is significantly larger than chlorine and bromine ones. The halogen substituent effect has a significant influence on changes of vibrational frequencies and Raman intensities.

Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases

The photophysics and photochemistry of DNA is of great importance due to the potential damage of the genetic code by UV light. Quantum mechanical studies have played a key role in interpretating the results of modern time-resolved pump-probe spectroscopy, and in elucidating the main photoactivated reactive paths. This review provides a concise, complete picture of the computational studies carried out, approximately, in the past decade. We start with an overview of the photophysics of the nucleobases in the gas phase and in solution. We discuss the proposed mechanisms for ultrafast decay to the ground state, that involve conical intersections, consider the role of triplet states, and analyze how the solvent modulates the photophysics. Then we move to larger systems, from dinucleotides to single- and double-stranded oligonucleotides. We focus on the possible role of charge transfer and delocalized or excitonic states in the photophysics of these systems and discuss the main photochemical paths. We finish with an outlook on the current challenges in the field and future directions of research.

Self-assembly of Natural and Unnatural Nucleobases at Surfaces and Interfaces

The self-assembly of small organic molecules interacting via non-covalent forces is a viable approach towards the construction of highly ordered nanostructured materials. Among various molecular components, natural and unnatural nucleobases can undergo non-covalent self-association to form supramolecular architectures with ad hoc structural motifs. Such structures, when decorated with appropriate electrically/optically active units, can be used as scaffolds to locate such units in pre-determined positions in 2D on a surface, thereby paving the way towards a wide range of applications, e.g., in optoelectronics. This review discusses some of the basic concepts of the supramolecular engineering of natural and unnatural nucleobases and derivatives thereof as well as self-assembly processes on conductive solid substrates, as investigated by scanning tunnelling microscopy in ultra-high vacuum and at the solid/liquid interface. By unravelling the structure and dynamics of these self-assembled architectures with a sub-nanometer resolution, a greater control over the formation of increasingly sophisticated functional systems is achieved. The ability to understand and predict how nucleobases interact, both among themselves as well as with other molecules, is extremely important, since it provides access to ever more complex DNA- and RNA-based nanostructures and nanomaterials as key components in nanomechanical devices.

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.

Isolated neutral peptides

This chapter examines the structural characterisation of isolated neutral amino-acids and peptides. After a presentation of the experimental and theoretical state-of-the-art in the field, a review of the major structures and shaping interactions is presented. Special focus is made on conformationally-resolved studies which enable one to go beyond simple structural characterisation; probing flexibility and excited-state photophysics are given as examples of promising future directions.

Understanding the ionic liquid [NC4111][NTf2] from individual building blocks: an IR-spectroscopic study

This study explores the interactions underlying the IR spectra of the ionic liquid [NC4111][NTf2] and its deuterated isotopomer [d9-NC4111][NTf2] by first isolating the spectra of charged ionic building blocks using mass-selective CIVP spectroscopy and then following the evolution of these bands upon sequential assembly of the ionic constituents. The spectra of the (1,1) and (2,2) neutral ion pairs are recorded using superfluid helium droplets as well as a solid neon matrix, while those of the larger charged aggregates are again obtained with CIVP. In general, the cluster spectra are similar to that of the bulk, with the (2,2) system displaying the closest resemblance. Analysis of the polarization-dependent band intensities of the neutral ion pairs in liquid droplets as a function of external electric field yields dipole moments of the neutral aggregates. This information allows a coarse assessment of the packing structure of the neutral pairs to be antiparallel at 0.37 K, in contrast to the parallel arrangement found for the assembly of small, high-dipole neutral molecules with large rotational constants (e.g., HCN). The role of an extra anion or cation attached to both the (1,1) and the (2,2) ion pairs to form the charged clusters is discussed in the context of an additional remote, more unfavorable binding site intrinsic to the nature of the charged IL clusters and as such not anticipated in the bulk phase. Whereas for the anion itself only the lowest energy trans conformer was observed, the higher clusters showed an additional population of the cis conformer. The interactions are found to be consistent with a minimal role of hydrogen bonding.

Automated quantum chemistry based molecular dynamics simulations of electron ionization induced fragmentations of the nucleobases Uracil, Thymine, Cytosine, and Guanine

The gas-phase decomposition pathways of electron ionization (EI)-induced radical cations of the nucleobases uracil, thymine, cytosine, and guanine are investigated by means of mixed quantum-classical molecular dynamics. No preconceived fragmentation channels are used in the calculations. The results compare well to a plethora of experimental and theoretical data for these important biomolecules. With our combined stochastic and dynamic approach, one can access in an unbiased way the energetically available decomposition mechanisms. Additionally, we are able to separate the EI mass spectra of different tautomers of cytosine and guanine. Our method (previously termed quantum chemistry electron ionization mass spectra) reproduces free nucleobase experimental mass spectra well and provides detailed mechanistic in-sight into high-energy unimolecular decomposition processes.

Experimental and first-principles investigation of the adsorption and entrapping of guanine with SiO2 clusters of sol–gel silicate material for understanding DNA photodamage

A sol–gel silicate matrix containing entrapped guanine was prepared. The SiO₂ matrix provides UVA protection by reducing the light penetration to the entrapped guanine molecules.

Electron Diffraction of Superfluid Helium Droplets

We present experimental results of electron diffraction of superfluid helium droplets and droplets doped with phthalocyanine gallium chloride and discuss the possibility of performing the same experiment with a laser aligned sample. The diffraction profile of pure droplets demonstrates dependence on the nozzle temperature, that is, on the average size of the droplets. Larger clusters demonstrate faster decay with increasing momentum transfer, whereas smaller clusters converge to isolated gas phase molecules at source temperatures of 18 K and higher. Electron diffraction of doped droplets shows similar modified molecular scattering intensity as that of the corresponding gas phase molecules. On the basis of fittings of the scattering profile, the number of remaining helium atoms of the doped droplets is estimated to be on the order of hundreds. This result offers guidance in assessing the possibility of electron diffraction from laser aligned molecules doped in superfluid helium droplets.

Atomic-scale investigation on the facilitation and inhibition of guanine tautomerization at Au(111) surface

Nucleobase tautomerization might induce mismatch of base pairing. Metals, involved in many important biophysical processes, have been theoretically proven to be capable of affecting tautomeric equilibria and stabilities of different nucleobase tautomers. However, direct real-space evidence on demonstrating different nucleobase tautomers and further revealing the effect of metals on their tautomerization at surfaces has not been reported to date. From the interplay of high-resolution STM imaging and DFT calculations, we show for the first time that tautomerization of guanine from G/9H to G/7H is facilitated on Au(111) by heating, whereas such tautomerization process is effectively inhibited by introducing Ni atoms due to its preferential coordination at the N7 site of G/9H. These findings may help to elucidate possible influence of metals on nucleobase tautomerization and provide from a molecular level some theoretical basis on metal-based drug design.

Coupled-cluster and density functional theory studies of the electronic 0–0 transitions of the DNA bases

The 0–0 transitions of the electronic excitation spectra of the lowest tautomers of the four nucleotide (DNA) bases have been studied using linear-response approximate coupled-cluster singles and doubles (CC2) calculations.

Experimental and first-principles study of guanine adsorption on ZnO clusters

Electronic structure for interaction of guanine with Zn₂O₂ cluster and the most preferred N1-site to form a stable G–Zn₂O₂ model.