Faster interprotein electron transfer in a [myoglobin, b⁵] complex with a redesigned 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.

Exploration of reaction rates of chlorine dioxide with tryptophan residue in oligopeptides and proteins

Chlorine dioxide (ClO₂), an alternative disinfectant to chlorine, has a superior ability to inactivate microorganisms, in which protein damage has been considered as the main inactivation mechanism. However, the reactivity of ClO₂ with amino acid residues in oligopeptides and proteins remains poorly investigated. In this research, we studied the reaction rate constants of ClO₂ with tryptophan residues in five heptapeptides and four proteins using stopped-flow or competition kinetic method. Each heptapeptide and protein contain only one tryptophan residue and the reactivity of tryptophan residue with ClO₂ was lower than that of free tryptophan (3.88 × 10⁴ (mol/L)^-1sec^-1 at pH 7.0). The neighboring amino acid residues affected the reaction rates through promoting inter-peptide aggregation, changing electron density, shifting pK_a values or inducing electron transfer via redox reactions. A single amino acid residue difference in oligopeptides can make the reaction rate constants differ by over 60% (e.g. 3.01 × 10⁴ (mol/L)^-1sec^-1 for DDDWNDD and 1.85 × 10⁴ (mol/L)^-1sec^-1 for DDDWDDD at pH 7.0 (D: aspartic acid, W: tryptophan, N: asparagine)). The reaction rates of tryptophan-containing oligopeptides were also highly pH-dependent with higher reactivity for deprotonated tryptophan than the neutral specie. Tryptophan residues in proteins spanned a 4-fold range reactivity toward ClO₂ (i.e. 0.84 × 10⁴ (mol/L)^-1sec^-1 for ribonuclease T1 and 3.21 × 10⁴ (mol/L)^-1sec^-1 for melittin at pH 7.0) with accessibility to the oxidant as the determinating factor. The local environment surrounding the tryptophan residue in proteins can also accelerate the reaction rates by increasing the electron density of the indole ring of tryptophan or inhibit the reaction rates by inducing electron transfer reactions. The results are of significance in advancing understanding of ClO₂ oxidative reactions with proteins and microbial inactivation mechanisms.

Electron Transfer Methods in Open Systems

Utilization of electron transfer methods for description of quantum transport is popular due to simplicity of the formulation and its ability to account for basic physics of electron exchange between the system and baths. At the same time, the necessity to go beyond simple golden rule-type expressions for rates was indicated in the literature and ad hoc formulations were proposed. Similarly, kinetic schemes for quantum transport beyond the usual second-order Lindblad/Redfield considerations were discussed. Here we utilize recently introduced the nonequilibrium Hubbard Green's function diagrammatic technique to analyze the construction of rates in open systems. We show that previous considerations for rates of second and fourth order can be obtained as a particular case of zero- and second-order Green's function diagrammatic series with bare diagrams. We discuss limitations of previous considerations, stress advantages of the Hubbard Green's function approach in constructing the rates, and indicate that standard dressing of the diagrams is a natural way to account for additional baths/degrees of freedom in the formulation of generalized expressions for the rates.

Expansion of Redox Chemistry in Designer Metalloenzymes

Many artificial enzymes that catalyze redox reactions have important energy, environmental, and medical applications. Native metalloenzymes use a set of redox-active amino acids and cofactors as redox centers, with a potential range between -700 and +800 mV versus standard hydrogen electrode (SHE, all reduction potentials are versus SHE). The redox potentials and the orientation of redox centers in native metalloproteins are optimal for their redox chemistry. However, the limited number and potential range of native redox centers challenge the design and optimization of novel redox chemistry in metalloenzymes. Artificial metalloenzymes use non-native redox centers and could go far beyond the natural range of redox potentials for novel redox chemistry. In addition to designing protein monomers, strategies for increasing the electron transfer rate in self-assembled protein complexes and protein-electrode or -nanomaterial interfaces will be discussed. Redox reactions in proteins occur on redox active amino acid residues (Tyr, Trp, Met, Cys, etc.) and cofactors (iron sulfur clusters, flavin, heme, etc.). The redox potential of these redox centers cover a ∼1.5 V range and is optimized for their specific functions. Despite recent progress, tuning the redox potential for amino acid residues or cofactors remains challenging. Many redox-active unnatural amino acids (UAAs) can be incorporated into protein via genetic codon expansion. Their redox potentials extend the range of physiologically relevant potentials. Indeed, installing new redox cofactors with fined-tuned redox potentials is essential for designing novel redox enzymes. By combining UAA and redox cofactor incorporation, we harnessed light energy to reduce CO₂ in a fluorescent protein, mimicking photosynthetic apparatus in nature. Manipulating the position and reduction potential of redox centers inside proteins is important for optimizing the electron transfer rate and the activity of artificial enzymes. Learning from the native electron transfer complex, protein-protein interactions can be enhanced by increasing the electrostatic interaction between proteins. An artificial oxidase showed close to native enzyme activity with optimized interaction with electron transfer partner and increased electron transfer efficiency. In addition to the de novo design of protein-protein interaction, protein self-assembly methods using scaffolds, such as proliferating cell nuclear antigen, to efficiently anchor enzymes and their redox partners. The self-assembly process enhances electron transfer efficiency and enzyme activity by bringing redox centers into close proximity of each other. In addition to protein self-assembly, protein-electrode or protein-nanomaterial self-assembly can also promote efficient electron transfer from inorganic materials to enzyme active sites. Such hybrid systems combine the efficiency of enzyme reactions and the robustness of electrodes or nanomaterials, often with advantageous catalytic activities. By combining these strategies, we can not only mimic some of nature's most fascinating reactions, such as photosynthesis and aerobic respiration, but also transcend nature toward environmental, energy, and health applications.

