Evaluating the roles of the heme a side chains in cytochrome c oxidase using designed heme proteins

Heme a is a redox cofactor unique to cytochrome c oxidases and vital to aerobic respiration. Heme a differs from the more common heme b by two chemical modifications, the C-8 formyl group and the C-2 hydroxyethylfarnesyl group. The effects of these porphyrin substituents on ferric and ferrous heme binding and electrochemistry were evaluated in a designed heme protein maquette. The maquette scaffold chosen, [Delta7-H3m](2), is a four-alpha-helix bundle that contains two bis(3-methyl-l-histidine) heme binding sites with known absolute ferric and ferrous heme b affinities. Hemes b, o, o+16, and heme a, those involved in the biosynthesis of heme a, were incorporated into the bis(3-methyl-l-histidine) heme binding sites in [Delta7-H3m](2). Spectroscopic analyses indicate that 2 equiv of each heme binds to [Delta7-H3m](2), as designed. Equilibrium binding studies of the hemes with the maquette demonstrate the tight affinity for hemes containing the C-2 hydroxyethylfarnesyl group in both the ferric and ferrous forms. Coupled with the measured equilibrium midpoint potentials, the data indicate that the hydroxyethylfarnesyl group stabilizes the binding of both ferrous and ferric heme by at least 6.3 kcal/mol via hydrophobic interactions. The data also demonstrate that the incorporation of the C-8 formyl substituent in heme a results in a 179 mV, or 4.1 kcal/mol, positive shift in the heme reduction potential relative to heme o due to the destabilization of ferric heme binding relative to ferrous heme binding. The two substituents appear to counterbalance each other to provide for tighter heme a affinity relative to heme b in both the ferrous and ferric forms by at least 6.3 and 2.1 kcal/mol, respectively. These results also provide a rationale for the reaction sequence observed in the biosynthesis of heme a.     

Differences in thermal stability between reduced and oxidized cytochrome b562 from Escherichia coli

The thermal stabilities of ferri- and ferrocytochrome b562 were examined. Thermally induced spectral changes, monitored by absorption and second-derivative spectroscopies, followed the dissociation of the heme moiety and the increased solvation of tyrosine residue(s) located in close proximity to the heme binding site. All observed thermal transitions were independent of the rate of temperature increase (0.5-2 degrees C/min), and the denatured protein exhibited partial to near-complete reversibility upon return to ambient temperature. The extent of renaturation of cytochrome b562 is dependent on the amount of time the unfolded conformer is exposed to temperatures above the transition temperature, Tm. All thermally induced spectra changes fit a simple two-state model, and the thermal transition was assumed to be reversible. The thermal transition for ferrocytochrome b562 yielded Tm and van't Hoff enthalpy (delta HvH) values of 81.0 degrees C and 137 kcal/mol, respectively. In contrast, Tm and delta HvH values obtained for the ferricytochrome were 66.7 degrees C and 110 kcal/mol, respectively. The estimated increase in the stabilization free energy at the Tm of ferricytochrome b562 following the one-electron reduction to the ferrous form, where delta delta G = delta Tm delta Sm [delta Sm = 324 cal/(K.mol), delta Tm = 14.3 degrees C] [Becktel, W. J., & Schellman, J. A. (1987) Biopolymers 26, 1859-1877], is 4.6 kcal/mol.     

Heme redox potential control in de novo designed four-alpha-helix bundle proteins

