Protein Dynamics in the Region of the Sixth Ligand Methionine Revealed by Studies of Imidazole Binding To Rhodobacter capsulatus Cytochrome c2 Hinge Mutants,

Structural and functional insights of GSU0105, a unique multiheme cytochrome from G. sulfurreducens

Geobacter sulfurreducens possesses over 100 cytochromes that assure an effective electron transfer to the cell exterior. The most abundant group of cytochromes in this microorganism is the PpcA family, composed of five periplasmic triheme cytochromes with high structural homology and identical heme coordination (His-His). GSU0105 is a periplasmic triheme cytochrome synthetized by G. sulfurreducens in Fe(III)-reducing conditions but is not present in cultures grown on fumarate. This cytochrome has a low sequence identity with the PpcA family cytochromes and a different heme coordination, based on the analysis of its amino acid sequence. In this work, amino acid sequence analysis, site-directed mutagenesis, and complementary biophysical techniques, including ultraviolet-visible, circular dichroism, electron paramagnetic resonance, and nuclear magnetic resonance spectroscopies, were used to characterize GSU0105. The cytochrome has a low percentage of secondary structural elements, with features of α-helices and β-sheets. Nuclear magnetic resonance shows that the protein contains three low-spin hemes (Fe(II), S = 0) in the reduced state. Electron paramagnetic resonance shows that, in the oxidized state, one of the hemes becomes high-spin (Fe(III), S = 5/2), whereas the two others remain low-spin (Fe(III), S = 1/2). The data obtained also indicate that the heme groups have distinct axial coordination. The apparent midpoint reduction potential of GSU0105 (-154 mV) is pH independent in the physiological range. However, the pH modulates the reduction potential of the heme that undergoes the low- to high-spin interconversion. The reduction potential values of cytochrome GSU0105 are more distinct compared to those of the PpcA family members, providing the protein with a larger functional working redox potential range. Overall, the results obtained, together with an amino acid sequence analysis of different multiheme cytochrome families, indicate that GSU0105 is a member of a new group of triheme cytochromes.

Comparing Properties of Common Bioinorganic Ligands with Switchable Variants of Cytochrome c

Ligand substitution at the metal center is common in catalysis and signal transduction of metalloproteins. Understanding the effects of particular ligands, as well as the polypeptide surrounding, is critical for uncovering mechanisms of these biological processes and exploiting them in the design of bioinspired catalysts and molecular devices. A series of switchable K79G/M80X/F82C (X = Met, His, or Lys) variants of cytochrome (cyt) c was employed to directly compare the stability of differently ligated proteins and activation barriers for Met, His, and Lys replacement at the ferric heme iron. Studies of these variants and their nonswitchable counterparts K79G/M80X have revealed stability trends Met < Lys < His and Lys < His < Met for the protein Fe^III-X and Fe^II-X species, respectively. The differences in the hydrogen-bonding interactions in folded proteins and in solvation of unbound X in the unfolded proteins explain these trends. Calculations of free energy of ligand dissociation in small heme model complexes reveal that the ease of the Fe^III-X bond breaking increases in the series amine < imidazole < thioether, mirroring trends in hardness of these ligands. Experimental rate constants for X dissociation in differently ligated cyt c variants are consistent with this sequence, but the differences between Met and His dissociation rates are attenuated because the former process is limited by the heme crevice opening. Analyses of activation parameters and comparisons to those for the Lys-to-Met ligand switch in the alkaline transition suggest that ligand dissociation is entropically driven in all the variants and accompanied by Lys protonation at neutral pH. The described thiolate redox-linked switches have offered a wealth of new information about interactions of different protein-derived ligands with the heme iron in cyt c model proteins, and we anticipate that the strategy of employing these switches could benefit studies of other redox metalloproteins and model complexes.

Apolar distal pocket mutants of yeast cytochrome c peroxidase: Binding of imidazole, 1-methylimidazole and 4-nitroimidazole to the triAla, triVal, and triLeu variants

