Alternative ground states enable pathway switching in biological electron transfer

Recent advances in tuning redox properties of electron transfer centers in metalloenzymes catalyzing oxygen reduction reaction and H2 oxidation important for fuel cells design

Current fuel-cell catalysts for oxygen reduction reaction (ORR) and H₂ oxidation use precious metals and, for ORR, require high overpotentials. In contrast, metalloenzymes perform their respective reaction at low overpotentials using earth-abundant metals, making metalloenzymes ideal candidates for inspiring electrocatalytic design. Critical to the success of these enzymes are redox-active metal centers surrounding the enzyme active sites that ensure fast electron transfer (ET) to or away from the active site, by tuning the catalytic potential of the reaction as observed in multicopper oxidases but also in dictating the catalytic bias of the reaction as realized in hydrogenases. This review summarizes recent advances in studying these ET centers in multicopper oxidases and heme-copper oxidases that perform ORR and hydrogenases in carrying out H₂ oxidation. Insights gained from understanding how the reduction potential of the ET centers effects reactivity at the active site in both the enzymes and their models are provided.

Bifurcation of excited state trajectories toward energy transfer or electron transfer directed by wave function symmetry

This work explores the concept that differential wave function overlap between excited states can be engineered within a molecular chromophore. The aim is to control excited state wave function symmetries, so that symmetry matches or mismatches result in differential orbital overlap and define low-energy trajectories or kinetic barriers within the excited state surface, that drive excited state population toward different reaction pathways. Two donor-acceptor assemblies were explored, where visible light absorption prepares excited states of different wave function symmetry. These states could be resolved using transient absorption spectroscopy, thanks to wave function symmetry-specific photoinduced optical transitions. One of these excited states undergoes energy transfer to the acceptor, while another undertakes a back-electron transfer to restate the ground state. This differential behavior is possible thanks to the presence of kinetic barriers that prevent excited state equilibration. This strategy can be exploited to avoid energy dissipation in energy conversion or photoredox catalytic schemes.

CuA-based chimeric T1 copper sites allow for independent modulation of reorganization energy and reduction potential

