Differential influence of dynamic processes on forward and reverse electron transfer across a protein-protein interface

Study on enhancement of hemoglobin antitoxic ability modified with chromium and ruthenium

Hemoglobin is essential for carrying oxygen (O₂) in the blood. However, its ability to bind excessively to carbon monoxide (CO) makes it susceptible to CO poisoning. To reduce the risk of CO poisoning, Cr-based heme and Ru-based heme were selected from among many transition metal-based hemes based on their characteristics of adsorption conformation, binding intensity, spin multiplicity, and electronic properties. The results showed that hemoglobin modified by Cr-based heme and Ru-based heme had strong anti-CO poisoning abilities. The Cr-based heme and Ru-based heme exhibited much stronger affinity for O₂ (-190.67 kJ/mol and -143.18 kJ/mol, respectively) than Fe-based heme (-44.60 kJ/mol). Moreover, Cr-based heme and Ru-based heme exhibited much weaker affinity for CO (-121.50 kJ/mol and -120.88 kJ/mol, respectively) than their affinity for O₂, suggesting that they were less likely to cause CO poisoning. The electronic structure analysis also supported this conclusion. Additionally, molecular dynamics analysis showed that hemoglobin modified by Cr-based heme and Ru-based heme was stable. Our findings offer a novel and effective strategy for enhancing the reconstructed hemoglobin's ability to bind O₂ and reduce its potential for CO poisoning.

Photoinduced hole hopping through tryptophans in proteins

Hole hopping through tryptophan/tyrosine chains enables rapid unidirectional charge transport over long distances. We have elucidated structural and dynamical factors controlling hopping speed and efficiency in two modified azurin constructs that include a rhenium(I) sensitizer, Re(His)(CO)₃(dmp)⁺, and one or two tryptophans (W₁, W₂). Experimental kinetics investigations showed that the two closely spaced (3 to 4 Å) intervening tryptophans dramatically accelerated long-range electron transfer (ET) from Cu^I to the photoexcited sensitizer. In our theoretical work, we found that time-dependent density-functional theory (TDDFT) quantum mechanics/molecular mechanics/molecular dynamics (QM/MM/MD) trajectories of low-lying triplet excited states of Re^I(His)(CO)₃(dmp)⁺-W₁(-W₂) exhibited crossings between sensitizer-localized (*Re) and charge-separated [Re^I(His)(CO)₃(dmp^•-)/(W₁ ^•+ or W₂ ^•+)] (CS1 or CS2) states. Our analysis revealed that the distances, angles, and mutual orientations of ET-active cofactors fluctuate in a relatively narrow range in which the cofactors are strongly coupled, enabling adiabatic ET. Water-dominated electrostatic field fluctuations bring *Re and CS1 states to a crossing where *Re(CO)₃(dmp)⁺←W₁ ET occurs, and CS1 becomes the lowest triplet state. ET is promoted by solvation dynamics around *Re(CO)₃(dmp)⁺(W₁); and CS1 is stabilized by Re(dmp^•-)/W₁ ^•+ electron/hole interaction and enhanced W₁ ^•+ solvation. The second hop, W₁ ^•+←W₂, is facilitated by water fluctuations near the W₁/W₂ unit, taking place when the electrostatic potential at W₂ drops well below that at W₁ ^•+ Insufficient solvation and reorganization around W₂ make W₁ ^•+←W₂ ET endergonic, shifting the equilibrium toward W₁ ^•+ and decreasing the charge-separation yield. We suggest that multiscale TDDFT/MM/MD is a suitable technique to model the simultaneous evolution of photogenerated excited-state manifolds.

Hole Hopping Across a Protein-Protein Interface

We have investigated photoinduced hole hopping in a Pseudomonas aeruginosa azurin mutant Re126WWCuI, where two adjacent tryptophan residues (W124 and W122) are inserted between the CuI center and a Re photosensitizer coordinated to a H126 imidazole (Re = ReI(H126)(CO)3(dmp)+, dmp = 4,7-dimethyl-1,10-phenanthroline). Optical excitation of this mutant in aqueous media (≤40 μM) triggers 70 ns electron transport over 23 Å, yielding a long-lived (120 μs) ReI(H126)(CO)3(dmp•-)WWCuII product. The Re126FWCuI mutant (F124, W122) is not redox-active under these conditions. Upon increasing the concentration to 0.2-2 mM, {Re126WWCuI}2 and {Re126FWCuI}2 are formed with the dmp ligand of the Re photooxidant of one molecule in close contact (3.8 Å) with the W122' indole on the neighboring chain. In addition, {Re126WWCuI}2 contains an interfacial tryptophan quadruplex of four indoles (3.3-3.7 Å apart). In both mutants, dimerization opens an intermolecular W122' → //*Re ET channel (// denotes the protein interface, *Re is the optically excited sensitizer). Excited-state relaxation and ET occur together in two steps (time constants of ∼600 ps and ∼8 ns) that lead to a charge-separated state containing a Re(H126)(CO)3(dmp•-)//(W122•+)' unit; then (CuI)' is oxidized intramolecularly (60-90 ns) by (W122•+)', forming ReI(H126)(CO)3(dmp•-)WWCuI//(CuII)'. The photocycle is closed by ∼1.6 μs ReI(H126)(CO)3(dmp•-) → //(CuII)' back ET that occurs over 12 Å, in contrast to the 23 Å, 120 μs step in Re126WWCuI. Importantly, dimerization makes Re126FWCuI photoreactive and, as in the case of {Re126WWCuI}2, channels the photoproduced "hole" to the molecule that was not initially photoexcited, thereby shortening the lifetime of ReI(H126)(CO)3(dmp•-)//CuII. Although two adjacent W124 and W122 indoles dramatically enhance CuI → *Re intramolecular multistep ET, the tryptophan quadruplex in {Re126WWCuI}2 does not accelerate intermolecular electron transport; instead, it acts as a hole storage and crossover unit between inter- and intramolecular ET pathways. Irradiation of {Re126WWCuII}2 or {Re126FWCuII}2 also triggers intermolecular W122' → //*Re ET, and the Re(H126)(CO)3(dmp•-)//(W122•+)' charge-separated state decays to the ground state by ∼50 ns ReI(H126)(CO)3(dmp•-)+ → //(W122•+)' intermolecular charge recombination. Our findings shed light on the factors that control interfacial hole/electron hopping in protein complexes and on the role of aromatic amino acids in accelerating long-range electron transport.

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.

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.

