A relationship between heme binding and protein stability in cytochrome b5

One small step for cytochrome P450 in its catalytic cycle, one giant leap for enzymology

The intermediates operating in the cytochrome P450 catalytic cycle have been investigated for more than half a century, fascinating many enzymologists. Each intermediate has its unique role to carry out diverse oxidations. Natural time course of the catalytic cycle is quite fast, hence, not all of the reactive intermediates could be isolated during physiological catalysis. Different high-valent iron intermediates have been proposed as primary oxidants: the candidates are compound 0 (Cpd 0, [FeOOH][Formula: see text]P450) and compound I (Cpd I, Fe(IV)[Formula: see text]O por[Formula: see text]P450). Among them, the role of Cpd I in hydroxylation is fairly well understood due the discovery of the peroxide shunt. This review endeavors to put the outstanding research efforts conducted to isolate and characterize the intermediates together. In addition to spectral features of each intermediate in the catalytic cycle, the oxidizing powers of Cpd 0 and Cpd I will be discussed along with most recent scientific findings.

Replacement of the Distal Histidine Reveals a Noncanonical Heme Binding Site in a 2-on-2 Hemoglobin

Heme ligation in hemoglobin is typically assumed by the "proximal" histidine. Hydrophobic contacts, ionic interactions, and the ligation bond secure the heme between two α-helices denoted E and F. Across the hemoglobin superfamily, several proteins also use a "distal" histidine, making the native state a bis-histidine complex. The group 1 truncated hemoglobin from Synechocystis sp. PCC 6803, GlbN, is one such bis-histidine protein. Ferric GlbN, in which the distal histidine (His46 or E10) has been replaced with a leucine, though expected to bind a water molecule and yield a high-spin iron complex at neutral pH, has low-spin spectral properties. Here, we applied nuclear magnetic resonance and electronic absorption spectroscopic methods to GlbN modified with heme and amino acid replacements to identify the distal ligand in H46L GlbN. We found that His117, a residue located in the C-terminal portion of the protein and on the proximal side of the heme, is responsible for the formation of an alternative bis-histidine complex. Simultaneous coordination by His70 and His117 situates the heme in a binding site different from the canonical site. This new holoprotein form is achieved with only local conformational changes. Heme affinity in the alternative site is weaker than in the normal site, likely because of strained coordination and a reduced number of specific heme-protein interactions. The observation of an unconventional heme binding site has important implications for the interpretation of mutagenesis results and globin homology modeling.

Distal regulation of heme binding of heme oxygenase-1 mediated by conformational fluctuations

Heme oxygenase-1 (HO-1) is an enzyme that catalyzes the oxidative degradation of heme. Since free heme is toxic to cells, rapid degradation of heme is important for maintaining cellular health. There have been useful mechanistic studies of the HO reaction based on crystal structures; however, how HO-1 recognizes heme is not completely understood because the crystal structure of heme-free rat HO-1 lacks electron densities for A-helix that ligates heme. In this study, we characterized conformational dynamics of HO-1 using NMR to elucidate the mechanism by which HO-1 recognizes heme. NMR relaxation experiments showed that the heme-binding site in heme-free HO-1 fluctuates in concert with a surface-exposed loop and transiently forms a partially unfolded structure. Because the fluctuating loop is located over 17 Å distal from the heme-binding site and its conformation is nearly identical among different crystal structures including catalytic intermediate states, the function of the loop has been unexamined. In the course of elucidating its function, we found interesting mutations in this loop that altered activity but caused little change to the conformation. The Phe79Ala mutation in the loop changed the conformational dynamics of the heme-binding site. Furthermore, the heme binding kinetics of the mutant was slower than that of the wild type. Hence, we concluded that the distal loop is involved in the regulation of the conformational change for heme binding through the conformational fluctuations. Similar to other enzymes, HO-1 effectively promotes its function using the identified distal sites, which might be potential targets for protein engineering.

SCMHBP: prediction and analysis of heme binding proteins using propensity scores of dipeptides

Background

Heme binding proteins (HBPs) are metalloproteins that contain a heme ligand (an iron-porphyrin complex) as the prosthetic group. Several computational methods have been proposed to predict heme binding residues and thereby to understand the interactions between heme and its host proteins. However, few in silico methods for identifying HBPs have been proposed.