A novel role for PsbO1 in photosynthetic electron transport as suggested by its light-triggered selective nitration in Arabidopsis thaliana

Exposure of Arabidopsis leaves to nitrogen dioxide (NO2) results in the selective nitration of specific proteins, such as PsbO1. The 9th tyrosine residue (9Tyr) of PsbO1 has been identified as the nitration site. This nitration is triggered by light and inhibited by photosynthetic electron transport inhibitors. During protein nitration, tyrosyl and NO2 radicals are formed concurrently and combine rapidly to form 3-nitrotyrosine. A selective oxidation mechanism for 9Tyr of PsbO1 is required. We postulated that, similar to 161Tyr of D1, 9Tyr of PsbO1 is selectively photo-oxidized by photosynthetic electron transport in response to illumination to a tyrosyl radical. In corroboration, after reappraising our oxygen evolution analysis, the nitration of PsbO1 proved responsible for decreased oxygen evolution from the thylakoid membranes. NO2 is reportedly taken into cells as nitrous acid, which dissociates to form NO2-. NO2- may be oxidized into NO2 by the oxygen-evolving complex. Light may synchronize this reaction with tyrosyl radical formation. These findings suggest a novel role for PsbO1 in photosynthetic electron transport.

Insights Into How Heme Reduction Potentials Modulate Enzymatic Activities of a Myoglobin-based Functional Oxidase