Attaining rational modulation of thermodynamic and kinetic redox parameters of metalloproteins is a key milestone towards the (re)design of proteins with new or improved redox functions. Here we report that implantation of ligand loops from natural T1 proteins into the scaffold of a CuA protein leads to a series of distorted T1-like sites that allow for independent modulation of reduction potentials (E°') and electron transfer reorganization energies (λ). On the one hand E°' values could be fine-tuned over 120 mV without affecting λ. On the other, λ values could be modulated by more than a factor of two while affecting E°' only by a few millivolts. These results are in sharp contrast to previous studies that used T1 cupredoxin folds, thus highlighting the importance of the protein scaffold in determining such parameters.

Unexpected electron spin density on the axial methionine ligand in CuA suggests its involvement in electron pathways

The CuA center is a paradigm for the study of long-range biological electron transfer. This metal center is an essential cofactor for terminal oxidases like cytochrome c oxidase, the enzymatic complex responsible for cellular respiration in eukaryotes and in most bacteria. CuA acts as an electron hub by transferring electrons from reduced cytochrome c to the catalytic site of the enzyme where dioxygen reduction takes place. Different electron transfer pathways have been proposed involving a weak axial methionine ligand residue, conserved in all CuA sites. This hypothesis has been challenged by theoretical calculations indicating the lack of electron spin density in this ligand. Here we report an NMR study with selectively labeled methionine in a native CuA. NMR spectroscopy discloses the presence of net electron spin density in the methionine axial ligand in the two alternative ground states of this metal center. Similar spin delocalization observed on two second sphere mutants further supports this evidence. These data provide a novel view of the electronic structure of CuA centers and support previously neglected electron transfer pathways.

Discrete Ligand Binding and Electron Transfer Properties of ba3-Cytochrome c Oxidase from Thermus thermophilus: Evolutionary Adaption to Low Oxygen and High Temperature Environments

Cytochrome c oxidase (C cO) couples the oxidation of cytochrome c to the reduction of molecular oxygen to water and links these electron transfers to proton translocation. The redox-driven C cO conserves part of the released free energy generating a proton motive force that leads to the synthesis of the main biological energy source ATP. Cytochrome ba₃ oxidase is a B-type oxidase from the extremely thermophilic eubacterium Thermus thermophilus with high O₂ affinity, expressed under elevated temperatures and limited oxygen supply and possessing discrete structural, ligand binding, and electron transfer properties. The origin and the cause of the peculiar, as compared to other C cOs, thermodynamic and kinetic properties remain unknown. Fourier transform infrared (FTIR) and time-resolved step-scan FTIR (TRS²-FTIR) spectroscopies have been employed to investigate the origin of the binding and electron transfer properties of cytochrome ba₃ oxidase in both the fully reduced (FR) and mixed valence (MV) forms. Several independent and not easily separated factors leading to increased thermostability and high O₂ affinity have been determined. These include (i) the increased hydrophobicity of the active center, (ii) the existence of a ligand input channel, (iii) the high affinity of Cu_B for exogenous ligands, (iv) the optimized electron transfer (ET) pathways, (v) the effective proton-input channel and water-exit pathway as well the proton-loading/exit sites, (vi) the specifically engineered protein structure, and (vii) the subtle thermodynamic and kinetic regulation. We correlate the unique ligand binding and electron transfer properties of cytochrome ba₃ oxidase with the existence of an adaption mechanism which is necessary for efficient function. These results suggest that a cascade of structural factors have been optimized by evolution, through protein architecture, to ensure the conversion of cytochrome ba₃ oxidase into a high O₂-affinity enzyme that functions effectively in its extreme native environment. The present results show that ba₃-cytochrome c oxidase uses a unique structural pattern of energy conversion that has taken into account all the extreme environmental factors that affect the function of the enzyme and is assembled in such a way that its exclusive functions are secured. Based on the available data of CcOs, we propose possible factors including the rigidity and nonpolar hydrophobic interactions that contribute to the behavior observed in cytochrome ba₃ oxidase.

Formation and Electronic Structure of an Atypical CuA Site

PmoD, a recently discovered protein from methane-oxidizing bacteria, forms a homodimer with a dicopper CuA center at the dimer interface. Although the optical and electron paramagnetic resonance (EPR) spectroscopic signatures of the PmoD CuA bear similarities to those of canonical CuA sites, there are also some puzzling differences. Here we have characterized the rapid formation (seconds) and slow decay (hours) of this homodimeric CuA site to two mononuclear Cu2+ sites, as well as its electronic and geometric structure, using stopped-flow optical and advanced paramagnetic resonance spectroscopies. PmoD CuA formation occurs rapidly and involves a short-lived intermediate with a λmax of 360 nm. Unlike other CuA sites, the PmoD CuA is unstable, decaying to two type 2 Cu2+ centers. Surprisingly, NMR data indicate that the PmoD CuA has a pure σu* ground state rather than the typical equilibrium between σu* and πu of all other CuA proteins. EPR, ENDOR, ESEEM, and HYSCORE data indicate the presence of two histidine and two cysteine ligands coordinating the CuA core in a highly symmetrical fashion. This report significantly expands the diversity and understanding of known CuA sites.

HYSCORE Insights into the Distribution of the Unpaired Spin Density in an Engineered CuA Site in Azurin and Its His120Gly Variant

A comparative study of the ¹H and ¹⁴N hyperfine interactions between the Cu_A site in an engineered Cu_A center in azurin (WT-Cu_AAz) and its His120Gly variant (H120G-Cu_AAz) using the two-dimensional ESEEM technique, HYSCORE, is reported. HYSCORE spectroscopy has clarified conflicting results in previous electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) studies and found clear differences between the two Cu_A azurins. Specifically, a hyperfine coupling A_N⊥ of 15.3 MHz was determined for the first time from the frequencies of double-quantum transitions of ¹⁴N histidine nitrogens coordinated to Cu_A in WT-Cu_AAz. In contrast, such coupling was not observed in the spectra of H120G-Cu_AAz, indicating at least a several megahertz increase in A_N⊥ for the coordinated nitrogen in this variant. In addition, ¹⁴N HYSCORE spectra of WT-Cu_AAz show interaction with only one type of weakly coupled nitrogen assigned to the remote N_ε atom of coordinated imidazole residues based on the quadrupole coupling constant ( e² Qq/4 h) of ∼0.4 MHz. The spectrum of H120G-Cu_AAz resolves additional features typical for backbone peptide nitrogens with larger e² Qq/4 h values of ∼0.7 MHz. Hyperfine couplings with these nitrogens vary between ∼0.4 and 0.7 MHz. In addition, the two resolved cross-peaks from C_β protons in H120G-Cu_AAz display only ∼1 MHz shifts relative to the corresponding peaks in WT-Cu_AAz. These new findings have provided the first experimental evidence of the previous density functional theory analysis that predicted changes in the delocalized electron spin population of ∼0.02-0.03 (i.e., ∼10%) on copper and sulfur atoms of the Cu_A center in H120 variants relative to WT-Cu_AAz and resolved contradicting results between EPR and ENDOR studies of the valence distribution in Cu_AAz and its variants.

