Toward Engineering the Stability and Hemin-Binding Properties of Microsomal Cytochromesb5into Rat Outer Mitochondrial Membrane Cytochromeb5: Examining the Influence of Residues 25 and 71†

Replacement of the heme axial lysine as a test of conformational adaptability in the truncated hemoglobin THB1

Amino acid replacement is a useful strategy to assess the roles of axial heme ligands in the function of native heme proteins. THB1, the protein product of the Chlamydomonas reinhardtii THB1 gene, is a group 1 truncated hemoglobin that uses a lysine residue in the E helix (Lys53, at position E10 by reference to myoglobin) as an iron ligand at neutral pH. Phylogenetic evidence shows that many homologous proteins have a histidine, methionine or arginine at the same position. In THB1, these amino acids would each be expected to convey distinct reactive properties if replacing the native lysine as an axial ligand. To explore the ability of the group 1 truncated Hb fold to support alternative ligation schemes and distal pocket conformations, the properties of the THB1 variants K53A as a control, K53H, K53M, and K53R were investigated by electronic absorption, EPR, and NMR spectroscopies. We found that His53 is capable of heme ligation in both the Fe(III) and Fe(II) states, that Met53 can coordinate only in the Fe(II) state, and that Arg53 stabilizes a hydroxide ligand in the Fe(III) state. The data illustrate that the group 1 truncated Hb fold can tolerate diverse rearrangement of the heme environment and has a strong tendency to use two protein side chains as iron ligands despite accompanying structural perturbations. Access to various redox pairs and different responses to pH make this protein an excellent test case for energetic and dynamic studies of heme ligation.

Heme-bound SiaA from Streptococcus pyogenes: Effects of mutations and oxidation state on protein stability

The protein SiaA (HtsA) is part of a heme uptake pathway in Streptococcus pyogenes. In this report, we present the heme binding of the alanine mutants of the axial histidine (H229A) and methionine (M79A) ligands, as well as a lysine (K61A) and cysteine (C58A) located near the heme propionates (based on homology modeling) and a control mutant (C47A). pH titrations gave pKa values ranging from 9.0 to 9.5, close to the value of 9.7 for WT SiaA. Resonance Raman spectra of the mutants suggested that the ferric heme environment may be distinct from the wild-type; spectra of the ferrous states were similar. The midpoint reduction potential of the K61A mutant was determined by spectroelectrochemical titration to be 61±3mV vs. SHE, similar to the wild-type protein (68±3mV). The addition of guanidine hydrochloride showed two processes for protein denaturation, consistent with heme loss from protein forms differing by the orientation of the heme in the binding pocket (the half-life for the slower process ranged from less than half a day to two days). The ease of protein unfolding was related to the strength of interaction of the residues with the heme. We hypothesize that kinetically facile but only partial unfolding, followed by a very slow approach to the completely unfolded state, may be a fundamental attribute of heme trafficking proteins. Small motions to release/transfer the heme accompanied by resistance to extensive unfolding may preserve the three dimensional form of the protein for further uptake and release.

Mechanistic basis of electron transfer to cytochromes p450 by natural redox partners and artificial donor constructs

Cytochromes P450 (P450s) are hemoproteins catalyzing oxidative biotransformation of a vast array of natural and xenobiotic compounds. Reducing equivalents required for dioxygen cleavage and substrate hydroxylation originate from different redox partners including diflavin reductases, flavodoxins, ferredoxins and phthalate dioxygenase reductase (PDR)-type proteins. Accordingly, circumstantial analysis of structural and physicochemical features governing donor-acceptor recognition and electron transfer poses an intriguing challenge. Thus, conformational flexibility reflected by togging between closed and open states of solvent exposed patches on the redox components was shown to be instrumental to steered electron transmission. Here, the membrane-interactive tails of the P450 enzymes and donor proteins were recognized to be crucial to proper orientation toward each other of surface sites on the redox modules steering functional coupling. Also, mobile electron shuttling may come into play. While charge-pairing mechanisms are of primary importance in attraction and complexation of the redox partners, hydrophobic and van der Waals cohesion forces play a minor role in docking events. Due to catalytic plasticity of P450 enzymes, there is considerable promise in biotechnological applications. Here, deeper insight into the mechanistic basis of the redox machinery will permit optimization of redox processes via directed evolution and DNA shuffling. Thus, creation of hybrid systems by fusion of the modified heme domain of P450s with proteinaceous electron carriers helps obviate the tedious reconstitution procedure and induces novel activities. Also, P450-based amperometric biosensors may open new vistas in pharmaceutical and clinical implementation and environmental monitoring.