Dissecting regulation mechanism of the FMN to heme interdomain electron transfer in nitric oxide synthases

Nitric oxide synthase (NOS), a flavo-hemoprotein, is responsible for biosynthesis of nitric oxide (NO) in mammals. Three NOS isoforms, iNOS, eNOS and nNOS (inducible, endothelial, and neuronal NOS), achieve their biological functions by tight control of interdomain electron transfer (IET) process through interdomain interactions. In particular, the FMN-heme IET is essential in coupling electron transfer in the reductase domain with NO synthesis in the heme domain by delivery of electrons required for O2 activation at the catalytic heme site. Emerging evidence indicates that calmodulin (CaM) activates NO synthesis in eNOS and nNOS by a conformational change of the FMN domain from its shielded electron-accepting (input) state to a new electron-donating (output) state, and that CaM is also required for proper alignment of the FMN and heme domains in the three NOS isoforms. In the absence of a structure of full-length NOS, an integrated approach of spectroscopic, rapid kinetic and mutagenesis methods is required to unravel regulation mechanism of the FMN-heme IET process. This is to investigate the roles of the FMN domain motions and the docking between the primary functional FMN and heme domains in regulating NOS activity. The recent developments in this area that are driven by the combined approach are the focuses of this review. A better understanding of the roles of interdomain FMN/heme interactions and CaM binding may serve as a basis for the rational design of new selective modulators of the NOS enzymes.

Distance-independent charge recombination kinetics in cytochrome c-cytochrome c peroxidase complexes: compensating changes in the electronic coupling and reorganization energies

Charge recombination rate constants vary no more than 3-fold for interprotein ET in the Zn-substituted wild type (WT) cytochrome c peroxidase (CcP):cytochrome c (Cc) complex and in complexes with four mutants of the Cc protein (i.e., F82S, F82W, F82Y, and F82I), despite large differences in the ET distance. Theoretical analysis indicates that charge recombination for all complexes involves a combination of tunneling and hopping via Trp191. For three of the five structures (WT and F82S(W)), the protein favors hopping more than that in the other two structures that have longer heme → ZnP distances (F82Y(I)). Experimentally observed biexponential ET kinetics is explained by the complex locking in alternative coupling pathways, where the acceptor hole state is either primarily localized on ZnP (slow phase) or on Trp191 (fast phase). The large conformational differences between the CcP:Cc interface for the F82Y(I) mutants compared to that the WT and F82S(W) complexes are predicted to change the reorganization energies for the CcP:Cc ET reactions because of changes in solvent exposure and interprotein ET distances. Since the recombination reaction is likely to occur in the inverted Marcus regime, an increased reorganization energy compensates the decreased role for hopping recombination (and the longer transfer distance) in the F82Y(I) mutants. Taken together, coupling pathway and reorganization energy effects for the five protein complexes explain the observed insensitivity of recombination kinetics to donor-acceptor distance and docking pose and also reveals how hopping through aromatic residues can accelerate long-range ET.

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

