A multigeneration analysis of cytochrome b(562) redox variants: evolutionary strategies for modulating redox potential revealed using a library approach

The structure of the diheme cytochrome c4 from Neisseria gonorrhoeae reveals multiple contributors to tuning reduction potentials

Cytochrome c4 (c4) is a diheme protein implicated as an electron donor to cbb3 oxidases in multiple pathogenic bacteria. Despite its prevalence, understanding of how specific structural features of c4 optimize its function is lacking. The human pathogen Neisseria gonorrhoeae (Ng) thrives in low oxygen environments owing to the activity of its cbb3 oxidase. Herein, we report characterization of Ng c4. Spectroelectrochemistry experiments of the wild-type (WT) protein have shown that the two Met/His-ligated hemes differ in potentials by ∼100 mV, and studies of the two His/His-ligated variants provided unambiguous assignment of heme A from the N-terminal domain of the protein as the high-potential heme. The crystal structure of the WT protein at 2.45 Å resolution has revealed that the two hemes differ in their solvent accessibility. In particular, interactions made by residues His57 and Ser59 in Loop1 near the axial ligand Met63 contribute to the tight enclosure of heme A, working together with the surface charge, to raise the reduction potential of the heme iron in this domain. The structure reveals a prominent positively-charged patch, which encompasses surfaces of both domains. In contrast to prior findings with c4 from Pseudomonas stutzeri, the interdomain interface of Ng c4 contributes minimally to the values of the heme iron potentials in the two domains. Analyses of the heme solvent accessibility, interface properties, and surface charges offer insights into the interplay of these structural elements in tuning redox properties of c4 and other multiheme proteins.

Reaction of Thiosulfate Dehydrogenase with a Substrate Mimic Induces Dissociation of the Cysteine Heme Ligand Giving Insights into the Mechanism of Oxidative Catalysis

Thiosulfate dehydrogenases are bacterial cytochromes that contribute to the oxidation of inorganic sulfur. The active sites of these enzymes contain low-spin c-type heme with Cys^–/His axial ligation. However, the reduction potentials of these hemes are several hundred mV more negative than that of the thiosulfate/tetrathionate couple (E_m, +198 mV), making it difficult to rationalize the thiosulfate oxidizing capability. Here, we describe the reaction of Campylobacter jejuni thiosulfate dehydrogenase (TsdA) with sulfite, an analogue of thiosulfate. The reaction leads to stoichiometric conversion of the active site Cys to cysteinyl sulfonate (C_α-CH₂-S-SO₃^–) such that the protein exists in a form closely resembling a proposed intermediate in the pathway for thiosulfate oxidation that carries a cysteinyl thiosulfate (C_α-CH₂-S-SSO₃^–). The active site heme in the stable sulfonated protein displays an E_m approximately 200 mV more positive than the Cys^–/His-ligated state. This can explain the thiosulfate oxidizing activity of the enzyme and allows us to propose a catalytic mechanism for thiosulfate oxidation. Substrate-driven release of the Cys heme ligand allows that side chain to provide the site of substrate binding and redox transformation; the neighboring heme then simply provides a site for electron relay to an appropriate partner. This chemistry is distinct from that displayed by the Cys-ligated hemes found in gas-sensing hemoproteins and in enzymes such as the cytochromes P450. Thus, a further class of thiolate-ligated hemes is proposed, as exemplified by the TsdA centers that have evolved to catalyze the controlled redox transformations of inorganic oxo anions of sulfur.

IonoBiology: The functional dynamics of the intracellular metallome, with lessons from bacteria

Metal ions are essential for life and represent the second most abundant constituent (after water) of any living cell. While the biological importance of inorganic ions has been appreciated for over a century, we are far from a comprehensive understanding of the functional roles that ions play in cells and organisms. In particular, recent advances are challenging the traditional view that cells maintain constant levels of ion concentrations (ion homeostasis). In fact, the ionic composition (metallome) of cells appears to be purposefully dynamic. The scientific journey that started over 60 years ago with the seminal work by Hodgkin and Huxley on action potentials in neurons is far from reaching its end. New evidence is uncovering how changes in ionic composition regulate unexpected cellular functions and physiology, especially in bacteria, thereby hinting at the evolutionary origins of the dynamic metallome. It is an exciting time for this field of biology, which we discuss and refer to here as IonoBiology.

