Ultrafast damage following radiation-induced oxidation of uracil in aqueous solution

Ultrafast fragmentation of highly-excited doubly-ionized deoxyribose: role of the liquid water environment

Ab initio molecular dynamics simulations are used to investigate the fragmentation dynamics following the double ionization of 2-deoxy-D-ribose (DR), a major component in the DNA chain. Different ionization scenarios are considered to provide a complete picture. First focusing on isolated DR2+, fragmentation patterns are determined for the ground electronic state, adding randomly distributed excitation energy to the nuclei. These patterns differ for the two isomers studied. To compare thermal and electronic excitation effects, Ehrenfest dynamics are also performed, allowing to remove the two electrons from selected molecular orbitals. Two intermediate-energy orbitals, localized on the carbon chain, were selected. The dissociation pattern corresponds to the most frequent pattern obtained when adding thermal excitation. On the contrary, targeting the four deepest orbitals, localized on the oxygen atoms, leads to selective ultrafast C-O and/or O-H bond dissociation. To probe the role of environment, a system consisting of a DR molecule embedded in liquid water is then studied. The two electrons are removed from either the DR or the water molecules directly linked to the sugar through hydrogen bonds. Although the dynamics onset is similar to that of isolated DR when removing the same deep orbitals localized on the sugar oxygen atoms, the subsequent fragmentation patterns differ. Sugar damage also occurs following the Coulomb explosion of neighboring H2O2+ molecules due to interaction with the emitted O or H atoms.

Pathways of two-body dissociation of BrCNq+ (q = 2, 3) induced by electron collision

Pathways of two-body fragmentation of BrCNq+ (q = 2, 3) have been explored by combined experimental and theoretical studies. In the experiment, the BrCN molecule is ionized by 1 keV electron impact and the created fragment ions are detected using an ion momentum imaging spectrometer. Six two-body fragmentation channels are identified. By measuring the momentum vectors of the fragment ions, the kinetic energy release (KER) distributions for these channels have been determined. Theoretically, the potential energy curves of BrCNq+ (q = 2, 3) as a function of Br-C and C-N internuclear distances are calculated by the complete active space self-consistent field method. By comparing the measured KER and theoretical predictions, pathways for the fragmentation channels are assigned. The relative branching ratios of the channels are also determined.

X-ray induced luminescence spectroscopy for DNA damaging intermediates aided by a monochromatic synchrotron radiation

PURPOSE

To identify the bonding sites of initial radiation interaction with DNA and to trace the following chemical reaction sequences on the pathway of damage induction, we carry out a spectroscopy XIL (X-ray induced luminescence) using soft X-ray synchrotron radiation. This is a nondestructive analysis of the excited intermediate species produced in a molecular mechanism on the damage induction pathway.

MATERIALS AND METHODS