The effects of various mechanisms of metalloporphyrin reduction potential modulation were investigated experimentally using a robust, well-characterized heme protein maquette, synthetic protein scaffold H10A24 [?CH(3)()CONH-CGGGELWKL.HEELLKK.FEELLKL.AEERLKK. L-CONH(2)()?(2)](2). Removal of the iron porphyrin macrocycle from the high dielectric aqueous environment and sequestration within the hydrophobic core of the H10A24 maquette raises the equilibrium reduction midpoint potential by 36-138 mV depending on the hydrophobicity of the metalloporphyrin structure. By incorporating various natural and synthetic metalloporphyrins into a single protein scaffold, we demonstrate a 300-mV range in reduction potential modulation due to the electron-donating/withdrawing character of the peripheral macrocycle substituents. Solution pH is used to modulate the metalloporphyrin reduction potential by 160 mV, regardless of the macrocycle architecture, by controlling the protonation state of the glutamate involved in partial charge compensation of the ferric heme. Attempts to control the reduction potential by inserting charged amino acids into the hydrophobic core at close proximity to the metalloporphyrin lead to varied success, with H10A24-L13E lowering the E(m8.5) by 40 mV, H10A24-E11Q raising it by 50 mV, and H10A24-L13R remaining surprisingly unaltered. Modifying the charge of the adjacent metalloporphyrin, +1 for iron(III) protoporphyrin IX or neutral for zinc(II) protoporphyrin IX resulted in a loss of 70 mV [Fe(III)PPIX](+) - [Fe(III)PPIX](+) interaction observed in maquettes. Using these factors in combination, we illustrate a 435-mV variation of the metalloporphyrin reduction midpoint potential in a simple heme maquette relative to the about 800-mV range observed for natural cytochromes. Comparison between the reduction potentials of the heme maquettes and other de novo designed heme proteins reveals global trends in the E(m) values of synthetic cytochromes.     

The thermodynamics of bovine and porcine insulin and proinsulin association determined by concentration difference spectroscopy

Difference spectroscopy was used to determine the equilibrium constants and thermodynamic parameters for the monomer-dimer association of bovine and porcine insulin and bovine proinsulin at pH 2.0 and 7.0. At pH 2 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine insulin were found to be -6.6 kcal/mol, -18 cal/mol-deg, and -12 kcal/mol, respectively. Porcine insulin behaved similarly to bovine insulin in its dimerization properties in that delta G degree 25, delta S degree, and delta H degree were found to be -6.8 kcal/mol, -14 cal/mol-deg, and -11 kcal/mol, respectively. At pH 7 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine insulin were found to be -7.2 kcal/mol, -16 cal/mol/deg, and -12 kcal/mol, respectively. At pH 7.0 delta G degree 25, delta S degree, and delta H degree for dimerization of porcine insulin were -6.7 kcal/mol, -11.6 cal/mol-deg, and -10 kcal/mol, respectively. The similarity in the thermodynamic parameters of both insulin species at the different pH's suggests that there are minimal structural changes at the monomer-monomer contact site over this pH range. The dimerization of both insulin species is under enthalpic control. This may suggest that the formation of the insulin dimer is not driven by hydrophobic bonding but, rather, is driven by the formation between subunits of four hydrogen bonds in an apolar environment. At pH 2 delta G degree 25, delta S degree, and delta H degree for dimerization of bovine proinsulin were found to be -5.3 kcal/mol, -26 cal/mol-deg, and -13 kcal/mol, respectively. At pH 7 delta G degree 25, delta S degree, and delta H degree for dimerization of proinsulin were -5.9 kcal/mol, -4.2 cal/mol-deg, and -7.2 kcal/mol, respectively. Although the presence of the C-peptide on proinsulin does not drastically affect the overall free energy change of dimer formation (as compared to insulin), the other thermodynamic parameters are rather drastically altered. This may be because of electrostatic interactions of groups on the C-peptide with groups on the B-chain which are near the subunit contact site in the insulin dimer.     

Hydrophobic modulation of heme properties in heme protein maquettes

