Distal pocket polarity in ligand binding to myoglobin: structural and functional characterization of a threonine68(E11) mutant

Sequencing of the Lumbricus terrestris genome reveals degeneracy in its erythrocruorin genes

The erythrocruorin of Lumbricus terrestris (LtEc) is a relatively large macromolecular assembly that consists of at least four different hemoglobin subunits (A, B, C, and D) and four linker subunits (L1, L2, L3, and L4). The complexity and stability of this large structure make LtEc an attractive hemoglobin-based oxygen carrier that could potentially be used as a substitute for donated red blood cells. However, the sequences of the LtEc subunit sequences must be determined before a scalable recombinant expression platform can be developed. The goal of this study was to sequence the L. terrestris genome to identify the complete sequences of the LtEc subunit genes. Our results revealed multiple homologous genes for each subunit (e.g., two homologous A globin genes; A1 and A2), with the exception of the L4 linker. Some of the homologous genes encoded identical peptide sequences (C1 and C2, L1a and L1b), while cDNA and mass spectrometry experiments revealed that some of the homologs are not expressed (e.g., A2). In contrast, multiple sequences for the B, D, L2, and L4 subunits were detected in LtEc samples. These observations reveal novel degeneracy in LtEc and other annelids, along with some new revisions to its previously published peptide sequences.

A Heme Pocket Aromatic Quadrupole Modulates Gas Binding to Cytochrome c'-β: Implications for NO Sensors

The structural basis by which gas-binding heme proteins control their interactions with NO, CO, and O₂, is fundamental to enzymology, biotechnology and human health. Cytochromes c´ (cyts c´) are a group of putative NO-binding heme proteins that fall into two families: the well characterised four alpha helix bundle fold (cyts c´-α) and an unrelated family with a largely beta sheet fold (cyts c´-β) resembling that of cytochromes P460. A recent structure of cyt c´-β from Methylococcus capsulatus Bath (McCP-β) revealed two heme pocket phenylalanine residues (Phe 32 and Phe 61) positioned near the distal gas binding site. This feature, dubbed the "Phe cap", is highly conserved within the sequences of other cyts c´-β, but is absent in their close homologues, the hydroxylamine oxidizing cytochromes P460, although some do contain a single Phe residue. Here we report an integrated structural, spectroscopic, and kinetic characterization of McCP-β complexes with diatomic gases, focusing on the interaction of the Phe cap with NO and CO. Significantly, crystallographic and resonance Raman data show that orientation of the electron rich aromatic ring face of Phe 32 towards distally-bound NO or CO is associated with weakened backbonding and higher off rates. Moreover, we propose that an aromatic quadrupole also contributes to the unusually weak backbonding reported for some heme-based gas sensors, including the mammalian NO-sensor, soluble guanylate cyclase (sGC). Collectively, this study sheds light on the influence of highly conserved distal Phe residues on heme-gas complexes of cytochrome c'-β, including the potential for aromatic quadrupoles to modulate NO and CO binding in other heme proteins.

Apoglobin Stability Is the Major Factor Governing both Cell-free and in Vivo Expression of Holomyoglobin

Expression levels in animal muscle tissues and in Escherichia coli vary widely for naturally occurring mammalian myoglobins (Mb). To explore this variation, we developed an in vitro transcription and wheat germ extract-based translation assay to examine quantitatively the factors that govern expression of holoMb. We constructed a library of naturally occurring Mbs from two terrestrial and four deep-diving aquatic mammals and three distal histidine mutants designed to enhance apoglobin stability but decrease hemin affinity. A strong linear correlation is observed between cell-free expression levels of holo-metMb variants and their corresponding apoglobin stabilities, which were measured independently by guanidine HCl-induced unfolding titrations using purified proteins. In contrast, there is little dependence of expression on hemin affinity. Our results confirm quantitatively that deep diving mammals have highly stable Mbs that express to higher levels in animal myocytes, E. coli, and the wheat germ cell-free system than Mbs from terrestrial mammals. Our theoretical analyses show that the rate of aggregation of unfolded apoMb is very large, and as a result, the key factor for high level expression of holoMb, and presumably other heme proteins, is an ultra high fraction of folded, native apoglobin that is capable of rapidly binding hemin. This fraction is determined by the overall equilibrium folding constant and not hemin affinity. These results also demonstrate that the cell-free transcription/translation system can be used as a high throughput platform to screen for apoglobin stability without the need to generate large amounts of protein for in vitro unfolding measurements.

How do heme-protein sensors exclude oxygen? Lessons learned from cytochrome c', Nostoc puntiforme heme nitric oxide/oxygen-binding domain, and soluble guanylyl cyclase

SIGNIFICANCE

Ligand selectivity for dioxygen (O(2)), carbon monoxide (CO), and nitric oxide (NO) is critical for signal transduction and is tailored specifically for each heme-protein sensor. Key NO sensors, such as soluble guanylyl cyclase (sGC), specifically recognized low levels of NO and achieve a total O(2) exclusion. Several mechanisms have been proposed to explain the O(2) insensitivity, including lack of a hydrogen bond donor and negative electrostatic fields to selectively destabilize bound O(2), distal steric hindrance of all bound ligands to the heme iron, and restriction of in-plane movements of the iron atom.

RECENT ADVANCES

Crystallographic structures of the gas sensors, Thermoanaerobacter tengcongensis heme-nitric oxide/oxygen-binding domain (Tt H-NOX(1)) or Nostoc puntiforme (Ns) H-NOX, and measurements of O(2) binding to site-specific mutants of Tt H-NOX and the truncated β subunit of sGC suggest the need for a H-bonding donor to facilitate O(2) binding.

CRITICAL ISSUES

However, the O(2) insensitivity of full length sGC with a site-specific replacement of isoleucine by a tyrosine on residue 145 and the very slow autooxidation of Ns H-NOX and cytochrome c' suggest that more complex mechanisms have evolved to exclude O(2) but retain high affinity NO binding. A combined graphical analysis of ligand binding data for libraries of heme sensors, globins, and model heme shows that the NO sensors dramatically inhibit the formation of six-coordinated NO, CO, and O(2) complexes by direct distal steric hindrance (cyt c'), proximal constraints of in-plane iron movement (sGC), or combinations of both following a sliding scale rule. High affinity NO binding in H-NOX proteins is achieved by multiple NO binding steps that produce a high affinity five-coordinate NO complex, a mechanism that also prevents NO dioxygenation.

FUTURE DIRECTIONS

Knowledge advanced by further extensive test of this "sliding scale rule" hypothesis should be valuable in guiding novel designs for heme based sensors.

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.

