Perturbations of the distal heme pocket in human myoglobin mutants probed by infrared spectroscopy of bound CO: correlation with ligand binding kinetics

A new inhibition mechanism in the multifunctional catalytic hemoglobin dehaloperoxidase as revealed by the DHP A(V59W) mutant: A spectroscopic and crystallographic study

As multifunctional catalytic hemoglobins, dehaloperoxidase isoenzymes A and B (DHP A and B) are among the most versatile hemoproteins in terms of activities displayed. The ability of DHP to bind over twenty different substrates in the distal pocket might appear to resemble the promiscuousness of monooxygenase enzymes, yet there are identifiable substrate-specific interactions that can steer the type of oxidation (O-atom vs. electron transfer) that occurs inside the DHP distal pocket. Here, we have investigated the DHP A(V59W) mutant in order to probe the limits of conformational flexibility in the distal pocket as it relates to the genesis of this substrate-dependent activity differentiation. The X-ray crystal structure of the metaquo DHP A(V59W) mutant (PDB 3K3U) and the V59W mutant in complex with fluoride [denoted as DHP A(V59W-F)] (PDB 7MNH) show significant mobility of the tryptophan in the distal pocket, with two parallel conformations having W59-N[Formula: see text] H-bonded to a heme-bound ligand (H2O or F[Formula: see text], and another conformation [observed only in DHP A(V59W-F)] that brings W59 sufficiently close to the heme as to preclude axial ligand binding. UV-vis and resonance Raman spectroscopic studies show that DHP A(V59W) is 5-coordinate high spin (5cHS) at pH 5 and 6-coordinate high spin (6cHS) at pH 7, whereas DHP A(V59W-F) is 6cHS from pH 5 to 7. Enzyme assays confirm robust peroxidase activity at pH 5, but complete loss of activity at pH 7. We find no evidence that tryptophan plays a role in the oxidation mechanism ([Formula: see text]. radical formation). Instead, the data reveal a new mechanism of DHP inhibition, namely a shift towards a non-reactive form by OH[Formula: see text] ligation to the heme-Fe that is strongly stabilized (presumably through H-bonding interactions) by the presence of W59 in the distal cavity.

Quantification of the Electrostatic Effect on Redox Potential by Positive Charges in a Catalyst Microenvironment

Charged functional groups in the secondary coordination sphere (SCS) of a heterogeneous nanoparticle or homogeneous electrocatalyst are of growing interest due to enhancements in reactivity that derive from specific interactions that stabilize substrate binding or charged intermediates. At the same time, accurate benchmarking of electrocatalyst systems most often depends on the development of linear free-energy scaling relationships. However, the thermodynamic axis in those kinetic-thermodynamic correlations is most often obtained by a direct electrochemical measurement of the catalyst redox potential and might be influenced by electrostatic effects of a charged SCS. In this report, we systematically probe positive charges in a SCS and their electrostatic contributions to the electrocatalyst redox potential. A series of 11 iron carbonyl clusters modified with charged and uncharged ligands was probed, and a linear correlation between the ν_CO absorption band energy and electrochemical redox potentials is observed except where the SCS is positively charged.

Effect of Occluded Ligand Migration on the Kinetics and Structural Dynamics of Homodimeric Hemoglobin

Small molecules such as molecular oxygen, nitric oxide, and carbon monoxide play important roles in life, and many proteins require the transport of small molecules to and from the bulk solvent for their function. Ligand migration within a protein molecule is expected to be closely related to the overall structural changes of the protein, but the detailed and quantitative connection remains elusive. For example, despite numerous studies, how occluded ligand migration affects the kinetics and structural dynamics of the R-T transition remains unclear. To shed light on this issue, we chose homodimeric hemoglobin (HbI) with the I114F mutation (I114F), which is known to interfere with ligand migration between the primary and secondary docking sites, and studied its kinetics and structural dynamics using time-resolved X-ray solution scattering. The kinetic analysis shows that I114F has three structurally distinct intermediates (I₁, I₂, and I₃) as in the wild type (WT), but its geminate CO recombination occurs directly from I₁ without the path via I₂ observed in WT. Moreover, the structural transitions, which involve ligand migration (the transitions from I₁ to I₂ and from I₃ to the initial state), are decelerated compared to WT. The structural analysis revealed that I114F involves generally smaller structural changes in all three intermediates compared to WT.

The reductive phase of Rhodobacter sphaeroides cytochrome c oxidase disentangled by CO ligation

Cytochrome c oxidase (CcO) is a membrane protein of the respiratory chain that catalytically reduces molecular oxygen (O2) to water while translocating protons across the membrane. The enzyme hosts two copper and two heme iron moieties (heme a/heme a3). The atomic details of the sequential steps that go along with this redox-driven proton translocation are a matter of debate. Particularly for the reductive phase of CcO that precedes oxygen binding experimental data are scarce. Here, we use CcO under anaerobic conditions where carbon monoxide (CO) is bound to heme a3 which in tandem with CuB forms the binuclear center (BNC). Fourier-transform infrared (FTIR) absorption spectroscopy is combined with electro-chemistry to probe different redox and protonation states populated by variation of the external electrostatic potential. With this approach, the redox behavior of heme a and the BNC could be separated and the corresponding redox potentials were determined. We also infer the protonation of one of the propionate side chains of heme a3 to correlate with the oxidation of heme a. Experimental changes in the local electric field surrounding CO bound to heme a3 are determined by their vibrational Stark effect and agree well with electrostatic computations. The comparison of experimental and computational results indicates that changes of the heme a3/CuB redox state are coupled to proton transfer towards heme a3. The latter supports the role of the heme a3 propionate D as proton loading site.

How does hemoglobin generate such diverse functionality of physiological relevance?

The absolute values of the O2-affinities (P50, Klow, and Khigh) of hemoglobin (Hb) are regulated neither by changes in the static T-/R-quaternary and associated tertiary structures nor the ligation states. They are pre-determined and regulated by the extrinsic environmental factors such as pH, buffers, and heterotropic effectors. The effect and role of O2 on Hb are reversibly to drive the structural allosteric equilibrium between the T(deoxy)- and R(oxy)-Hb toward R(oxy)-Hb (the structural allostery). R(oxy)-Hb has a higher O2-affinity (Khigh) relative to that (Klow) of the T(deoxy)-Hb (Khigh>Klow) under any fixed environmental conditions. The apparent O2-affinity of Hb is high, as the globin matrix interferes with the dissociation process of O2, forcing the dissociated O2 geminately to re-bind to the heme Fe. This artificially increases [oxy-Hb] and concomitantly decreases [deoxy-Hb], leading to the apparent increases of the O2-affinity of Hb. The effector-linked high-frequency thermal fluctuations of the globin matrix act as a gating mechanism to modulate such physical, energetic, and kinetic barriers to enhance the dissociation process of O2, resulted in increases in [deoxy-Hb] and concomitant decrease in [oxy-Hb], leading to apparent reductions of the O2-affinity of Hb (the entropic allostery). The heme in Hb is simply a low-affinity O2-trap, the coordination structure of which is not altered by static T-/R-quaternary and associated tertiary structural changes of Hb. Thus, heterotrophic effectors are the signal molecule, which acts as a functional link between these two allosteries and generates the diverse functionality of Hb of physiological relevance. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Photolysis of Hi-CO Nitrogenase - Observation of a Plethora of Distinct CO Species using Infrared Spectroscopy

Fourier transform infrared spectroscopy (FT-IR) was used to study the photochemistry of CO-inhibited Azotobacter vinelandii nitrogenase using visible light at cryogenic temperatures. The FT-IR difference spectrum of photolyzed hi-CO at 4 K comprises negative bands at 1973 cm^-1 and 1679 cm^-1 together with positive bands at 1711 cm^-1, 2135 and 2123 cm^-1. The negative bands are assigned to a hi-CO state that comprises 2 metal-bound CO ligands, one terminally bound, and one bridged and/or protonated species. The positive band at 1711 cm^-1 is assigned to a lo-CO product with a single bridged and/or protonated metal-CO group. We term these species 'Hi-1' and 'Lo-1' respectively. The high-energy bands are assigned to a liberated CO trapped in the protein pocket. Warming results in CO recombination, and the temperature dependence of the recombination rate yields an activation energy of 4 kJ mol^-1. Two α-H195 variant enzymes yielded additional signals. Asparagine substitution, α-H195N, gives a spectrum containing 2 negative 'Hi-2' bands at 1936 and 1858 cm^-1 with a positive 'Lo-2' band at 1780 cm^-1, while glutamine substitution, α-H195Q, produces a complex spectrum that includes a third CO species, with negative 'Hi-3' bands at 1938 and 1911 cm^-1 and a positive feature 'Lo-3' band at 1921 cm^-1. These species can be assigned to a combination of terminal, bridged, and possibly protonated CO groups bound to the FeMo-cofactor active site. The proposed structures are discussed in terms of both CO inhibition and the mechanism nitrogenase catalysis. Given the intractability of observing nitrogenase intermediates by crystallographic methods, IR-monitored photolysis appears to be a promising and information-rich probe of nitrogenase structure and chemistry.

Engineering of the heme pocket of an H-NOX domain for direct cyanide detection and quantification

A new cyanide sensing system, the Heme-Nitric oxide and/or OXygen binding domain (H-NOX domain) from Thermoanaerobacter tengcongensis (Tt H-NOX), has been investigated. With straightforward absorbance-based detection, we have achieved a cyanide detection limit of 0.5 microM (approximately 10 ppb) with an upper detection range that is adjustable with protein concentration. We find a linear correlation of multiple spectroscopic features with cyanide concentration. These spectroscopic features include the Soret band maximum and absorbance changes in both the Soret and alpha/beta band regions of the spectrum. Multiple probes for cyanide detection makes sensing with Tt H-NOX unique compared to other cyanide sensing methods. Furthermore, using site-directed mutagenesis, we have rationally engineered the heme pocket of Tt H-NOX to improve its cyanide sensing properties. Using a mutant that alters the heme structure of Tt H-NOX (P115A) we were able to introduce colorimetric detection of cyanide. Through substituting phenylalanine 78 with a smaller (valine, F78V) or a larger residue (tyrosine, F78Y), we demonstrate a correlation with distal pocket steric crowding and affinity for cyanide. In particular, F78V Tt H-NOX shows a significant increase in CN(-) binding affinity and selectivity. Thus, we demonstrate the ability to fine-tune the affinity and specificity of Tt H-NOX for cyanide, suggesting that Tt H-NOX can be readily tailored into a practical cyanide sensor.

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.

Probing the role of hydration in the unfolding transitions of carbonmonoxy myoglobin and apomyoglobin

We show that the equilibrium unfolding transition of horse carbonmonoxy myoglobin monitored by the stretching vibration of the CO ligand, a local environmental probe, is very sharp and, thus, quite different from those measured by global conformational reporters. In addition, the denatured protein exhibits an A(0)-like CO band. We hypothesize that this sharp transition reports penetration of water into the heme pocket of the protein. Parallel experiments on horse apomyoglobin, wherein an environment-sensitive fluorescent probe, nile red, was used, also reveals a similar putative hydration event. Given the importance of dehydration in protein folding and also the recent debate over the interpretation of probe-dependent unfolding transitions, these results have strong implications on the mechanism of protein folding.

DFT analysis of axial and equatorial effects on heme-CO vibrational modes: applications to CooA and H-NOX heme sensor proteins

Determinants of the Fe-CO and C-O stretching frequencies in (imidazole)heme-CO adducts have been investigated via density functional theory (DFT) analysis, in connection with puzzling characteristics of the heme sensor protein CooA and of the H-NOX (Heme-Nitric Oxide and/or OXygen binding) family of proteins, including soluble guanylate cyclase (sGC). The computations show that two mechanisms of Fe-histidine bond weakening have opposite effects on the nuFeC/nuCO pattern. Mechanical tension is expected to raise nuFeC with little change in nuCO whereas the weakening of H-bond donation from the imidazole ligand has the opposite effect. Data on CooA indicate imidazole H-bond weakening associated with heme displacement, as part of the activation mechanism. The computations also reveal that protein-induced distortion of the porphyrin ring, a prominent structural feature of the H-NOX protein TtTar4H (Thermoanaerobacter tengcongensis Tar4 protein heme domain), has surprisingly little effect on nuFeC or nuCO. However, another structural feature, strong H-bonding to the propionates, is suggested to account for the weakened back bonding that is evident in sGC. TtTar4H-CO itself has an elevated nuFeC, which is successfully modeled as a compression effect, resulting from steric crowding in the distal pocket. nuFeC/nuCO data, in conjunction with modeling, can provide valuable insight into mechanisms for heme-protein modulation.

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of the dimeric hemoglobin of Scapharca inaequivalvis (HbI) after photodissociation. By combining these studies with X-ray crystallography, three transient ligand docking sites were identified: a primary docking site B in close vicinity to the heme iron, and two secondary docking sites C and D corresponding to the Xe4 and Xe2 cavities of myoglobin. To assess the relevance of these findings for physiological binding, we also performed flash photolysis experiments on HbICO at room temperature and equilibrium binding studies with dioxygen. Our results show that the Xe4 and Xe2 cavities serve as transient docking sites for unbound ligands in the protein, but not as way stations on the entry/exit pathway. For HbI, the so-called histidine gate mechanism proposed for other globins appears as a plausible entry/exit route as well.

A study of the horseradish peroxidase catalytic site by FTIR spectroscopy

Vibrational changes in the catalytic site of horseradish peroxidase were investigated by FTIR (Fourier-transform infrared) spectroscopy in the 1000-2500 cm(-1) range. Difference spectra were generated by photolysis of the haemII-CO compound at different pH/pD values. The spectra report on the fine structure around the catalytic site and show vibrational changes of protein backbone, amino acid residues and cofactors. Assignments of the FTIR vibrations can be made based upon known crystal structures, comparisons with absorption frequencies and extinction coefficients of model amino acids and cofactors, effects of H2O/2H2O exchange and changes of pH/pD. Concomitant with the photolysis of the CO ligand are changes due to haem and protein vibrations, predominant among which are arginine and histidine residue vibrations.

Ultrafast time-resolved IR studies of protein-ligand interactions

Time-resolved mid-IR spectroscopy combines molecular sensitivity with ultrafast capability to incisively probe protein-ligand interactions in model heme proteins. Highly conserved residues near the heme binding site fashion a ligand-docking site that mediates the transport of ligands to and from the binding site. We employ polarization anisotropy measurements to probe the orientation and orientational distribution of CO when bound to and docked near the active binding site, as well as the dynamics of ligand trapping in the primary docking site. In addition, we use more conventional transient absorption methods to probe the dynamics of ligand escape from this site, as well as the ultrafast dynamics of NO geminate recombination with the active binding site. The systems investigated include myoglobin, hemoglobin, and microperoxidase.

Conformational substates and dynamic properties of carbonmonoxy hemoglobin

Heme pocket dynamics of human carbonmonoxy hemoglobin (HbCO) is studied by Fourier transform infrared spectroscopy. The CO stretching band at various temperatures in the interval 300-10 K is analyzed in terms of three taxonomic A substates; however, in HbCO the band attributed to the A(1) taxonomic substate accounts for approximately 90% of the total intensity in the pH range 8.8-4.5. Two different regimes as a function of temperature are observed: below 160 K, the peak frequency and the bandwidth of the A(1) band have constant values whereas, above this temperature, a linear temperature dependence is observed, suggesting the occurrence of transitions between statistical substates within the A(1) taxonomic substate in this protein. The relationship between the heme pocket dynamics (as monitored by the thermal behavior of the CO stretching band), the overall dynamic properties of the protein matrix (as monitored by the thermal behavior of Amide II and Amide I' bands) and the glass transition of the solvent (as monitored by the thermal behavior of the bending band of water) is also investigated. From this analysis, we derive the picture of a very soft heme pocket of hemoglobin characterized by rather large anharmonic terms and strongly coupled to the dynamic properties of the solvent.

An electrostatic model for the frequency shifts in the carbonmonoxy stretching band of myoglobin: correlation of hydrogen bonding and the stark tuning rate

The effect of internal and applied external electric fields on the vibrational stretching frequency for bound CO (nu(CO)) in myoglobin mutants was studied using density functional theory. Geometry optimization and frequency calculations were carried out for an imidazole-iron-porphine-carbonmonoxy adduct with various small molecule hydrogen-bonding groups. Over 70 vibrational frequency calculations of different model geometries and hydrogen-bonding groups were compared to derive overall trends in the C-O stretching frequency (nu(CO)) in terms of the C-O bond length and Mulliken charge. Simple linear functions were derived to predict the Stark tuning rate using an approach analogous to the vibronic theory of activation.(1) Potential energy calculations show that the strongest interaction occurs for C-H or N-H hydrogen bonding nearly perpendicular to the Fe-C-O bond axis. The calculated frequencies are compared to the structural data available from 18 myoglobin crystal structures, supporting the hypothesis that the vast majority of hydrogen-bonding interactions with CO occur from the side, rather than the end, of the bound CO ligand. The nu(CO) frequency shifts agree well with experimental frequency shifts for multiple bands, known as A states, and site-directed mutations in the distal pocket of myoglobin. The model calculations quantitatively explain electrostatic effects in terms of specific hydrogen-bonding interactions with bound CO in heme proteins.

Formation and characterization of carbon monoxide adducts of iron "twin coronet" porphyrins. Extremely low CO affinity and a strong negative polar effect on bound CO

The carbon monoxide (CO) adducts of iron "twin coronet" porphyrins (TCPs) are characterized by UV-vis, resonance Raman (RR), IR, and 13C NMR spectroscopies. A superstructured porphyrin, designated as TCP, was used as a common framework for the four different types of iron complexes. TCP bears two binaphthalene bridges on each side and creates two hydrophobic pockets surrounded by the bulky aromatic rings. In the CO-binding cavities, the hydroxyl groups are oriented toward the center above the heme. The iron complexes investigated are as follows: TCP (which is without a covalently linked axial ligand), TCP-PY (which has a linked pyridine ligand), and TCP-TB and TCP-TG (both of which have a linked thiolate ligand). These complexes were synthesized as ferric forms and identified by the various spectroscopic methods. The UV-vis spectra of TCP-CO and TCP-PY-CO exhibit lambda(max) at 432, 546 and 428, 541 nm, respectively. On the other hand, the CO adducts of TCP-TB and TCP-TG show typical hyperporphyrin spectra for a thiolate-ligated iron(II) porphyrin-CO complex. In the RR spectra, the nu(Fe-CO) bands were observed at 506, 489 cm(-1) (TCP), 465 cm(-1) (TCP-PY), 458, 437 cm(-1) (TCP-TG) and 429 cm(-1) (TCP-TB). Compared with the reported nu(Fe-CO) frequencies of hemoproteins and their model systems, these observed values are unusually low. Further, abnormally high nu(C-O) bands are observed at 1990 cm(-1) (TCP-CO) and 2008 cm(-1) (TCP-PY-CO) in IR spectra. The lower nu(Fe-CO) and the higher nu(C-O) frequencies can be ascribed to the strong negative polar effect caused by the vicinal hydroxyl groups in the cavity. This prediction is further supported by the observation of significant 13C shieldings exhibited by TCP-CO (delta = 202.6 ppm) and TCP-PY-CO (delta = 202.3 ppm), in comparison to hemoproteins and other heme models. The CO affinity of TCP-PY (P1/2CO = 0.017 Torr at 25 C) is unusually lower than other heme models. The unique behavior of these CO adducts is discussed in context of the TCP structures.

Effect of a charge relay on the vibrational frequencies of carbonmonoxy iron porphine adducts: the coupling of changes in axial ligand bond strength and porphine core size

The effect of a charge relay involving Asp-His-Fe in peroxidase enzymes is explored using density functional theory (DFT) calculations of vibrational spectra and potential energy surfaces of carbonmonoxy model systems. The series of models consists of a carbonmonoxy iron porphine molecule with a trans imidazole ligand hydrogen-bonded to six different partners at the Ndelta position. Calculations on the oxy system and on models of the Asp-His-Ser catalytic triad of serine proteases were also performed to obtain an understanding of how the redistribution of charge in these systems may contribute to enzymatic function. The goal of the study is to relate the experimental frequencies in resonance Raman and Fourier transform infrared studies to bonding that is important for the function of heme enzymes. Calculations of both axial and in-plane modes exhibit trends that agree with experimental data. Comparisons of the charge distribution on the different models show that polarization of iron carbonomonoxy bonds consistent with the mechanism for peroxidase function leads to a frequency reduction in the C-O stretching mode nuCO. The combination of axial trans sigma-bonding and pi-bonding effects that include expansion of the porphine core result in little change in the Fe-C stretching frequency nuFe-CO in the series of molecules studied with different Ndelta-H hydrogen bonding. A particular role for the core size is discussed that demonstrates the applicability of trends observed in vibrational spectroscopy of hemes to the charge relay mechanism and other axial ligation effects. The bonding interactions described account for the increase in electron density on bound diatomic ligands, which is required for peroxidase function.

Connection between the taxonomic substates and protonation of histidines 64 and 97 in carbonmonoxy myoglobin

Infrared spectra of heme-bound CO in sperm whale carbonmonoxy myoglobin and two mutants (H64L and H97F) were studied in the pH range from 4.2 to 9.5. Comparison of the native protein with the mutants shows that the observed pH effects can be traced to protonations of two histidine residues, H64 and H97, near the active site. Their imidazole sidechains experience simple, uncoupled Henderson-Hasselbalch type protonations, giving rise to four different protonation states. Because two of the protonation states are linked by a pH-independent equilibrium, the overall pH dependence of the spectra is described by a linear combination of three independent components. Global analysis, based on singular value decomposition and matrix least-squares algorithms enabled us to extract the pK values of the two histidines and the three basis spectra of the protonating species. The basis spectra were decomposed into the taxonomic substates A(0), A(1), and A(3), previously introduced in a heuristic way to analyze CO stretch spectra in heme proteins at fixed pH (see for instance, Biophys. J. 71:1563-1573). Moreover, an additional, weakly populated substate, called A(x), was identified. Protonation of H97 gives rise to a blue shift of the individual infrared lines by about 2 cm(-1), so that the A substates actually appear in pairs, such as A(0) and A(0)(+). The blue shift can be explained by reduced backbonding from the heme iron to the CO. Protonation of the distal histidine, H64, leads to a change of the infrared absorption from the A(1) or A(3) substate lines to A(0). This behavior can be explained by a conformational change upon protonation that moves the imidazole sidechain of H64 away from the CO into the high-dielectric solvent environment, which avoids the energetically unfavorable situation of an uncompensated electric charge in the apolar, low-dielectric protein interior. Our results suggest that protonation reactions serve as an important mechanism to create taxonomic substates in proteins.

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.

Structural and spectroscopic studies of azide complexes of horse heart myoglobin and the His-64-->Thr variant

The high-resolution X-ray crystallographic structures of horse heart azidometmyoglobin complexes of the wild-type protein and the His-64-->Thr variant have been determined to 2.0 and 1.8 A respectively. Azide binds to wild-type metmyoglobin in a bent configuration with an Fe-N-1-N-3 angle of 119 degrees and is oriented into the distal crevice in the direction of Ile-107. The proximity of the His-64 NE2 atom to the N-1 atom of the bound azide indicates stabilization of the ligand by the His-64 side chain through hydrogen bonding. In addition, structural characterization of wild-type horse heart azidometmyoglobin establishes that the only structural change induced by ligand binding is a small movement of the Leu-29 side chain away from the azide ligand. EPR and Fourier transform infrared spectroscopy were used to characterize the myoglobin azide complexes further. EPR spectroscopy revealed that, in contrast with wild-type azidometmyoglobin, two slightly different low-spin species are formed by azide bound to the His-64-->Thr variant both in solution and in a polycrystalline sample. One of these low-spin species has a greater relative intensity, with g values very similar to those of the azide complex of the wild-type protein. These EPR results together with structural information on this variant indicate the presence of two distinct conformations of bound azide, with one form predominating. The major conformation is comparable to that formed by wild-type myoglobin in which azide is oriented into the distal crevice. In the minor conformation the azide is oriented towards the exterior of the protein.

Theoretical study of the electrostatic and steric effects on the spectroscopic characteristics of the metal-ligand unit of heme proteins. 2. C-O vibrational frequencies, 17O isotropic chemical shifts, and nuclear quadrupole coupling constants

The quantum chemical calculations, vibronic theory of activation, and London-Pople approach are used to study the dependence of the C-O vibrational frequency, 17O isotropic chemical shift, and nuclear quadrupole coupling constant on the distortion of the porphyrin ring and geometry of the CO coordination, changes in the iron-carbon and iron-imidazole distances, magnitude of the iron displacement out of the porphyrin plane, and presence of the charged groups in the heme environment. It is shown that only the electrostatic interactions can cause the variation of all these parameters experimentally observed in different heme proteins, and the heme distortions could modulate this variation. The correlations between the theoretically calculated parameters are shown to be close to the experimentally observed ones. The study of the effect of the electric field of the distal histidine shows that the presence of the four C-O vibrational bands in the infrared absorption spectra of the carbon monoxide complexes of different myoglobins and hemoglobins can be caused by the different orientations of the different tautomeric forms of the distal histidine. The dependence of the 17O isotropic chemical shift and nuclear quadrupole coupling constant on pH and the distal histidine substitution can be also explained from the same point of view.

Effect of distal cavity mutations on the binding and activation of oxygen by ferrous horseradish peroxidase

Mutations have been introduced at residues Arg-38 or His-42 in horseradish peroxidase isoenzyme C (HRPC) in order to probe the role of these key distal residues in the reaction of ferrous HRPC with dioxygen. The association and dissociation rate constants for dioxygen binding to His-42 --> Leu, His-42 --> Arg, Arg-38 --> Leu, Arg-38 --> Lys, Arg-38 --> Ser, and Arg-38 --> Gly variants have been measured using stopped-flow spectrophotometry. Replacement of His-42 by Leu or Arg increases the oxygen binding rate constant by less than an order of magnitude, whereas changing the polar distal Arg-38 causes increases of more than 2 orders. These results demonstrate that His-42 and Arg-38 impede the binding of dioxygen to ferrous HRPC, presumably by steric and/or electrostatic interactions in the distal heme cavity. Recombinant HRPC oxyperoxidase reverted slowly to the ferric state with no spectrophotometrically detectable intermediates and with an apparent first-order rate constant of 9.0 x 10(-3) s(-1), which is essentially the same as that for the native, glycosylated enzyme. This reaction was accelerated when His-42 was replaced by Leu or Arg (k(decay) = 0.10 and 0.07 s(-1), respectively) presumably due to the loss of the hydrogen bond between the His-42 imidazole and the bound dioxygen. Substitution of Arg-38 by Leu, Lys, or Gly also produced a less stable oxyperoxidase (k(decay) = 0.22, 0.20, and 0.58 s(-1), respectively). However, with the Arg-38 --> Ser variant, a transient intermediate, proposed to be a ferric-superoxide complex, was detected by rapid-scan stopped-flow spectrophotometry during the conversion of oxyperoxidase to the ferric state. This variant also exhibits an unusually high affinity for dioxygen. It is proposed that Arg-38 interacts with the bound dioxygen to promote superoxide character, thereby stabilizing the oxyperoxidase state and making the binding of dioxygen to ferrous HRPC essentially irreversible. We conclude that Arg-38 and His-42 not only promote the heterolytic cleavage of bound hydrogen peroxide to form compound I but also decrease the lability of the ferrous enzyme-dioxygen complex in order to suppress the formation of the inactive ferrous state.

Theoretical study of the distal-side steric and electrostatic effects on the vibrational characteristics of the FeCO unit of the carbonylheme proteins and their models

The vibronic theory of activation and quantum chemical intermediate neglect of differential overlap (INDO) calculations are used to study the activation of carbon monoxide (change of the C-O bond index and force field constant) by the imidazole complex with heme in dependence on the distortion of the porphyrin ring, geometry of the CO coordination, iron-carbon and iron-imidazole distances, iron displacement out of the porphyrin plane, and presence of the charged groups in the heme environment. It is shown that the main contribution to the CO activation stems from the change in the sigma donation from the 5 sigma CO orbital to iron, and back-bonding from the iron to the 2 pi orbital of CO. It follows from the results that none of the studied distortions can explain, by itself, the wide variation of the C-O vibrational frequency in the experimentally studied model compounds and heme proteins. To study the dependence of the properties of the FeCO unit on the presence of charged groups in the heme environment, the latter are simulated by the homogeneous electric field and point charges of different magnitude and location. The results show that charged groups can strongly affect the strength of the C-O bond and its vibrational frequency. It is found that the charges located on the distal side of the heme plane can affect the Fe-C and C-O bond indexes (and, consequently, the Fe-C and C-O vibrational frequencies), both in the same and in opposite directions, depending on their position. The theoretical results allow us to understand the peculiarities of the effect of charged groups on the properties of the FeCO unit both in heme proteins and in their model compounds.

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.

Origin of the pH-dependent spectroscopic properties of pentacoordinate metmyoglobin variants

The pH dependence of the electronic and EPR spectra of two variants of horse heart myoglobin (Mb) in which the distal His64 ligand has been replaced by either Thr or Ile has been studied. Both of these variants exhibit spectroscopic changes with pH that are indicative of a transition between two ferric high-spin forms that occurs with a pKa of 9.49 for the His64Thr variant and 9.26 for the His64Ile variant and that is distinctly different from the pH-dependent spectroscopic changes related to titration of the distal aquo ligand of wild-type Mb. The electronic and EPR spectra of both variants at all values of pH studied are consistent with the presence of a pentacoordinate heme iron center. For the His64Thr variant, a high-resolution (1.9 A) structure determination establishes the lack of the distal aquo ligand and demonstrates an out-of-plane movement of the ferric iron toward the proximal histidine together with a decrease of the Fe-His bond length. Investigation of this pH-linked equilibrium by EPR spectroscopy reveals rhombically split high-spin signals at both pH 7 and 11 with a greater degree of rhombicity exhibited by the alkaline species. We propose that the pH-linked spectroscopic transition exhibited by these distal histidine variants results from the deprotonation of the proximal His93 residue to produce imidazolate ligation at alkaline pH.

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.

Distal residue-CO interaction in carbonmonoxy myoglobins: a molecular dynamics study of three distal mutants

Six 90-ps molecular dynamics trajectories, two for each of three distal mutants of sperm whale carbonmonoxy myoglobin, are reported; solvent waters within 16 A of the active site have been included. In both His64GIn trajectories, the distal side chain remains part of the heme pocket, forming a "closed" conformation similar to that of the wild type 64N delta H tautomer. Despite a connectivity more closely resembling the N epsilon H histidine tautomer, close interactions with the carbonyl ligand similar to those observed for the wild type 64N epsilon H tautomer are prevented in this mutant by repulsive interactions between the carbonyl O and the 64O epsilon. The aliphatic distal side chain of the His64Leu mutant shows little interaction with the carbonyl ligand in either His64Leu trajectory. Solvent water molecules move into and out of the active site in the His64Gly mutant trajectories; during all the other carbonmonoxy myoglobin trajectories, including the wild type distal tautomers considered in an earlier work, solvent molecules rarely encroach closer than 6 A of the active site. These results are consistent with a recent structural interpretation of the wild type infrared spectrum, and the current reinterpretation that the distal-ligand interaction in carbonmonoxy myoglobin is largely electrostatic, not steric, in nature.

Spectroscopic study of Ser92 mutants of human myoglobin: hydrogen bonding effect of Ser92 to proximal His93 on structure and property of myoglobin

Neutron diffraction studies have demonstrated that the hydroxyl group oxygen of Ser92(F7) is hydrogen bonded to the proximal His93(48) N epsilon H proton in myoglobin (Mb) [Cheng, X., & Shoenborn, B. P. (1991) J. Mol. Biol. 220, 381-399]. In order to examine the importance of this hydrogen bond, Ser92 was replaced with Ala and Asp in human Mb. By comparing the optical, 1H-NMR, resonance Raman, and IR spectra of Mb(S92A) in several spin and oxidation states with those of wild-type Mb, it was found that the mutation causes a structural change on the heme proximal side but not on the distal side. Comparison of the NMR spectra of the cyanomet form of Mb(S92A) and Mb(WT) suggests that the imidazole plane of His93 rotates somewhat around the Fe-N delta (His93) bond upon loss of the hydrogen bond between His93 and Ser92. The 2D 1H-NMR measurements of the CO complexes show that mutation of Ser92 to Ala changes the relative position of the His97 imidazole group to the heme plane, but the change is not so drastic as reported in the crystal data of Ser92 mutant of pig Mb [Smerdon et al. (1993) Biochemistry 32, 5132-5138]. On the other hand, ligand (CO, O2) binding is only slightly affected by this mutation. From these results, we conclude that the Ser92-His93 hydrogen bond maintains the protein structure of the proximal heme pocket, but it does not strongly affect the electronic structure of the heme as well as of the His93 imidazole ring.(ABSTRACT TRUNCATED AT 250 WORDS)

FTIR analysis of the interaction of azide with horse heart myoglobin variants

The interaction of azide with variants of horse heart myoglobin (Mb) has been characterized by Fourier transform infrared (FTIR), electron paramagnetic resonance (EPR), and UV-VIS absorption spectroscopy and by molecular modeling calculations. Distal histidine variants (His64Thr, His64Ile, His64Lys) and charged surface variants (Val67Arg, Lys45Glu, Lys45Glu/Lys63Glu) were included in this study. All variants, with the exception of Val67Arg, have a lower azide affinity than the wild-type protein. Analysis of the temperature dependence of the FTIR spectra (277-313 K) revealed that the wild-type protein and all variants exhibit a high-spin/low-spin equilibrium. Introduction of positively charged amino acid residues shifts nu max for the low-spin form to higher energy while negatively charged residues shifted this maximum to lower energy. The low azide binding affinity exhibited by the His64Thr and His64Ile variants is accompanied by a shift of the nu max for the low-spin infrared band to lower energy and by a significant increase in the corresponding half-bandwidths. This observation indicates greater mobility of the bound azide ligand in these variants. The His64Lys variant exhibits two infrared bands attributable to low-spin forms that are assigned to two different conformations of the lysyl residue. In one conformation, the lysine is proposed to form a hydrogen bond with the bound azide similar to that proposed to occur between the distal histidine and bound azide, and in the other conformation no interaction occurs.

Anatomy and dynamics of a ligand-binding pathway in myoglobin: the roles of residues 45, 60, 64, and 68

In order for diatomic ligands to enter and exit myoglobin, there must be substantial displacements of amino acid side chains from their positions in the static X-ray structure. One pathway, involving Arg/Lys45, His64, and Val68, has been studied in greatest detail. In an earlier study (Lambright et al., 1989) we reported the surprising result that mutation of the surface residue Lys45 to arginine lowers the inner barrier to CO rebinding. Until then, it had been thought that this barrier primarily involves interior distal pocket residues such as His64 and Val68. In this report, we present a detailed study of the CO rebinding kinetics in aqueous solution of a series of single- and double-site mutants of human myoglobin at positions 64, 68, 45, and 60. On the basis of the observed kinetics, we propose that the effect of surface residue 45 on the inner barrier can be explained by a chain of interactions between surface and pocket residues. Very large, and in some cases unexpected, changes are observed in the kinetics of recombination and in the partitioning between geminate and bimolecular recombination.

Site-directed mutants of the cytochrome bo ubiquinol oxidase of Escherichia coli: amino acid substitutions for two histidines that are putative CuB ligands

The bo-type ubiquinol oxidase of Escherichia coli is a member of the superfamily of structurally related heme-copper respiratory oxidases. The members of this family, which also includes the aa3-type cytochrome c oxidases, contain at least two heme prosthetic groups, a six-coordinate low-spin heme, and a high-spin heme. The high-spin heme is magnetically coupled to a copper, CuB, forming a binuclear center which is the site of oxygen reduction to water. Vectorial proton translocation across the membrane bilayer appears to be another common feature of this superfamily of oxidases. It has been proposed previously that the two adjacent histidines in putative transmembrane helix VII (H333 and H334 in the E. coli sequence) of the largest subunit of the heme-copper oxidases are ligands to CuB. Previously reported mutagenesis studies of the E. coli bo-type oxidase and the aa3-type oxidase of Rhodobacter sphaeroides supported this model, as substitutions at these two positions produced nonfunctional enzymes but did not perturb the visible spectra of the two heme groups. In this work, six different amino acids, including potential copper-liganding residues, were substituted for H333 and H334 of the E. coli oxidase. All of the mutations resulted in inactive, but assembled, oxidase with both of the heme components present. However, cryogenic Fourier transform infrared (FTIR) spectroscopy of the CO adducts revealed that dramatic changes occur at the binuclear center as a result of each mutation and that CuB appears to be absent.(ABSTRACT TRUNCATED AT 250 WORDS)

Dynamics of protein relaxation in site-specific mutants of human myoglobin

We have recently reported spectroscopic evidence for structural relaxation of myoglobin (Mb) following photodissociation of MbCO [Lambright, D. G., Balasubramanian, S., & Boxer, S. G. (1991) Chem. Phys. 158, 249-260]. In this paper we report measurements for a series of single amino acid mutants of human myoglobin on the distal side of the heme pocket (positions 45, 64, and 68) in order to examine specific structural determinants involved in this conformational relaxation and to determine the nature of the coupling between relaxation and the functional process of ligand binding. The kinetics of ligand binding and conformational relaxation were monitored by transient absorption spectroscopy in the Soret spectral region, and the results are analyzed using a four-state ligand binding model. Two principal results emerge: (1) amino acid substitutions in the distal heme pocket affect the kinetics of the nonequilibrium conformational relaxation and (2) the rate of ligand escape from the protein matrix is not significantly perturbed by the distal heme pocket mutations.