Accommodating a nonconservative internal mutation by water-mediated hydrogen bonding between β-sheet strands: a comparison of human and rat type B (mitochondrial) cytochrome b5

Mammalian type B (mitochondrial) b(5) cytochromes exhibit greater amino acid sequence diversity than their type A (microsomal) counterparts, as exemplified by the type B proteins from human (hCYB5B) and rat (rCYB5B). The comparison of X-ray crystal structures of hCYB5B and rCYB5B reported herein reveals a striking difference in packing involving the five-strand β-sheet, which can be attributed to fully buried residue 21 in strand β4. The greater bulk of Leu21 in hCYB5B in comparison to that of Thr21 in rCYB5B results in a substantial displacement of the first two residues in β5, and consequent loss of two of the three hydrogen bonds between β5 and β4. Hydrogen bonding between the residues is instead mediated by two well-ordered, fully buried water molecules. In a 10 ns molecular dynamics simulation, one of the buried water molecules in the hCYB5B structure exchanged readily with solvent via intermediates having three water molecules sandwiched between β4 and β5. When the buried water molecules were removed prior to a second 10 ns simulation, β4 and β5 formed persistent hydrogen bonds identical to those in rCYB5B, but the Leu21 side chain was forced to adopt a rarely observed conformation. Despite the apparently greater ease of access of water to the interior of hCYB5B than of rCYB5B suggested by these observations, the two proteins exhibit virtually identical stability, dynamic, and redox properties. The results provide new insight into the factors stabilizing the cytochrome b(5) fold.

Effective approach for calculations of absolute stability of proteins using focused dielectric constants

The ability to predict the absolute stability of proteins based on their corresponding sequence and structure is a problem of great fundamental and practical importance. In this work, we report an extensive, refinement and validation of our recent approach (Roca et al., FEBS Lett 2007;581:2065-2071) for predicting absolute values of protein stability DeltaG(fold). This approach employs the semimacroscopic protein dipole Langevin dipole method in its linear response approximation version (PDLD/S-LRA) while using the best fitted values of the dielectric constants epsilon'(p) and epsilon'(eff) for the self energy and charge-charge interactions, respectively. The method is validated on a diverse set of 45 proteins. It is found that the best fitted values of both dielectric constants are around 40. However, the self energy of internal residues and the charge-charge interactions of Lys have to be treated with care, using a somewhat lower values of epsilon'(p) and epsilon'(eff). The predictions of DeltaG(fold) reported here, have an average error of only 1.8 kcal/mole compared to the observed values, making our method very promising for estimating protein stability. It also provides valuable insight into the complex electrostatic phenomena taking place in folded proteins.

Structural characterization of the hemophore HasAp from Pseudomonas aeruginosa: NMR spectroscopy reveals protein-protein interactions between Holo-HasAp and hemoglobin

Pseudomonas aeruginosa secretes a 205 residue long hemophore (full-length HasAp) that is subsequently cleaved at the C'-terminal domain to produce mainly a 184 residue long truncated HasAp that scavenges heme [Letoffé, S., Redeker, V., and Wandersman, C. (1998) Mol. Microbiol. 28, 1223-1234]. HasAp has been characterized by X-ray crystallography and in solution by NMR spectroscopy. The X-ray crystal structure of truncated HasAp revealed a polypeptide alphabeta fold and a ferriheme coordinated axially by His32 and Tyr75, with the side chain of His83 poised to accept a hydrogen bond from the Tyr75 phenolic acid group. NMR investigations conducted with full-length HasAp showed that the carboxyl-terminal tail (21 residues) is disordered and conformationally flexible. NMR spectroscopic investigations aimed at studying a complex between apo-HasAp and human methemoglobin were stymied by the rapid heme capture by the hemophore. In an effort to circumvent this problem NMR spectroscopy was used to monitor the titration of 15N-labeled holo-HasAp with hemoglobin. These studies allowed identification of a specific area on the surface of truncated HasAp, encompassing the axial ligand His32 loop that serves as a transient site of interaction with hemoglobin. These findings are discussed in the context of a putative encounter complex between apo-HasAp and hemoglobin that leads to efficient hemoglobin-heme capture by the hemophore. Similar experiments conducted with full-length 15N-labeled HasAp and hemoglobin revealed a transient interaction site in full-length HasAp similar to that observed in the truncated hemophore. The spectral perturbations observed while investigating these interactions, however, are weaker than those observed for the interactions between hemoglobin and truncated HasAp, suggesting that the disordered tail in the full-length HasAp must be proteolyzed in the extracellular milieu to make HasAp a more efficient hemophore.

