Structural basis of electron transfer modulation in the purple CuA center

Capturing the Binuclear Copper State of Peptidylglycine Monooxygenase Using a Peptidyl-Homocysteine Lure

Peptidylglycine monooxygenase is a copper-dependent enzyme that catalyzes C-alpha hydroxylation of glycine extended pro-peptides, a critical post-translational step in peptide hormone processing. The canonical mechanism posits that dioxygen binds at the mononuclear M-center to generate a Cu(II)-superoxo species capable of H atom abstraction from the peptidyl substrate, followed by long-range electron tunneling from the CuH center. Recent crystallographic and biochemical data have challenged this mechanism, suggesting instead that an "open-to-closed" transition brings the copper centers closer, allowing reactivity within a binuclear intermediate. Here we present the first direct observation of an enzyme-bound binuclear copper species, captured by the use of an Ala-Ala-Phe-hCys inhibitor complex. This molecule reacts with the fully reduced enzyme to form a thiolate-bridged binuclear species characterized by EXAFS of the WT and its M314H variant and with the oxidized enzyme to form a novel mixed valence entity characterized by UV/vis and EPR. Mechanistic implications are discussed.

Direct Methane Oxidation by Copper- and Iron-Dependent Methane Monooxygenases

Methane is a potent greenhouse gas that contributes significantly to climate change and is primarily regulated in Nature by methanotrophic bacteria, which consume methane gas as their source of energy and carbon, first by oxidizing it to methanol. The direct oxidation of methane to methanol is a chemically difficult transformation, accomplished in methanotrophs by complex methane monooxygenase (MMO) enzyme systems. These enzymes use iron or copper metallocofactors and have been the subject of detailed investigation. While the structure, function, and active site architecture of the copper-dependent particulate methane monooxygenase (pMMO) have been investigated extensively, its putative quaternary interactions, regulation, requisite cofactors, and mechanism remain enigmatic. The iron-dependent soluble methane monooxygenase (sMMO) has been characterized biochemically, structurally, spectroscopically, and, for the most part, mechanistically. Here, we review the history of MMO research, focusing on recent developments and providing an outlook for future directions of the field. Engineered biological catalysis systems and bioinspired synthetic catalysts may continue to emerge along with a deeper understanding of the molecular mechanisms of biological methane oxidation. Harnessing the power of these enzymes will necessitate combined efforts in biochemistry, structural biology, inorganic chemistry, microbiology, computational biology, and engineering.

Biochemical characterization of Dimocarpus longan polyphenol oxidase provides insights into its catalytic efficiency

The "dragon-eye" fruits produced by the tropical longan tree are rich in nutrients and antioxidants. They suffer from post-harvest enzymatic browning, a process for which mainly the polyphenol oxidase (PPO) family of enzymes is responsible. In this study, two cDNAs encoding the PPO have been cloned from leaves of Dimocarpus longan (Dl), heterologously expressed in Escherichia coli and purified by affinity chromatography. The prepro-DlPPO1 contains two signal peptides at its N-terminal end that facilitate transportation of the protein into the chloroplast stroma and to the thylakoid lumen. Removal of the two signal peptides from prepro-DlPPO1 yields pro-DlPPO1. The prepro-DlPPO1 exhibited higher thermal tolerance than pro-DlPPO1 (unfolding at 65 °C vs. 40 °C), suggesting that the signal peptide may stabilize the fold of DlPPO1. DlPPO1 can be classified as a tyrosinase because it accepts both monophenolic and diphenolic substrates. The pro-DlPPO1 exhibited the highest specificity towards the natural diphenol (-)-epicatechin (k_cat/K_M of 800 ± 120 s^-1 mM^-1), which is higher than for 4-methylcatechol (590 ± 99 s^-1 mM^-1), pyrogallol (70 ± 9.7 s^-1 mM^-1) and caffeic acid (4.3 ± 0.72 s^-1 mM^-1). The kinetic efficiencies of prepro-DlPPO1 are 23, 36, 1.7 and 4.7-fold lower, respectively, than those observed with pro-DlPPO1 for the four aforementioned diphenolic substrates. Additionally, docking studies showed that (-)-epicatechin has a lower binding energy than any other investigated substrate. Both kinetic and in-silico studies strongly suggest that (-)-epicatechin is a good substrate of DlPPO1 and ascertain the affinity of PPOs towards specific flavonoid compounds.