Fine Tuning of Functional Features of the CuA Site by Loop-Directed Mutagenesis

Here we report the spectroscopic and electrochemical characterization of three novel chimeric Cu_A proteins in which either one or the three loops surrounding the metal ions in the Thermus thermophilus protein have been replaced by homologous human and plant sequences while preserving the set of coordinating amino acids. These conservative modifications mimic basic differences between Cu_A sites from different organisms and allow for fine tuning the energy gap between alternative electronic ground states of Cu_A.. This results in a systematic modulation of thermodynamic and kinetic electron transfer (ET) parameters and in the selection of one of two possible redox-active molecular orbitals, which differ in the ET reorganization energy by a factor of 2. Moreover, the ET mechanism is found to be frictionally controlled, and the modifications introduced into the different chimeras do not affect the frictional activation parameter.

Dramatic Electronic Perturbations of CuA Centers via Subtle Geometric Changes

Cu_A is a binuclear copper site acting as electron entry port in terminal heme-copper oxidases. In the oxidized form, Cu_A is a mixed valence pair whose electronic structure can be described using a potential energy surface with two minima, σ_u* and π_u, that are variably populated at room temperature. We report that mutations in the first and second coordination spheres of the binuclear metallocofactor can be combined in an additive manner to tune the energy gap and, thus, the relative populations of the two lowest-lying states. A series of designed mutants span σ_u*/π_u energy gaps ranging from 900 to 13 cm^-1. The smallest gap corresponds to a variant with an effectively degenerate ground state. All engineered sites preserve the mixed-valence character of this metal center and the electron transfer functionality. An increase of the Cu-Cu distance less than 0.06 Å modifies the σ_u*/π_u energy gap by almost 2 orders of magnitude, with longer distances eliciting a larger population of the π_u state. This scenario offers a stark contrast to synthetic systems, as model compounds require a lengthening of 0.5 Å in the Cu-Cu distance to stabilize the π_u state. These findings show that the tight control of the protein environment allows drastic perturbations in the electronic structure of Cu_A sites with minor geometric changes.

DEPC modification of the CuA protein from Thermus thermophilus

The Cu_A center is the initial electron acceptor in cytochrome c oxidase, and it consists of two copper ions bridged by two cysteines and ligated by two histidines, a methionine, and a carbonyl in the peptide backbone of a nearby glutamine. The two ligating histidines are of particular interest as they may influence the electronic and redox properties of the metal center. To test for the presence of reactive ligating histidines, a portion of cytochrome c oxidase from the bacteria Thermus thermophilus that contains the Cu_A site (the TtCu_A protein) was treated with the chemical modifier diethyl pyrocarbonate (DEPC) and the reaction followed through UV-visible, circular dichroism, and electron paramagnetic resonance spectroscopies at pH 5.0-9.0. A mutant protein (H40A/H117A) with the non-ligating histidines removed was similarly tested. Introduction of an electron-withdrawing DEPC-modification onto the ligating histidine 157 of TtCu_A increased the reduction potential by over 70 mV, as assessed by cyclic voltammetry. Results from both proteins indicate that DEPC reacts with one of the two ligating histidines, modification of a ligating histidine raises the reduction potential of the Cu_A site, and formation of the DEPC adduct is reversible at room temperature. The existence of the reactive ligating histidine suggests that this residue may play a role in modulating the electronic and redox properties of TtCu_A through kinetically-controlled proton exchange with the solvent. Lack of reactivity by the metalloproteins Sco and azurin, both of which contain a mononuclear copper center, indicate that reactivity toward DEPC is not a characteristic of all ligating histidines.

