Perturbation of the Fe−O2 Bond by Nearby Residues in Heme Pocket: Observation of νFe-O2 Raman Bands for Oxymyoglobin Mutants

Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies

After a general introduction to the features and mechanisms of cytochrome c oxidases (CcOs) in mitochondria and aerobic bacteria, we present DFT calculated physical and spectroscopic properties for the catalytic reaction cycle compared with experimental observations in bacterial ba3 type CcO, also with comparisons/contrasts to aa3 type CcOs. The Dinuclear Complex (DNC) is the active catalytic reaction center, containing a heme a3 Fe center and a near lying Cu center (called CuB) where by successive reduction and protonation, molecular O2 is transformed to two H2O molecules, and protons are pumped from an inner region across the membrane to an outer region by transit through the CcO integral membrane protein. Structures, energies and vibrational frequencies for Fe-O and O-O modes are calculated by DFT over the catalytic cycle. The calculated DFT frequencies in the DNC of CcO are compared with measured frequencies from Resonance Raman spectroscopy to clarify the composition, geometry, and electronic structures of different intermediates through the reaction cycle, and to trace reaction pathways. X-ray structures of the resting oxidized state are analyzed with reference to the known experimental reaction chemistry and using DFT calculated structures in fitting observed electron density maps. Our calculations lead to a new proposed reaction pathway for coupling the PR → F → OH (ferryl-oxo → ferric-hydroxo) pathway to proton pumping by a water shift mechanism. Through this arc of the catalytic cycle, major shifts in pKa's of the special tyrosine and a histidine near the upper water pool activate proton transfer. Additional mechanisms for proton pumping are explored, and the role of the CuB+ (cuprous state) in controlling access to the dinuclear reaction site is proposed.

A Water Molecule Residing in the Fea33+···CuB2+ Dinuclear Center of the Resting Oxidized as-Isolated Cytochrome c Oxidase: A Density Functional Study

Although the dinuclear center (DNC) of the resting oxidized "as-isolated" cytochrome c oxidase (CcO) is not a catalytically active state, its detailed structure, especially the nature of the bridging species between the Fea33+ and CuB2+ metal sites, is still both relevant and unsolved. Recent crystallographic work has shown an extended electron density for a peroxide type dioxygen species (O1-O2) bridging the Fea3 and CuB centers. In this paper, our density functional theory (DFT) calculations show that the observed peroxide type electron density between the two metal centers is most likely a mistaken analysis due to overlap of the electron density of a water molecule located at different positions between apparent O1 and O2 sites in DNCs of different CcO molecules with almost the same energy. Because the diffraction pattern and the resulting electron density map represent the effective long-range order averaged over many molecules and unit cells in the X-ray structure, this averaging can lead to an apparent observed superposition of different water positions between the Fea33+ and CuB2+ metal sites.

DFT Fea3-O/O-O Vibrational Frequency Calculations over Catalytic Reaction Cycle States in the Dinuclear Center of Cytochrome c Oxidase

Density functional vibrational frequency calculations have been performed on eight geometry optimized cytochrome c oxidase (CcO) dinuclear center (DNC) reaction cycle intermediates and on the oxymyoglobin (oxyMb) active site. The calculated Fe-O and O-O stretching modes and their frequency shifts along the reaction cycle have been compared with the available resonance Raman (rR) measurements. The calculations support the proposal that in state A[Fea33+-O2-•···CuB+] of CcO, O2 binds with Fea32+ in a similar bent end-on geometry to that in oxyMb. The calculations show that the observed 20 cm-1 shift of the Fea3-O stretching mode from the PR to F state is caused by the protonation of the OH- ligand on CuB2+ (PR[Fea34+═O2-···HO--CuB2+] → F[Fea34+═O2-···H2O-CuB2+]), and that the H2O ligand is still on the CuB2+ site in the rR identified F[Fea34+═O2-···H2O-CuB2+] state. Further, the observed rR band at 356 cm-1 between states PR and F is likely an O-Fea3-porphyrin bending mode. The observed 450 cm-1 low Fea3-O frequency mode for the OH active oxidized state has been reproduced by our calculations on a nearly symmetrically bridged Fea33+-OH-CuB2+ structure with a relatively long Fea3-O distance near 2 Å. Based on Badger's rule, the calculated Fea3-O distances correlate well with the calculated νFe-O-2/3 (νFe-O is the Fea3-O stretching frequency) with correlation coefficient R = 0.973.

Hydrogen atom abstraction by synthetic heme ferric superoxide and hydroperoxide species

To date, artificial dioxygen adducts of heme have not been demonstrated to be able to oxidize organic substrates in sharp contrast to their non-heme analogues and naturally occurring enzymes like heme dioxygenases. To address this apparent anomaly, an iron porphyrin complex is synthesized which includes a pendant quinol group. The corresponding dioxygen bound iron porphyrin species is demonstrated to perform hydrogen atom transfer (HAT) from the quinol group appended to the porphyrin ligand. The resultant ferric peroxide, formed by the first HAT, performs a 2nd HAT generating a ferryl species (FeIV[double bond, length as m-dash]O) and resulting in the 2e-/2H+ oxidation of the pendant hydroquinone to quinone.

