Electronic structure of the ground and excited states of the Cu(A) site by NMR spectroscopy

13C Direct Detected NMR for Challenging Systems

Thanks to recent improvements in NMR spectrometer hardware and pulse sequence design, modern ¹³C NMR has become a useful tool for biomolecular applications. The complete assignment of a protein can be accomplished by using ¹³C detected multinuclear experiments and it can provide unique information relevant for the study of a variety of different biomolecules including paramagnetic proteins and intrinsically disordered proteins. A wide range of NMR observables can be measured, concurring to the structural and dynamic characterization of a protein in isolation, as part of a larger complex, or even inside a living cell. We present the different properties of ¹³C with respect to ¹H, which provide the rationale for the experiments developed and their application, the technical aspects that need to be faced, and the many experimental variants designed to address different cases. Application areas where these experiments successfully complement proton NMR are also described.

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.

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.

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.

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.

Engineering a bifunctional copper site in the cupredoxin fold by loop-directed mutagenesis

Copper sites in proteins are designed to perform either electron transfer or redox catalysis. Type 1 and CuA sites are electron transfer hubs bound to a rigid protein fold that prevents binding of exogenous ligands and side reactions. Here we report the engineering of two Type 1 sites by loop-directed mutagenesis within a CuA scaffold with unique electronic structures and functional features. A copper-thioether axial bond shorter than the copper-thiolate bond is responsible for the electronic structure features, in contrast to all other natural or chimeric sites where the copper thiolate bond is short. These sites display highly unusual features, such as: (1) a high reduction potential despite a strong interaction with the axial ligand, which we attribute to changes in the hydrogen bond network and (2) the ability to bind exogenous ligands such as imidazole and azide. This strategy widens the possibility of using natural protein scaffolds with functional features not present in nature.

Walking the seven lines: binuclear copper A in cytochrome c oxidase and nitrous oxide reductase

The enzymes nitrous oxide reductase (N2OR) and cytochrome c oxidase (COX) are constituents of important biological processes. N2OR is the terminal reductase in a respiratory chain converting N₂O to N₂ in denitrifying bacteria; COX is the terminal oxidase of the aerobic respiratory chain of certain bacteria and eukaryotic organisms transforming O₂ to H₂O accompanied by proton pumping. Different spectroscopies including magnetic resonance techniques, were applied to show that N2OR has a mixed-valent Cys-bridged [Cu^1.5+(CyS)₂Cu^1.5+] copper site, and that such a binuclear center, called CuA, does also exist in COX. A sequence motif shared between the CuA center of N2OR and the subunit II of COX raises the issue of a putative evolutionary relationship of the two enzymes. The suggestion of a binuclear CuA in COX, with one unpaired electron delocalized between two equivalent Cu nuclei, was difficult to accept originally, even though regarded as a clever solution to many experimental observations. This minireview in honor of Helmut Sigel traces several of the critical steps forward in understanding the nature of CuA in N2OR and COX, and discusses its unique electronic features to some extent including the contributions made by the development of methodology and the discovery of a novel multi-copper enzyme. Left: X-band (9.130 GHz) and C-band (4.530 GHz, 1st harmonic display of experimental spectrum) EPR spectra of bovine heart cytochrome c oxidase, recorded at 20K. Right: Ribbon presentation of the CuA domain in cytochrome c oxidase and nitrous oxide reductase.

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.

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.

Spin-state-selective methods in solution- and solid-state biomolecular 13C NMR

Spin-state-selective methods to achieve homonuclear decoupling in the direct acquisition dimension of (13)C detected NMR experiments have been one of the key contributors to converting (13)C detected NMR experiments into really useful tools for studying biomolecules. We discuss here in detail the various methods that have been proposed, summarize the large array of new experiments that have been developed and present applications to different kinds of proteins in different aggregation states.

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.

The IR-¹⁵N-HSQC-AP experiment: a new tool for NMR spectroscopy of paramagnetic molecules

A crucial factor for the understanding of structure-function relationships in metalloproteins is the identification of NMR signals from residues surrounding the metal cofactor. When the latter is paramagnetic, the NMR information in the proximity of the metal center may be scarce, because fast nuclear relaxation quenches signal intensity and coherence transfer efficiency. To identify residues at a short distance from a paramagnetic center, we developed a modified version of the ¹⁵N-HSQC experiment where (1) an inversion recovery filter is added prior to HSQC, (2) the INEPT period has been optimized according to fast relaxation of interested spins, (3) the inverse INEPT has been eliminated and signals acquired as antiphase doublets. The experiment has been successfully tested on a human [Fe₂S₂] protein which is involved in the biogenesis of iron-sulfur proteins. Thirteen HN resonances, unobserved with conventional HSQC experiments, could be identified. The structural arrangement of the protein scaffold in the proximity of the Fe/S cluster is fundamental to comprehend the molecular processes responsible for the transfer of Fe/S groups in the iron-sulfur protein assembly machineries.

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.

