Emilie Cauet

Histidine-aromatic interactions in proteins and protein-ligand complexes: quantum chemical study of X-ray and model structures

His-aromatic complexes, with the His located above the aromatic plane, are stabilized by π-π, δ(+)-π and/or cation-π interactions according to whether the His is neutral or protonated and the partners are in stacked or T-shape conformations. Here we attempt to probe the relative strength of these interactions as a function of the geometry and protonation state, in gas phase, in water and protein-like environments (acetone, THF and CCl4), by means of quantum chemistry calculations performed up to second order of the Møller-Plesset pertubation theory. Two sets of conformations are considered for that purpose. The first set contains 89 interactions between His and Phe, Tyr, Trp, or Ade, observed in X-ray structures of proteins and protein-ligand complexes. The second set contains model structures obtained by moving an imidazolium/imidazole moiety above a benzene ring or an adenine moiety. We found that the protonated complexes are much more stable than the neutral ones in gas phase. This higher stability is due to the electrostatic contributions, the electron correlation contributions being equally important in the two forms. Thus, π-π and δ(+)-π interactions present essentially favorable electron correlation energy terms, whereas cation-π interactions feature in addition favorable electrostatic energies. The protonated complexes remain more stable than the neutral ones in protein-like environments, but the difference is drastically reduced. Furthermore, the T-shape conformation is undoubtedly more favorable than the stacked one in gas phase. This advantage decreases in the solvents, and the stacked conformation becomes even slightly more favorable in water. The frequent occurrence of His-aromatic interactions in catalytic sites, at protein-DNA or protein-ligand interfaces and in 3D domain swapping proteins emphasize their importance in biological processes.

Probing Electron-Induced Bond Cleavage at the Single-Molecule Level Using DNA Origami Templates

Low-energy electrons (LEEs) play an important role in nanolithography, atmospheric chemistry, and DNA radiation damage. Previously, the cleavage of specific chemical bonds triggered by LEEs has been demonstrated in a variety of small organic molecules such as halogenated benzenes and DNA nucleobases. Here we present a strategy that allows for the first time to visualize the electron-induced dissociation of single chemical bonds within complex, but well-defined self-assembled DNA nanostructures. We employ atomic force microscopy to image and quantify LEE-induced bond dissociations within specifically designed oligonucleotide targets that are attached to DNA origami templates. In this way, we use a highly selective approach to compare the efficiency of the electron-induced dissociation of a single disulfide bond with the more complex cleavage of the DNA backbone within a TT dinucleotide sequence. This novel technique enables the fast and parallel determination of DNA strand break yields with unprecedented control over the DNAs primary and secondary structure. Thus the detailed investigation of DNA radiation damage in its most natural environment, e.g., DNA nucleosomes constituting the chromatin, now becomes feasible.

Vertical ionization potentials of nucleobases in a fully solvated DNA environment.

Vertical ionization potentials (IPs) of nucleobases embedded in a fully solvated DNA fragment (12-mer B-DNA fragment + 22 sodium counterions + 5760 water molecules equilibrated to 298 K) have been calculated using a combined quantum mechanical molecular mechanics (QM/MM) approach. Calculations of the vertical IP of the anion Cl(-) are reported that support the accuracy of the application of a QM/MM method to this problem. It is shown that the pi nucleotide HOMO origin for the emitted electron is localized on the base by the hydration structure surrounding the DNA in a way similar to that recently observed for pyrimidine nucleotides in aqueous solutions (Slavicek, P.; et al. J. Am. Chem. Soc. 2009, 131, 6460). In a first step, a high level of theory, CCSD(T)/aug-cc-pVDZ, was used to calculate the vertical IP of each of the four single bases isolated in the QM region while the remaining DNA fragment, counterions, and water solvent molecules were included in the MM region. The calculated vertical IPs show a large positive shift of 3.2-3.3 eV compared to the corresponding gas-phase values. This shift is similar for all four DNA bases. The origin of the large increase in vertical IPs of nucleobases is found to be the long-range electrostatic interactions with the solvation structure outside the DNA helix. Thermal fluctuations in the fluid can result in IP changes of roughly 1 eV on a picosecond time scale. IPs of pi-stacked and H-bonded clusters of DNA bases were also calculated using the same QM/MM model but with a lower level of theory, B3LYP/6-31G(d=0.2). An IP shift of 4.02 eV relative to the gas phase is found for a four-base-pair B-DNA duplex configuration. The primary goal of this work was to estimate the influence of long-range solvation interactions on the ionization properties of DNA bases rather than provide highly precise IP evaluations. The QM/MM model presented in this work provides an attractive method to treat the difficult problem of incorporating a detailed long-range structural model of physiological conditions into investigations of the electronic processes in DNA.