Catalysis and Electron Transfer in De Novo Designed Metalloproteins

One of the hallmark advances in our understanding of metalloprotein function is showcased in our ability to design new, non-native, catalytically active protein scaffolds. This review highlights progress and milestone achievements in the field of de novo metalloprotein design focused on reports from the past decade with special emphasis on de novo designs couched within common subfields of bioinorganic study: heme binding proteins, monometal- and dimetal-containing catalytic sites, and metal-containing electron transfer sites. Within each subfield, we highlight several of what we have identified as significant and important contributions to either our understanding of that subfield or de novo metalloprotein design as a discipline. These reports are placed in context both historically and scientifically. General suggestions for future directions that we feel will be important to advance our understanding or accelerate discovery are discussed.

Designing Artificial Metalloenzymes by Tuning of the Environment beyond the Primary Coordination Sphere

Metalloenzymes catalyze a variety of reactions using a limited number of natural amino acids and metallocofactors. Therefore, the environment beyond the primary coordination sphere must play an important role in both conferring and tuning their phenomenal catalytic properties, enabling active sites with otherwise similar primary coordination environments to perform a diverse array of biological functions. However, since the interactions beyond the primary coordination sphere are numerous and weak, it has been difficult to pinpoint structural features responsible for the tuning of activities of native enzymes. Designing artificial metalloenzymes (ArMs) offers an excellent basis to elucidate the roles of these interactions and to further develop practical biological catalysts. In this review, we highlight how the secondary coordination spheres of ArMs influence metal binding and catalysis, with particular focus on the use of native protein scaffolds as templates for the design of ArMs by either rational design aided by computational modeling, directed evolution, or a combination of both approaches. In describing successes in designing heme, nonheme Fe, and Cu metalloenzymes, heteronuclear metalloenzymes containing heme, and those ArMs containing other metal centers (including those with non-native metal ions and metallocofactors), we have summarized insights gained on how careful controls of the interactions in the secondary coordination sphere, including hydrophobic and hydrogen bonding interactions, allow the generation and tuning of these respective systems to approach, rival, and, in a few cases, exceed those of native enzymes. We have also provided an outlook on the remaining challenges in the field and future directions that will allow for a deeper understanding of the secondary coordination sphere a deeper understanding of the secondary coordintion sphere to be gained, and in turn to guide the design of a broader and more efficient variety of ArMs.

Acquisition of ionic copper by the bacterial outer membrane protein OprC through a novel binding site

Copper, while toxic in excess, is an essential micronutrient in all kingdoms of life due to its essential role in the structure and function of many proteins. Proteins mediating ionic copper import have been characterised in detail for eukaryotes, but much less so for prokaryotes. In particular, it is still unclear whether and how gram-negative bacteria acquire ionic copper. Here, we show that Pseudomonas aeruginosa OprC is an outer membrane, TonB-dependent transporter that is conserved in many Proteobacteria and which mediates acquisition of both reduced and oxidised ionic copper via an unprecedented CxxxM-HxM metal binding site. Crystal structures of wild-type and mutant OprC variants with silver and copper suggest that acquisition of Cu(I) occurs via a surface-exposed “methionine track” leading towards the principal metal binding site. Together with whole-cell copper quantitation and quantitative proteomics in a murine lung infection model, our data identify OprC as an abundant component of bacterial copper biology that may enable copper acquisition under a wide range of conditions.

How do Gram-negative bacteria acquire copper? This study shows that the outer membrane protein OprC from Pseudomonas aeruginosa is abundant during infection and mediates highly selective acquisition of both copper redox states via an extracellular "methionine track" and an unprecedented near-irreversible binding site.

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.

Understanding and Modulating Metalloenzymes with Unnatural Amino Acids, Non-Native Metal Ions, and Non-Native Metallocofactors

