Ab Initio Study of the Electronic Spectrum of Formamide with Explicit Solvent

Photochemistry of 2-thiooxazole: a plausible prebiotic precursor to RNA nucleotides

Potentially prebiotic chemical reactions leading to RNA nucleotides involve periods of UV irradiation, which are necessary to promote selectivity and destroy biologially irrelevant side products. Nevertheless, UV light has only been applied to promote specific stages of prebiotic reactions and its effect on complete prebiotic reaction sequences has not been extensively studied. Here, we report on an experimental and computational investigation of the photostability of 2-thiooxazole (2-TO), a potential precursor of pyrimidine and 8-oxopurine nucleotides on early Earth. Our UV-irradiation experiments resulted in rapid decomposition of 2-TO into unidentified small molecule photoproducts. We further clarify the underlying photochemistry by means of accurate ab initio calculations and surface hopping molecular dynamics simulations. Overall, the computational results show efficient rupture of the aromatic ring upon the photoexcitation of 2-TO via breaking of the C-O bond. Consequently, the initial stage of the divergent prebiotic synthesis of pyrimidine and 8-oxopurine nucleotides would require periodic shielding from UV light either with sun screening chromophores or through a planetary scenario that would protect 2-TO until it is transformed into a more stable intermediate compound, e.g. oxazolidinone thione.

Excited state dynamics of protonated dopamine: hydration and conformation effects

Electronic and vibrational spectroscopy in a cryogenic ion trap has been applied to protonated dopamine water clusters and assigned with the help of quantum chemistry calculations performed in the ground and electronic excited states. A dramatic hydration effect is observed when dopamine is solvated by three water molecules. The broad electronic spectra recorded for the bare and small water clusters containing protonated dopamine turn to sharp, well-resolved vibronic transitions in the 1-3 complex. This reflects the change induced by hydration in the photodynamics of protonated dopamine which is initially controlled by an excited state proton transfer (ESPT) reaction from the ammonium group toward the catechol ring. Interestingly, conformer selectivity is revealed in the 1-3 complex which shows two low lying energy conformers for which the ESPT reaction is prevented or not depending on the H-bond network formed between the dopamine and water molecules.

IR-VUV spectroscopy of pyridine dimers, trimers and pyridine-ammonia complexes in a supersonic jet

The infrared spectra of the C-H stretching vibrations of (pyridine)m, m = 1-3, and the N-H stretching vibrations of (pyridine)m-(NH3)n, m = 1, 2; n = 1-4, complexes were investigated by infrared (IR)-vacuum ultraviolet (VUV) spectroscopy under jet-cooled conditions. The ionization potential (IP0) of the pyridine monomer was determined to be 74 546 cm-1 (9.242 eV), while its complexes showed only smooth curves of the ionization thresholds at ∼9 eV, indicating large structural changes in the ionic form. The pyridine monomer exhibits five main features with several satellite bands in the C-H stretching region at 3000-3200 cm-1. Anharmonic calculations including Fermi-resonance were carried out to analyze the candidates of the overtone and combination bands which can couple to the C-H stretching fundamentals. For (pyridine)2 and (pyridine)3, most C-H bands are blue-shifted by 3-5 cm-1 from those of the monomer. The structures revealed by random searching algorithms with density functional methods indicate that the π-stacked structure is most stable for (pyridine)2, while (pyridine)3 prefers the structures stabilized by dipole-dipole and C-Hπ interactions. For the (pyridine)m-(NH3)n complexes, the mass spectrum exhibited a wide range distribution of the complexes. The observed IR spectra in the N-H stretching vibrations of the complexes showed four main bands in the 3200-3450 cm-1 region. These features are very similar to those of (NH3)n complexes, and the bands are assigned to the anti-symmetric N-H stretching band (ν3), the symmetric N-H stretching (ν1) band, and the first overtone bands of the N-H bending vibrations (2ν4). The anharmonic calculations including the Fermi-resonance between ν1 and 2ν4 well reproduced the observed spectra.

Sequential electron transfer governs the UV-induced self-repair of DNA photolesions