Trema and parasponia hemoglobins reveal convergent evolution of oxygen transport in plants

All plants contain hemoglobins that fall into distinct phylogenetic classes. The subset of plants that carry out symbiotic nitrogen fixation expresses hemoglobins that scavenge and transport oxygen to bacterial symbiotes within root nodules. These "symbiotic" oxygen transport hemoglobins are distinct in structure and function from the nonoxygen transport ("nonsymbiotic") Hbs found in all plants. Hemoglobins found in two closely related plants present a paradox concerning hemoglobin structure and function. Parasponia andersonii is a nitrogen-fixing plant that expresses a symbiotic hemoglobin (ParaHb) characteristic of oxygen transport hemoglobins in having a pentacoordinate ferrous heme iron, moderate oxygen affinity, and a relatively rapid oxygen dissociation rate constant. A close relative that does not fix nitrogen, Trema tomentosa, expresses hemoglobin (TremaHb) sharing 93% amino acid identity to ParaHb, but its phylogeny predicts a typical nonsymbiotic hemoglobin with a hexacoordinate heme iron, high oxygen affinity, and slow oxygen dissociation rate constant. Here we characterize heme coordination and oxygen binding in TremaHb and ParaHb to investigate whether or not two hemoglobins with such high sequence similarity are actually so different in functional behavior. Our results indicate that the two proteins resemble nonsymbiotic hemoglobins in the ferric oxidation state and symbiotic hemoglobins in the ferrous oxidation state. They differ from each other only in oxygen affinity and oxygen dissociation rate constants, two factors key to their different functions. These results demonstrate distinct mechanisms for convergent evolution of oxygen transport in different phylogenetic classes of plant hemoglobins.

Optical detection of disordered water within a protein cavity

Internal water molecules are important to protein structure and function, but positional disorder and low occupancies can obscure their detection by X-ray crystallography. Here, we show that water can be detected within the distal cavities of myoglobin mutants by subtle changes in the absorbance spectrum of pentacoordinate heme, even when the presence of solvent is not readily observed in the corresponding crystal structures. A well-defined, noncoordinated water molecule hydrogen bonded to the distal histidine (His64) is seen within the distal heme pocket in the crystal structure of wild type (wt) deoxymyoglobin. Displacement of this water decreases the rate of ligand entry into wt Mb, and we have shown previously that the entry of this water is readily detected optically after laser photolysis of MbCO complexes. However, for L29F and V68L Mb no discrete positions for solvent molecules are seen in the electron density maps of the crystal structures even though His64 is still present and slow rates of ligand binding indicative of internal water are observed. In contrast, time-resolved perturbations of the visible absorption bands of L29F and V68L deoxyMb generated after laser photolysis detect the entry and significant occupancy of water within the distal pockets of these variants. Thus, the spectral perturbation of pentacoordinate heme offers a potentially robust system for measuring nonspecific hydration of the active sites of heme proteins.

A comparative study on axial coordination and ligand binding in ferric mini myoglobin and horse heart myoglobin

The absorption and resonance Raman spectra and the azide binding kinetics of ferric horse heart myoglobin (Mb) and mini myoglobin (a chemically truncated form of horse heart Mb containing residues 32-139) have been compared. The steady-state spectra show that an additional six-coordinated low-spin form (not present in entire horse heart Mb, which is purely six-coordinated high spin) predominates in mini Mb. The distal histidine is possibly the sixth ligand in this species. The presence of two species corresponds to a kinetic biphasicity for mini Mb that is not observed for horse heart Mb. Azide binds to horse heart Mb much more slowly than to sperm whale Mb. This difference may result from a sterically hindered distal pocket in horse heart Mb. In both cases, the rate constants level off at high azide concentrations, implying the existence of a rate-limiting step (likely referable to the dissociation of the axial sixth ligand). The faster rate constant of mini Mb is similar to that of sperm whale Mb, whereas the slower one is similar to that of entire horse heart Mb.

Diatomic ligand discrimination by the heme oxygenases from Neisseria meningitidis and Pseudomonas aeruginosa

Heme oxygenases have an increased binding affinity for O2 relative to CO. Such discrimination is critical to the function of HO enzymes because one of the main products of heme catabolism is CO. Kinetic studies of mammalian and bacterial HO proteins reveal a significant decrease in the dissociation rate of O2 relative to other heme proteins such as myoglobin. Here we report the kinetic rate constants for the binding of O2 and CO by the heme oxygenases from Neisseria meningitidis (nmHO) and Pseudomonas aeruginosa (paHO). A combination of stopped-flow kinetic and laser flash photolysis experiments reveal that nmHO and paHO both maintain a similar degree of ligand discrimination as mammalian HO-1 and the HO from Corynebacterium diphtheriae. However, in addition to the observed decrease in dissociation rate for O2 by both nmHO and paHO, kinetic analyses show an increase in dissociation rate for CO by these two enzymes. The crystal structures of nmHO and paHO both contain significant differences from the mammalian HO-1 and bacterial C. diphtheriae HO structures, which suggests a structural basis for ligand discrimination in nmHO and paHO.

Mechanisms of heme protein-mediated lipid oxidation using hemoglobin and myoglobin variants in raw and heated washed muscle

The hemoglobin variant rHb 0.1, which possesses a decreased ability to form subunits, stimulated lipid oxidation in washed fish muscle less effectively as compared to wild-type hemoglobin (rHb 0.0). This could be due to the lower hemin affinity and more rapid autoxidation rate of subunits as compared to tetramers. To differentiate between hemin affinity and autoxidation effects, ferrous V68T Mb was compared to ferrous wild-type myoglobin (WT Mb). WT Mb has a more rapid hemin loss rate (25-fold) than does V68T, while V68T autoxidized more rapidly than did WT Mb (60-fold). Ferrous WT Mb promoted TBARS and lipid peroxide formation more rapidly than did ferrous V68T (p < 0.01). This indicated hemin loss rate was more critical in determining onset of lipid oxidation as compared to autoxidation rate. Hemin alone was capable of stimulating lipid oxidation. Albumin enhanced the ability of hemin to promote lipid oxidation. MetMb promoted lipid oxidation more effectively than did ferrous Mb, which could be due to the lower hemin affinity of metMb as compared to that of ferrous Mb. EDTA, an iron chelator, had no effect on the rate or extent of lipid oxidation mediated by Mb in the cooked system. Variants with a 975-fold range of hemin affinities promoted lipid oxidation with equivalent efficacy in cooked washed cod contrary to results in uncooked washed cod. The cooking temperatures apparently denature the globin and release hemin reactant to such an extent that the impact of hemin affinity on lipid oxidation observed in the raw state is negated in the cooked state. These studies collectively suggest released hemin is of primary importance in promoting lipid oxidation in raw and cooked washed fish muscle.

