Heme axial methionine fluxion in Pseudomonas aeruginosa Asn64Gln cytochrome c551

Unusual Cytochrome c552 from Thioalkalivibrio paradoxus: Solution NMR Structure and Interaction with Thiocyanate Dehydrogenase

The search of a putative physiological electron acceptor for thiocyanate dehydrogenase (TcDH) newly discovered in the thiocyanate-oxidizing bacteria Thioalkalivibrio paradoxus revealed an unusually large, single-heme cytochrome c (CytC552), which was co-purified with TcDH from the periplasm. Recombinant CytC552, produced in Escherichia coli as a mature protein without a signal peptide, has spectral properties similar to the endogenous protein and serves as an in vitro electron acceptor in the TcDH-catalyzed reaction. The CytC552 structure determined by NMR spectroscopy reveals significant differences compared to those of the typical class I bacterial cytochromes c: a high solvent accessible surface area for the heme group and so-called “intrinsically disordered” nature of the histidine-rich N- and C-terminal regions. Comparison of the signal splitting in the heteronuclear NMR spectra of oxidized, reduced, and TcDH-bound CytC552 reveals the heme axial methionine fluxionality. The TcDH binding site on the CytC552 surface was mapped using NMR chemical shift perturbations. Putative TcDH-CytC552 complexes were reconstructed by the information-driven docking approach and used for the analysis of effective electron transfer pathways. The best pathway includes the electron hopping through His528 and Tyr164 of TcDH, and His83 of CytC552 to the heme group in accordance with pH-dependence of TcDH activity with CytC552.

Methionine ligand lability of homologous monoheme cytochromes c

Direct electrochemical analysis of adsorbed bacterial monoheme cytochromes c has revealed a phenomenological loss of the axial methionine when examined using pyrolytic "edge-plane" graphite (EPG) electrodes. While prior findings have reported that the Met-loss state may be quantitatively understood using the cytochrome c from Hydrogenobacter thermophilus as a model system, here we demonstrate that the formation of the Met-loss state upon EPG electrodes can be observed for a range of cytochrome orthologs. Through an electrochemical comparison of the wild-type proteins from organisms of varying growth temperature optima, we establish that Met-ligand losses at graphite surfaces have similar energetics to the "foldons" for known protein folding pathways. Furthermore, a downward shift in reduction potential to approximately -100 mV vs standard hydrogen electrode was observed, similar to that of the alkaline transition found in mitochondrial cytochromes c. Pourbaix diagrams for the Met-loss forms of each cytochrome, considered here in comparison to mutants where the Met-ligand has been substituted to His or Ala, suggest that the nature of the Met-loss state is distinct from either a His-/aquo- or a bis-His-ligated heme center, yet more closely matches the pKa values found for bis-His-ligated hemes., We find the propensity for adoption of the Met-loss state in bacterial monoheme cytochromes c scales with their overall thermal stability, though not with the specific stability of the Fe-Met bond.

NMR investigations of nitrophorin 2 belt side chain effects on heme orientation and seating of native N-terminus NP2 and NP2(D1A)

Nitrophorin 2 (NP2), one of the four NO-storing and NO-releasing proteins found in the saliva of the blood-sucking bug Rhodnius prolixus, has a more ruffled heme and a high preference for a particular heme orientation (B) compared with nitrophorin 1 and nitrophorin 4, which show not a preference (A to B ratio of approximately 1:1), suggesting that it fits more tightly in the β-barrel protein. In this work we have prepared a series of "belt" mutants of NP2(D1A) and (ΔM0)NP2 aimed at reducing the size of aromatic or other residues that surround the heme, and investigated them as the high-spin aqua and low-spin N-methylimidazole complexes. The belt mutants included Y38A, Y38F, F42A, F66A, Y85A, Y85F, Y104A, I120T, and a triple mutant of NP2(D1A), the F42L, L106F, I120T mutant. Although I120 has been mainly considered to be a distal pocket residue, CδH3 of I120 lies directly above the heme 3-methyl, at 2.67 Å, of heme orientation B, or the 2-vinyl of A, and it thus plays a role as a belt mutant, a role that turns out to be extremely important in creating the strong favoring of the B heme orientation [A to B ratio of 1:14 for NP2(D1A) or 1:12 for (ΔM0)NP2]. The results show that the 1D (1)H NMR spectra of the high-spin forms are quite sensitive to changes in the shape of the heme binding cavity. The single mutation I120T eliminates the favorability of the B heme orientation by producing a heme A to B orientation ratio of 1:1, whereas the single mutation F42A reverses the heme orientation from an A to B ratio of 1:14 seen for NP2(D1A) to 10:1 for NP2(D1A,F42A). The most extreme ratio was found for the triple mutant of NP2(D1A), NP2(D1A,F42L,L105F,I120T), in which the A to B ratio is approximately 25:1, a ΔG change of about -3.5 kcal/mol or -14.1 kJ/mol with respect to NP2(D1A). The seating of the heme is modified as well in that mutant and in several others, by rotations of the heme by up to 4° from the seating observed in NP2(D1A), in order to relieve steric interactions between a vinyl β-carbon and a protein side chain, or to fill a cavity created by replacing a large protein side chain by a much smaller one; the latter was observed for all tyrosine to alanine mutants. These relatively small changes in seating have a measurable effect on the NMR spectra of the mutants, but are indeed minor in terms of overall seating and reactivity of the NP2(D1A) protein. The (1)H NMR resonances of the hemin substituents of the low-spin N-methylimidazole complexes of NP2(D1A,F42L,L105F,I120T) as well as NP2(D1A,I120T), NP2(D1A,Y104A), and NP2(D1A,F42A) have been assigned using natural abundance (1)H{(13)C} heteronuclear multiple quantum correlation and (1)H-(1)H nuclear Overhauser effect spectroscopy spectra.