Metalloproteins set the gold standard for performing important functions, including catalyzing demanding reactions under mild conditions. Designing artificial metalloenzymes (ArMs) to catalyze abiological reactions has been a major endeavor for many years, but most ArM activities are far below those of native enzymes, making them unsuitable for most pratical applications. A critical step to advance the field is to fundamentally understand what it takes to not only confer but also fine-tune ArM activities so they match those of native enzymes. Indeed, only once we can freely modulate ArM activity to rival (or surpass!) natural enzymes can the potential of ArMs be fully realized. A key to unlocking ArM potential is the observation that one metal primary coordination sphere can display a range of functions and levels of activity, leading to the realization that secondary coordination sphere (SCS) interactions are critically important. However, SCS interactions are numerous, long-range, and weak, making them very difficult to reproduce in ArMs. Furthermore, natural enzymes are tied to a small set of biologically available functional moieties from canonical amino acids and physiologically available metal ions and metallocofactors, severely limiting the chemical space available to probe and tune ArMs. In this Account, we summarize the use of unnatural amino acids (UAAs) and non-native metal ions and metallocofactors by our group and our collaborators to probe and modulate ArM functions. We incorporated isostructural UAAs in a type 1 copper (T1Cu) protein azurin to provide conclusive evidence that axial ligand hydrophobicity is a major determinant of T1Cu redunction potential ( E°'). Closely related work from other groups are also discussed. We also probed the role of protein backbone interactions that cannot be altered by standard mutagenesis by replacing the peptide bond with an ester linkage. We used insight gained from these studies to tune the E°' of azurin across the entire physiological range, the broadest range ever achieved in a single metalloprotein. Introducing UAA analogues of Tyr into ArM models of heme-copper oxidase (HCO) revealed a linear relationship between p K_a, E°', and activity. We also substituted non-native hemes and non-native metal ions for their native equivalents in these models to resolve several issues that were intractable in native HCOs and the closely related nitric oxide reductases, such as their roles in modulating substrate affinity, electron transfer rate, and activity. We incorporated abiological cofactors such as ferrocene and Mn(salen) into azurin and myoglobin, respectively, to stabilize these inorganic and organometallic compounds in water, confer abiological functions, tune their E°' and activity through SCS interactions, and show that the approach to metallocofactor anchoring and orientation can tune enantioselectivity and alter function. Replacing Cu in azurin with non-native Fe or Ni can impart novel activities, such as superoxide reduction and C-C bond formation. While progress was made, we have identified only a small fraction of the interactions that can be generally applied to ArMs to fine-tune their functions. Because SCS interactions are subtle and heavily interconnected, it has been difficult to characterize their effects quantitatively. It is vital to develop spectroscopic and computational techniques to detect and quantify their effects in both resting states and catalytic intermediates.

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.

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.

Characterization of a long overlooked copper protein from methane- and ammonia-oxidizing bacteria

Methane-oxidizing microbes catalyze the oxidation of the greenhouse gas methane using the copper-dependent enzyme particulate methane monooxygenase (pMMO). Isolated pMMO exhibits lower activity than whole cells, however, suggesting that additional components may be required. A pMMO homolog, ammonia monooxygenase (AMO), converts ammonia to hydroxylamine in ammonia-oxidizing bacteria (AOB) which produce another potent greenhouse gas, nitrous oxide. Here we show that PmoD, a protein encoded within many pmo operons that is homologous to the AmoD proteins encoded within AOB amo operons, forms a copper center that exhibits the features of a well-defined Cu_A site using a previously unobserved ligand set derived from a cupredoxin homodimer. PmoD is critical for copper-dependent growth on methane, and genetic analyses strongly support a role directly related to pMMO and AMO. These findings identify a copper-binding protein that may represent a missing link in the function of enzymes critical to the global carbon and nitrogen cycles.

Methane- and ammonia-oxidizing bacteria use the integral membrane, copper-dependent enzymes particulate methane monooxygenase (pMMO) and ammonia monooxygenase (AMO) to oxidize methane and ammonia. Here the authors structurally characterize the copper-binding protein PmoD, which contains an unusual CuA site and their genetic analyses strongly support a pMMO and AMO related function of PmoD.

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.

Probing the role of the backbone carbonyl interaction with the CuA center in azurin by replacing the peptide bond with an ester linkage

The role of a backbone carbonyl interaction with an engineered CuA center in azurin was investigated by developing a method of synthesis and incorporation of a depsipeptide where one of the amide bonds in azurin is replaced by an ester bond using expressed protein ligation. Studies by electronic absorption and electron paramagnetic resonance spectroscopic techniques indicate that, while the substitution does not significantly alter the geometry of the site, it weakens the axial interaction to the CuA center and strengthens the Cu-Cu bond, as evidenced by the blue shift of the near-IR absorption that has been assigned to the Cu-Cu ψ → ψ* transition. Interestingly, the changes in the electronic structure from the replacement did not result in a change in the reduction potential of the CuA center, suggesting that the diamond core structure of Cu2SCys2 is resistant to variations in axial interactions.

Binuclear Cu(A) Formation in Biosynthetic Models of Cu(A) in Azurin Proceeds via a Novel Cu(Cys)2His Mononuclear Copper Intermediate

Cu(A) is a binuclear electron transfer (ET) center found in cytochrome c oxidases (CcOs), nitrous oxide reductases (N₂ORs), and nitric oxide reductase (NOR). In these proteins, the Cu(A) centers facilitate efficient ET (kET > 10⁴s⁻¹) under low thermodynamic driving forces (10-90 mV). While the structure and functional properties of Cu(A) are well understood, a detailed mechanism of the incorporation of copper into the protein and the identity of the intermediates formed during the Cu(A) maturation process are still lacking. Previous studies of the Cu(A) assembly mechanism in vitro using a biosynthetic model Cu(A) center in azurin (Cu(A)Az) identified a novel intermediate X (Ix) during reconstitution of the binuclear site. However, because of the instability of Ix and the coexistence of other Cu centers, such as Cu(A)' and type 1 copper centers, the identity of this intermediate could not be established. Here, we report the mechanism of Cu(A) assembly using variants of Glu114XCuAAz (X = Gly, Ala, Leu, or Gln), the backbone carbonyl of which acts as a ligand to the Cu(A) site, with a major focus on characterization of the novel intermediate Ix. We show that Cu(A) assembly in these variants proceeds through several types of Cu centers, such as mononuclear red type 2 Cu, the novel intermediate Ix, and blue type 1 Cu. Our results show that the backbone flexibility of the Glu114 residue is an important factor in determining the rates of T2Cu → Ix formation, suggesting that Cu(A) formation is facilitated by swinging of the ligand loop, which internalizes the T2Cu capture complex to the protein interior. The kinetic data further suggest that the nature of the Glu114 side chain influences the time scales on which these intermediates are formed, the wavelengths of the absorption peaks, and how cleanly one intermediate is converted to another. Through careful understanding of these mechanisms and optimization of the conditions, we have obtained Ix in ∼80-85% population in these variants, which allowed us to employ ultraviolet-visible, electron paramagnetic resonance, and extended X-ray absorption fine structure spectroscopic techniques to identify the Ix as a mononuclear Cu(Cys)(2)(His) complex. Because some of the intermediates have been proposed to be involved in the assembly of native Cu(A), these results shed light on the structural features of the important intermediates and mechanism of Cu(A) formation.

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.

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.

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.

Computational approaches for rational design of proteins with novel functionalities

Proteins are the most multifaceted macromolecules in living systems and have various important functions, including structural, catalytic, sensory, and regulatory functions. Rational design of enzymes is a great challenge to our understanding of protein structure and physical chemistry and has numerous potential applications. Protein design algorithms have been applied to design or engineer proteins that fold, fold faster, catalyze, catalyze faster, signal, and adopt preferred conformational states. The field of de novo protein design, although only a few decades old, is beginning to produce exciting results. Developments in this field are already having a significant impact on biotechnology and chemical biology. The application of powerful computational methods for functional protein designing has recently succeeded at engineering target activities. Here, we review recently reported de novo functional proteins that were developed using various protein design approaches, including rational design, computational optimization, and selection from combinatorial libraries, highlighting recent advances and successes.

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.

Design of functional metalloproteins

Metalloproteins catalyse some of the most complex and important processes in nature, such as photosynthesis and water oxidation. An ultimate test of our knowledge of how metalloproteins work is to design new metalloproteins. Doing so not only can reveal hidden structural features that may be missing from studies of native metalloproteins and their variants, but also can result in new metalloenzymes for biotechnological and pharmaceutical applications. Although it is much more challenging to design metalloproteins than non-metalloproteins, much progress has been made in this area, particularly in functional design, owing to recent advances in areas such as computational and structural biology.

Continuous and smooth potential energy surface for conductorlike screening solvation model using fixed points with variable areas

Rigorously continuous and smooth potential energy surfaces, as well as exact analytic gradients, are obtained for a conductorlike screening solvation model (CPCM, a variant of the general COSMO) with Hartree-Fock (RHF, ROHF, UHF, and MCSCF) and density functional theory (R-DFT, RO-DFT, and U-DFT) methods using a new tessellation scheme, fixed points with variable areas (FIXPVA). In FIXPVA, spheres centered at atoms are used to define the molecular cavity and surface. The surface of each sphere is divided into 60, 240, or 960 tesserae, which have positions fixed relative to the sphere center and areas scaled by switching functions of their distances to neighboring spheres. Analytic derivatives of the positions and areas of the surface tesserae with respect to atomic coordinates can be obtained and used to evaluate the solvation energy gradients. Due to the accurate analytic gradients and smooth potential energy surface, geometry optimization processes using these methods are stable and convergent.

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

The electronic properties of Thermus thermophilus Cu(A) in the oxidized form were studied by (1)H and (13)C NMR spectroscopy. All of the (1)H and (13)C resonances from cysteine and imidazole ligands were observed and assigned in a sequence-specific fashion. The detection of net electron spin density on a peptide moiety is attributed to the presence of a H-bond to a coordinating sulfur atom. This hydrogen bond is conserved in all natural Cu(A) variants and plays an important role for maintaining the electronic structure of the metal site, rendering the two Cys ligands nonequivalent. The anomalous temperature dependence of the chemical shifts is explained by the presence of a low-lying excited state located about 600 cm(-1) above the ground state. The room-temperature shifts can be described as the thermal average of a sigma(u)* ground state and a pi(u) excited state. These results provide a detailed description of the electronic structure of the Cu(A) site at atomic resolution in solution at physiologically relevant temperature.

Perturbations to the geometric and electronic structure of the CuA site: factors that influence delocalization and their contributions to electron transfer

Using a combination of electronic spectroscopies and DFT calculations, the effect of pH perturbation on the geometric and electronic structure of the CuA site has been defined. Descriptions are developed for high pH (pH = 7) and low pH (pH = 4) forms of CuA azurin and its H120A mutant which address the discrepancies concerning the extent of delocalization indicated by multifrequency EPR and ENDOR data (J. Am. Chem. Soc. 2005, 127, 7274; Biophys. J. 2002, 82, 2758). Our resonance Raman and MCD spectra demonstrate that the low pH and H120A mutant forms are essentially identical and are the perturbed forms of the completely delocalized high pH CuA site. However, in going from high pH to low pH, a seven-line hyperfine coupling pattern associated with complete delocalization of the electron (S = 1/2) over two Cu coppers (I(Cu) = 3/2) changes into a four-line pattern reflecting apparent localization. DFT calculations show that the unpaired electron is delocalized in the low pH form and reveal that its four-line hyperfine pattern results from the large EPR spectral effects of approximately 1% 4s orbital contribution of one Cu to the ground-state spin wave function upon protonative loss of its His ligand. The contribution of the Cu-Cu interaction to electron delocalization in this low symmetry protein site is evaluated, and the possible functional significance of the pH-dependent transition in regulating proton-coupled electron transfer in cytochrome c oxidase is discussed.

The role of ligand-containing loops at copper sites in proteins

Many approaches are being used to engineer metalloproteins, with most of these informed by, and aiming to further elucidate, the basic structural requirements for biological metal centers. Cupredoxins are type 1 (T1) copper-containing electron transfer (ET) proteins with a -barrel fold that is thought to constrain metal site structure. The T1 copper ion is bound by ligands mainly originating from a single active site loop whose length and structure varies. This Highlight article will focus on protein engineering studies which have investigated the role of the metal-binding loop for active site integrity and functionality. Scaffold differences are present within the cupredoxin family and their influence has also been assessed. Given the widespread occurrence of -barrel domains in nature, and the array of metal sites in proteins composed of loop regions, the studies described on this model system have implications for a variety of metalloproteins.

Reorganization Energy of the CuA Center in Purple Azurin: Impact of the Mixed Valence-to-Trapped Valence State Transition

Mixed valence (MV) coordination compounds play important roles in redox reactions in chemistry and biology. Details of the contribution of a mixed valence state to protein electron transfer (ET) reactivity such as reorganization energy, however, have not been experimentally defined. Herein we report measurements of reorganization energies of a binuclear CuA center engineered into Pseudomonas aeruginosa azurin that exhibits a reversible transition between a totally delocalized MV state at pH 8.0 and a trapped valence (TV) state at pH 4.0. The reorganization energy of a His120Ala variant of CuA azurin that displays a TV state at both the above pH values has also been determined. We found that the MV-to-TV state transition increases the reorganization energy by 0.18 eV, providing evidence that the MV state of the CuA center has lower reorganization energy than its TV counterpart. We have also shown that lowering the pH from 8.0 to 4.0 results in a similar (approximately 0.4 eV) decrease in reorganization energy for both blue (type 1) and purple (CuA) azurins, even though the reorganization energies of the two different copper centers are different at a given pH. These results suggest that the MV state plays only a secondary role in modulation of the ET reactivity via the reorganization energy, as compared to that of the driving force.

Respiratory transformation of nitrous oxide (N2O) to dinitrogen by Bacteria and Archaea

N2O is a potent greenhouse gas and stratospheric reactant that has been steadily on the rise since the beginning of industrialization. It is an obligatory inorganic metabolite of denitrifying bacteria, and some production of N2O is also found in nitrifying and methanotrophic bacteria. We focus this review on the respiratory aspect of N2O transformation catalysed by the multicopper enzyme nitrous oxide reductase (N2OR) that provides the bacterial cell with an electron sink for anaerobic growth. Two types of Cu centres discovered in N2OR were both novel structures among the Cu proteins: the mixed-valent dinuclear Cu(A) species at the electron entry site of the enzyme, and the tetranuclear Cu(Z) centre as the first catalytically active Cu-sulfur complex known. Several accessory proteins function as Cu chaperone and ABC transporter systems for the biogenesis of the catalytic centre. We describe here the paradigm of Z-type N2OR, whose characteristics have been studied in most detail in the genera Pseudomonas and Paracoccus. Sequenced bacterial genomes now provide an invaluable additional source of information. New strains harbouring nos genes and capability of N2O utilization are being uncovered. This reveals previously unknown relationships and allows pattern recognition and predictions. The core nos genes, nosZDFYL, share a common phylogeny. Most principal taxonomic lineages follow the same biochemical and genetic pattern and share the Z-type enzyme. A modified N2OR is found in Wolinella succinogenes, and circumstantial evidence also indicates for certain Archaea another type of N2OR. The current picture supports the view of evolution of N2O respiration prior to the separation of the domains Bacteria and Archaea. Lateral nos gene transfer from an epsilon-proteobacterium as donor is suggested for Magnetospirillum magnetotacticum and Dechloromonas aromatica. In a few cases, nos gene clusters are plasmid borne. Inorganic N2O metabolism is associated with a diversity of physiological traits and biochemically challenging metabolic modes or habitats, including halorespiration, diazotrophy, symbiosis, pathogenicity, psychrophily, thermophily, extreme halophily and the marine habitat down to the greatest depth. Components for N2O respiration cover topologically the periplasm and the inner and outer membranes. The Sec and Tat translocons share the task of exporting Nos components to their functional sites. Electron donation to N2OR follows pathways with modifications depending on the host organism. A short chronology of the field is also presented.

Engineered lanthanide-binding metallohomeodomains: designing folded chimeras by modular turn substitution

A series of chimeric metallohomeodomains are described, engineered by rational design of a flexible Ca/Ln binding site into a DNA-binding scaffold. A modular turn-substitution approach was used to create proteins that both bind DNA and lanthanide ions, while retaining the secondary structure of the full homeodomain (determined by circular dichroism [CD]). Four similar metallohomeodomains were designed (C1-C4), their structural stability predicted by molecular dynamics (MD) simulation of loop-mutations into the known homeodomain structure, and each designed protein cloned, expressed, and purified using standard molecular biology techniques. Two of the four loop insertions resulted in folded, metal- and DNA-binding proteins (EuC2 Kd = 2.1 +/- 0.4 microM; EuC4 Kd = 3.2 +/- 1.0 microM). These results show the successful incorporation of a metal site into a full protein domain, without compromising long-range structure. This is an important achievement in biomolecular design, as it provides a critical starting point for exploring metallonuclease function and substrate accessibility in a well-organized chimeric protein domain (rather than only in small HTH peptide systems).

Biosynthetic Inorganic Chemistry

Inorganic chemistry and biology can benefit greatly from each other. Although synthetic and physical inorganic chemistry have been greatly successful in clarifying the role of metal ions in biological systems, the time may now be right to utilize biological systems to advance coordination chemistry. One such example is the use of small, stable, easy-to-make, and well-characterized proteins as ligands to synthesize novel inorganic compounds. This biosynthetic inorganic chemistry is possible thanks to a number of developments in biology. This review summarizes the progress in the synthesis of close models of complex metalloproteins, followed by a description of recent advances in using the approach for making novel compounds that are unprecedented in either inorganic chemistry or biology. The focus is mainly on synthetic "tricks" learned from biology, as well as novel structures and insights obtained. The advantages and disadvantages of this biosynthetic approach are discussed.

Spectroscopic characterizations of bridging cysteine ligand variants of an engineered Cu2(Scys)2 CuA azurin

Bridging cysteine ligands of the Cu(A) center in an engineered Cu(A) azurin were replaced with serine, and the variants (Cys116Ser and Cys112Ser Cu(A) azurin) were characterized by mass spectrometry, as well as UV-vis and electron paramagnetic resonance (EPR) spectroscopic techniques. The replacements resulted in dramatically perturbed spectroscopic properties, indicating that the cysteines play a critical role in maintaining the structural integrity of the Cu center. The replacements at different cysteine residues resulted in different perturbations, even though the two cysteines are geometrically symmetrical in the primary coordination sphere with respect to the two copper ions. The Cys112Ser variant contains two distinct type 2 copper centers, while the Cys116Ser variant has one type 1 copper center with slight tetragonal distortion. Both the UV-vis and EPR spectra of the Cys116Ser variant change with pH, and the pK(a) of the transition is 6.0. A type 1 copper EPR spectrum with A(||) = 26 G was obtained at pH 7.0, while a type 2 copper EPR spectrum with A(||) = 140 G was found at pH 5.0. Interestingly, lowering the temperature from 290 to 85 K resulted in conversion of the Cys116Ser variant from a type 1 copper center to a type 2 copper center, suggesting rearrangement of the ligand around the copper or binding of an exogenous ligand at low temperature. This difference in mutation effects at different cysteines may be due to different constraints exerted on the two cysteines by hydrogen-bonding patterns in the ligand loop.

Axial methionine has much less influence on reduction potentials in a CuA center than in a blue copper center

The role of the highly conserved axial methionine of the purple CuA center in an engineered CuA azurin on modulating the reduction potentials of the copper center was investigated by a systematic replacement of the methionine with glutamate, aspartate, and leucine. In contrast to the same substitutions in the structurally related blue copper azurin, much smaller changes in reduction potential were observed in the CuA azurin upon replacing the methionine ligand with negatively charged Glu (-8 mV) and Asp (-5 mV) and more hydrophobic Leu (+16 mV). These findings are important in understanding the different roles of the two cupredoxins. The diamond core Cu2S2(Cys) structure of the CuA is much more resistant to variations of axial ligand interactions than the distorted tetrahedral structure of the blue copper protein. This difference may translate into a much wider range of reduction potentials (>1000 mV) for blue copper proteins that transfer electrons to a variety of partners in many different biological systems and a much narrower range of reduction potentials (<40 mV) for CuA proteins where a small difference in reduction potentials between the CuA and its redox partners is required.

pH-dependent transition between delocalized and trapped valence states of a CuA center and its possible role in proton-coupled electron transfer

A pH-dependent transition between delocalized and trapped mixed valence states of an engineered CuA center in azurin has been investigated by UV-visible absorption and electron paramagnetic resonance spectroscopic techniques. At pH 7.0, the CuA azurin displays a typical delocalized mixed valence dinuclear [Cu(1.5)....Cu(1.5)] spectra with optical absorptions at 485, 530, and 760 nm, and with a seven-line EPR hyperfine. Upon lowering of the pH from 7.0 to 4.0, the absorption at 760 nm shifted to lower energy toward 810 nm, and a four-line EPR hyperfine, typical of a trapped valence, was observed. The pH-dependent transition is reversible because increasing the pH restores all delocalized spectral features. Lowering the pH resulted in not only a trapped valence state, but also a dramatically increased reduction potential of the Cu center (from 160 mV to 340 mV). Mutation of the titratable residues around the metal-binding site ruled out Glu-114 and identified the C-terminal histidine ligand (His-120) as a site of protonation, because the His120Ala mutation abolished the above pH-dependent transition. The corresponding histidine in cytochrome c oxidases is along a major electron transfer pathway from CuA center to heme a. Because the protonation of this histidine can result in an increased reduction potential that will prevent electron flow from the CuA to heme a, the CuA and the histidine may play an important role in regulating proton-coupled electron transfer.

Spin Distribution and the Location of Protons in Paramagnetic Proteins

▪ Abstract  Two current frontiers in EPR research are high-field ([Formula: see text]) electron paramagnetic resonance (EPR) and high-field electron-nuclear double resonance (ENDOR). This review focuses on recent advances in high-field ENDOR and its applications to the study of proteins containing native paramagnetic sites. It concentrates on two aspects; the first concerns the determination of the location of protons and is related to the site geometry, and the second focuses on the spin density distribution within the site, which is inherent to the electronic structure. Both spin density and proton locations can be derived from ligand hyperfine couplings determined by ENDOR measurements. A brief description of the experimental methods is presented along with a discussion of the advantages and disadvantages of high-field ENDOR compared with conventional X-band (∼ 9.5 GHz) experiments. Specific examples of both protein single crystals and frozen solutions are then presented. These include the determination of the coordinates of water ligand protons in the Mn(II) site of concanavalin A, the detection of hydrogen bonds in a quinone radical in the bacterial photosynthetic reaction center as well as in the tyrosyl radical in ribonuclease reductase, and the study of the spin distribution in copper proteins. The copper proteins discussed are the type I copper of azurin and the binuclear CuA center in a number of proteins. The last part of the review presents a brief discussion of the interpretation of hyperfine couplings using quantum chemical calculations, primarily density functional theory (DFT) methods. Such methods are becoming an integral part of the data analysis tools, as they can facilitate signal assignment and provide the ultimate relation between the experimental hyperfine couplings and the electronic wave function.

Spectroscopic evidence for interactions between hexacyanoiron(II/III) and an engineered purple CuA azurin

Interactions between hexacyanoiron(II/III) and a dinuclear, mixed valence Cu(A) center in engineered Cu(A) azurin have been investigated by UV-visible (UV-vis) and electron paramagnetic resonance (EPR) spectroscopic techniques. Addition of ferricyanide (hexacyanoiron(III)) to the Cu(A) azurin resulted in a new absorption band around 500 nm in the UV-vis and an isotropic line at g = 2.16 in the EPR spectra. Control experiments, including additions of Cu(II)SO(4) or Cu(I)(CH(3)CN)(4)PF(6) to ferricyanide or ferrocyanide, as well as gel filtration purification of the ferricyanide-Cu(A) azurin adduct indicate complex formation between cupric ion and ferrocyanide ion in the protein. Solvent or small molecule accessibility, metal oxidation state and the presence of more than one metal ion are potential factors important for the complex formation. These findings must be taken into consideration when using ferricyanide or ferrocyanide as redox agents for studying Cu(A) centers in proteins.

Determination of reduction potential of an engineered CuA azurin by cyclic voltammetry and spectrochemical titrations

The reduction potentials of an engineered CuA azurin in its native and thermally denatured states have been determined using cyclic voltammetry and spectrochemical titrations. Using a 4,4'-dipyridyl disulfide modified gold electrode, the reduction potentials of native and thermally denatured CuA azurin are the same within the experimental error (422 +/- 5 and 425 +/- 5 mV vs. NHE, respectively, in 50 mM ammonium acetate buffer, pH 5.1, 300 mM NaCl, 25 degrees C), indicating that the potential is that of a nonnative state. In contrast, using a didodecyldimethylammonium bromide (DDAB) film-pyrolytic graphite edge (PGE) electrode, the reduction potentials of native and thermally denatured CuA azurin have been determined to be 271 +/- 7 mV (50 mM ammonium acetate buffer, pH 5.1, 4 degrees C) and 420 +/- 1 mV (50 mM ammonium acetate buffer, pH 5.1, 25 degrees C), respectively. Spectroscopic redox titration using [Ru(NH3)5Py]2+ resulted in a reduction potential (254+/-4 mV) (50 mM ammonium acetate buffer, pH 5.1, 4 degrees C) similar to the value obtained using the DDAB film-PGE electrochemical method. Complete reoxidation of [Ru(NH3)5Py]2+-reduced CuA azurin is also consistent with the conclusion that this spectrochemical titration method using [Ru(NH3)5Py]2+ measures the reduction potential of native CuA azurin.