Formation of a transient [myoglobin (Mb), cytochrome b(5) (cyt b(5))] complex is required for the reductive repair of inactive ferri-Mb to its functional ferro-Mb state. The [Mb, cyt b(5)] complex exhibits dynamic docking (DD), with its cyt b(5) partner in rapid exchange at multiple sites on the Mb surface. A triple mutant (Mb(3M)) was designed as part of efforts to shift the electron-transfer process to the simple docking (SD) regime, in which reactive binding occurs at a restricted, reactive region on the Mb surface that dominates the docked ensemble. An electrostatically-guided brownian dynamics (BD) docking protocol was used to generate an initial ensemble of reactive configurations of the complex between unrelaxed partners. This ensemble samples a broad and diverse array of heme-heme distances and orientations. These configurations seeded all-atom constrained molecular dynamics simulations (MD) to generate relaxed complexes for the calculation of electron tunneling matrix elements (T(DA)) through tunneling-pathway analysis. This procedure for generating an ensemble of relaxed complexes combines the ability of BD calculations to sample the large variety of available conformations and interprotein distances, with the ability of MD to generate the atomic level information, especially regarding the structure of water molecules at the protein-protein interface, that defines electron-tunneling pathways. We used the calculated T(DA) values to compute ET rates for the [Mb(wt), cyt b(5)] complex and for the complex with a mutant that has a binding free energy strengthened by three D/E → K charge-reversal mutations, [Mb(3M), cyt b(5)]. The calculated rate constants are in agreement with the measured values, and the mutant complex ensemble has many more geometries with higher T(DA) values than does the wild-type Mb complex. Interestingly, water plays a double role in this electron-transfer system, lowering the tunneling barrier as well as inducing protein interface remodeling that screens the repulsion between the negatively-charged propionates of the two hemes.

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.

The reversible opening of water channels in cytochrome c modulates the heme iron reduction potential

Dynamic protein-solvent interactions are fundamental for life processes, but their investigation is still experimentally very demanding. Molecular dynamics simulations up to hundreds of nanoseconds can bring to light unexpected events even for extensively studied biomolecules. This paper reports a combined computational/experimental approach that reveals the reversible opening of two distinct fluctuating cavities in Saccharomyces cerevisiae iso-1-cytochrome c. Both channels allow water access to the heme center. By means of a mixed quantum mechanics/molecular dynamics (QM/MD) theoretical approach, the perturbed matrix method (PMM), that allows to reach long simulation times, changes in the reduction potential of the heme Fe(3+)/Fe(2+) couple induced by the opening of each cavity are calculated. Shifts of the reduction potential upon changes in the hydration of the heme propionates are observed. These variations are relatively small but significant and could therefore represent a tool developed by cytochrome c for the solvent driven, fine-tuning of its redox functionality.

Inter-monomer electron transfer is too slow to compete with monomeric turnover in bc(1) complex

The homodimeric bc(1) complexes are membrane proteins essential in respiration and photosynthesis. The ~11Å distance between the two b(L)-hemes of the dimer opens the possibility of electron transfer between them, but contradictory reports make such inter-monomer electron transfer controversial. We have constructed in Rhodobacter sphaeroides a heterodimeric expression system similar to those used before, in which the bc(1) complex can be mutated differentially in the two copies of cyt b to test for inter-monomer electron transfer, but found that genetic recombination by cross-over then occurs to produce wild-type homodimer. Selection pressure under photosynthetic growth always favored the homodimer over heterodimeric variants enforcing inter-monomer electron transfer, showing that the latter are not competitive. These results, together with kinetic analysis of myxothiazol titrations, demonstrate that inter-monomer electron transfer does not occur at rates competitive with monomeric turnover. We examine the results from other groups interpreted as demonstrating rapid inter-monomer electron transfer, conclude that similar mechanisms are likely to be in play, and suggest that such claims might need to be re-examined.

Role of an isoform-specific serine residue in FMN-heme electron transfer in inducible nitric oxide synthase

In the crystal structure of a calmodulin (CaM)-bound FMN domain of human inducible nitric oxide synthase (NOS), the CaM-binding region together with CaM forms a hinge, and pivots on an R536(NOS)/E47(CaM) pair (Xia et al. J Biol Chem 284:30708-30717, 2009). Notably, isoform-specific human inducible NOS S562 and C563 residues form hydrogen bonds with the R536 residue through their backbone oxygens. In this study, we investigated the roles of the S562 and C563 residues in the NOS FMN-heme interdomain electron transfer (IET), the rates of which can be used to probe the interdomain FMN/heme alignment. Human inducible NOS S562K and C563R mutants of an oxygenase/FMN (oxyFMN) construct were made by introducing charged residues at these sites as found in human neuronal NOS and endothelial NOS isoforms, respectively. The IET rate constant of the S562K mutant is notably decreased by one third, and its flavin fluorescence intensity per micromole per liter is diminished by approximately 24 %. These results suggest that a positive charge at position 562 destabilizes the hydrogen-bond-mediated NOS/CaM alignment, resulting in slower FMN-heme IET in the mutant. On the other hand, the IET rate constant of the C563R mutant is similar to that of the wild-type, indicating that the mutational effect is site-specific. Moreover, the human inducible NOS oxyFMN R536E mutant was constructed to disrupt the bridging CaM/NOS interaction, and its FMN-heme IET rate was decreased by 96 %. These results demonstrated a new role of the isoform-specific serine residue of the key CaM/FMN(NOS) bridging site in regulating the FMN-heme IET (possibly by tuning the alignment of the FMN and heme domains).

Mechanism of Nitric Oxide Synthase Regulation: Electron Transfer and Interdomain Interactions