Redox state dependence of axial ligand dynamics in Nitrosomonas europaea cytochrome c552

Analysis of NMR spectra reveals that the heme axial Met ligand orientation and dynamics in Nitrosomonas europaea cytochrome c552 (Ne cyt c) are dependent on the heme redox state. In the oxidized state, the heme axial Met is fluxional, interconverting between two conformers related to each other by inversion through the Met δS atom. In the reduced state, there is no evidence of fluxionality, with the Met occupying one conformation similar to that seen in the homologous Pseudomonas aeruginosa cytochrome c551. Comparison of the observed and calculated pseudocontact shifts for oxidized Ne cyt c using the reduced protein structure as a reference structure reveals a redox-dependent change in the structure of the loop bearing the axial Met (loop 3). Analysis of nuclear Overhauser effects (NOEs) and existing structural data provides further support for the redox state dependence of the loop 3 structure. Implications for electron transfer function are discussed.

Modulation of ligand-field parameters by heme ruffling in cytochromes c revealed by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectra of variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht c-552) and Pseudomonas aeruginosa cytochrome c(551) (Pa c-551) are analyzed to determine the effect of heme ruffling on ligand-field parameters. Mutations introduced at positions 13 and 22 in Ht c-552 were previously demonstrated to influence hydrogen bonding in the proximal heme pocket and to tune reduction potential (E(m)) over a range of 80 mV [Michel, L. V.; Ye, T.; Bowman, S. E. J.; Levin, B. D.; Hahn, M. A.; Russell, B. S.; Elliott, S. J.; Bren, K. L. Biochemistry 2007, 46, 11753-11760]. These mutations are shown here to also increase heme ruffling as E(m) decreases. The primary effect on electronic structure of increasing heme ruffling is found to be a decrease in the axial ligand-field term Δ/λ, which is proposed to arise from an increase in the energy of the d(xy) orbital. Mutations at position 7, previously demonstrated to influence heme ruffling in Pa c-551 and Ht c-552, are utilized to test this correlation between molecular and electronic structure. In conclusion, the structure of the proximal heme pocket of cytochromes c is shown to play a role in determining heme conformation and electronic structure.

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.

Review: studies of ferric heme proteins with highly anisotropic/highly axial low spin (S = 1/2) electron paramagnetic resonance signals with bis-histidine and histidine-methionine axial iron coordination

Six-coordinated heme groups are involved in a large variety of electron transfer reactions because of their ability to exist in both the ferrous (Fe(2+)) and ferric (Fe(3+)) state without any large differences in structure. Our studies on hemes coordinated by two histidines (bis-His) and hemes coordinated by histidine and methionine (His-Met) will be reviewed. In both of these coordination environments, the heme core can exhibit ferric low spin (electron paramagnetic resonance EPR) signals with large g(max) values (also called Type I, highly anisotropic low spin, or highly axial low spin, HALS species) as well as rhombic EPR (Type II) signals. In bis-His coordinated hemes rhombic and HALS envelopes are related to the orientation of the His groups with respect to each other such that (i) parallel His planes results in a rhombic signal and (ii) perpendicular His planes results in a HALS signal. Correlation between the structure of the heme and its ligands for heme with His-Met axial ligation and ligand-field parameters, as derived from a large series of cytochrome c variants, show, however, that for such a combination of axial ligands there is no clear-cut difference between the large g(max) and the "small g-anisotropy" cases as a result of the relative Met-His arrangements. Nonetheless, a new linear correlation links the average shift delta of the heme methyl groups with the g(max) values.