Imidazole binding to three apolar distal heme pocket mutants of yeast cytochrome c peroxidase (CcP) has been investigated between pH4 and 8. The three CcP variants have Arg-48, Trp-51, and His-52 mutated to either all alanine, CcP(triAla), all valine, CcP(triVal), or all leucine residues, CcP(triLeu). The imidazole binding curves for all three mutants are biphasic indicating that each of the mutants exists in at least two conformational states with different affinities for imidazole. At pH7, the high-affinity conformations of the three CcP mutants bind imidazole between 3.8 and 4.7 orders of magnitude stronger than that of wild-type CcP while the low-affinity conformations have binding affinities about 2.5 orders of magnitude larger than wild-type CcP. Imidazole binding to the three CcP mutants is pH dependent with the strongest binding observed at high pH. Apparent pK(a) values for the transition in binding vary between 5.6 and 7.5 for the high-affinity conformations and between 6.2 and 6.8 for the low-affinity conformations of the CcP triple mutants. The kinetics of imidazole binding are also biphasic. The fast phase of imidazole binding to CcP(triAla) and CcP(triLeu) is linearly dependent on the imidazole concentration while the slow phase is independent of imidazole concentration. Both phases of imidazole binding to CcP(triVal) have a hyperbolic dependence on the imidazole concentration. The apparent association rate constants vary between 30 and 170 M(-1)s(-1) while the apparent dissociation rate constants vary between 0.05 and 0.43 s(-1). The CcP triple mutants have higher binding affinities for 1-methylimidazole and 4-nitroimidazole than does wild-type CcP.

Binding of imidazole, 1-methylimidazole and 4-nitroimidazole to yeast cytochrome c peroxidase (CcP) and the distal histidine mutant, CcP(H52L)

Imidazole, 1-methylimidazole and 4-nitroimidazole bind to yeast cytochrome c peroxidase (yCcP) with apparent equilibrium dissociation constants (KD(app)) of 3.3±0.4, 0.85±0.11, and ~0.2M, respectively, at pH7. This is the weakest imidazole binding to a heme protein reported to date and it is about 120 times weaker than imidazole binding to metmyoglobin. Spectroscopic changes associated with imidazole and 1-methylimidazole binding to yCcP suggest partial ionization of bound imidazole to imidazolate. The pKa for ionization of bound imidazole is estimated to be 7.4±0.2, about 7 units lower than that of free imidazole and about 3 units lower than imidazole bound to metmyoglobin. Equilibrium binding of imidazole to CcP(H52L) is biphasic with low- and high-affinity phases having KD(app) values of 9.5±4.5 and 0.13±0.04M, respectively. CcP(H52L) binding of 1-methylimidazole is monophasic with an affinity similar to those of yCcP and rCcP. Binding of 1-methylimidazole to rCcP is associated with two kinetic phases, the initial binding complete within 10s, followed by a process that is consistent with 1-methylimidazole binding to a cavity created by movement of Trp-191 from the interior of the protein to the surface. Both the equilibrium binding and kinetics of 1-methylimidazole binding to yCcP are pH dependent. yCcP has a four-fold increase in 1-methylimidazole binding affinity on decreasing the pH from 7.5 to 4.0, an observation that is unique among the many studies on binding of imidazole and imidazole derivatives to heme proteins.

Kinetic and equilibrium studies of acrylonitrile binding to cytochrome c peroxidase and oxidation of acrylonitrile by cytochrome c peroxidase compound I

Ferric heme proteins bind weakly basic ligands and the binding affinity is often pH dependent due to protonation of the ligand as well as the protein. In an effort to find a small, neutral ligand without significant acid/base properties to probe ligand binding reactions in ferric heme proteins we were led to consider the organonitriles. Although organonitriles are known to bind to transition metals, we have been unable to find any prior studies of nitrile binding to heme proteins. In this communication we report on the equilibrium and kinetic properties of acrylonitrile binding to cytochrome c peroxidase (CcP) as well as the oxidation of acrylonitrile by CcP compound I. Acrylonitrile binding to CcP is independent of pH between pH 4 and 8. The association and dissociation rate constants are 0.32±0.16 M(-1) s(-1) and 0.34±0.15 s(-1), respectively, and the independently measured equilibrium dissociation constant for the complex is 1.1±0.2 M. We have demonstrated for the first time that acrylonitrile can bind to a ferric heme protein. The binding mechanism appears to be a simple, one-step association of the ligand with the heme iron. We have also demonstrated that CcP can catalyze the oxidation of acrylonitrile, most likely to 2-cyanoethylene oxide in a "peroxygenase"-type reaction, with rates that are similar to rat liver microsomal cytochrome P450-catalyzed oxidation of acrylonitrile in the monooxygenase reaction. CcP compound I oxidizes acrylonitrile with a maximum turnover number of 0.61 min(-1) at pH 6.0.