Bridging Theory and Experiment to Address Structural Properties of Truncated Haemoglobins: Insights from Thermobifida fusca HbO

In this chapter, we will discuss the paradigmatic case of Thermobifida fusca (Tf-trHb) HbO in its ferrous and ferric states and its behaviour towards a battery of possible ligands. This choice was dictated by the fact that it has been one of the most extensively studied truncated haemoglobins, both in terms of spectroscopic and molecular dynamics studies. Tf-trHb typifies the structural properties of group II trHbs, as the active site is characterized by a highly polar distal environment in which TrpG8, TyrCD1, and TyrB10 provide three potential H-bond donors in the distal cavity capable of stabilizing the incoming ligands. The role of these residues in key topological positions, and their interplay with the iron-bound ligands, has been addressed in studies carried out on the CO, F(-), OH(-), CN(-), and HS(-) adducts formed with the wild-type protein and a combinatorial set of mutants, in which the distal polar residues, TrpG8, TyrCD1, and TyrB10, have been singly, doubly, or triply replaced by a Phe residue. In this context, such a complete analysis provides an excellent benchmark for the investigation of the relationship between protein structure and function, allowing one to translate physicochemical properties of the active site into the observed functional behaviour. Tf-trHb will be compared with other members of the group II trHbs and, more generally, with members of the other trHb subgroups.

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

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

NMR, IR, Mössbauer and quantum chemical investigations of metalloporphyrins and metalloproteins

We review contributions made towards the elucidation of CO and O2binding geometries in respiratory proteins. Nuclear magnetic resonance, infrared spectroscopy, Mössbauer spectroscopy, X-ray crystallography and quantum chemistry have all been used to investigate the Fe –ligand interactions. Early experimental results showed linear correlations between17O chemical shifts and the infrared stretching frequency (νCO) of the CO ligand in carbonmonoxyheme proteins and between the17O chemical shift and the13CO shift. These correlations led to early theoretical investigations of the vibrational frequency of carbon monoxide and of the13C and17O NMR chemical shifts in the presence of uniform and non-uniform electric fields. Early success in modeling these spectroscopic observables then led to the use of computational methods, in conjunction with experiment, to evaluate ligand-binding geometries in heme proteins. Density functional theory results are described which predict57Fe chemical shifts and Mössbauer electric field gradient tensors,17O NMR isotropic chemical shifts, chemical shift tensors and nuclear quadrupole coupling constants (e2qQ/h) as well as13C isotropic chemical shifts and chemical shift tensors in organometallic clusters, heme model metalloporphyrins and in metalloproteins. A principal result is that CO in most heme proteins has an essentially linear and untilted geometry (τ = 4 °, β = 7 °) which is in extremely good agreement with a recently published X-ray synchrotron structure. CO / O2discrimination is thus attributable to polar interactions with the distal histidine residue, rather than major Fe–C–O geometric distortions.

Spectroscopic features of cytochrome P450 reaction intermediates

Cytochromes P450 constitute a broad class of heme monooxygenase enzymes with more than 11,500 isozymes which have been identified in organisms from all biological kingdoms [1]. These enzymes are responsible for catalyzing dozens chemical oxidative transformations such as hydroxylation, epoxidation, N-demethylation, etc., with very broad range of substrates [2,3]. Historically these enzymes received their name from 'pigment 450' due to the unusual position of the Soret band in UV-vis absorption spectra of the reduced CO-saturated state [4,5]. Despite detailed biochemical characterization of many isozymes, as well as later discoveries of other 'P450-like heme enzymes' such as nitric oxide synthase and chloroperoxidase, the phenomenological term 'cytochrome P450' is still commonly used as indicating an essential spectroscopic feature of the functionally active protein which is now known to be due to the presence of a thiolate ligand to the heme iron [6]. Heme proteins with an imidazole ligand such as myoglobin and hemoglobin as well as an inactive form of P450 are characterized by Soret maxima at 420nm [7]. This historical perspective highlights the importance of spectroscopic methods for biochemical studies in general, and especially for heme enzymes, where the presence of the heme iron and porphyrin macrocycle provides rich variety of specific spectroscopic markers available for monitoring chemical transformations and transitions between active intermediates of catalytic cycle.

Formation of an end-on ferric peroxo intermediate upon one-electron reduction of a ferric superoxo heme

The low-spin end-on ferric peroxo heme intermediate has been proposed as an alternative reactive intermediate involved in the catalytic cycles of enzymes such as nitric oxide synthase and cytochrome P450. This transient heme intermediate has never been captured using synthetic heme models. We demonstrate herein our success in the solution preparation of such an end-on ferric peroxo intermediate derived from a heme model, which features both a group hanging over the porphyrin macrocycle and a covalently appended axial imidazole ligand, through one-electron reduction of its ferric superoxo precursor. The obtained ferric peroxo intermediate was further transformed into the corresponding ferric hydroperoxo species upon protonation. This heme model compound provides a convenient system for sequential preparation of the important and biologically relevant superoxo/peroxo/hydroperoxo heme intermediates through an oxygenation/one-electron reduction/protonation process similar to the mechanisms used by enzyme systems.

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.

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.

