Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study.

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or ß hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe-histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or ß chains, prior to the quaternary R-T and T-R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe-His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem. 2002, 98, 149). Experimentally, the νFeHis evolution is faster for ß than for α chains, and pump-probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 µs for ß but not for α hemes, an interval previously shown to be the first step in the R-T transition. In the α chain νFeHis dropped sharply at 20 µs, the final step in the R-T transition. The time courses are fully consistent with recent computational mapping of the R-T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.

Determination of ligand pathways in globins: apolar tunnels versus polar gates.

Although molecular dynamics simulations suggest multiple interior pathways for O(2) entry into and exit from globins, most experiments indicate well defined single pathways. In 2001, we highlighted the effects of large-to-small amino acid replacements on rates for ligand entry and exit onto the three-dimensional structure of sperm whale myoglobin. The resultant map argued strongly for ligand movement through a short channel from the heme iron to solvent that is gated by the distal histidine (His-64(E7)) near the solvent edge of the porphyrin ring. In this work, we have applied the same mutagenesis mapping strategy to the neuronal mini-hemoglobin from Cerebratulus lacteus (CerHb), which has a large internal tunnel from the heme iron to the C-terminal ends of the E and H helices, a direction that is 180° opposite to the E7 channel. Detailed comparisons of the new CerHb map with expanded results for Mb show unambiguously that the dominant (>90%) ligand pathway in CerHb is through the internal tunnel, and the major (>75%) ligand pathway in Mb is through the E7 gate. These results demonstrate that: 1) mutagenesis mapping can identify internal pathways when they exist; 2) molecular dynamics simulations need to be refined to address discrepancies with experimental observations; and 3) alternative pathways have evolved in globins to meet specific physiological demands.

Ligand migration in the apolar tunnel of Cerebratulus lacteus mini-hemoglobin.

The large apolar tunnel traversing the mini-hemoglobin from Cerebratulus lacteus (CerHb) has been examined by x-ray crystallography, ligand binding kinetics, and molecular dynamic simulations. The addition of 10 atm of xenon causes loss of diffraction in wild-type (wt) CerHbO(2) crystals, but Leu-86(G12)Ala CerHbO(2), which has an increased tunnel volume, stably accommodates two discrete xenon atoms: one adjacent to Leu-86(G12) and another near Ala-55(E18). Molecular dynamics simulations of ligand migration in wt CerHb show a low energy pathway through the apolar tunnel when Leu or Ala, but not Phe or Trp, is present at the 86(G12) position. The addition of 10-15 atm of xenon to solutions of wt CerHbCO and L86A CerHbCO causes 2-3-fold increases in the fraction of geminate ligand recombination, indicating that the bound xenon blocks CO escape. This idea was confirmed by L86F and L86W mutations, which cause even larger increases in the fraction of geminate CO rebinding, 2-5-fold decreases in the bimolecular rate constants for ligand entry, and large increases in the computed energy barriers for ligand movement through the apolar tunnel. Both the addition of xenon to the L86A mutant and oxidation of wt CerHb heme iron cause the appearance of an out Gln-44(E7) conformer, in which the amide side chain points out toward the solvent and appears to lower the barrier for ligand escape through the E7 gate. However, the observed kinetics suggest little entry and escape (≤ 25%) through the E7 pathway, presumably because the in Gln-44(E7) conformer is thermodynamically favored.

The stretching frequencies of bound alkyl isocyanides indicate two distinct ligand orientations within the distal pocket of myoglobin.

The FTIR spectra for alkyl isocyanides (CNRs) change from a single nu(CN) band centered at approximately 2175 cm(-1) to two peaks at approximately 2075 and approximately 2125 cm(-1) upon binding to sperm whale myoglobin (Mb). The low- and high-frequency peaks have been assigned to in and out conformations, respectively. In the in conformation, the ligand is pointing toward the protein interior, and the distal His64(E7) is in a closed position, donates a H-bond to the bound isocyano group, enhances back-bonding, and lowers the C-N bond order. In the out conformation, the ligand side chain points toward solvent through a channel opened by outward rotation of His64. Loss of positive polarity near the binding site causes an increase in C-N bond order. Support for this interpretation is threefold: (1) similar shifts to lower frequency occur for MbCO complexes when H-bond donation from His64(E7) occurs; (2) only one peak at approximately 2125 cm(-1), indicative of an apolar environment, is observed for CNRs bound to H64A or H64L Mb mutants or to chelated protoheme in soap micelles; and (3) the fraction of in conformation based on FTIR spectra correlates strongly with the fraction of geminate recombination after nanosecond laser photolysis. The in alkyl side chain conformation causes the photodissociated ligand to be "stuck" in the distal pocket, promoting internal rebinding, whereas the out conformation inhibits geminate recombination because part of the ligand is already in an open E7 channel, poised for rapid escape.

