Ionization and Electron Attachment for Nucleobases in Water

Charge Transfer Mechanism in Guanine-Based Self-Assembled Monolayers on a Gold Surface

In this work, we have theoretically determined the one-electron oxidation potentials and charge transfer mechanisms in complex systems based on a self-assembled monolayer of guanine molecules adsorbed on a gold surface through different organic linkers. Classical molecular dynamics simulations were carried out to sample the conformational space of both the neutral and the cationic species. Thus, the redox potentials were determined for the ensembles of geometries through multiscale quantum-mechanics/molecular-mechanics/continuum solvation model calculations in the framework of the Marcus theory and in combination with an additive scheme previously developed. In this context, conformational sampling, description of the environment, and effects caused by the linker have been considered. Applying this methodology, we unravel the phenomena of electric current transport by evaluating the different stages in which charge transfer could occur. The results revealed how the positive charge migrates from the organic layer to the gold surface. Specifically, the transport mechanism seems to take place mainly along a single ligand and driven with the help of the electrostatic interactions of the surrounding molecules. Aside, several self-assembled monolayers with different linkers have been analyzed to understand how the nature of that moiety can tune the redox properties and the efficiency of the transport. We have found that the conjugation between the guanine and the linker, at the same time conjugated to the gold surface, gives rise to a more efficient transport. In conclusion, the established computational protocol sheds light on the mechanism behind charge transport in electrochemical DNA-based biosensor nanodevices.

Toward a numerically efficient description of bulk-solvated anionic states

We investigate the vertical electron attachment energy (VAE) of 1-methyl-4-nitroimidazole, a model radiosensitizer, employing quantum mechanics/molecular mechanics (QM/MM) and QM/polarized continuum (QM/PCM) solvation models. We considered the solvent-excluded surface (QM/PCM-SES) and Van der Waals (QM/PCM-VDW) cavities within the PCM framework, the electrostatic embedding QM/MM (EE-QM/MM) model, and the self-consistent sequential QM/MM polarizable electrostatic embedding (scPEE-S-QM/MM) model. Due to slow VAE convergence concerning the number of QM solvent molecules, full QM calculations prove inefficient. Ensemble averages in these calculations do not align with VAEs computed for the representative solute-solvent configuration. QM/MM and QM/PCM calculations show agreement with each other for sufficiently large QM regions, although the QM/PCM-VDW model exhibits artifacts linked to the cavity. QM/MM models demonstrate good agreement between ensemble averages and VAEs calculated with the representative configuration. Notably, the VAE computed with the scPEE-S-QM/MM model achieves faster convergence concerning the number of QM water molecules compared to the EE-QM/MM model, attributed to enhanced efficiency from MM charge polarization in the scPEE-S-QM/MM approach. This emphasizes the importance of QM/classical models with accurate solute-solvent and solvent-solvent mutual polarization for obtaining converged VAEs at a reasonable computational cost. The full-QM approach is very inefficient, while the microsolvation model is inaccurate. Computational savings in QM/MM models result from electrostatic embedding and the representative configuration, with the scPEE-S-QM/MM approach emerging as an efficient tool for describing bulk-solvated anions within the QM/MM framework. Its potential extends to improving transient anion state descriptions in biomolecules and radiosensitizers, especially given the frequent employment of microsolvation models.

Modeling One-Electron Oxidation Potentials and Hole Delocalization in Double-Stranded DNA by Multilayer and Dynamic Approaches

The number of innovative applications for DNA nowadays is growing quickly. Its use as a nanowire or electrochemical biosensor leads to the need for a deep understanding of the charge-transfer process along the strand, as well as its redox properties. These features are computationally simulated and analyzed in detail throughout this work by combining molecular dynamics, multilayer schemes, and the Marcus theory. One-electron oxidation potential and hole delocalization have been analyzed for six DNA double strands that cover all possible binary combinations of nucleotides. The results have revealed that the one-electron oxidation potential decreases with respect to the single-stranded DNA, giving evidence that the greater rigidity of a double helix induces an increase in the capacity of storing the positive charge generated upon oxidation. In addition, the hole is mainly stored in nucleobases with large reducer character, i.e., purines, especially when those are arranged in a stacked configuration in the same strand. From the computational point of view, the sampling needed to describe biological systems implies a significant computational cost. Here, we show that a small number of representative conformations generated by clustering analysis provides accurate results when compared with those obtained from sampling, reducing considerably the computational cost.