Structure and ligand selection of hemoglobin II from Lucina pectinata

Lucina pectinata ctenidia harbor three heme proteins: sulfide-reactive hemoglobin I (HbI(Lp)) and the oxygen transporting hemoglobins II and III (HbII(Lp) and HbIII(Lp)) that remain unaffected by the presence of H(2)S. The mechanisms used by these three proteins for their function, including ligand control, remain unknown. The crystal structure of oxygen-bound HbII(Lp) shows a dimeric oxyHbII(Lp) where oxygen is tightly anchored to the heme through hydrogen bonds with Tyr(30)(B10) and Gln(65)(E7). The heme group is buried farther within HbII(Lp) than in HbI(Lp). The proximal His(97)(F8) is hydrogen bonded to a water molecule, which interacts electrostatically with a propionate group, resulting in a Fe-His vibration at 211 cm(-1). The combined effects of the HbII(Lp) small heme pocket, the hydrogen bonding network, the His(97) trans-effect, and the orientation of the oxygen molecule confer stability to the oxy-HbII(Lp) complex. Oxidation of HbI(Lp) Phe(B10) --> Tyr and HbII(Lp) only occurs when the pH is decreased from pH 7.5 to 5.0. Structural and resonance Raman spectroscopy studies suggest that HbII(Lp) oxygen binding and transport to the host bacteria may be regulated by the dynamic displacements of the Gln(65)(E7) and Tyr(30)(B10) pair toward the heme to protect it from changes in the heme oxidation state from Fe(II) to Fe(III).

The roles of Tyr(CD1) and Trp(G8) in Mycobacterium tuberculosis truncated hemoglobin O in ligand binding and on the heme distal site architecture

The crystal structure of the cyano-met form of Mt-trHbO revealed two unusual distal residues Y(CD1) and W(G8) forming a hydrogen-bond network with the heme-bound ligand [Milani, M., et al. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 5766-5771]. W(G8) is an invariant residue in group II and group III trHbs and has no counterpart in other globins. A previous study reported that changing Y(CD1) for a Phe causes a significant increase in the O2 combination rate, but almost no change in the O2 dissociation rate [Ouellet, H., et al. (2003) Biochemistry 42, 5764-5774]. Here we investigated the role of the W(G8) in ligand binding by using resonance Raman spectroscopy, stopped-flow spectrophotometry, and X-ray crystallography. For this purpose, W(G8) was changed, by site-directed mutagenesis, to a Phe in both the wild-type protein and the mutant Y(CD1)F to create the single mutant W(G8)F and the double mutant Y(CD1)F/W(G8)F, respectively. Resonance Raman results suggest that W(G8) interacts with the heme-bound O2 and CO, as evidenced by the increase of the Fe-O2 stretching mode from 559 to 564 cm-1 and by the lower frequency of the Fe-CO stretching modes (514 and 497 cm-1) compared to that of the wild-type protein. Mutation of W(G8) to Phe indicates that this residue controls ligand binding, as evidenced by a dramatic increase of the combination rates of both O2 and CO. Also, the rate of O2 dissociation showed a 90-1000-fold increase in the W(G8)F and Y(CD1)F/W(G8)F mutants, that is in sharp contrast with the values obtained for the other distal mutants Y(B10)F and Y(CD1)F [Ouellet, H., et al. (2003) Biochemistry 42, 5764-5774]. Taken together, these data indicate a pivotal role for the W(G8) residue in O2 binding and stabilization.

Structural and Functional Properties of a Truncated Hemoglobin from a Food-borne Pathogen Campylobacter jejuni

Campylobacter jejuni contains two hemoglobins, Cgb and Ctb. Cgb has been suggested to perform an NO detoxification reaction to protect the bacterium against NO attack. On the other hand, the physiological function of Ctb, a class III truncated hemoglobin, remains unclear. By using CO as a structural probe, resonance Raman data show that the distal heme pocket of Ctb exhibits a positive electrostatic potential. In addition, two ligand-related vibrational modes, nu(Fe-O(2)) and nu(O-O), were identified in the oxy derivative, with frequencies at 542 and 1132 cm(-1), respectively, suggesting the presence of an intertwined H-bonding network surrounding the heme-bound ligand, which accounts for its unusually high oxygen affinity (222 microm(-1)). Mutagenesis studies of various distal mutants suggest that the heme-bound dioxygen is stabilized by H-bonds donated from the Tyr(B10) and Trp(G8) residues, which are highly conserved in the class III truncated hemoglobins; furthermore, an additional H-bond donated from the His(E7) to the Tyr(B10) further regulates these H-bonding interactions by restricting the conformational freedom of the phenolic side chain of the Tyr(B10). Taken together, the data suggest that it is the intricate balance of the H-bonding interactions that determines the unique ligand binding properties of Ctb. The extremely high oxygen affinity of Ctb makes it unlikely to function as an oxygen transporter; on the other hand, the distal heme environment of Ctb is surprisingly similar to that of cytochrome c peroxidase, suggesting a role of Ctb in performing a peroxidase or P450-type of oxygen chemistry.

