Evaluation of Electron-Withdrawing Group Effects on Heme Binding in Designed Proteins: Implications for Heme a in Cytochrome c Oxidase

Finding the E-channel proton loading sites by calculating the ensemble of protonation microstates

The aerobic electron transfer chain builds a proton gradient by proton coupled electron transfer reactions through a series of proteins. Complex I is the first enzyme in the sequence. Here transfer of two electrons from NADH to quinone yields four protons pumped from the membrane N- (negative, higher pH) side to the P- (positive, lower pH) side. Protons move through three linear antiporter paths, with a few amino acids and waters providing the route; and through the E-channel, a complex of competing paths, with clusters of interconnected protonatable residues. Proton loading sites (PLS) transiently bind protons as they are transported from N- to P-compartments. PLS can be individual residues or extended clusters of residues. The program MCCE uses Monte Carlos sampling to analyze the E-channel proton binding in equilibrium with individual Molecular Dynamics snapshots from trajectories of Thermus thermuphillus Complex I in the apo, quinone and quinol bound states. At pH 7, the five E-channel subunits (Nqo4, Nqo7, Nqo8, Nqo10, and Nqo11) take >25,000 protonation microstates, each with different residues protonated. The microstate explosion is tamed by analyzing interconnected clusters of residues along the proton transfer paths. A proton is bound and released from a cluster of five coupled residues on the protein N-side and to six coupled residues in the protein center. Loaded microstates bind protons to sites closer to the P-side in the forward pumping direction. MCCE microstate analysis identifies strongly coupled proton binding amongst individual residues in the two PLS clusters.

Biosynthesis and trafficking of heme o and heme a: new structural insights and their implications for reaction mechanisms and prenylated heme transfer

Aerobic respiration is a key energy-producing pathway in many prokaryotes and virtually all eukaryotes. The final step of aerobic respiration is most commonly catalyzed by heme-copper oxidases embedded in the cytoplasmic or mitochondrial membrane. The majority of these terminal oxidases contain a prenylated heme (typically heme a or occasionally heme o) in the active site. In addition, many heme-copper oxidases, including mitochondrial cytochrome c oxidases, possess a second heme a cofactor. Despite the critical role of heme a in the electron transport chain, the details of the mechanism by which heme b, the prototypical cellular heme, is converted to heme o and then to heme a remain poorly understood. Recent structural investigations, however, have helped clarify some elements of heme a biosynthesis. In this review, we discuss the insight gained from these advances. In particular, we present a new structural model of heme o synthase (HOS) based on distance restraints from inferred coevolutionary relationships and refined by molecular dynamics simulations that are in good agreement with the experimentally determined structures of HOS homologs. We also analyze the two structures of heme a synthase (HAS) that have recently been solved by other groups. For both HOS and HAS, we discuss the proposed catalytic mechanisms and highlight how new insights into the heme-binding site locations shed light on previously obtained biochemical data. Finally, we explore the implications of the new structural data in the broader context of heme trafficking in the heme a biosynthetic pathway and heme-copper oxidase assembly.

Manganese and Cobalt in the Nonheme-Metal-Binding Site of a Biosynthetic Model of Heme-Copper Oxidase Superfamily Confer Oxidase Activity through Redox-Inactive Mechanism

The presence of a nonheme metal, such as copper and iron, in the heme-copper oxidase (HCO) superfamily is critical to the enzymatic activity of reducing O2 to H2O, but the exact mechanism the nonheme metal ion uses to confer and fine-tune the activity remains to be understood. We herein report that manganese and cobalt can bind to the same nonheme site and confer HCO activity in a heme-nonheme biosynthetic model in myoglobin. While the initial rates of O2 reduction by the Mn, Fe, and Co derivatives are similar, the percentages of reactive oxygen species (ROS) formation are 7%, 4%, and 1% and the total turnovers are 5.1 ± 1.1, 13.4 ± 0.7, and 82.5 ± 2.5, respectively. These results correlate with the trends of nonheme-metal-binding dissociation constants (35, 22, and 9 μM) closely, suggesting that tighter metal binding can prevent ROS release from the active site, lessen damage to the protein, and produce higher total turnover numbers. Detailed spectroscopic, electrochemical, and computational studies found no evidence of redox cycling of manganese or cobalt in the enzymatic reactions and suggest that structural and electronic effects related to the presence of different nonheme metals lead to the observed differences in reactivity. This study of the roles of nonheme metal ions beyond the Cu and Fe found in native enzymes has provided deeper insights into nature's choice of metal ion and reaction mechanism and allows for finer control of the enzymatic activity, which is a basis for the design of efficient catalysts for the oxygen reduction reaction in fuel cells.

Equilibrium Studies of Designed Metalloproteins

Complete thermodynamic descriptions of the interactions of cofactors with proteins via equilibrium studies are challenging, but are essential to the evaluation of designed metalloproteins. While decades of studies on protein-protein interaction thermodynamics provide a strong underpinning to the successful computational design of novel protein folds and de novo proteins with enzymatic activity, the corresponding paucity of data on metal-protein interaction thermodynamics limits the success of computational metalloprotein design efforts. By evaluating the thermodynamics of metal-protein interactions via equilibrium binding studies, protein unfolding free energy determinations, proton competition equilibria, and electrochemistry, a more robust basis for the computational design of metalloproteins may be provided. Our laboratory has shown that such studies provide detailed insight into the assembly and stability of designed metalloproteins, allow for parsing apart the free energy contributions of metal-ligand interactions from those of porphyrin-protein interactions in hemeproteins, and even reveal their mechanisms of proton-coupled electron transfer. Here, we highlight studies that reveal the complex interplay between the various equilibria that underlie metalloprotein assembly and stability and the utility of making these detailed measurements.

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

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