Heme ligation and redox chemistry in two bacterial thiosulfate dehydrogenase (TsdA) enzymes

Thiosulfate dehydrogenases (TsdAs) are bidirectional bacterial di-heme enzymes that catalyze the interconversion of tetrathionate and thiosulfate at measurable rates in both directions. In contrast to our knowledge of TsdA activities, information on the redox properties in the absence of substrates is rather scant. To address this deficit, we combined magnetic CD (MCD) spectroscopy and protein film electrochemistry (PFE) in a study to resolve heme ligation and redox chemistry in two representative TsdAs. We examined the TsdAs from Campylobacter jejuni, a microaerobic human pathogen, and from the purple sulfur bacterium Allochromatium vinosum. In these organisms, the enzyme functions as a tetrathionate reductase and a thiosulfate oxidase, respectively. The active site Heme 1 in both enzymes has His/Cys ligation in the ferric and ferrous states and the midpoint potentials (E_m) of the corresponding redox transformations are similar, −185 mV versus standard hydrogen electrode (SHE). However, fundamental differences are observed in the properties of the second, electron transferring, Heme 2. In C. jejuni, TsdA Heme 2 has His/Met ligation and an E_m of +172 mV. In A. vinosum TsdA, Heme 2 reduction triggers a switch from His/Lys ligation (E_m, −129 mV) to His/Met (E_m, +266 mV), but the rates of interconversion are such that His/Lys ligation would be retained during turnover. In summary, our findings have unambiguously assigned E_m values to defined axial ligand sets in TsdAs, specified the rates of Heme 2 ligand exchange in the A. vinosum enzyme, and provided information relevant to describing their catalytic mechanism(s).

De novo design of symmetric ferredoxins that shuttle electrons in vivo

A symmetric origin for bacterial ferredoxins was first proposed over 50 y ago, yet, to date, no functional symmetric molecule has been constructed. It is hypothesized that extant proteins have drifted from their symmetric roots via gene duplication followed by mutations. Phylogenetic analyses of extant ferredoxins support the independent evolution of N- and C-terminal sequences, thereby allowing consensus-based design of symmetric 4Fe-4S molecules. All designs bind two [4Fe-4S] clusters and exhibit strongly reducing midpoint potentials ranging from -405 to -515 mV. One of these constructs efficiently shuttles electrons through a designed metabolic pathway in Escherichia coli These finding establish that ferredoxins consisting of a symmetric core can be used as a platform to design novel electron transfer carriers for in vivo applications. Outer-shell asymmetry increases sequence space without compromising electron transfer functionality.

Redox tuning of cytochrome b562 through facile metal porphyrin substitution

The biologically and nanotechnologically important heme protein cytochrome b(562) was reconstructed with zinc and copper porphyrins, leading to significant changes in the spectral, redox and electron transfer properties. The Cu form shifts the redox potential by +300 mV and exhibits high electron transfer, while the Zn form is redox inert.

Single-molecule mapping of long-range electron transport for a cytochrome b(562) variant

Cytochrome b(562) was engineered to introduce a cysteine residue at a surface-exposed position to facilitate direct self-assembly on a Au(111) surface. The confined protein exhibited reversible and fast electron exchange with a gold substrate over a distance of 20 Å between the heme redox center and the gold surface, a clear indication that a long-range electron-transfer pathway is established. Electrochemical scanning tunneling microscopy was used to map electron transport features of the protein at the single-molecule level. Tunneling resonance was directly imaged and apparent molecular conductance was measured, which both show strong redox-gated effects. This study has addressed the first case of heme proteins and offered new perspectives in single-molecule bioelectronics.

AMD biofilms: using model communities to study microbial evolution and ecological complexity in nature