Tuning of Enthalpic/Entropic Parameters of a Protein Redox Center through Manipulation of the Electronic Partition Function

Manipulation of the partition function (Q) of the redox center Cu_A from cytochrome c oxidase is attained by tuning the accessibility of a low lying alternative electronic ground state and by perturbation of the electrostatic potential through point mutations, loop engineering and pH variation. We report clear correlations of the entropic and enthalpic contributions to redox potentials with Q and with the identity and hydrophobicity of the weak axial ligand, respectively.

Design of a single protein that spans the entire 2-V range of physiological redox potentials

The reduction potential (E°') is a critical parameter in determining the efficiency of most biological and chemical reactions. Biology employs three classes of metalloproteins to cover the majority of the 2-V range of physiological E°'s. An ultimate test of our understanding of E°' is to find out the minimal number of proteins and their variants that can cover this entire range and the structural features responsible for the extreme E°'. We report herein the design of the protein azurin to cover a range from +970 mV to -954 mV vs. standard hydrogen electrode (SHE) by mutating only five residues and using two metal ions. Spectroscopic methods have revealed geometric parameters important for the high E°'. The knowledge gained and the resulting water-soluble redox agents with predictable E°'s, in the same scaffold with the same surface properties, will find wide applications in chemical, biochemical, biophysical, and biotechnological fields.

Reversible Switching of Redox-Active Molecular Orbitals and Electron Transfer Pathways in Cu(A) Sites of Cytochrome c Oxidase

The Cu(A) site of cytochrome c oxidase is a redox hub that participates in rapid electron transfer at low driving forces with two redox cofactors in nearly perpendicular orientations. Spectroscopic and electrochemical characterizations performed on first and second-sphere mutants have allowed us to experimentally detect the reversible switching between two alternative electronic states that confer different directionalities to the redox reaction. Specifically, the M160H variant of a native Cu(A) shows a reversible pH transition that allows to functionally probe both states in the same protein species. Alternation between states exerts a dramatic impact on the kinetic redox parameters, thereby suggesting this effect as the mechanism underlying the efficiency and directionality of Cu(A) electron transfer in vivo. These findings may also prove useful for the development of molecular electronics.

Long-Range Electron Transfer in Engineered Azurins Exhibits Marcus Inverted Region Behavior

The Marcus theory of electron transfer (ET) predicts that while the ET rate constants increase with rising driving force until it equals a reaction's reorganization energy, at higher driving force the ET rate decreases, having reached the Marcus inverted region. While experimental evidence of the inverted region has been reported for organic and inorganic ET reactions as well as for proteins conjugated with ancillary redox moieties, evidence of the inverted region in a "protein-only" system has remained elusive. We herein provide such evidence in a series of nonderivatized proteins. These results may facilitate the design of ET centers for future applications such as advanced energy conversions.

The role of protein dynamics and thermal fluctuations in regulating cytochrome c/cytochrome c oxidase electron transfer

In this overview we present recent combined electrochemical, spectroelectrochemical, spectroscopic and computational studies from our group on the electron transfer reactions of cytochrome c and of the primary electron acceptor of cytochrome c oxidase, the CuA site, in biomimetic complexes. Based on these results, we discuss how protein dynamics and thermal fluctuations may impact on protein ET reactions, comment on the possible physiological relevance of these results, and finally propose a regulatory mechanism that may operate in the Cyt/CcO electron transfer reaction in vivo. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Molecular dynamics simulations of apocupredoxins: insights into the formation and stabilization of copper sites under entatic control

Cupredoxins perform copper-mediated long-range electron transfer (ET) in biological systems. Their copper-binding sites have evolved to force copper ions into ET-competent systems with decreased reorganization energy, increased reduction potential, and a distinct electronic structure compared with those of non-ET-competent copper complexes. The entatic or rack-induced state hypothesis explains these special properties in terms of the strain that the protein matrix exerts on the metal ions. This idea is supported by X-ray structures of apocupredoxins displaying "closed" arrangements of the copper ligands like those observed in the holoproteins; however, it implies completely buried copper-binding atoms, conflicting with the notion that they must be exposed for copper loading. On the other hand, a recent work based on NMR showed that the copper-binding regions of apocupredoxins are flexible in solution. We have explored five cupredoxins in their "closed" apo forms through molecular dynamics simulations. We observed that prearranged ligand conformations are not stable as the X-ray data suggest, although they do form part of the dynamic landscape of the apoproteins. This translates into variable flexibility of the copper-binding regions within a rigid fold, accompanied by fluctuations of the hydrogen bonds around the copper ligands. Major conformations with solvent-exposed copper-binding atoms could allow initial binding of the copper ions. An eventual subsequent incursion to the closed state would result in binding of the remaining ligands, trapping the closed conformation thanks to the additional binding energy and the fastening of noncovalent interactions that make up the rack.

Native Cu(A) redox sites are largely resilient to pH variations within a physiological range

Previous studies on engineered CuA centres have shown that one of the histidine ligands is protonated and dissociated from the metal site at physiological pH values, thus suggesting a role in regulating proton-coupled electron transfer of cytochrome c oxidases in vivo. Here we report that for native CuA such protonation does not take place at physiologically relevant pH values and, furthermore, no significant changes in the spectroscopic and redox properties of the metal site occur at low pH.

Redox-state sensing by hydrogen bonds in the CuA center of cytochrome c oxidase

Cytochrome c oxidases (CcO) couple electron transfer to active proton translocation through a gated mechanism that minimizes energy losses by preventing protons from flowing backwards or leaking. Such a complex mechanism requires that information about the redox and protonation states of the different centers be transmitted between different parts of the oxidase. Here we report a network of residues located around the electron entry point of CcO, the CuA site in subunit II, that experience collective pH equilibria around neutral pH. This network starts at the occluded side of the CuA site and extends to the interface between subunits I and II of the CcO, where the proton exit is located and through which electrons flow into subunit I. One of the residues in this network is directly involved in a hydrogen bond to one of the CuA ligands, whose strength is highly sensitive to the redox state of the metal center. We propose that this interaction mediates the transmission of redox changes from ET centers to other functional regions of the oxidase, and possibly also in other similar machineries, as part of their gating and regulatory mechanisms.

Axial interactions in the mixed-valent CuA active site and role of the axial methionine in electron transfer

Within Cu-containing electron transfer active sites, the role of the axial ligand in type 1 sites is well defined, yet its role in the binuclear mixed-valent CuA sites is less clear. Recently, the mutation of the axial Met to Leu in a CuA site engineered into azurin (CuA Az) was found to have a limited effect on E(0) relative to this mutation in blue copper (BC). Detailed low-temperature absorption and magnetic circular dichroism, resonance Raman, and electron paramagnetic resonance studies on CuA Az (WT) and its M123X (X = Q, L, H) axial ligand variants indicated stronger axial ligation in M123L/H. Spectroscopically validated density functional theory calculations show that the smaller ΔE(0) is attributed to H2O coordination to the Cu center in the M123L mutant in CuA but not in the equivalent BC variant. The comparable stabilization energy of the oxidized over the reduced state in CuA and BC (CuA ∼ 180 mV; BC ∼ 250 mV) indicates that the S(Met) influences E(0) similarly in both. Electron delocalization over two Cu centers in CuA was found to minimize the Jahn-Teller distortion induced by the axial Met ligand and lower the inner-sphere reorganization energy. The Cu-S(Met) bond in oxidized CuA is weak (5.2 kcal/mol) but energetically similar to that of BC, which demonstrates that the protein matrix also serves an entatic role in keeping the Met bound to the active site to tune down E(0) while maintaining a low reorganization energy required for rapid electron transfer under physiological conditions.