One-Electron Oxidation Potentials and Hole Delocalization in Heterogeneous Single-Stranded DNA

The study of DNA processes is essential to understand not only its intrinsic biological functions but also its role in many innovative applications. The use of DNA as a nanowire or electrochemical biosensor leads to the need for a deep investigation of the charge transfer process along the strand as well as of the redox properties. In this contribution, the one-electron oxidation potential and the charge delocalization of the hole formed after oxidation are computationally investigated for different heterogeneous single-stranded DNA strands. We have established a two-step protocol: (i) molecular dynamics simulations in the frame of quantum mechanics/molecular mechanics (QM/MM) were performed to sample the conformational space; (ii) energetic properties were then obtained within a QM1/QM2/continuum approach in combination with the Marcus theory over an ensemble of selected geometries. The results reveal that the one-electron oxidation potential in the heterogeneous strands can be seen as a linear combination of that property within the homogeneous strands. In addition, the hole delocalization between different nucleobases is, in general, small, supporting the conclusion of a hopping mechanism for charge transport along the strands. However, charge delocalization becomes more important, and so does the tunneling mechanism contribution, when the reducing power of the nucleobases forming the strand is similar. Moreover, charge delocalization is slightly enhanced when there is a correlation between pairs of some of the interbase coordinates of the strand: twist/shift, twist/slide, shift/slide, and rise/tilt. However, the internal structure of the strand is not the predominant factor for hole delocalization but the specific sequence of nucleotides that compose the strand.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Intramolecular and intermolecular hole delocalization rules the reducer character of isolated nucleobases and homogeneous single-stranded DNA

The use of DNA strands as nanowires or electrochemical biosensors requires a deep understanding of charge transfer processes along the strand, as well as of the redox properties. These properties are computationally assessed in detail throughout this study. By applying molecular dynamics and hybrid QM/continuum and QM/QM/continuum schemes, the vertical ionization energies, adiabatic ionization energies, vertical attachment energies, one-electron oxidation potentials, and delocalization of the hole generated upon oxidation have been determined for nucleobases in their free form and as part of a pure single-stranded DNA. We show that the reducer ability of the isolated nucleobases is explained by the intramolecular delocalization of the positively charged hole, while the enhancement of the reducer character when going from aqueous solution to the strand correlates very well with the intermolecular hole delocalization. Our simulations suggest that the redox properties of DNA strands can be tuned by playing with the balance between intramolecular and intermolecular charge delocalization.

An Efficient Multilayer Approach to Model DNA-Based Nanobiosensors

In this work, we present a full computational protocol to successfully obtain the one-electron reduction potential of nanobiosensors based on a self-assembled monolayer of DNA nucleobases linked to a gold substrate. The model is able to account for conformational sampling and environmental effects at a quantum mechanical (QM) level efficiently, by combining molecular mechanics (MM) molecular dynamics and multilayer QM/MM/continuum calculations within the framework of Marcus theory. The theoretical model shows that a guanine-based biosensor is more prone to be oxidized than the isolated nucleobase in water due to the electrostatic interactions between the assembled guanine molecules. In addition, the redox properties of the biosensor can be tuned by modifying the nature of the linker that anchor the nucleobases to the metal support.

Ultrafast Intersystem Crossing Dynamics of 6-Selenoguanine in Water