We have investigated the properties of the two hemes bound to histidine in the H10 positions of the uniquely structured apo form of the heme binding four-helix bundle protein maquette [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2) for the amino acids at positions 6 (I), 13 (F) and 24 (H), respectively. The primary structure of each alpha-helix, alpha-SH, in [I(6)F(13)H(24)](2) is Ac-CGGGEI(6)WKL.H(10)EEF(13)LKK.FEELLKL.H(24)EERLKK.L-CONH(2). In our nomenclature, [I(6)F(13)H(24)] represents the disulfide-bridged di-alpha-helical homodimer of this sequence, i.e., (alpha-SS-alpha), and [I(6)F(13)H(24)](2) represents the dimeric four helix bundle composed of two di-alpha-helical subunits, i.e., (alpha-SS-alpha)(2). We replaced the histidines at positions H24 in [I(6)F(13)H(24)](2) with hydrophobic amino acids incompetent for heme ligation. These maquette variants, [I(6)F(13)I(24)](2), [I(6)F(13)A(24)](2), and [I(6)F(13)F(24)](2), are distinguished from the tetraheme binding parent peptide, [I(6)F(13)H(24)](2), by a reduction in the heme:four-helix bundle stoichiometry from 4:1 to 2:1. Iterative redesign has identified phenylalanine as the optimal amino acid replacement for H24 in the context of apo state conformational specificity. Furthermore, the novel second generation diheme [I(6)F(13)F(24)](2) maquette was related to the first generation diheme [H10A24](2) prototype, [L(6)L(13)A(24)](2) in the present nomenclature, via a sequential path in sequence space to evaluate the effects of conservative hydrophobic amino acid changes on heme properties. Each of the disulfide-linked dipeptides studied was highly helical (>77% as determined from circular dichroism spectroscopy), self-associates in solution to form a dimer (as determined by size exclusion chromatography), is thermodynamically stable (-DeltaG(H)2(O) >18 kcal/mol), and possesses conformational specificity that NMR data indicate can vary from multistructured to single structured. Each peptide binds one heme with a dissociation constant, K(d1) value, tighter than 65 nM forming a series of monoheme maquettes. Addition of a second equivalent of heme results in heme binding with a K(d2) in the range of 35-800 nM forming the diheme maquette state. Single conservative amino acid changes between peptide sequences are responsible for up to 10-fold changes in K(d) values. The equilibrium reduction midpoint potential (E(m7.5)) determined in the monoheme state ranges from -156 to -210 mV vs SHE and in the diheme state ranges from -144 to -288 mV. An observed heme-heme electrostatic interaction (>70 mV) in the diheme state indicates a syn global topology of the di-alpha-helical monomers. The heme affinity and electrochemistry of the three H24 variants studied identify the tight binding sites (K(d1) and K(d2) values <200 nM) having the lower reduction midpoint potentials (E(m7.5) values of -155 and -260 mV) with the H10 bound hemes in the parent tetraheme state of [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2). The results of this study illustrate that conservative hydrophobic amino acid changes near the heme binding site can modulate the E(m) by up to +/-50 mV and the K(d) by an order of magnitude. Furthermore, the effects of multiple single amino acid changes on E(m) and K(d) do not appear to be additive.     

Protein denaturation with guanidine hydrochloride or urea provides a different estimate of stability depending on the contributions of electrostatic interactions.

The objective of this study was to address the question of whether or not urea and guanidine hydrochloride (GdnHCl) give the same estimates of the stability of a particular protein. We previously suspected that the estimates of protein stability from GdnHCl and urea denaturation data might differ depending on the electrostatic interactions stabilizing the proteins. Therefore, 4 coiled-coil analogs were designed, where the number of intrachain and interchain electrostatic attractions (A) were systematically changed to repulsions (R): 20A, 15A5R, 10A10R, and 20R. The GdnHCl denaturation data showed that the 4 coiled-coil analogs, which had electrostatic interactions ranging from 20 attractions to 20 repulsions, had very similar [GdnHCl]1/2 values (average of congruent to 3.5 M) and, as well, their delta delta Gu values were very close to 0 (0.2 kcal/mol). In contrast, urea denaturation showed that the [urea]1/2 values proportionately decreased with the stepwise change from 20 electrostatic attractions to 20 repulsions (20A, 7.4 M; 15A5R, 5.4 M; 10A10R, 3.2 M; and 20R, 1.4 M), and the delta delta Gu values correspondingly increased with the increasing differences in electrostatic interactions (20A-15A5R, 1.5 kcal/mol; 20A-10A10R, 3.7 kcal/mol; and 20A-20R, 5.8 kcal/mol). These results indicate that the ionic nature of GdnHCl masks electrostatic interactions in these model proteins, a phenomenon that was absent when the unchanged urea was used. Thus, GdnHCl and urea denaturations may give vastly different estimates of protein stability, depending on how important electrostatic interactions are to the protein.     

Membrane binding of the colicin E1 channel: activity requires an electrostatic interaction of intermediate magnitude.

