Elise Dumont

Insights into the chemical meanings of the reaction electronic flux

Abstract The negative derivative of the chemical potential with respect to the reaction coordinate is called reaction electronic flux and has recently focused a wide interest to better understand chemical reactions at molecular level. After much consideration, it is now well accepted that positive REF values are associated with spontaneous processes, while negative REF ones translate unspontaneous phenomena. These characteristics of the REF are based on a thermodynamic analogy and have been shown right through computational results. In this paper, we develop two analytical expressions of the REF in both the canonical and the grand canonical ensembles. The connection between both equations is established. They are then analyzed, and some arguments are put forward to support the alleged characteristic of the REF and its ability to properly discriminate spontaneous from unspontaneous phenomena.

Interaction of Palmatine with DNA: An Environmentally Controlled Phototherapy Drug

Palmatine is one of the four main protoberberine alkaloids and is largely employed in pharmacy and medicine as a versatile drug with considerable biological activities. More recently, palmatine has been proposed as a promising DNA phototherapy drug, notably due to its ability to produce in situ singlet oxygen only when interacting with DNA. The fine mechanisms of palmatine-DNA interactions as well as its complicated photophysics are not yet fully understood. In this paper, we identify via molecular dynamic techniques two stable interaction modes between palmatine and B-DNA, namely insertion and minor groove binding, whose structural and electronic bases are analyzed and rationalized. These two competitive modes share the same UV-vis signature and estimated binding free energies, and thus they may indeed coexist. By using hybrid quantum mechanics/molecular mechanics protocols coupled to molecular dynamics, we analyze palmatine excited state properties in water solution and in interaction with DNA. The environmentally controlled production of singlet oxygen is thus rationalized in terms of the competition between local and charge-transfer excited states.

Resolving the Benzophenone DNA-Photosensitization Mechanism at QM/MM Level.

Benzophenone, the parent of the diarylketone family, is a versatile compound commonly used as a UV blocker. It may also trigger triplet-based DNA photosensitization. Therefore, benzophenone is involved in DNA photodamage induction. In the absence of experimentally resolved structure, the mechanism of DNA damage production remains elusive. Employing a hybrid quantum mechanics/molecular mechanics approach, here we address the spin transfer mechanism between this drug and proximal thymine, that is, the DNA nucleobase most prone to suffer triplet damages.

Understanding DNA under oxidative stress and sensitization: the role of molecular modeling

DNA is constantly exposed to damaging threats coming from oxidative stress, i.e., from the presence of free radicals and reactive oxygen species. Sensitization from exogenous and endogenous compounds that strongly enhance the frequency of light-induced lesions also plays an important role. The experimental determination of DNA lesions, though a difficult subject, is somehow well established and allows to elucidate even extremely rare DNA lesions. In parallel, molecular modeling has become fundamental to clearly understand the fine mechanisms related to DNA defects induction. Indeed, it offers an unprecedented possibility to get access to an atomistic or even electronic resolution. Ab initio molecular dynamics may also describe the time-evolution of the molecular system and its reactivity. Yet the modeling of DNA (photo-)reactions does necessitate elaborate multi-scale methodologies to tackle a damage induction reactivity that takes place in a complex environment. The double-stranded DNA environment is first characterized by a very high flexibility, but also a strongly inhomogeneous electrostatic embedding. Additionally, one aims at capturing more subtle effects, such as the sequence selectivity which is of critical important for DNA damage. The structure and dynamics of the DNA/sensitizers complexes, as well as the photo-induced electron- and energy-transfer phenomena taking place upon sensitization, should be carefully modeled. Finally the factors inducing different repair ratios for different lesions should also be rationalized. In this review we will critically analyze the different computational strategies used to model DNA lesions. A clear picture of the complex interplay between reactivity and structural factors will be sketched. The use of proper multi-scale modeling leads to the in-depth comprehension of DNA lesions mechanisms and also to the rational design of new chemo-therapeutic agents.

What singles out the G[8-5]C intrastrand DNA cross-link? Mechanistic and structural insights from quantum mechanics/molecular mechanics simulations.

Naturally occurring intrastrand oxidative cross-link lesions have proven to be a potent source of endogenous DNA damage. Among the variety of lesions that can be formed and have been identified, G[8-5]C damage (in which the C8 atom of a guanine is covalently bonded to the C5 atom of a nearby cytosine belonging to the same strand) occurs with a low incidence yet takes on special importance because of its high mutagenicity. Hybrid Car-Parrinello molecular dynamics simulations, rooted in density functional theory and coupled to molecular mechanics, have been performed to shed light on the cyclization process. The activation free energy of the reacting subsystem embedded in a solvated dodecamer is estimated to be ∼12.4 kcal/mol, which is ∼3 kcal/mol higher than the value for the prototypical G[8-5m]T lesion inferred employing the same theoretical framework [Garrec, J., Patel, C., Rothlisberger, U., and Dumont, E. (2012) J. Am. Chem. Soc.134, 2111-2119]. This study also situates the G[8-5m]mC lesion at an intermediate activation free energy (∼10.5 kcal/mol). The order of reactivity in DNA (T(•) > mC(•) > C(•)) is reversed compared to that in the reacting subsystems in the gas phase (C(•) > mC(•) > T(•)), stressing the crucial role of the solvated B-helix environment. The results of our simulations also characterize a more severe distortion for G[8-5]C than for methylene-bridged intrastrand cross-links.

