Effect of the Thermal Fluctuations of the Photophysics of GC and CG DNA Steps: A Computational Dynamical Study

Here we refine and assess two computational procedures aimed to include the effect of thermal fluctuations on the electronic spectra and the ultrafast excited state dynamics of multichromophore systems, focusing on DNA duplexes. Our approach is based on a fragment diabatization procedure that, from a given Quantum Mechanical (QM) reference method, can provide the parameters (energy and coupling) of the reference diabatic states on the basis of the isolated fragments, either for a purely electronic excitonic Hamiltonian (FrDEx) or a linear vibronic coupling Hamiltonian (FrD-LVC). After having defined the most cost-effective procedure for DNA duplexes on two smaller fragments, FrDEx is used to simulate the absorption and Electronic Circular Dichroism (ECD) spectra of (GC)₅ sequences, including the coupling with the Charge Transfer (CT) states, on a number of structures extracted from classical Molecular Dynamics (MD) simulations. The computed spectra are close to the reference TD-DFT calculations and fully consistent with the experimental ones. We then couple MD simulations and FrD-LVC to simulate the interplay between local excitations and CT transitions, both intrastrand and interstrand, in GC and CG steps when included in a oligoGC or in oligoAT DNA sequence. We predict that for both sequences a substantial part of the photoexcited population on G and C decays, within 50-100 fs, to the corresponding intrastrand CT states. This transfer is more effective for GC steps that, on average, are more closely stacked than CG ones.

High-Energy Long-Lived Emitting Mixed Excitons in Homopolymeric Adenine-Thymine DNA Duplexes

The publication deals with polymeric pA●pT and oligomeric A₂₀●T₂₀ DNA duplexes whose fluorescence is studied by time-correlated single photon counting. It is shown that their emission on the nanosecond timescale is largely dominated by high-energy components peaking at a wavelength shorter than 305 nm. Because of their anisotropy (0.02) and their sensitivity to base stacking, modulated by the duplex size and the ionic strength of the solution, these components are attributed to mixed ππ*/charge transfer excitons. As high-energy long-lived excited states may be responsible for photochemical reactions, their identification via theoretical studies is an important challenge.

Why Do Most Aromatics Fail to Support Hole Hopping in the Cytochrome c Peroxidase-Cytochrome c Complex?

Electron transport through aromatic species (especially tryptophan and tyrosine) plays a central role in water splitting, redox signaling, oxidative damage protection, and bioenergetics. The cytochrome c peroxidase (CcP)-cytochrome c (Cc) complex (CcP:Cc) is used widely to study interprotein electron transfer (ET) mechanisms. Tryptophan 191 (Trp191) of CcP supports hole hopping charge recombination in the CcP:Cc complex. Experimental studies find that when Trp191 is substituted by tyrosine, phenylalanine, or redox-active aniline derivatives bound in the W191G cavity, enzymatic activity and charge recombination rates both decrease. Theoretical analysis of these CcP:Cc complexes finds that the ET kinetics depend strongly on the chemistry of the modified Trp site. The computed electronic couplings in the W191F and W191G species are orders of magnitude smaller than in the native protein, due largely to the absence of a hopping intermediate and the large tunneling distance. Small molecules bound in the W191G cavity are weakly coupled electronically to the Cc heme, and the structural disorder of the guest molecule in the binding pocket may contribute further to the lack of enzymatic activity. The couplings in W191Y are not substantially weakened compared to the native species, but the redox potential difference for tyrosine vs tryptophan oxidation accounts for the slower rate in the Tyr mutant. Thus, theoretical analysis explains why only the native Trp supports rapid hole hopping in the CcP:Cc complex. Favorable free energies and electronic couplings are essential for establishing an efficient hole hopping relay in this protein-protein complex.

Exciton transfer free energy from Car-Parrinello molecular dynamics

A computational approach is presented which allows the calculation of free energies profiles for exciton transfer processes within the framework of ab initio molecular dynamics (AIMD) simulations, sampling both the electronic and the nuclear degrees of freedom. To achieve this, restraining potentials are imposed on the centres of maximally localized Wannier orbitals. The resulting quantum-mechanical orbital forces are derived analytically and implemented in an AIMD program. In analogy to classical umbrella sampling techniques, these restraints are used to control an exciton transfer by incrementally moving the Wannier centres corresponding to the electron-hole pair along a suitable reaction coordinate. The new method is applied to study exciton transfer between two stacked penta(3-methylthiophene) molecules as a function of intermolecular distance. From the resulting free energy profiles, exciton transfer rates and diffusion constants are estimated, which prove to be in line with experimental results.

Computational modeling of photoexcitation in DNA single and double strands

The photoexcitation of DNA strands triggers extremely complex photoinduced processes, which cannot be understood solely on the basis of the behavior of the nucleobase building blocks. Decisive factors in DNA oligomers and polymers include collective electronic effects, excitonic coupling, hydrogen-bonding interactions, local steric hindrance, charge transfer, and environmental and solvent effects. This chapter surveys recent theoretical and computational efforts to model real-world excited-state DNA strands using a variety of established and emerging theoretical methods. One central issue is the role of localized vs delocalized excitations and the extent to which they determine the nature and the temporal evolution of the initial photoexcitation in DNA strands.

Photodynamic behavior of electronic coupling in a N-methylformamide dimer

The role of the bridging water molecules has been studied during the excited state photodynamics of a N-methylformamide dimer in complex with water molecules employing the complete active space self-consistent field (CASSCF) and CAS perturbation theory (CASPT2) methods.