1
|
Migliore A, Messina A. Controlling the charge-transfer dynamics of two-level systems around avoided crossings. J Chem Phys 2024; 160:084112. [PMID: 38415830 DOI: 10.1063/5.0188749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 02/01/2024] [Indexed: 02/29/2024] Open
Abstract
Two-level quantum systems are fundamental physical models that continue to attract growing interest due to their crucial role as a building block of quantum technologies. The exact analytical solution of the dynamics of these systems is central to control theory and its applications, such as that to quantum computing. In this study, we reconsider the two-state charge transfer problem by extending and using a methodology developed to study (pseudo)spin systems in quantum electrodynamics contexts. This approach allows us to build a time evolution operator for the charge transfer system and to show new opportunities for the coherent control of the system dynamics, with a particular emphasis on the critical dynamic region around the transition state coordinate, where the avoided crossing of the energy levels occurs. We identify and propose possible experimental implementations of a class of rotations of the charge donor (or acceptor) that endow the electronic coupling matrix element with a time-dependent phase that can be employed to realize controllable coherent dynamics of the system across the avoided level crossing. The analogy of these rotations to reference frame rotations in generalized semiclassical Rabi models is discussed. We also show that the physical rotations in the charge-transfer systems can be performed so as to implement quantum gates relevant to quantum computing. From an exquisitely physical-mathematical viewpoint, our approach brings to light situations in which the time-dependent state of the system can be obtained without resorting to the special functions appearing in the Landau-Zener approach.
Collapse
Affiliation(s)
- Agostino Migliore
- Department of Chemical Sciences, University of Padova, Via Marzolo 1, 35131 Padova, Italy
| | - Antonino Messina
- Dipartimento di Matematica e Informatica, Università degli Studi di Palermo, Via Archirafi, 34, I-90123 Palermo, Italy
| |
Collapse
|
2
|
Ziogos OG, Kubas A, Futera Z, Xie W, Elstner M, Blumberger J. HAB79: A new molecular dataset for benchmarking DFT and DFTB electronic couplings against high-level ab initio calculations. J Chem Phys 2021; 155:234115. [PMID: 34937363 DOI: 10.1063/5.0076010] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
A new molecular dataset called HAB79 is introduced to provide ab initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is composed of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular, sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings, whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob's ladder, giving mean relative unsigned errors (MRUEs) of 27.7% (26.3%) for the generalized gradient approximation (GGA) functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP, and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost-effective POD in combination with a GGA functional and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude. The ab initio reference data reported here are expected to be useful for benchmarking other DFT or semi-empirical approaches for electronic coupling calculations.
Collapse
Affiliation(s)
- Orestis George Ziogos
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Adam Kubas
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland
| | - Zdenek Futera
- Faculty of Science, University of South Bohemia, 370 05 Ceske Budejovice, Czech Republic
| | - Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
3
|
Fluctuating hydrogen-bond networks govern anomalous electron transfer kinetics in a blue copper protein. Proc Natl Acad Sci U S A 2018; 115:6129-6134. [PMID: 29844178 PMCID: PMC6004490 DOI: 10.1073/pnas.1805719115] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
We combine experimental and computational methods to address the anomalous kinetics of long-range electron transfer (ET) in mutants of Pseudomonas aeruginosa azurin. ET rates and driving forces for wild type (WT) and three N47X mutants (X = L, S, and D) of Ru(2,2'-bipyridine)2 (imidazole)(His83) azurin are reported. An enhanced ET rate for the N47L mutant suggests either an increase of the donor-acceptor (DA) electronic coupling or a decrease in the reorganization energy for the reaction. The underlying atomistic features are investigated using a recently developed nonadiabatic molecular dynamics method to simulate ET in each of the azurin mutants, revealing unexpected aspects of DA electronic coupling. In particular, WT azurin and all studied mutants exhibit more DA compression during ET (>2 Å) than previously recognized. Moreover, it is found that DA compression involves an extended network of hydrogen bonds, the fluctuations of which gate the ET reaction, such that DA compression is facilitated by transiently rupturing hydrogen bonds. It is found that the N47L mutant intrinsically disrupts this hydrogen-bond network, enabling particularly facile DA compression. This work, which reveals the surprisingly fluctional nature of ET in azurin, suggests that hydrogen-bond networks can modulate the efficiency of long-range biological ET.
Collapse
|
4
|
Computational evidence support the hypothesis of neuroglobin also acting as an electron transfer species. J Biol Inorg Chem 2017; 22:615-623. [DOI: 10.1007/s00775-017-1455-2] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Accepted: 03/28/2017] [Indexed: 12/31/2022]
|
5
|
Gillet N, Berstis L, Wu X, Gajdos F, Heck A, de la Lande A, Blumberger J, Elstner M. Electronic Coupling Calculations for Bridge-Mediated Charge Transfer Using Constrained Density Functional Theory (CDFT) and Effective Hamiltonian Approaches at the Density Functional Theory (DFT) and Fragment-Orbital Density Functional Tight Binding (FODFTB) Level. J Chem Theory Comput 2016; 12:4793-4805. [DOI: 10.1021/acs.jctc.6b00564] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Natacha Gillet
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse
12, 76131 Karlsruhe, Germany
| | - Laura Berstis
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| | - Xiaojing Wu
- Laboratoire
de Chimie-Physique, Université Paris Sud, CNRS, Université Paris Saclay, Campus d’Orsay. 15, avenue Jean Perrin, 91405 Cedex Orsay, France
| | - Fruzsina Gajdos
- Department
of Physics and Astronomy, University College London, Gower Street, London WCIE 6BT, United Kingdom
| | - Alexander Heck
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse
12, 76131 Karlsruhe, Germany
| | - Aurélien de la Lande
- Laboratoire
de Chimie-Physique, Université Paris Sud, CNRS, Université Paris Saclay, Campus d’Orsay. 15, avenue Jean Perrin, 91405 Cedex Orsay, France
| | - Jochen Blumberger
- Department
of Physics and Astronomy, University College London, Gower Street, London WCIE 6BT, United Kingdom
| | - Marcus Elstner
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse
12, 76131 Karlsruhe, Germany
| |
Collapse
|
6
|
Park JW, Rhee YM. Diabatic Population Matrix Formalism for Performing Molecular Mechanics Style Simulations with Multiple Electronic States. J Chem Theory Comput 2015; 10:5238-53. [PMID: 26583208 DOI: 10.1021/ct5006856] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An accurate description of nonbonded interactions is important in investigating dynamics of molecular systems. In many situations, fixed point charge models are successfully applied to explaining various chemical phenomena. However, these models with conventional formulations will not be appropriate in elucidating the detailed dynamics during nonadiabatic events. This is mainly because the chemical properties of any molecule, especially its electronic populations, significantly change with respect to molecular distortions in the vicinity of the surface crossing. To overcome this issue in molecular simulations yet within the framework of the fixed point charge model, we define a diabatic electronic population matrix and substitute it for the conventional adiabatic partial charges. We show that this matrix can be readily utilized toward attaining more reliable descriptions of Coulombic interactions, in combination with the interpolation formalism for obtaining the intramolecular interaction potential. We demonstrate how the mixed formalism with the diabatic charges and the interpolation can be applied to molecular simulations by conducting adiabatic and nonadiabatic molecular dynamics trajectory calculations of the green fluorescent protein chromophore anion in aqueous environment.
Collapse
Affiliation(s)
- Jae Woo Park
- Center for Self-Assembly and Complexity, Institute for Basic Science (IBS) , Pohang 790-784, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 790-784, Korea
| | - Young Min Rhee
- Center for Self-Assembly and Complexity, Institute for Basic Science (IBS) , Pohang 790-784, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 790-784, Korea
| |
Collapse
|
7
|
Blumberger J. Recent Advances in the Theory and Molecular Simulation of Biological Electron Transfer Reactions. Chem Rev 2015; 115:11191-238. [DOI: 10.1021/acs.chemrev.5b00298] [Citation(s) in RCA: 238] [Impact Index Per Article: 26.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jochen Blumberger
- Department of Physics and
Astronomy, University College London, Gower Street, London WC1E 6BT, U.K
| |
Collapse
|
8
|
Kubas A, Gajdos F, Heck A, Oberhofer H, Elstner M, Blumberger J. Electronic couplings for molecular charge transfer: benchmarking CDFT, FODFT and FODFTB against high-level ab initio calculations. II. Phys Chem Chem Phys 2015; 17:14342-54. [PMID: 25573447 DOI: 10.1039/c4cp04749d] [Citation(s) in RCA: 102] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
A new database (HAB7-) of electronic coupling matrix elements (Hab) for electron transfer in seven medium-sized negatively charged π-conjugated organic dimers is introduced. Reference data are obtained with spin-component scaled approximate coupled cluster method (SCS-CC2) and large basis sets. Assessed DFT-based approaches include constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), self-consistent charge density functional tight-binding (FODFTB) and the recently described analytic overlap method (AOM). This complements the previously reported HAB11 database where only cationic dimers were considered. The CDFT method in combination with a functional based on PBE and including 50% of exact exchange (HFX) was found to provide best estimates, with a mean relative unsigned error (MRUE) of 8.2%. CDFT couplings systematically increase with decreasing fraction of HFX as a consequence of increasing delocalisation of the SOMO orbital. The FODFT method is found to be very robust underestimating electronic couplings by 28%. The FODFTB and AOM methods, although orders of magnitude more efficient in terms of computational effort than the DFT approaches, perform well with reasonably small errors of 54% and 29%, respectively, translating in errors in the non-adiabatic electron transfer rate of a factor of 2.4 and 1.7, respectively. We discuss carefully various sources of errors and the scope and limitations of all assessed methods taking into account the results obtained for both HAB7- and HAB11 databases.
Collapse
Affiliation(s)
- Adam Kubas
- University College London, Department of Physics and Astronomy, Gower Street, London WC1E 6BT, UK.
| | | | | | | | | | | |
Collapse
|
9
|
de la Lande A, Gillet N, Chen S, Salahub DR. Progress and challenges in simulating and understanding electron transfer in proteins. Arch Biochem Biophys 2015; 582:28-41. [PMID: 26116376 DOI: 10.1016/j.abb.2015.06.016] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2015] [Revised: 06/15/2015] [Accepted: 06/22/2015] [Indexed: 11/19/2022]
Abstract
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.
Collapse
Affiliation(s)
- Aurélien de la Lande
- Laboratoire de Chimie Physique, UMR 8000, CNRS, Université Paris Sud. 15, av. Jean Perrin, 91405 Orsay, France.
| | - Natacha Gillet
- Laboratoire de Chimie Physique, UMR 8000, CNRS, Université Paris Sud. 15, av. Jean Perrin, 91405 Orsay, France
| | - Shufeng Chen
- Laboratoire de Chimie Physique, UMR 8000, CNRS, Université Paris Sud. 15, av. Jean Perrin, 91405 Orsay, France
| | - Dennis R Salahub
- Department of Chemistry, CMS - Centre for Molecular Simulation and IQST - Institute for Quantum Science and Technology, University of Calgary, 2500 University Drive N.W., Calgary, Alberta T2N 1N4, Canada.
| |
Collapse
|
10
|
Kubas A, Hoffmann F, Heck A, Oberhofer H, Elstner M, Blumberger J. Electronic couplings for molecular charge transfer: benchmarking CDFT, FODFT, and FODFTB against high-level ab initio calculations. J Chem Phys 2014; 140:104105. [PMID: 24628150 DOI: 10.1063/1.4867077] [Citation(s) in RCA: 127] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
We introduce a database (HAB11) of electronic coupling matrix elements (H(ab)) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H(ab) values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.
Collapse
Affiliation(s)
- Adam Kubas
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Felix Hoffmann
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Alexander Heck
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 6, 76131 Karlsruhe, Germany
| | - Harald Oberhofer
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 6, 76131 Karlsruhe, Germany
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
11
|
Migliore A, Polizzi NF, Therien M, Beratan DN. Biochemistry and theory of proton-coupled electron transfer. Chem Rev 2014; 114:3381-465. [PMID: 24684625 PMCID: PMC4317057 DOI: 10.1021/cr4006654] [Citation(s) in RCA: 350] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2013] [Indexed: 02/01/2023]
Affiliation(s)
- Agostino Migliore
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Nicholas F. Polizzi
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - Michael
J. Therien
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| | - David N. Beratan
- Department
of Chemistry, Department of Biochemistry, and Department of Physics, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
12
|
Brancolini G, Migliore A, Corni S, Fuentes-Cabrera M, Luque FJ, Di Felice R. Dynamical treatment of charge transfer through duplex nucleic acids containing modified adenines. ACS NANO 2013; 7:9396-406. [PMID: 24060008 PMCID: PMC3903158 DOI: 10.1021/nn404165y] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
We address the issue of whether chemical alterations of nucleobases are an effective tool to modulate charge transfer through DNA molecules. Our investigation uses a multilevel computational approach based on classical molecular dynamics and quantum chemistry. We find yet another piece of evidence that structural fluctuations are a key factor to determine the electronic structure of double-stranded DNA. We argue that the electronic structure and charge transfer ability of flexible polymers is the result of a complex intertwining of various structural, dynamical and chemical factors. Chemical intuition may be used to design molecular wires, but this is not the sole component in the complex charge transfer mechanism through DNA.
Collapse
Affiliation(s)
- Giorgia Brancolini
- CNR Institute of Nanoscience, S3 Center, Via Campi 213/A, 41125 Modena, Italy
- (GB); (RDF). Phone: +39-059-205-5320. Fax: +39-059-205-5651
| | - Agostino Migliore
- School of Chemistry, Tel Aviv University, 69978 Tel Aviv, Israel
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Stefano Corni
- CNR Institute of Nanoscience, S3 Center, Via Campi 213/A, 41125 Modena, Italy
| | - Miguel Fuentes-Cabrera
- Center for Nanophase Materials Sciences, and Computer Science and Mathematics Division, Oak Ridge National Laboratory, P O Box 2008, Oak Ridge, Tennessee 37831 6494, USA
| | - F. Javier Luque
- Department de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avenida Diagonal 643, Barcelona 08028, Spain
| | - Rosa Di Felice
- CNR Institute of Nanoscience, S3 Center, Via Campi 213/A, 41125 Modena, Italy
- Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089 USA
- (GB); (RDF). Phone: +39-059-205-5320. Fax: +39-059-205-5651
| |
Collapse
|
13
|
Chen Z, Mo Y. Electron Transfer in Electrophilic Aromatic Nitration and Nitrosation: Computational Evidence for the Marcus Inverted Region. J Chem Theory Comput 2013; 9:4428-35. [DOI: 10.1021/ct400618k] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zhenhua Chen
- Department
of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, United States
| | - Yirong Mo
- Department
of Chemistry, Western Michigan University, Kalamazoo, Michigan 49008, United States
| |
Collapse
|
14
|
Oberhofer H, Blumberger J. Revisiting electronic couplings and incoherent hopping models for electron transport in crystalline C60 at ambient temperatures. Phys Chem Chem Phys 2012; 14:13846-52. [PMID: 22858858 DOI: 10.1039/c2cp41348e] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
We assess the validity of incoherent hopping models that have previously been used to describe electron transport in crystalline C(60) at room temperature. To this end we present new density functional theory based calculations of the electron transfer parameter defining these models. Specifically, we report electronic coupling matrix elements for several ten thousand configurations that are thermally accessible to the C(60) molecules through rotational diffusion around their lattice sites. We find that the root-mean-square fluctuations of the electronic coupling matrix element (11 meV) are almost as large as the average value (14 meV) and that the distribution is well approximated by a Gaussian. Importantly, due to the small reorganisation energy of the C(60) dimer (≈0.1 eV), the ET is almost activationless for the majority of configurations. Yet, for a small but significant fraction of orientations the coupling is so strong compared to reorganisation energy that no charge-localised state exists, a situation that is aggravated if zero-point motion of the nuclei is taken into account. The present calculations indicate that standard hopping models do not provide a sound description of electron transport in C(60), which might be the case for many other organics as well, and that approaches are needed that solve the electron dynamics directly.
Collapse
Affiliation(s)
- Harald Oberhofer
- Theoretical Chemistry, TU Munich, Lichtenbergstr. 4, D-85747 Garching, Germany
| | | |
Collapse
|
15
|
Nazmutdinov RR, Bronshtein MD, Zinkicheva TT, Chi Q, Zhang J, Ulstrup J. Modeling and computations of the intramolecular electron transfer process in the two-heme protein cytochrome c(4). Phys Chem Chem Phys 2012; 14:5953-65. [PMID: 22430606 DOI: 10.1039/c2cp24084j] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The di-heme protein Pseudomonas stutzeri cytochrome c(4) (cyt c(4)) has emerged as a useful model for studying long-range protein electron transfer (ET). Recent experimental observations have shown a dramatically different pattern of intramolecular ET between the two heme groups in different local environments. Intramolecular ET in homogeneous solution is too slow (>10 s) to be detected but fast (ms-μs) intramolecular ET in an electrochemical environment has recently been achieved by controlling the molecular orientation of the protein assembled on a gold electrode surface. In this work we have performed computational modeling of the intramolecular ET process by a combination of density functional theory (DFT) and quantum mechanical charge transfer theory to disclose reasons for this difference. We first address the electronic structures of the model heme core with histidine and methionine axial ligands in both low- and high-spin states by structure-optimized DFT. The computations enable estimating the intramolecular reorganization energy of the ET process for different combinations of low- and high-spin heme couples. Environmental reorganization free energies, work terms ("gating") and driving force were determined using dielectric continuum models. We then calculated the electronic transmission coefficient of the intramolecular ET rate using perturbation theory combined with the electronic wave functions determined by the DFT calculations for different heme group orientations and Fe-Fe separations. The reactivity of low- and high-spin heme groups was notably different. The ET rate is exceedingly low for the crystallographic equilibrium orientation but increases by several orders of magnitude for thermally accessible non-equilibrium configurations. Deprotonation of the propionate carboxyl group was also found to enhance the ET rate significantly. The results are discussed in relation to the observed surface immobilization effect and support the notion of conformationally gated ET.
Collapse
Affiliation(s)
- Renat R Nazmutdinov
- Kazan National Research Technological University, 420015 Kazan, Republic Tatarstan, Russian Federation
| | | | | | | | | | | |
Collapse
|
16
|
Mo Y, Song L, Lin Y, Liu M, Cao Z, Wu W. Block-Localized Wavefunction (BLW) Based Two-State Approach for Charge Transfers between Phenyl Rings. J Chem Theory Comput 2012; 8:800-5. [DOI: 10.1021/ct2008318] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Yirong Mo
- Department of Chemistry, Western
Michigan University, Kalamazoo, Michigan 49008, United States
| | - Lingchun Song
- Department of Chemistry, Western
Michigan University, Kalamazoo, Michigan 49008, United States
| | - Yuchun Lin
- Department of Chemistry, Western
Michigan University, Kalamazoo, Michigan 49008, United States
| | - Minghong Liu
- Department of Chemistry, Western
Michigan University, Kalamazoo, Michigan 49008, United States
| | - Zexing Cao
- The State Key Laboratory of
Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory
of Theoretical and Computational Chemistry, and College of Chemistry
and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005,
China
| | - Wei Wu
- The State Key Laboratory of
Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory
of Theoretical and Computational Chemistry, and College of Chemistry
and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005,
China
| |
Collapse
|
17
|
Affiliation(s)
- Benjamin Kaduk
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA
| | | | | |
Collapse
|
18
|
Migliore A. Nonorthogonality Problem and Effective Electronic Coupling Calculation: Application to Charge Transfer in π-Stacks Relevant to Biochemistry and Molecular Electronics. J Chem Theory Comput 2011; 7:1712-25. [PMID: 26596435 DOI: 10.1021/ct200192d] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recently proposed method for the calculation of the effective electronic coupling (or charge-transfer integral) in a two-state system is discussed and related to other methods in the literature. The theoretical expression of the coupling is exact within the two-state model and applies to the general case where the charge transfer (CT) process involves nonorthogonal initial and final diabatic (localized) states. In this work, it is shown how this effective electronic coupling is also the one to be used in a suitable extension of Rabi's formula to the nonorthogonal representation of two-state dynamical problems. The formula for the transfer integral is inspected in the regime of long-range CT and applied to CT reactions in redox molecular systems of interest to biochemistry and/or to molecular electronics: the guanine-thymine stack from regular B-DNA, the polyaromatic perylenediimide stack, and the quinol-semiquinone couple. The calculations are performed within the framework of the Density Functional Theory (DFT), using hybrid exchange-correlation (XC) density functionals, which also allowed investigation of the appropriateness of such hybrid-DFT methods for computing electronic couplings. The use of the recently developed M06-2X and M06-HF density functionals in appropriate ways is supported by the results of this work.
Collapse
|
19
|
Oberhofer H, Blumberger J. Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set. J Chem Phys 2011; 133:244105. [PMID: 21197974 DOI: 10.1063/1.3507878] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a plane wave basis set implementation for the calculation of electronic coupling matrix elements of electron transfer reactions within the framework of constrained density functional theory (CDFT). Following the work of Wu and Van Voorhis [J. Chem. Phys. 125, 164105 (2006)], the diabatic wavefunctions are approximated by the Kohn-Sham determinants obtained from CDFT calculations, and the coupling matrix element calculated by an efficient integration scheme. Our results for intermolecular electron transfer in small systems agree very well with high-level ab initio calculations based on generalized Mulliken-Hush theory, and with previous local basis set CDFT calculations. The effect of thermal fluctuations on the coupling matrix element is demonstrated for intramolecular electron transfer in the tetrathiafulvalene-diquinone (Q-TTF-Q(-)) anion. Sampling the electronic coupling along density functional based molecular dynamics trajectories, we find that thermal fluctuations, in particular the slow bending motion of the molecule, can lead to changes in the instantaneous electron transfer rate by more than an order of magnitude. The thermal average, (<|H(ab)|(2)>)(1/2)=6.7 mH, is significantly higher than the value obtained for the minimum energy structure, |H(ab)|=3.8 mH. While CDFT in combination with generalized gradient approximation (GGA) functionals describes the intermolecular electron transfer in the studied systems well, exact exchange is required for Q-TTF-Q(-) in order to obtain coupling matrix elements in agreement with experiment (3.9 mH). The implementation presented opens up the possibility to compute electronic coupling matrix elements for extended systems where donor, acceptor, and the environment are treated at the quantum mechanical (QM) level.
Collapse
Affiliation(s)
- Harald Oberhofer
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | | |
Collapse
|
20
|
de la Lande A, Salahub DR. Derivation of interpretative models for long range electron transfer from constrained density functional theory. ACTA ACUST UNITED AC 2010. [DOI: 10.1016/j.theochem.2009.11.012] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
21
|
Migliore A. Full-electron calculation of effective electronic couplings and excitation energies of charge transfer states: Application to hole transfer in DNA pi-stacks. J Chem Phys 2010; 131:114113. [PMID: 19778106 PMCID: PMC2766402 DOI: 10.1063/1.3232007] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work I develop and apply a theoretical method for calculating effective electronic couplings (or transfer integrals) between redox sites involved in hole or electron transfer reactions. The resulting methodology is a refinement and a generalization of a recently developed approach for transfer integral evaluation. In fact, it holds for any overlap between the charge-localized states used to represent charge transfer (CT) processes in the two-state model. The presented theoretical and computational analyses show that the prototype approach is recovered for sufficiently small overlaps. The method does not involve any empirical parameter. It allows a complete multielectron description, therefore including electronic relaxation effects. Furthermore, its theoretical formulation holds at any value of the given reaction coordinate and yields a formula for the evaluation of the vertical excitation energy (i.e., the energy difference between the adiabatic ground and first-excited electronic states) that rests on the same physical quantities used in transfer integral calculation. In this paper the theoretical approach is applied to CT in B-DNA base dimers within the framework of Density Functional Theory (DFT), although it can be implemented in other computational schemes. The results of this work, as compared with previous Hartree-Fock (HF) and post-HF evaluations, support the applicability of the current implementation of the method to larger pi-stacked arrays, where post-HF approaches are computationally unfeasible.
Collapse
Affiliation(s)
- Agostino Migliore
- Department of Chemistry and Center for Molecular Modeling, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, USA
| |
Collapse
|
22
|
Vázquez-Mayagoitia Á, Huertas O, Brancolini G, Migliore A, Sumpter BG, Orozco M, Luque FJ, Di Felice R, Fuentes-Cabrera M. Ab initio Study of the Structural, Tautomeric, Pairing, and Electronic Properties of Seleno-Derivatives of Thymine. J Phys Chem B 2009; 113:14465-72. [DOI: 10.1021/jp9057077] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Álvaro Vázquez-Mayagoitia
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Oscar Huertas
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Giorgia Brancolini
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Agostino Migliore
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Bobby G. Sumpter
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Modesto Orozco
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - F. Javier Luque
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Rosa Di Felice
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| | - Miguel Fuentes-Cabrera
- Chemistry Department, University of Tennessee, 1416 Circle Drive, 552 Dabney-Buehler Hall, Knoxville, Tennessee 37996-1600, Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultat de Farmàcia, Universitat de Barcelona, Avgda Diagonal 643, Barcelona, 08028, Spain, National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Via Campi 213/A, 41125 Modena, Italy, Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street,
| |
Collapse
|
23
|
Migliore A, Corni S, Varsano D, Klein ML, Di Felice R. First principles effective electronic couplings for hole transfer in natural and size-expanded DNA. J Phys Chem B 2009; 113:9402-15. [PMID: 19537767 DOI: 10.1021/jp904295q] [Citation(s) in RCA: 59] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Hole transfer processes between base pairs in natural DNA and size-expanded DNA (xDNA) are studied and compared, by means of an accurate first principles evaluation of the effective electronic couplings (also known as transfer integrals), in order to assess the effect of the base augmentation on the efficiency of charge transport through double-stranded DNA. According to our results, the size expansion increases the average electronic coupling, and thus the CT rate, with potential implications in molecular biology and in the implementation of molecular nanoelectronics. Our analysis shows that the effect of the nucleobase expansion on the charge-transfer (CT) rate is sensitive to the sequence of base pairs. Furthermore, we find that conformational variability is an important factor for the modulation of the CT rate. From a theoretical point of view, this work offers a contribution to the CT chemistry in pi-stacked arrays. Indeed, we compare our methodology against other standard computational frameworks that have been adopted to tackle the problem of CT in DNA, and unravel basic principles that should be accounted for in selecting an appropriate theoretical level.
Collapse
Affiliation(s)
- Agostino Migliore
- Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104-6323, USA.
| | | | | | | | | |
Collapse
|
24
|
Guallar V, Wallrapp F. Mapping protein electron transfer pathways with QM/MM methods. J R Soc Interface 2009; 5 Suppl 3:S233-9. [PMID: 18445553 DOI: 10.1098/rsif.2008.0061.focus] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Mixed quantum mechanics/molecular mechanics (QM/MM) methods offer a valuable computational tool for understanding the electron transfer pathway in protein-substrate interactions and protein-protein complexes. These hybrid methods are capable of solving the Schrödinger equation on a small subset of the protein, the quantum region, describing its electronic structure under the polarization effects of the remainder of the protein. By selectively turning on and off different residues in the quantum region, we are able to obtain the electron pathway for short- and large-range interactions. Here, we summarize recent studies involving the protein-substrate interaction in cytochrome P450 camphor, ascorbate peroxidase and cytochrome c peroxidase, and propose a novel approach for the long-range protein-protein electron transfer. The results on ascorbate peroxidase and cytochrome c peroxidase reveal the importance of the propionate groups in the electron transfer pathway. The long-range protein-protein electron transfer has been studied on the cytochrome c peroxidase-cytochrome c complex. The results indicate the importance of Phe82 and Cys81 on cytochrome c, and of Asn196, Ala194, Ala176 and His175 on cytochrome c peroxidase.
Collapse
Affiliation(s)
- Victor Guallar
- Life Science Department, Barcelona Supercomputing Center, Jordi Girona 29, Barcelona, Spain.
| | | |
Collapse
|
25
|
Migliore A, Sit PHL, Klein ML. Evaluation of Electronic Coupling in Transition-Metal Systems Using DFT: Application to the Hexa-Aquo Ferric−Ferrous Redox Couple. J Chem Theory Comput 2009; 5:307-23. [DOI: 10.1021/ct800340v] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Agostino Migliore
- Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323
| | - Patrick H.-L. Sit
- Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323
| | - Michael L. Klein
- Center for Molecular Modeling and Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323
| |
Collapse
|
26
|
Wallrapp F, Masone D, Guallar V. Electron Transfer in the P450cam/PDX Complex. The QM/MM e-Pathway. J Phys Chem A 2008; 112:12989-94. [DOI: 10.1021/jp803538u] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Frank Wallrapp
- ICREA Research Professor, Life Science Department, Barcelona Supercomputing Center, Jordi Girona, 29, 08034 Barcelona, Spain
| | - Diego Masone
- ICREA Research Professor, Life Science Department, Barcelona Supercomputing Center, Jordi Girona, 29, 08034 Barcelona, Spain
| | - Victor Guallar
- ICREA Research Professor, Life Science Department, Barcelona Supercomputing Center, Jordi Girona, 29, 08034 Barcelona, Spain
| |
Collapse
|
27
|
de la Lande A, Martí S, Parisel O, Moliner V. Long Distance Electron-Transfer Mechanism in Peptidylglycine α-Hydroxylating Monooxygenase: A Perfect Fitting for a Water Bridge. J Am Chem Soc 2007; 129:11700-7. [PMID: 17764178 DOI: 10.1021/ja070329l] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The active sites of copper enzymes have been the subject of many theoretical and experimental investigations from a number of years. Such studies have embraced topics devoted to the modeling of the first coordination sphere at the metallic cations up to the development of biomimetic, or bioinspired, catalytic systems. At least from the theoretical viewpoint, fewer efforts have been dedicated to elucidate how the two copper cations act concertedly in noncoupled dicopper enzymes such as peptidylglycine alpha-hydroxylating monooxygenase (PHM) and dopamine beta-monooxygenase (DbetaM). In these metalloenzymes, an electronic transfer is assumed between the two distant copper cations (11 A). Recent experimental results suggest that this transfer occurs through water molecules, a phenomenon which has been theoretically evidenced to be of high efficiency in the case of cytochrome b5 (Science, 2005, 310, 1311). In the present contribution dedicated to PHM, we overpass the common theoretical approaches dedicated to the electronic and geometrical structures of sites CuM or CuH restricted to their first coordination spheres and aim at directly comparing theoretical results to the experimentally measured activity of the PHM enzyme. To achieve this goal, molecular dynamics simulations were performed on wild-type and various mutants of PHM. More precisely, we provide an estimate of the electron-transfer efficiency between the CuM and CuH sites by means of such molecular dynamics simulations coupled to Marcus theory joined to the Beratan model to approximate the required coupling matrix elements. The theoretical results are compared to the kinetics measurements performed on wild and mutated PHM. The present work, the dynamic aspects of which are essential, accounts for the experimental results issued from mutagenesis. It supports the conclusion that an electronic transfer can occur between two copper(I) sites along a bridge involving a set of hydrogen and chemical bonds. Residue Gln170 is evidenced to be the keystone of this water-mediated pathway.
Collapse
Affiliation(s)
- Aurélien de la Lande
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castelló, Spain
| | | | | | | |
Collapse
|
28
|
Venkat AS, Corni S, Di Felice R. Electronic coupling between azurin and gold at different protein/substrate orientations. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2007; 3:1431-7. [PMID: 17591735 DOI: 10.1002/smll.200700001] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
By means of constrained classical molecular dynamics simulations, we have computed the structure of azurin deposited on a Au(111) surface at different possible orientations and the azimuthal forces acting on the protein at each sampled conformation. We have then evaluated the effect of the angular variation on the speed of electron tunneling between the protein redox site and the metal surface. We find that the azurin/gold electronic coupling has a strong dependence on the molecular orientation and is greatly enhanced by inclining the protein to lie as flat as possible on the surface. We discuss the implications of our results for scanning probe microscopy experiments in which tunneling currents are measured while the protein is subjected to mechanical forces exerted by the tip of the instrument.
Collapse
Affiliation(s)
- Anurag Setty Venkat
- National Research Center on nanoStructures and bioSystems at Surfaces S3 of INFM-CNR, Modena, Italy
| | | | | |
Collapse
|
29
|
Naraharisetty SRG, Kurochkin DV, Rubtsov IV. C–D Modes as structural reporters via dual-frequency 2DIR spectroscopy. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.02.020] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
30
|
Migliore A, Corni S, Felice RD, Molinari E. Water-Mediated Electron Transfer between Protein Redox Centers. J Phys Chem B 2007; 111:3774-81. [PMID: 17388538 DOI: 10.1021/jp068773i] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Recent experimental and theoretical investigations show that water molecules between or near redox partners can significantly affect their electron-transfer (ET) properties. Here we study the effects of intervening water molecules on the electron self-exchange reaction of azurin (Az), by performing a conformational sampling on the water medium and by using a newly developed ab initio method to calculate transfer integrals between molecular redox sites. We show that the insertion of water molecules at the interface between the copper active sites of Az dimers slightly increases the overall ET rate, while some favorable water conformations can considerably enhance the ET kinetics. These features are traced back to the interplay of two competing factors: the electrostatic interaction between the water and protein subsystems (mainly opposing the ET process for the water arrangements drawn from MD simulations) and the effectiveness of water in mediating ET coupling pathways. Such an interplay provides a physical basis for the found absence of correlation between the electronic couplings derived through ab initio electronic structure calculations and the related quantities obtained through the Empirical Pathways (EP) method. In fact, the latter does not account for electrostatic effects on the transfer integrals. Thus, we conclude that the water-mediated electron tunneling is not controlled by the geometry of a single physical pathway. We discuss the results in terms of the interplay between different ET pathways controlled by the conformational changes of one of the water molecules via its electrostatic influence. Finally, we examine the dynamical effects of the interfacial water and check the validity of the Condon approximation.
Collapse
Affiliation(s)
- Agostino Migliore
- National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Modena, Italy.
| | | | | | | |
Collapse
|
31
|
Migliore A, Corni S, Di Felice R, Molinari E. Water Effects on Electron Transfer in Azurin Dimers. J Phys Chem B 2006; 110:23796-800. [PMID: 17125342 DOI: 10.1021/jp064690q] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recent experimental and theoretical analyses indicate that water molecules between or near redox partners can significantly affect their electron-transfer (ET) properties. Here, we study the effects of intervening water molecules on the electron self-exchange reaction of azurin (Az) by using a newly developed ab-initio method to calculate transfer integrals between molecular sites. We show that the insertion of water molecules in the gap between the copper active sites of Az dimers slows down the exponential decay of the ET rates with the copper-to-copper distance. Depending on the distance between the redox sites, water can enhance or suppress the electron-transfer kinetics. We show that this behavior can be ascribed to the simultaneous action of two competing effects: the electrostatic interaction of water with the protein subsystem and its ability to mediate ET coupling pathways.
Collapse
Affiliation(s)
- Agostino Migliore
- National Center on nanoStructures and bioSystems at Surfaces (S3) of INFM-CNR, Modena, Italy.
| | | | | | | |
Collapse
|
32
|
Wu Q, Van Voorhis T. Extracting electron transfer coupling elements from constrained density functional theory. J Chem Phys 2006; 125:164105. [PMID: 17092061 DOI: 10.1063/1.2360263] [Citation(s) in RCA: 210] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Constrained density functional theory (DFT) is a useful tool for studying electron transfer (ET) reactions. It can straightforwardly construct the charge-localized diabatic states and give a direct measure of the inner-sphere reorganization energy. In this work, a method is presented for calculating the electronic coupling matrix element (Hab) based on constrained DFT. This method completely avoids the use of ground-state DFT energies because they are known to irrationally predict fractional electron transfer in many cases. Instead it makes use of the constrained DFT energies and the Kohn-Sham wave functions for the diabatic states in a careful way. Test calculations on the Zn2+ and the benzene-Cl atom systems show that the new prescription yields reasonable agreement with the standard generalized Mulliken-Hush method. We then proceed to produce the diabatic and adiabatic potential energy curves along the reaction pathway for intervalence ET in the tetrathiafulvalene-diquinone (Q-TTF-Q) anion. While the unconstrained DFT curve has no reaction barrier and gives Hab approximately 17 kcal/mol, which qualitatively disagrees with experimental results, the Hab calculated from constrained DFT is about 3 kcalmol and the generated ground state has a barrier height of 1.70 kcal/mol, successfully predicting (Q-TTF-Q)- to be a class II mixed-valence compound.
Collapse
Affiliation(s)
- Qin Wu
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA.
| | | |
Collapse
|