Multi-heme cytochromes--new structures, new chemistry

Heme is one of the most pervasive cofactors in nature and the c-type cytochromes represent one of the largest families of heme-containing proteins. Recent progress in bacterial genomic analysis has revealed a vast range of genes encoding novel c-type cytochromes that contain multiple numbers of heme cofactors. The genome sequence of Geobacter sulfurreducens, for example, includes some one hundred genes encoding c-type cytochromes, with around seventy of these containing two, or more, heme groups and with one protein containing an astonishing twenty seven heme groups. This wealth of cytochromes is of great significance in the respiratory flexibility shown by bacteria such as Geobacter. In addition, we are now discovering that many of these multi-heme cytochromes have associated enzymatic activities and in some cases this is revealing new chemistries. The purpose of this perspective is to describe recent progress in the structural and functional analyses of these new multi-heme cytochromes. To illustrate this we have chosen to focus on three of these cytochromes which exhibit catalytic activities; nitrite reductase, hydroxylamine oxidoreductase and tetrathionate reductase. In addition we consider the multi-heme cytochromes from Geobacter and Desulfovibrio species. Finally, we consider and contrast the repeating structural modules found in these multi-heme cytochromes.     

Kinetics and Mechanism of Mineral Respiration: How Iron Hemes Synchronize Electron Transfer Rates

Anaerobic microorganisms of the Geobacter genus are effective electron sources for the synthesis of nanoparticles, for bioremediation of polluted water, and for the production of electricity in fuel cells. In multistep reactions, electrons are transferred via iron/heme cofactors of c‐type cytochromes from the inner cell membrane to extracellular metal ions, which are bound to outer membrane cytochromes. We measured electron production and electron flux rates to 5×10⁵ e s^?1 per G. sulfurreducens. Remarkably, these rates are independent of the oxidants, and follow zero order kinetics. It turned out that the microorganisms regulate electron flux rates by increasing their Fe²⁺/Fe³⁺ ratios in the multiheme cytochromes whenever the activity of the extracellular metal oxidants is diminished. By this mechanism the respiration remains constant even when oxidizing conditions are changing. This homeostasis is a vital condition for living systems, and makes G. sulfurreducens a versatile electron source.     

Heme axial methionine fluxion in Pseudomonas aeruginosa Asn64Gln cytochrome c551

Heme axial methionine ligands in ferricytochromes c552 from Hydrogenobacter thermophilus (HT) and Nitrosomonas europaea, both members of the cyt c8 family, display fluxional behavior. The ligand motion, proposed to be inversion at sulfur, results in an unusually small range of hyperfine shifts for heme substituents in these proteins. Herein, heme axial Met fluxion is induced in a structurally homologous cytochrome c551 from Pseudomonas aeruginosa (PA) by substituting heme pocket residue Asn64 with Gln. The mutant, PA-N64Q, displays a highly compressed range of heme substituent hyperfine shifts, temperature-dependent heme methyl resonance line broadening, low rhombic magnetic anisotropy, and a magnetic axes orientation consistent with Met orientational averaging. Analysis of NMR properties of PA-N64Q demonstrates that the heme pocket of the mutant resembles that of HT. This result confirms the importance of peripheral interactions and, in particular, residue 64 in determining axial Met orientation and heme electronic structure in proteins in the cyt c8 family.     

Changes in the heme ligation during folding of a Geobacter sulfurreducens sensor GSU0935

Ligand binding and substitution reactions are important for metalloprotein folding and function. The heme sensor of a methyl-accepting chemotaxis GSU0935 is a c-type cytochrome from the bacterium Geobacter sulfurreducens. The heme domain switches one of its axial ligands from H(2)O to a low-spin ligand, presumably Met, upon reduction. The study analyzes the stability and folding kinetics of the ferric domain. Guanidine hydrochloride denaturation yields the low-spin heme species arising from coordination of the ferric heme by non-native His residues. The population of the low-spin species further increases and then declines during protein refolding. Kinetics and mutational effects suggest that His54, from the N-terminal region of the domain, is the transient ligand to the heme. The capture and release of a non-native ligand within the compact partially-folded structures illustrates the flexibility of the heme environment in GSU0935, which may relate to the domain sensor function.     

Folding, conformational changes, and dynamics of cytochromes C probed by NMR spectroscopy