Alkyl isocyanides serve as transition state analogues for ligand entry and exit in myoglobin.

Alkyl isocyanides (CNRs) identify pathways for diatomic ligand movement into and out of Mb, with their side chains acting as transition state analogues. The bound alkyl groups point either into the back of the distal pocket (in conformation, nu(CN) approximately 2070-2090 cm(-1)), which allows hydrogen bond donation from His64(E7) to the isocyano group, or toward solvent through an open His(E7) channel (out conformation, nu(CN) approximately 2110-2130 cm(-1)), which prevents polar interactions with the isocyano atoms. Fractions of the in conformer (F(in)) were measured by FTIR spectroscopy for methyl through n-pentyl isocyanide bound to a series of 20 different distal pocket mutants of sperm whale myoglobin and found to be governed by the ease of rotation of the His(E7) side chain, distal pocket volume and steric interactions, and, for the longer isocyanides, the unfavorable hydrophobic effect of placing their terminal carbon atoms into the solvent phase in the out conformation. There are strong correlations between the fraction of in conformer, F(in), for long-chain MbCNR complexes measured by FTIR spectroscopy, the fraction of geminate recombination of photodissociated O(2), and the bimolecular rates of O(2) entry into the distal pocket. These correlations indicate that alkyl isocyanides serve as transition state analogues for the movement of O(2) into and out of the binding pocket of Mb.

Straight-chain alkyl isocyanides open the distal histidine gate in crystal structures of myoglobin .

Crystal structures of methyl, ethyl, propyl, and butyl isocyanide bound to sperm whale myoglobin (Mb) reveal two major conformations. In the in conformer, His(E7) is in a "closed" position, forcing the ligand alkyl chain to point inward. In the out conformer, His(E7) is in an "open" position, allowing the ligand side chain to point outward. A progressive increase in the population of the out conformer is observed with increasing ligand length in P2(1) crystals of native Mb at pH 7.0. This switch from in to out with increasing ligand size also occurs in solution as measured by the decrease in the relative intensity of the low-frequency ( approximately 2075 cm(-1)) versus high-frequency ( approximately 2125 cm(-1)) isocyano bands. In contrast, all four isocyanides in P6 crystals of wild-type recombinant Mb occupy the in conformation. However, mutating either His64 to Ala, creating a "hole" to solvent, or Phe46 to Val, freeing rotation of His64, causes bound butyl isocyanide to point completely outward in P6 crystals. Thus, the unfavorable hindrance caused with crowding a large alkyl side chain into the distal pocket appears to be roughly equal to that for pushing open the His(E7) gate and is easily affected by crystal packing. This structural conclusion supports the "side path" kinetic mechanism for O(2) release, in which the dissociated ligand first moves toward the protein interior and then encounters steric resistance, which is roughly equal to that for escaping to solvent through the His(E7) channel.

The structure and NO binding properties of the nitrophorin-like heme-binding protein from Arabidopsis thaliana gene locus At1g79260.1.

The protein from Arabidopsis thaliana gene locus At1g79260.1 is comprised of 166-residues and is of previously unknown function. Initial structural studies by the Center for Eukaryotic Structural Genomics (CESG) suggested that this protein might bind heme, and consequently, the crystal structures of apo and heme-bound forms were solved to near atomic resolution of 1.32 A and 1.36 A, respectively. The rate of hemin loss from the protein was measured to be 3.6 x 10(-5) s(-1), demonstrating that it binds heme specifically and with high affinity. The protein forms a compact 10-stranded beta-barrel that is structurally similar to the lipocalins and fatty acid binding proteins (FABPs). One group of lipocalins, the nitrophorins (NP), are heme proteins involved in nitric oxide (NO) transport and show both sequence and structural similarity to the protein from At1g79260.1 and two human homologues, all of which contain a proximal histidine capable of coordinating a heme iron. Rapid-mixing and laser photolysis techniques were used to determine the rate constants for carbon monoxide (CO) binding to the ferrous form of the protein (k'(CO) = 0.23 microM(-1) s(-1), k(CO) = 0.050 s(-1)) and NO binding to the ferric form (k'(NO) = 1.2 microM(-1) s(-1), k(NO) = 73 s(-1)). Based on both structural and functional similarity to the nitrophorins, we have named the protein nitrobindin and hypothesized that it plays a role in NO transport. However, one of the two human homologs of nitrobindin contains a THAP domain, implying a possible role in apoptosis. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Review: correlations between oxygen affinity and sequence classifications of plant hemoglobins.