Enhancing the thermal stability of mitochondrial cytochrome b5 by introducing a structural motif characteristic of the less stable microsomal isoform

Outer mitochondrial membrane cytochrome b5 (OM b5) is the most thermostable cytochrome b5 isoform presently known. Herein, we show that OM b5 thermal stability is substantially enhanced by swapping an apparently invariant motif in its heme-independent folding core with the corresponding motif characteristic of its less stable evolutionary relative, microsomal cytochrome b5 (Mc b5). The motif swap involved replacing two residues, Arg15 with His and Glu20 with Ser, thereby introducing a Glu11-His15-Ser20 H-bonding triad on the protein surface along with a His15/Trp22 pi-stacking interaction. The ferric and ferrous forms of the OM b5 R15H/E20S double mutant have thermal denaturation midpoints (Tm values) of approximately 93 degrees C and approximately 104 degrees C, respectively. A 15 degrees C increase in apoprotein Tm plays a key role in the holoprotein thermal stability enhancement, and is achieved by one of the most common natural mechanisms for stabilization of thermophilic versus mesophilic proteins: raising the unfolding free energy along the entire stability curve.

Comparison of cytochromes b5 from insects and vertebrates

We report a 1.55 A X-ray crystal structure of the heme-binding domain of cytochrome b(5) from Musca domestica (house fly; HF b(5)), and compare it with previously published structures of the heme-binding domains of bovine microsomal cytochrome b(5) (bMc b(5)) and rat outer mitochondrial membrane cytochrome b(5) (rOM b(5)). The structural comparison was done in the context of amino acid sequences of all known homologues of the proteins under study. We show that insect b(5)s contain an extended hydrophobic patch at the base of the heme binding pocket, similar to the one previously shown to stabilize mammalian OM b(5)s relative to their Mc counterparts. The hydrophobic patch in insects includes a residue with a bulky hydrophobic side chain at position 71 (Met). Replacing Met71 in HF b(5) with Ser, the corresponding residue in all known mammalian Mc b(5)s, is found to substantially destabilize the holoprotein. The destabilization is a consequence of two related factors: (1) a large decrease in apoprotein stability and (2) extension of conformational disruption in the apoprotein beyond the empty heme binding pocket (core 1) and into the heme-independent folding core (core 2). Analogous changes have previously been shown to accompany replacement of Leu71 in rOM b(5) with Ser. That the stabilizing role of Met71 in HF b(5) is manifested primarily in the apo state is highlighted by the fact that its crystallographic Calpha B factor is modestly larger than that of Ser71 in bMc b(5), indicating that it slightly destabilizes local polypeptide conformation when heme is in its binding pocket. Finally, we show that the final unit of secondary structure in the cytochrome b(5) heme-binding domain, a 3(10) helix known as alpha6, differs substantially in length and packing interactions not only for different protein isoforms but also for given isoforms from different species.

Insight into heme protein redox potential control and functional aspects of six-coordinate ligand-sensing heme proteins from studies of synthetic heme peptides