Roles of the Heme Distal Residues of FixL in O2 Sensing: A Single Convergent Structure of the Heme Moiety Is Relevant to the Downregulation of Kinase Activity

FixL is a heme-based O(2) sensor, in which the autophosphorylation is regulated by the binding of exogenous ligands such as O(2) and CN(-). In this study, mutants of the heme distal Arg200, Arg208, Ile209, Ile210, and Arg214 residues of SmFixL were characterized biochemically and physicochemically, because it has been suggested that they are significant residues in ligand-linked kinase regulation. Measurements of the autoxidation rate, affinities, and kinetics of ligand binding revealed that all of the above residues are involved in stabilization of the O(2)-heme complex of FixL. However, Arg214 was found to be the only residue that is directly relevant to the ligand-dependent regulation of kinase activity. Although the wild type and R214K and R214Q mutants exhibited normal kinase regulation, R214A, R214M, R214H, and R214Y did not. (13)C and (15)N NMR analyses for (13)C(15)N(-) bound to the truncated heme domains of the Arg214 mutants indicated that, in the wild type and the foregoing two mutants, the heme moiety is present in a single conformation, but in the latter four, the conformations fluctuate possibly because of the lack of an interaction between the iron-bound ligand and residue 214. It is likely that such a rigid conformation of the ligand-bound form is important for the downregulation of histidine kinase activity. Furthermore, a comparison of the NMR data between the wild type and R214K and R214Q mutants suggests that a strong electrostatic interaction between residue 214 and the iron-bound ligand is not necessarily required for the single convergent structure and eventually for the downregulation of FixL.

Mechanism for transduction of the ligand-binding signal in heme-based gas sensory proteins revealed by resonance Raman spectroscopy

Gene analysis has revealed a variety of new heme-containing gas sensory proteins in organisms ranging from bacteria to mammals. These proteins are composed of sensor, communication, and functional domains. The sensor domain contains a heme that binds effector molecules such as NO, O2, or CO. Ligand binding by the sensor domain modulates the physiological role of the protein, such as DNA binding in the case of transcriptional factors or the catalytic reaction rate in the case of enzymes. This Account summarizes resonance Raman (RR) studies, including static and time-resolved measurements, which have enabled elucidation of the mechanisms by which binding of specific target molecule by the sensor domain is transduced to alteration of the functional domain. These studies have shown that signals can be conveyed from the heme to the functional domain via three different pathways: (i) a distal pathway, (ii) a proximal pathway, and (iii) a heme peripheral pathway.

Resonance Raman study of deoxy and ligated (O2 and CO) mesoheme IX-reconstituted myoglobin, hemoglobin and its alpha and beta subunits

In this work, we corrected the resonance Raman (RR) results presented earlier for deoxy mesoheme IX-reconstituted hemoglobin (mesoHb) alpha and beta subunits implied that mesohemes in these subunits undergo substantial structural changes upon formation of a hemoglobin tetramer (Biochemistry 29 (1990) 5087). We show that these data were probably due to the improper handling of the deoxy mesoheme subunit preparation. Additionally, we discuss the RR spectra of deoxy, oxy, and CO species of mesoheme IX-reconstituted myoglobin (mesoMb) and alpha and beta deoxy meso hemoglobin subunits, including their analogues with deuterium-substituted mesoheme IX in all methyl groups (d(12)). Based on the obtained data, we propose a complete RR band assignment for all of the investigated molecules. The most pronounced changes are observed for the gamma(7) mode (out-of-plane movement of methane carbon atoms) associated with the interaction of the ethyl groups with the globin. We also show that in mesoheme IX-reconstituted proteins, the O(2) molecule binds stronger than in the case of native species. This is manifested by the up-shift of nu(Fe-O(2)).

Role of Arginine 220 in the Oxygen Sensor FixL from Bradyrhizobium japonicum

In the heme-based oxygen sensor protein FixL, conformational changes induced by oxygen binding to the heme sensor domain regulate the activity of a neighboring histidine kinase, eventually restricting expression of specific genes to hypoxic conditions. The conserved arginine 220 residue is suggested to play a key role in the signal transduction mechanism. To obtain detailed insights into the role of this residue, we replaced Arg(220) by histidine (R220H), glutamine (R220Q), glutamate (R220E), and isoleucine (R220I) in the heme domain FixLH from Bradyrhizobium japonicum. These mutations resulted in dramatic changes in the O(2) affinity with K(d) values in the order R220I < R220Q < wild type < R220H. For the R220H and R220Q mutants, residue 220 interacts with the bound O(2) or CO ligands, as seen by resonance Raman spectroscopy. For the oxy-adducts, this H-bond modifies the pi acidity of the O(2) ligand, and its strength is correlated with the back-bonding-sensitive nu(4) frequency, the k(off) value for O(2) dissociation, and heme core-size conformational changes. This effect is especially strong for the wild-type protein where Arg(220) is, in addition, positively charged. These observations strongly suggest that neither strong ligand fixation nor the displacement of residue 220 into the heme distal pocket are solely responsible for the reported heme conformational changes associated with kinase activity regulation, but that a significant decrease of the heme pi(*) electron density because of strong back-bonding toward the oxygen ligand also plays a key role.