Nitric oxide synthase (NOS), a flavo-hemoprotein, tightly regulates nitric oxide (NO) synthesis and thereby its dual biological activities as a key signaling molecule for vasodilatation and neurotransmission at low concentrations, and also as a defensive cytotoxin at higher concentrations. Three NOS isoforms, iNOS, eNOS and nNOS (inducible, endothelial, and neuronal NOS), achieve their key biological functions by tight regulation of interdomain electron transfer (IET) process via interdomain interactions. In particular, the FMN-heme IET is essential in coupling electron transfer in the reductase domain with NO synthesis in the heme domain by delivery of electrons required for O(2) activation at the catalytic heme site. Compelling evidence indicates that calmodulin (CaM) activates NO synthesis in eNOS and nNOS through a conformational change of the FMN domain from its shielded electron-accepting (input) state to a new electron-donating (output) state, and that CaM is also required for proper alignment of the domains. Another exciting recent development in NOS enzymology is the discovery of importance of the the FMN domain motions in modulating reactivity and structure of the catalytic heme active site (in addition to the primary role of controlling the IET processes). In the absence of a structure of full-length NOS, an integrated approach of spectroscopic (e.g. pulsed EPR, MCD, resonance Raman), rapid kinetics (laser flash photolysis and stopped flow) and mutagenesis methods is critical to unravel the molecular details of the interdomain FMN/heme interactions. This is to investigate the roles of dynamic conformational changes of the FMN domain and the docking between the primary functional FMN and heme domains in regulating NOS activity. The recent developments in understanding of mechanisms of the NOS regulation that are driven by the combined approach are the focuses of this review. An improved understanding of the role of interdomain FMN/heme interaction and CaM binding may serve as the basis for the design of new selective inhibitors of NOS isoforms.

Intraprotein electron transfer between the FMN and heme domains in endothelial nitric oxide synthase holoenzyme

Intraprotein electron transfer (IET) from flavin mononucleotide (FMN) to heme is an essential step in nitric oxide (NO) synthesis by NO synthase (NOS). The IET kinetics in neuronal and inducible NOS (nNOS and iNOS) holoenzymes have been previously determined in our laboratories by laser flash photolysis [reviewed in: C.J. Feng, G. Tollin, Dalton Trans., (2009) 6692-6700]. Here we report the kinetics of the IET in a bovine endothelial NOS (eNOS) holoenzyme in the presence and absence of added calmodulin (CaM). The IET rate constant in the presence of CaM is estimated to be ~4.3s(-1). No IET was observed in the absence of CaM, indicating that CaM is the primary factor in controlling the FMN-heme IET in the eNOS enzyme. The IET rate constant value for the eNOS holoenzyme is approximately 10 times smaller than those obtained for the iNOS and CaM-bound nNOS holoenzymes. Possible mechanisms underlying the difference in IET kinetics among the NOS isoforms are discussed. Because the rate-limiting step in the IET process in these enzymes is the conformational change from input state to output state, a slower conformational change (than in the other isoforms) is most likely to cause the slower IET in eNOS.

Photoinitiated singlet and triplet electron transfer across a redesigned [myoglobin, cytochrome b5] interface

We describe a strategy by which reactive binding of a weakly bound, 'dynamically docked (DD)' complex without a known structure can be strengthened electrostatically through optimized placement of surface charges, and discuss its use in modulating complex formation between myoglobin (Mb) and cytochrome b(5) (b(5)). The strategy employs paired Brownian dynamics (BD) simulations, one which monitors overall binding, the other reactive binding, to examine [X --> K] mutations on the surface of the partners, with a focus on single and multiple [D/E --> K] charge reversal mutations. This procedure has been applied to the [Mb, b(5)] complex, indicating mutations of Mb residues D44, D60, and E85 to be the most promising, with combinations of these showing a nonlinear enhancement of reactive binding. A novel method of displaying BD profiles shows that the 'hits' of b(5) on the surfaces of Mb(WT), Mb(D44K/D60K), and Mb(D44K/D60K/E85K) progressively coalesce into two 'clusters': a 'diffuse' cluster of hits that are distributed over the Mb surface and have negligible electrostatic binding energy and a 'reactive' cluster of hits with considerable stability that are localized near its heme edge, with short Fe-Fe distances favorable to electron transfer (ET). Thus, binding and reactivity progressively become correlated by the mutations. This finding relates to recent proposals that complex formation is a two-step process, proceeding through the formation of a weakly bound encounter complex to a well-defined bound complex. The design procedure has been tested through measurements of photoinitiated ET between the Zn-substituted forms of Mb(WT), Mb(D44K/D60K), and Mb(D44K/D60K/E85K) and Fe(3+)b(5). Both mutants convert the complex from the DD regime exhibited by Mb(WT), in which the transient complex is in fast kinetic exchange with its partners, k(off) >> k(et), to the slow-exchange regime, k(et) >> k(off), and both mutants exhibit rapid intracomplex ET from the triplet excited state to Fe(3+)b(5) (rate constant, k(et) approximately 10(6) s(-1)). The affinity constants of the mutant Mbs cannot be derived through conventional analysis procedures because intracomplex singlet ET quenching causes the triplet-ground absorbance difference to progressively decrease during a titration, but this effect has been incorporated into a new procedure for computing binding constants. Most importantly, these measurements reveal the presence of fast photoinduced singlet ET across the protein-protein interface, (1)k(et) approximately 2 x 10(8) s(-1).