Cyclobutane pyrimidine dimers (CpDs) are among the most common DNA lesions occurring due to the interaction with ultraviolet light. While photolyases have been well known as external factors repairing CpDs, the intrinsic self-repairing capabilities of the GAT 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 T DNA sequence were discovered only recently and are still largely obscure. Here, we elucidate the mechanistic details of this self-repair process by means of MD simulations and QM/MM computations involving the algebraic diagrammatic construction to the second order [ADC(2)] method. We show that local UV-excitation of guanine may be followed by up to three subsequent electron transfers, which may eventually enable efficient CpD ring opening when the negative charge resides on the T 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 T dimer. Consequently, the molecular mechanism of GAT 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 T self-repair can be envisaged as sequential electron transfer (SET) occurring downhill along the slope of the S₁ potential energy surface. Even though the general features of the SET mechanism are retained in both of the studied stacked conformers, our optimizations of different S₁/S₀ state crossings revealed minor differences which could influence their self-repair efficiencies. We expect that such assessment of the availability and efficiency of the SET process in other DNA oligomers could hint towards other sequences exhibiting similar photochemical properties. Such explorations will be particularly fascinating in the context of the origins of biomolecules on Earth, owing to the lack of external repairing factors in the Archean age.

Curve crossing in a manifold of coupled electronic states: direct quantum dynamics simulations of formamide

Quantum dynamics simulations are an important tool to evaluate molecular behaviour including the, often key, quantum nature of the system. In this paper we present an algorithm that is able to simulate the time evolution of a molecule after photo-excitation into a manifold of states. The direct dynamics variational multi-configurational Gaussian (DD-vMCG) method circumvents the computational bottleneck problems of traditional grid-based methods by computing the potential energy functions on-the-fly, i.e. only where required. Unlike other commonly used direct dynamics methods, DD-vMCG is fully quantum mechanical. Here, the method is combined with a novel on-the-fly diabatisation scheme to simulate the short-time dynamics of the key molecule formamide and its acid analogue formimidic acid. This is a challenging test system due to the nature and large number of excited states, and eight coupled states are included in the calculations. It is shown that the method is able to provide unbiased information on the product channels open after excitation at different energies and demonstrates the potential to be a practical scheme, limited mainly by the quality of the quantum chemistry used to describe the excited states.

A prebiotically plausible synthesis of pyrimidine β-ribonucleosides and their phosphate derivatives involving photoanomerization

Previous research has identified ribose aminooxazoline as a potential intermediate in the prebiotic synthesis of the pyrimidine nucleotides with remarkable properties. It crystallizes spontaneously from reaction mixtures, with an enhanced enantiomeric excess if initially enantioenriched, which suggests that reservoirs of this compound might have accumulated on the early Earth in an optically pure form. Ribose aminooxazoline can be converted efficiently into α-ribocytidine by way of 2,2'-anhydroribocytidine, although anomerization to β-ribocytidine by ultraviolet irradiation is extremely inefficient. Our previous work demonstrated the synthesis of pyrimidine β-ribonucleotides, but at the cost of ignoring ribose aminooxazoline, using arabinose aminooxazoline instead. Here we describe a long-sought route through ribose aminooxazoline to the pyrimidine β-ribonucleosides and their phosphate derivatives that involves an extraordinarily efficient photoanomerization of α-2-thioribocytidine. In addition to the canonical nucleosides, our synthesis accesses β-2-thioribouridine, a modified nucleoside found in transfer RNA that enables both faster and more-accurate nucleic acid template-copying chemistry.

Oxocomplexes of Mo(VI) and W(VI) with 8-hydroxyquinoline-5-sulfonate in solution: structural studies and the effect of the metal ion on the photophysical behaviour

Multinuclear ((1)H, (13)C, (95)Mo and (183)W) NMR spectroscopy, combined with DFT calculations, provides detailed information on the complexation between the Mo(VI) and W(VI) oxoions and 8-hydroxyquinoline-5-sulfonate (8-HQS) in aqueous solution. Over the concentration region studied, Mo(VI) and W(VI) oxoions form three homologous complexes with 8-HQS in water in the pH range 2-8. Two of these, detected at pH < 6, are mononuclear 1 : 2 (metal : ligand) isomers, with the metal centre (MO2(2+)) coordinated to two 8-HQS ligands. An additional complex, dominant at slightly higher pH values (5-8) for solutions with a 1 : 1 metal : ligand molar ratio, has a binuclear M2O5(2+) centre coordinated to two 8-HQS ligands. The two metal atoms are bridged by three oxygen atoms, two coming from 8-HQS, together with the M-O-M bridge of the bimetallic centre. We show that the long-range exchange corrected BOP functional with local response dispersion (LCBOPLRD), together with explicit solvent molecules, leads to geometries that readily converge to equilibrium structures having realistic bridging O8-HQS-M bonds. Previous attempts to calculate the structures of such binuclear complexes using DFT with the B3LYP functional have failed due to difficulties in treating the weak interaction in these bridged structures. We believe that the LCBOPLRD method may be of more general application in theoretical studies in related binuclear metal complexes. UV/visible absorption and luminescence spectra of all the complexes have also been recorded. The complex between Mo(vi) and 8-HQS is only weakly luminescent, in contrast to what has been observed with this ligand and many other metal ions. We suggest that this is due to the presence of low-lying ligand-to-metal charge transfer (LMCT) states close to the emitting ligand-based level which quench the emission. However, with W(VI), DFT calculations show that the LMCT states are now much higher in energy than the ligand based levels, leading to a marked increase in fluorescence.