In vitro channel activity of the C-terminal colicin E1 channel polypeptide under conditions of variable electrostatic interaction with synthetic lipid membranes showed distinct maxima with respect to pH and membrane surface potential. The membrane binding energy was determined from fluorescence quenching of the intrinsic tryptophans of the channel polypeptide by liposomes containing N-trinitrophenyl-phosphatidylethanolamine. Maximum in vitro colicin channel activity correlated with an intermediate magnitude of the electrostatic interaction. For conditions associated with maximum activity (40% anionic lipid, I = 0.12 M, pH 4.0), the free energy of binding was delta G approximately -9 kcal/mol, with nonelectrostatic and electrostatic components, delta Gnel approximately -5 kcal/mol and delta Gel approximately -4 kcal/mol, and an effective binding charge of +7 at pH 4.0. Binding of the channel polypeptide to negative membranes at pH 8 is minimal, whereas initial binding at pH 4 followed by a shift to pH 8 causes only 3-10% reversal of binding, implying that it is kinetically trapped, probably by a hydrophobic interaction. It was inferred that membrane binding and insertion involves an initial electrostatic interaction responsible for concentration and binding to the membrane surface. This is followed by insertion into the bilayer driven by hydrophobic forces, which are countered in the case of excessive electrostatic binding.     

Binding properties of androgen receptors. Evidence for identical receptors in rat testis, epididymis, and prostate.

Androgen receptors in crude and partially purified 105,000 X g supernatant fractions from rat testis, epididymis, and prostate were studied in vitro using a charcoal adsorption assay and sucrose gradient centrifugation. Androgen metabolism was eliminated during receptor purification allowing determination of the kinetics of [3H]-androgen-receptor complex formation. In all three tissues, receptors were found to have essentially identical capabilities to bind androgen, with the affinity for [3H] dihydrotestosterone being somewhat higher than for [3H] testosterone. Equilibrium dissociation constants for [3H] dihydrotestosterone and [3H] testosterone (KD = 2 to 5 X 10(-10) M) were estimated from independently determined rates of association (ka congruent to 6 X 10(7) M-1 h-1 for [3H] dihydrotestosterone and 2 X 10(8) M-1 h-1 for [3H] testosterone) and dissociation (t 1/2 congruent to 40 hr for [3H] dihydrotestosterone and 15 h [3H] testosterone). Evaluation of the effect of temperature on androgen receptor binding of [3H]testosterone allowed estimation of several thermodynamic parameters, including activation energies of association and dissociation (delta H congruent to 14 kcal/mol), the apparent free energy (delta G congruent to -12 kcal/mol), enthalpy (delta H congruent to -2.5 kcal/mol), and entropy (delta S congruent to 35 cal col-1 K-1). Optimum receptor binding occurred at a pH of 8. Receptor stability was greatly enhanced when bound with androgen. Receptor specificity for testosterone and dihydrotestosterone was demonstrated by competitive binding assays. The potent synthetic androgen, 7 alpha, 17 alpha-dimethyl-19-nortestosterone, inhibited binding of [3H] testosterone or [3H] dihydrotesterone nearly as well as testosterone and dihydrotestosterone while larger amounts of 5 alpha-androstane-3alpha, 17 beta-diol and nonandrogenic steroids were required. Sedimentation coefficients of androgen receptors in all unfractionated supernatants were 4 and 5 to 8 S. Differences in sedimentation coefficients were observed following (NH4)2SO4 precipitation which did not influence the binding properties of the receptors. These results, together with measurements of3alpha/beta-hydroxysteroid oxidoreductase activity in vitro, suggest that organ differences in receptor binding of [3H] dihydrotestosterone and [3H] testosterone in vivo result from relative differences in intracellular concentrations of these androgens rather than from differences in receptor affinities.     

Nonclassical hydrophobic effect in membrane binding equilibria.

