Dynamics of the excited-state hydrogen transfer in a (dG)·(dC) homopolymer: intrinsic photostability of DNA

Electronic spectroscopy of gemcitabine and derivatives for possible dual-action photodynamic therapy applications

In this paper, we explore the molecular basis of combining photodynamic therapy (PDT), a light-triggered targeted anticancer therapy, with the traditional chemotherapeutic properties of the well-known cytotoxic agent gemcitabine. A photosensitizer prerequisite is significant absorption of biocompatible light in the visible/near IR range, ideally between 600 and 1000 nm. We use highly accurate multiconfigurational CASSCF/MS-CASPT2/MM and TD-DFT methodologies to determine the absorption properties of a series of gemcitabine derivatives with the goal of red-shifting the UV absorption band toward the visible region and facilitating triplet state population. The choice of the substitutions and, thus, the rational design is based on important biochemical criteria and on derivatives whose synthesis is reported in the literature. The modifications tackled in this paper consist of: (i) substitution of the oxygen atom at O2 position with heavier atoms (O → S and O → Se) to red shift the absorption band and increase the spin-orbit coupling, (ii) addition of a lipophilic chain at the N7 position to enhance transport into cancer cells and slow down gemcitabine metabolism, and (iii) attachment of aromatic systems at C5 position to enhance red shift further. Results indicate that the combination of these three chemical modifications markedly shifts the absorption spectrum toward the 500 nm region and beyond and drastically increases spin-orbit coupling values, two key PDT requirements. The obtained theoretical predictions encourage biological studies to further develop this anticancer approach.

Solvent assisted excited-state deactivation pathways in isolated 2,7-diazaindole-S1-3 (S = Water and Ammonia) complexes

The role of solvent molecules in the deactivation of photo-excited 2,7-diazaindole (DAI) - (H₂O)_1-3 and DAI - (NH₃)_1-3 complexes were computationally investigated. An excited-state proton transfer (ESPT) path from the solvent to the DAI molecule was followed using the TD-DFT-D4 (B3LYP) level of theory. The computed potential energy profile of ESPT process has shown intersection between ππ* and nπ* states facilitated via relative stabilization of the nπ* state with decreasing N(7)-H_b bond length. The ESPT process, starting from the DAI-S_n (ππ*) state, crosses through a barrier ranging from 27 to 36 kJmol^-1 for water complexes and 26-30 kJmol^-1 for ammonia complexes. The energy of the excited state was rapidly decreased with a shorter N(7)-H_b bond length. Subsequently, a significant trend of finding a second intersection between the ground and the excited state was observed for all the complexes. The results firmly suggested a significant deactivation channel of excited azaindole derivatives. In the present system, two competing channels, ESPT and ESHT, were found to be energetically accessible. The energy barriers associated with the ESPT barriers for DAI-(H₂O)_1-3 complexes are similar to the ESHT barrier, depicting equal dominance of both processes. The increased basicity of the N(7) atom in the excited state resulted a facile ESPT process from the water to N(7) of the DAI molecule. However, DAI-(NH₃)_1-3 complexes show clear preference for ESHT over ESPT process owing to its higher gas-phase pK_a value making it a poor proton donor. The above systems can be used as a model to computationally and experimentally investigate the competing radiative and deactivation pathways of photo-excited solvated complexes of N-H-bearing bio-relevant molecules via proton and hydrogen transfer reactions.

Spatial and Temporal Resolution of the Oxygen-Independent Photoinduced DNA Interstrand Cross-Linking by a Nitroimidazole Derivative