Results

This work proposes a scoring card method (SCM) based method (named SCMHBP) for predicting and analyzing HBPs from sequences. A balanced dataset of 747 HBPs (selected using a Gene Ontology term GO:0020037) and 747 non-HBPs (selected from 91,414 putative non-HBPs) with an identity of 25% was firstly established. Consequently, a set of scores that quantified the propensity of amino acids and dipeptides to be HBPs is estimated using SCM to maximize the predictive accuracy of SCMHBP. Finally, the informative physicochemical properties of 20 amino acids are identified by utilizing the estimated propensity scores to be used to categorize HBPs. The training and mean test accuracies of SCMHBP applied to three independent test datasets are 85.90% and 71.57%, respectively. SCMHBP performs well relative to comparison with such methods as support vector machine (SVM), decision tree J48, and Bayes classifiers. The putative non-HBPs with high sequence propensity scores are potential HBPs, which can be further validated by experimental confirmation. The propensity scores of individual amino acids and dipeptides are examined to elucidate the interactions between heme and its host proteins. The following characteristics of HBPs are derived from the propensity scores: 1) aromatic side chains are important to the effectiveness of specific HBP functions; 2) a hydrophobic environment is important in the interaction between heme and binding sites; and 3) the whole HBP has low flexibility whereas the heme binding residues are relatively flexible.

Conclusions

SCMHBP yields knowledge that improves our understanding of HBPs rather than merely improves the prediction accuracy in predicting HBPs.

A missense mutation in the human cytochrome b5 gene causes 46,XY disorder of sex development due to true isolated 17,20 lyase deficiency

CONTEXT

Isolated 17,20 lyase deficiency is commonly defined by apparently normal 17α-hydroxylase activity but severely reduced 17,20 lyase activity of the bifunctional enzyme cytochrome P450 (CYP) enzyme 17A1 (CYP17A1), resulting in sex steroid deficiency but normal glucocorticoid and mineralocorticoid reserve. Cytochrome b5 (CYB5A) is thought to selectively enhance 17,20 lyase activity by facilitating the allosteric interaction of CYP17A1 with its electron donor P450 oxidoreductase (POR).

OBJECTIVE

We investigated a large consanguineous family including three siblings with 46,XY disorder of sex development (DSD) presenting with isolated 17,20 lyase deficiency.

DESIGN

We investigated the clinical and biochemical phenotype, conducted genetic analyses, and functionally characterized the identified CYB5A mutation in cell-based CYP17A1 coexpression assays.

RESULTS

All three siblings presented with 46,XY DSD, sex steroid deficiency, normal mineralocorticoids and glucocorticoids, and a urine steroid metabolome suggestive of isolated 17,20 lyase deficiency. CYP17A1 and POR sequences were normal, but we detected a homozygous CYB5A missense mutation (g.28,400A→T; p.H44L). Functional in vitro analysis revealed normal CYP17A1 17α-hydroxylase activity but severely impaired 17,20 lyase activity. In silico analysis suggested the disruption of CYB5A heme binding by p.H44L.

CONCLUSION

We have identified the first human CYB5A missense mutation as the cause of isolated 17,20 lyase deficiency in three individuals with 46,XY DSD. Detailed review of previously reported cases with apparently isolated 17,20 lyase deficiency due to mutant CYP17A1 and POR reveals impaired 17α-hydroxylase activity as assessed by steroid metabolome analysis and short cosyntropin testing. This suggests that truly isolated 17,20 lyase deficiency is observed only in individuals with inactivating CYB5A mutations.

The Bradyrhizobium japonicum frcB gene encodes a diheme ferric reductase