Studies with myoglobin variants indicate that released hemin is the primary promoter of lipid oxidation in washed fish muscle

Variants of sperm whale myoglobin (Mb) were used to assess the mechanism of heme protein-mediated lipid oxidation in washed cod muscle. A myoglobin variant with high hemin affinity (V68T) was an exceptionally poor promoter of lipid oxidation, while a Mb variant with low hemin affinity (H97A) was a potent promoter of lipid oxidation. V68T releases hemin slowly due to the ability of threonine to hydrogen bond with coordinated water and the distal histidine within the heme crevice. H97A rapidly releases hemin because the relatively small alanine residue creates a channel for water to easily enter the heme crevice which weakens the covalent linkage of hemin to the proximal histidine. A variant sensitive to heme degradation (L29F/H64Q) was a weaker promoter of lipid oxidation compared to wild-type Mb. This suggests that degrading the heme ring and releasing iron decreased the ability of Mb to promote lipid oxidation. Free radicals resulting from hemin-mediated decomposition of lipid hydroperoxides have the capacity to propagate lipid oxidation and degrade hemin catalyst. This may explain why heme proteins behave as reactants rather than "catalysts" of lipid oxidation in washed cod. Collectively these studies strongly suggest that released hemin is the critical entity that drives heme protein-mediated lipid oxidation in washed fish muscle.

On the Relationship of Coral Allene Oxide Synthase to Catalase

A heme domain of coral allene oxide synthase (cAOS) catalyzes the formation of allene oxide from fatty acid hydroperoxide. Although cAOS has a similar heme active site to that of catalase, cAOS is completely lacking in catalase activity. A close look at the hydrogen-bonding possibilities around the distal His in cAOS suggested that the imidazole ring is rotated by 180 degrees relative to that of catalase because of the hydrogen bond between Thr-66 and the distal His-67. This could contribute to the functional differences between cAOS and catalase, and to examine this possibility, we mutated Thr-66 in cAOS to Val, the corresponding residue in catalase. In contrast to the complete absence of catalase activity in wild type (WT) cAOS, T66V had a modest catalase activity. On the other hand, the mutation suppressed the native enzymatic activity of the formation of allene oxide to 14% of that of WT cAOS. In the resonance Raman spectrum, whereas WT cAOS has only a 6-coordinate/high spin heme, T66V has a 5-coordinate/high spin heme as a minor species. Because catalase adopts a 5-coordinate/high spin structure, probably the 5-coordinate/high spin portion of T66V showed the catalase activity. Furthermore, in accord with the fact that the CN affinity of catalase is higher than that of WT cAOS, the CN affinity of T66V was 8-fold higher than that of WT cAOS, indicating that the mutation could mimic the heme active site in catalase. We, therefore, propose that the hydrogen bond between Thr-66 and distal His-67 could modulate the orientation of distal His, thereby regulating the enzymatic activity in cAOS.

Water and ligand entry in myoglobin: assessing the speed and extent of heme pocket hydration after CO photodissociation

A previously undescribed spectrokinetic assay for the entry of water into the distal heme pocket of wild-type and mutant myoglobins is presented. Nanosecond photolysis difference spectra were measured in the visible bands of sperm whale myoglobin as a function of distal pocket mutation and temperature. A small blue shift in the 560-nm deoxy absorption peak marked water entry several hundred nanoseconds after CO photodissociation. The observed rate suggests that water entry is rate-limited by the escape of internal dissociated CO. The heme pocket hydration and geminate recombination yields were found to be the primary factors controlling the overall bimolecular association rate constants for CO binding to the mutants studied. The kinetic analysis provides estimates of 84%, 60%, 40%, 0%, and 99% for the steady-state hydrations of wild-type, H64Q, H64A, H64L, and V68F deoxymyoglobin, respectively. The second-order rate constants for CO and H(2)O entry into the empty distal pocket of myoglobin are markedly different, 8 x 10(7) and 2 x 10(5) M(-1).s(-1), respectively, suggesting that hydrophobic partitioning of the apolar gas from the aqueous phase into the relatively apolar protein interior lowers the free energy barrier for CO entry.

Human myoglobin recognition of oxygen: dynamics of the energy landscape

Femtosecond to nanosecond dynamics of O(2) rebinding to human WT myoglobin and its mutants, V68F and I107F, have been studied by using transient absorption. The results are compared with NO rebinding. Even though the immediate environment around the heme binding site is changed by the mutations, the picosecond geminate rebinding of oxygen is at most minimally affected. On the other hand, the V68F (E11) mutation causes drastic differences in rebinding on the nanosecond time scale, whereas the effect of the I107F (G8) mutation remains relatively small within our 10-ns time window. Unlike traditional homogeneous kinetics and molecular dynamics collisional simulations, we propose a "bifurcation model" for populations of directed and undirected dynamics on the ultrafast time scale, reflecting the distribution of initial protein conformations. The major mutation effect occurs on the time scale on which global protein conformational change is possible, consistent with transitions between the conformations of directed and undirected population playing a role in the O(2) binding. We discuss the relevance of these findings to the bimolecular function of the protein.

Viscosity-dependent Relaxation Significantly Modulates the Kinetics of CO Recombination in the Truncated Hemoglobin TrHbN from Mycobacterium tuberculosis

Kinetic traces were generated for the nanosecond and slower rebinding of photodissociated CO to trHbN in solution and in porous sol-gel matrices as a function of viscosity, conformation, and mutation. TrHbN is one of the two truncated hemoglobins from Mycobacterium tuberculosis. The kinetic traces were analyzed in terms of three distinct phases. These three phases are ascribed to rebinding: (i) from the distal heme pocket, (ii) from the adjacent apolar tunnel prior to conformational relaxation, and (iii) from the apolar tunnel subsequent to conformational relaxation. The fractional content of each of these phases was shown to be a function of the viscosity and, in the case of the sol-gel-encapsulated samples, sample preparation history. The observed kinetic patterns support a model consisting of the following elements: (i) the viscosity and conformation-sensitive dynamics of the Tyr(B10) side chain facilitate diffusion of the dissociated ligand from the distal heme pocket into the adjacent tunnel; (ii) the distal heme pocket architecture determines ligand access from the tunnel back to the heme iron; (iii) the distal heme pocket architecture is governed by a ligand-dependent hydrogen bonding network that limits the range of accessible side chain positions; and (iv) the apolar tunnel linking the heme site to the solvent biases the competition between water and ligand for occupancy of the vacated polar distal heme pocket greatly toward the nonpolar ligand. Implications of these finding with respect to biological function are discussed.