Similar to virtually all components of natural environments, microbial systems are inherently complex and dynamic. Advances in cultivation-independent molecular methods have provided a route to study microbial consortia in their natural surroundings and to begin resolving the community structure, dominant metabolic processes and inter-organism interactions. However, the utility of these methods generally scales inversely with community complexity. By applying genomics-enabled methods to the study of natural microbial communities with reduced levels of species richness, a relatively comprehensive understanding of the metabolic networks and evolutionary processes within these communities can be attained. In such well-defined model systems, it is also possible to link emergent ecological patterns to their molecular and evolutionary underpinnings, facilitating construction of predictive ecosystem models. In this study, we review over a decade of research on one such system-acid mine drainage biofilm communities. We discuss the value and limitations of tractable model microbial communities in developing molecular methods for microbial ecology and in uncovering principles that may explain behavior in more complex systems.

Analysis of the electrochemistry of hemes with E(m)s spanning 800 mV

The free energy of heme reduction in different proteins is found to vary over more than 18 kcal/mol. It is a challenge to determine how proteins manage to achieve this enormous range of E(m)s with a single type of redox cofactor. Proteins containing 141 unique hemes of a-, b-, and c-type, with bis-His, His-Met, and aquo-His ligation were calculated using Multi-Conformation Continuum Electrostatics (MCCE). The experimental E(m)s range over 800 mV from -350 mV in cytochrome c(3) to 450 mV in cytochrome c peroxidase (vs. SHE). The quantitative analysis of the factors that modulate heme electrochemistry includes the interactions of the heme with its ligands, the solvent, the protein backbone, and sidechains. MCCE calculated E(m)s are in good agreement with measured values. Using no free parameters the slope of the line comparing calculated and experimental E(m)s is 0.73 (R(2) = 0.90), showing the method accounts for 73% of the observed E(m) range. Adding a +160 mV correction to the His-Met c-type hemes yields a slope of 0.97 (R(2) = 0.93). With the correction 65% of the hemes have an absolute error smaller than 60 mV and 92% are within 120 mV. The overview of heme proteins with known structures and E(m)s shows both the lowest and highest potential hemes are c-type, whereas the b-type hemes are found in the middle E(m) range. In solution, bis-His ligation lowers the E(m) by approximately 205 mV relative to hemes with His-Met ligands. The bis-His, aquo-His, and His-Met ligated b-type hemes all cluster about E(m)s which are approximately 200 mV more positive in protein than in water. In contrast, the low potential bis-His c-type hemes are shifted little from in solution, whereas the high potential His-Met c-type hemes are raised by approximately 300 mV from solution. The analysis shows that no single type of interaction can be identified as the most important in setting heme electrochemistry in proteins. For example, the loss of solvation (reaction field) energy, which raises the E(m), has been suggested to be a major factor in tuning in situ E(m)s. However, the calculated solvation energy vs. experimental E(m) shows a slope of 0.2 and R(2) of 0.5 thus correlates weakly with E(m)s. All other individual interactions show even less correlation with E(m). However the sum of these terms does reproduce the range of observed E(m)s. Therefore, different proteins use different aspects of their structures to modulate the in situ heme electrochemistry. This study also shows that the calculated E(m)s are relatively insensitive to different heme partial charges and to the protein dielectric constant used in the simulation.

Design of a novel heme protein with a non-heme globin scaffold

A binding site for iron protoporphyrin IX (heme) was designed and embedded in a photosynthetic non-heme protein, phycocyanin, which forms a globin-like backbone structure, called a globin fold, but lacks sequence similarity to the globin family containing myoglobins and hemoglobins. Based on the structural alignment of the phycocyanin and myoglobin molecules, the proximal and distal His residues were repositioned in the phycocyanin sequence for heme ligation. The heme-binding pocket was created around the His residues by several residue replacements in the phycocyanin core. The synthesized phycocyanin variant, designated as HPY, bound one heme per protein molecule and showed spectroscopic features characteristic of six-coordinated heme proteins. The heme-binding HPY exhibited redox activity with an electrochemical midpoint potential of -130 mV against the standard hydrogen electrode, which was approximately 200 mV lower than the potential of natural myoglobins but 50 mV higher than the typical values of designed heme proteins with four-helix bundle or globin scaffolds. HPY also displayed native-like folding properties, in contrast to these designed heme proteins. However, the bound ferrous heme of HPY was quickly reoxidized by air and did not stably bind O(2), unlike the natural globins. The present results demonstrated that the globin fold of a non-globin protein is suitable for binding heme but is not sufficient for the reversible O(2) binding and myoglobin functions. The comparison of HPY with the natural globins may yield new insights into the essential features for realizing the natural heme protein functions.