1H and 13C NMR spectroscopic studies of the ferriheme resonances of three low-spin complexes of wild-type nitrophorin 2 and nitrophorin 2(V24E) as a function of pH

The ferriheme resonances of the low-spin (S = 1/2) complexes of wild-type (wt) nitrophorin 2 (NP2) and its heme pocket mutant NP2(V24E) with imidazole (ImH), histamine (Hm), and cyanide (CN⁻) as the sixth ligand have been investigated by NMR spectroscopy as a function of pH (4.0–7.5). For the three wt NP2 complexes, the ratio of the two possible heme orientational isomers, A and B, remains almost unchanged (ratio of A:B approximately 1:6 to 1:5) over this wide pH range. However, strong chemical exchange cross peaks appear in the nuclear Overhauser effect spectroscopy/exchange spectroscopy (NOESY/EXSY) spectra for the heme methyl resonances at low pH (pH* 4.0–5.5), which indicate chemical exchange between two species. We have shown these to be two different exogenous ImH or Hm orientations that are denoted B and B′, with the ImH plane nearly parallel and perpendicular to the ImH plane of the protein-provided His57, respectively. The wt NP2–CN complex also shows EXSY cross peaks due to chemical exchange, which is shown to be a result of interchange between two ruffling distortions of the heme. The same ruffling distortion interchange is also responsible for the ImH and Hm chemical exchange. For the three NP2(V24E) ligand complexes, no EXSY cross peaks are observed, but the A:B ratios change dramatically with pH. The fact that heme favors the A orientation highly for NP2(V24E) at low pH as compared with wt NP2 is believed to be due to the steric effect of the V24E mutation. The existence of the B′ species at lower pH for wt NP2 complexes and the increase in A heme orientation at lower pH for NP2(V24E) are believed to be a result of a change in structure near Glu53 when it is protonated at low pH. ¹H{¹³C} heteronuclear multiple quantum coherence (HMQC) spectra are very helpful for the assignment of heme and nearby protein side chain resonances.

Modulation of the ligand-field anisotropy in a series of ferric low-spin cytochrome c mutants derived from Pseudomonas aeruginosa cytochrome c-551 and Nitrosomonas europaea cytochrome c-552: a nuclear magnetic resonance and electron paramagnetic resonance study

Cytochromes of the c type with histidine-methionine (His-Met) heme axial ligation play important roles in electron-transfer reactions and in enzymes. In this work, two series of cytochrome c mutants derived from Pseudomonas aeruginosa (Pa c-551) and from the ammonia-oxidizing bacterium Nitrosomonas europaea (Ne c-552) were engineered and overexpressed. In these proteins, point mutations were induced in a key residue (Asn64) near the Met axial ligand; these mutations have a considerable impact both on heme ligand-field strength and on the Met orientation and dynamics (fluxionality), as judged by low-temperature electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectra. Ne c-552 has a ferric low-spin (S = 1/2) EPR signal characterized by large g anisotropy with g(max) resonance at 3.34; a similar large g(max) value EPR signal is found in the mitochondrial complex III cytochrome c1. In Ne c-552, deletion of Asn64 (NeN64Delta) changes the heme ligand field from more axial to rhombic (small g anisotropy and g(max) at 3.13) and furthermore hinders the Met fluxionality present in the wild-type protein. In Pa c-551 (g(max) at 3.20), replacement of Asn64 with valine (PaN64V) induces a decrease in the axial strain (g(max) at 3.05) and changes the Met configuration. Another set of mutants prepared by insertion (ins) and/or deletion (Delta) of a valine residue adjacent to Asn64, resulting in modifications in the length of the axial Met-donating loop (NeV65Delta, NeG50N/V65Delta, PaN50G/V65ins), did not result in appreciable alterations of the originally weak (Ne c-552) or very weak (Pa c-551) axial field but had an impact on Met orientation, fluxionality, and relaxation dynamics. Comparison of the electronic fingerprints in the overexpressed proteins and their mutants reveals a linear relationship between axial strain and average paramagnetic heme methyl shifts, irrespective of Met orientation or dynamics. Thus, for these His-Met axially coordinated Fe(III), the large g(max) value EPR signal does not represent a special case as is observed for bis-His axially coordinated Fe(III) with the two His planes perpendicular to each other.

