Interfacial hydration, dynamics and electron transfer: multi-scale ET modeling of the transient [myoglobin, cytochrome b5] complex

Tuning Radical Relay Residues by Proton Management Rescues Protein Electron Hopping

Transient tyrosine and tryptophan radicals play key roles in the electron transfer (ET) reactions of photosystem (PS) II, ribonucleotide reductase (RNR), photolyase, and many other proteins. However, Tyr and Trp are not functionally interchangeable, and the factors controlling their reactivity are often unclear. Cytochrome c peroxidase (CcP) employs a Trp191^•+ radical to oxidize reduced cytochrome c (Cc). Although a Tyr191 replacement also forms a stable radical, it does not support rapid ET from Cc. Here we probe the redox properties of CcP Y191 by non-natural amino acid substitution, altering the ET driving force and manipulating the protic environment of Y191. Higher potential fluorotyrosine residues increase ET rates marginally, but only addition of a hydrogen bond donor to Tyr191^• (via Leu232His or Glu) substantially alters activity by increasing the ET rate by nearly 30-fold. ESR and ESEEM spectroscopies, crystallography, and pH-dependent ET kinetics provide strong evidence for hydrogen bond formation to Y191^• by His232/Glu232. Rate measurements and rapid freeze quench ESR spectroscopy further reveal differences in radical propagation and Cc oxidation that support an increased Y191^• formal potential of ∼200 mV in the presence of E232. Hence, Y191 inactivity results from a potential drop owing to Y191^•+ deprotonation. Incorporation of a well-positioned base to accept and donate back a hydrogen bond upshifts the Tyr^• potential into a range where it can effectively oxidize Cc. These findings have implications for the Y_Z/Y_D radicals of PS II, hole-hopping in RNR and cryptochrome, and engineering proteins for long-range ET reactions.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

Quantum chemical calculations of tryptophan → heme electron and excitation energy transfer rates in myoglobin

The development of optical multidimensional spectroscopic techniques has opened up new possibilities for the study of biological processes. Recently, ultrafast two-dimensional ultraviolet spectroscopy experiments have determined the rates of tryptophan → heme electron transfer and excitation energy transfer for the two tryptophan residues in myoglobin (Consani et al., Science, 2013, 339, 1586). Here, we show that accurate prediction of these rates can be achieved using Marcus theory in conjunction with time-dependent density functional theory. Key intermediate residues between the donor and acceptor are identified, and in particular the residues Val68 and Ile75 play a critical role in calculations of the electron coupling matrix elements. Our calculations demonstrate how small changes in structure can have a large effect on the rates, and show that the different rates of electron transfer are dictated by the distance between the heme and tryptophan residues, while for excitation energy transfer the orientation of the tryptophan residues relative to the heme is important. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Charge-Disproportionation Symmetry Breaking Creates a Heterodimeric Myoglobin Complex with Enhanced Affinity and Rapid Intracomplex Electron Transfer

We report rapid photoinitiated intracomplex electron transfer (ET) within a "charge-disproportionated" myoglobin (Mb) dimer with greatly enhanced affinity. Two mutually supportive Brownian Dynamics (BD) interface redesign strategies, one a new "heme-filtering" approach, were employed to "break the symmetry" of a Mb homodimer by pairing Mb constructs with complementary highly positive and highly negative net surface charges, introduced through D/E → K and K → E mutations, respectively. BD simulations using a previously developed positive mutant, Mb(+6) = Mb(D44K/D60K/E85K), led to construction of the complementary negative mutant Mb(-6) = Mb(K45E, K63E, K95E). Simulations predict the pair will form a well-defined complex comprising a tight ensemble of conformations with nearly parallel hemes, at a metal-metal distance ∼18-19 Å. Upon expression and X-ray characterization of the partners, BD predictions were verified through ET photocycle measurements enabled by Zn-deuteroporphyrin substitution, forming the [ZnMb(-6), Fe(3+)Mb(+6)] complex. Triplet ET quenching shows charge disproportionation increases the binding constant by no less than ∼5 orders of magnitude relative to wild-type Mb values. All progress curves for charge separation (CS) and charge recombination (CR) are reproduced by a generalized kinetic model for the interprotein ET photocycle. The intracomplex ET rate constants for both CS and CR are increased by over 5 orders of magnitude, and their viscosity independence is indicative of true interprotein ET, rather than dynamic gating as seen in previous studies. The complex displays an unprecedented timecourse for CR of the CS intermediate I. After a laser flash, I forms through photoinduced CS, accumulates to a maximum concentration, then dies away through CR. However, before completely disappearing, I reappears without another flash and reaches a second maximum before disappearing completely.