Computational design and elaboration of a de novo heterotetrameric alpha-helical protein that selectively binds an emissive abiological (porphinato)zinc chromophore

The first example of a computationally de novo designed protein that binds an emissive abiological chromophore is presented, in which a sophisticated level of cofactor discrimination is pre-engineered. This heterotetrameric, C(2)-symmetric bundle, A(His):B(Thr), uniquely binds (5,15-di[(4-carboxymethyleneoxy)phenyl]porphinato)zinc [(DPP)Zn] via histidine coordination and complementary noncovalent interactions. The A(2)B(2) heterotetrameric protein reflects ligand-directed elements of both positive and negative design, including hydrogen bonds to second-shell ligands. Experimental support for the appropriate formulation of [(DPP)Zn:A(His):B(Thr)](2) is provided by UV/visible and circular dichroism spectroscopies, size exclusion chromatography, and analytical ultracentrifugation. Time-resolved transient absorption and fluorescence spectroscopic data reveal classic excited-state singlet and triplet PZn photophysics for the A(His):B(Thr):(DPP)Zn protein (k(fluorescence) = 4 x 10(8) s(-1); tau(triplet) = 5 ms). The A(2)B(2) apoprotein has immeasurably low binding affinities for related [porphinato]metal chromophores that include a (DPP)Fe(III) cofactor and the zinc metal ion hemin derivative [(PPIX)Zn], underscoring the exquisite active-site binding discrimination realized in this computationally designed protein. Importantly, elements of design in the A(His):B(Thr) protein ensure that interactions within the tetra-alpha-helical bundle are such that only the heterotetramer is stable in solution; corresponding homomeric bundles present unfavorable ligand-binding environments and thus preclude protein structural rearrangements that could lead to binding of (porphinato)iron cofactors.

Visualization of the encounter ensemble of the transient electron transfer complex of cytochrome c and cytochrome c peroxidase

Recent studies have provided experimental evidence for the existence of an encounter complex, a transient intermediate in the formation of protein complexes. We use paramagnetic relaxation enhancement NMR spectroscopy in combination with Monte Carlo simulations to characterize and visualize the ensemble of encounter orientations in the short-lived electron transfer complex of yeast cytochrome c (Cc) and cytochrome c peroxidase (CcP). The complete conformational space sampled by the protein molecules during the dynamic part of the interaction was mapped experimentally. The encounter complex was described by an electrostatic ensemble of orientations based on Monte Carlo calculations, considering the protein structures in atomic detail. We demonstrate that this visualization of the encounter complex, in combination with the specific complex, is in excellent agreement with the experimental data. Our results indicate that Cc samples only about 15% of the surface area of CcP, surrounding the specific binding interface. The encounter complex is populated for 30% of the time, representing a mere 0.5 kcal/mol difference in the free energies between the two states. This delicate balance is interpreted to be a consequence of the conflicting requirements of fast electron transfer and high turnover of the complex.

Molecular basis for directional electron transfer

Biological macromolecules involved in electron transfer reactions display chains of closely packed redox cofactors when long distances must be bridged. This is a consequence of the need to maintain a rate of transfer compatible with metabolic activity in the framework of the exponential decay of electron tunneling with distance. In this work intermolecular electron transfer was studied in kinetic experiments performed with the small tetraheme cytochrome from Shewanella oneidensis MR-1 and from Shewanella frigidimarina NCIMB400 using non-physiological redox partners. This choice allowed the effect of specific recognition and docking to be eliminated from the measured rates. The results were analyzed with a kinetic model that uses the extensive thermodynamic characterization of these proteins reported in the literature to discriminate the kinetic contribution of each heme to the overall rate of electron transfer. This analysis shows that, in this redox chain that spans 23 A, the kinetic properties of the individual hemes establish a functional specificity for each redox center. This functional specificity combined with the thermodynamic properties of these soluble proteins ensures directional electron flow within the cytochrome even outside of the context of a functioning respiratory chain.

Electrostatic redesign of the [myoglobin, cytochrome b5] interface to create a well-defined docked complex with rapid interprotein electron transfer

Cyt b(5) is the electron-carrier "repair" protein that reduces met-Mb and met-Hb to their O(2)-carrying ferroheme forms. Studies of electron transfer (ET) between Mb and cyt b(5) revealed that they react on a "Dynamic Docking" (DD) energy landscape on which binding and reactivity are uncoupled: binding is weak and involves an ensemble of nearly isoenergetic configurations, only a few of which are reactive; those few contribute negligibly to binding. We set the task of redesigning the surface of Mb so that its reaction with cyt b(5) instead would occur on a conventional "simple docking" (SD) energy landscape, on which a complex exhibits a well-defined (set of) reactive binding configuration(s), with binding and reactivity thus no longer being decoupled. We prepared a myoglobin (Mb) triple mutant (D44K/D60K/E85K; Mb(+6)) substituted with Zn-deuteroporphyrin and monitored cytochrome b(5) (cyt b(5)) binding and electron transfer (ET) quenching of the (3)ZnMb(+6) triplet state. In contrast, to Mb(WT), the three charge reversals around the "front-face" heme edge of Mb(+6) have directed cyt b(5) to a surface area of Mb adjacent to its heme, created a well-defined, most-stable structure that supports good ET pathways, and apparently coupled binding and ET: both K(a) and k(et) are increased by the same factor of approximately 2 x 10(2), creating a complex that exhibits a large ET rate constant, k(et) = 10(6 1) s(-1), and is in slow exchange (k(off) << k(et)). In short, these mutations indeed appear to have created the sought-for conversion from DD to simple docking (SD) energy landscapes.