Rationalizing the photochemistry of nucleobases where an oxygen is replaced by a heavier atom is essential for applications that exploit near-unity triplet quantum yields. Herein, we report on the ultrafast excited-state deactivation mechanism of 6-selenoguanine (6SeGua) in water by combining nonadiabatic trajectory surface-hopping dynamics with an electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. We find that the predominant relaxation mechanism after irradiation starts on the bright singlet S₂ state that converts internally to the dark S₁ state, from which the population is transferred to the triplet T₂ state via intersystem crossing and finally to the lowest T₁ state. This S₂ → S₁ → T₂ → T₁ deactivation pathway is similar to that observed for the lighter 6-thioguanine (6tGua) analogue, but counterintuitively, the T₁ lifetime of the heavier 6SeGua is shorter than that of 6tGua. This fact is explained by the smaller activation barrier to reach the T₁/S₀ crossing point and the larger spin-orbit couplings of 6SeGua compared to 6tGua. From the dynamical simulations, we also calculate transient absorption spectra (TAS), which provide two time constants (τ₁ = 131 fs and τ₂ = 191 fs) that are in excellent agreement with the experimentally reported value (τ_exp = 130 ± 50 fs) (Farrel et al. J. Am. Chem. Soc. 2018, 140, 11214). Intersystem crossing itself is calculated to occur with a time scale of 452 ± 38 fs, highlighting that the TAS is the result of a complex average of signals coming from different nonradiative processes and not intersystem crossing alone.

Low-Energy Electron Damage to Plasmid DNA in Thin Films: Dependence on Substrates, Surface Density, Charging, Environment, and Uniformity

The interaction of low-energy electrons (LEEs) with DNA plays a significant role in the mechanisms leading to biological damage induced by ionizing radiation, particularly in radiotherapy, and its sensitization by chemotherapeutic drugs and nanoparticles. Plasmids constitute the form of DNA found in mitochondria and appear as a suitable model of genomic DNA. In a search for the best LEE targets, damage was induced to plasmids, in thin films in vacuum, by 6, 10, and 100 eV electrons under single collision conditions. The yields of single- and double-strand breaks, other cluster damage, isolated base lesions, and crosslinks were measured by electrophoresis and enzyme treatment. The films were deposited on oriented graphite or polycrystalline tantalum, with or without DNA autoassembly via diaminopropane (Dap) intercalation. Yields were correlated with the influence of vacuum, film uniformity, surface density, substrates, and the DNA environment. Aided by surface potential measurements and scanning electron microscopy and atomic force microscopy images, the lyophilized Dap-DNA films were found to be the most practical high-quality targets. These studies pave the way to the fabrication of LEE target-films composed of plasmids intercalated with biomolecules that could mimic the cellular environment; for example, as a first step, by replacing Dap with an amino acid.

Computation of Oxidation Potentials of Solvated Nucleobases by Static and Dynamic Multilayer Approaches

The determination of the redox properties of nucleobases is of paramount importance to get insight into the charge-transfer processes in which they are involved, such as those occurring in DNA-inspired biosensors. Although many theoretical and experimental studies have been conducted, the value of the one-electron oxidation potentials of nucleobases is not well-defined. Moreover, the most appropriate theoretical protocol to model the redox properties has not been established yet. In this work, we have implemented and evaluated different static and dynamic approaches to compute the one-electron oxidation potentials of solvated nucleobases. In the static framework, two thermodynamic cycles have been tested to assess their accuracy against the direct determination of oxidation potentials from the adiabatic ionization energies. Then, the introduction of vibrational sampling, the effect of implicit and explicit solvation models, and the application of the Marcus theory have been analyzed through dynamic methods. The results revealed that the static direct determination provides more accurate results than thermodynamic cycles. Moreover, the effect of sampling has not shown to be relevant, and the results are improved within the dynamic framework when the Marcus theory is applied, especially in explicit solvent, with respect to the direct approach. Finally, the presence of different tautomers in water does not affect significantly the one-electron oxidation potentials.

Damage Induced to DNA and Its Constituents by 0-3 eV UV Photoelectrons†