We introduce aqueous samples of UMP (uridine-5'-monophosphate) in the vacuum by the use of a liquid micro-jet technique. The luminescence in the region of UV-VIS (from visible to ultraviolet) radiation induced after the absorption of monochromatic soft X-ray by aqueous UMP is measured with sweeping the soft X-ray energy in the region of 370-560 eV.

RESULTS

The enhanced XIL intensities for aqueous UMP in the region of soft X-ray of 410-530 eV (in "water window" region) are obtained. The enhancement of XIL intensities in the UV-VIS region, relative to the water control, is explained by the excitation and ionization of a K-shell electron of nitrogen atoms in the uracil moiety. The enhanced XIL intensities do not match the structure of XANES (X-ray absorption near-edge structure) of the aqueous UMP. This suggests that the XIL intensities reflect the quantum yields of luminescence, or the quantum yields for conversion by UMP of an absorbed X-ray into UV-VIS radiation. In this paper, spectra of luminescence are shown to be resolved by combining low pass filters. The filtered luminescence spectra are obtained at the center of gravity (λc) of the band pass wavelength regions at λc = 270nm, 295 nm, 340 nm, 385 nm, 450 nm, and 525 nm., which show a trend similar to the fluorescence of nucleobases induced by ultraviolet radiation.

CONCLUSION

It is concluded that the origin of the observed XIL is the hydrated uracil moiety in aqueous UMP, decomposition of which is suppressed by the migration of excess charge and internal energy after the double ionization due to Auger decay.

Suppression of X-ray-Induced Radiation Damage to Biomolecules in Aqueous Environments by Immediate Intermolecular Decay of Inner-Shell Vacancies

The predominant reason for the damaging power of high-energy radiation is multiple ionization of a molecule, either direct or via the decay of highly excited intermediates, as, e.g., in the case of X-ray irradiation. Consequently, the molecule is irreparably damaged by the subsequent fragmentation in a Coulomb explosion. In an aqueous environment, however, it has been observed that irradiated molecules may be saved from fragmentation presumably by charge and energy dissipation mechanisms. Here, we show that the protective effect of the environment sets in even earlier than hitherto expected, namely immediately after single inner-shell ionization. By combining coincidence measurements of the fragmentation of X-ray-irradiated microsolvated pyrimidine molecules with theoretical calculations, we identify direct intermolecular electronic decay as the protective mechanism, outrunning the usually dominant Auger decay. Our results demonstrate that such processes play a key role in charge delocalization and have to be considered in investigations and models on high-energy radiation damage in realistic environments.

First-Principles Simulations of Biological Molecules Subjected to Ionizing Radiation

Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.

Hydration of Nucleobase as Probed by Electron Emission of Uridine-5′-Mono-Phosphate (UMP) in Aqueous Solution Induced by Nitrogen K-Shell Ionization

To identify the precise early radiation processes of DNA lesions, we measure electron kinetic energy spectra emitted from uridine-5′ monophosphate (UMP) in aqueous solution for the photoionization of the N 1s orbital electron and for the following Auger effect using a monochromatic soft X-ray synchrotron radiation at energies above the nitrogen K-shell ionization threshold. The change of photoelectron spectra for UMP in aqueous solutions at different proton concentrations (pH = 7.5 and 11.3) is ascribed to the chemical shift of the N3 nitrogen atom in uracil moiety of canonical and deprotonated forms. The lowest double ionization potentials for aqueous UMP at different pH obtained from the Auger electron spectra following the N 1s photoionization values show the electrostatic aqueous interaction of uracil moiety of canonical (neutral) and deprotonated (negatively charged) forms with hydrated water molecules.

The physical stage of radiolysis of solvated DNA by high-energy-transfer particles: insights from new first principles simulations

Electron dynamics simulations based on density functional theory are carried out on nanometric molecular systems to decipher the primary processes following irradiation of bio-macromolecules by high energy transfer charged particles.

Ab Initio Molecular Dynamics Simulations to Interpret the Molecular Fragmentation Induced in Deoxyribose by Synchrotron Soft X-Rays

It has been suggested that core ionization in DNA atoms could induce complex, irreparable damage. Synchrotron soft X-rays have been used to probe the damage induced by such events in thin films of DNA components. In a complementary approach, we investigate the fragmentation dynamics following a carbon or oxygen K-shell ionization in 2-deoxy-D-ribose (DR), a major component in the DNA chain. Core ionization of the sugars hydration layer is also studied. To that aim, we use state-of-the-art ab initio Density Functional Theory-based Molecular Dynamics (MD) simulations. The ultrafast dissociation dynamics of the core ionized molecule, prior Auger decay, is modeled for about 10 fs. We show that the core-ionization of oxygen atoms within DR or its hydration layer may induce proton transfers towards nearby molecules, before Auger decay. In a second step, we model an Auger effect occurring either at the beginning or at the end of the core–hole dynamics. Two electrons are removed from the deepest valence molecular orbitals localized on the initially core-ionized oxygen atom (O*), and this electronic state is propagated by means of Ehrenfest MD. We show an ultrafast dissociation of the DR2+ molecule C-O* bonds, which, in most cases, seems independent of the time at which Auger decay occurs.

Fragmentation dynamics of nitrogen trifluoride induced by electron collision

The fragmentation dynamics of nitrogen trifluoride (NF₃) in collisions with a 500 eV electron is studied by using a momentum imaging spectrometer. The kinetic energy releases of two-body, three-body, and four-body fragmentation channels of NF₃ ^q+ (q = 2, 3) are investigated. The fragmentation dynamics of three-body, as well as four-body, dissociation channels is analyzed by the Dalitz plot and the Newton diagram. It is found that for all of the dissociation channels, the fragment including N atom (ion) always shares significant momenta, regardless of whether it is charged. For F atom, however, it is always emitted with negligible momenta.

Proton Collision on Deoxyribose Originating from Doubly Ionized Water Molecule Dissociation

In this work, we studied the fragmentation dynamics of 2-deoxy-d-ribose (DR) in solution that arises from the double ionization of a water molecule in its primary hydration shell. This process was modeled in the framework of ab initio molecular dynamics. The charge unbalanced in the solvent molecules produces a Coulomb explosion with the consequent release of protons with kinetic energy in the few electronvolts range, which collide with the surrounding molecules in solution inducing further chemical reactions. In particular, we observe proton collisions with the solute molecule DR, which leads to a complete ring opening. In DNA, damage to the DR moiety may lead to DNA strand breaking. This mechanism can be understood as one of the possible steps in the radiation-induced fragmentation of DNA chains.

Assessment of Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) in Radiation Chemistry: Ionized Water Dimer