Kinetic-dynamic model for conformational control of an electron transfer photocycle: mixed-metal hemoglobin hybrids

It is becoming increasingly clear that the transfer of an electron across a protein-protein interface is coupled to the dynamics of conformational conversion between and within ensembles of interface conformations. Electron transfer (ET) reactions in conformationally mobile systems provide a "clock" against which the rapidity of a dynamic process may be measured, and we here report a simple kinetic (master equation) model that self-consistently incorporates conformational dynamics into an ET photocycle comprised of a photoinitiated "forward" step and thermal return to ground. This kinetic/dynamic (KD) model assumes an ET complex exists as multiple interconverting conformations which partition into an ET-optimized (reactive; R) population and a less-reactive population ( S). We take the members of each population to be equivalent by constraining them to have the same conformational energy, the same average rate constant for conversion to members of the other population, and the same rate constants for forward and back ET. The result is a mapping of a complicated energy surface onto the simple "gating", two-well surface, but with rate constants that are defined microscopically. This model successfully describes the changes in the ET photocycle within the "predocked" mixed-metal hemoglobin (Hb) hybrid, [alpha(Zn), beta(Fe3+N 3 (-))], as conformational kinetics are modulated by variations in viscosity (eta = 1-15 cP; 20 degrees C). The description reveals how the conformational "routes" by which a hybrid progresses through a photocycle differ in different dynamic regimes. Even at eta = 1 cP, the populations are not in fast exchange, and ET involves a complex interplay between conformational and ET processes; at intermediate viscosities the hybrid exhibits "differential dynamics" in which the forward and back ET processes involve different initial ensembles of configurational substates; by eta = 15 cP, the slow-exchange limit is approached. Even at low viscosity, the ET-coupled motions are fairly slow, with rate constants of <10 (3) s (-1). Current ideas about Hb function lead to the testable hypothesis that ET in the hybrid may be coupled to allosteric fluctuations of the two [alpha 1, beta 2] dimers of Hb.

The Q-cycle reviewed: How well does a monomeric mechanism of the bc(1) complex account for the function of a dimeric complex?

Recent progress in understanding the Q-cycle mechanism of the bc(1) complex is reviewed. The data strongly support a mechanism in which the Q(o)-site operates through a reaction in which the first electron transfer from ubiquinol to the oxidized iron-sulfur protein is the rate-determining step for the overall process. The reaction involves a proton-coupled electron transfer down a hydrogen bond between the ubiquinol and a histidine ligand of the [2Fe-2S] cluster, in which the unfavorable protonic configuration contributes a substantial part of the activation barrier. The reaction is endergonic, and the products are an unstable ubisemiquinone at the Q(o)-site, and the reduced iron-sulfur protein, the extrinsic mobile domain of which is now free to dissociate and move away from the site to deliver an electron to cyt c(1) and liberate the H(+). When oxidation of the semiquinone is prevented, it participates in bypass reactions, including superoxide generation if O(2) is available. When the b-heme chain is available as an acceptor, the semiquinone is oxidized in a process in which the proton is passed to the glutamate of the conserved -PEWY- sequence, and the semiquinone anion passes its electron to heme b(L) to form the product ubiquinone. The rate is rapid compared to the limiting reaction, and would require movement of the semiquinone closer to heme b(L) to enhance the rate constant. The acceptor reactions at the Q(i)-site are still controversial, but likely involve a "two-electron gate" in which a stable semiquinone stores an electron. Possible mechanisms to explain the cyt b(150) phenomenon are discussed, and the information from pulsed-EPR studies about the structure of the intermediate state is reviewed. The mechanism discussed is applicable to a monomeric bc(1) complex. We discuss evidence in the literature that has been interpreted as shown that the dimeric structure participates in a more complicated mechanism involving electron transfer across the dimer interface. We show from myxothiazol titrations and mutational analysis of Tyr-199, which is at the interface between monomers, that no such inter-monomer electron transfer is detected at the level of the b(L) hemes. We show from analysis of strains with mutations at Asn-221 that there are coulombic interactions between the b-hemes in a monomer. The data can also be interpreted as showing similar coulombic interaction across the dimer interface, and we discuss mechanistic implications.

Comparison of intra- vs intermolecular long-range electron transfer in crystals of ruthenium-modified azurin