Plants express three phylogenetic classes of hemoglobins (Hb) based on sequence analyses. Class 1 and 2 Hbs are full-length globins with the classical eight helix Mb-like fold, whereas Class 3 plant Hbs resemble the truncated globins found in bacteria. With the exception of the specialized leghemoglobins, the physiological functions of these plant hemoglobins remain unknown. We have reviewed and, in some cases, measured new oxygen binding properties of a large number of Class 1 and 2 plant nonsymbiotic Hbs and leghemoglobins. We found that sequence classification correlates with distinct extents of hexacoordination with the distal histidine and markedly different overall oxygen affinities and association and dissociation rate constants. These results suggest strong selective pressure for the evolution of distinct physiological functions. The leghemoglobins evolved from the Class 2 globins and show no hexacoordination, very high rates of O(2) binding ( approximately 250 muM(-1) s(-1)), moderately high rates of O(2) dissociation ( approximately 5-15 s(-1)), and high oxygen affinity (K(d) or P(50) approximately 50 nM). These properties both facilitate O(2) diffusion to respiring N(2) fixing bacteria and reduce O(2) tension in the root nodules of legumes. The Class 1 plant Hbs show weak hexacoordination (K(HisE7) approximately 2), moderate rates of O(2) binding ( approximately 25 muM(-1) s(-1)), very small rates of O(2) dissociation ( approximately 0.16 s(-1)), and remarkably high O(2) affinities (P(50) approximately 2 nM), suggesting a function involving O(2) and nitric oxide (NO) scavenging. The Class 2 Hbs exhibit strong hexacoordination (K(HisE7) approximately 100), low rates of O(2) binding ( approximately 1 muM(-1) s(-1)), moderately low O(2) dissociation rate constants ( approximately 1 s(-1)), and moderate, Mb-like O(2) affinities (P(50) approximately 340 nM), perhaps suggesting a sensing role for sustained low, micromolar levels of oxygen.

The mechanism of direct heme transfer from the streptococcal cell surface protein Shp to HtsA of the HtsABC transporter.

The heme-binding proteins Shp and HtsA are part of the heme acquisition machinery found in Streptococcus pyogenes. The hexacoordinate heme (Fe(II)-protoporphyrin IX) or hemochrome form of holoShp (hemoShp) is stable in air in Tris-HCl buffer, pH 8.0, binds to apoHtsA with a K(d) of 120 +/- 18 microm, and transfers its heme to apoHtsA with a rate constant of 28 +/- 6s(-1) at 25 degrees C, pH 8.0. The hemoHtsA product then autoxidizes to the hexacoordinate hemin (Fe(III)-protoporphyrin IX) or hemichrome form (hemiHtsA) with an apparent rate constant of 0.017 +/- 0.002 s(-1). HemiShp also rapidly transfers hemin to apoHtsA through a hemiShp.apoHtsA complex (K(d) = 48 +/- 7 microM) at a rate approximately 40,000 times greater than the rate of simple hemin dissociation from hemiShp into solvent (k(transfer) = 43 +/- 3s(-1) versus k(-hemin) = 0.0003 +/- 0.00006 s(-1)). The rate constants for hemin binding to and dissociation from HtsA (k'(hemin) approximately 80 microm(-1) s(-1), k(-hemin) = 0.0026 +/- 0.0002 s(-1)) are 50- and 10-fold greater than the corresponding rate constants for Shp (k(hemin) approximately 1.6 microM(-1) s(-1), k(-hemin) = 0.0003 s(-1)), which implies that HtsA has a more accessible active site. However, the affinity of apoHtsA for hemin (k(hemin) approximately 31,000 microm(-1)) is roughly 5-fold greater than that of apoShp (k(hemin) approximately 5,300 microM(-1)), accounting for the net transfer from Shp to HstA. These results support a direct, rapid, and affinity-driven mechanism of heme and hemin transfer from the cell surface receptor Shp to the ATP-binding cassette transporter system.

Tyrosine B10 inhibits stabilization of bound carbon monoxide and oxygen in soybean leghemoglobin.

Detailed comparisons of the carbon monoxide FTIR spectra and ligand-binding properties of a library of E7, E11, and B10 mutants indicate significant differences in the role of electrostatic interactions in the distal pockets of wild-type sperm whale myoglobin and soybean leghemoglobin. In myoglobin, strong hydrogen bonds from several closely related conformations of the distal histidine (His(E7)) side chain preferentially stabilize bound oxygen. In leghemoglobin, the imidazole side chain of His(E7) is confined to a single conformation, which only weakly hydrogen bonds to bound ligands. The phenol side chain of Tyr(B10) appears to "fix" the position of His(E7), probably by donating a hydrogen bond to the Ndelta atom of the imidazole side chain. The proximal pocket of leghemoglobin is designed to favor strong coordination bonds between the heme iron and axial ligands. Thus, high oxygen affinity in leghemoglobin is established by a favorable staggered geometry of the proximal histidine. The interaction between His(E7) and Tyr(B10) prevents overstabilization of bound oxygen. If hydrogen bonding from His(E7) were as strong as it is in mammalian myoglobin, the resultant ultrahigh affinity of leghemoglobin would prevent oxygen transport in root nodules.