Proton coupled electron transfer and redox active tyrosines in Photosystem II

In this article, progress in understanding proton coupled electron transfer (PCET) in Photosystem II is reviewed. Changes in acidity/basicity may accompany oxidation/reduction reactions in biological catalysis. Alterations in the proton transfer pathway can then be used to alter the rates of the electron transfer reactions. Studies of the bioenergetic complexes have played a central role in advancing our understanding of PCET. Because oxidation of the tyrosine results in deprotonation of the phenolic oxygen, redox active tyrosines are involved in PCET reactions in several enzymes. This review focuses on PCET involving the redox active tyrosines in Photosystem II. Photosystem II catalyzes the light-driven oxidation of water and reduction of plastoquinone. Photosystem II provides a paradigm for the study of redox active tyrosines, because this photosynthetic reaction center contains two tyrosines with different roles in catalysis. The tyrosines, YZ and YD, exhibit differences in kinetics and midpoint potentials, and these differences may be due to noncovalent interactions with the protein environment. Here, studies of YD and YZ and relevant model compounds are described.

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.

Ligand preference and orientation in b- and c-type heme-binding proteins

Hemes are often incorporated into designed proteins. The importance of the heme ligand type and its orientation is still a matter of debate. Here, heme ligands and ligand orientation were investigated using a nonredundant (87 structures) and a redundant (1503 structures) set of structures to compare and contrast design features of natural b- and c-type heme-binding proteins. Histidine is the most common ligand. Marked differences in ligation motifs between b- and c-type hemes are higher occurrence of His-Met in c-type heme binding motifs (16.4% vs. 1.4%) and higher occurrence of exchangeable, small molecules in b-type heme binding motifs (67.6% vs. 9.9%). Histidine ligands that are part of the c-type CXXCH heme-binding motif show a distinct asymmetric distribution of orientation. They tend to point between either the heme propionates or between the NA and NB heme nitrogens. Molecular mechanics calculations show that this asymmetry is due to the bonded constraints of the covalent attachment between the heme and the protein. In contrast, the orientations of b-type hemes histidine ligands are found evenly distributed with no preference. Observed histidine heme ligand orientations show no dominating influence of electrostatic interactions between the heme propionates and the ligands. Furthermore, ligands in bis-His hemes are found more frequently perpendicular rather than parallel to each other. These correlations support energetic constraints on ligands that can be used in designing proteins.

Development of a heme protein structure-electrochemical function database

Proteins containing heme, iron(protoporphyrin IX) and its variants, continue to be one of the most-studied classes of biomolecules due to their diverse range of biological functions. The literature is abundant with reports of structural and functional characterization of individual heme proteins which demonstrate that heme protein reduction potential values, E_m, span the range from –550 mV to +450 mV versus SHE. In order to unite these data for the purposes of global analysis, a new web-based resource of heme protein structure–function relationships is presented: the Heme Protein Database (HPD). This database is the first of its kind to combine heme protein structural classifications including protein fold, heme type and heme axial ligands, with heme protein reduction potential values in a web-searchable format. The HPD is located at http://heme.chem.columbia.edu/heme.php. The data illustrate that heme protein E_m values are modulated over a 300 mV range by the type of global protein fold, a 600 mV range by the type of porphyrin and an 800 mV range by the axial ligands. Thus, the 1 V range observed in heme protein reduction potential values in biological systems arises from subtle combinations of these various factors.

Proton-coupled electron transfer in a biomimetic peptide as a model of enzyme regulatory mechanisms

Proton-coupled electron-transfer reactions are central to enzymatic mechanism in many proteins. In several enzymes, essential electron-transfer reactions involve oxidation and reduction of tyrosine side chains. For these redox-active tyrosines, proton transfer couples with electron transfer, because the phenolic pKA of the tyrosine is altered by changes in the tyrosine redox state. To develop an experimentally tractable peptide system in which the effect of proton and electron coupling can be investigated, we have designed a novel amino acid sequence that contains one tyrosine residue. The tyrosine can be oxidized by ultraviolet photolysis or electrochemical methods and has a potential cross-strand interaction with a histidine residue. NMR spectroscopy shows that the peptide forms a beta-hairpin with several interstrand dipolar contacts between the histidine and tyrosine side chains. The effect of the cross-strand interaction was probed by electron paramagnetic resonance and electrochemistry. The data are consistent with an increase in histidine pKA when the tyrosine is oxidized; the effect of this thermodynamic coupling is to increase tyrosyl radical yield at low pH. The coupling mechanism is attributed to an interstrand pi-cation interaction, which stabilizes the tyrosyl radical. A similar interaction between histidine and tyrosine in enzymes provides a regulatory mechanism for enzymatic electron-transfer reactions.

Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations

A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.