"Bind and Crawl" Association Mechanism of Leishmania major Peroxidase and Cytochrome c Revealed by Brownian and Molecular Dynamics Simulations

Leishmania major, the parasitic causative agent of leishmaniasis, produces a heme peroxidase (LmP), which catalyzes the peroxidation of mitochondrial cytochrome c (LmCytc) for protection from reactive oxygen species produced by the host. The association of LmP and LmCytc, which is known from kinetics measurements to be very fast (∼10(8) M(-1) s(-1)), does not involve major conformational changes and has been suggested to be dominated by electrostatic interactions. We used Brownian dynamics simulations to investigate the mechanism of formation of the LmP-LmCytc complex. Our simulations confirm the importance of electrostatic interactions involving the negatively charged D211 residue at the LmP active site, and reveal a previously unrecognized role in complex formation for negatively charged residues in helix A of LmP. The crystal structure of the D211N mutant of LmP reported herein is essentially identical to that of wild-type LmP, reinforcing the notion that it is the loss of charge at the active site, and not a change in structure, that reduces the association rate of the D211N variant of LmP. The Brownian dynamics simulations further show that complex formation occurs via a "bind and crawl" mechanism, in which LmCytc first docks to a location on helix A that is far from the active site, forming an initial encounter complex, and then moves along helix A to the active site. An atomistic molecular dynamics simulation confirms the helix A binding site, and steady state activity assays and stopped-flow kinetics measurements confirm the role of helix A charges in the association mechanism.

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

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

Electron-transfer acceleration investigated by time resolved infrared spectroscopy

Ultrafast electron transfer (ET) processes are important primary steps in natural and artificial photosynthesis, as well as in molecular electronic/photonic devices. In biological systems, ET often occurs surprisingly fast over long distances of several tens of angströms. Laser-pulse irradiation is conveniently used to generate strongly oxidizing (or reducing) excited states whose reactions are then studied by time-resolved spectroscopic techniques. While photoluminescence decay and UV-vis absorption supply precise kinetics data, time-resolved infrared absorption (TRIR) and Raman-based spectroscopies have the advantage of providing additional structural information and monitoring vibrational energy flows and dissipation, as well as medium relaxation, that accompany ultrafast ET. We will discuss three cases of photoinduced ET involving the Re(I)(CO)3(N,N) moiety (N,N = polypyridine) that occur much faster than would be expected from ET theories. [Re(4-N-methylpyridinium-pyridine)(CO)3(N,N)](2+) represents a case of excited-state picosecond ET between two different ligands that remains ultrafast even in slow-relaxing solvents, beating the adiabatic limit. This is caused by vibrational/solvational excitation of the precursor state and participation of high-frequency quantum modes in barrier crossing. The case of Re-tryptophan assemblies demonstrates that excited-state Trp → *Re(II) ET is accelerated from nanoseconds to picoseconds when the Re(I)(CO)3(N,N) chromophore is appended to a protein, close to a tryptophan residue. TRIR in combination with DFT calculations and structural studies reveals an interaction between the N,N ligand and the tryptophan indole. It results in partial electronic delocalization in the precursor excited state and likely contributes to the ultrafast ET rate. Long-lived vibrational/solvational excitation of the protein Re(I)(CO)3(N,N)···Trp moiety, documented by dynamic IR band shifts, could be another accelerating factor. The last discussed process, back-ET in a porphyrin-Re(I)(CO)3(N,N) dyad, demonstrates that formation of a hot product accelerates highly exergonic ET in the Marcus inverted region. Overall, it follows that ET can be accelerated by enhancing the electronic interaction and by vibrational excitation of the reacting system and its medium, stressing the importance of quantum nuclear dynamics in ET reactivity. These effects are experimentally accessible by time-resolved vibrational spectroscopies (IR, Raman) in combination with quantum chemical calculations. It is suggested that structural dynamics play different mechanistic roles in light-triggered ET involving electronically excited donors or acceptors than in ground-state processes. While TRIR spectroscopy is well suitable to elucidate ET processes on a molecular-level, transient 2D-IR techniques combining optical and two IR (or terahertz) laser pulses present future opportunities for investigating, driving, and controlling ET.

