Solvent Reorganization Energies in A-DNA, B-DNA, and Rhodamine 6G−DNA Complexes from Molecular Dynamics Simulations with a Polarizable Force Field

Transient and Enduring Electronic Resonances Drive Coherent Long Distance Charge Transport in Molecular Wires

It is shown that the yields of oxidative damage observed in double-stranded DNA oligomers consisting of two guanines separated by adenine-thymine (A:T) _n bridges of various lengths are reliably accounted for by a multistep mechanism, in which transient and nontransient electronic resonances induce charge transport and solvent relaxation stabilizes the hole transfer products. The proposed multistep mechanism leads to results in excellent agreement with the observed yield ratios for both the short and the long distance regime; the almost distance independence of yield ratios for longer bridges ( n ≥ 3) is the consequence of the significant energy decrease of the electronic levels of the bridge, which, as the bridge length increases, become quasi-degenerate with those of the acceptor and donor groups (enduring resonance). These results provide significant guidelines for the design of novel DNA sequences to be employed in organic electronics.

How abasic sites impact hole transfer dynamics in GC-rich DNA sequences

Hole transfer dynamics through GC-rich DNA duplexes containing abasic sites is strongly modulated by the nature of the unpaired nucleobase.

Progress and challenges in simulating and understanding electron transfer in proteins

This Review presents an overview of the most common numerical simulation approaches for the investigation of electron transfer (ET) in proteins. We try to highlight the merits of the different approaches but also the current limitations and challenges. The article is organized into three sections. Section 2 deals with direct simulation algorithms of charge migration in proteins. Section 3 summarizes the methods for testing the applicability of the Marcus theory for ET in proteins and for evaluating key thermodynamic quantities entering the reaction rates (reorganization energies and driving force). Recent studies interrogating the validity of the theory due to the presence of non-ergodic effects or of non-linear responses are also described. Section 4 focuses on the tunneling aspects of electron transfer. How can the electronic coupling between charge transfer states be evaluated by quantum chemistry approaches and rationalized? What interesting physics regarding the impact of protein dynamics on tunneling can be addressed? We will illustrate the different sections with examples taken from the literature to show what types of system are currently manageable with current methodologies. We also take care to recall what has been learned on the biophysics of ET within proteins thanks to the advent of atomistic simulations.

Photoinduced Electron Transfer in DNA: Charge Shift Dynamics Between 8-Oxo-Guanine Anion and Adenine

Femtosecond time-resolved IR spectroscopy is used to investigate the excited-state dynamics of a dinucleotide containing an 8-oxoguanine anion at the 5'-end and neutral adenine at the 3'-end. UV excitation of the dinucleotide transfers an electron from deprotonated 8-oxoguanine to its π-stacked neighbor adenine in less than 1 ps, generating a neutral 8-oxoguanine radical and an adenine radical anion. These species are identified by the excellent agreement between the experimental and calculated IR difference spectra. The quantum efficiency of this ultrafast charge shift reaction approaches unity. Back electron transfer from the adenine radical anion to the 8-oxguanine neutral radical occurs in 9 ps, or approximately 6 times faster than between the adenine radical anion and the 8-oxoguanine radical cation (Zhang, Y. et al. Proc. Natl. Acad. Sci. U.S.A. 2014, 111, 11612-11617). The large asymmetry in forward and back electron transfer rates is fully rationalized by semiclassical nonadiabatic electron transfer theory. Forward electron transfer is ultrafast because the driving force is nearly equal to the reorganization energy, which is estimated to lie between 1 and 2 eV. Back electron transfer is highly exergonic and takes place much more slowly in the Marcus inverted region.

Electrostatics interactions in classical simulations

Electrostatic interactions are crucial for both the accuracy and performance of atomistic biomolecular simulations. In this chapter we review well-established methods and current developments aiming at efficiency and accuracy. Specifically, we review the classical Ewald summations, particle-particle particle-method particle-method Ewald algorithms, multigrid, fast multipole, and local methods. We also highlight some recent developments targeting more accurate, yet classical, representation of the molecular charge distribution.

Polarizable force fields

This chapter provides an overview of the most common methods for including an explicit description of electronic polarization in molecular mechanics force fields: the induced point dipole, shell, and fluctuating charge models. The importance of including polarization effects in biomolecular simulations is discussed, and some of the most important achievements in the development of polarizable biomolecular force fields to date are highlighted.

A comparative study of the binding of QSY 21 and Rhodamine 6G fluorescence probes to DNA: structure and dynamics

Molecular dynamics (MD) simulations and ab initio quantum chemical calculations were employed to investigate the structure, dynamics and interactions of the QSY 21 nonfluorescent quencher and the fluorescence dye Rhodamine 6G bound to a B-DNA decamer. For QSY 21, two binding motifs were observed. In the first motif, the central xanthene ring is stacked on one base of the adjacent cytosine-guanine DNA base pair, whereas one of the 2,3-dihydro-1-indolyl aromatic side rings is stacked on the other base. In the second motif, the QSY 21 stacking interaction with the DNA base pair is mediated only by one of the side rings. Several transitions between the motifs are observed during a MD simulation. The ab initio calculations show that none of these motifs is energetically preferred. Two binding motifs were found also for Rhodamine 6G, with the xanthene ring stacked predominantly either on the cytosine or on the guanine. These results suggest that the side rings of QSY 21 play a crucial role in its stacking on the DNA and indicate novel binding mode absent in the case of Rhodamine 6G, which lacks aromatic side rings.

A detailed binding free energy study of 2:1 ligand-DNA complex formation by experiment and simulation

In 2004, we used NMR to solve the structure of the minor groove binder thiazotropsin A bound in a 2:1 complex to the DNA duplex, d(CGACTAGTCG)2. In this current work, we have combined theory and experiment to confirm the binding thermodynamics of this system. Molecular dynamics simulations that use polarizable or non-polarizable force fields with single and separate trajectory approaches have been used to explore complexation at the molecular level. We have shown that the binding process invokes large conformational changes in both the receptor and ligand, which is reflected by large adaptation energies. This is compensated for by the net binding free energy, which is enthalpy driven and entropically opposed. Such a conformational change upon binding directly impacts on how the process must be simulated in order to yield accurate results. Our MM-PBSA binding calculations from snapshots obtained from MD simulations of the polarizable force field using separate trajectories yield an absolute binding free energy (-15.4 kcal mol(-1)) very close to that determined by isothermal titration calorimetry (-10.2 kcal mol(-1)). Analysis of the major energy components reveals that favorable non-bonded van der Waals and electrostatic interactions contribute predominantly to the enthalpy term, whilst the unfavorable entropy appears to be driven by stabilization of the complex and the associated loss of conformational freedom. Our results have led to a deeper understanding of the nature of side-by-side minor groove ligand binding, which has significant implications for structure-based ligand development.

Polarization effects in molecular mechanical force fields

The focus here is on incorporating electronic polarization into classical molecular mechanical force fields used for macromolecular simulations. First, we briefly examine currently used molecular mechanical force fields and the current status of intermolecular forces as viewed by quantum mechanical approaches. Next, we demonstrate how some components of quantum mechanical energy are effectively incorporated into classical molecular mechanical force fields. Finally, we assess the modeling methods of one such energy component-polarization energy-and present an overview of polarizable force fields and their current applications. Incorporating polarization effects into current force fields paves the way to developing potentially more accurate, though more complex, parameterizations that can be used for more realistic molecular simulations.