Designing photosystem II: molecular engineering of photo-catalytic proteins

Biological photosynthesis utilizes membrane-bound pigment/protein complexes to convert light into chemical energy through a series of electron-transfer events. In the unique photosystem II (PSII) complex these electron-transfer events result in the oxidation of water to molecular oxygen. PSII is an extremely complex enzyme and in order to exploit its unique ability to convert sunlight into chemical energy it will be necessary to make a minimal model. Here we will briefly describe how PSII functions and identify those aspects that are essential in order to catalyze the oxidation of water into O(2), and review previous attempts to design simple photo-catalytic proteins and summarize our current research exploiting the E. coli bacterioferritin protein as a scaffold into which multiple cofactors can be bound, to oxidize a manganese metal center upon illumination. Through the reverse engineering of PSII and light driven water splitting reactions it may be possible to provide a blueprint for catalysts that can produce clean green fuel for human energy needs.

Mutagenic evidence for the optimal control of evolutionary dynamics

Elucidating the fitness measures optimized during the evolution of complex biological systems is a major challenge in evolutionary theory. We present experimental evidence and an analytical framework demonstrating how biochemical networks exploit optimal control strategies in their evolutionary dynamics. Optimal control theory explains a striking pattern of extremization in the redox potentials of electron transport proteins, assuming only that their fitness measure is a control objective functional with bounded controls.

Engineering model proteins for Photosystem II function

Our knowledge of Photosystem II and the molecular mechanism of oxygen production are rapidly advancing. The time is now ripe to exploit this knowledge and use it as a blueprint for the development of light-driven catalysts, ultimately for the splitting of water into O2 and H2. In this article, we outline the background and our approach to this technological application through the reverse engineering of Photosystem II into model proteins.

Modulation of heme redox potential in the cytochrome c6 family

Cytochrome c6A is a unique dithio-cytochrome of green algae and plants. It has a very similar core structure to that of bacterial and algal cytochromes c6 but is unable to fulfill the same function of transferring electrons from cytochrome f to photosystem I. A key feature is that its heme midpoint potential is more than 200 mV below that of cytochrome c6 despite having His and Met as axial heme-iron ligands. To identify the molecular origins of the difference in potential, the structure of cytochrome c6 from the cyanobacterium Phormidium laminosum has been determined by X-ray crystallography and compared with the known structure of cytochrome c6A. One salient difference of the heme pockets is that a highly conserved Gln (Q51) in cytochrome c6 is replaced by Val (V52) in c6A. Using protein film voltammetry, we found that swapping these residues raised the c6A potential by +109 mV and decreased that of c6 by almost the same extent, -100 mV. X-ray crystallography of the V52Q protein showed that the Gln residue adopts the same configuration relative to the heme as in cytochrome c6 and we propose that this stereochemistry destabilizes the oxidized form of the heme. Consequently, replacement of Gln by Val was probably a key step in the evolution of cytochrome c6A from cytochrome c6, inhibiting reduction by the cytochrome b6f complex and facilitating establishment of a new function.

Structure and redox properties of the haem centre in the C357M mutant of cytochrome P450cam