Electron self-exchange of cytochrome c measured via13C detected protonless NMR

The use of protonless 13C′–13C′ EXSY (COCO-EXSY) is proposed here to measure electron self-exchange rates. The experiment is compared to the commonly employed 1H and 15N EXSY experiments using as a reference system human cytochrome c. In COCO-EXSY, the exchange peaks are stronger than in the other experiments with respect to the self peaks and their intensity is less dependent on the choice of the EXSY mixing time. The use of 13C directed detection may be essential for all those cases where T2 relaxation is detrimental, as in the case of proteins containing highly paramagnetic metal centers, or rotating slowly in solution, or where the amide signals are difficult to detect due to chemical or conformational exchange. The proposed experiment has a general applicability and can be used to monitor exchange phenomena different from electron self-exchange.

Electron spin density on the axial His ligand of high-spin and low-spin nitrophorin 2 probed by heteronuclear NMR spectroscopy

The electronic structure of heme proteins is exquisitely tuned by the interaction of the iron center with the axial ligands. NMR studies of paramagnetic heme systems have been focused on the heme signals, but signals from the axial ligands have been rather difficult to detect and assign. We report an extensive assignment of the (1)H, (13)C and (15)N resonances of the axial His ligand in the NO-carrying protein nitrophorin 2 (NP2) in the paramagnetic high-spin and low-spin forms, as well as in the diamagnetic NO complex. We find that the high-spin protein has σ spin delocalization to all atoms in the axial His57, which decreases in size as the number of bonds between Fe(III) and the atom in question increases, except that within the His57 imidazole ring the contact shifts are a balance between positive σ and negative π contributions. In contrast, the low-spin protein has π spin delocalization to all atoms of the imidazole ring. Our strategy, adequately combined with a selective residue labeling scheme, represents a straightforward characterization of the electron spin density in heme axial ligands.

Alternative ground states enable pathway switching in biological electron transfer

Electron transfer is the simplest chemical reaction and constitutes the basis of a large variety of biological processes, such as photosynthesis and cellular respiration. Nature has evolved specific proteins and cofactors for these functions. The mechanisms optimizing biological electron transfer have been matter of intense debate, such as the role of the protein milieu between donor and acceptor sites. Here we propose a mechanism regulating long-range electron transfer in proteins. Specifically, we report a spectroscopic, electrochemical, and theoretical study on WT and single-mutant Cu(A) redox centers from Thermus thermophilus, which shows that thermal fluctuations may populate two alternative ground-state electronic wave functions optimized for electron entry and exit, respectively, through two different and nearly perpendicular pathways. These findings suggest a unique role for alternative or "invisible" electronic ground states in directional electron transfer. Moreover, it is shown that this energy gap and, therefore, the equilibrium between ground states can be fine-tuned by minor perturbations, suggesting alternative ways through which protein-protein interactions and membrane potential may optimize and regulate electron-proton energy transduction.

Analysis of copper-ligand bond lengths in X-ray structures of different types of copper sites in proteins

An updated picture of the ligand sets and copper-ligand atom bond lengths in proteins is presented which takes advantage of (i) the approximately twofold increase in the number of entries for copper-containing proteins in the PDB since the last study of this kind, especially benefiting from the recent incorporation of the structures of proteins involved in copper homeostasis, and (ii) a preliminary classification of copper sites based on their structural, electronic and functional features. This classification allowed the calculation of reliable target copper-ligand distances for several bonds that were not available in previous work and that are in good agreement with EXAFS data and the known chemistry of these sites. The analysis presented here further disclosed an artifactual dependence of the average of the reported Cu-NHis bond lengths on structure resolution, highlighting the importance of taking this into account when computing target distances even from high-resolution structures. Finally, a relationship between the two Cu-O distances in bidentate carboxylates is disclosed, similar to that reported previously for other metal ions.

A theoretical investigation of the functional role of the axial methionine ligand of the Cu(A) site in cytochrome c oxidase