DNA damage is ubiquitous in nature and is at the basis of emergent treatments such as photodynamic therapy, which is based on the activation of highly oxidative reactive oxygen species by photosensitizing O₂. However, hypoxia observed in solid tumors imposes the necessity to devise oxygen-independent modes of action able to induce DNA damage under a low oxygen concentration. The complexity of these DNA damage mechanisms in realistic environments grows exponentially when taking into account light absorption and subsequent excited-state population, photochemical and (photo)-redox reactions, the multiple species involved in different electronic states, noncovalent interactions, multiple reaction steps, and the large number of DNA reactive sites. This work tackles all the intricate reactivity of a photosensitizer based on a nitroimidazole derivative reacting toward DNA in solution under UV light exposition. This is performed through a combination of ground- and excited-state quantum chemistry, classical molecular dynamics, and hybrid QM/MM simulations to rationalize in detail the formation of DNA interstrand cross-links (ICLs) exerted by the noncanonical noncovalent photosensitizer. Unprecedented spatial and temporal resolution of these phenomena is achieved, revealing that the ICL is sequence-specific and that the fastest reactions take place at AT, GC, and GT steps involving either the opposite nucleobases or adjacent Watson-Crick base pairs. The N7 and O6 positions of guanine, the N7 and N3 sites of adenine, the N4 position of cytosine, and the O2 atom of thymine are deemed as the most nucleophile sites and are positively identified to participate in the ICL productions. This work provides a multiscale computational protocol to study DNA reactivity with noncovalent photosensitizers, and contributes to the understanding of therapies based on photoinduced DNA damage at molecular and electronic levels. In addition, we believe the depth understanding of these processes should assist the design of new photosensitizers considering their molecular size, electronic properties, and the observed regioselectivity toward nucleic acids.

The Behavior of Triplet Thymine in a Model B-DNA Strand. Energetics and Spin Density Localization Revealed by ab initio Molecular Dynamics Simulations†

Among the naturally occurring nucleobases, thymine presents the lowest triplet state, hence it represents a hotspot for energy transfer and photosensitization. In turn, the population of the triplet state may lead to thymine dimerization and hence to the production of toxic DNA lesions and has been the subject of intensive theoretical and experimental investigations. Relying on QM/MM molecular dynamics simulations, we have sought to situate the energy of the lowest triplet state of thymine embedded in a B-DNA environment. The energy gap varies between 305 and 329 kJ mol^-1 when a single thymine is treated at the quantum chemistry level, depending on its position in the model double-stranded 16-bp oligonucleotide. The energy of triplet state decreases up to 300 kJ mol^-1 , due to polarization effects, when we consider coupled stacked nucleobases up to the inclusion of four nucleobases. Our value lies in good agreement with the energy inferred experimentally by Miranda and coworkers (270 kJ mol^-1 ), and our theoretical exploration opens the door to investigations toward other more complex and biologically relevant environments, such as thymines embedded in nucleosome core particles. Our investigations also provide a reference for further studies using semi-empirical approaches such as density functional-based tight-binding, allowing to further rationalize sequence effects.

Fragment Diabatization Linear Vibronic Coupling Model for Quantum Dynamics of Multichromophoric Systems: Population of the Charge-Transfer State in the Photoexcited Guanine-Cytosine Pair

We introduce a method (FrD-LVC) based on a fragment diabatization (FrD) for the parametrization of a linear vibronic coupling (LVC) model suitable for studying the photophysics of multichromophore systems. In combination with effective quantum dynamics (QD) propagations with multilayer multiconfigurational time-dependent Hartree (ML-MCTDH), the FrD-LVC approach gives access to the study of the competition between intrachromophore decays, like those at conical intersections, and interchromophore processes, like exciton localization/delocalization and the involvement of charge-transfer (CT) states. We used FrD-LVC parametrized with time-dependent density functional theory (TD-DFT) calculations, adopting either CAM-B3LYP or ωB97X-D functionals, to study the ultrafast photoexcited QD of a guanine-cytosine (GC) hydrogen-bonded pair, within a Watson-Crick arrangement, considering up to 12 coupled diabatic electronic states and the effect of all of the 99 vibrational coordinates. The bright excited states localized on C and, especially, on G are predicted to be strongly coupled to the G → C CT state, which is efficiently and quickly populated after an excitation to any of the four lowest energy bright local excited states. Our QD simulations show that more than 80% of the excited population on G and ∼50% of that on C decay to this CT state in less than 50 fs. We investigate the role of vibronic effects in the population of the CT state and show that it depends mainly on its large reorganization energy so that it can occur even when it is significantly less stable than the bright states in the Franck-Condon region. At the same time, we document that the formation of the GC pair almost suppresses the involvement of dark nπ* excited states in the photoactivated dynamics.