Resonance Raman Investigation of the Specific Sensing Mechanism of a Target Molecule by Gas Sensory Proteins

Specific sensing of gas molecules such as CO, NO, and O2 is a unique function of gas sensory hemoproteins, while hemoproteins carry out a wide variety of functions such as oxygen storage/transport, electron transfer, and catalysis as enzymes. It is important in gas sensory proteins that the heme domain not only recognizes its target molecule but also discriminates against other gases having similar molecular structures. Coordination of a target molecule to the heme is assumed to alter the protein conformation in the vicinity of heme, and the conformation change is propagated to the effector domain where substrate turnover, DNA binding, or interaction with a signal transduction protein is performed differently than the binding of other gases. To understand the appearance of such a specificity, we focus our attention on the ligand-protein interactions in the distal side of heme. In practice, the metal-ligand vibrations as well as internal modes of ligand and heme are measured with resonance Raman spectroscopy for wild-type and some mutant proteins with full-length or limited sensory regions. On the basis of such observations together with the knowledge currently available, we discuss the mechanism of specific sensing of a diatomic molecule in gas sensory proteins.

Structure and function of the Gondwanian hemoglobin of Pseudaphritis urvillii, a primitive notothenioid fish of temperate latitudes

The suborder Notothenioidei dominates the Antarctic ichthyofauna. The non-Antarctic monotypic family Pseudaphritidae is one of the most primitive families. The characterization of the oxygen-transport system of euryhaline Pseudaphritis urvillii is herewith reported. Similar to most Antarctic notothenioids, this temperate species has a single major hemoglobin (Hb 1, over 95% of the total). Hb 1 has strong Bohr and Root effects. It shows two very uncommon features in oxygen binding: At high pH values, the oxygen affinity is exceptionally high compared to other notothenioids, and subunit cooperativity is modulated by pH in an unusual way, namely the curve of the Hill coefficient is bell-shaped, with values approaching 1 at both extremes of pH. Molecular modeling, electronic absorption and resonance Raman spectra have been used to characterize the heme environment of Hb 1 in an attempt to explain these features, particularly in view of some potentially important nonconservative replacements found in the primary structure. Compared to human HbA, no major changes were found in the structure of the proximal cavity of the alpha-chain of Hb 1, although an altered distal histidyl and heme position was identified in the models of the beta-chain, possibly facilitated by a more open heme pocket due to reduced steric constraints on the vinyl substituent groups. This conformation may lead to the hemichrome form identified by spectroscopy in the Met state, which likely fulfils a potentially important physiological role.

Distal heme pocket conformers of carbonmonoxy derivatives of Ascaris hemoglobin: evidence of conformational trapping in porous sol-gel matrices

We report the ligand dependence of the conformer distribution in the distal heme pocket of Ascaris suum hemoglobin (Hb) studied by resonance Raman spectroscopy. The heme-bound CO is used as a spectroscopic antenna to probe the original distribution of conformers in the dioxygen derivative of Ascaris Hb, by utilizing sol-gel encapsulation. The first step is to encapsulate the dioxygen derivative in the porous sol-gel and let the gel age, thus trapping the equilibrium conformational distribution of Ascaris dioxygen Hb. In the second step, the dioxygen ligand is replaced by CO. The sol-gel environment impedes any large scale movements, drastically slowing down the conformational relaxation triggered by the ligation change, essentially "locking in" the initial quaternary and even tertiary structure of the protein. Studying the Fe-CO frequencies of the latter sample allows evaluation of the distribution of the distal heme pocket conformers that was originally associated with the dioxygen derivative. Extending the study to the Ascaris mutants allows for examination of the effect of specific residues in the distal pocket on the conformational distribution. The choice of mutants was largely based on the anticipated variation in hydrogen bonding patterns. The results show that the sol-gel encapsulation can slow or prevent re-equilibration within the distal heme pocket of Ascaris Hb and that the distribution of distal heme pocket conformers for the CO derivative of Ascaris Hb in the sol-gel is highly dependent on the history of the sample. Additionally, we report a detailed study of the CO complex of the mutants in solution for assignment of the various heme pocket conformers, and we present a comparison of the sol-gel data with solution data. The results support a picture in which the dioxygen derivative biases the population strongly toward a tightly packed configuration that favors the network of strong hydrogen bonding interactions, and suggest that Ascaris Hb is uniquely designed for dioxygen capture.

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.

Structural dynamics of myoglobin: spectroscopic and structural characterization of ligand docking sites in myoglobin mutant L29W