Tyrosine triad at the interface between the Rieske iron-sulfur protein, cytochrome c1 and cytochrome c2 in the bc1 complex of Rhodobacter capsulatus

A triad of tyrosine residues (Y152-154) in the cytochrome c(1) subunit (C1) of the Rhodobacter capsulatus cytochrome bc(1) complex (BC1) is ideally positioned to interact with cytochrome c(2) (C2). Mutational analysis of these three tyrosines showed that, of the three, Y154 is the most important, since its mutation to alanine resulted in significantly reduced levels, destabilization, and inactivation of BC1. A second-site revertant of this mutant that regained photosynthetic capacity was found to have acquired two further mutations-A181T and A200V. The Y152Q mutation did not change the spectral or electrochemical properties of C1, and showed wild-type enzymatic C2 reduction rates, indicating that this mutation did not introduce major structural changes in C1 nor affect overall activity. Mutations Y153Q and Y153A, on the other hand, clearly affect the redox properties of C1 (e.g. by lowering the midpoint potential as much as 117 mV in Y153Q) and the activity by 90% and 50%, respectively. A more conservative Y153F mutant on the other hand, behaves similarly to wild-type. This underscores the importance of an aromatic residue at position Y153, presumably to maintain close packing with P184, which modeling indicates is likely to stabilize the sixth heme ligand conformation.

Methionine ligand lability in bacterial monoheme cytochromes c: an electrochemical study

The direct electrochemical analysis of adsorbed redox active proteins has proven to be a powerful technique in biophysical chemistry, frequently making use of the electrode material pyrolytic "edge-plane" graphite. However, many heme-bearing proteins such as cytochromes c have been also examined systematically at alkanethiol-modified gold surfaces, and previously we reported the characterization of the redox properties of a series of bacterial cytochromes c in a side-by-side comparison of carbon and gold electrode materials. In our prior findings, we reported an unanticipated, low potential (E(m) ∼ -100 mV vs SHE) redox couple that could be analogously observed when a variety of monoheme cytochromes c are adsorbed onto carbon-based electrodes. Here we demonstrate that our prior phenomological data can be understood quantitatively in the loss of the methionine ligand of the heme iron, using the cytochrome c from Hydrogenbacter thermophilum as a model system. Through the comparison of wild-type protein with M61H and M61A mutants, in direct electrochemical analyses conducted as a function of temperature and exogenous ligand concentration, we are able to show that Met-ligated cytochromes c have a propensity to lose their Met ligand at graphite surfaces, and that energetics of this process (6.3 ± 0.2 kJ/mol) is similar to the energies associated with "foldons" of known protein folding pathways.

Binding of imidazole to the heme of cytochrome c1 and inhibition of the bc1 complex from Rhodobacter sphaeroides: II. Kinetics and mechanism of binding