NMR spectroscopy has become a vital tool for studies of protein conformational changes and dynamics. Oxidized Fe(III)cytochromes c are a particularly attractive target for NMR analysis because their paramagnetism (S = (1)/(2)) leads to high (1)H chemical shift dispersion, even for unfolded or otherwise disordered states. In addition, analysis of shifts induced by the hyperfine interaction reveals details of the structure of the heme and its ligands for native and nonnative protein conformational states. The use of NMR spectroscopy to investigate the folding and dynamics of paramagnetic cytochromes c is reviewed here. Studies of nonnative conformations formed by denaturation and by anomalous in vivo maturation (heme attachment) are facilitated by the paramagnetic, low-spin nature of native and nonnative forms of cytochromes c. Investigation of the dynamics of folded cytochromes c also are aided by their paramagnetism. As an example of this analysis, the expression in Escherichia coli of cytochrome c(552) from Nitrosomonas europaea is reported here, along with analysis of its unusual heme hyperfine shifts. The results are suggestive of heme axial methionine fluxion in N. europaea ferricytochrome c(552). The application of NMR spectroscopy to investigate paramagnetic cytochrome c folding and dynamics has advanced our understanding of the structure and dynamics of both native and nonnative states of heme proteins.     

The magnetic properties of myoglobin as studied by NMR spectroscopy

Deoxymyoglobin has been investigated by NMR spectroscopy to determine the magnetic anisotropy through pseudocontact shifts and the total magnetic susceptibility through Evans measurements. The magnetic anisotropy values were found to be Deltachi(ax)=-2.03+/-0.08 x 10(-32) m(3) and Deltachi(rh)=-1.02+/-0.09 x 10(-32) m(3). The negative value of the axial susceptibility anisotropy originates from the z tensor axis lying in the heme plane, unlike all other heme systems investigated so far. This magnetic axis is almost exactly orthogonal to the axial histidine plane. The other two axes lie essentially in the histidine plane, the closest to the heme normal being tilted by about 36 degrees from it, towards pyrrole A on the side of the proximal histidine. From the comparison with cytochrome c' it clearly appears that the position of the one axis lying in the heme plane is related to the axial histidine orientation. Irrespective of the directions, the magnetic anisotropy is smaller than that of the analogous reduced cytochrome c' and of the order of that of low-spin iron(III). The magnetic anisotropy of the system permits the measurement of residual dipolar couplings, which, together with pseudocontact shifts, prove that the solution structure is very similar to that in the crystalline state. Magnetic measurements, at variance with previous data, demonstrate that there is an orbital contribution to the magnetic moment, micro(eff)=5.5 micro(B). Finally, from the magnetic anisotropy data, the hyperfine shifts of iron ligands could be separated in pseudocontact and contact components, and hints are provided to understand the spin-delocalisation mechanism in S=2 systems by keeping in mind the delocalisation patterns in low-spin S=1/2 and high-spin S= 5/2 iron(III) systems.     

Catalytic cycles for the reduction of [UO2]2+ by cytochrome c7 proteins proposed from DFT calculations

The mechanism of the reduction of the hydrated uranyl cation, [UO2](2+), by the cytochromes G. sulfurreducens and D. acetoxidans has been studied using density functional theory calculations. We propose that the initial electron transfer step from the heme is to a cation-cation complex in the case of D. acetoxidans, but for G. sulfurreducens, it is to a single uranyl cation, which then forms a U(V)-U(VI) complex with a second uranyl cation. For both enzymes, the subsequent catalytic pathways are very similar. A U(V)-U(V) complex is formed, which then undergoes disproportionation via two successive protonation steps of one uranyl group, to give a U(VI)-U(IV) complex which dissociates to individual U(VI) and U(IV) species, the former being bound at the enzyme active site. Intermediate structures along the catalytic pathway are consistent with EXAFS data.     

Anisotropic nuclear spin interactions for the morphology analysis of proteins in solution by NMR spectroscopy.

