Biological Charge Transfer in Redox Regulation and Signaling
by
Ruijie Darius Teo
Department of Chemistry Duke University
Date:______Approved:
______David Beratan, Advisor
______Patrick Charbonneau
______Agostino Migliore
______Kenichi Yokoyama
Dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Department of Chemistry in the Graduate School of Duke University
2020
i
v
ABSTRACT
Biological Charge Transfer in Redox Regulation and Signaling
by
Ruijie Darius Teo
Department of Chemistry Duke University
Date:______Approved:
______David Beratan, Advisor
______Patrick Charbonneau
______Agostino Migliore
______Kenichi Yokoyama
An abstract of a dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in the Department of Chemistry in the Graduate School of Duke University
2020
i
v
Copyright by Ruijie Darius Teo 2020
Abstract
Biological signaling via DNA-mediated charge transfer between high-potential
[4Fe4S]2+/3+ clusters is widely discussed in the literature. Recently, it was proposed that for DNA replication on the lagging strand, primer handover from primase to polymerase α is facilitated by DNA-mediated charge transfer between the [4Fe4S] clusters housed in the respective C-terminal domains of the proteins. Using a theoretical-computational approach, I established that redox signaling between the clusters in primase and polymerase α cannot be accomplished solely by DNA-mediated charge transport, due to the unidirectionality of charge transfer between the [4Fe4S] cluster and the nucleic acid. I extended the study by developing an open-source electron hopping pathway search code to characterize hole hopping pathways in proteins and nucleic acids. I used this module to analyze protective hole escape routes in cytochrome p450, cytochrome c oxidase, and benzylsuccinate synthase. Next, I used the module to analyze molecular dynamics snapshots of a mutant primase, where the Y345C mutation
(found in gastric tumors) attenuates charge transfer between the [4Fe4S] cluster and nucleic acid, which in turn, could disrupt the signaling process between primase and polymerase α. In another protein-nucleic acid system, I found that charge transfer in the p53-DNA complex plays an important role for p53 to differentiate Gadd45 DNA and p21 DNA in metabolic pathway regulation. Using density functional theory calculations
iv
on molecular dynamics snapshots, I found that hole transfer (HT) from Gadd45 DNA to the proximal cysteine residue in the DNA-binding domain of p53 is preferred over HT from p21 DNA to cysteine. This preference ensures that the p21 DNA remains bound to the transcription factor p53 which induces the transcription of the gene under cellular oxidative stress. This dissertation concludes with a study that demonstrates similar electron conductivities between an artificial nucleic acid, 2'-deoxy-2'-fluoro- arabinonucleic acid (2’F-ANA), and DNA. Compared to DNA, 2’F-ANA offers the additional benefit of chemical stability with respect to hydrolysis and nuclease degradation, thereby promoting its use as a sensor in biological systems and cellular environments.
v
Contents
Abstract ...... iv
List of Tables ...... x
List of Figures ...... xiii
Acknowledgements ...... xvi
1. Introduction ...... 1
1.1 Charge Transfer in Proteins ...... 2
1.2 Charge Transfer in DNA ...... 4
1.3 Charge Transfer in Protein-DNA Systems ...... 5
1.4 Aims and Objectives ...... 7
2. Classical Marcus Theory ...... 8
2.1 Charge Transfer Mechanisms ...... 11
2.1.1 Tunneling ...... 12
2.1.2 Superexchange ...... 13
2.1.3 Flickering Resonance ...... 14
2.1.4 Hopping ...... 16
3. Electronic Couplings ...... 17
3.1 Empirical/Semi-Empirical Models ...... 17
3.1.1 Hopfield Model ...... 17
3.1.2 Pathway Tunneling Model ...... 18
3.1.3 Average Packing Density Model ...... 19
vi
3.2 Electronic Structure Models ...... 20
3.2.1 Generalized Mulliken-Hush Method ...... 20
3.2.2 Effective Hamiltonians ...... 21
3.2.2.1 Localized Orbitals ...... 21
3.2.3 Constrained Density Functional Theory (CDFT) ...... 22
3.2.4 Orbital Splitting ...... 25
4. Reorganization Energy ...... 26
4.1 Marcus’ Two-Sphere Model ...... 26
4.1.1 Inner-Sphere Reorganization Energy ...... 27
4.1.2 Outer-Sphere Reorganization Energy ...... 29
5. Molecular Dynamics ...... 30
5.1 MD Force Fields ...... 31
5.1.1 Seminario Method for Calculating Force Constants ...... 33
5.1.2 Atomic Charges for Coulomb Interactions ...... 34
5.2 Integration Algorithms for System Evolution ...... 35
5.3 Solvent Representation ...... 37
5.4 Periodic Boundary Conditions ...... 39
5.5 Temperature/Pressure Coupling Algorithms ...... 40
6. Kinetic Models and Master Equations for Charge Transfer Dynamics ...... 43
7. Charge Transfer in the Human Primosome ...... 46
7.1 Methodology ...... 49
7.1.1 Docking and Modeling ...... 49
vii
7.1.2 Redox Potentials ...... 51
7.1.3 Reorganization Energies ...... 52
7.1.4 Electronic Couplings ...... 53
7.1.5 Kinetic Model ...... 56
7.2 Unidirectionality ...... 57
7.3 Concluding Remarks ...... 63
8. Mapping Hole Hopping Pathways in Proteins ...... 67
8.1 Description of EHPath ...... 70
8.2 Cytochrome p450 ...... 78
8.3 Cytochrome c Peroxidase ...... 83
8.4 Benzylsuccinate Synthase ...... 85
8.5 Concluding Remarks ...... 87
9. Role of Y345C Mutation on Primase-RNA/DNA Hole Hopping Pathways ...... 89
9.1 Introduction ...... 89
9.2 Broken-Symmetry DFT ...... 90
9.3 Methodology and Results ...... 92
9.3.1 Generation of Truncated Models ...... 92
9.3.2 DFT Calculations ...... 92
9.3.3 Structural Comparison ...... 96
9.3.4 MD Simulations ...... 97
9.3.5 Binding Free Energies ...... 101
9.3.6 Hole Hopping Pathways ...... 102
viii
10. Hole Transfer at the p53-DNA Interface ...... 106
10.1 Methodology ...... 111
10.2 Results and Analysis ...... 113
10.3 Concluding Remarks ...... 118
11. Hole Transfer in 2’-Deoxy-2’-Fluoro-Arabinonucleic Acid ...... 120
11.1 Methodology ...... 122
11.1.1 MD Simulation Setup ...... 123
11.1.2 Generation of Truncated Models ...... 124
11.1.3 Free Energy Change ∆G° ...... 125
11.1.4 Reorganization Energies ...... 125
11.2 Results and Analysis ...... 126
11.2.1 Effective Electronic Couplings ...... 126
11.2.2 Reorganization Energies ...... 128
11.2.3 Charge Transport ...... 130
11.3 Concluding Remarks ...... 133
12. Conclusions ...... 135
Appendix ...... 137
References ...... 153
Biography ...... 177
ix
List of Tables
Table 1: Reorganization energies (λX) associated with the redox group X in electron self- exchange...... 50
Table 2: Values of VEA and E0red of Tyr, Trp and Met...... 52
DFT Table 3: DFT (M11/6-311g**) values of the electronic couplings (VIF ); and DFT electronic couplings “dressed” with the effects of the protein medium...... 54
Table 4: Semiempirical estimates of electronic couplings (VIF); DFT (M11/6-311g**) values of the electronic couplings; and “dressed” DFT electronic couplings...... 55
Table 5: ET paths, steps, parameters, and rate constants for CT mechanisms potentially at play in the p58c-RNA/DNA complex...... 58
Table 6: Fastest 5 hole hopping routes in P450BM3 (the heme is the hole donor, while Y305 or Y334 is the terminal hole acceptor)...... 82
Table 7: Mean residence time τM of the hole for the 5 fastest hole hopping pathways in Ccp1 with selected terminal HT sites identified in ref. 159...... 84
Table 8: Mean residence times τM (eq. 76) and τM,approx (eq. 77) of the hole in the 5 fastest hole-hopping escape pathways of BSS...... 86
Table 9: RMSD comparison of the six optimized [4Fe4S]3+ structures (optimized with the Cosmo model) with the crystal structure...... 97
Table 10: Top hole hopping pathways between the [4Fe4S] cluster in wild-type or Y345C primase, and the nucleic acid...... 102
x
Table 11: Top hole hopping pathways between the [4Fe4S] cluster in wild-type or Y345C primase, and the nucleic acid...... 103
Table 12: RDA of D-A pairs from the top hopping pathways listed in Table 10...... 104
Table 13: RDA of D-A pairs from the top hopping pathways listed in Table 11...... 105
Table 14: Values of mean-square electronic coupling,
Table 15: Ranges of distances (in Å) spanned by the indicated bp-Cys pairs, over the selected MD snapshots, at the DNA contacts with the top and bottom proteins...... 115
Table 16: Mean-square electronic coupling
Table 17: Reorganization energy (λDA) and hole-transfer rate (kDA) values for the indicated base-pair dimers in DNA, using the S1, S2 and S3 dielectric constant sets...... 129
Table 18: Reorganization energy (λDA) and hole-transfer rate (kDA) values for the indicated base-pair dimers in 2’-FANA, using the S1, S2 and S3 dielectric constant sets.130
Table 19: Mean travel time (τ) spent by the hole to traverse the path from TA to a charge drain in contact with CG...... 131
Table 20: Mean travel time (τ) spent by the hole to traverse the path from CG to a charge drain in contact with TA. The notation is the same as in Table 19...... 131
Table A1: Interatomic distances Req (Å) for the six optimized geometries 13+-63+ of [4Fe4S]3+ and the corresponding average values...... 137
xi
Table A2: Bond stretching force constants Kr (kcal mol−1 Å−2) for the six optimized geometries 13+-63+ of [4Fe4S]3+ and the corresponding average values...... 138
Table A3: Angles qeq (°) for the six optimized geometries 13+-63+ of [4Fe4S]3+ and the corresponding average values...... 139
Table A4: Angle bending force constants Kq (kcal mol−1 rad−2) for the six optimized geometries 13+-63+ of [4Fe4S]3+ and the corresponding average values...... 141
Table A5: Table of the interatomic distances Req (Å) for the three optimized geometries 12+-32+ of [4Fe4S]2+ and the corresponding average values...... 143
Table A6: Table of the bond stretching force constants Kr (kcal mol−1 Å−2) for the three optimized geometries 12+-32+ of [4Fe4S]2+ and the corresponding average values...... 144
Table A7: Table of the angles qeq (°) for the three optimized geometries 12+-32+ of [4Fe4S]2+ and the corresponding average values...... 145
Table A8: Table of the angle bending force constants Kq (kcal mol−1 rad−2) for the three optimized geometries 12+-32+ of [4Fe4S]2+ and the corresponding average values...... 147
Table A9: RESP charges for atoms corresponding to the optimized geometries 13+-63+ of [4Fe4S]3+ and the corresponding average values...... 149
Table A10: Table of the RESP charges for atoms corresponding to the three optimized geometries 12+-32+ of [4Fe4S]2+ and the corresponding average values...... 151
xii
List of Figures
Figure 1: Marcus parabolas for CT (in this case, excess electron transfer) from D- (electron on donor) and A...... 8
Figure 2: (A) Superexchange, (B) flickering resonance, and (C) hopping CT mechanisms for a D-1-2-3-A chain comprising of three bridge sites...... 12
Figure 3: Bonded (solid lines) and non-bonded interactions (dashed line) represented by atoms 1-5...... 32
Figure 4: (A) 3-site water model, (B) 4-site model, (C) 5-site model, and (D) 6-site model. Color code: pink (O), green (H), blue (lone pair), orange (dummy atom)...... 38
Figure 5: Kinetic model for a hopping network consisting of N+1 sites, where the charge starts at site 1 and eventually arrives at the trap site N+1 shown in orange...... 44
Figure 6: Human primosome-nucleic acid complex...... 47
Figure 7: Proposed mechanism of primer handoff driven by DNA charge transport chemistry...... 48
Figure 8: Strongest tunneling pathways between the [4Fe4S] cluster in primase and purine nucleobases in anchored nucleic acid...... 57
Figure 9: CT steps and rate constants in hole hopping between the [4Fe4S] cluster and the nucleobases in the p58c-RNA/DNA complex...... 59
Figure 10: Fastest ET routes from the iron-sulfur cluster in primase to DNA/RNA (i.e., hole transfer in the opposite direction)...... 61
Figure 11: CT-mediated [4Fe4S] protein signaling and pertinent redox potentials...... 62
xiii
Figure 12. Side-by-side comparison of p58C structures with correctly folded [PDB code 5F0Q] and misfolded (PDB code 3L9Q) substrate-binding regions...... 65
Figure 13. Structure of a directed graph representing hopping routes in a protein...... 71
Figure 14: Part of the catalytic cycle of cytochrome P450...... 79
Figure 15: Probable routes for hole hopping from the active site of P450BM3 to the protein surface...... 80
Figure 16: Redox-active residues (C, M, W, Y) in Ccp1 (PDB 1ZBY184)...... 83
Figure 17: Locations of the BSS redox-active residues involved in the five most rapid hole hopping routes from G829• to the terminal hole acceptor [4Fe4S]2+ (Table 8)...... 86
Figure 18: a-h) Optimized geometries 13+-63+ of the [4Fe4S]3+ cluster corresponding to the six spin layer assignments, and g) average geometry of the [4Fe4S]3+ cluster...... 94
Figure 19: RMSD (Å) of the wild-type primase-RNA/DNA complex (excluding water and counterions) across 120 ns using the A3+ average force field parameters...... 99
Figure 20: RMSD (Å) of the wild-type primase-RNA/DNA complex (excluding water and counterions) across 130 ns using the 33+ force field parameters...... 99
Figure 21: RMSD (Å) of the mutated Y345C primase-RNA/DNA complex (excluding water and counterions) across 125 ns using the A3+ force field parameters...... 100
Figure 22: RMSD (Å) of the mutated Y345C primase-RNA/DNA complex (excluding water and counterions) across 110 ns using the 33+ force field parameters...... 100
Figure 23: RMSD (Å) of the wild-type primase-RNA/DNA complex (excluding water and counterions) across 120 ns using the A2+ average force field parameters...... 101
xiv
Figure 24: Protein-DNA models used in the MD simulation of the p53 protein complexed with the Gadd45 and p21 half-site sequences...... 108
Figure 25: RMSD (without hydrogens) along the MD production run for the protein- DNA complexes containing the Gadd45 (blue) and p21 (orange) DNA models...... 112
Figure 26: HT models for the contacts of p53 with (a) Gadd45 and (b) p21 DNA sequences...... 113
Figure 27: (a) VCG-GC, (b) VGC-A1T, (c) VA1T-A2T, and (d) VA2T-TA versus the MD simulation time for DNA and 2’F-ANA...... 126
Figure 28: Instantaneous value of the Marcus expression for the outer-sphere reorganization energy λoDA(t) vs. the MD simulation time (ns) for the D-A dimers...... 128
xv
Acknowledgements
Firstly, I am indebted to my advisor, Professor David Beratan, for his tremendous support, encouragement, and mentorship throughout my graduate education. His kindness, humility, and generosity are some of the traits I would always remember and strive to emulate. Next, I am deeply grateful to my mentor, Professor
Agostino Migliore, for his constant guidance, thoughtfulness, and camaraderie. His passion for research is certainly infectious, and I have become a much better researcher and thinker under his tutelage.
I would like to specially thank my other committee members - Professor Patrick
Charbonneau and Professor Kenichi Yokoyama - for their generous time, advice, and support. I would also like to thank Meg Avery, Michael Conti, Professor Jiyong Hong,
Professor Michael Fitzgerald, and Professor Katherine Franz for their help with administrative matters and fellowship applications. In addition, I am grateful for the interaction with several faculty members and researchers that I had the fortune of meeting as a graduate student, including Professor Stephen Craig, Professor Emily
Derbyshire, Professor Steve Haase, Professor Richard MacPhail, Professor Henry Pfister,
Professor Michael Therien, Dr. Ronald Venters, and Professor Weitao Yang.
I had the privilege of mentoring several students at Duke (Xiaochen Du, Daniel
Koceja, Elizabeth Smithwick, Kiriko Terai, and Hector Torres), and I would like to thank
xvi
them for their valuable research contributions and the opportunity to hone my mentoring skills. I wish them the best in their academic endeavors.
I would also like to thank the past and present members of the Beratan group for enriching my experience. Many thanks to Dr. Shuming Bai, Dr. Tomasz Janowski,
Benjamin Rousseau, Professor Luis Rego, Dr. Xuyan Ru, Jesus Valdiviezo, Jonathon
Yuly, Professor Peng Zhang, Zhendian Zhang, and Dr. Ellie Zheng for making this possible. As busy as graduate school may be, I would like to thank Brandon Bowser,
Heather Folliard, Tony and Faye Hilger, Jiachen Li, Jeffrey Lin, Justin Ma, Melyssa
Minto, Robert Tennant, Ruobing Wang, Yujia Zhai, and Hong Zhou for making the process enjoyable.
In addition, I am grateful to Blue Waters, the Duke Graduate School, and the
National Institute of Health for supporting my research. The Blue Waters staff (special shoutouts to Dr. Victor Anisimov, Dr. Greg Bauer, Dr. Maxim Belkin, Dr. Scott Lathrop,
Noni Ledford, Kjellrun Olson, and Susan Vinson) and the Blue Waters Fellows are an incredible bunch of people to be around with.
Above all, I would like express my immense gratitude to my family and my incredible wife. This dissertation is only possible through her selflessness and sacrificial love. I am also extremely grateful for the joy my precious son brings to me after a long day of research. Daddy loves you so very much. I am truly indebted to both of you.
xvii
1. Introduction
The overarching aim of my research is to understand the intricate relationship among charge transfer (CT), redox regulation, and signaling in biological systems. A CT event occurs when a donor group releases an electron or a hole to the acceptor group.
CT processes govern cellular redox signaling processes, which occur, for example, in response to a change in reactive oxygen species (ROS) or reactive nitrogen species (RNS) levels of the cellular system. The mitochondria, for example, initiate plasma membrane repair in injured cells through their production of ROS facilitated by increased calcium uptake.1 Redox signaling processes are paramount in the regulation of major cellular pathways, while structural differences in proteins and DNA can alter CT pathways, dysregulate cellular signaling, and result in the over proliferation of cells, and human diseases like cancer.
Using theoretical chemistry, modeling, and simulation, I have investigated CT in a wide variety of protein/nucleic acid systems, including the human primase-RNA/DNA complex, p53-DNA complex, azurin, 2’fluoroarabinonucleic acid (2’F-ANA), cytochrome p450, cytochrome c peroxidase (Ccp1), and benzylsuccinate synthase (BSS). As CT in these biological systems helps facilitate the systems’ responses to elevated levels of ROS and
RNS, detailed analyses and characterizations of the kinetics and thermodynamics of CT pathways in these systems are crucial in establishing the role of CT in cellular response.
This chapter aims to introduce and frame CT mechanisms in proteins, DNA, protein-DNA
1
systems, and provide a detailed description of the aims and objectives of this dissertation research.
1.1 Charge Transfer in Proteins
Biological CT is ubiquitous in nature. A classic example is the photosynthetic photosystem II that catalyzes the light-driven oxidation of water. At the core of photosystem II lies the oxygen evolving complex (OEC), an Mn4Ca cofactor that is electronically coupled to the chlorophyll complex P680 through a redox-active tyrosine- histidine pair (TyrZ-D1His190).2 This pair acts as an electron relay between P680 and the
OEC components by mediating four sequential oxidations of the OEC from the S0 state to the S4 state through proton-coupled electron transfer (PCET).2-4 When the S4 state is attained, the oxidation of water, the generation of dioxygen and the regeneration of the S0 state results.
CT in proteins can occur either through single-step tunneling (superexchange, if virtual states of the bridge (B) are present between the charge donor (D) and acceptor (A)) or multi-step hopping. In both cases, the electronic coupling between D and A is relatively weak (see section 1.4 for a more detailed discussion about these mechanisms). This weak coupling imposes a limit of ~20 Å for the D-A distance in order for a single-step tunneling event to be feasible on a biological millisecond timescale.6 Charge hopping, on the other hand, can be viewed as an incoherent multi-step tunneling process where sequential tunneling reactions can deliver charges to distant sites (> 25 Å). Many studies have
2
revealed that electrons are transferred over large distances in enzymes like Complex I
(NADH dehydrogenase) of the mitochondrial respiratory chain5 and in synthetically modified proteins like Ru-cytochrome (cyt) c.6 It has also been proposed that redox-active amino acids like tyrosine and tryptophan can protect proteins from harmful reactive nitrogen and oxygen species by acting as electron or hole conduits to divert highly oxidizing holes from the active site to the protein surface for scavenging.7 The arrangement of these redox-active aromatic residues modulates their electronic coupling interaction energies and could be a key factor in the evolution of hole-hopping chains in proteins.8
Protein CT can be broadly divided into two categories – 1) inter-protein CT, and
2) intra-protein CT. A classic example of inter-protein CT occurs between cytochrome c
(Cc) and cytochrome c oxidase (CcO). A conformational change at the protein interface and a possible redox-dependent gating mechanism9 was attributed to the much lower ET rate from Cc to cytochrome c oxidase CcO (10–102 s−1)10 compared to the rate estimated from electrochemically triggered redox reaction of a surface-modified Cc bound to self- assembled monolayers on a gold electrode (∼107 s−1).11,12 A classic example of intra-protein
CT is in the purple bacterium Rhodobacter sphaeroides R26, where the primary electron donor comprises a pair of symmetrically arranged bacteriochlorophyll a cofactors (PL and
PM) coordinated by different histidine residues (His-L168 for PL and His-M202 for PM).
Other cofactors present include two accessory BChls, two BPhes a, two ubiquinones, and
3
a nonheme iron that are all arranged a near twofold symmetry.13 Interestingly, the electron does not travel along both pathways but instead ET occurs solely along one single branch.
1.2 Charge Transfer in DNA
CT in DNA has been extensively studied. DNA has robust structural, biological, and electronic properties that allow it to function as a charge mediator,14 where the helical nature of DNA and the ordered π-stacking of nucleobases allow charge propagation over nm-scale distances.15 Long-range CT is attractive due to its applications in nanotechnology and biology.16 Optimizing charge transport dynamics through DNA has thus been an intense area of research through the exploration of different hole donors, length of base-pair sequences, as well as the type of sequences.17-24 In biology, extensive research has been focused on DNA-mediated CT in
DNA repair and replication proteins containing redox-active residues like iron-sulfur clusters and cysteine (such as DNA primase and p53) where charges can be shuttled between DNA and these residues. It has been hypothesized that interprotein communication can be established using CT.25-28 In nanotechnology, DNA is, for example, explored in self-assembling 3D architectures. Previous studies, for example, have demonstrated the structural robustness of DNA in a torsionally-stressed ‘tensegrity triangle’ motif.29,30 Engineering such architectures is useful in programmable DNA nanoelectronics based on long-range CT.31
4
DNA CT often involves the migration of a positive charge (a hole) across the
HOMOs of the GC and AT nucleobase pairs. These HOMOs tend to localize on the purines (G or A) due to their lower oxidation potentials (and higher HOMO energies) than the pyrimidines.32 An established mechanism of long-range DNA CT is the hopping mechanism,24,33 where the charge ‘hops’ among the purine bridging sites via multistep tunneling. In contrast to standard coherent tunneling (often used to describe short-range DNA CT), recent theoretical models based on coherent transport, such as ‘deep-hole transfer’34 and transient coherence or flickering resonance (FR),35 were introduced (see section 1.4.4) to describe CT in the hopping regime; the development of these models is motivated by the observed small tunneling decay exponent β (< 1.0 Å-1) at longer distances.36 These coherent mechanisms exist since small energy gaps (< 0.1 eV) between a D and A nucleobase pair can be transiently eliminated by thermal fluctuations (i.e., � ≈ �� � ~ 0.1 eV) at room temperature) of the energy levels. This allows transient delocalization of the charge from D to multiple Bs.
1.3 Charge Transfer in Protein-DNA Systems
Barton’s group has devoted much effort to studying DNA-mediated CT between high potential [4Fe4S]2+/3+-containing DNA replication and repair proteins.16 DNA- mediated CT between [4Fe4S]2+/3+ clusters, which is hypothesized to form the basis for redox signaling between [4Fe4S] cluster-containing proteins,37-39 while one study even
5
suggests that these [4Fe4S] cofactors may have been a key to the origins of life.40 Once hypothesized to play a structural role, recent experiments suggest that these iron-sulfur clusters are involved in the regulation of protein activity.41,42 Switching between the reduced and oxidized redox states of a high-potential [4Fe4S] cluster via DNA-mediated
CT may synchronize the binding and unbinding events required for repair and replication proteins. These events are governed by electrostatics, by which a protein containing a [4Fe4S]3+ cluster binds more tightly to the negatively-charged DNA than a
[4Fe4S]2+-containing protein.43 This increased DNA binding shifts the [4Fe4S]2+/3+ couple to ca. +80 mV vs NHE.26,44,45 From the results of protein activity assays,46,47 [4Fe4S]2+/3+ cluster pairs are proposed to act as points of communication between proteins.16,41
In base excision repair enzymes like Endonuclease III and MutY, it is proposed that mismatched base pairs are recognized by detecting changes in base stacking through DNA CT.44,46,48 Once the [4Fe4S] cluster in one protein becomes oxidized to the
3+ state, it binds more tightly to the DNA and a second protein will attempt to transfer an electron to the first protein. If the double helix is well formed (does not contain base pair mismatches) and the bases are properly stacked, the electron is thought to transfer from a [4Fe4S]2+ cluster in the first protein to the [4Fe4S] 3+ cluster in the second protein through the intervening DNA and protein medium. However, if the base stacking is disrupted by a mismatched base pair or by an amino acid mutation, DNA-mediated CT would be interrupted between these two proteins. Using this method, DNA damage
6
could be located in a more efficient fashion than checking every base pair for defects.
Likewise, DinG, an R-loop repairing helicase, is thought to co-operate with the base excision repair enzymes to detect R-loops in the DNA through a similar mechanism.47
1.4 Aims and Objectives
In order to examine the feasibility of CT in the biological systems presented in this dissertation, the standard methodology of first simulating the biological system of interest using molecular dynamics (MD) , followed by using electronic structure methods like density functional theory (DFT) to calculate Marcus parameters for the CT rate constants (see Chapter 2), is utilized. Chapters 2-6 introduce pertinent theoretical and computational methods that enable these studies. The main objectives for each biological system studied (Chapters 7-11) are outlined as follows – 1) the examination of the directionality of CT between the [4Fe4S]2+/3+-containing primase and the bound
RNA/DNA duplex, 2) development of a computer code that maps hole hopping pathways in proteins, 3) impact of Y345C mutation on CT between primase and nucleic acid duplex, 4) modulation of CT at the p53-nucleic acid binding interface with differing
DNA substrates, and 5) comparison of charge conductivity of an artificial nucleic acid
(2’F-ANA) with DNA.
7
2. Classical Marcus Theory
CT (excess electron transfer, EET, or hole transfer, HT) reactions between a donor and acceptor can be represented by means of two parabolic curves representing the
(free) energies of the initial and final CT states as a function of a nuclear reaction coordinate (see Figure 1). The vertical energy gap ∆�, which is the energy difference between the two adiabatic states (red curves, Figure 1) at a certain nuclear configuration, can represent the nuclear reaction coordinate, although other suitable choices of the reaction coordinate can also be made. Taking the probability distribution � of ∆� to be
Gaussian, this gives rise to two parabolic free energy curves via the Landau free energy
49 equation �(∆�) = −� � ln � (∆�) + constant.
G1
Gb
G0
Figure 1: Marcus parabolas for CT (in this case, excess electron transfer) from D- (electron on donor) and A. The reaction coordinate is the vertical energy gap ΔE. The driving force ΔGo is the free energy difference between the two minima ΔEamin and ΔEbmin. This figure is adapted from ref. 49.
Now, the free energies � (∆�) and � (∆�) for the respective diabatic states a and b relate to the force constant � via the parabolic equation 8
1 � (∆�) = � ∆� − ∆� (1) 2
1 � (∆�) = � ∆� + ∆� (2) 2