Probing the reactivity of singlet oxygen with purines

The reaction of singlet molecular oxygen with purine DNA bases is investigated by computational means. We support the formation of a transient endoperoxide for guanine and by classical molecular dynamics simulations we demonstrate that the formation of this adduct does not affect the B-helicity. We thus identify the guanine endoperoxide as a key intermediate, confirming a low-temperature nuclear magnetic resonance proof of its existence, and we delineate its degradation pathway, tracing back the preferential formation of 8-oxoguanine versus spiro-derivates in B-DNA. Finally, the latter oxidized 8-oxodGuo product exhibits an almost barrierless reaction profile, and hence is found, coherently with experience, to be much more reactive than guanine itself. On the contrary, in agreement with experimental observations, singlet-oxygen reactivity onto adenine is kinetically blocked by a higher energy transition state.

Improved DFT Description of Intrastrand Cross-Link Formation by Inclusion of London Dispersion Corrections

The formation of covalent linkages between two vicinal nucleotides has been proved experimentally to constitute a particularly deleterious class of DNA lesions. These tandem lesions by essence present a competitive chemistry. The density functional theory with dispersion (DFT-D) method is shown to dramatically improve the theoretical description of the formation of a prototypical intrastrand cross-link, when compared to pure or hybrid GGA functionals which strongly deviate from the π-π self-stacking mode, as dinucleotides are artificially stabilized by the formation of unrealistic intramolecular hydrogen bonds (HBs). Inclusion of London dispersion correction restores a more realistic picture of the reactant structure and also of geometries and energies along the reaction profile. This paves the way toward a robust insilico screening of intrastrand cross-link DNA defects.

Probing deactivation pathways of DNA nucleobases by two-dimensional electronic spectroscopy: first principles simulations

The SOS//QM/MM [Rivalta et al., Int. J. Quant. Chem., 2014, 114, 85] method consists of an arsenal of computational tools allowing accurate simulation of one-dimensional (1D) and bi-dimensional (2D) electronic spectra of monomeric and dimeric systems with unprecedented details and accuracy. Prominent features like doubly excited local and excimer states, accessible in multi-photon processes, as well as charge-transfer states arise naturally through the fully quantum-mechanical description of the aggregates. In this contribution the SOS//QM/MM approach is extended to simulate time-resolved 2D spectra that can be used to characterize ultrafast excited state relaxation dynamics with atomistic details. We demonstrate how critical structures on the excited state potential energy surface, obtained through state-of-the-art quantum chemical computations, can be used as snapshots of the excited state relaxation dynamics to generate spectral fingerprints for different de-excitation channels. The approach is based on high-level multi-configurational wavefunction methods combined with non-linear response theory and incorporates the effects of the solvent/environment through hybrid quantum mechanics/molecular mechanics (QM/MM) techniques. Specifically, the protocol makes use of the second-order Perturbation Theory (CASPT2) on top of Complete Active Space Self Consistent Field (CASSCF) strategy to compute the high-lying excited states that can be accessed in different 2D experimental setups. As an example, the photophysics of the stacked adenine-adenine dimer in a double-stranded DNA is modeled through 2D near-ultraviolet (NUV) spectroscopy.

Correlation of bistranded clustered abasic DNA lesion processing with structural and dynamic DNA helix distortion

Clustered apurinic/apyrimidinic (AP; abasic) DNA lesions produced by ionizing radiation are by far more cytotoxic than isolated AP lesion entities. The structure and dynamics of a series of seven 23-bp oligonucleotides featuring simple bistranded clustered damage sites, comprising of two AP sites, zero, one, three or five bases 3′ or 5′ apart from each other, were investigated through 400 ns explicit solvent molecular dynamics simulations. They provide representative structures of synthetically engineered multiply damage sites-containing oligonucleotides whose repair was investigated experimentally (Nucl. Acids Res. 2004, 32:5609-5620; Nucl. Acids Res. 2002, 30: 2800–2808). The inspection of extrahelical positioning of the AP sites, bulge and non Watson–Crick hydrogen bonding corroborates the experimental measurements of repair efficiencies by bacterial or human AP endonucleases Nfo and APE1, respectively. This study provides unprecedented knowledge into the structure and dynamics of clustered abasic DNA lesions, notably rationalizing the non-symmetry with respect to 3′ to 5′ position. In addition, it provides strong mechanistic insights and basis for future studies on the effects of clustered DNA damage on the recognition and processing of these lesions by bacterial or human DNA repair enzymes specialized in the processing of such lesions.

Photophysics of Acetophenone Interacting with DNA: Why the Road to Photosensitization is Open

Deoxyribonucleic acid photosensitization, i.e. the photoinduced electron‐ or energy‐transfer of chromophores interacting with DNA, is a crucial phenomenon that triggers important DNA lesions such as pyrimidine dimerization, even upon absorption of relatively low‐energy radiation. Oxidative lesions may also be produced via the photoinduced production of reactive oxygen species. Aromatic ketones, and acetophenone in particular, are well known for their sensitization effects. In this contribution we model the structural and dynamical properties of the acetophenone/DNA aggregates as well as their spectroscopic and photophysical properties using high‐level hybrid quantum mechanics/molecular mechanics methods. We show that the key steps of the photochemistry of acetophenone in gas phase are conserved in the macromolecular environment and thus an ultrafast singlet–triplet conversion of acetophenone is expected prior to the transfer to DNA.