We describe detailed studies of peptide-sandwiched mesohemes PSMA and PSMW, which comprise two histidine (His)-containing peptides covalently attached to the propionate groups of iron mesoporphyrin II. Some of the energy produced by ligation of the His side chains to Fe in the PSMs is invested in inducing helical conformations in the peptides. Replacing an alanine residue in each peptide of PSMA with tryptophan (Trp) to give PSMW generates additional energy via Trp side chain-porphyrin interactions, which enhances the peptide helicity and stability of the His-ligated state. The structural change strengthened His-FeIII ligation to a greater extent than His-FeII ligation, leading to a 56-mV negative shift in the midpoint reduction potential at pH 8 (Em,8 value). This is intriguing because converting PSMA to PSMW decreased heme solvent exposure, which would normally be expected to stabilize FeII relative to FeIII. This and other results presented herein suggest that differences in stability may be at least as important as differences in porphyrin solvent exposure in governing redox potentials of heme protein variants having identical heme ligation motifs. Support for this possibility is provided by the results of studies from our laboratories comparing the microsomal and mitochondrial isoforms of mammalian cytochrome b5. Our studies of the PSMs also revealed that reduction of FeIII to FeII reversed the relative affinities of the first and second His ligands for Fe (K2III > K1III; K2II < K1II). We propose that this is a consequence of conformational mobility of the peptide components, coupled with the much greater ease with which FeII can be pulled from the mean plane of a porphyrin. An interesting consequence of this phenomenon, which we refer to as "dynamic strain", is that an exogenous ligand can compete with one of the His ligands in an FeII-PSM, a reaction accompanied by peptide helix unwinding. In this regard, the PSMs are better models of neuroglobin, CooA, and other six-coordinate ligand-sensing heme proteins than of stably bis(His)-ligated electron-transfer heme proteins such as cytochrome b5. Exclusive binding of exogenous ligands by the FeII form of PSMA led to positive shifts in its Em,8 value, which increases with increasing ligand strength. The possible relevance of this observation to the function of six-coordinate ligand-sensing heme proteins is discussed.

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.

Intracellular sorting of the tail-anchored protein cytochrome b5 in plants: a comparative study using different isoforms from rabbit and Arabidopsis

Tail-anchored (TA) proteins are bound to membranes by a hydrophobic sequence located very close to the C-terminus, followed by a short luminal polar region. Their active domains are exposed to the cytosol. TA proteins are synthesized on free cytosolic ribosomes and are found on the surface of every subcellular compartment, where they play various roles. The basic mechanisms of sorting and targeting of TA proteins to the correct membrane are poorly characterized. In mammalian cells, the net charge of the luminal region determines the sorting to the correct target membrane, a positive charge leading to mitochondria and negative or null charge to the endoplasmic reticulum (ER). Here sorting signals of TA proteins were studied in plant cells and compared with those of mammalian proteins, using in vitro translation-translocation and in vivo expression in tobacco protoplasts or leaves. It is shown that rabbit cytochrome b5 (cyt b5) with a negative charge is faithfully sorted to the plant ER, whereas a change to a positive charge leads to chloroplast targeting (instead of to mitochondria as observed in mammalian cells). The subcellular location of two cyt b5 isoforms from Arabidopsis thaliana (At1g26340 and At5g48810, both with positive net charge) was then determined. At5g48810 is targeted to the ER, and At1g26340 to the chloroplast envelope. The results show that the plant ER, unlike the mammalian ER, can accommodate cytochromes with opposite C-terminal net charge, and plant cells have a specific and as yet uncharacterized mechanism to sort TA proteins with the same positive C-terminal charge to different membranes.

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.

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.

Insertion of the cytochrome b5 heme-binding loop into an SH3 domain. Effects on structure and stability, and clues about the cytochrome's architecture

Under native conditions, apocytochrome b(5) exhibits a stable core and a disordered heme-binding region that refolds upon association with the cofactor. The termini of this flexible region are in close proximity, suggesting that loop closure may contribute to the thermodynamic properties of the apocytochrome. A chimeric protein containing 43 residues encompassing the cytochrome loop was constructed using the cyanobacterial photosystem I accessory protein E (PsaE) from Synechococcus sp. PCC 7002 as a structured scaffold. PsaE has the topology of an SH3 domain, and the insertion was engineered to replace its 14-residue CD loop. NMR and optical spectroscopies showed that the hybrid protein (named EbE1) was folded under native conditions and that it retained the characteristics of an SH3 domain. NMR spectroscopy revealed that structural and dynamic differences were confined near the site of loop insertion. Variable-temperature 1D NMR spectra of EbE1 confirmed the presence of a kinetic unfolding barrier. Thermal and chemical denaturations of PsaE and EbE1 demonstrated cooperative, two-state transitions; the stability of the PsaE scaffold was found only moderately compromised by the insertion, with a DeltaT(m) of 8.3 degrees C, a DeltaC(m) of 1.5 M urea, and a DeltaDeltaG degrees of 4.2 kJ/mole. The data implied that the penalty for constraining the ends of the inserted region was lower than the approximately 6.4 kJ/mole calculated for a self-avoiding chain. Extrapolation of these results to cytochrome b(5) suggested that the intrinsic stability of the folded portion of the apoprotein reflected only a small detrimental contribution from the large heme-binding domain.

