Fatemeh Fateminasab, Mohammad Aarabi, Aurelien de La Lande, Reza Omidyan. J. Mol. Liquids. 2021, 344, 117961. doi.org/10.1016/j.molliq.2021.117961.
In the present study, we have employed a high-level density functional theory (DFT) model to investigate both implicit and explicit effects of solvation (i.e., the bulk solvation effect as well as microhydration) in addition to the effect of distal Histidine on the binding of CO to Ferrous and Ferric-Heme. It has been predicted that the distal N-methylimidazole (MI) in addition to microsolvation, as a simple mode of environment, leads to a significant stabilization on the binding of CO to Heme-FeII (by ∼ 4 kcal mol−1). This is while no clear alteration has been predicted for the implicit solvation model. For the Ferric heme, and in contrast to the Ferrous system, strong destabilization from the environment has been predicted for binding CO to the Ferric heme analog.
J. D Samaniego-Rojas, L. I. Hernández-Segura, L. López-Sosa, R. I Delgado-Venegas, B. Gomez, J.-C. Lambry, A. de la Lande, T. Mineva, J. Alejandre, B. A Zúñiga-Gutiérrez, R. Flores-Moreno, P. Calaminici, G. Geudtner, A. M Köster.Andreas in Multiscale Dynamics Simulations: Nano and Nano-bio Systems in Complex Environments.Edited by D. Wei and D. R. Salahub. 2021. Royal Chemical Society. doi.org/10.1039/9781839164668-00001
This chapter describes the theoretical background of the quantum mechanical/molecular mechanical (QM/MM) implementation in deMon2k within the framework of auxiliary density functional theory (ADFT). It aims to give the reader an overview of the current state of the art of this QM/MM implementation and perspectives for its future development. To this end, we first derive the ADFT working equations for the QM and QM/MM energy and gradient expressions. Based on the joint QM/MM gradient expression, we present algorithms for QM/MM structure optimizations, transition-state searches and molecular dynamics simulations. The use of auxiliary density perturbation theory (ADPT) in the framework of QM/MM is discussed using illustrative implementations including analytic second-order ADFT energy derivatives, nuclear magnetic resonance chemical shift calculations and excited state calculations using time-dependent ADFT. The chapter closes with the description of a transformation program used to generate deMon2k QM/MM inputs.
A. de la Lande, S. Denisov, M. Mostafavi. Phys. Chem. Chem. Phys., 2021, 23, 21148-21162. doi.org/10.1039/D0CP06482C. Link to HAL
The early mechanisms by which ionizing rays damage biological structures by so-called direct effects are largely elusive. In a recent picosecond pulse radiolysis study of concentrated uridine monophosphate solutions [J. Ma, S. A. Denisov, J.-L. Marignier, P. Pernot, A. Adhikary, S. Seki and M. Mostafavi, J. Phys. Chem. Lett., 2018, 9, 5105], unexpected results were found regarding the oxidation of the nucleobase. The signature of the oxidized nucleobase could not be detected 5 ps after the electron pulse, but only the oxidized phosphate, raising intriguing questions about the identity of charge-transfer mechanisms that could explain the absence of U+. We address here this question by means of advanced first-principles atomistic simulations of solvated uridine monophosphate, combining Density Functional Theory (DFT) with polarizable embedding schemes. We contrast three very distinct mechanisms of charge transfer covering the atto-, femto- and pico-second timescales. We first investigate the ionization mechanism and subsequent hole/charge migrations on a timescale of attoseconds to a few femtoseconds under the frozen nuclei approximation. We then consider a nuclear-driven phosphate-to-oxidized-nucleobase electron transfer, showing that it is an uncompetitive reaction channel on the sub-picosecond timescale, despite its high exothermicity and significant electronic coupling. Finally, we show that non-adiabatic charge transfer is enabled by femtosecond nuclear relaxation after ionization. We show that electronic decoherence and the electronic coupling strength are the key parameters that determine the hopping probabilities. Our results provide important insight into the interplay between electronics and nuclear motions in the early stages of the multiscale responses of biological matter subjected to ionizing radiation.
This review summarizes the state-of-the-art knowledge of heme ligation in the gas phase. The unique aspect of the gas phase approach is to allow a step-by-step ligation of heme and thus enables the analysis of the properties of -four, -five and -six coordinate hemes in vacuo, under conditions directly comparable with quantum calculations. This approach also allows the characterization of situations uncommon in Nature, completing the coordination spectrum of hemes: four coordinate heme and protonated heme, an intermediate between ferrous and ferric heme. Therefore, a complete set of systems is described for the ferrous and ferric cases and there is no discontinuity between the two oxidation states of iron, so that the same mechanisms are at work, donation and back donation of different strengths depending upon the ligand. The similarity of ligation properties in ferrous and ferric hemes is consistent with calculations of the electron density at the Fe atom level, rather independent of the formal oxidation state in contrast with the porphyrin cycle. Hemes spin states have been reviewed, for they identify the electronic distribution of the metal. In ligated ferrous and ferric hemes, we find that binding energy measurements combined with spectroscopy describe their properties most effectively.
We recently adapted the Auxiliary DFT framework as implemented in deMon2k to the simulation of time-dependent problems via the Runge and Gross equations. Our implementation of the so-called Real-Time-Time-Dependent ADFT (RT-TD-ADFT) fully benefits from the algorithms available in deMon2k to carry out variational density fitting, notably the MINRES algorithm recently proposed for self-consistent-field calculations. We test here MINRES for the first time in the context of RT-TD-ADFT. We report extensive benchmarks calculations to assess the reliability of the ADFT framework. These encompass the construction of absorption spectra in the gas phase and in solvent, the calculation of electronic stopping power curves, the irradiation of zeolites by swift ions and the investigation of charge migrations with attosecond time resolution. All our results are very encouraging. We show that even small auxiliary basis sets are sufficient to obtain results almost undisguisable from those obtained with large and flexible auxiliary bases. Overall, we establish the reliability of RT-TD-ADFT to simulate electronics dynamics in large or very large molecular systems.
n this article, we investigate the mechanisms of DNA ionization upon irradiation by 0.5 meV alpha particles. We focus on the sugar-phosphate group and its hydration shell. In radiation chemistry, the term quasi-direct effect refers the physical and chemical responses taking place after irradiation of solvent molecules pertaining to the solvation shells of solutes. The molecular mechanisms accounting for the quasi-direct effect are actually largely elusive, especially for those prevailing in the early timescales (< 10–12 s). We report Real-Time Time-Dependent Auxiliary Density Functional Theory simulations carried out within the framework of hybrid QM/MM scheme (Quantum Mechanics/Molecular Mechanics) with polarizable and non-polarizable embedding. Ten water molecules from the solvation shell of DNA backbone are independently irradiated. We find that during the first femtoseconds after irradiation, the holes formed on the irradiated water remain at their sites of formation. Electrostatic induction within the environment does not significantly impact charge migrations. We address the hypothesis that charge migration driven by electron correlation is responsible for an ultrafast H2O+ to DNA charge transfers, which would account for a quasi-direct effect. We find that pure charge migration at fixed nuclear positions is not responsible for the quasi-direct effect when considering sugar-phosphate solvation shells.
X. Wu, J. Hénin, L. Baciou, M. Baaden, F. Cailliez, A. de La Lande. Frontiers in Chemistry, 2021, 9. doi.org/10.3389/fchem.2021.650651.
NOX5 is a member of the NADPH oxidase family which is dedicated to the production of reactive oxygen species. The molecular mechanisms governing transmembrane electron transfer (ET) that permits to shuttle electrons over the biological membrane have remained elusive for a long time. Using computer simulations, we report conformational dynamics of NOX5 embedded within a realistic membrane environment. We assess the stability of the protein within the membrane and monitor the existence of cavities that could accommodate dioxygen molecules. We investigate the heme-to-heme electron transfer. We find a reaction free energy of a few tenths of eV (ca. −0.3 eV) and a reorganization free energy of around 1.1 eV (0.8 eV after including electrostatic induction corrections). The former indicates thermodynamically favorable ET, while the latter falls in the expected values for transmembrane inter-heme ET. We estimate the electronic coupling to fall in the range of the μeV. We identify electron tunneling pathways showing that not only the W378 residue is playing a central role, but also F348. Finally, we reveal the existence of two connected O2−binding pockets near the outer heme with fast exchange between the two sites on the nanosecond timescale. We show that when the terminal heme is reduced, O2 binds closer to it, affording a more efficient tunneling pathway than when the terminal heme is oxidized, thereby providing an efficient mechanism to catalyze superoxide production in the final step. Overall, our study reveals some key molecular mechanisms permitting reactive oxygen species production by NOX5 and paves the road for further investigation of ET processes in the wide family of NADPH oxidases by computer simulations.
Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.
Aurelio Alvarez-Ibarra, Angela Parise, Karim Hasnaoui, Aurélien de La Lande. Phys. Chem. Chem. Phys. 2020, 22, 7747-7758. doi.org/10.1039/D0CP00165A.
Selected by the editor as a PCCP hot paper.
The primary processes that occur following direct irradiation of bio-macromolecules by ionizing radiation determine the multiscale responses that lead to biomolecular lesions. The so-called physical stage loosely describes processes of energy deposition and molecular ionization/excitation but remains largely elusive. We propose a new approach based on first principles density functional theory to simulate energy deposition in large and heterogeneous biomolecules by high-energy-transfer particles. Unlike traditional Monte Carlo approaches, our methodology does not rely on pre-parametrized sets of cross-sections, but captures excitation, ionization and low energy electron emission at the heart of complex biostructures. It furthermore gives access to valuable insights on ultrafast charge and hole dynamics on the femtosecond time scale. With this new tool, we reveal the mechanisms of ionization by swift ions in microscopic DNA models and solvated DNA comprising almost 750 atoms treated at the DFT level of description. We reveal a so-called ebb-and-flow ionization mechanism in which polarization of the irradiated moieties appears as a key feature. We also investigate where secondary electrons produced by irradiation localize on chemical moieties composing DNA. We compare irradiation of solvated DNA by light (H+, and He2+) vs. heavier (C6+) ions, highlighting the much higher probability of double ionization with the latter. Our methodology constitutes a stepping stone towards a greater understanding of the chemical stage and more generally towards the multiscale modelling of radiation damage in biology using first principles.
Multicomponent Density Functional Theory (MDFT) is a promising methodology to incorporate nuclear quantum effects, such as zero-point energy or tunneling, or to simulate other types of particles such as muons or positrons using particle densities as basic quantities. As for standard electronic DFT, a still ongoing challenge is to achieve the most efficient implementations. We introduce a multicomponent DFT implementation within the framework of auxiliary DFT, focusing on molecular systems comprising electrons and quantum protons. We introduce a dual variational procedure to determine auxiliary electron and proton densities which leads to a succession of approximate energy expressions. Electronic and protonic fitted densities are employed in (i) electron-electron, proton-proton, and electron-proton classical Coulomb interactions and (ii) electron exchange-correlation, proton-proton exchange, and electron-proton correlation (EPC) potentials. If needed, exact exchange among electrons or among protons is computed by the variational fitting of the corresponding Fock potential. The implementation is carried out in deMon2k. We test various electron proton correlation functionals on proton affinities. We find that auxiliary densities can be safely used in electron-electron, proton-proton, and electron-proton classical Coulomb interactions, as well as in EPC, albeit with some precautions related to the choice of the electronic auxiliary basis set that must be flexible enough. Computational tests reported indicate that introduction of density fitting in MDFT is clearly advantageous in terms of computational effort with good scaling properties with respect to the number of electrons and protons treated at the DFT level.
