1
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
2
|
Narth C, Gillet N, Cailliez F, Lévy B, de la Lande A. Electron transfer, decoherence, and protein dynamics: insights from atomistic simulations. Acc Chem Res 2015; 48:1090-7. [PMID: 25730126 DOI: 10.1021/ar5002796] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Christophe Narth
- Laboratoire
de Chimie Théorique, CNRS UMR 7616, Université Pierre et Marie Curie, case courrier 137. 4, Place Jussieu, 75252 Cedex 05 Paris, France
| | - Natacha Gillet
- Laboratoire
de Chimie-Physique, CNRS UMR 8000, Université Paris Sud, Bâtiment
349 - Campus d’Orsay. 15, avenue Jean Perrin, 91405 Cedex Orsay, France
| | - Fabien Cailliez
- Laboratoire
de Chimie-Physique, CNRS UMR 8000, Université Paris Sud, Bâtiment
349 - Campus d’Orsay. 15, avenue Jean Perrin, 91405 Cedex Orsay, France
| | - Bernard Lévy
- Laboratoire
de Chimie-Physique, CNRS UMR 8000, Université Paris Sud, Bâtiment
349 - Campus d’Orsay. 15, avenue Jean Perrin, 91405 Cedex Orsay, France
| | - Aurélien de la Lande
- Laboratoire
de Chimie-Physique, CNRS UMR 8000, Université Paris Sud, Bâtiment
349 - Campus d’Orsay. 15, avenue Jean Perrin, 91405 Cedex Orsay, France
| |
Collapse
|
3
|
Heck A, Woiczikowski PB, Kubař T, Giese B, Elstner M, Steinbrecher TB. Charge transfer in model peptides: obtaining Marcus parameters from molecular simulation. J Phys Chem B 2012; 116:2284-93. [PMID: 22260641 DOI: 10.1021/jp2086297] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Charge transfer within and between biomolecules remains a highly active field of biophysics. Due to the complexities of real systems, model compounds are a useful alternative to study the mechanistic fundamentals of charge transfer. In recent years, such model experiments have been underpinned by molecular simulation methods as well. In this work, we study electron hole transfer in helical model peptides by means of molecular dynamics simulations. A theoretical framework to extract Marcus parameters of charge transfer from simulations is presented. We find that the peptides form stable helical structures with sequence dependent small deviations from ideal PPII helices. We identify direct exposure of charged side chains to solvent as a cause of high reorganization energies, significantly larger than typical for electron transfer in proteins. This, together with small direct couplings, makes long-range superexchange electron transport in this system very slow. In good agreement with experiment, direct transfer between the terminal amino acid side chains can be dicounted in favor of a two-step hopping process if appropriate bridging groups exist.
Collapse
Affiliation(s)
- Alexander Heck
- Department for Theoretical Chemical Biology, Institute for Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | | | | | | | | | | |
Collapse
|
4
|
Dutta S, Cook RJ, Houtman JC, Kohen A, Cheatum CM. Characterization of azido-NAD+ to assess its potential as a two-dimensional infrared probe of enzyme dynamics. Anal Biochem 2010; 407:241-6. [PMID: 20705046 DOI: 10.1016/j.ab.2010.08.008] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2010] [Accepted: 08/05/2010] [Indexed: 11/23/2022]
Abstract
Enzyme active-site dynamics at femtosecond to picosecond time scales are of great biochemical importance, but remain relatively unexplored due to the lack of appropriate analytical methods. Two-dimensional infrared (2D IR) spectroscopy is one of the few methods that can examine chemical biological motions at this time scale, but all the IR probes used so far were specific to a few unique enzymes. The lack of IR probes of broader specificity is a major limitation to further 2D IR studies of enzyme dynamics. Here we describe the synthesis of a general IR probe for nicotinamide-dependent enzymes. This azido analog of the ubiquitous cofactor nicotinamide adenine dinucleotide is found to be stable and bind to several dehydrogenases with dissociation constants similar to that for the native cofactor. The infrared absorption spectra of this probe bound to several enzymes indicate that it has significant potential as a 2D IR probe to investigate femtosecond dynamics of nicotinamide-dependent enzymes.
Collapse
|
5
|
Abstract
We develop nonorthogonal projectors, called Löwdin projectors, to construct an effective donor-acceptor system composed of localized donor (D) and acceptor (A) states of a long-distance electron transfer problem. When these states have a nonvanishing overlap with the bridge states these projectors are non-Hermitian and there are various possible effective two-level systems that can be built. We show how these can be constructed directly from the Schrödinger or Dyson equation projected onto the D-A subspace of the Hilbert space and explore these equations to determine the connection between Hamiltonian and Green function partitioning. We illustrate the use of these effective two-level systems in estimating the electron transfer rate in the context of a simple electron transfer model.
Collapse
Affiliation(s)
- Paulo C P de Andrade
- Departamento de Física, Universidade Federal do Paraná, 81531-990, Curitiba-PR, Brazil.
| | | |
Collapse
|
6
|
|
7
|
Skourtis SS, Balabin IA, Kawatsu T, Beratan DN. Protein dynamics and electron transfer: electronic decoherence and non-Condon effects. Proc Natl Acad Sci U S A 2005; 102:3552-7. [PMID: 15738409 PMCID: PMC553344 DOI: 10.1073/pnas.0409047102] [Citation(s) in RCA: 155] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
We compute the autocorrelation function of the donor-acceptor tunneling matrix element <T(DA)(t)T(DA)(0)> for six Ru-azurin derivatives. Comparison of this decay time to the decay time of the time-dependent Franck-Condon factor {computed by Rossky and coworkers [Lockwood, D. M., Cheng, Y.-K. & Rossky, P. J. (2001) Chem. Phys. Lett. 345, 159-165]} reveals the extent to which non-Condon effects influence the electron-transfer rate. <T(DA)(t)T(DA)(0)> is studied as a function of donor-acceptor distance, tunneling pathway structure, tunneling energy, and temperature to explore the structural and dynamical origins of non-Condon effects. For azurin, the correlation function is remarkably insensitive to tunneling pathway structure. The decay time is only slightly shorter than it is for solvent-mediated electron transfer in small organic molecules and originates, largely, from fluctuations of valence angles rather than bond lengths.
Collapse
Affiliation(s)
- Spiros S Skourtis
- Department of Physics, University of Cyprus, P.O. Box 20537, Nicosia 1678, Cyprus.
| | | | | | | |
Collapse
|
8
|
Aquino AJ, Beroza P, Reagan J, Onuchic J. Estimating the effect of protein dynamics on electron transfer to the special pair in the photosynthetic reaction center. Chem Phys Lett 1997; 275:181-7. [DOI: 10.1016/s0009-2614(97)00742-2] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
9
|
Aquino A, Beroza P, Beratan D, Onuchic J. Docking and electron transfer between cytochrome c2 and the photosynthetic reaction center. Chem Phys 1995; 197:277-88. [DOI: 10.1016/0301-0104(95)00099-a] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
10
|
|
11
|
Blanchard JP, Blackman CF. Clarification and application of an ion parametric resonance model for magnetic field interactions with biological systems. Bioelectromagnetics 1994; 15:217-38. [PMID: 8074738 DOI: 10.1002/bem.2250150306] [Citation(s) in RCA: 158] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Theoretical models proposed to date have been unable to clearly predict biological results from exposure to low-intensity electric and magnetic fields (EMF). Recently a predictive ionic resonance model was proposed by Lednev, based on an earlier atomic spectroscopy theory described by Podgoretskii and Podgoretskii and Khrustalev. The ion parametric resonance (IPR) model developed in this paper corrects mathematical errors in the earlier Lednev model and extends that model to give explicit predictions of biological responses to parallel AC and DC magnetic fields caused by field-induced changes in combinations of ions within the biological system. Distinct response forms predicted by the IPR model depend explicitly on the experimentally controlled variables: magnetic flux densities of the AC and DC magnetic fields (Bac and Bdc, respectively); AC frequency (fac); and, implicitly, charge to mass ratio of target-ions. After clarifying the IPR model and extending it to combinations of different resonant ions, this paper proposes a basic set of experiments to test the IPR model directly which do not rely on the choice of a particular specimen or endpoint. While the fundamental bases of the model are supported by a variety of other studies, the IPR model is necessarily heuristic when applied to biological systems, because it is based on the premise that the magnitude and form of magnetic field interactions with unhydrated resonant ions in critical biological structures alter ion-associated biological activities that may in turn be correlated with observable effects in living systems.
Collapse
Affiliation(s)
- J P Blanchard
- Rochtel Corporation, San Francisco, California 94119-3965
| | | |
Collapse
|
12
|
|
13
|
|