Iron utilization by bacteria in aerobic environments involves uptake as a ferric chelate from the environment, followed by reduction to the ferrous form. Ferric iron reduction is poorly understood in most bacterial species. Here, we identified Bradyrhizobium japonicum frcB (bll3557) as a gene adjacent to, and coregulated with, the pyoR gene (blr3555) encoding the outer membrane receptor for transport of a ferric pyoverdine. FrcB is a membrane-bound, diheme protein, characteristic of eukaryotic ferric reductases. Heme was essential for FrcB stability, as were conserved histidine residues in the protein that likely coordinate the heme moieties. Expression of the frcB gene in Escherichia coli conferred ferric reductase activity on those cells. Furthermore, reduced heme in purified FrcB was oxidized by ferric iron in vitro. B. japonicum cells showed inducible ferric reductase activity in iron-limited cells that was diminished in an frcB mutant. Steady-state levels of frcB mRNA were strongly induced under iron-limiting conditions, but transcript levels were low and unresponsive to iron in an irr mutant lacking the global iron response transcriptional regulator Irr. Thus, Irr positively controls the frcB gene. FrcB belongs to a family of previously uncharacterized proteins found in many proteobacteria and some cyanobacteria. This suggests that membrane-bound, heme-containing ferric reductase proteins are not confined to eukaryotes but may be common in bacteria.

Preparation of a biologically active apo-cytochrome b5 via heterologous expression in Escherichia coli

Cytochrome b(5) (b(5)) has been shown to modulate many cytochrome P450 (CYP)-dependent reactions. In order to elucidate the mechanism of such modulations, it is necessary to evaluate not only the effect of native b(5) on CYP-catalyzed reactions, but also that of the apo-cytochrome b(5) (apo-b(5)). Therefore, the apo-b(5) protein was prepared using a heterologous expression in Escherichia coli. The gene for rabbit b(5) was constructed from synthetic oligonucleotides using polymerase chain reaction (PCR), cloned into pUC19 plasmid and amplified in DH5 alpha cells. The gene sequence was verified by DNA sequencing. The sequence coding b(5) was cleaved from pUC19 by NdeI and XhoI restriction endonucleases and subcloned to the expression vector pET22b. This vector was used to transform E. coli BL-21 (DE3) Gold cells by heat shock. Expression of b(5) was induced with isopropyl beta-D-1-thiogalactopyranoside (IPTG). The b(5) protein, produced predominantly in its apo-form, was purified from isolated membranes of E. coli cells by chromatography on a column of DEAE-Sepharose. Using such procedures, the homogenous preparation of apo-b(5) protein was obtained. Oxidized and reduced forms of the apo-b(5) reconstituted with heme exhibit the same absorbance spectra as native b(5). The prepared recombinant apo-b(5) reconstituted with heme can be reduced by NADPH:CYP reductase. The reconstituted apo-b(5) is also fully biologically active, exhibiting the comparable stimulation effect on the CYP3A4 enzymatic activity towards oxidation of 1-phenylazo-2-hydroxynaphthalene (Sudan I) as native rabbit and human b(5).

Structural and thermodynamic encoding in the sequence of rat microsomal cytochrome b(5)

The water-soluble domain of rat microsomal cytochrome b(5) is a convenient protein with which to inspect the connection between amino acid sequence and thermodynamic properties. In the absence of its single heme cofactor, cytochrome b(5) contains a partially folded stretch of 30 residues. This region is recognized as prone to disorder by programs that analyze primary structures for such intrinsic features. The cytochrome was subjected to amino acid replacements in the folded core (I12A), in the portion that refolds only when in contact with the heme group (N57P), and in both (F35H/H39A/L46Y). Despite the difficulties associated with measuring thermodynamic quantities for the heme-bound species, it was possible to rationalize the energetic consequences of both types of replacements and test a simple equation relating apoprotein and holoprotein stability. In addition, a phenomenological relationship between the change in T(m) (the temperature at the midpoint of the thermal transition) and the change in thermodynamic stability determined by chemical denaturation was observed that could be used to extend the interpretation of incomplete holoprotein stability data. Structural information was obtained by nuclear magnetic resonance spectroscopy toward an atomic-level analysis of the effects.

Heme binding properties of heterologously expressed spinach cytochromeb6: Implications for transmembraneb-type cytochrome formation