Thr-E11 regulates O2 affinity in Cerebratulus lacteus mini-hemoglobin

The mini-hemoglobin from Cerebratulus lacteus (CerHb) belongs to a class of globins containing the polar Tyr-B10/Gln-E7 amino acid pair that normally causes low rates of O2 dissociation and ultra-high O2 affinity, which suggest O2 sensing or NO scavenging functions. CerHb, however, has high rates of O2 dissociation (kO2 = 200-600 s(-1)) and moderate O2 affinity (KO2) approximately 1 microm(-1)) as a result of a third polar amino acid in its active site, Thr-E11. When Thr-E11 is replaced by Val, kO2 decreases 1000-fold and KO2 increases 130-fold at pH 7.0, 20 degrees C. The mutation also shifts the stretching frequencies of both heme-bound and photodissociated CO, indicating marked changes of the electrostatic field at the active site. The crystal structure of Thr-E11 --> Val CerHbO2 at 1.70 A resolution is almost identical to that of the wild-type protein (root mean square deviation of 0.12 A). The dramatic functional and spectral effects of the Thr-E11 --> Val mutation are due exclusively to changes in the hydrogen bonding network in the active site. Replacing Thr-E11 with Val "frees" the Tyr-B10 hydroxyl group to rotate toward and donate a strong hydrogen bond to the heme-bound ligand, causing a selective increase in O2 affinity, a decrease of the rate coefficient for O2 dissociation, a 40 cm(-1) decrease in nuCO of heme-bound CO, and an increase in ligand migration toward more remote intermediate sites.

Crystal structure of the dioxygen-bound heme oxygenase from Corynebacterium diphtheriae: implications for heme oxygenase function

HmuO, a heme oxygenase of Corynebacterium diphtheriae, catalyzes degradation of heme using the same mechanism as the mammalian enzyme. The oxy form of HmuO, the precursor of the catalytically active ferric hydroperoxo species, has been characterized by ligand binding kinetics, resonance Raman spectroscopy, and x-ray crystallography. The oxygen association and dissociation rate constants are 5 microm(-1) s(-1) and 0.22 s(-1), respectively, yielding an O(2) affinity of 21 microm(-1), which is approximately 20 times greater than that of mammalian myoglobins. However, the affinity of HmuO for CO is only 3-4-fold greater than that for mammalian myoglobins, implying the presence of strong hydrogen bonding interactions in the distal pocket of HmuO that preferentially favor O(2) binding. Resonance Raman spectra show that the Fe-O(2) vibrations are tightly coupled to porphyrin vibrations, indicating the highly bent Fe-O-O geometry that is characteristic of the oxy forms of heme oxygenases. In the crystal structure of the oxy form the Fe-O-O angle is 110 degrees, the O-O bond is pointed toward the heme alpha-meso-carbon by direct steric interactions with Gly-135 and Gly-139, and hydrogen bonds occur between the bound O(2) and the amide nitrogen of Gly-139 and a distal pocket water molecule, which is a part of an extended hydrogen bonding network that provides the solvent protons required for oxygen activation. In addition, the O-O bond is orthogonal to the plane of the proximal imidazole side chain, which facilitates hydroxylation of the porphyrin alpha-meso-carbon by preventing premature O-O bond cleavage.

Structural dynamics of myoglobin: ligand migration and binding in valine 68 mutants

We have combined Fourier transform infrared/temperature derivative (FTIR-TDS) spectroscopy at cryogenic temperatures and flash photolysis at ambient temperature to examine the effects of polar and bulky amino acid replacements of the highly conserved distal valine 68 in sperm whale myoglobin. In FTIR-TDS experiments, the CO ligand can serve as an internal voltmeter that monitors the local electrostatic field not only at the active site but also at intermediate ligand docking sites. Mutations of residue 68 alter size, shape, and electric field of the distal pocket, especially in the vicinity of the primary docking site (state B). As a consequence, the infrared bands associated with the ligand at site B are shifted. The effect is most pronounced in mutants with large aromatic side chains. Polar side chains (threonine or serine) have only little effect on the peak frequencies. Ligands that migrate toward more remote sites C and D give rise to IR bands with altered frequencies. TDS experiments separate the photoproducts according to their recombination temperatures. The rates and extent of ligand migration among internal cavities at cryogenic temperatures can be used to interpret geminate and bimolecular O2 and CO recombination at room temperature. The kinetics of geminate recombination can be explained by steric arguments alone, whereas both the polarity and size of the position 68 side chain play major roles in regulating bimolecular ligand binding from the solvent.

Distal heme pocket regulation of ligand binding and stability in soybean leghemoglobin

Leghemoglobins facilitate diffusion of oxygen through root tissue to a bacterial terminal oxidase in much the same way that myoglobin transports oxygen from blood to muscle cell mitochondria. Leghemoglobin serves an additional role as an oxygen scavenger to prevent inhibition of nitrogen fixation. For this purpose, the oxygen affinity of soybean leghemoglobin is 20-fold greater than myoglobin, resulting from an 8-fold faster association rate constant combined with a 3-fold slower dissociation rate constant. Although the biochemical mechanism used by myoglobin to bind oxygen has been described in elegant detail, an explanation for the difference in affinity between these two structurally similar proteins is not obvious. The present work demonstrates that, despite their similar structures, leghemoglobin uses methods different from myoglobin to regulate ligand affinity. Oxygen and carbon monoxide binding to a comprehensive set of leghemoglobin distal heme pocket mutant proteins in comparison to their myoglobin counterparts has revealed some of these mechanisms. The "distal histidine" provides a crucial hydrogen bond to stabilize oxygen in myoglobin but has little effect on bound oxygen in leghemoglobin and is retained mainly for reasons of protein stability and prevention of heme loss. Furthermore, soybean leghemoglobin uses an unusual combination of HisE7 and TyrB10 to sustain a weak stabilizing interaction with bound oxygen. Thus, the leghemoglobin distal heme pocket provides a much lower barrier to oxygen association than occurs in myoglobin and oxygen dissociation is regulated from the proximal heme pocket.

Ligand migration in human myoglobin: steric effects of isoleucine 107(G8) on O(2) and CO binding