Physicochemical properties of the ternary complexes of Pt(ii) with uracil and small peptide moieties: an experimental and computational study

Interactions of Pt(Ura)(GL) with DNA.

Photodynamic behavior of electronic coupling in a N-methylformamide dimer

The role of the bridging water molecules has been studied during the excited state photodynamics of a N-methylformamide dimer in complex with water molecules employing the complete active space self-consistent field (CASSCF) and CAS perturbation theory (CASPT2) methods.

Quantum mechanical investigations on the role of C-terminal residue in influencing the structural features of dipeptides containing N-terminal proline

This study investigates the influence of the side chain moiety of C-terminal residue on the structural and molecular properties of seven dipeptides having proline at their N-terminal positions. The C-terminal component of the dipeptides is varied with seven different combinations viz. Ala, Leu, Asp, Thr, Asn, Arg and Sec. The calculations are carried out using B3LYP/6-311++G(d,p) level of theory in gas and implicit aqueous phase. Effects of explicit aqueous environment on the dipeptide structures are also investigated for two systems. The results furnished by this DFT study provide valuable information regarding the role of the side chain groups of C-terminal residues in determining the structural features of the amide planes, values of the ψ and ф dihedrals, geometry about the α-carbon atoms, theoretical IR spectra as well as the number and type of intramolecular H-bond interactions existing in the dipeptides, and extend a fine corroboration to the earlier theoretical and experimental observations. In aqueous phase the dipeptide geometries exhibit larger values of total dipole moments, greater HOMO-LUMO energy gaps and enhanced thermodynamic stability than those in gas phase. The explicit water molecules are found to modify the geometrical parameters related to the amide planes and vibrational spectra of the dipeptides.

Solvent induced shifts in the UV spectrum of amides

Solvent effects on the electronic spectra of formamide and trans-N-methylacetamide are studied using four different levels of theory: singly excited configuration interaction (CIS), equations of motion coupled-cluster theory with singles and doubles (EOM-CCSD), completely renormalized coupled-cluster theory with singles and doubles with perturbative triple excitations (CR-EOM-CCSD(T)), and time-dependent density functional theory (TDDFT), employing small clusters of water molecules. The simulated electronic spectrum is obtained via molecular dynamics simulations with 100 waters modeled with the effective fragment potential method and exhibits a blue-shift and red-shift, respectively, for the n → π* and πnb → π* vertical excitation energies, in good agreement with the experimental electronic spectra of amides.

NMR, DFT and luminescence studies of the complexation of Al(III) with 8-hydroxyquinoline-5-sulfonate

Multinuclear ((1)H, (13)C and (27)Al) magnetic resonance spectroscopy (1D and 2D), DFT calculations and fluorescence have been used to study the complexation of 8-hydroxyquinoline-5-sulfonate (8-HQS) with Al(III). The study combines the high sensitivity of luminescence techniques, the selectivity of multinuclear NMR spectroscopy with the structural details accessible through DFT calculations, and aims to provide a detailed understanding of the complexation between the Al(3+) ion and 8-HQS. A full speciation study has been performed and over the concentration region studied, the Al(3+) ion forms complexes with 8-HQS in an aqueous solution in the pH range 2-6. At higher pH, the extensive hydrolysis of the metal limits complexation. Using Job's method, three complexes were detected, with 1 : 1, 1 : 2 and 1 : 3 (metal : ligand) stoichiometries. These results are in agreement with those previously reported using potentiometric and electrochemical techniques. The geometries of the complexes are proposed based on the combination of NMR results with optimized DFT calculations. All the complexes in aqueous solutions at 25 °C are mononuclear species, and have an approximately octahedral geometry with the metal coordinated to one molecule of 8-HQS and four molecules of water (1 : 1 complex), two molecules of 8-HQS and two molecules of water mutually cis (1 : 2 complex), and to three molecules of 8-HQS in non-symmetrical arrangement (mer-isomer), for the 1 : 3 (metal : ligand) complex. On binding to Al(III), 8-HQS shows a more marked fluorescence than the weakly fluorescent free ligand. In addition, as previously noted, there are marked changes in the absorption spectra, which support the use of 8-HQS as a sensitive optical sensor to detect Al(3+) metal ions in surface waters and biological fluids. These complexes also show potential for applications in organic light emitting diodes (OLEDs).