ATP binding and aspartate protonation enhance photoinduced electron transfer in plant cryptochrome

Cryptochromes are flavoproteins encountered in most vegetal and animal species. They play a role of blue-light receptors in plants and in invertebrates. The putative resting state of the FAD cofactor in these proteins is its fully oxidized form, FADox. Upon blue-light excitation, the isoalloxazine ring (ISO) may undergo an ultrafast reduction by a nearby tryptophan residue W400. This primary reduction triggers a cascade of electron and proton transfers, ultimately leading to the formation of the FADH° radical. A recent experimental study has shown that the yield of FADH° formation in Arabidopsis cryptochrome can be strongly modulated by ATP binding and by pH, affecting the protonation state of D396 (proton donor to FAD°(-)). Here we provide a detailed molecular analysis of these effects by means of combined classical molecular dynamics simulations and time-dependent density functional theory calculations. When ATP is present and D396 protonated, FAD remains in close contact with W400, thereby enhancing electron transfer (ET) from W400 to ISO*. In contrast, deprotonation of D396 and absence of ATP introduce flexibility to the photoactive site prior to FAD excitation, with the consequence of increased ISO-W400 distance and diminished tunneling rate by almost two orders of magnitude. We show that under these conditions, ET from the adenine moiety of FAD becomes a competitive relaxation pathway. Overall, our data suggest that the observed effects of ATP and pH on the FAD photoreduction find their roots in the earliest stage of the photoreduction process; i.e., ATP binding and the protonation state of D396 determine the preferred pathway of ISO* relaxation.

Symmetrized photoinitiated electron flow within the [myoglobin:cytochrome b₅] complex on singlet and triplet time scales: energetics vs dynamics

We report here that photoinitiated electron flow involving a metal-substituted (M = Mg, Zn) myoglobin (Mb) and its physiological partner protein, cytochrome b5 (cyt b5) can be "symmetrized": the [Mb:cyt b5] complex stabilized by three D/E → K mutations on Mb (D44K/D60K/E85K, denoted MMb) exhibits both oxidative and reductive ET quenching of both the singlet and triplet photoexcited MMb states, the direction of flow being determined by the oxidation state of the cyt b5 partner. The first-excited singlet state of MMb ((1)MMb) undergoes ns-time scale reductive ET quenching by Fe(2+)cyt b5 as well as ns-time scale oxidative ET quenching by Fe(3+)cyt b5, both processes involving an ensemble of structures that do not interconvert on this time scale. Despite a large disparity in driving force favoring photooxidation of (1)MMb relative to photoreduction (δ(-ΔG(0)) ≈ 0.4 eV, M = Mg; ≈ 0.2 eV, M = Zn), for each M the average rate constants for the two reactions are the same within error, (1)k(f) > 10(8) s(-1). This surprising observation is explained by considering the driving-force dependence of the Franck-Condon factor in the Marcus equation. The triplet state of the myoglobin ((3)MMb) created by intersystem crossing from (1)MMb likewise undergoes reductive ET quenching by Fe(2+)cyt b5 as well as oxidative ET quenching by Fe(3+)cyt b5. As with singlet ET, the rate constants for oxidative ET quenching and reductive ET quenching on the triplet time scale are the same within error, (3)k(f) ≈ 10(5) s(-1), but here the equivalence is attributable to gating by intracomplex conversion among a conformational ensemble.

Tryptophan-accelerated electron flow across a protein-protein interface