The complex physical and chemical interactions between DNA and 0-3 eV electrons released by UV photoionization can lead to the formation of various lesions such as base modifications and cleavage, crosslinks and single strand breaks. Furthermore, in the presence of platinum chemotherapeutic agents, these electrons can cause clustered lesions, including double strand breaks. We explain the mechanisms responsible for these damages via the production 0-3 eV electrons by UVC radiation, and by UV photons of any wavelengths, when they are produced by photoemission from nanoparticles lying within about 10 nm from DNA. We review experimental evidence showing that a single 0-3 eV electron can produce these damages. The foreseen benefits UV-irradiation of nanoparticles targeted to the cell nucleus are mentioned in the context of cancer therapy, as well as the potential hazards to human health when they are present in cells.

Low energy (6-18 eV) electron scattering from condensed thymidine (dT) III: absolute electronic excitation cross sections

Absolute cross sections (CSs) for electronic excitation by low-energy electron (LEE) scattering, from condensed thymidine (dT) in the 6-18 eV incident energy range, were measured by high-resolution electron energy loss spectroscopy (HREELS). Various electron energy loss (EEL) spectra were acquired using 1 ML of dT condensed on a multilayer film of Ar held at about 20 K under ultra-high vacuum (∼1 × 10-11 Torr). dT is one of the most complex DNA constituents to be studied by HREELS and these spectra provide the first LEE energy-loss data for electronic excitation of a nucleoside. CSs for transitions to the states 13A', 13A'', 23A', 21A', 33A', 23A'', 43A', 33A'', 53A' and 51A' of dT were extracted from the EEL spectra. These states correlate to those previously measured for the thymine moiety. Two broad resonances are observed in the energy dependence of the CSs at around 8 and 10 eV; these energies are close to those found in earlier gas- and solid-phase studies on the interaction of LEEs with dT, thymine and related molecules. A quantitative comparison between the electronic CSs of dT and those of thymine and tetrahydrofuran indicates that no variation is induced in the electronic CSs of thymine upon chemically binding to a deoxyribose group.

Change of Initial Yield of a Hydrated Electron with Uridine Monophosphate Concentration Is Related to the Excitation Photon Energy in Transient Absorption Spectroscopy

The initial yield of a hydrated electron (e_aq^-) in a solution under laser pulse irradiation was investigated by pump-probe transient absorption spectroscopy. The initial quantum yield of e_aq^- varies with the concentration of uridine monophosphate (UMP). The variation of the concentration of e_aq^- is often used to study the prehydrated electron (e_pre^-) and e_aq^- attachment to UMP. The results of 320 and 260 nm excitations were compared. It was found that with the increase of UMP concentration, the initial yield of e_aq^- increases at 320 nm excitation, but decreases at 260 nm excitation. The further analysis indicates that some of the e_pre^- attachments to UMP before solvation at 260 nm excitation result in the decrease of the e_aq^- yield. In addition, the absorption of UMP to 260 nm also causes the decrease of the e_aq^- yield. After the excitation at 320 nm, the phosphate group of UMP can release electrons more easily than that of water molecules by two-photon absorption, and therefore the e_aq^- yield increases. With the increase of UMP concentration, the decay rate of e_aq^- increases because e_aq^- is captured by UMP. The change of excitation photon does not affect the reaction rate of e_aq^- attachment to UMP. The longer lifetime of e_aq^- obtained at 260 nm excitation than 320 nm excitation is induced by the larger e_aq^- escape probability at 260 nm excitation. Our results show that the femtosecond pulse pump-probe transient absorption spectroscopy method should be cautiously used because of its complexity in studying the e_pre^- attachment to nucleotides in an aqueous solution.

Electron and hole interactions with P, Z, and P:Z and the formation of mutagenic products by proton transfer reactions

Z would act as an electron acceptor and P would capture a hole in the unnatural DNA. The latter process would produce mutagenic products via a proton transfer reaction.

Notable effect of water on excess electron attachment to aqueous DNA deoxyribonucleosides

As excess electrons are vertically attached to aqueous deoxyribonucleosides, ∼50% of excess electrons would be delocalized over the water molecules.