Submolecular unfolding units of Pseudomonas aeruginosa cytochrome c-551

Hydrogen exchange rates for backbone amide protons of oxidized Pseudomonas aeruginosa cytochrome c-551 (P. aeruginosa cytochrome c) have been measured in the presence of low concentrations of the denaturant guanidine hydrochloride. Analysis of the data has allowed identification of submolecular unfolding units known as foldons. The highest-energy foldon bears similarity to the proposed folding intermediate for P. aeruginosa cytochrome c. Parallels are seen to the foldons of the structurally homologous horse cytochrome c, although the heme axial methionine-bearing loop has greater local stability in P. aeruginosa cytochrome c, in accord with previous folding studies. Regions of low local stability are observed to correspond with regions that interact with redox partners, providing a link between foldon properties and function.

Heme attachment motif mobility tunes cytochrome c redox potential

Hydrogen exchange (HX) rates and midpoint potentials (Em) of variants of cytochrome c from Pseudomonas aeruginosa (Pa cyt c551) and Hydrogenobacter thermophilus (Ht cyt c552) have been characterized in an effort to develop an understanding of the impact of properties of the Cys-X-X-Cys-His pentapeptide c-heme attachment (CXXCH) motif on heme redox potential. Despite structural conservation of the CXXCH motif, Ht cyt c552 exhibits a low level of protection from HX for amide protons within this motif relative to Pa cyt c551. Site-directed mutants have been prepared to determine the structural basis for and functional implications of these variations on HX behavior. The double mutant Ht-M13V/K22M displays suppressed HX within the CXXCH motif as well as a decreased Em (by 81 mV), whereas the corresponding double mutant of Pa cyt c551 (V13M/M22K) exhibits enhanced HX within the CXXCH pentapeptide and a modest increase in Em (by 30 mV). The changes in Em correlate with changes in axial His chemical shifts in the ferric proteins reflecting the extent of histidinate character. Thus, the mobility of the CXXCH pentapeptide is found to impact the His-Fe(III) interaction and therefore the heme redox potential.

Heme crevice disorder after sixth ligand displacement in the cytochrome c-551 family

1H NMR and visible absorption spectroscopy were used to monitor sixth ligand methionine displacement reactions in four members of the ferricytochrome c-551 family from Pseudomonas aeruginosa, Pseudomonas stutzeri, Pseudomonas stutzeri substrain ZoBell, and Nitrosomonas europae. Potassium cyanide displaces the methionine ligand with very modest changes in the visible spectra, but profound changes in the NMR spectra. The initial product formed kinetically, designated complex I, changes with time and/or heating to a more thermodynamically favored product termed complex II. Spectra indicate that both I and II are actually a family of closely related conformational isomers. Low temperature NMR spectra of complex II indicate that some of the isomers are in chemical exchange on the NMR time scale. High pH also displaces the methionine ligand in a manner similar to the well-known alkaline transition of mitochondrial cytochrome c. However, the reaction occurs at higher pH values and over a narrower pH range for the c-551 family, and the transition pH range is different for the different proteins studied. The final alkaline forms also show peak widths and a number of peaks indicative of multiple conformational isomers.

The heme environment of mouse neuroglobin: histidine imidazole plane orientations obtained from solution NMR and EPR spectroscopy as compared with X-ray crystallography

The 1H NMR chemical shifts of the heme methyl groups of the ferriheme complex of metneuroglobin (Du et al. in J. Am. Chem. Soc. 125:8080-8081, 2003) predict orientations of the axial histidine ligands (Shokhirev and Walker in J. Biol. Inorg. Chem. 3:581-594, 1998) that are not consistent with the X-ray data (Vallone et al. in Proteins Struct. Funct. Bioinf. 56:85-94, 2004), and the EPR spectrum (Vinck et al. in J. Am. Chem. Soc. 126:4516-4517, 2004) is only marginally consistent with these data. The reasons for these inconsistencies appear to be rooted in the high degree of aqueous solution exposure of the heme group and the fact that there are no strong hydrogen-bond acceptors for the histidine imidazole N-H protons provided by the protein. Similar inconsistencies may exist for other water-soluble heme proteins, and 1H NMR spectroscopy provides a simple means to verify whether the solution structure of the heme center is the same as or different from that in the crystalline state.