To investigate the ligand pathway in myoglobin, some mutant myoglobins, in which one of the amino acid residues constituting a putative ligand-docking site, Ile107, is replaced by Ala, Val, Leu, or Phe, were prepared and their structural and ligand binding properties were characterized. The kinetic barrier for the ligand entry to protein inside was lowered by decreasing the side-chain volume at position 107, indicating that the bulky side chain interferes with the formation of the activation state for the ligand migration and the free space near position 107 would be filled with the ligand in the activation state. Another prominent effect of the reduced side-chain volume at position 107 is to stabilize the ligand-binding intermediate state. Because the stabilization can be ascribed to decrease of the positive enthalpy, the enlarged free space near position 107 would relieve unfavorable steric interactions between the ligand and nearby amino acid residues. The side-chain volume at position 107, therefore, is crucial for the kinetic barrier for the ligand migration and free energy of the ligand-binding intermediate state, which allows us to propose that some photodissociated O(2) moves toward position 107 to be trapped and then expelled to the solvent.

Mapping the pathways for O2 entry into and exit from myoglobin

The effects of mutagenesis on geminate and bimolecular O2 rebinding to 90 mutants at 27 different positions were used to map pathways for ligand movement into and out of sperm whale myoglobin. By analogy to a baseball glove, the protein "catches" and then "holds" incoming ligand molecules long enough to allow bond formation with the iron atom. Opening of the glove occurs by outward movements of the distal histidine (His(64)), and the ligands are trapped in the interior "webbing" of the distal pocket, in the space surrounded by Ile(28), Leu(29), Leu(32), Val(68), and Ile(107). The size of this pocket is a major determinant of the rate of ligand entry into the protein. Immediately after photo- or thermal dissociation, O2 moves away from the iron into this interior pocket. The majority of the dissociated ligands return to the active site and either rebind to the iron atom or escape through the His(64) gate. A fraction of the ligands migrate further away from the heme group into cavities that have been defined as Xe binding sites 4 and 1; however, most of these ligands also return to the distal pocket, and net escape through the interior of wild-type myoglobin is <20-25%.

EPR studies on the photoinduced intermediates of NO complexes in recombinant ferric-Mb trapped at low temperatures

The nitrosyl complex of ferric myoglobin is EPR-silent. Upon photolysis at low temperatures, the photoinduced intermediates trapped in the distal heme cavity exhibit new EPR spectra due to the interaction between the photodissociated NO (S=1/2) and the ferric high spin heme (S=5/2). In order to elucidate the effect of distal E7 (His64) and E11 (Val68) mutations upon the electronic structure of the metal center, its immediate environment, and its interaction with the photodissociated NO, EPR spectra of the photoproducts of the NO complexes of recombinant ferric Mb mutants were measured at 5 K. EPR spectra of the photoproducts were closely related to the size and/or the polarity of the distal pocket residues. The distal pocket of the E7 mutants seemed to be sterically crowded, even decreasing the side chain volume or changing its hydrophobicity by replacing amino acid at position 64. We have found that the mobility of the photodissociated NO molecule in the distal heme pocket was strongly governed by the nature of the amino acid residue at E11 position.

Multiple geminate ligand recombinations in human hemoglobin

The geminate ligand recombination reactions of photolyzed carbonmonoxyhemoglobin were studied in a nanosecond double-excitation-pulse time-resolved absorption experiment. The second laser pulse, delayed by intervals as long as 400 ns after the first, provided a measure of the geminate kinetics by rephotolyzing ligands that have recombined during the delay time. The peak-to-trough magnitude of the Soret band photolysis difference spectrum measured as a function of the delay between excitation pulses showed that the room temperature kinetics of geminate recombination in adult human hemoglobin are best described by two exponential processes, with lifetimes of 36 and 162 ns. The relative amounts of bimolecular recombination to T- and R-state hemoglobins and the temperature dependence of the submicrosecond kinetics between 283 and 323 K are also consistent with biexponential kinetics for geminate recombination. These results are discussed in terms of two models: geminate recombination kinetics modulated by concurrent protein relaxation and heterogeneous kinetics arising from alpha and beta chain differences.

Altered ligand rebinding kinetics due to distal-side effects in hemoglobin chico (Lysbeta66(E10) --> thr)

Hb Chico is an unusual human hemoglobin variant that has lowered oxygen affinity, but unaltered cooperativity and anion sensitivity. Previous studies showed these features to be associated with distal-side heme pocket alterations that confer increased structural rigidity on the molecule and that increase water content in the beta-chain heme pocket. We report here that the extent of nanosecond geminate rebinding of oxygen to the variant and its isolated beta-chains is appreciably decreased. Structural alterations in this variant decrease its oxygen recombination rates without significantly altering rates of migration out of the heme pocket. Data analysis indicates that one or more barriers that impede rebinding of oxygen from docking sites in the heme pocket are increased, with less consequence for CO rebinding. Resonance Raman spectra show no significant alterations in spectral regions sensitive to interactions between the heme iron and the proximal histidine residue, confirming that the functional differences in the variant are due to distal-side heme pocket alterations. These effects are discussed in the context of a schematic representation of heme pocket wells and barriers that could aid the design of novel hemoglobins with altered ligand affinity without loss of the normal allosteric responses that facilitate unloading of oxygen to respiring tissues.

Myoglobin mutants giving the largest geminate yield in CO rebinding in the nanosecond time domain

We have measured the rebinding of carbon monoxide (CO) to some distal mutants of myoglobin (Mb) in the time range from 10(-8) to 10(-1) s by flash photolysis, in which the photodissociated CO rebinds to the heme iron without escaping to the solvent water from the protein matrix. We have found that the double mutants [His64-->Val/Val68-->Thr (H64V/V68T) and His64-->Val/Val68-->Ser (H64V/V68S)] have an extremely large geminate yield (70-80%) in water at 5 degreesC, in contrast to the 7% of the geminate yield of wild-type Mb. The CO geminate yields for these two mutants are the largest in those of Mb mutants reported so far, showing that the two mutants have a unique heme environment that favors CO geminate rebinding. Comparing the crystal structures and 1H-NMR and vibrational spectral data of H64V/V68T and H64V/V68S with those of other mutants, we discuss factors that may control the nanosecond geminate CO rebinding and CO migration in the protein matrix.

Properties of human hemoglobins with increased polarity in the alpha- or beta-heme pocket. Carbonmonoxy derivatives