The effects of site-specific mutation of the axial cysteine (C357M) to a methionine residue in cytochrome P450cam on the enzyme's coordination geometry and redox potential have been investigated. The absorption spectra of the haem centre in the C357M mutant of the enzyme showed close similarity to those of cytochrome c both in the oxidised and reduced forms. A well-defined absorption peak at 695 nm, similar to that seen in the case of cytochrome c and characteristic of methionine ligation to the ferric haem, was observed. The results indicated that the haem of C357M cytochrome P450cam is possibly axially coordinated to a methionine and a histidine, analogously to cytochrome c. The circular dichroism spectra in the visible and the far-UV regions suggested that the tertiary structure of the haem cavity in the C357M mutant cytochrome P450cam was distinctly different from that in the wild-type enzyme or in cytochrome c, although the secondary structure of the mutant remained identical to that of the wild-type cytochrome P450cam. Comparison of the natures of the CD spectra in the 400 nm and 695 nm regions of the C357M mutant of cytochrome P450cam with those of horse cytochrome c suggested (R) chirality at the sulfur atom of the iron-bound methionine residue in the mutant. The redox potential of the haem centre, estimated by redox titration of the C357M mutant, was found to be +260 mV, which is much higher than that in the wild-type enzyme and similar to the redox potential of cytochrome c. This supported the concept that axial ligation of the haem plays the major role in tuning the redox potential of the haem centre in haem proteins.

Multielectron Redox Chemistry of Iron Porphyrinogens

Iron octamethylporphyrinogens were prepared and structurally characterized in three different oxidation states in the absence of axial ligands and with sodium or tetrafluoroborate as the only counterions. Under these conditions, the iron- and ligand-based redox chemistry of iron porphyrinogens can be defined. The iron center is easily oxidized by a single electron (E(1/2) = -0.57 V vs NHE in CH(3)CN) when confined within the fully reduced macrocycle. The porphyrinogen ligand also undergoes oxidation but in a single four-electron step (E(p) = +0.77 V vs NHE in CH(3)CN); one of the ligand-based electrons is intercepted for the reduction of Fe(III) to Fe(II) to result in an overall three-electron oxidation process. The oxidation equivalents in the macrocycle are stored in C(alpha)-C(alpha) bonds of spirocyclopropane rings, formed between adjacent pyrroles. EPR, magnetic and Mossbauer measurements, and DFT computations of the redox states of the iron porphyrinogens reveal that the reduced ligand gives rise to iron in intermediate spin states, whereas the fully oxidized ligand possesses a weaker sigma-donor framework, giving rise to high-spin iron. Taken together, the results reported herein establish a metal-macrocycle cooperativity that engenders a multielectron chemistry for iron porphyrinogens that is unavailable to heme cofactors.

Electric-field-induced redox potential shifts of tetraheme cytochromes c3 immobilized on self-assembled monolayers: surface-enhanced resonance Raman spectroscopy and simulation studies

The tetraheme protein cytochrome c(3) (Cyt-c(3)) from Desulfovibrio gigas, immobilized on a self-assembled monolayer (SAM) of 11-mercaptoundecanoic acid, is studied by theoretical and spectroscopic methods. Molecular dynamics simulations indicate that the protein docks to the negatively charged SAM via its lysine-rich domain around the exposed heme IV. Complex formation is associated with only little protein structural perturbations. This finding is in line with the resonance Raman and surface-enhanced resonance Raman (SERR) spectroscopic results that indicate essentially the same heme pocket structures for the protein in solution and adsorbed on SAM-coated Ag electrodes. Electron- and proton-binding equilibrium calculations reveal substantial negative shifts of the redox potentials compared to the protein in solution. The magnitude of these shifts decreases in the order heme IV (-161 mV) > heme III (-73 mV) > heme II (-57 mV) > heme I (-26 mV), resulting in a change of the order of reduction. These shifts originate from the distance-dependent electrostatic interactions between the SAM headgroups and the individual hemes, leading to a stabilization of the oxidized forms. The results of the potential-dependent SERR spectroscopic analyses are consistent with the theoretical predictions and afford redox potential shifts of -160 mV (heme IV), -90 mV (heme III), -70 mV (heme II), and +20 mV (heme I) relative to the experimental redox potentials for Cyt-c(3) in solution. SERR spectroscopic experiments reveal electric-field-induced changes of the redox potentials also for the structurally very similar Cyt-c(3) from Desulfovibrio vulgaris, although the shifts are somewhat smaller compared to Cyt-c(3) from D. gigas. This study suggests that electric-field-induced redox potential shifts may also occur upon binding to biomembranes or partner proteins and thus may affect biological electron transfer processes.

De novo proteins from designed combinatorial libraries