In vivo and in vitro requirements for the formation of cytochrome b(6) were examined to analyze the mechanisms of transmembrane b-type cytochrome formation. After heterologous expression of spinach cytochrome b(6), formation of the holo-cytochrome was observed within the E. coli inner membrane. The transmembrane orientation of cytochrome b(6) appeared not to be critical for heme binding and holo-cytochrome formation. Furthermore, in vitro reconstitution of cytochrome b(6) was possible under oxidizing as well as under reducing conditions. Taken together these observations strongly indicate that transmembrane b-type cytochromes can spontaneously assemble in vitro as well as in a membrane.

Structural and thermodynamic consequences of b heme binding for monomeric apoglobins and other apoproteins

The binding of a cofactor to a protein matrix often involves a reorganization of the polypeptide structure. b Hemoproteins provide multiple examples of this behavior. In this minireview, selected monomeric and single b heme proteins endowed with distinct topological properties are inspected for the extent of induced refolding upon heme binding. To complement the data reported in the literature, original results are presented on a two-on-two globin of cyanobacterial origin (Synechococcus sp. PCC 7002 GlbN) and on the heme-containing module of FixL, an oxygen-sensing protein with the mixed alpha/beta topology of PAS domains. GlbN had a stable apoprotein that was further stabilized and locally refolded by heme binding; in contrast, apoFixLH presented features of a molten globule. Sequence analyses (helicity, disorder, and polarity) and solvent accessibility calculations were performed to identify trends in the architecture of b hemoproteins. In several cases, the primary structure appeared biased toward a partially disordered binding pocket in the absence of the cofactor.

A dynamic N-capping motif in cytochrome b5: evidence for a pH-controlled conformational switch

Apocytochrome b5 is a marginally stable protein exhibiting under native conditions a slow conformational exchange in its C-terminal region. The affected elements of secondary structure include a 3(10)-helix containing at its N-terminus a histidine Ncap and a subsequent proline. Participation of the neutral histidine side-chain in backbone amide capping lowers the imidazole pKa. To explore the nature of the conformational exchange in the protein and determine whether it is related to cis-trans isomerization of the His-Pro bond, three octapeptides encompassing the helix were synthesized and studied by NMR spectroscopy. One corresponded to the wild-type sequence, the second was the D-histidine epimer, and the third contained an alanine in place of the proline. It was found that the rates of cis-trans interconversion in the proline-containing peptides were slower than the rates of the conformational exchange in the protein. In addition, the wild-type peptide hinted at a predisposition for Ncap formation when in the trans configuration. Analysis of the pH response of the peptides and protein suggested that at pH near neutral, the conformational exchange detected in the protein involved only species with a trans His-Pro bond and could be approximated with a three-state model by which the terminal helix sampled a locally unfolded state. This state, which contained an uncapped histidine with a normal pKa, partitioned into neutral and protonated populations according to pH. The intrinsic conformational bias of the wild-type peptide and the pH-driven equilibria illustrated how a 3(10)-element could serve as a nucleation site for structural rearrangement.

Influence of point mutations on the flexibility of cytochrome b5: molecular dynamics simulations of holoproteins

Two membrane-bound isoforms of cytochrome b5 have been identified in mammals, one associated with the outer mitochondrial membrane (OM b5) and the other with the endoplasmic reticulum (microsomal, or Mc b5). The soluble heme binding domains of OM and Mc b5 have highly similar three-dimensional structures but differ significantly in physical properties, with OM b5 exhibiting higher stability due to stronger heme association. In this study, we present results of 8.5-ns length molecular dynamics simulations for rat Mc b5, bovine Mc b5, and rat OM b5, as well as for two rat OM b5 mutants that were anticipated to exhibit properties intermediate between those of rat OM b5 and the two Mc proteins: the A18S/I32L/L47R triple mutant (OM3M) and the A18S/I25L/I32L/L47R/L71S quintuple mutant (OM5M). Analysis of the structure, fluctuations, and interactions showed that the five b5 variants used in this study differed in organization of their molecular surfaces and heme binding cores in a way that could be used to explain certain experimentally observed physical differences. Overall, our simulations provided qualitative microscopic explanations of many of the differences in physical properties between OM and Mc b5 and two mutants in terms of localized changes in structure and flexibility. They also reveal that opening of a surface cleft between hydrophobic cores 1 and 2 in bovine Mc b5, observed in two previously reported simulations (E. M. Storch and V. Daggett, Biochemistry, 1995, Vol. 34, pp. 9682-9693; A. Altuve, Biochemistry, 2001, Vol. 40, pp. 9469-9483), probably resulted from removal of crystal contacts and likely does not occur on the nanosecond time scale. Finally, the MD simulations of OM5M b5 verify that stability and dynamic properties of cytochrome b5 are remarkably resistant to mutations that dramatically alter the stability and structure of the apoprotein.