The spectroscopic, conformational, and functional properties of mutant carbonmonoxy hemoglobins in which either the beta-globin Val67(E11) or the alpha-globin Val62(E11) is replaced by threonine have been investigated. The thermal evolution of the Soret absorption band and the stretching frequency of the bound CO were used to probe the stereodynamic properties of the heme pocket. The functional properties were investigated by kinetic measurements. The spectroscopic and functional data were related to the conformational properties through molecular analysis. The effects of this nonpolar-to-polar isosteric mutation are: (i) increase of heme pocket anharmonic motions, (ii) stabilization of the A0 conformer in the IR spectrum, (iii) increased CO dissociation rates. The spectroscopic data indicate that for the carbonmonoxy derivatives, the Val --> Thr mutation has a larger conformational effect on the beta-subunits than on the alpha-subunits. This is at variance with the deoxy derivatives where the conformational modification was larger in the heme pocket of the alpha-subunit (Cupane, A., Leone, M., Militello, V., Friedman, R. K., Koley, A. P., Vasquez, G. P., Brinigar, W. S., Karavitis, M., and Fronticelli, C. (1997) J. Biol. Chem. 272, 26271-26278). These effects are attributed to a different electrostatic interaction between Ogamma of Thr(E11) and the bound CO molecule. Molecular analysis indicates a more favorable interaction of the bound CO with Thr Ogamma in the beta-subunit heme pocket.

The effects of heme pocket hydrophobicity on the ligand binding dynamics in myoglobin as studied with leucine 29 mutants

To examine the effects of heme pocket hydrophobicity on the ligand binding in myoglobin, some artificial mutants of human myoglobin have been prepared, in which less hydrophobic amino acid residue (Ala, Gly, Ser) is located at the Leu29 (10th residue of the B helix) position. CO rebinding rates for the mutants were markedly decelerated, while the 1H, and 15N NMR spectra of the mutants show that the structural changes around the heme iron for these mutants are rather small. The kinetic and structural properties of the mutants indicate that the ligand binding rate depends on the hydrophobicity inside the heme cavity for these mutants in addition to the volume of the side chain at the 29-position. On the basis of the IR stretching frequency of liganded CO, invasion of water molecules into the heme pocket in the mutants is suggested, which would be induced by the decrease in the hydrophobicity due to the amino acid substitution. A slight red shift of the position of the Soret peak for the serine mutant L29S also supports the reduced hydrophobicity inside the heme cavity. We can concluded that, together with the kinetic properties of the mutants, the hydrophobicity of the heme pocket is one of the key factors in regulating the ligand binding to the heme iron.

Modification of alpha-chain or beta-chain heme pocket polarity by Val(E11) --> thr substitution has different effects on the steric, dynamic, and functional properties of human recombinant hemoglobin. Deoxy derivatives

The dynamic and functional properties of mutant deoxyhemoglobins in which either the beta-globin Val67(E11) or the alpha-globin Val62(E11) is replaced by threonine have been investigated through the thermal evolution of the Soret absorption band in the temperature range 300 to 20 K and through the kinetics of CO rebinding after flash photolysis at room temperature. The conformational properties of the modified alpha chain and beta chain distal heme pockets were also studied through x-ray crystallography and molecular modeling. The data obtained with the various techniques consistently indicate that the polar isosteric mutation in the distal side of the alpha chain heme pocket has a larger effect on the investigated properties than the analogous mutation on the beta chain. We attribute the observed differences to the presence of a water molecule in the distal heme pocket of the modified alpha chains, interacting with the hydroxyl of the threonine side chain. This is indicated by molecular modeling which showed that the water molecule present in the alpha chain distal heme pocket can bridge by H bonding between Thr62(E11) and His58(E7) without introducing any unfavorable steric interactions. Consistent with the dynamic and functional data, the presence of a water molecule in the distal heme pocket of the modified beta chains is not observed by x-ray crystallography.

Reaction of variant sperm-whale myoglobins with hydrogen peroxide: the effects of mutating a histidine residue in the haem distal pocket

The reaction of hydrogen peroxide with a number of variants of sperm-whale myoglobin in which the distal pocket histidine residue (His64) had been mutated was studied with a combination of stopped-flow spectroscopy and freeze-quench EPR. The rate of the initial bimolecular reaction with hydrogen peroxide in all the proteins studied was found to depend on the polarity of the amino acid side chain at position 64. In wild-type myoglobin there were no significant optical changes subsequent to this reaction, suggesting the rapid formation of the well-characterized oxyferryl species. This conclusion was supported by freeze-quench EPR data, which were consistent with the pattern of reactivity previously reported [King and Winfield (1963) J. Biol. Chem. 238, 1520-1528]. In those myoglobins bearing a mutation at position 64, the initial bimolecular reaction with hydrogen peroxide yielded an intermediate species that subsequently decayed via a second hydrogen peroxide-dependent step leading to modification or destruction of the haem. In the mutant His64-->Gln the calculated electronic absorption spectrum of the intermediate was not that of an oxyferryl species but seemed to be that of a low-spin ferric haem. Freeze-quench EPR studies of this mutant and the apolar mutant (His64-->Val) revealed the accumulation of a novel intermediate after the first hydrogen peroxide-dependent reaction. The unusual EPR characteristics of this species are provisionally assigned to a low-spin ferric haem with bound peroxide as the distal ligand. These results are interpreted in terms of a reaction scheme in which the polarity of the distal pocket governs the rate of binding of hydrogen peroxide to the haem iron and the residue at position 64 governs both the rate of heterolytic oxygen scission and the stability of the oxyferryl product.

A chemometric analysis of the resonance Raman spectra of mutant carbonmonoxy-myoglobins reveals the effects of polarity

Resonance Raman spectra of 10 carbonmonoxy-myoglobins have been obtained, including sperm whale native, pig wild-type, and the mutants H64L, H64A, V68T, V68N, H64V/V68T, F43W, F46V, and L29F. This series was chosen in order to study the effect of ligand binding pocket polarity on the positions of the v(Fe-CO) and delta (Fe-C-O) bands. Spectra of both 12CO and 13CO isotopic forms have been obtained and a detailed analysis has facilitated the identification of both the ligand-specific bands and six underlying porphyrin bands which are insensitive to this isotopic substitution. Along with a band-fitting analysis of infrared spectra, these resonance Raman data provide a comprehensive evaluation of the vibrations of the FeCO unit. The band positions of the ligand-specific modes are found to depend on the structure of the ligand binding pocket, arising from the strength of back-bonding within the FeCO unit, and clear correlations exist between the v(Fe-CO), delta (Fe-C-O), and v(C-O) band positions which characterize this synergic bonding. The results are consistent with the proposal that the vibration band positions are determined primarily by the electrostatic potential at the ligand. Five discrete band sets are observed for this set of mutants, suggesting that 5 discrete conformations occur.

Protein engineering strategies for designing more stable hemoglobin-based blood substitutes