We report a new metallolabeled blue copper protein, Re126W122Cu(I) Pseudomonas aeruginosa azurin, which has three redox sites at well-defined distances in the protein fold: Re(I)(CO)3(4,7-dimethyl-1,10-phenanthroline) covalently bound at H126, a Cu center, and an indole side chain W122 situated between the Re and Cu sites (Re-W122(indole) = 13.1 Å, dmp-W122(indole) = 10.0 Å, Re-Cu = 25.6 Å). Near-UV excitation of the Re chromophore leads to prompt Cu(I) oxidation (<50 ns), followed by slow back ET to regenerate Cu(I) and ground-state Re(I) with biexponential kinetics, 220 ns and 6 μs. From spectroscopic measurements of kinetics and relative ET yields at different concentrations, it is likely that the photoinduced ET reactions occur in protein dimers, (Re126W122Cu(I))2 and that the forward ET is accelerated by intermolecular electron hopping through the interfacial tryptophan: *Re//←W122←Cu(I), where // denotes a protein-protein interface. Solution mass spectrometry confirms a broad oligomer distribution with prevalent monomers and dimers, and the crystal structure of the Cu(II) form shows two Re126W122Cu(II) molecules oriented such that redox cofactors Re(dmp) and W122-indole on different protein molecules are located at the interface at much shorter intermolecular distances (Re-W122(indole) = 6.9 Å, dmp-W122(indole) = 3.5 Å, and Re-Cu = 14.0 Å) than within single protein folds. Whereas forward ET is accelerated by hopping through W122, BET is retarded by a space jump at the interface that lacks specific interactions or water molecules. These findings on interfacial electron hopping in (Re126W122Cu(I))2 shed new light on optimal redox-unit placements required for functional long-range charge separation in protein complexes.

Evolving the [myoglobin, cytochrome b(5)] complex from dynamic toward simple docking: charging the electron transfer reactive patch

We describe photoinitiated electron transfer (ET) from a suite of Zn-substituted myoglobin (Mb) variants to cytochrome b(5) (b(5)). An electrostatic interface redesign strategy has led to the introduction of positive charges into the vicinity of the heme edge through D/E → K charge-reversal mutation combinations at "hot spot" residues (D44, D60, and E85), augmented by the elimination of negative charges from Mb or b(5) by neutralization of heme propionates. These variations create an unprecedentedly large range in the product of the ET partners' total charges (-5 < -q(Mb)q(b(5)) < 40). The binding affinity (K(a)) increases 1000-fold as -q(Mb)q(b(5)) increases through this range and exhibits a surprisingly simple, exponential dependence on -q(Mb)q(b(5)). This is explained in terms of electrostatic interactions between a "charged reactive patch" (crp) on each partner's surface, defined as a compact region around the heme edge that (i) contains the total protein charge of each variant and (ii) encompasses a major fraction of the "reactive region" (Rr) comprising surface atoms with large matrix elements for electron tunneling to the heme. As -q(Mb)q(b(5)) increases, the complex undergoes a transition from fast to slow-exchange dynamics on the triplet ET time scale, with a correlated progression in the rate constants for intracomplex (k(et)) and bimolecular (k(2)) ET. This progression is analyzed by integrating the crp and Rr descriptions of ET into the textbook steady-state treatment of reversible binding between partners that undergo intracomplex ET and found to encompass the full range of behaviors predicted by the model. The generality of this approach is demonstrated by its application to the extensive body of data for the ET complex between the photosynthetic reaction center and cytochrome c(2). Deviations from this model also are discussed.

Temperature invariance of the nitrogenase electron transfer mechanism

Earlier studies of electron transfer (ET) from the nitrogenase Fe protein to the MoFe protein concluded that the mechanism for ET changed during cooling from 25 to 5 °C, based on the observation that the rate constant for Fe protein to MoFe protein ET decreases strongly, with a nonlinear Arrhenius plot. They further indicated that the ET was reversible, with complete ET at ambient temperature but with an equilibrium constant near unity at 5 °C. These studies were conducted with buffers having a strong temperature coefficient. We have examined the temperature variation in the kinetics of oxidation of the Fe protein by the MoFe protein at a constant pH of 7.4 fixed by the buffer 3-(N-morpholino)propanesulfonic acid (MOPS), which has a very small temperature coefficient. Using MOPS, we also observe temperature-dependent ET rate constants, with nonlinear Arrhenius plots, but we find that ET is gated across the temperature range by a conformational change that involves the binding of numerous water molecules, consistent with an unchanging ET mechanism. Furthermore, there is no solvent kinetic isotope effect throughout the temperature range studied, again consistent with an unchanging mechanism. In addition, the nonlinear Arrhenius plots are explained by the change in heat capacity caused by the binding of waters in an invariant gating ET mechanism. Together, these observations contradict the idea of a change in ET mechanism with cooling. Finally, the extent of ET at constant pH does not change significantly with temperature, in contrast to the previously proposed change in ET equilibrium.