Proximal influences in two-on-two globins: effect of the Ala69Ser replacement on Synechocystis sp. PCC 6803 hemoglobin

The cyanobacterium Synechocystis sp. PCC 6803 (S6803) expresses a two-on-two globin in which His46 (distal side) and His70 (proximal) function as heme iron axial ligands. His46 can be displaced by O2, CO, and CN-, among others, whereas His70 is not labile under native conditions. The residue preceding the proximal histidine has been implicated in controlling globin axial ligand reactivity; the details of the mechanism, however, are not well understood, and little information exists for bis-histidyl hexacoordinate proteins. In many vertebrate hemoglobins and in the Synechocystis protein, the position is occupied by an alanine, whereas, in myoglobins, it is a serine involved in an intricate hydrogen-bond network. We examined the role of Ala69 in S6803 hemoglobin through the effects of an Ala --> Ser replacement. The substitution resulted in minor structural perturbations, but the response of the holoprotein to temperature-, urea-, and acid-induced denaturation was measurably affected. Enhanced three-state behavior was manifested in the decoupling of heme binding and secondary-structure formation. Urea-gradient gel experiments revealed that the stability of the apoprotein was unchanged by the replacement and that a slight alteration of the folding kinetics occurred in the holoproteins. Cyanide-binding experiments were performed to assess trans effects. The apparent rate constant for association decreased 2-fold upon Ala69Ser replacement. This deceleration was attributed to a change in the lifetime of a state containing a decoordinated His46. The results demonstrated that, as in vertebrate globins and leghemoglobin, proximal influences operate to determine fundamental dynamic and thermodynamic properties of the protein.

Assembly of a transmembrane b-type cytochrome is mainly driven by transmembrane helix interactions

Folding, assembly and stability of alpha-helical membrane proteins is still not very well understood. Several of these membrane proteins contain cofactors, which are essential for their function and which can be involved in protein assembly and/or stabilization. The effect of heme binding on the assembly and stability of the transmembrane b-type cytochrome b'559 was studied by fluorescence resonance energy transfer. Cytochrome b'559 consists of two monomers of a 44 amino acid long polypeptide, which contains one transmembrane domain. The synthesis of two variants of the b'559 monomer, each carrying a specific fluorescent dye, allowed monitoring helix-helix interactions in micelles by resonance energy transfer. The measurements demonstrate that the transmembrane peptides dimerize in detergent in the absence and presence of the heme cofactor. Cofactor binding only marginally enhances dimerization and, apparently, the redox state of the heme group has no effect on dimerization.

Enhancing the stability of microsomal cytochrome b5: a rational approach informed by comparative studies with the outer mitochondrial membrane isoform

The outer mitochondrial membrane isoform of mammalian cytochrome b5 (OM b5) is much less prone to lose heme than the microsomal isoform (Mc b5), with a conserved difference at position 71 (leucine versus serine) playing a major role. We replaced Ser71 in Mc b5 with Leu, with the prediction that it would retard heme loss by diminishing polypeptide expansion accompanying rupture of the histidine to iron bonds. The strategy was partially successful in that it slowed dissociation of heme from its less stable orientation in bMc b5 (B). Heme dissociation from orientation A was accelerated to a similar extent, however, apparently owing to increased binding pocket dynamic mobility related to steric strain. A second mutation (L32I) guided by results of previous comparative studies of Mc and OM b5s diminished the steric strain, but much greater relief was achieved by replacing heme with iron deuteroporphyrin IX (FeDPIX). Indeed, the stability of the Mc(S71L) b5 FeDPIX complex is similar to that of the FeDPIX complex of OM b5. The results suggest that maximizing heme binding pocket compactness in the apo state is a useful general strategy for increasing the stability of engineered or designed proteins.