The kinetics of imidazole (Im) and N-methylimidazole (MeIm) binding to oxidized cytochrome (cyt) c(1) of detergent-solubilized bc(1) complex from Rhodobacter sphaeroides are described. The rate of formation of the cyt c(1)-Im complex exhibited three separated regions of dependence on the concentration of imidazole: (i) below 8 mM Im, the rate increased with concentration in a parabolic manner; (ii) above 20 mM, the rate leveled off, indicating a rate-limiting conformational step with lifetime approximately 1 s; and (iii) at Im concentrations above 100 mM, the rate substantially increased again, also parabolically. In contrast, binding of MeIm followed a simple hyperbolic concentration dependence. The temperature dependences of the binding and release kinetics of Im and MeIm were also measured and revealed very large activation parameters for all reactions. The complex concentration dependence of the Im binding rate is not consistent with the popular model for soluble c-type cytochromes in which exogenous ligand binding is preceded by spontaneous opening of the heme cleft, which becomes rate-limiting at high ligand concentrations. Instead, binding of ligand to the heme is explained by a model in which an initial and superficial binding facilitates access to the heme by disruption of hydrogen-bonded structures in the heme domain. For imidazole, two separate pathways of heme access are indicated by the distinct kinetics at low and high concentration. The structural basis for ligand entry to the heme cleft is discussed.

Binding of imidazole to the heme of cytochrome c1 and inhibition of the bc1 complex from Rhodobacter sphaeroides: I. Equilibrium and modeling studies

We have used imidazole (Im) and N-methylimidazole (MeIm) as probes of the heme-binding cavity of membrane-bound cytochrome (cyt) c(1) in detergent-solubilized bc(1) complex from Rhodobacter sphaeroides. Imidazole binding to cyt c(1) substantially lowers the midpoint potential of the heme and fully inhibits bc(1) complex activity. Temperature dependences showed that binding of Im (K(d) approximately 330 microM, 25 degrees C, pH 8) is enthalpically driven (DeltaH(0) = -56 kJ/mol, DeltaS(0) = -121 J/mol/K), whereas binding of MeIm is 30 times weaker (K(d) approximately 9.3 mM) and is entropically driven (DeltaH(0) = 47 kJ/mol, DeltaS(0)(o) = 197 J/mol/K). The large enthalpic and entropic contributions suggest significant structural and solvation changes in cyt c(1) triggered by ligand binding. Comparison of these results with those obtained previously for soluble cyts c and c(2) suggested that Im binding to cyt c(1) is assisted by formation of hydrogen bonds within the heme cleft. This was strongly supported by molecular dynamics simulations of Im adducts of cyts c, c(2), and c(1), which showed hydrogen bonds formed between the N(delta)H of Im and the cyt c(1) protein, or with a water molecule sequestered with the ligand in the heme cleft.

Plasmon waveguide resonance spectroscopic evidence for differential binding of oxidized and reduced Rhodobacter capsulatus cytochrome c2 to the cytochrome bc1 complex mediated by the conformation of the Rieske iron-sulfur protein

The dissociation constants for the binding of Rhodobacter capsulatus cytochrome c2 and its K93P mutant to the cytochrome bc1 complex embedded in a phospholipid bilayer were measured by plasmon waveguide resonance spectroscopy in the presence and absence of the inhibitor stigmatellin. The reduced form of cytochrome c2 strongly binds to reduced cytochrome bc1 (Kd = 0.02 microM) but binds much more weakly to the oxidized form (Kd = 3.1 microM). In contrast, oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a biphasic fashion with Kd values of 0.11 and 0.58 microM. Such a biphasic interaction is consistent with binding to two separate sites or conformations of oxidized cytochrome c2 and/or cytochrome bc1. However, in the presence of stigmatellin, we find that oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a monophasic fashion with high affinity (Kd = 0.06 microM) and reduced cytochrome c2 binds less strongly (Kd = 0.11 microM) but approximately 30-fold more tightly than in the absence of stigmatellin. Structural studies with cytochrome bc1, with and without the inhibitor stigmatellin, have led to the proposal that the Rieske protein is mobile, moving between the cytochrome b and cytochrome c1 components during turnover. In one conformation, the Rieske protein binds near the heme of cytochrome c1, while the cytochrome c2 binding site is also near the cytochrome c1 heme but on the opposite side from the Rieske site, where cytochrome c2 cannot directly interact with Rieske. However, the inhibitor, stigmatellin, freezes the Rieske protein iron-sulfur cluster in a conformation proximal to cytochrome b and distal to cytochrome c1. We conclude from this that the dual conformation of the Rieske protein is primarily responsible for biphasic binding of oxidized cytochrome c2 to cytochrome c1. This optimizes turnover by maximizing binding of the substrate, oxidized cytochrome c2, when the iron-sulfur cluster is proximal to cytochrome b and minimizing binding of the product, reduced cytochrome c2, when it is proximal to cytochrome c1.