The enthalpy of transfer of four different amphiphilic molecules from the aqueous phase to the lipid membrane was determined by titration calorimetry. The four molecules investigated were the potential-sensitive dye 2-(p-toluidinyl)naphthalene-6-sulfonate (TNS), the membrane conductivity inducing anion tetraphenylborate (TPB), the Ca2+ channel blocker amlodipine [B?uerle, H. D., & Seelig, J. (1991) Biochemistry 30, 7203-7211], and the positively charged local anesthetic dibucaine. All four amphiphiles penetrate into the hydrophobic part of the membrane, and their binding constants, after correcting for electrostatic effects, range between 600 M-1 for dibucaine and 60,000 M-1 for tetraphenylborate. The corresponding changes in free energy were about -6 to -9 kcal/mol. Binding of the amphiphiles to membrane vesicles composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine was accompanied by exothermic heats of reaction for all four molecules. For TNS, TPB, and amlodipine, the enthalpies of transfer were almost identical and corresponded to delta H approximately -9 kcal/mol, essentially accounting for the total free energy change. Thus, the binding of these charged amphiphiles to the hydrophobic membrane was driven by enthalpy. This is in contrast to the classical hydrophobic effect, where the transfer is considered to be entropy driven. For dibucaine, the enthalpy of transfer was smaller with delta H approximately -2 kcal/mol but was still about one-third of the total free energy change. All enthalpies of transfer exhibited a distinct temperature dependence with molar heat capacities delta Cp of -30 to -100 cal mol-1K-1 for the transfer from water to the membrane.(ABSTRACT TRUNCATED AT 250 WORDS)     

Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.     

The Gibbs-Donnan near-equilibrium system of heart

The gradients of the major inorganic ions across the plasma membrane of heart were examined to determine the factors controlling the extent and direction of the changes induced during injury, certain diseases, and electrolyte disturbances. The ionic environment was altered by changing only the concentration of inorganic phosphate, [sigma Pi]o, from 0 to 1.2 to 5 mM in the Krebs-Henseleit buffer perfusing working rat hearts. Raising [sigma Pi]o from 1.2 to 5 mM resulted in a decrease in total Mg2+ content and calculated free cytosolic [Mg2+] from 0.44 to 0.04 mM, conversion of 4 mmol of MgATP2- to ATP4- and a decrease in measured intracellular [Cl-]i from 41 to 16 mM. At all levels of [sigma Pi]o, both the [Na+]i and [K+]i were invariant at about 3 mM and 130 mM, respectively, as was the energy of hydrolysis of the terminal phosphate bond of sigma ATP, delta GATP Hydr, of -13.2 kcal/mol. The relationship maintained between the ions on both sides of the plasma membrane by the 3Na+/2K(+)transporting ATPase (EC 3.6.1.37) and an open K+ channel was: (formula; see text) The energy of the gradients of the other inorganic ions across the plasma membrane, delta G[ion]o/i, exhibited three distinct quanta of energy derived from the prime quantum of delta GATP Hydr of -13.2 kcal/mol. The second quantum was about one-third of delta GATP Hydr or +/- 4.4 kcal/mol and comprised the delta G[Na+]o/i, delta G[Mg2+]o/i, and delta G[HPO42-]o/i. These results indicated near-equilibrium was achieved by the reactants of the 3Na+/2K(+)-ATPase, the K+ channel, the Na(+)-Pi co-transporter, and a postulated net Mg2+/H2PO4- exchanger. The third quantum was one-third of delta G[Na+]o/i or about +/- 1.5 kcal/mol and comprised delta G[H+]o/i, delta G[HCO3-]o/i, and delta G[Cl-]o/i. The delta G[K+]o/i was 0, indicating near-equilibrium between the chemical energy of [K+]o/i and the E across the plasma membrane of -83 mV. It is concluded that the gradients of the major inorganic ions across the plasma membrane and the potential across that membrane constitute a Gibbs-Donnan equilibrium system catalyzed by transport enzymes sharing common substrates. The chemical and electrical energies of those gradients are equal in magnitude and opposite in sign to the chemical energy of ATP hydrolysis.     

Thermodynamics of follitropin binding to solubilized calf testis receptor