Heme-copper oxidase (HCO) is a class of respiratory enzymes that use a heme-copper center to catalyze O₂ reduction to H₂ O. While heme reduction potential (E°') of different HCO types has been found to vary >500 mV, its impact on HCO activity remains poorly understood. Here, we use a set of myoglobin-based functional HCO models to investigate the mechanism by which heme E°' modulates oxidase activity. Rapid stopped-flow kinetic measurements show that increasing heme E°' by ca. 210 mV results in increases in electron transfer (ET) rates by 30-fold, rate of O₂ binding by 12-fold, O₂ dissociation by 35-fold, while decreasing O₂ affinity by 3-fold. Theoretical calculations reveal that E°' modulation has significant implications on electronic charge of both heme iron and O₂ , resulting in increased O₂ dissociation and reduced O₂ affinity at high E°' values. Overall, this work suggests that fine-tuning E°' in HCOs and other heme enzymes can modulate their substrate affinity, ET rate and enzymatic activity.

Glutamate 350 Plays an Essential Role in Conformational Gating of Long-Range Radical Transport in Escherichia coli Class Ia Ribonucleotide Reductase

Escherichia coli class Ia ribonucleotide reductase (RNR) is composed of two subunits that form an active α2β2 complex. The nucleoside diphosphate substrates (NDP) are reduced in α2, 35 Å from the essential diferric-tyrosyl radical (Y₁₂₂^•) cofactor in β2. The Y₁₂₂^•-mediated oxidation of C₄₃₉ in α2 occurs by a pathway (Y₁₂₂ ⇆ [W₄₈] ⇆ Y₃₅₆ in β2 to Y₇₃₁ ⇆ Y₇₃₀ ⇆ C₄₃₉ in α2) across the α/β interface. The absence of an α2β2 structure precludes insight into the location of Y₃₅₆ and Y₇₃₁ at the subunit interface. The proximity in the primary sequence of the conserved E₃₅₀ to Y₃₅₆ in β2 suggested its importance in catalysis and/or conformational gating. To study its function, pH-rate profiles of wild-type β2/α2 and mutants in which 3,5-difluorotyrosine (F₂Y) replaces residue 356, 731, or both are reported in the presence of E₃₅₀ or E₃₅₀X (X = A, D, or Q) mutants. With E₃₅₀, activity is maintained at the pH extremes, suggesting that protonated and deprotonated states of F₂Y₃₅₆ and F₂Y₇₃₁ are active and that radical transport (RT) can occur across the interface by proton-coupled electron transfer at low pH or electron transfer at high pH. With E₃₅₀X mutants, all RNRs were inactive, suggesting that E₃₅₀ could be a proton acceptor during oxidation of the interface Ys. To determine if E₃₅₀ plays a role in conformational gating, the strong oxidants, NO₂Y₁₂₂^•-β2 and 2,3,5-F₃Y₁₂₂^•-β2, were reacted with α2, CDP, and ATP in E₃₅₀ and E₃₅₀X backgrounds and the reactions were monitored for pathway radicals by rapid freeze-quench electron paramagnetic resonance spectroscopy. Pathway radicals are generated only when E₃₅₀ is present, supporting its essential role in gating the conformational change(s) that initiates RT and masking its role as a proton acceptor.

Improving artificial metalloenzymes' activity by optimizing electron transfer

This feature article discusses the strategies to optimize electron transfer efficiency, towards enhancing the activity of artificial metalloenzymes.

Using Biosynthetic Models of Heme-Copper Oxidase and Nitric Oxide Reductase in Myoglobin to Elucidate Structural Features Responsible for Enzymatic Activities

In biology, a heme-Cu center in heme-copper oxidases (HCOs) is used to catalyze the four-electron reduction of oxygen to water, while a heme-nonheme diiron center in nitric oxide reductases (NORs) is employed to catalyze the two-electron reduction of nitric oxide to nitrous oxide. Although much progress has been made in biochemical and biophysical studies of HCOs and NORs, structural features responsible for similarities and differences within the two enzymatic systems remain to be understood. Here, we discuss the progress made in the design and characterization of myoglobin-based enzyme models of HCOs and NORs. In particular, we focus on use of these models to understand the structure-function relations between HCOs and NORs, including the role of nonheme metals, conserved amino acids in the active site, heme types and hydrogen-bonding network in tuning enzymatic activities and total turnovers. Insights gained from these studies are summarized and future directions are proposed.

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.

Design and fine-tuning redox potentials of metalloproteins involved in electron transfer in bioenergetics

Redox potentials are a major contributor in controlling the electron transfer (ET) rates and thus regulating the ET processes in the bioenergetics. To maximize the efficiency of the ET process, one needs to master the art of tuning the redox potential, especially in metalloproteins, as they represent major classes of ET proteins. In this review, we first describe the importance of tuning the redox potential of ET centers and its role in regulating the ET in bioenergetic processes including photosynthesis and respiration. The main focus of this review is to summarize recent work in designing the ET centers, namely cupredoxins, cytochromes, and iron-sulfur proteins, and examples in design of protein networks involved these ET centers. We then discuss the factors that affect redox potentials of these ET centers including metal ion, the ligands to metal center and interactions beyond the primary ligand, especially non-covalent secondary coordination sphere interactions. We provide examples of strategies to fine-tune the redox potential using both natural and unnatural amino acids and native and nonnative cofactors. Several case studies are used to illustrate recent successes in this area. Outlooks for future endeavors are also provided. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Structural basis of interprotein electron transfer in bacterial sulfite oxidation

Interprotein electron transfer underpins the essential processes of life and relies on the formation of specific, yet transient protein-protein interactions. In biological systems, the detoxification of sulfite is catalyzed by the sulfite-oxidizing enzymes (SOEs), which interact with an electron acceptor for catalytic turnover. Here, we report the structural and functional analyses of the SOE SorT from Sinorhizobium meliloti and its cognate electron acceptor SorU. Kinetic and thermodynamic analyses of the SorT/SorU interaction show the complex is dynamic in solution, and that the proteins interact with Kd = 13.5 ± 0.8 μM. The crystal structures of the oxidized SorT and SorU, both in isolation and in complex, reveal the interface to be remarkably electrostatic, with an unusually large number of direct hydrogen bonding interactions. The assembly of the complex is accompanied by an adjustment in the structure of SorU, and conformational sampling provides a mechanism for dissociation of the SorT/SorU assembly.

Structure and Function of Transient Encounters of Redox Proteins

Many biomolecular interactions proceed via lowly populated, transient intermediates. Believed to facilitate formation of a productive complex, these short-lived species are inaccessible to conventional biophysical and structural techniques and, until recently, could only be studied by theoretical simulations. Recent development of experimental approaches sensitive to the presence of minor species--in particular paramagnetic relaxation enhancement (PRE) NMR spectroscopy--has enabled direct visualization and detailed characterization of such lowly populated states. Collectively referred to as an encounter complex, the binding intermediates are particularly important in transient protein interactions, such as those orchestrating signaling cascades or energy-generating electron transfer (ET) chains. Here I discuss encounter complexes of redox proteins mediating biological ET reactions, which are essential for many vital cellular activities including oxidative phosphorylation and photosynthesis. In particular, this Account focuses on the complex of cytochrome c (Cc) and cytochrome c peroxidase (CcP), which is a paradigm of biomolecular ET and an attractive system for studying protein binding and enzymatic catalysis. The Cc-CcP complex formation proceeds via an encounter state, consisting of multiple protein-protein orientations sampled in the search of the dominant, functionally active bound form and exhibiting a broad spatial distribution, in striking agreement with earlier theoretical simulations. At low ionic strength, CcP binds another Cc molecule to form a weak ternary complex, initially inferred from kinetics experiments and postulated to account for the measured ET activity. Despite strenuous efforts, the ternary complex could not be observed directly and remained eagerly sought for the past two decades. Very recently, we have solved its structure in solution and shown that it consists of two binding forms: the dominant, ET-inactive geometry and an ensemble of lowly populated species with short separations between Cc and CcP cofactors, which summarily account for the measured ET rate. Unlike most protein complexes, which require accurate alignment of the binding surfaces in a single, well-defined orientation to carry out their function, redox proteins can form multiple productive complexes. As fast ET will occur any time the redox centers of the binding partners are close enough to ensure efficient electron tunneling across the interface, many protein-protein orientations are expected to be ET active. The present analysis confirms that the low-occupancy states can support the functional ET activity and contribute to the stability of redox protein complexes. As illustrated here, boundaries between the dominant and the encounter forms become blurred for many dynamic ET systems, which are more aptly described by ensembles of functionally and structurally heterogeneous bound forms.

"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.

Evidence for Functionally Relevant Encounter Complexes in Nitrogenase Catalysis

Nitrogenase is the only enzyme that can convert atmospheric dinitrogen (N2) into biologically usable ammonia (NH3). To achieve this multielectron redox process, the nitrogenase component proteins, MoFe-protein (MoFeP) and Fe-protein (FeP), repeatedly associate and dissociate in an ATP-dependent manner, where one electron is transferred from FeP to MoFeP per association. Here, we provide experimental evidence that encounter complexes between FeP and MoFeP play a functional role in nitrogenase catalysis. The encounter complexes are stabilized by electrostatic interactions involving a positively charged patch on the β-subunit of MoFeP. Three single mutations (βAsn399Glu, βLys400Glu, and βArg401Glu) in this patch were generated in Azotobacter vinelandii MoFeP. All of the resulting variants displayed decreases in specific catalytic activity, with the βK400E mutation showing the largest effect. As simulated by the Thorneley-Lowe kinetic scheme, this single mutation lowered the rate constant for FeP-MoFeP association 5-fold. We also found that the βK400E mutation did not affect the coupling of ATP hydrolysis with electron transfer (ET) between FeP and MoFeP. These data suggest a mechanism where FeP initially forms encounter complexes on the MoFeP β-subunit surface en route to the ATP-activated, ET-competent complex over the αβ-interface.

A Designed Metalloenzyme Achieving the Catalytic Rate of a Native Enzyme

Terminal oxidases catalyze four-electron reduction of oxygen to water, and the energy harvested is utilized to drive the synthesis of adenosine triphosphate. While much effort has been made to design a catalyst mimicking the function of terminal oxidases, most biomimetic catalysts have much lower activity than native oxidases. Herein we report a designed oxidase in myoglobin with an O₂ reduction rate (52 s^–1) comparable to that of a native cytochrome (cyt) cbb₃ oxidase (50 s^–1) under identical conditions. We achieved this goal by engineering more favorable electrostatic interactions between a functional oxidase model designed in sperm whale myoglobin and its native redox partner, cyt b₅, resulting in a 400-fold electron transfer (ET) rate enhancement. Achieving high activity equivalent to that of native enzymes in a designed metalloenzyme offers deeper insight into the roles of tunable processes such as ET in oxidase activity and enzymatic function and may extend into applications such as more efficient oxygen reduction reaction catalysts for biofuel cells.

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.

Bidirectional Photoinduced Electron Transfer in Ruthenium(II)-Tris-bipyridyl-Modified PpcA, a Multi-heme c-Type Cytochrome from Geobacter sulfurreducens

PpcA, a tri-heme cytochrome c7 from Geobacter sulfurreducens, was investigated as a model for photosensitizer-initiated electron transfer within a multi-heme "molecular wire" protein architecture. Escherichia coli expression of PpcA was found to be tolerant of cysteine site-directed mutagenesis, demonstrated by the successful expression of natively folded proteins bearing cysteine mutations at a series of sites selected to vary characteristically with respect to the three -CXXCH- heme binding domains. The introduced cysteines readily reacted with Ru(II)-(2,2'-bpy)2(4-bromomethyl-4'-methyl-2,2'-bipyridine) to form covalently linked constructs that support both photo-oxidative and photo-reductive quenching of the photosensitizer excited state, depending upon the initial heme redox state. Excited-state electron-transfer times were found to vary from 6 × 10(-12) to 4 × 10(-8) s, correlated with the distance and pathways for electron transfer. The fastest rate is more than 10(3)-fold faster than previously reported for photosensitizer-redox protein constructs using amino acid residue linking. Clear evidence for inter-heme electron transfer within the multi-heme protein is not detected within the lifetimes of the charge-separated states. These results demonstrate an opportunity to develop multi-heme c-cytochromes for investigation of electron transfer in protein "molecular wires" and to serve as frameworks for metalloprotein designs that support multiple-electron-transfer redox chemistry.

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.

Systematic tuning of heme redox potentials and its effects on O2 reduction rates in a designed oxidase in myoglobin

Cytochrome c Oxidase (CcO) is known to catalyze the reduction of O2 to H2O efficiently with a much lower overpotential than most other O2 reduction catalysts. However, methods by which the enzyme fine-tunes the reduction potential (E°) of its active site and the corresponding influence on the O2 reduction activity are not well understood. In this work, we report systematic tuning of the heme E° in a functional model of CcO in myoglobin containing three histidines and one tyrosine in the distal pocket of heme. By removing hydrogen-bonding interactions between Ser92 and the proximal His ligand and a heme propionate, and increasing hydrophobicity of the heme pocket through Ser92Ala mutation, we have increased the heme E° from 95 ± 2 to 123 ± 3 mV. Additionally, replacing the native heme b in the CcO mimic with heme a analogs, diacetyl, monoformyl, and diformyl hemes, that posses electron-withdrawing groups, resulted in higher E° values of 175 ± 5, 210 ± 6, and 320 ± 10 mV, respectively. Furthermore, O2 consumption studies on these CcO mimics revealed a strong enhancement in O2 reduction rates with increasing heme E°. Such methods of tuning the heme E° through a combination of secondary sphere mutations and heme substitutions can be applied to tune E° of other heme proteins, allowing for comprehensive investigations of the relationship between E° and enzymatic activity.

Long-range electron transfer with myoglobin immobilized at Au/mixed-SAM junctions: mechanistic impact of the strong protein confinement

Horse muscle myoglobin (Mb) was tightly immobilized at Au-deposited ~15-Å-thick mixed-type (1:1) alkanethiol SAMs, HS-(CH₂)₁₁-COOH/HS-(CH₂)₁₁-OH, and placed in contact with buffered H₂O or D₂O solutions. Fast-scan cyclic voltammetry (CV) and a Marcus-equation-based analysis were applied to determine unimolecular standard rate constants and reorganization free energies for electron transfer (ET), under variable-temperature (15-55 °C) and -pressure (0.01-150 MPa) conditions. The CV signal was surprisingly stable and reproducible even after multiple temperature and pressure cycles. The data analysis revealed the following values: standard rate constant, 33 s⁻¹ (25 °C, 0.01 MPa, H₂O); reorganization free energy, 0.5 ± 0.1 eV (throughout); activation enthalpy, 12 ± 3 kJ mol⁻¹; activation volume, -3.1 ± 0.2 cm³ mol⁻¹; and pH-dependent solvent kinetic isotope effect (k(H)⁰/k(D)⁰), 0.7-1.4. Furthermore, the values for the rate constant and reorganization free energy are very similar to those previously found for cytochrome c electrostatically immobilized at the monocomponent Au/HS-(CH₂)₁₁-COOH junction. In vivo, Mb apparently forms a natural electrostatic complex with cytochrome b₅ (cyt-b₅) through the "dynamic" (loose) docking pattern, allowing for a slow ET that is intrinsically coupled to the water's removal from the "defective" heme iron (altogether shaping the biological repair mechanism for Mb's "met" form). In contrary, our experiments rather mimic the case of a "simple" (tight) docking of the redesigned (mutant) Mb with cyt-b₅ (Nocek et al. J. Am. Chem. Soc. 2010, 132, 6165-6175). According to our analysis, in this configuration, Mb's distal pocket (linked to the "ligand channel") seems to be arrested within the restricted configuration, allowing the rate-determining reversible ET process to be coupled only to the inner-sphere reorganization (minimal elongation/shortening of an Fe-OH₂ bond) rather than the pronounced detachment (rebinding) of water and, hence, to be much faster.

Femtosecond dynamics of short-range protein electron transfer in flavodoxin

Intraprotein electron transfer (ET) in flavoproteins is important for understanding the correlation of their redox, configuration, and reactivity at the active site. Here, we used oxidized flavodoxin as a model system and report our complete characterization of a photoinduced redox cycle from the initial charge separation in 135-340 fs to subsequent charge recombination in 0.95-1.6 ps and to the final cooling relaxation of the product(s) in 2.5-4.3 ps. With 11 mutations at the active site, we observed that these ultrafast ET dynamics, much faster than active-site relaxation, mainly depend on the reduction potentials of the electron donors with minor changes caused by mutations, reflecting a highly localized ET reaction between the stacked donor and acceptor at a van der Waals distance and leading to a gas-phase type of bimolecular ET reaction confined in the active-site nanospace. Significantly, these ultrafast ET reactions ensure our direct observation of vibrationally excited reaction product(s), suggesting that the back ET barrier is effectively reduced because of the decrease in the total free energy in the Marcus inverted region, leading to the accelerated charge recombination. Such vibrationally coupled charge recombination should be a general feature of flavoproteins with similar configurations and interactions between the cofactor flavin and neighboring aromatic residues.

Formation of transient protein complexes

The encounter complex of two proteins is a dynamic intermediate state that guides proteins to their binding site, thus enhancing the rate of complex formation. It is particularly useful for complexes that must balance a biological requirement for high turnover with the need for specific binding, such as electron transfer complexes. Here, we describe the current methods for studying and visualizing encounter complexes. We discuss recent developments in mapping the energy landscapes, the role of hydrophobic interactions during encounter complex formation and the discovery of futile encounter complexes. These studies have not only provided insight into encounter complexes of electron transfer proteins, but also opened up new questions and approaches for studying encounter complexes in other weakly associated proteins.

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.

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.

The key role of water in the dioxygenase function of Escherichia coli flavohemoglobin

Flavohemoglobins (FHbs) are members of the globin superfamily, widely distributed among prokaryotes and eukaryotes that have been shown to carry out nitric oxide dioxygenase (NOD) activity. In prokaryotes, such as Escherichia coli, NOD activity is a defence mechanism against the NO release by the macrophages of the hosts' immune system during infection. Because of that, FHbs have been studied thoroughly and several drugs have been developed in an effort to fight infectious processes. Nevertheless, the protein's structural determinants involved in the NOD activity are still poorly understood. In this context, the aim of the present work is to unravel the molecular basis of FHbs structural dynamics-to-function relationship using state of the art computer simulation tools. In an effort to fulfill this goal, we studied three key processes that determine NOD activity, namely i) ligand migration into the active site ii) stabilization of the coordinated oxygen and iii) intra-protein electron transfer (ET). Our results allowed us to determine key factors related to all three processes like the presence of a long hydrophobic tunnel for ligand migration, the presence of a water mediated hydrogen bond to stabilize the coordinated oxygen and therefore achieve a high affinity, and the best possible ET paths between the FAD and the heme, where water molecules play an important role. Taken together the presented results close an important gap in our understanding of the wide and diverse globin structural-functional relationships.

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.

Dynamics in transient complexes of redox proteins

Recent studies have provided experimental information about the initial stage of protein complex formation, the encounter complex. This stage is particularly important in the weak and transient complexes formed between electron transfer proteins and their partners. These studies are discussed and the role of the encounter complex is interpreted in terms of the specific requirements that the biological function puts on these complexes.

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.