Combinatorial libraries of de novo amino acid sequences can provide a rich source of diversity for the discovery of novel proteins with interesting and important activities. Randomly generated sequences, however, rarely fold into well-ordered proteinlike structures. To enhance the quality of a library, features of rational design must be used to focus sequence diversity into those regions of sequence space that are most likely to yield folded structures. This review describes how focused libraries can be constructed by designing the binary pattern of polar and nonpolar amino acids to favor proteins that contain abundant secondary structure, while simultaneously burying hydrophobic side chains and exposing hydrophilic side chains to solvent. The "binary code" for protein design was used to construct several libraries of de novo proteins, including both alpha-helical and beta-sheet structures. The recently determined solution structure of a binary patterned four-helix bundle is well ordered, thereby demonstrating that sequences that have neither been selected by evolution (in vivo or in vitro) nor designed by computer can form nativelike proteins. Examples are presented demonstrating how binary patterned libraries have successfully produced well-ordered structures, cofactor binding, catalytic activity, self-assembled monolayers, amyloid-like nanofibrils, and protein-based biomaterials.

Distance dependence of interactions between charged centres in proteins with common structural features

Data collected for interactions among redox centres, and interactions between redox centres and acid-base residues in a family of small multihaem cytochromes are analysed. The distance dependent attenuation of the interactions between non-surface charges, for separations that range from 8 to 23 angstroms, can be described by a simple function derived from the Debye-Huckel formalism, fit to 9.5 and 7.6 as values for the relative dielectric constant and Debye length, respectively. However, there is considerable scatter in the data despite the structural similarities among the proteins, which is discussed in the framework of using such simple models in predicting properties of novel proteins.

Midpoint reduction potentials and heme binding stoichiometries of de novo proteins from designed combinatorial libraries

We previously reported the de novo design of combinatorial libraries of proteins targeted to fold into four-helix bundles. The sequences of these proteins were designed using a binary code strategy in which each position in the linear sequence is designated as either polar or nonpolar, but the exact identity of the amino acid at each position is varied combinatorially. We subsequently reported that approximately half of these binary coded proteins were capable of binding heme. These de novo heme-binding proteins showed CO binding characteristics similar to natural heme proteins, and several were active as peroxidases. Here we analyze the midpoint reduction potentials and heme binding stoichiometries of several of these de novo heme proteins. All the proteins bound heme with a 1:1 stoichiometry. The reduction potentials ranged from -112 to -176 mV. We suggest that this represents an estimate of the default range of potentials for heme proteins that have neither been prejudiced by rational design nor selected by evolution.

Effect of redox potential of the heme on the peroxidase activity of cytochrome b562

Measurements of peroxidase activities of two site-specific mutants and wild type cytochrome b562 suggest that the enzymatic activity correlates with the redox potential of the metal center. A lower value of the Fe(3+)/Fe(2+) redox potential seems to be important for promoting peroxidase activity of the hemeprotein possibly by stabilization of the high-valent redox intermediate involved in the catalytic function. The results provide an approach towards rational tuning of enzyme function when 'grafted' into a new protein environment.

Designing redox metalloproteins from bottom-up and top-down perspectives

The design of redox-active metalloproteins has been approached from two different directions. The de novo design approach has recently reached an important stage, at which structural information on several different designed metalloproteins has been obtained. This new information highlights the real challenge of this approach. The alternative approach involving re-engineering of evolved proteins has also made significant advances recently.

Tuning of the redox potential of the primary electron donor in reaction centres of purple bacteria: effects of amino acid polarity and position

Mutation of residues His L168 and Phe M197 in the reaction centre from Rhodobacter sphaeroides has an unusually strong effect on the mid-point redox potential (E(m)) of the pair of bacteriochlorophylls that form the primary donor of electrons, tuning E(m) over a range of nearly 250 mV. This effect is correlated to the accompanying change in the permanent dipole of the L168 or M197 residue, suggesting it is mediated by changes in charge-dipole interactions. Comparisons with mutations made at a variety of other positions show that this correlation is particular to this residue pair, perhaps reflecting their proximity to the ring I regions of the dimer bacteriochlorophylls that form the overlap region between these molecules.