Thermodynamic parameters of follitropin binding to solubilized testicular receptors were measured in order to assess the forces involved in the binding reaction. Reversibility of follitropin binding to solubilized receptor decreased only 20% over the temperature range 4-24 degrees C, whereas earlier studies indicated reversibility of binding to membrane-bound receptor decreased by more than 40% over the same range [Anderson, T. T., Curatolo, L. M., & Reichert, L. E., Jr. (1983) Mol. Cell. Endocrinol. 33, 37-52]. Thermodynamic analysis of follitropin binding to solubilized receptors showed that the hydrophobic effect was important in the binding reaction. The mean values, at 25 degrees C, for delta H and delta S were -31.8 kcal/mol and -66.0 cal mol-1 K-1, respectively, and delta Cp was -3.0 kcal mol-1 K-1. This is an unusually large heat capacity for protein-protein association reactions, indicating an enhanced role for the hydrophobic effect with the solubilized (compared to membrane-bound) receptor. Since glycerol was necessary to stabilize the solubilized receptor, we determined whether glycerol affected the thermodynamic parameters measured for the binding reaction. Control experiments, performed with membrane-bound receptor in the presence or absence of glycerol, indicated that delta Cp actually decreased upon addition of glycerol (-0.8 kcal mol-1 K-1 in the presence of glycerol compared to -2.3 kcal mol-1 K-1 in the absence of glycerol). Thus, the large negative delta Cp observed for the soluble receptor was a result of its removal from the membrane and was not due to the presence of glycerol.(ABSTRACT TRUNCATED AT 250 WORDS)     

Dynamic analysis of efficient heme rotation in myoglobin by NMR spectroscopy.

Sperm whale myoglobin was reconstituted with 1,4,5,8-tetramethylhemin. The hyperfine-shifted proton NMR signals from the prosthetic group exhibit remarkable pattern changes around 15 degrees C, while the globin resonances are normal to obey the Curie law. The NMR anomaly specifically observed for the heme signals suggests a slow to rapid rotational transition of the hemin about the iron-histidine bond. The temperature-dependent pattern changes were quantitatively analyzed by a dynamic NMR method. Two sets of analyses with the heme-methyl and pyrrole-proton lines consistently afforded delta H not equal to = 16.3 kcal/mol, delta S not equal to = 14.0 e.u., delta G not equal to = 12.1 kcal/mol at 298 K, and a frequency of 90 degrees heme rotation 5600 s-1 at 20 degrees C. The relatively large activation entropy suggests that structural rearrangements at the direct heme vicinity are involved and that efficient heme rotation is accomplished by a number of fluctuative local heme-globin contacts within a conserved crevice structure.     

Reduction of horse heart ferricytochrome c by bovine liver ferrocytochrome b5. Experimental and theoretical analysis.

The reduction of horse heart ferricytochrome c by the tryptic fragment of bovine liver cytochrome b5 and its dimethyl ester heme (DME)-substituted derivative has been studied as a function of ionic strength, pH, and temperature under solution conditions where the reaction is bimolecular. The rate constant for ferricytochrome c reduction by native ferrocytochrome b5 is 1.8 (+/- 0.2) x 10(7) M-1 s-1 (25 degrees C) with delta H++ = 7.5 (+/- 0.2) kcal/mol and delta S++ = -0.3 (+/- 0.6) eu (pH 7.0, I = 0.348 M). Under the same solution conditions, the reduction of ferricytochrome c by DME-ferrocytochrome b5 proceeds with a rate constant of 1.7 (+/- 0.1) x 10(7) M-1 s-1 with delta H++ = 7.9 (+/- 0.4) kcal/mol and delta S++ = 1 (+/- 1) eu. The rate constants for both reactions are strongly dependent on ionic strength. A detailed electrostatic analysis of the proteins has been performed. Two relatively simple Brownian dynamics simulation models predict rate constants for the reaction between the two native proteins that demonstrate a dependence on ionic strength similar to that observed experimentally. In one of these models, the proteins are treated as spheres with reactive surface patches that are defined by a 5 degrees cone generated about the dipole vector calculated for each protein and aligned with the presumed electron-transfer site near the partially exposed heme edge. The second model replaces the reactive patch assumption with an exponential distance dependence for the probability of reaction that permits estimation of a value for the distance-dependence factor alpha. Calculations with this latter model in combination with the aligned dipole assumption provide a reasonable approximation to the observed ionic strength dependence for the reaction and are consistent with a value of alpha = 1.2 A-1.     