Selective metal-ion incorporation and ligand substitution are employed to control whether electrons tunnel over intra- or intermolecular separations in crystals of P. aeruginosa azurin modified with Ru-polypyridine complexes. Cu(1+)-to-Ru3+ electron transfer (ET) across a specific protein-protein interface in the crystal lattice has a time constant 5-10 times longer than ET between the same donor and acceptor within a single protein (tauET = 5 vs 0.5-1.0 micros). Slower intermolecular ET agrees well with a longer distance between redox centers across the inter-protein (18.9 A) compared to the intra-protein separation (17.0 A) and indicates that the closest donor/acceptor pair dominates crystal ET. Lowering the crystal pH accelerates inter-protein ET (tauET = 1.0 micros) but not intra-protein ET. Faster inter-protein ET likely results from a pH-induced peptide bond flip that perturbs hydrogen bonding in the path between Ru and Cu centers on adjacent molecules.

Distal charge transport in peptides

Biological systems often transport charges and reactive processes over substantial distances. Traditional models of chemical kinetics generally do not describe such extreme distal processes. In this Review, an atomistic model for a distal transport of information, which was specifically developed for peptides, is considered. Chemical reactivity is taken as the result of distal effects based on two-step bifunctional kinetics involving unique, very rapid motional properties of peptides in the subpicosecond regime. The bifunctional model suggests highly efficient transport of charge and reactivity in an isolated peptide over a substantial distance; conversely, a very low efficiency in a water environment was found. The model suggests ultrafast transport of charge and reactivity over substantial molecular distances in a peptide environment. Many such domains can be active in a protein.

Sulfite oxidizing enzymes

Sulfite oxidizing enzymes are essential mononuclear molybdenum (Mo) proteins involved in sulfur metabolism of animals, plants and bacteria. There are three such enzymes presently known: (1) sulfite oxidase (SO) in animals, (2) SO in plants, and (3) sulfite dehydrogenase (SDH) in bacteria. X-ray crystal structures of enzymes from all three sources (chicken SO, Arabidopsis thaliana SO, and Starkeya novella SDH) show nearly identical square pyramidal coordination around the Mo atom, even though the overall structures of the proteins and the presence of additional cofactors vary. This structural information provides a molecular basis for studying the role of specific amino acids in catalysis. Animal SO catalyzes the final step in the degradation of sulfur-containing amino acids and is critical in detoxifying excess sulfite. Human SO deficiency is a fatal genetic disorder that leads to early death, and impaired SO activity is implicated in sulfite neurotoxicity. Animal SO and bacterial SDH contain both Mo and heme domains, whereas plant SO only has the Mo domain. Intraprotein electron transfer (IET) between the Mo and Fe centers in animal SO and bacterial SDH is a key step in the catalysis, which can be studied by laser flash photolysis in the presence of deazariboflavin. IET studies on animal SO and bacterial SDH clearly demonstrate the similarities and differences between these two types of sulfite oxidizing enzymes. Conformational change is involved in the IET of animal SO, in which electrostatic interactions may play a major role in guiding the docking of the heme domain to the Mo domain prior to electron transfer. In contrast, IET measurements for SDH demonstrate that IET occurs directly through the protein medium, which is distinctly different from that in animal SO. Point mutations in human SO can result in significantly impaired IET or no IET, thus rationalizing their fatal effects. The recent developments in our understanding of sulfite oxidizing enzyme mechanisms that are driven by a combination of molecular biology, rapid kinetics, pulsed electron paramagnetic resonance (EPR), and computational techniques are the subject of this review.

Coupling coherence distinguishes structure sensitivity in protein electron transfer

Quantum mechanical analysis of electron tunneling in nine thermally fluctuating cytochrome b562 derivatives reveals two distinct protein-mediated coupling limits. A structure-insensitive regime arises for redox partners coupled through dynamically averaged multiple-coupling pathways (in seven of the nine derivatives) where heme-edge coupling leads to the multiple-pathway regime. A structure-dependent limit governs redox partners coupled through a dominant pathway (in two of the nine derivatives) where axial-ligand coupling generates the single-pathway limit and slower rates. This two-regime paradigm provides a unified description of electron transfer rates in 26 ruthenium-modified heme and blue-copper proteins, as well as in numerous photosynthetic proteins.

Role of protein frame and solvent for the redox properties of azurin from Pseudomonas aeruginosa

We have coupled hybrid quantum mechanics (density functional theory; Car-Parrinello)/molecular mechanics molecular dynamics simulations to a grand-canonical scheme, to calculate the in situ redox potential of the Cu(2+) + e(-) --> Cu(+) half reaction in azurin from Pseudomonas aeruginosa. An accurate description at atomistic level of the environment surrounding the metal-binding site and finite-temperature fluctuations of the protein structure are both essential for a correct quantitative description of the electronic properties of this system. We report a redox potential shift with respect to copper in water of 0.2 eV (experimental 0.16 eV) and a reorganization free energy lambda = 0.76 eV (experimental 0.6-0.8 eV). The electrostatic field of the protein plays a crucial role in fine tuning the redox potential and determining the structure of the solvent. The inner-sphere contribution to the reorganization energy is negligible. The overall small value is mainly due to solvent rearrangement at the protein surface.