We have studied CO binding to the heme and CO migration among protein internal cavities after photodissociation in sperm whale carbonmonoxy myoglobin (MbCO) mutant L29W using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) and kinetic experiments at cryogenic temperatures. Photoproduct intermediates, characterized by CO at particular locations in the protein, were selectively enhanced by applying special laser illumination protocols. These studies were performed on the L29W mutant protein and a series of double mutants constructed so that bulky amino acid side chains block passageways between cavities or fill these sites. Binding of xenon was also employed as an alternative means of occluding cavities. All mutants exhibit two conformations, A(I) and A(II), with distinctly different photoproduct states and ligand binding properties. These differences arise mainly from different positions of the W29 and H64 side chains in the distal heme pocket [Ostermann, A., et al. (2000) Nature 404, 205-208]. The detailed knowledge of the interplay between protein structure, protein dynamics, and ligand migration at cryogenic temperatures allowed us to develop a dynamic model that explains the slow CO and O(2) bimolecular association observed after flash photolysis at ambient temperature.

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.

Stationary and time-resolved resonance Raman spectra of His77 and Met95 mutants of the isolated heme domain of a direct oxygen sensor from Escherichia coli

The heme environments of Met(95) and His(77) mutants of the isolated heme-bound PAS domain (Escherichia coli DOS PAS) of a direct oxygen sensing protein from E. coli (E. coli DOS) were investigated with resonance Raman (RR) spectroscopy and compared with the wild type (WT) enzyme. The RR spectra of both the reduced and oxidized WT enzyme were characteristic of six-coordinate low spin heme complexes from pH 4 to 10. The time-resolved RR spectra of the photodissociated CO-WT complex had an iron-His stretching band (nu(Fe-His)) at 214 cm(-1), and the nu(Fe-CO) versus nu(CO) plot of CO-WT E. coli DOS PAS fell on the line of His-coordinated heme proteins. The photodissociated CO-H77A mutant complex did not yield the nu(Fe-His) band but gave a nu(Fe-Im) band in the presence of imidazole. The RR spectrum of the oxidized M95A mutant was that of a six-coordinate low spin complex (i.e. the same as that of the WT enzyme), whereas the reduced mutant appeared to contain a five-coordinate heme complex. Taken together, we suggest that the heme of the reduced WT enzyme is coordinated by His(77) and Met(95), and that Met(95) is displaced by CO and O(2). Presumably, the protein conformational change that occurs upon exchange of an unknown ligand for Met(95) following heme reduction may lead to activation of the phosphodiesterase domain of E. coli DOS.

Resonance Raman and ligand binding studies of the oxygen-sensing signal transducer protein HemAT from Bacillus subtilis

HemAT-Bs is a heme-containing signal transducer protein responsible for aerotaxis of Bacillus subtilis. The recombinant HemAT-Bs expressed in Escherichia coli was purified as the oxy form in which oxygen was bound to the ferrous heme. Oxygen binding and dissociation rate constants were determined to be k(on) = 32 microm(-1) s(-1) and k(off) = 23 s(-1), respectively, revealing that HemAT-Bs has a moderate oxygen affinity similar to that of sperm whale myoglobin (Mb). The rate constant for autoxidation at 37 degrees C was 0.06 h(-1), which is also close to that of Mb. Although the electronic absorption spectra of HemAT-Bs were similar to those of Mb, HemAT-Bs showed some unique characteristics in its resonance Raman spectra. Oxygen-bound HemAT-Bs gave the nu(Fe-O(2)) band at a noticeably low frequency (560 cm(-1)), which suggests a unique hydrogen bonding between a distal amino acid residue and the proximal atom of the bound oxygen molecule. Deoxy HemAT-Bs gave the nu(Fe-His) band at a higher frequency (225 cm(-1)) than those of ordinary His-coordinated deoxy heme proteins. CO-bound HemAT-Bs gave the nu(Fe-CO) and nu(C-O) bands at 494 and 1964 cm(-1), respectively, which fall on the same nu(C-O) versus nu(Fe-CO) correlation line as that of Mb. Based on these results, the structural and functional properties of HemAT-Bs are discussed.

A comparative resonance Raman analysis of heme-binding PAS domains: heme iron coordination structures of the BjFixL, AxPDEA1, EcDos, and MtDos proteins

The heme-PAS is a specialized domain with which a broad class of signal-transducing heme proteins detect physiological heme ligands. Such domains exhibit a wide range of ligand binding parameters, yet they are all expected to feature an alpha-beta heme binding fold and a predominantly hydrophobic heme distal pocket without a distal histidine. We have compared, for the first time, the resonance Raman spectra of several heme-PASs: the heme-binding domains of Bradyrhizobium japonicum FixL, Escherichia coli Dos, Acetobacter xylinum PDEA1, and Methanobacterium thermoautotrophicum Dos. In all cases, the nu(Fe)-(CO) and nu(C-O) values of the carbonmonoxy forms were consistent with coordination of the heme iron to histidine on the proximal side and binding of the CO without electrostatic interaction with the heme distal pocket. EcDos was unusual in having predominantly hexacoordinate heme iron in the deoxy and met forms. Despite an evident lack of CO interaction with the EcDos heme pocket, relatively low Fe-O(2) (562 cm(-1)) and N-O (1576 cm(-1)) stretching frequencies indicated that strong polar interactions with that heme distal pocket are possible for highly bent ligands such as O(2) or NO. None of the newly studied NO adducts exhibited evidence of the Fe-His rupture and pentacoordination previously noted for Sinorhizobium meliloti FixL. A low Fe-His stretching frequency, formerly interpreted as a strained Fe-His bond, and the slow association of O(2) with S. meliloti FixL failed to correlate with the newly studied proteins having low association rate or low equilibrium association constants for binding of O(2). We conclude that although heme-PASs share some features, they represent distinct signal transduction mechanisms.