Ionization in the condensed phase and molecular clusters leads to a complicated chain of processes with coupled electron-nuclear dynamics. It is difficult to describe such dynamics with conventional nonadiabatic molecular dynamics schemes since the number of states swiftly increases as the molecular system grows. It is therefore attractive to use a direct electron and nuclear propagation such as the real-time time-dependent density functional theory (RT-TDDFT). Here we report a RT-TDDFT benchmark study on simulations of singly and doubly ionized states of a water monomer and dimer as a prototype for more complex processes in a condensed phase. We employed the RT-TDDFT based Ehrenfest molecular dynamics with a generalized gradient approximate (GGA) functional and compared it with wave-function-based surface hopping (SH) simulations. We found that the initial dynamics of a singly HOMO ionized water dimer is similar for both the RT-TDDFT/GGA and the SH simulations but leads to completely different reaction channels on a longer time scale. This failure is attributed to the self-interaction error in the GGA functionals and it can be avoided by using hybrid functionals with large fraction of exact exchange (represented here by the BHandHLYP functional). The simulations of doubly ionized states are reasonably described already at the GGA level. This suggests that the RT-TDDFT/GGA method could describe processes following the autoionization processes such as Auger emission, while its applicability to more complex processes such as intermolecular Coulombic decay remains limited.

Quantum modeling of ultrafast photoinduced charge separation

Phenomena involving electron transfer are ubiquitous in nature, photosynthesis and enzymes or protein activity being prominent examples. Their deep understanding thus represents a mandatory scientific goal. Moreover, controlling the separation of photogenerated charges is a crucial prerequisite in many applicative contexts, including quantum electronics, photo-electrochemical water splitting, photocatalytic dye degradation, and energy conversion. In particular, photoinduced charge separation is the pivotal step driving the storage of sun light into electrical or chemical energy. If properly mastered, these processes may also allow us to achieve a better command of information storage at the nanoscale, as required for the development of molecular electronics, optical switching, or quantum technologies, amongst others. In this Topical Review we survey recent progress in the understanding of ultrafast charge separation from photoexcited states. We report the state-of-the-art of the observation and theoretical description of charge separation phenomena in the ultrafast regime mainly focusing on molecular- and nano-sized solar energy conversion systems. In particular, we examine different proposed mechanisms driving ultrafast charge dynamics, with particular regard to the role of quantum coherence and electron-nuclear coupling, and link experimental observations to theoretical approaches based either on model Hamiltonians or on first principles simulations.

Reactivity of prehydrated electrons toward nucleobases and nucleotides in aqueous solution

DNA damage induced via dissociative attachment by low-energy electrons (0 to 20 eV) is well studied in both gas and condensed phases. However, the reactivity of ultrashort-lived prehydrated electrons ([Formula: see text]) with DNA components in a biologically relevant environment has not been fully explored to date. The electron transfer processes of [Formula: see text] to the DNA nucleobases G, A, C, and T and to nucleosides/nucleotides were investigated by using 7-ps electron pulse radiolysis coupled with pump-probe transient absorption spectroscopy in aqueous solutions. In contrast to previous results, obtained by using femtosecond laser pump-probe spectroscopy, we show that G and A cannot scavenge [Formula: see text] at concentrations of ≤50 mM. Observation of a substantial decrease of the initial yield of hydrated electrons ([Formula: see text]) and formation of nucleobase/nucleotide anion radicals at increasing nucleobase/nucleotide concentrations present direct evidence for the earliest step in reductive DNA damage by ionizing radiation. Our results show that [Formula: see text] is more reactive with pyrimidine than purine nucleobases/nucleotides with a reactivity order of T > C > A > G. In addition, analyses of transient signals show that the signal due to formation of the resulting anion radical directly correlates with the loss of the initial [Formula: see text] signal. Therefore, our results do not agree with the previously proposed dissociation of transient negative ions in nucleobase/nucleotide solutions within the timescale of these experiments. Moreover, in a molecularly crowded medium (for example, in the presence of 6 M phosphate), the scavenging efficiency of [Formula: see text] by G is significantly enhanced. This finding implies that reductive DNA damage by ionizing radiation depends on the microenvironment around [Formula: see text].