Sequence dependence of electron-induced DNA strand breakage revealed by DNA nanoarrays

The electronic structure of DNA is determined by its nucleotide sequence, which is for instance exploited in molecular electronics. Here we demonstrate that also the DNA strand breakage induced by low-energy electrons (18 eV) depends on the nucleotide sequence. To determine the absolute cross sections for electron induced single strand breaks in specific 13 mer oligonucleotides we used atomic force microscopy analysis of DNA origami based DNA nanoarrays. We investigated the DNA sequences 5′-TT(XYX)3TT with X = A, G, C and Y = T, BrU 5-bromouracil and found absolute strand break cross sections between 2.66 · 10−14 cm2 and 7.06 · 10−14 cm2. The highest cross section was found for 5′-TT(ATA)3TT and 5′-TT(ABrUA)3TT, respectively. BrU is a radiosensitizer, which was discussed to be used in cancer radiation therapy. The replacement of T by BrU into the investigated DNA sequences leads to a slight increase of the absolute strand break cross sections resulting in sequence-dependent enhancement factors between 1.14 and 1.66. Nevertheless, the variation of strand break cross sections due to the specific nucleotide sequence is considerably higher. Thus, the present results suggest the development of targeted radiosensitizers for cancer radiation therapy.

Cisplatin cytotoxicity: A theoretical study of induced mutations

We investigate possible mutations in the genetic code induced by cisplatin with an approach combining molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Specifically, the impact of platination on the natural tautomeric equilibrium in guanine-cytosine (GC) base pairs is assessed to disclose the possible role played by non-canonical forms in anti-tumour activity. To obtain valuable predictions, the main interactions present in a real DNA environment, namely hydration and stacking, are simultaneously taken into account. According to our results, the Pt-DNA adduct promotes a single proton transfer reaction in GC in the DNA sequence AG[combining low line]G[combining low line]C. Such rare tautomers might play an important role in the cisplatin biological activity since they meet the stability requirements necessary to promote a permanent mutation.

Near-Quantitative Agreement of Model-Free DFT-MD Predictions with XAFS Observations of the Hydration Structure of Highly Charged Transition-Metal Ions.

First-principles dynamics simulations (DFT, PBE96, and PBE0) and electron scattering calculations (FEFF9) provide near-quantitative agreement with new and existing XAFS measurements for a series of transition-metal ions interacting with their hydration shells via complex mechanisms (high spin, covalency, charge transfer, etc.). This analysis does not require either the development of empirical interparticle interaction potentials or structural models of hydration. However, it provides consistent parameter-free analysis and improved agreement with the higher-R scattering region (first- and second-shell structure, symmetry, dynamic disorder, and multiple scattering) for this comprehensive series of ions. DFT+GGA MD methods provide a high level of agreement. However, improvements are observed when exact exchange is included. Higher accuracy in the pseudopotential description of the atomic potential, including core polarization and reducing core radii, was necessary for very detailed agreement. The first-principles nature of this approach supports its application to more complex systems.

Benchmarking DFT and TD-DFT Functionals for the Ground and Excited States of Hydrogen-Rich Peptide Radicals

We assess the pros and cons of a large panel of DFT exchange-correlation functionals for the prediction of the electronic structure of hydrogen-rich peptide radicals formed after electron attachment on a protonated peptide. Indeed, despite its importance in the understanding of the chemical changes associated with the reduction step, the question of the attachment site of an electron and, more generally, of the reduced species formed in the gas phase through electron-induced dissociation (ExD) processes in mass spectrometry is still a matter of debate. For hydrogen-rich peptide radicals in which several positive groups and low-lying π* orbitals can capture the incoming electron in ExD, inclusion of full Hartree-Fock exchange at long-range interelectronic distance is a prerequisite for an accurate description of the electronic states, thereby excluding several popular exchange-correlation functionals, e.g., B3LYP, M06-2X, or CAM-B3LYP. However, we show that this condition is not sufficient by comparing the results obtained with asymptotically correct range-separated hybrids (M11, LC-BLYP, LC-BPW91, ωB97, ωB97X, and ωB97X-D) and with reference CASSCF-MRCI and EOM-CCSD calculations. The attenuation parameter ω significantly tunes the spin density distribution and the excited states vertical energies. The investigated model structures, ranging from methylammonium to hexapeptide, allow us to obtain a description of the nature and energy of the electronic states, depending on (i) the presence of hydrogen bond(s) around the cationic site(s), (ii) the presence of π* molecular orbitals (MOs), and (iii) the selected DFT approach. It turns out that, in the present framework, LC-BLYP and ωB97 yields the most accurate results.