Local stability of Rhodobacter capsulatus cytochrome c2 probed by solution phase hydrogen/deuterium exchange and mass spectrometry

The hydrogen/deuterium exchange kinetics of Rhodobacter capsulatus cytochrome c2 have been determined using mass spectrometry. As expected, the relative domain stability was generally similar to that of the cytochrome c2 structural homolog, horse heart cytochrome c, but we were able to find evidence to support the presence of a second, small beta-sheet not found in the horse cytochrome, which stabilizes a structural region dominated by Omega loops. Importantly, we find that the so-called hinge region, comprised of 15 amino acids, which include the methionine sixth heme ligand (M96), is destabilized on oxidation, and this destabilization is propagated to a portion of the second Omega loop, most likely through perturbation of two hydrogen bonds that couple these two domains in the three dimensional structure. The mutation of a lysine at position 93 to proline amplifies the destabilization observed on oxidation of the wild-type cytochrome c2 and results in further destabilization observed in regions 52-60, 75-82, and 83-97. This suggests that hydrogen bond interactions involving two bound waters, the T94 hydroxyl, the front heme propionate and the Y75 hydroxyl, are significantly compromised upon mutation. In summary, these observations are consistent with the approximately 20-fold increase in the movement of the hinge away from the heme face in the oxidized cytochrome c2 as determined by ligand binding kinetics. Thus, H/D exchange kinetics can be used to identify relatively subtle structural features and at least in some cases facilitate the understanding of the structural basis of the dynamic properties of proteins.

The effect of replacing the axial methionine ligand with a lysine residue in cytochrome c-550 from Paracoccus versutus assessed by X-ray crystallography and unfolding

The structure of cytochrome c-550 from the nonphotosynthetic bacteria Paraccocus versutus has been solved by X-ray crystallography to 1.90 A resolution, and reveals a high structural homology to other bacterial cytochromes c(2). The effect of replacing the axial heme-iron methionine ligand with a lysine residue on protein structure and unfolding has been assessed using the M100K variant. From X-ray structures at 1.95 and 1.55 A resolution it became clear that the amino group of the lysine side chain coordinates to the heme-iron. Structural differences compared to the wild-type protein are confined to the lysine ligand loop connecting helices four and five. In the heme cavity an additional water molecule is found which participates in an H-bonding interaction with the lysine ligand. Under cryo-conditions extra electron density in the lysine ligand loop is revealed, leading to residues K97 to T101 being modeled with a double main-chain conformation. Upon unfolding, dissociation of the lysine ligand from the heme-iron is shown to be pH dependent, with NMR data consistent with the occurrence of a ligand exchange mechanism similar to that seen for the wild-type protein.

The Effects of Ligand Exchange and Mobility on the Peroxidase Activity of a Bacterial Cytochrome c upon Unfolding

The effect on the heme environment upon unfolding Paracoccus versutus ferricytochrome c-550 and two site-directed variants, K99E and H118Q, has been assessed through a combination of peroxidase activity increase and one-dimensional NMR spectroscopy. At pH 4.5, the data are consistent with a low- to high-spin heme transition, with the K99E mutation resulting in a protein with increased peroxidase activity in the absence of or at low concentrations of denaturant. Furthermore, the mobility of the polypeptide chain at pH 4.5 for the wild-type protein has been monitored in the absence and presence of denaturant through heteronuclear NMR experiments. The results are discussed in terms of local stability differences between bacterial and mitochondrial cytochromes c that are inferred from peroxidase activity assays. At pH 7.0, a mixture of misligated heme states arising from protein-based ligands assigned to lysine and histidine is detected. At low denaturant concentrations, these partially unfolded misligated heme forms inhibit the peroxidase activity. Data from the K99E mutation at pH 7.0 indicate that K99 is not involved in heme misligation, whereas histidine coordination is proven by the data from the H118Q variant.