M3C: A Computational Approach To Describe Statistical Fragmentation of Excited Molecules and Clusters

The Microcanonical Metropolis Monte Carlo method, based on a random sampling of the density of states, is revisited for the study of molecular fragmentation in the gas phase (isolated molecules, atomic and molecular clusters, complex biomolecules, etc.). A random walk or uniform random sampling in the configurational space (atomic positions) and a uniform random sampling of the relative orientation, vibrational energy, and chemical composition of the fragments is used to estimate the density of states of the system, which is continuously updated as the random sampling populates individual states. The validity and usefulness of the method is demonstrated by applying it to evaluate the caloric curve of a weakly bound rare gas cluster (Ar₁₃), to interpret the fragmentation of highly excited small neutral and singly positively charged carbon clusters (C_n, n = 5,7,9 and C_n⁺, n = 4,5) and to simulate the mass spectrum of the acetylene molecule (C₂H₂).

Nonadiabatic molecular dynamics simulations: synergies between theory and experiments

Recent developments in nonadiabatic dynamics enabled ab inito simulations of complex ultrafast processes in the condensed phase. These advances have opened new avenues in the study of many photophysical and photochemical reactions triggered by the absorption of electromagnetic radiation. In particular, theoretical investigations can be combined with the most sophisticated femtosecond experimental techniques to guide the interpretation of measured time-resolved observables. At the same time, the availability of experimental data at high (spatial and time) resolution offers a unique opportunity for the benchmarking and the improvement of those theoretical models used to describe complex molecular systems in their natural environment. The established synergy between theory and experiments can produce a better understanding of new ultrafast physical and chemical processes at atomistic scale resolution. Furthermore, reliable ab inito molecular dynamics simulations can already be successfully employed as predictive tools to guide new experiments as well as the design of novel and better performing materials. In this paper, I will give a concise account on the state of the art of molecular dynamics simulations of complex molecular systems in their excited states. The principal aim of this approach is the description of a given system of interest under the most realistic ambient conditions including all environmental effects that influence experiments, for instance, the interaction with the solvent and with external time-dependent electric fields, temperature, and pressure. To this end, time-dependent density functional theory (TDDFT) is among the most efficient and accurate methods for the representation of the electronic dynamics, while trajectory surface hopping gives a valuable representation of the nuclear quantum dynamics in the excited states (including nonadiabatic effects). Concerning the environment and its effects on the dynamics, the quantum mechanics/molecular mechanics (QM/MM) approach has the advantage of providing an atomistic (even though approximated) description of the solvent molecules, which is crucial for the characterization of all ultrafast relaxation phenomena that depend on the geometrical arrangement at the interface between a molecule and the solvent, for example, the hydrogen bond network. After a short description of the TDDFT-based implementation of Ehrenfest and trajectory surface hopping dynamics, I will present applications in different domains of molecular chemistry and physics: the analysis and the understanding of (time-resolved) X-ray absorption spectra, the interpretation of the ultrafast relaxation dynamics of photoexcited dyes in solution, and the design of specific laser pulses (capable of inducing desired chemical reactions) using local control theory.

Investigation of the fragmentation of core-ionised deoxyribose: a study as a function of the tautomeric form

Fragmentation dynamics following the core ionisation of isolated 2-deoxy-d-ribose by DFT-MD simulations.

Looking at radiation damage on prebiotic building blocks

A number of complex organic molecules have been detected in the interstellar medium, as well as in meteorites or comets. Among them, some exobiologic-relevant molecules have attracted particular interest. In the hypothesis of an exogen transport of prebiotic building blocks at the origin of life, the survival of such species and particularly their resistance to the solar UV radiation or cosmic rays is a key issue. For that purpose, we have performed a theoretical approach of the charge transfer dynamics induced by collision of protons with nucleobases and the 2-deoxy-d-ribose sugar moiety in a wide collision energy range. Calculations have been carried out by means of ab initio quantum chemistry molecular methods and compared to previous theoretical results using carbon projectile ions. Qualitative trends can be exhibited on DNA or RNA building blocks damage, which may concern studies on prebiotic species under spatial radiation.