Determining the relative orientation of domains within a protein is an important problem in structural biology, which has been difficult to address by either X-ray crystallography or NMR. The structure of a multidomain protein in a crystal lattice can be altered by crystal packing forces, resulting in different domain arrangements from those in solution. On the other hand, conventional NMR primarily provides short-range structural information, including proton-proton distances derived from nuclear Overhauser effects (NOEs) and torsion angles through vicinal spin couplings. Thus, NMR cannot always determine the precise interdomain arrangements due to the sparsely observed spin interactions between domains. However, the weak alignment of proteins in solution has enabled a new NMR technique to determine the domain arrangement based on the different structural information, which defines the orientation of a structural unit in protein against the magnetic field. This technique relies on the anisotropic nuclear spin interactions that only occur for a molecule in a weakly aligned state. In this review, the basics of the new NMR approach are described with focusing on its application to domain orientation analysis. We also describe our recently established NMR approach using the same spin interactions, which expands the domain arrangement analysis to higher-molecular weight proteins over 100 kDa.     

Assignment of the ferriheme resonances of the low-spin complexes of nitrophorins 1 and 4 by (1)H and (13)C NMR spectroscopy: comparison to structural data obtained from X-ray crystallography

In this work we report the assignment of the majority of the ferriheme resonances of low-spin nitrophorins (NP) 1 and 4 and compare them to those of NP2, published previously. It is found that the structure of the ferriheme complexes of NP1 and NP4, in terms of the orientation of the ligand(s), can be determined with good accuracy by NMR techniques in the low-spin forms and that angle plots proposed previously (Shokhirev, N. V.; Walker, F. A. J. Biol. Inorg. Chem. 1998, 3, 581-594) describe the angle of the effective nodal plane of the axial ligands in solution. The effective nodal plane of low-spin NP1, NP4, and NP2 complexes is in all cases of imidazole and histamine complexes quite similar to the average of the His-59 or -57 and the exogenous ligand angles seen in the X-ray crystal structures. For the cyanide complexes of the nitrophorins, however, the effective nodal plane of the axial ligand does not coincide with the actual histidine-imidazole plane orientation. This appears to be a result of the contribution of an additional source of asymmetry, the orientation of one of the zero-ruffling lines of the heme. Probably this effect exists for the imidazole and histamine complexes as well, but because the effect of asymmetry that occurs from planar exogenous axial ligands is much larger than the effect of heme ruffling the effect of the zero-ruffling line can only be detected for the cyanide complexes, where the only ligand plane is that of the proximal histidine. The three-dimensional structures of the three NP-CN complexes, including that of NP2-CN reported herein, confirm the high degree of ruffling of these complexes. There is an equilibrium between the two heme orientations (A and B) that depends on the heme cavity shape and changes somewhat with exogenous axial ligand. The A:B ratio can be much more accurately measured by NMR spectroscopy than by X-ray crystallography.     

Manipulation of the director in bicellar mesophases by sample spinning: a new tool for NMR spectroscopy.

It is shown that bicellar nematic liquid-crystalline phases can be oriented with the director (the normal to the bicellar plane) at an arbitrary angle to the applied magnetic field by sample rotation around one axis (variable-angle sample spinning) or around two axes successively (switched-angle spinning). This promises to open novel possibilities for NMR studies of bicelles and proteins incorporated into bicelles or dissolved in a solution containing bicelles, including the correlation of several orientations in a two-dimensional NMR experiment.     

Assignment of the ferriheme resonances of the high-spin forms of nitrophorins 1 and 4 by 1H NMR spectroscopy: comparison to structural data obtained from X-ray crystallography

In this work, we report the assignment of the majority of the ferriheme resonances of high-spin nitrophorins (NPs) 1 and 4 and compare them to those of NP2, published previously. It is found that the structures of the ferriheme complexes of NP1 and NP4, in terms of the orientation of the histidine imidazole ligand, can be described with good accuracy by NMR techniques and that the angle plot proposed previously for the high-spin form of the NPs (Shokhireva, T. Kh.; Shokhirev, N. V.; Walker, F. A. Biochemistry 2003, 42, 679-693) describes the angle of the effective nodal plane of the axial histidine imidazole in solution. There is an equilibrium between the two heme orientations (A and B), which depends on the heme cavity shape, which can be altered by mutation of amino acids with side chains (phenyl vs tyrosyl) near the potential position where a heme vinyl group would be in one of the isomers. The A:B ratio can be much more accurately measured by NMR spectroscopy than by X-ray crystallography.     

1H and 13C NMR investigation of the influence of nonligated residue contacts on the heme electronic structure in cyanometmyoglobin complexes reconstituted with centro- and pseudocentrosymmetric hemins