Modeling the infrared and circular dichroism spectroscopy of a bridged cyclic diamide

Density functional theory based molecular dynamics simulations are used to study the structure, infrared (IR) spectroscopy, circular dichroism (CD) spectroscopy, and coupling between the amide I vibrations of a bridged cyclic diamide in the gas phase and in aqueous solution. IR spectra computed via the dipole moment time correlation function show a large red-shift of 30 cm(-1) in the amide I vibration in solution compared to the gas phase, and are in good agreement with experiment. Conformationally averaged CD spectra computed using the CIS(D) method are highly sensitive to the structures used, and structures sampled in the aqueous phase simulation are required to obtain qualitatively correct CD spectra. Analysis of the coupling between the amide I modes shows that in the aqueous phase there is an increased localization of the vibrations on the individual peptide groups and a reduction in the mode coupling parameter compared to the gas phase. Overall, the results illustrate the significance of incorporating molecular dynamics in the simulation of IR and CD spectra.

Lewis acid-base adducts: a quantitative Raman analysis of formamide and dimethylsulfoxide mixtures

Raman spectra of pure liquid dimethylsulfoxide (DMSO) and of binary mixtures of formamide (FA) and DMSO in different compositions were obtained. The vibrations involving the SO functional group in the band envelope at ca. 1050 cm(-1) of pure liquid DMSO are assigned to monomers, dimers and higher aggregates of DMSO. The appearance of a new band at 1024 cm(-1), whose intensity shows large dependence on the FA concentration, is assigned to a FA-DMSO adduct. This has been possible due to the two H-bond donor sites of FA and the strong donor character of DMSO that become the environment propitious for the donor-acceptor reaction. Quantitative analysis performed in the SO stretching region in the binary mixtures gives a 1:1 stoichiometry in this adduct in the limit of infinite dilution. This proportion is in full agreement with our previous determination for the FA-ACN adduct. The experimental evidence of the 1:1 FA-DMSO adduct is presented for the first time using Raman spectroscopy. The results described here open new possibilities to study the acid-base reactions nature of FA adducts.

Two-state model based on the block-localized wave function method

The block-localized wave function (BLW) method is a variant of ab initio valence bond method but retains the efficiency of molecular orbital methods. It can derive the wave function for a diabatic (resonance) state self-consistently and is available at the Hartree-Fock (HF) and density functional theory (DFT) levels. In this work we present a two-state model based on the BLW method. Although numerous empirical and semiempirical two-state models, such as the Marcus-Hush two-state model, have been proposed to describe a chemical reaction process, the advantage of this BLW-based two-state model is that no empirical parameter is required. Important quantities such as the electronic coupling energy, structural weights of two diabatic states, and excitation energy can be uniquely derived from the energies of two diabatic states and the adiabatic state at the same HF or DFT level. Two simple examples of formamide and thioformamide in the gas phase and aqueous solution were presented and discussed. The solvation of formamide and thioformamide was studied with the combined ab initio quantum mechanical and molecular mechanical Monte Carlo simulations, together with the BLW-DFT calculations and analyses. Due to the favorable solute-solvent electrostatic interaction, the contribution of the ionic resonance structure to the ground state of formamide and thioformamide significantly increases, and for thioformamide the ionic form is even more stable than the covalent form. Thus, thioformamide in aqueous solution is essentially ionic rather than covalent. Although our two-state model in general underestimates the electronic excitation energies, it can predict relative solvatochromic shifts well. For instance, the intense pi-->pi* transition for formamide upon solvation undergoes a redshift of 0.3 eV, compared with the experimental data (0.40-0.5 eV).