Over the past five years our laboratory has been using rational, comparative, and random combinatorial mutagenesis strategies to optimize the alpha and beta subunits of recombinant human hemoglobin (Hb) for efficient O2 transport, greater stability, and minimum interference with vascular activity. In each approach, mammalian myoglobin (Mb) has been used as a prototype to develop experimental methodologies and to study the stereochemical mechanisms that govern O2 affinity, discrimination against CO, rates of ligand binding, auto- and chemically induced oxidation, resistance to hemin loss, and stability to globin denaturation. Multiple replacements in the distal portion of the heme pocket have been designed rationally to lower oxygen affinity and at the same time inhibit oxidative side reactions. The P50 values are adjusted by altering electrostatic and steric interactions between the bound ligand and residues at the Leu(B10), His(E7), and Va(E11) positions. Large apolar residues (Leu, Phe, Trp) at the B10 and E11 positions inhibit NO-induced and autooxidation in both myoglobin and hemoglobin by excluding oxidants and proton donors from the immediate vicinity of the bound ligand. Similar strategies appear to have evolved in a number of animal myoglobins and hemoglobins which have unusual amino acids at the E7, B10, and E11 positions. Random combinatorial mutagenesis techniques have been developed to insert new amino acid combinations near the bound ligand in sperm whale Mb. The objective is to obtain "unnatural" distal pocket structures that enhance O2 transport and resistance to oxidation by alternative mechanisms.

Binding of nitric oxide and carbon monoxide to soluble guanylate cyclase as observed with Resonance raman spectroscopy

Resonance Raman spectra have been recorded for the ferrous heme of soluble guanylate cyclase (sGC), the only receptor known thus far for .NO. On the basis of the frequencies of porphyrin core sensitive vibrations in the high frequency region of the Raman spectrum, we conclude that the ferrous heme is five-coordinate, high spin, when no exogenous ligands are present. We assign a prominent vibration that occurs at 204 cm-1 in the reduced enzyme to the heme Fe(2+)-proximal histidine stretching vibration. In the .NO bound form of the enzyme, the heme Fe2+ retains a five-coordinate geometry. Assuming that .NO binds to the distal side of the heme, this observation indicates that the weak Fe-His bond breaks when .NO binds. Two isotope-sensitive vibrations are observed in the .NO bound enzyme, one at 1677 cm-1, attributed to the N-O stretching vibration, and one at 525 cm-1, attributed to the Fe-NO stretching vibration. When CO is bound to the ferrous heme, the heme ligation is six-coordinate. From this, we conclude that the Fe-His bond is still intact and that, if cleavage of the Fe-proximal ligand bond is necessary for complete activation of sGC, then CO should only weakly activate the enzyme, which has been shown to be the case. In the carbonmonoxy enzyme, the Fe-CO stretching vibration is observed at 472 cm-1 and the Fe-C-O bending vibration is detected at 562 cm-1. These frequencies are the lowest yet observed for the Fe-CO stretching and Fe-C-O bending modes in heme proteins or model systems with imidazole as the proximal ligand and suggest that there is significant negative polarity in the distal pocket. The negative polarity and the low frequency of the Fe-His stretching vibration may account for the very low O2 affinity of sGC.

A double mutant of sperm whale myoglobin mimics the structure and function of elephant myoglobin

The functional, spectral, and structural properties of elephant myoglobin and the L29F/H64Q mutant of sperm whale myoglobin have been compared in detail by conventional kinetic techniques, infrared and resonance Raman spectroscopy, 1H NMR, and x-ray crystallography. There is a striking correspondence between the properties of the naturally occurring elephant protein and those of the sperm whale double mutant, both of which are quite distinct from those of native sperm whale myoglobin and the single H64Q mutant. These results and the recent crystal structure determination by Bisig et al. (Bisig, D. A., Di Iorio, E. E., Diederichs, K., Winterhalter, K. H., and Piontek, K. (1995) J. Biol. Chem. 270, 20754-20762) confirm that a Phe residue is present at position 29 (B10) in elephant myoglobin, and not a Leu residue as is reported in the published amino acid sequence. The single Gln64(E7) substitution lowers oxygen affinity approximately 5-fold and increases the rate of autooxidation 3-fold. These unfavorable effects are reversed by the Phe29(B10) replacement in both elephant myoglobin and the sperm whale double mutant. The latter, genetically engineered protein was originally constructed to be a blood substitute prototype with moderately low O2 affinity, large rate constants, and increased resistance to autooxidation. Thus, the same distal pocket combination that we designed rationally on the basis of proposed mechanisms for ligand binding and autooxidation is also found in nature.

Phe-46(CD4) orients the distal histidine for hydrogen bonding to bound ligands in sperm whale myoglobin

The role of Phe-46(CD4) in modulating the functional properties of sperm whale myoglobin was investigated by replacing this residue with Leu, Ile, Val, Ala, Trp, Tyr, and Glu. This highly conserved amino acid almost makes direct contact with the distal histidine and has been postulated to affect ligand binding. The overall association rate constants for CO, O2, and NO binding were little affected by decreasing the size of residue 46 step-wise from Phe to Leu to Val to Ala. In contrast, the rates of CO, O2, and NO dissociation increased 4-, 10-, and 25-fold, respectively, for the same series of mutants, causing large decreases in the affinity of myoglobin for all three diatomic gases. The rates of autooxidation at 37 degrees C, pH 7.0 increased dramatically from approximately 0.1-0.3 h-1 for wild-type, Tyr-46, and Trp-46 myoglobins to 1.5, 5.2, 4.9, and 5.0 h-1 for the Leu-46, Ile-46, Val-46 and Ala-46 mutants, respectively. Rates of NO and O2 geminate recombination were measured using 35 ps and 9 ns laser excitation pulses. Decreasing the size of residue 46 causes significant decreases in the extent of both picosecond and nanosecond rebinding processes. High resolution structures of Leu-46 and Val-46 metmyoglobins, Val-46 CO-myoglobin, and Val-46 deoxymyoglobin were determined by X-ray crystallography. When Phe-46 is replaced by Val, the loss of internal packing volume is compensated by (1) contraction of the CD corner toward the core of the protein, (2) movement of the E-helix toward the mutation site, (3) greater exposure of the distal pocket to intruding solvent molecules, and (4) large disorder in the position of the side chain of the distal histidine (His-64). In wild-type myoglobin, the van der Waals contact between C zeta of Phe-46 and C beta of His-64 appears to restrict rotation of the imidazole side chain. Insertion of Val at position 46 relieves this steric restriction, allowing the imidazole side chain to rotate about the C alpha - C beta bond toward the surface of the globin and about the C beta - C gamma bond toward the space previously occupied by the native Phe-46 side chain. This movement disrupts hydrogen bonding with bound ligands, causing significant decreases in affinity, and opens the distal pocket to solvent water molecules, causing marked increases in the rate of autooxidation.(ABSTRACT TRUNCATED AT 400 WORDS)