The functional roles of the amino acid residues of the Cu(A) site in bovine cytochrome c oxidase (CcO) were investigated by utilizing hybrid quantum mechanics (QM)/molecular mechanics (MM) calculations. The energy levels of the molecular orbitals (MOs) involving Cu d(zx) orbitals unexpectedly increased, as compared with those found previously with a simplified model system lacking the axial Met residue (i.e., Cu(2)S(2)N(2)). This elevation of MO energies stemmed from the formation of the anti-bonding orbitals, which are generated by hybridization between the d(zx) orbitals of Cu ions and the p-orbitals of the S and O atoms of the axial ligands. To clarify the roles of the axial Met ligand, the inner-sphere reorganization energies of the Cu(A) site were computed, with the Met residue assigned to either the QM or MM region. The reorganization energy slightly increased when the Met residue was excluded from the QM region. The existing experimental data and the present structural modeling study also suggested that the axial Met residue moderately increased the redox potential of the Cu(A) site. Thus, the role of the Met may be to regulate the electron transfer rate through the fine modulation of the electronic structure of the Cu(A) "platform", created by two Cys/His residues coordinated to the Cu ions. This regulation would provide the optimum redox potential/reorganization energy of the Cu(A) site, and thereby facilitate the subsequent cooperative reactions, such as the proton pump and the enzymatic activity, of CcO. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

NMR study of the exchange coupling in the trinuclear cluster of the multicopper oxidase Fet3p

Fet3p from Saccharomyces cerevisiae is a multicopper oxidase (MCO) which oxidizes Fe(2+) to Fe(3+). The electronic structure of the different copper centers in this family of enzymes has been extensively studied and discussed for years with a particular focus on the exchange coupling regime in the trinuclear cluster (TNC). Using NMR spectroscopy we have quantified the exchange coupling constant in the type 3 center in a fully metalated oxidase; this value in Fet3p is significantly higher than that reported for proteins containing isolated type 3 centers as tyrosinase. We also provide evidence of exchange coupling between the type 2 and the type 3 Cu(2+) ions, which supports the crystallographic evidence of dioxygen binding to the TNC. This work provides the foundation for the application of NMR to these complex systems.

Transforming a blue copper into a red copper protein: engineering cysteine and homocysteine into the axial position of azurin using site-directed mutagenesis and expressed protein ligation

Interactions of the axial ligand with its blue copper center are known to be important in tuning spectroscopic and redox properties of cupredoxins. While conversion of the blue copper center with a weak axial ligand to a green copper center containing a medium strength axial ligand has been demonstrated in cupredoxins, converting the blue copper center to a red copper center with a strong axial ligand has not been reported. Here we show that replacing Met121 in azurin from Pseudomonas aeruginosa with Cys caused an increased ratio (R(L)) of absorption at 447 nm over that at 621 nm. Whereas no axial Cu-S(Cys121) interaction in Met121Cys was detectable by extended X-ray absorption fine structure (EXAFS) spectroscopy at pH 5, similar to what was observed in native azurin with Met121 as the axial ligand, the Cu-S(Cys121) interaction at 2.74 A is clearly visible at higher pH. Despite the higher R(L) and stronger axial Cys121 interaction with Cu(II) ion, the Met121Cys variant remains largely a type 1 copper protein at low pH (with hyperfine coupling constant A( parallel) = 54 x 10(-4) cm(-1) at pH 4 and 5), or distorted type 1 or green copper protein at high pH (A(parallel) = 87 x 10(-4) cm(-1) at pH 8 and 9), attributable to the relatively long distance between the axial ligand and copper and the constraint placed by the protein scaffold. To shorten the distance between axial ligand and copper, we replaced Met121 with a nonproteinogenic amino acid homocysteine that contains an extra methylene group, resulting in a variant whose spectra (R(L)= 1.5, and A(parallel) = 180 x 10(-4) cm(-1)) and Cu-S(Cys) distance (2.22 A) are very similar to those of the red copper protein nitrosocyanin. Replacing Met121 with Cys or homocysteine resulted in lowering of the reduction potential from 222 mV in the native azurin to 95 +/- 3 mV for Met121Cys azurin and 113 +/- 6 mV for Met121Hcy azurin at pH 7. The results strongly support the "coupled distortion" model that helps explain axial ligand tuning of spectroscopic properties in cupredoxins, and demonstrate the power of using unnatural amino acids to address critical chemical biological questions.

Speeding up (13)C direct detection biomolecular NMR spectroscopy

After the exploitation of (1)H polarization as a starting source for (13)C direct detection experiments, pulse sequences are designed which exploit the accelerated (1)H longitudinal relaxation to expedite (13)C direct detection experiments. We show here that 2D experiments based on (13)C direct detection on a 0.5 mM water sample of ubiquitin can be recorded in a few minutes and 3D experiments in a few hours. We also show that fast methods like nonuniform sampling can be easily implemented. As overall experimental time has always been a counter indication for the use of (13)C direct detection experiments, this research opens new avenues for the application of (13)C NMR to biological molecules.