The 1H and 13C chemical shifts for the heme methyls of low-spin, ferric sperm whale cyanometmyoglobin reconstituted with a variety of centrosymmetric and pseudocentrosymmetric hemins have been recorded and analyzed to shed light on the nature of heme-protein contacts, other than that of the axial His, that modulate the rhombic perturbation to the heme's in-plane electronic asymmetry. The very similar 1H dipolar shifts for heme pocket residues in all complexes yield essentially the same magnetic axes as in wild type, and the resultant dipolar shifts allow the direct determination of the heme methyl proton and 13C contact shifts in all complexes. It is demonstrated that, even when the magnetic axes and anisotropies are known, the intrinsic uncertainties in the orientational parameters lead to a sufficiently large uncertainty in dipolar shift that the methyl proton contact shifts are inherently significantly less reliable indicators of the unpaired electron spin distribution than the methyl 13C contact shifts. The pattern of the noninversion symmetry in 13C contact shifts in the centro- or pseudocentrosymmetric hemes is shown to correlate with the positions of aromatic rings of Phe43(CD1) and His97(FG3) parallel to, and in contact with, the heme. These results indicate that such pi-pi interactions significantly perturb the in-plane asymmetry of the heme pi spin distribution and cannot be ignored in a quantitative interpretation of the heme methyl 13C contact shifts in terms of the axial His orientation in b-type hemoproteins.     

A justification for using NMR model-free methods when investigating the solution structures of rhombic paramagnetic lanthanide complexes

The detailed analysis of the 1H NMR hyperfine shifts according to the model-free methods shows that the semi-rigid monometallic complexes [Ln(L)(NO3)3] (Ln = Eu-Yb) are isostructural in solution. The associated separation of contact and pseudo-contact contributions to the hyperfine NMR shifts in each rhombic lanthanide complex at room temperature provides paramagnetic susceptibility tensors whose principal magnetic axes match the expected symmetry requirements. Moreover, both axial (Delta chi(ax)) and rhombic (Delta chi(rh)) paramagnetic anisotropies display satisfactory linear dependence on Bleaney's factors, a correlation predicted by the approximate high-temperature expansion of the magnetic susceptibility limited to T(-2). Consequently, the simple, and chemically attracting NMR model-free methods are not limited to axial systems, and can be safely used for the investigation of the solution structures of any lanthanide complexes. Molecular-based structural criteria for the reliable estimation of paramagnetic susceptibility tensors by NMR are discussed, together with the assignment of the labels of the crystal-field and magnetic axes within Bleaney's approach.     

Determining the orientation and localization of membrane-bound peptides

Many naturally occurring bioactive peptides bind to biological membranes. Studying and elucidating the mode of interaction is often an essential step to understand their molecular and biological functions. To obtain the complete orientation and immersion depth of such compounds in the membrane or a membrane-mimetic system, a number of methods are available, which are separated in this review into four main classes: solution NMR, solid-state NMR, EPR and other methods. Solution NMR methods include the Nuclear Overhauser Effect (NOE) between peptide and membrane signals, residual dipolar couplings and the use of paramagnetic probes, either within the membrane-mimetic or in the solvent. The vast array of solid state NMR methods to study membrane-bound peptide orientation and localization includes the anisotropic chemical shift, PISA wheels, dipolar waves, the GALA, MAOS and REDOR methods and again the use of paramagnetic additives on relaxation rates. Paramagnetic additives, with their effect on spectral linewidths, have also been used in EPR spectroscopy. Additionally, the orientation of a peptide within a membrane can be obtained by the anisotropic hyperfine tensor of a rigidly attached nitroxide label. Besides these magnetic resonance techniques a series of other methods to probe the orientation of peptides in membranes has been developed, consisting of fluorescence-, infrared- and oriented circular dichroism spectroscopy, colorimetry, interface-sensitive X-ray and neutron scattering and Quartz crystal microbalance.     

Oxidation of heme to beta- and delta-biliverdin by Pseudomonas aeruginosa heme oxygenase as a consequence of an unusual seating of the heme