Os(II) Oligothienyl Complexes as a Hypoxia-Active Photosensitizer Class for Photodynamic Therapy

Hypoxia presents a challenge to anticancer therapy, reducing the efficacy of many available treatments. Photodynamic therapy is particularly susceptible to hypoxia, given that its mechanism relies on oxygen. Herein, we introduce two new osmium-based polypyridyl photosensitizers that are active in hypoxia. The lead compounds emerged from a systematic study of two Os(II) polypyridyl families derived from 2,2'-bipyridine (bpy) or 4,4'-dimethyl-2,2'-bipyridine (dmb) as coligands combined with imidazo[4,5-f][1,10]phenanthroline ligands tethered to n = 0-4 thiophenes (IP-nT). The compounds were characterized and investigated for their spectroscopic and (photo)biological activities. The two hypoxia-active Os(II) photosensitizers had n = 4 thiophenes, with the bpy analogue 1-4T being the most potent. In normoxia, 1-4T had low nanomolar activity (half-maximal effective concentration (EC50) = 1-13 nM) with phototherapeutic indices (PI) ranging from 5500 to 55 000 with red and visible light, respectively. A sub-micromolar potency was maintained even in hypoxia (1% O2), with light EC50 and PI values of 732-812 nM and 68-76, respectively -currently among the largest PIs for hypoxic photoactivity. This high degree of activity coincided with a low-energy, long-lived (0.98-3.6 μs) mixed-character intraligand charge-transfer (3ILCT)/ligand-to-ligand charge-transfer (3LLCT) state only accessible in quaterthiophene complexes 1-4T and 2-4T. The coligand identity strongly influenced the photophysical and photobiological results in this study, whereby the bpy coligand led to longer lifetimes (3.6 μs) and more potent photo-cytotoxicity relative to those of dmb. The unactivated compounds were relatively nontoxic both in vitro and in vivo. The maximum tolerated dose for 1-4T and 2-4T in mice was greater than or equal to 200 mg kg-1, an excellent starting point for future in vivo validation.

Molecular Basis of SARS-CoV-2 Infection and Rational Design of Potential Antiviral Agents: Modeling and Simulation Approaches

The emergence in late 2019 of the coronavirus SARS-CoV-2 has resulted in the breakthrough of the COVID-19 pandemic that is presently affecting a growing number of countries. The development of the pandemic has also prompted an unprecedented effort of the scientific community to understand the molecular bases of the virus infection and to propose rational drug design strategies able to alleviate the serious COVID-19 morbidity. In this context, a strong synergy between the structural biophysics and molecular modeling and simulation communities has emerged, resolving at the atomistic level the crucial protein apparatus of the virus and revealing the dynamic aspects of key viral processes. In this Review, we focus on how in silico studies have contributed to the understanding of the SARS-CoV-2 infection mechanism and the proposal of novel and original agents to inhibit the viral key functioning. This Review deals with the SARS-CoV-2 spike protein, including the mode of action that this structural protein uses to entry human cells, as well as with nonstructural viral proteins, focusing the attention on the most studied proteases and also proposing alternative mechanisms involving some of its domains, such as the SARS unique domain. We demonstrate that molecular modeling and simulation represent an effective approach to gather information on key biological processes and thus guide rational molecular design strategies.

Photoinduced DNA Lesions in Dormant Bacteria: The Peculiar Route Leading to Spore Photoproducts Characterized by Multiscale Molecular Dynamics*

