The physical stage of radiolysis of solvated DNA by high-energy-transfer particles: insights from new first principles simulations

Electron dynamics simulations based on density functional theory are carried out on nanometric molecular systems to decipher the primary processes following irradiation of bio-macromolecules by high energy transfer charged particles.

The surprisingly high ligation energy of CO to ruthenium porphyrins

A combined theoretical and experimental approach has been used to investigate the binding energy of a ruthenium metalloporphyrin ligated with CO, ruthenium tetraphenylporphyrin [RuII TPP], in the RuII oxidation degree. Measurements performed with VUV ionization using the DESIRS beamline at Synchrotron SOLEIL led to adiabatic ionization energies of [RuII TPP] and its complex with CO, [RuII TPP-CO], of 6.48 ± 0.03 eV and 6.60 ± 0.03 eV, respectively, while the ion dissociation threshold of [RuII TPP-CO]+ is measured to be 8.36 ± 0.03 eV using the ground-state neutral complex. These experimental data are used to derive the binding energies of the CO ligand in neutral and cationic complexes (1.88 ± 0.06 eV and 1.76 ± 0.06 eV, respectively) using a Born-Haber cycle. Density functional theory calculations, in very satisfactory agreement with the experimental results, help to get insights into the metal-ligand bond. Notably, the high ligation energies can be rationalized in terms of the ruthenium orbital structure, which is singular compared to that of the iron atom. Thus, beyond indications of a strengthening of the Ru-CO bond due to the decrease in the CO vibrational frequency in the complex as compared to the Fe-CO bond, high-level calculations are essential to accurately describe the metal ligand (CO) bond and show that the Ru-CO bond energy is strongly affected by the splitting of triplet and singlet spin states in uncomplexed [Ru TPP].

Theoretical estimation of redox potential of biological quinone cofactors

Redox potentials are essential to understand biological cofactor reactivity and to predict their behavior in biological media. Experimental determination of redox potential in biological system is often difficult due to complexity of biological media but computational approaches can be used to estimate them. Nevertheless, the quality of the computational methodology remains a key issue to validate the results. Instead of looking to the best absolute results, we present here the calibration of theoretical redox potential for quinone derivatives in water coupling QM + MM or QM/MM scheme. Our approach allows using low computational cost theoretical level, ideal for long simulations in biological systems, and determination of the uncertainties linked to the calculations. © 2017 Wiley Periodicals, Inc.

Electron transfer, decoherence, and protein dynamics: insights from atomistic simulations

Electron transfer in biological systems drives the processes of life. From cellular respiration to photosynthesis and enzymatic catalysis, electron transfers (ET) are chemical processes on which essential biological functions rely. Over the last 40 years, scientists have sought understanding of how these essential processes function in biology. One important breakthrough was the discovery that Marcus theory (MT) of electron transfer is applicable to biological systems. Chemists have experimentally collected both the reorganization energies (λ) and the driving forces (ΔG°), two parameters of Marcus theory, for a large variety of ET processes in proteins. At the same time, theoretical chemists have developed computational approaches that rely on molecular dynamics and quantum chemistry calculations to access numerical estimates of λ and ΔG°. Yet another crucial piece in determining the rate of an electron transfer is the electronic coupling between the initial and final electronic wave functions. This is an important prefactor in the nonadiabatic rate expression, since it reflects the probability that an electron tunnels from the electron donor to the acceptor through the intervening medium. The fact that a protein matrix supports electron tunneling much more efficiently than vacuum is now well documented, both experimentally and theoretically. Meanwhile, many chemists have provided examples of the rich physical chemistry that can be induced by protein dynamics. This Account describes our studies of the dynamical effects on electron tunneling. We present our analysis of two examples of natural biological systems through MD simulations and tunneling pathway analyses. Through these examples, we show that protein dynamics sustain efficient tunneling. Second, we introduce two time scales: τcoh and τFC. The former characterizes how fast the electronic coupling varies with nuclear vibrations (which cause dephasing). The latter reflects the time taken by the system to leave the crossing region. In the framework of open quantum systems, τFC is a short time approximation of the characteristic decoherence time of the electronic subsystem in interaction with its nuclear environment. The comparison of the respective values of τcoh and τFC allows us to probe the occurrence of non-Condon effects. We use ab initio MD simulations to analyze how decoherence appears in several biological cofactors. We conclude that we cannot account for its order of magnitude by considering only the atoms or bonds directly concerned with the transfer. Decoherence results from contributions from all atoms of the system appearing with a time delay that increases with the distance from the primarily concerned atoms or bonds. The delay and magnitude of the contributions depend on the chemical nature of the system. Finally, we present recent developments based on constrained DFT for efficient and accurate evaluations of the electronic coupling in ab initio MD simulations. These are promising methods to study the subtle fluctuations of the electronic coupling and the mechanisms of electronic decoherence in biological systems.

Enhanced semiempirical QM methods for biomolecular interactions

Recent successes and failures of the application of 'enhanced' semiempirical QM (SQM) methods are reviewed in the light of the benefits and backdraws of adding dispersion (D) and hydrogen-bond (H) correction terms. We find that the accuracy of SQM-DH methods for non-covalent interactions is very often reported to be comparable to dispersion-corrected density functional theory (DFT-D), while computation times are about three orders of magnitude lower. SQM-DH methods thus open up a possibility to simulate realistically large model systems for problems both in life and materials science with comparably high accuracy.