Alteration of axial coordination by protein engineering in myoglobin. Bisimidazole ligation in the His64-->Val/Val68-->His double mutant

Pig and human myoglobin have been engineered to reverse the positions of the distal histidine and valine (i.e. His64(E7)-->Val and Val68(E11)-->His). Spectroscopic and ligand binding properties have been measured for human and pig H64V/V68H myoglobin, and the structure of the pig H64V/V68H double mutant has been determined to 2.07-A resolution by x-ray crystallography. The crystal structure shows that the N epsilon of His68 is located 2.3 A away from the heme iron, resulting in the formation of a hexacoordinate species. The imidazole plane of His68 is tilted relative to the heme normal; moreover it is not parallel to that of His93, in agreement with our previous proposal (Qin, J., La Mar, G. N., Dou, Y., Admiraal, S. J., and Ikeda-Saito, M. (1994) J. Biol. Chem. 269, 1083-1090). At cryogenic temperatures, the heme iron is in a low spin state, which exhibits a highly anisotropic EPR spectrum (g1 = 3.34, g2 = 2.0, and g3 < 1), quite different from that of the imidazole complex of metmyoglobin. The mean iron-nitrogen distance is 2.01 A for the low spin ferric state as determined by x-ray spectroscopy. The ferrous form of H64V/V68H myoglobin shows an optical spectrum that is similar to that of b-type cytochromes and consistent with the hexacoordinate bisimidazole hemin structure determined by the x-ray crystallography. The double mutation lowers the ferric/ferrous couple midpoint potential from +54 mV of the wild-type protein to -128 mV. Ferrous H64V/V68H myoglobin binds CO and NO to form stable complexes, but its reaction with O2 results in a rapid autooxidation to the ferric species. All of these results demonstrate that the three-dimensional positions of His64 and Val68 in the wild-type myoglobin are as important as the chemical nature of the side chains in facilitating reversible O2 binding and inhibiting autooxidation.

Stability of myoglobin: a model for the folding of heme proteins

Factors governing the stability of sperm whale, pig, and human metmyoglobin were examined by (1) measuring guanidinium chloride induced unfolding of apoglobins containing 22 replacements at positions 29(B10), 43(CD1), 64(E7), 68(E11), and 107(G8), (2) determining the rates of hemin loss from the recombinant holoproteins, and (3) estimating constitutive expression levels of the corresponding genes in Escherichia coli TB-1 cells. The denaturant titrations were analyzed in terms of a two-step unfolding reaction, N(native apoprotein)-->I(intermediate)-->U(unfolded), in which the intermediate is visualized by an increase in tryptophan fluorescence emission. Two key conclusions were reached. First, high rates of hemin loss are not necessarily correlated with unstable globin structures and vice versa. In general, both rates of hemin loss and the equilibrium constants for apoprotein unfolding must be determined in order to understand the overall stability of heme proteins and to predict the efficiency of their expression. Second, polar residues in the distal pocket cause marked decreases in the overall stability of apomyoglobin. Removal of hemin from V68N and L29N sperm whale myoglobins produces the molten globular I state at pH 7, 25 degrees C, without addition of denaturant. In contrast, the H64L and H64F mutations produce apoproteins which are 10-30 times more stable than wild-type apoglobin. The latter results show that protein stability is sacrificed in order to have the distal histidine (H64) present to increase O2 affinity and inhibit autooxidation.

Engineering Ascaris hemoglobin oxygen affinity in sperm whale myoglobin: role of tyrosine B10

The contribution to oxygen stabilization of a tyrosine residue in topological position (B10) has been studied in sperm whale myoglobin by simultaneous replacement of residues at positions (B10), (E7) and (E10) as suggested by analysis of the sequence of high oxygen affinity hemoglobins, such as that of the nematode Ascaris suum. Kinetic and equilibrium experiments with the gaseous ligands oxygen and carbon monoxide show that indeed the introduction of tyrosine (B10), together with replacement of the distal histidine (E7) with glutamine, is associated with a large decrease in the oxygen dissociation rate constant. Our results are consistent with the possible formation in the distal pocket of two hydrogen bonds with the iron-bound oxygen.

Correlation between the steric bulk of the distal E7 and E11 residues and the tilt of the FeCN unit in cyanometmyoglobin as determined by NMR from the orientation of the magnetic axes in single and double point mutants

The amino acids in the heme pocket of sperm whale myoglobin single E11 and double E7 and E11 point mutants in the metcyano form have been assigned by NMR methods to assess the role of steric bulk in modulating ligand tilt. The five mutants investigated are the single mutants His64(E7)-->Gly (H[E7]G), Val68(E11)-->Ile (V[E11]I), and Val68(E11)-->Ala (V[E11]A) and the double mutants His64-(E7)-->Gly:Val68(E11)-->Ile (H,V[E7,E11]G,I) and His64(E7)-->Gly:Val68(E11)-->Ala (H,V[E7,E11]G,A). The dipolar (NOESY) contacts on the proximal side of the heme confirm a conserved molecular structure for all of the mutants. The proximal residue coordinates, together with the dipolar shifts for proximal side residues, quantitatively yield the orientations of the magnetic susceptibility tensors, whose major axis corresponds to the orientation of the ligand. It is observed that upon reduction of the steric bulk in the V[E11]A mutant, the tilt of the ligand is significantly reduced (approximately 8 degrees) from that in the wild type (WT) (approximately 16 degrees), with little change in the direction of tilt. In the case of increased steric bulk at position 68 in the V[E11]I mutant, it is observed that the extent and direction of the tilt are essentially the same as in WT, and it is shown that this is due to the fact that Ile68 is oriented in the pocket with its C delta H3 directed away from the iron.(ABSTRACT TRUNCATED AT 250 WORDS)

Discovery of new ligand binding pathways in myoglobin by random mutagenesis

A random library of single amino acid mutants of myoglobin was generated using a highly efficient, single-base-misincorporation random mutagenesis method to discover new ligand-binding pathways in myoglobin. A surprisingly large fraction of the library exhibits ligand-binding kinetics that are substantially different from the wild-type protein. In addition to residues 45, 64 and 68, which comprise the best studied ligand-binding pathway single mutations of several other clusters of residues far away from that pathway are discovered which profoundly affect the ligand-binding kinetics. These results provide a new approach to explore the relationship between the fluctuations in protein structure and function.