Radical cations of the nucleic bases and radiation damage to DNA: ab initio study

Abstract This review summarizes the contribution of high level quantum chemical calculations to the investigation of some elementary reactive processes related to the radiation damage to DNA. It is focused on the biomimetic species that govern these processes at the molecular level. These species are the DNA bases, isolated or embedded in base clusters. Their cations, formed by ionization in their ground and first excited electronic states, are at the center of the present work. We present a synthetic and critical overview of the computational methods used to predict accurate ionization potentials, to correctly describe the non-bonding interactions (stacking, H-bonding and cation- π ) stabilizing the studied biomimetic clusters, to characterize their excited states and to investigate the topology of the corresponding potential energy surfaces (minima, transition states, avoided crossings, conical intersections, reaction paths). All these aspects are illustrated by the recent literature and by our own research work, namely on the electron transfer occurring within a stacked dimer of guanines.

Accurate ground-state potential energy surfaces of the C2H2–Kr and C2H2–Xe van der Waals complexes

Accurate ab initio intermolecular potential energy surfaces (IPES) have been obtained for the first time for the ground electronic state of the C2H2–Kr and C2H2–Xe van der Waals complexes. Extensive tests, including complete basis set and all-electron scalar relativistic results, support their calculation at the CCSD(T) level of theory, using small-core relativistic pseudopotentials for the rare-gas atoms and aug-cc-pVQZ basis sets extended with a set of 3s3p2d1f1g mid-bond functions. All results are corrected for the basis set superposition error. The importance of the scalar relativistic and rare-gas outer-core (n–1)d correlation effects is investigated. The calculated IPES, adjusted to analytical functions, are characterized by global minima corresponding to skew T-shaped geometries, in which the Jacobi vector positioning the rare-gas atom with respect to the center of mass of the C2H2 moiety corresponds to distances of 4.064 and 4.229 Å, and angles of 65.22° and 68.67° for C2H2–Kr and C2H2–Xe, respectively. The interaction energy of both complexes is estimated to be −151.88 (1.817 kJ mol−1) and −182.76 cm−1 (2.186 kJ mol−1), respectively. The evolution of the topology of the IPES as a function of the rare-gas atom, from He to Xe, is also discussed.

Unique hole-trapping property of the human telomere sequence

The hypothetical protection of genes from oxidative damage provided by the G-rich telomeric overhangs located at the end of chromosomes, which consist, in humans, of single strands of TTAGGG sequence repeats, is investigated here. First principle Moller-Plesset perturbation theory calculations reveal that the TTAGGG human telomere sequence is particularly prone to oxidation and can act as a profound hole trap as deep as a sequence of five consecutive guanines. In addition, we show that the sequence dependence is very important and that modifications in the human telomeric sequence can induce crucial changes in the electronic structure of the sequence, with concomitant increase of the ionization energy. These theoretical results provide, for the first time, quantitative data indicating a high and unique efficiency of the human telomeric sequence as a trap in long-range hole migration which will aid in the design of subsequent experiments.Abstract The hypothetical protection of genes from oxidative damage provided by the G-rich telomeric overhangs located at the end of chromosomes, which consist, in humans, of single strands of TTAGGG sequence repeats, is investigated here. First principle Møller-Plesset perturbation theory calculations reveal that the TTAGGG human telomere sequence is particularly prone to oxidation and can act as a profound hole trap as deep as a sequence of five consecutive guanines. In addition, we show that the sequence dependence is very important and that modifications in the human telomeric sequence can induce crucial changes in the electronic structure of the sequence, with concomitant increase of the ionization energy. These theoretical results provide, for the first time, quantitative data indicating a high and unique efficiency of the human telomeric sequence as a trap in long-range hole migration which will aid in the design of subsequent experiments.