The origin of the unusual regioselectivity of heme oxygenation, i.e. the oxidation of heme to delta-biliverdin (70%) and beta-biliverdin (30%), that is exhibited by heme oxygenase from Pseudomonas aeruginosa (pa-HO) has been studied by (1)H NMR, (13)C NMR, and resonance Raman spectroscopies. Whereas resonance Raman indicates that the heme-iron ligation in pa-HO is homologous to that observed in previously studied alpha-hydroxylating heme oxygenases, the NMR spectroscopic studies suggest that the heme in this enzyme is seated in a manner that is distinct from that observed for all other alpha-hydroxylating heme oxygenase enzymes for which a structure is known. In pa-HO, the heme is rotated in-plane approximately 110 degrees, so the delta-meso-carbon of the major orientational isomer is located within the HO-fold in the place where the alpha-hydroxylating enzymes typically place the alpha-meso-carbon. The unusual heme seating displayed by pa-HO places the heme propionates so that these groups point in the direction of the solvent-exposed heme edge and appears to originate in large part from the absence of stabilizing interactions between the polypeptide and the heme propionates, which are typically found in alpha-hydroxylating heme oxygenase enzymes. These interactions typically involve Lys-16 and Tyr-112, in Neisseriae meningitidis HO, and Lys-16 and Tyr-134, in human and rat HO-1. The corresponding residues in pa-HO are Asn-19 and Phe-117, respectively. In agreement with this hypothesis, we found that the Asn-19 Lys/Phe-117 Tyr double mutant of pa-HO exists as a mixture of molecules exhibiting two distinct heme seatings; one seating is identical to that exhibited by wild-type pa-HO, whereas the alternative seating is very similar to that typical of alpha-hydroxylating heme oxygenase enzymes and is related to the wild-type seating by approximately 110 degrees in-plane rotation of the heme. Furthermore, each of these heme seatings in the pa-HO double mutant gives rise to a subset of two heme isomeric orientations that are related to each other by 180 degrees rotation about the alpha-gamma-meso-axis. The coexistence of these molecules in solution, in the proportions suggested by the corresponding area under the peaks in the (1)H NMR spectrum, explains the unusual regioselectivity of heme oxygenation observed with the double mutant, which we found produces alpha- (55%), delta- (35%), and beta-biliverdin (10%). Alpha-biliverdin is obtained by oxidation of the heme seated similar to that of alpha-hydroxylating enzymes, whereas beta- and delta-biliverdin are formed from the oxidation of heme seated as in wild-type pa-HO.     

Recent developments in the <Superscript>13</Superscript>C NMR spectroscopic analysis of paramagnetic hemes and heme proteins

Despite the wealth of information that has been obtained from the study of paramagnetic hemes and heme proteins by ¹H NMR spectroscopy, there are certain limitations imposed by the nature of paramagnetically affected resonances that are difficult to overcome. Although it has long been recognized that ¹³C NMR spectroscopy is likely to be a powerful complementary technique to overcome some of these limitations, the low sensitivity and low natural abundance of ¹³C nuclei has resulted in a lag in the application of ¹³C NMR spectroscopy to the study of paramagnetic hemes and heme proteins. The tremendous advances in methodology and instrumentation witnessed in the NMR field, coupled to the advent of recombinant DNA methods that have made possible the preparation and purification of significant quantities of proteins, and the biosynthesis of ¹³C-labeled heme, have contributed to an increased interest in the study of paramagnetic heme active sites by ¹³C NMR spectroscopy. As a consequence, ¹³C NMR spectroscopy is emerging as a powerful tool to study heme electronic structure and structure–function relationships in heme-containing proteins. In this report we strive to summarize some of the recent developments in the analysis of paramagnetic hemes and heme-containing proteins by ¹³C NMR spectroscopy.     

HYSCORE spectroscopy in the cytochrome b(559) of the photosystem II reaction center

A HYSCORE investigation of the heme center in the cytochrome b(559) is presented. To assign the observed signals to specific nuclei, bis-imidazol coordinated heme compounds that model the iron environment in cytochrome b(559) are also studied. In the model compounds selective isotopic substitution of nitrogen atoms has been performed. The HYSCORE spectra allow us to obtain the hyperfine and quadrupolar coupling tensors of heme and imidazol bonding nitrogen atoms. The results can be interpreted in terms of the structure and the electronic distribution of the active center. The hyperfine tensors indicate that the unpaired electron is confined in a nonbonding iron orbital with a negligible nitrogen p orbital contribution. Quadrupolar coupling tensors suggest that the orientation of the semioccupied orbital is driven by the orientation of the two parallel imidazol rings of the axial histidine side chains. The results are discussed in terms of the structure-function relationship of cytochromes.     

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.