Solution structure and dynamics of the complex between cytochrome c and cytochrome c peroxidase determined by paramagnetic NMR

The physiological complex of yeast cytochrome c peroxidase and iso-1-cytochrome c is a paradigm for biological electron transfer. Using paramagnetic NMR spectroscopy, we have determined the conformation of the protein complex in solution, which is shown to be very similar to that observed in the crystal structure [Pelletier H, Kraut J (1992) Science 258:1748-1755]. Our results support the view that this transient electron transfer complex is dynamic. The solution structure represents the dominant protein-protein orientation, which, according to our estimates, is occupied for >70% of the lifetime of the complex, with the rest of the time spent in the dynamic encounter state. Based on the observed paramagnetic effects, we have delineated the conformational space sampled by the protein molecules during the dynamic part of the interaction, providing experimental support for the theoretical predictions of the classical Brownian dynamics study [Northrup SH, Boles JO, Reynolds JCL (1988) Science 241:67-70]. Our findings corroborate the dynamic behavior of this complex and offer an insight into the mechanism of the protein complex formation in solution.

Sugar-derived glasses support thermal and photo-initiated electron transfer processes over macroscopic distances

Trehalose-derived glasses are shown to support long range electron transfer reactions between spatially well separated donors and protein acceptors. The results indicate that these matrices can be used not only to greatly stabilize protein structures but also to facilitate both thermal and photo-initiated hemeprotein reduction over large macroscopic distances. To date the promise of exciting new protein-based technologies that can harness the exceptional tunability of protein functionality has been significantly thwarted by both intrinsic instability and stringent solvent/environment requirements for the expression of functional properties. The presented results raise the prospect of overcoming these limitations with respect to incorporating redox active proteins into solid state devices such as tunable batteries, switches, and solar cells. The findings also have implications for formulations intended to enhance long term storage of biomaterials, new protein-based synthetic strategies, and biophysical studies of functional intermediates trapped under nonequilibrium conditions. In addition, the study shows that certain sugars such as glucose or tagatose, when added to redox-inactive glassy matrices, can be used as a source of thermal electrons that can be harvested by suitable redox active proteins, raising the prospect of using common sugars as an electron source in solid state thermal fuel cells.

Electron-transfer mechanism of the triplet state quenching of aluminium tetrasulfonated phthalocyanine by cytochrome c

The mechanism of electron-transfer from aluminium tetrasulfonated phthalocyanine triplet state to cytochrome c was investigated in this work. This reaction successfully quenches the dye triplet state due to the formation of complexes between the solute and the protein at the active site. The electron-transfer rate constant is around 3x10(7) s(-1), and is in accordance with previous results for the singlet excited state quenching [C.A.T. Laia, S.M.B. Costa, D. Phillips, A. Beeby. Electron-transfer kinetics in sulfonated aluminum phthalocyanines/cytochrome c complexes, J. Phys. Chem. B 108 (2004) 7506-7514.] in the framework of the Marcus theory, with a reorganization energy equal to 0.94 eV. The complex formation is diffusion controlled, but heterogeneities of the protein surface charge distribution lead to quenching rate constants smaller than predicted on a hard-spheres model with electrostatic interactions. Also the binding equilibrium constant is strongly affected by this phenomenon. Ionic strength plays an important role on the complex formation, but its effect on the unimolecular electron-transfer rate constant is negligible within experimental error.

Long-range protein electron transfer observed at the single-molecule level: In situ mapping of redox-gated tunneling resonance

A biomimetic long-range electron transfer (ET) system consisting of the blue copper protein azurin, a tunneling barrier bridge, and a gold single-crystal electrode was designed on the basis of molecular wiring self-assembly principles. This system is sufficiently stable and sensitive in a quasi-biological environment, suitable for detailed observations of long-range protein interfacial ET at the nanoscale and single-molecule levels. Because azurin is located at clearly identifiable fixed sites in well controlled orientation, the ET configuration parallels biological ET. The ET is nonadiabatic, and the rate constants display tunneling features with distance-decay factors of 0.83 and 0.91 A(-1) in H(2)O and D(2)O, respectively. Redox-gated tunneling resonance is observed in situ at the single-molecule level by using electrochemical scanning tunneling microscopy, exhibiting an asymmetric dependence on the redox potential. Maximum resonance appears around the equilibrium redox potential of azurin with an on/off current ratio of approximately 9. Simulation analyses, based on a two-step interfacial ET model for the scanning tunneling microscopy redox process, were performed and provide quantitative information for rational understanding of the ET mechanism.

Long-range electron transfer

Recent investigations have shed much light on the nuclear and electronic factors that control the rates of long-range electron tunneling through molecules in aqueous and organic glasses as well as through bonds in donor-bridge-acceptor complexes. Couplings through covalent and hydrogen bonds are much stronger than those across van der Waals gaps, and these differences in coupling between bonded and nonbonded atoms account for the dependence of tunneling rates on the structure of the media between redox sites in Ru-modified proteins and protein-protein complexes.