Unique ligand-protein interactions in a new truncated hemoglobin from Mycobacterium tuberculosis

A new truncated hemoglobin (HbO) from Mycobacterium tuberculosis has been expressed and purified. Sequence alignment of HbO with other hemoglobins suggests that the proximal F8 residue is histidine and the distal E7 and the B10 positions are occupied by alanine and tyrosine, respectively. The highly conserved residue at the CD1 position, surprisingly, is tyrosine, making HbO the first exception in the hemoglobin family that does not contain phenylalanine at this position. Resonance Raman data suggest that a strong hydrogen bonding network, involving the B10 Tyr and the CD1 Tyr, stabilizes the heme-bound O2 and CO as evidenced by the relatively low frequency of the Fe-O2 stretching mode (559 cm(-1)) and the high frequency of the Fe-CO stretching mode (527 cm(-1)). The presence of this hydrogen bonding network is supported by mutagenesis studies with the B10 tyrosine or the CD1 tyrosine mutated to phenylalanine. Taken together, these data demonstrate a rigid and polar distal pocket in HbO, which is significantly different from that of HbN, the other hemoglobin from M. tuberculosis. The distinct features in the heme active site structures and the temporal expression patterns of HbO and HbN suggest that these two hemoglobins may have very different physiological functions.

Controlling ligand binding in myoglobin by mutagenesis

A quadruple mutant of sperm whale myoglobin was constructed to mimic the structure found in Ascaris suum hemoglobin. The replacements include His(E7)-->Gln, Leu(B10)-->Tyr, Thr(E10)--> Arg, and Ile(G8)-->Phe. Single, double, and triple mutants were characterized to dissect out the effects of the individual substitutions. The crystal structures of the deoxy and oxy forms of the quadruple mutant were determined and compared with that of native Ascaris hemoglobin. Tyr(B10) myoglobin displays low O(2) affinity, high dissociation rate constants, and heterogeneous kinetic behavior, suggesting unfavorable steric interactions between the B10 phenol side chain and His(E7). In contrast, all mutants containing the Tyr(B10)/Gln(E7) pair show high O(2) affinity, low dissociation rate constants, and simple, monophasic kinetic behavior. Replacement of Ile(107) with Phe enhances nanosecond geminate recombination singly and in combination with the Tyr(B10)/Gln(E7)/Arg(E10) mutation by limiting access to the Xe4 site. These kinetic results and comparisons with native Ascaris hemoglobin demonstrate the importance of distal pocket cavities in governing the kinetics of ligand binding. The approximately 150-fold higher O(2) affinity of Ascaris hemoglobin compared with that for Tyr(B10)/Gln(E7)-containing myoglobin mutants appears to be the result of favorable proximal effects in the Ascaris protein, due to a staggered orientation of His(F8), the lack of a hydrogen bonding lattice between the F4, F7, and F8 residues, and the presence of a large polar Trp(G5) residue in the interior portion of the proximal heme pocket.

The effect of ligand dynamics on heme electronic transition band III in myoglobin

Band III is a near-infrared electronic transition at ~13,000 cm(-1) in heme proteins that has been studied extensively as a marker of protein conformational relaxation after photodissociation of the heme-bound ligand. To examine the influence of the heme pocket structure and ligand dynamics on band III, we have studied carbon monoxide recombination in a variety of myoglobin mutants after photolysis at 3 K using Fourier transform infrared temperature-derivative spectroscopy with monitoring in three spectral ranges, (1) band III, the mid-infrared region of (2) the heme-bound CO, and (3) the photodissociated CO. Here we present data on mutant myoglobins V68F and L29W, which both exhibit pronounced ligand movements at low temperature. From spectral and kinetic analyses in the mid-infrared, a small number of photoproduct populations can be distinguished, differing in their distal heme pocket conformations and/or CO locations. We have decomposed band III into its individual photoproduct contributions. Each photoproduct state exhibits a different "kinetic hole-burning" (KHB) effect, a coupling of the activation enthalpy for rebinding to the position of band III. The analysis reveals that the heme pocket structure and the photodissociated CO markedly affect the band III transition. A strong kinetic hole-burning effect results only when the CO ligand resides in the docking site on top of the heme group. Migration of CO away from the heme group leads to an overall blue shift of band III. Consequently, band III can be used as a sensitive tool to study ligand dynamics after photodissociation in heme proteins.

The heme environment of mouse neuroglobin. Evidence for the presence of two conformations of the heme pocket