Some bacterial species enter a dormant state in the form of spores to resist to unfavorable external conditions. Spores are resistant to a wide series of stress agents, including UV radiation, and can last for tens to hundreds of years. Due to the suspension of biological functions, such as DNA repair, they accumulate DNA damage upon exposure to UV radiation. Differently from active organisms, the most common DNA photoproducts in spores are not cyclobutane pyrimidine dimers, but rather the so-called spore photoproducts. This noncanonical photochemistry results from the dry state of DNA and its binding to small, acid-soluble proteins that drastically modify the structure and photoreactivity of the nucleic acid. Herein, multiscale molecular dynamics simulations, including extended classical molecular dynamics and quantum mechanics/molecular mechanics based dynamics, are used to elucidate the coupling of electronic and structural factors that lead to this photochemical outcome. In particular, the well-described impact of the peculiar DNA environment found in spores on the favored formation of the spore photoproduct, given the small free energy barrier found for this path, is rationalized. Meanwhile, the specific organization of spore DNA precludes the photochemical path that leads to cyclobutane pyrimidine dimer formation.

Excimer Intermediates en Route to Long-Lived Charge-Transfer States in Single-Stranded Adenine DNA as Revealed by Nonadiabatic Dynamics

The ultrafast time evolution of a single-stranded adenine DNA is studied using a hybrid multiscale quantum mechanics/molecular mechanics (QM/MM) scheme coupled to nonadiabatic surface hopping dynamics. As a model, we use (dA)₂₀ where a stacked adenine tetramer is treated quantum chemically. The dynamical simulations combined with on-the-fly quantitative wave function analysis evidence the nature of the long-lived electronically excited states formed upon absorption of UV light. After a rapid decrease of the initially excited excitons, relaxation to monomer-like states and excimers occurs within 100 fs. The former monomeric states then relax into additional excimer states en route to forming stabilized charge-transfer states on a longer timescale of hundreds of femtoseconds. The different electronic-state characters is reflected on the spatial separation between the adenines: excimers and charge-transfer states show a much smaller spatial separation than the monomer-like states and the initially formed excitons.

On the Intrinsically Low Quantum Yields of Pyrimidine DNA Photodamages: Evaluating the Reactivity of the Corresponding Minimum Energy Crossing Points

The low quantum yield of photoformation of cyclobutane pyrimidine dimers and pyrimidine-pyrimidone (6-4) adducts in DNA bases is usually associated with the presence of more favorable nonreactive decay paths and with the unlikeliness of exciting the system in a favorable conformation. Here, we prove that the ability of the reactive conical intersection to bring the system either back to the absorbing conformation or to the photoproduct must be considered as a fundamental factor in the low quantum yields of the mentioned photodamage. In support of the proposed model, the one order of magnitude difference in the quantum yield of formation of the cyclobutane thymine dimer with respect to the thymine-thymine (6-4) adduct is rationalized here by comparing the reactive ability of the seam of intersections leading respectively to the cyclobutane thymine dimer and the oxetane precursor of the thymine-thymine (6-4) adduct at the CASPT2 level of theory.

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

First-principles characterization of the singlet excited state manifold in DNA/RNA nucleobases

TD-DFT characterization of the high-energy singlet excited state manifold of the canonical DNA/RNA nucleobasesin vacuumis assessed against RASPT2 reference computations for reliable simulations of linear and non-linear electronic spectra.

Enhanced reactivity of the pyrimidine peroxyl radical towards the C–H bond in duplex DNA – a theoretical study

The peroxyl radical exhibits a much stronger reactivity towards C1′–H1′ in duplex DNA with respect to single-stranded DNA.

UV-induced hydrogen transfer in DNA base pairs promoted by dark nπ* states

Formation of an excited-state complex enables ultrafast photorelaxation of dark nπ* states in GC and HC base pairs.

Experimental and theoretical studies on thymine photodimerization mediated by oxidatively generated DNA lesions and epigenetic intermediates

Combined spectroscopic and computational studies reveal that, in spite of their structural similarities, 5-formyluracil and 5-formylcytosine photosensitize cyclobutane thymine dimers through two different types of mechanisms.

Can microsolvation effects be estimated from vacuum computations? A case-study of alcohol decomposition at the H2O/Pt(111) interface

Activation and reaction energies of alcohol decomposition at Pt(111) are barely modified by a PCM, in contrast to adding a single water molecule, whose effect can be predicted based on vacuum computations.