Thermodynamics of phospholipid bilayer assembly

Thermodynamic properties of bilayer assembly have been obtained from measurements of the solubility of the sodium salt of dimyristoylphosphatidylglycerol (DMPG) in water. The standard free energy of bilayer assembly delta G degree a is shown to be RT 1n Xs + zF psi 0 where Xs is the mole fraction of dissolved lipid, F is the Faraday constant, z is the valence of the counterion (Na+), and psi 0 is the electrical double-layer potential of the ionized bilayer. The function d 1n Xs/dT was found to be discontinuous at 24 degrees C, the gel-liquid-crystal transition temperature (Tm) for DMPG. This function was unaffected when solubilities were measured in 0.001 M NaCl solutions; thus, psi 0 is constant in the experimental temperature interval (4-40 degrees C). Using a value of psi 0 = -180 mV [Eisenberg et al. (1979) Biochemistry 18, 5213-5223], and the temperature dependence of delta G degrees a, values for delta H degrees a and delta S degree a at 24 degrees C were calculated for the gel and liquid-crystal states of DMPG. For the gel, delta H degrees a and T delta S a are -26.2 and 12.7 kcal/mol, respectively; for the liquid-crystal, delta H degrees a and T delta S degrees a are -19.2 and -5.7 kcal/mol, respectively. The calculated value for the latent heat of the gel-liquid-crystal transition is 7 kcal/mol, in agreement with calorimetric measurements.     

Ligand binding equilibrium and kinetic measurements on the dimeric myoglobin of Busycon canaliculatum and the comparative ligand binding of diverse non-cooperative heme proteins

The rate constants and delta H degrees for the non-cooperative dimeric Busycon myoglobin are: oxygen, k' = 4.75 X 10(7) M-1 sec-1, k = 71 sec-1, and CO, l'= 3.46 X 10(5) M-1 sec-1, l = 0.0052 sec-1 at 20 degrees C, pH 7, delta H degrees = -3 kcal/mol for O2 and CO.2. Log-log plots of k vs K for oxygen and of l' vs L for CO binding for numerous non-cooperative hemoglobins and myoglobins point to a large steric influence of the protein on heme ligation reactions. Many of the proteins behave as "R" state for one ligand, but "T" for the other.     

Azide binding to the cytochrome c ferric heme octapeptide. A model for anion binding to the active site of high spin ferric heme proteins.

Equilibrium constants for the binding of azide to the ferric heme c octapeptide in 50% ethylene glycol 50% buffer were measured spectrophotometrically. The equilibrium constant for azide binding at 20 degrees C and pH* 7.4 is 29.2, which is approximately 3 to 4 orders of magnitude lower than that observed for azide binding to various ferric hemeproteins. The equilibrium constant was indepent of pH* in the range from 7 to 8. Equilibrium constants at several temperatures exhibited an apparent van't Hoff relationship yielding thermodynamic values of delta H0 = -26,100 J/mol (-6240 cal/mol) and delta S0 = -61.5 J/0K mol (-14.7 e.u.). Comparison of these values to the values for the heme proteins enables one to explain the differences in equiliberium constants in terms of differences in the polarity of the heme environments. The results are consistent with the concept that the oxygen affinity of heme complexes increases with the polarity of the heme environment. The data also suggest that an increase in the polarity of the heme environment should result in a corresponding increase in the susceptibility of ferrous heme dioxygen complexes toward oxidation by the dioxygen.     

The redox couple of the cytochrome c cyanide complex: the contribution of heme iron ligation to the structural stability, chemical reactivity, and physiological behavior of horse cytochrome c