Stabilizing roles of residual structure in the empty heme binding pockets and unfolded states of microsomal and mitochondrial apocytochrome b5

The microsomal (Mc) and mitochondrial (OM) isoforms of mammalian cytochrome b5 are the products of different genes, which likely arose via duplication of a primordial gene and subsequent functional divergence. Despite sharing essentially identical folds, heme-polypeptide interactions are stronger in OM b5s than in Mc b5s due to the presence of two conserved patches of hydrophobic amino acid side chains in the OM heme binding pockets. This is of fundamental interest in terms of understanding heme protein structure-function relationships, because stronger heme-polypeptide interactions in OM b5s in comparison to Mc b5s may represent a key source of their more negative reduction potentials. Herein we provide evidence that interactions amongst the amino acid side chains contributing to the hydrophobic patches in rat OM (rOM) b5 persist when heme is removed, rendering the empty heme binding pocket of rOM apo-b5 more compact and less conformationally dynamic than that in bovine Mc (bMc) apo-b5. This may contribute to the stronger heme binding by OM apo-b5 by reducing the entropic penalty associated with polypeptide folding. We also show that when bMc apo-b5 unfolds it adopts a structure that is more compact and contains greater nonrandom secondary structure content than unfolded rOM apo-b5. We propose that a more robust beta-sheet in Mc apo-b5s compensates for the absence of the hydrophobic packing interactions that stabilize the heme binding pocket in OM apo-b5s.

Mammalian mitochondrial and microsomal cytochromes b(5) exhibit divergent structural and biophysical characteristics

The only outer mitochondrial membrane cytochrome b(5) examined to date, from rat (rOM b(5)), exhibits greater stability than known mammalian microsomal (Mc) isoforms, as well as a much higher kinetic barrier for hemin dissociation and a more negative reduction potential. A BlastP search of available databases using the protein sequence of rOM b(5) as template revealed entries for analogous proteins from human (hOM b(5)) and mouse (mOM b(5)). We prepared a synthetic gene coding for the heme-binding domain of hOM b(5), and expressed the protein to high levels. The hOM protein exhibits stability, hemin-binding, and redox properties similar to those of rOM b(5), suggesting that they are characteristic of the OM b(5) subfamily. The divergence in properties between the OM and Mc b(5) isoforms in mammals can be attributed, at least in part, to the presence of two extended hydrophobic patches in the former. The biophysical properties characteristic of the OM proteins may be important in facilitating the two functions proposed for them so far, reduction of ascorbate radical and stimulation of androgen synthesis.

House fly cytochrome b5 exhibits kinetically trapped hemin and selectivity in hemin binding

We report that cytochrome b(5) (cyt b(5)) from Musca domestica (house fly) is more thermally stable than all other microsomal (Mc) cytochromes b(5) that have been examined to date. It also exhibits a much higher barrier to equilibration of the two isomeric forms of the protein, which differ by a 180 degrees rotation about the alpha-gamma-meso axis of hemin (ferric heme). In fact, hemin is kinetically trapped in a nearly statistical 1.2:1 ratio of rotational forms in freshly expressed protein. The equilibrium ratio (5.5:1) is established only upon incubation at temperatures above 37 degrees C. House fly Mc cyt b(5) is only the second b-hemoprotein that has been shown to exhibit kinetically trapped hemin at room temperature or above, the first being cyt b(5) from the outer membrane of rat liver mitochondria (rat OM cyt b(5)). Finally, we show that the small excess of one orientational isomer over the other in freshly expressed protein results from selective binding of hemin by the apoprotein, a phenomenon that has not heretofore been established for any apocyt b(5).