Neuroglobin (Ngb) is a newly discovered oxygen-binding heme protein that is primarily expressed in the brain of humans and other vertebrates. To characterize the structure/function relationships of this new heme protein, we have used resonance Raman spectroscopy to determine the structure of the heme environment in Ngb from mice. In the Fe(2+)CO complex, two conformations of the Fe-CO unit are present, one of which arises from an open conformation of the heme pocket in which the CO is not interacting with any nearby residue, and the other arises from a closed conformation where a positively charged residue near the CO group stabilizes the complex. For the Fe(2+)O(2) complex, we detect a single nu(Fe-OO) stretching mode at a frequency similar to that of oxymyoglobins and oxyhemoglobins of vertebrates (571 cm(-1)). Based on the Fe-C-O frequencies of the closed conformation of Ngb, a highly polar distal environment is indicated from which the O(2) off-rate is predicted to be lower than that of Mb. In the absence of exogenous ligands, a heme pocket residue coordinates to the heme iron, forming a six-coordinate complex, thereby predicting a low on-rate for exogenous ligands. These structural properties of the heme pocket of Ngb are discussed with respect to its proposed in vivo oxygen delivery function.

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%.

Roles of Ile209 and Ile210 on the heme pocket structure and regulation of histidine kinase activity of oxygen sensor FixL from Rhizobium meliloti

FixL is a sensor histidine kinase having a heme-containing domain as an O(2) sensing site. In the study presented here, Ile209 and Ile210 located near the heme iron of the heme domain of Rhizobium meliloti FixL (RmFixL) were mutated, and the mutational effects on the regulation of the kinase activity and the heme pocket structure were examined by the autophosphorylation assay and UV-visible absorption and resonance Raman (RR) spectroscopies. The mutation of these residues disrupted the regulation of the kinase activity by the sensor (heme) domain, indicating that Ile209 and Ile210 play important roles in the signal transduction between the heme and the kinase domains. By measurement of the resonance Raman and optical absorption spectra of Ile209 and Ile210 mutants in several oxidation, spin, and ligation states, it was found that both residues are highly flexible, and their side chains sterically interact with the O(2) ligand, when it binds to the heme iron. On the basis of the results, we propose an O(2) sensing mechanism of RmFixL; the kinase activity is regulated via conformational changes of Ile209 and Ile210 induced by the O(2) binding to the sensory center.

Sensory mechanism of oxygen sensor FixL from Rhizobium meliloti: crystallographic, mutagenesis and resonance Raman spectroscopic studies

FixL of Rhizobium meliloti (RmFixL) is a sensor histidine kinase of the two-component system, which regulates the expression of the genes related to nitrogen fixation in the root nodule in response to the O(2) levels. The crystal structure of the sensor domain of FixL (RmFixLH), which contains a heme (Fe-porphyrin) as a sensing site, was determined at 1.4 A resolution. Based on the structural and spectroscopic analyses, we propose the O(2) sensing mechanism that differs from the case proposed in BjFixLH as follows; conformational changes in the F/G loop, which are induced by steric repulsion between the bent-bound O(2) and the Ile209 side-chain, would be transmitted to the histidine kinase domain. Interaction between the iron-bound O(2) and Ile209 was also observed in the resonance Raman spectra of RmFixLH as evidenced by the fact that the Fe-O(2) and Fe-CN stretching frequencies were shifted from 575 to 570 cm(-1) (Fe-O(2)), and 504 to 499 cm(-1), respectively, as the result of the replacement of Ile209 with an Ala residue. In the I209A mutant of RmFixL, the O(2) sensing activity was destroyed, thus confirming our proposed mechanism.

Properties and reactivity of myoglobin reconstituted with chemically modified protohemin complexes

The synthetic complexes protohemin-6(7)-L-arginyl-L-alanine (HM-RA) and protohemin-6(7)-L-histidine methyl ester (HM-H) were prepared by condensation of suitably protected Arg-Ala or His residues with protohemin IX. HM-RA and HM-H were used for reconstitution of apomyoglobin from horse heart, yielding the Mb-RA and Mb-H derivatives, respectively, of the protein. The spectral, binding and catalytic properties of Mb-RA and Mb-H are significantly different from those of Mb. As shown by MM and MD calculations, these differences are determined by some local structural changes around the heme which are generated by increased mobility of a key peptide segment (Phe43-Lys47), containing the residue (Lys45) that in native Mb interacts with one of the porphyrin carboxylate groups. In the reconstituted Mbs this carboxylate group is bound to the Arg-Ala or His residue and is no longer available for electrostatic interaction with Lys45. The mobility of the peptide segment near the active site allows the distal histidine to come to a closer contact with the heme, and in fact Mb-RA and Mb-H exist as an equilibrium between a high-spin form and a major low-spin, six-coordinated form containing a bis-imidazole ligated heme. The two forms are clearly distinguishable in the NMR spectra, that also show that each of them consists of a mixture of the two most stable isomers resulting from cofactor reconstitution, as also anticipated by MM and MD calculations. Exogenous ligands such as cyanide, azide, or hydrogen peroxide can displace the bound distal histidine, but their affinity is reduced. On the other hand, mobilization of the peptide chain around the heme in the reconstituted Mbs increases the accessibility of large donor molecules at the heme periphery, with respect to native Mb, where a rigid backbone limits access to the distal pocket. The increased active site accessibility of Mb-RA and Mb-H facilitates the binding and electron transfer of phenolic substrates in peroxidase-type oxidations catalyzed by the reconstituted proteins in the presence of hydrogen peroxide.