Contrary to most heme proteins, ferrous cytochrome c does not bind ligands such as cyanide and CO. In order to quantify this observation, the redox potential of the ferric/ferrous cytochrome c-cyanide redox couple was determined for the first time by cyclic voltammetry. Its E0' was -240 mV versus SHE, equivalent to -23.2 kJ/mol. The entropy of reaction for the reduction of the cyanide complex was also determined. From a thermodynamic cycle that included this new value for the cyt c cyanide complex E0', the binding constant of cyanide to the reduced protein was estimated to be 4.7 x 10(-3) L M(-1) or 13.4 kJ/mol (3.2 kcal/mol), which is 48.1 kJ/mol (11.5 kcal/mol) less favorable than the binding of cyanide to ferricytochrome c. For coordination of cyanide to ferrocytochrome c, the entropy change was earlier experimentally evaluated as 92.4 J mol(-1) K(-1) (22.1 e.u.) at 25 K, and the enthalpy change for the same net reaction was calculated to be 41.0 kJ/mol (9.8 kcal/mol). By taking these results into account, it was discovered that the major obstacle to cyanide coordination to ferrocytochrome c is enthalpic, due to the greater compactness of the reduced molecule or, alternatively, to a lower rate of conformational fluctuation caused by solvation, electrostatic, and structural factors. The biophysical consequences of the large difference in the stabilities of the closed crevice structures are discussed.     

Crystal structures of ferrous and CO-, CN(-)-, and NO-bound forms of rat heme oxygenase-1 (HO-1) in complex with heme: structural implications for discrimination between CO and O2 in HO-1

Heme oxygenase (HO) catalyzes heme degradation by utilizing O(2) and reducing equivalents to produce biliverdin IX alpha, iron, and CO. To avoid product inhibition, the heme[bond]HO complex (heme[bond]HO) is structured to markedly increase its affinity for O(2) while suppressing its affinity for CO. We determined the crystal structures of rat ferrous heme[bond]HO and heme[bond]HO bound to CO, CN(-), and NO at 2.3, 1.8, 2.0, and 1.7 A resolution, respectively. The heme pocket of ferrous heme-HO has the same conformation as that of the previously determined ferric form, but no ligand is visible on the distal side of the ferrous heme. Fe[bond]CO and Fe[bond]CN(-) are tilted, whereas the Fe[bond]NO is bent. The structure of heme[bond]HO bound to NO is identical to that bound to N(3)(-), which is also bent as in the case of O(2). Notably, in the CO- and CN(-)-bound forms, the heme and its ligands shift toward the alpha-meso carbon, and the distal F-helix shifts in the opposite direction. These shifts allow CO or CN(-) to bind in a tilted fashion without a collision between the distal ligand and Gly139 O and cause disruption of one salt bridge between the heme and basic residue. The structural identity of the ferrous and ferric states of heme[bond]HO indicates that these shifts are not produced on reduction of heme iron. Neither such conformational changes nor a heme shift occurs on NO or N(3)(-) binding. Heme[bond]HO therefore recognizes CO and O(2) by their binding geometries. The marked reduction in the ratio of affinities of CO to O(2) for heme[bond]HO achieved by an increase in O(2) affinity [Migita, C. T., Matera, K. M., Ikeda-Saito, M., Olson, J. S., Fujii, H., Yoshimura, T., Zhou, H., and Yoshida, T. (1998) J. Biol. Chem. 273, 945-949] is explained by hydrogen bonding and polar interactions that are favorable for O(2) binding, as well as by characteristic structural changes in the CO-bound form.     

Characterization by NMR of the heme-myoglobin adduct formed during the reductive metabolism of BrCCl3. Covalent bonding of the proximal histidine to the ring I vinyl group

The reductive debromination of BrCCl3 by ferrous deoxymyoglobin leads to the covalent bonding of the prosthetic heme to the protein. We have previously shown, by the use of peptide mapping and mass spectrometry, that histidine residue 93 is covalently bound to the heme moiety. In the present study the structure of the heme adduct was more completely determined by 1H and 13C NMR techniques. We have found that the ring I vinyl group of the prosthetic heme was altered by the addition of a histidine imidazole nitrogen to the alpha-carbon and a CCl2 moiety to the beta-carbon. The electronic absorption spectra of the oxidized and reduced states of the altered heme-protein indicated that the heme-iron exists in a bis-histidine-ligated form. Analysis of the crystal structure of native myoglobin suggested that for the altered heme-protein, histidine residues 97 and 64 are ligated to the heme-iron and that residue 97 has replaced the native proximal histidine residue 93. These movements, in effect a "histidine shuffle" at the active site, may be responsible for the enhanced reducing activity of the altered protein.