The leghemoglobin proximal heme pocket directs oxygen dissociation and stabilizes bound heme

Iron-histidine bonding in bishistidyl hemoproteins-A local vibrational mode study

We investigated the intrinsic strength of distal and proximal FeN bonds for both ferric and ferrous oxidation states of bishistidyl hemoproteins from bacteria, animals, human, and plants, including two cytoglobins, ten hemoglobins, two myoglobins, six neuroglobins, and six phytoglobins. As a qualified measure of bond strength, we used local vibrational force constants ka (FeN) based on local mode theory developed in our group. All calculations were performed with a hybrid QM/MM ansatz. Starting geometries were taken from available x-ray structures. ka (FeN) values were correlated with FeN bond lengths and covalent bond character. We also investigated the stiffness of the axial NFeN bond angle. Our results highlight that protein effects are sensitively reflected in ka (FeN), allowing one to compare trends in diverse protein groups. Moreover, ka (NFeN) is a perfect tool to monitor changes in the axial heme framework caused by different protein environments as well as different Fe oxidation states.

Best Practices on QM/MM Simulations of Biological Systems

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

A novel single sensor hemoglobin domain from the thermophilic cyanobacteria Thermosynechococcus elongatus BP-1 exhibits higher pH but lower thermal stability compared to globins from mesophilic organisms

Thermosynechococcus elongatus-BP1 belongs to the class of photoautotrophic cyanobacterial organisms. The presence of chlorophyll a, carotenoids, and phycocyanobilin are the characteristics that categorize T. elongatus as a photosynthetic organism. Here, we report the structural and spectroscopic characteristics of novel hemoglobin (Hb) Synel Hb from T.elongatus, synonymous with Thermosynechococcus vestitus BP-1. The X-ray crystal structure (2.15 Å) of Synel Hb suggests the presence of a globin domain with a pre-A helix similar to the sensor domain (S) family of Hbs. The rich hydrophobic core accommodates heme in a penta-coordinated state and readily binds an extraneous ligand(imidazole). The absorption and circular dichroic spectral analysis of Synel Hb reiteratedthat the heme is in Fe^III+ state with a predominantly α-helical structure similar to myoglobin. Synel Hb displays higher resistance to structural perturbations induced via external stresses like pH and guanidium hydrochloride, which is comparable to Synechocystis Hb. However, Synel Hb exhibited lower thermal stability compared to mesophilic hemoglobins. Overall, the data is suggestive of the structural sturdiness of Synel Hb, which probably corroborates its origin in extreme thermophilic conditions. The stable globin provides scope for further investigation and may lead to new insights with scope for engineering stability in hemoglobin-based oxygen carriers.

Soluble guanylyl cyclase: Molecular basis for ligand selectivity and action in vitro and in vivo

Nitric oxide (NO), carbon monoxide (CO), oxygen (O₂), hydrogen sulfide (H₂S) are gaseous molecules that play important roles in the physiology and pathophysiology of eukaryotes. Tissue concentrations of these physiologically relevant gases vary remarkable from nM range for NO to high μM range of O₂. Various hemoproteins play a significant role in sensing and transducing cellular signals encoded by gaseous molecules or in transporting them. Soluble guanylyl cyclase (sGC) is a hemoprotein that plays vital roles in a wide range of physiological functions and combines the functions of gaseous sensor and signal transducer. sGC uniquely evolved to sense low non-toxic levels of NO and respond to elevated NO levels by increasing its catalytic ability to generate the secondary signaling messenger cyclic guanosine monophosphate (cGMP). This review discusses sGC's gaseous ligand selectivity and the molecular basis for sGC function as high-affinity and selectivity NO receptor. The effects of other gaseous molecules and small molecules of cellular origin on sGC's function are also discussed.

Stability and Folding of the Unusually Stable Hemoglobin from Synechocystis is Subtly Optimized and Dependent on the Key Heme Pocket Residues

AIMS

The aim of our study is to understand the biophysical traits that govern the stability and folding of Synechocystis hemoglobin, a unique cyanobacterial globin that displays unusual traits not observed in any of the other globins discovered so far.

BACKGROUND

For the past few decades, classical hemoglobins such as vertebrate hemoglobin and myoglobin have been extensively studied to unravel the stability and folding mechanisms of hemoglobins. However, the expanding wealth of hemoglobins identified in all life forms with novel properties, like heme coordination chemistry and globin fold, have added complexity and challenges to the understanding of hemoglobin stability, which has not been adequately addressed. Here, we explored the unique truncated and hexacoordinate hemoglobin from the freshwater cyanobacterium Synechocystis sp. PCC 6803 known as "Synechocystis hemoglobin (SynHb)". The "three histidines" linkages to heme are novel to this cyanobacterial hemoglobin.

OBJECTIVE

Mutational studies were employed to decipher the residues within the heme pocket that dictate the stability and folding of SynHb.

METHODS

Site-directed mutants of SynHb were generated and analyzed using a repertoire of spectroscopic and calorimetric tools.

RESULTS

The results revealed that the heme was stably associated to the protein under all denaturing conditions with His117 playing the anchoring role. The studies also highlighted the possibility of existence of a "molten globule" like intermediate at acidic pH in this exceptionally thermostable globin. His117 and other key residues in the heme pocket play an indispensable role in imparting significant polypeptide stability.

CONCLUSION

Synechocystis hemoglobin presents an important model system for investigations of protein folding and stability in general. The heme pocket residues influenced the folding and stability of SynHb in a very subtle and specific manner and may have been optimized to make this Hb the most stable known as of date. Other: The knowledge gained hereby about the influence of heme pocket amino acid side chains on stability and expression is currently being utilized to improve the stability of recombinant human Hbs for efficient use as oxygen delivery vehicles.

A new paradigm for gaseous ligand selectivity of hemoproteins highlighted by soluble guanylate cyclase

Nitric oxide (NO), carbon monoxide (CO), and oxygen (O₂) are important physiological messengers whose concentrations vary in a remarkable range, [NO] typically from nM to several μM while [O₂] reaching to hundreds of μM. One of the machineries evolved in living organisms for gas sensing is sensor hemoproteins whose conformational change upon gas binding triggers downstream response cascades. The recently proposed "sliding scale rule" hypothesis provides a general interpretation for gaseous ligand selectivity of hemoproteins, identifying five factors that govern gaseous ligand selectivity. Hemoproteins have intrinsic selectivity for the three gases due to a neutral proximal histidine ligand while proximal strain of heme and distal steric hindrance indiscriminately adjust the affinity of these three gases for heme. On the other hand, multiple-step NO binding and distal hydrogen bond donor(s) specifically enhance affinity for NO and O₂, respectively. The "sliding scale rule" hypothesis provides clear interpretation for dramatic selectivity for NO over O₂ in soluble guanylate cyclase (sGC) which is an important example of sensor hemoproteins and plays vital roles in a wide range of physiological functions. The "sliding scale rule" hypothesis has so far been validated by all experimental data and it may guide future designs for heme-based gas sensors.

Comparative analysis of the heme iron electronic structure and stereochemistry in tetrameric rabbit hemoglobin and monomeric soybean leghemoglobin a using Mössbauer spectroscopy with a high velocity resolution

A comparative study of tetrameric rabbit hemoglobin and monomeric soybean leghemoglobin a in the oxy- and deoxy-forms was carried out using 57Fe Mössbauer spectroscopy with a high velocity resolution in order to analyze the heme iron electronic structure and stereochemistry in relation to the Mössbauer hyperfine parameters. The Mössbauer spectra of tetrameric rabbit hemoglobin in both forms were fitted using two quadrupole doublets related to the 57Fe in ɑ- and β-subunits. In contrast, the Mössbauer spectra of monomeric soybean leghemoglobin a were fitted using: (i) two quadrupole doublets for the oxy-form related to two conformational states of the distal His E7 imidazole ring and different hydrogen bonding of oxygen molecule in the oxy-form and (ii) using three quadrupole doublets for deoxy-form related to three conformational states of the proximal His F8 imidazole ring. Small variations of Mössbauer hyperfine parameters related to small differences in the heme iron electronic structure and stereochemistry in tetrameric rabbit hemoglobin and monomeric soybean leghemoglobin a are discussed.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

Phe28B10 Induces Channel-Forming Cytotoxic Amyloid Fibrillation in Human Neuroglobin, the Brain-Specific Hemoglobin

Since its discovery, neuroglobin (Ngb), a neuron-specific oxygen binding hemoglobin, distinct from the classical myoglobin and blood hemoglobin, has attracted attention as an endogenous neuroprotectant. Recent reports suggest that Ngb protects neurons from brain stroke, ischemic stress-induced degeneration, and other brain disorders. Proteins with a specific role in neuroprotection are often associated with neurodegeneration, as well, depending on the cellular environment or specific cellular triggers that tilt the balance one way or the other. This investigation explored the potential role of Ngb in amyloid fibril-related neuronal disorder. Ngb was capable of amyloid formation in vitro at neutral pH and ambient temperature, in both apo and holo forms, albeit at a slower rate in the holo form, unlike other hemoglobins that exhibit such behavior exclusively in the apo states. Elevated temperature enhanced the rate of fibril formation significantly. The B-helix, which is known to play a major role in Ngb ligand binding kinetics, was found to be amyloidogenic with the Phe28^B10 amino acid side chain as the key inducer of fibrillation. The Ngb amyloid fibril was also significantly cytotoxic to neuroblastoma cell lines, compared to those obtained from reference hemoglobins. The Ngb fibril probably promoted toxicity by inducing channel formation in the cell membrane, as investigated here using synthetic lipid bilayer membranes and the propidium iodide uptake assay. These findings imply that Ngb plays a role in neurodegenerative disorders in vivo, for which there seems to be indirect evidence by association. Ngb thus presents a novel prospect for understanding amyloid-related brain disorders beyond the limited set of proteins currently investigated for such diseases.

Structural and Functional Significance of the N- and C-Terminal Appendages in Arabidopsis Truncated Hemoglobin

Plant hemoglobins constitute three distinct groups: symbiotic, nonsymbiotic, and truncated hemoglobins. Structural investigation of symbiotic and nonsymbiotic (class I) hemoglobins revealed the presence of a vertebrate-like 3/3 globin fold in these proteins. In contrast, plant truncated hemoglobins are similar to bacterial truncated hemoglobins with a putative 2/2 α-helical globin fold. While multiple structures have been reported for plant hemoglobins of the first two categories, for plant truncated globins only one structure has been reported of late. Here, we report yet another crystal structure of the truncated hemoglobin from Arabidopsis thaliana (AHb3) with two water molecules in the heme pocket, of which one is distinctly coordinated to the heme iron, unlike the only available crystal structure of AHb3 with a hydroxyl ligand. AHb3 was monomeric in its crystallographic asymmetric unit; however, dimer was evident in the crystallographic symmetry, and the globin indeed existed as a stable dimer in solution. The tertiary structure of the protein exhibited a bacterial-like 2/2 α-helical globin fold with an additional N-terminal α-helical extension and disordered C-termini. To address the role of these extended termini in AHb3, which is yet unknown, N- and C-terminal deletion mutants were created and characterized and molecular dynamics simulations performed. The C-terminal deletion had an insignificant effect on most properties but perturbed the dimeric equilibrium of AHb3 and significantly influenced azide binding kinetics in the ferric state. These results along with the disordered nature of the C-terminus indicated its putative role in intramolecular or intermolecular interactions probably regulating protein-ligand and protein-protein interactions. While the N-terminal deletion did not change the overall globin fold, stability, or ligand binding kinetics, it seemed to have influenced coordination at the heme iron, the hydration status of the active site, and the quaternary structure of AHb3. Evidence indicated that the N-terminus is the predominant factor regulating the quaternary interaction appropriate to physiological requirements, dynamics of the side chains in the heme pocket, and tunnel organization in the protein matrix.

Molecular dynamics simulations indicate that tyrosineB10 limits motions of distal histidine to regulate CO binding in soybean leghemoglobin

Myoglobin (Mb) uses strong electrostatic interaction in its distal heme pocket to regulate ligand binding. The mechanism of regulation of ligand binding in soybean leghemoglobin a (Lba) has been enigmatic and more so due to the absence of gaseous ligand bound atomic resolution three-dimensional structure of the plant globin. While the 20-fold higher oxygen affinity of Lba compared with Mb is required for its dual physiological function, the mechanism by which this high affinity is achieved is only emerging. Extensive mutational analysis combined with kinetic and CO-FT-IR spectroscopic investigation led to the hypothesis that Lba depended on weakened electrostatic interaction between distal HisE7 and bound ligand achieved by invoking B10Tyr, which itself hydrogen bonds with HisE7 thus restricting it in a single conformation detrimental to Mb-like strong electrostatic interaction. Such theory has been re-assessed here using CO-Lba in silico model and molecular dynamics simulation. The investigation supports the presence of at least two major conformations of HisE7 in Lba brought about by imidazole ring flip, one of which makes hydrogen bonds effectively with B10Tyr affecting the former's ability to stabilize bound ligand, while the other does not. However, HisE7 in Lba has limited conformational freedom unlike high frequency of imidazole ring flips observed in Mb and in TyrB10Leu mutant of Lba. Thus, it appears that TyrB10 limits the conformational freedom of distal His in Lba, tuning down ligand dissociation rate constant by reducing the strength of hydrogen bonding to bound ligand, which the freedom of distal His of Mb allows.

Significantly enhanced heme retention ability of myoglobin engineered to mimic the third covalent linkage by nonaxial histidine to heme (vinyl) in synechocystis hemoglobin

Heme proteins, which reversibly bind oxygen and display a particular fold originally identified in myoglobin (Mb), characterize the "hemoglobin (Hb) superfamily." The long known and widely investigated Hb superfamily, however, has been enriched by the discovery and investigation of new classes and members. Truncated Hbs typify such novel classes and exhibit a distinct two-on-two α-helical fold. The truncated Hb from the freshwater cyanobacterium Synechocystis exhibits hexacoordinate heme chemistry and bears an unusual covalent bond between the nonaxial His(117) and a heme porphyrin 2-vinyl atom, which remains tightly associated with the globin unlike any other. It seems to be the most stable Hb known to date, and His(117) is the dominant force holding the heme. Mutations of amino acid residues in the vicinity did not influence this covalent linkage. Introduction of a nonaxial His into sperm whale Mb at the topologically equivalent position and in close proximity to vinyl group significantly increased the heme stability of this prototype globin. Reversed phase chromatography, electrospray ionization-MS, and MALDI-TOF analyses confirmed the presence of covalent linkage in Mb I107H. The Mb mutant with the engineered covalent linkage was stable to denaturants and exhibited ligand binding and auto-oxidation rates similar to the wild type protein. This indeed is a novel finding and provides a new perspective to the evolution of Hbs. The successful attempt at engineering heme stability holds promise for the production of stable Hb-based blood substitute.

Influence of heme environment structure on dioxygen affinity for the dual function Amphitrite ornata hemoglobin/dehaloperoxidase. Insights into the evolutional structure-function adaptations

Sea worm, Amphitrite ornata, has evolved its globin (an O(2) carrier) also to serves as a dehaloperoxidase (DHP) to detoxify haloaromatic pollutants generated by competing species. A previous mutagenesis study by our groups on both DHP and sperm whale myoglobin (SW Mb) revealed some structural factors that influence the dehaloperoxidase activities (significantly lower for Mb) of both proteins. Using an isocyanide/O(2) partition constant measurement method in this study, we have examined the effects of these structural factors on the O(2) equilibrium constants (KO2) of DHP, SW Mb, and their mutants. A clear trend of decreasing O(2) affinity and increasing catalytic activity along with the increase in the distal His N(ε)-heme iron distance is observed. An H93K/T95H Mb double mutant mimicking the DHP proximal His positioning exhibited markedly enhanced O(2) affinity, confirming the essential effect of proximal His rotation on the globin function of DHP. For DHP, the L100F, T56G and M86E variants showed the effects of distal volume, distal His flexibility and proximal electronic push, respectively, on the O(2) affinity. This study provides insights into how DHP has evolved its heme environment to gain significantly enhanced peroxidase capability without compromising its primary function as an O(2) carrier.

Small ligand-globin interactions: reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

X-ray crystallographic structural characteristics of Arabidopsis hemoglobin I and their functional implications

Genome of the model dicot flowering plant, Arabidopsis thaliana, a popular tool for understanding molecular biology of plant physiology, encodes all three classes of plant hemoglobins that differ in their sequence, ligand binding and spectral properties. As such these globins are of considerable attention. Crystal structures of few members of plant class I nonsymbiotic hemoglobin have been described earlier. Here we report the crystal structure of Arabidopsis class I hemoglobin (AHb1) to 2.2Ǻ and compare its key features with the structures of similar nonsymbiotic hemoglobin from other species. Crystal structure of AHb1 is homologous to the related members with similar globin fold and heme pocket architecture. The structure is homodimeric in the asymmetric unit with both distal and proximal histidines coordinating to the heme iron atom. Residues lining the dimeric interface are also conserved in AHb1 with the exception of additional electrostatic interaction between H112 and E113 of each subunit and that involving Y119 through two water molecules. In addition, differences in heme pocket non-covalent interactions, a novel Ser residue at F7 position, Xe binding site variability, internal cavity topology differences, CD loop conformation and stability and other such properties might explain kinetic variability in AHb1. Detailed cavity analysis of AHb1 showed the presence of a novel long tunnel connecting the distal pockets of both the monomers. Presence of such tunnel, along with conformational heterogeneity observed in the two chains, might suggest cooperative ligand binding and support its role in NO scavenging. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

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

SIGNIFICANCE

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

RECENT ADVANCES

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

CRITICAL ISSUES

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

FUTURE DIRECTIONS

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

A singular value decomposition approach for kinetic analysis of reactions of HNO with myoglobin

The reactions of several horse heart myoglobin species with nitrosyl hydride, HNO, derived from Angeli's salt (AS) and Piloty's acid (PA) have been followed by UV-visible, (1)H NMR and EPR spectroscopies. Spectral analysis of myoglobin-derived speciation during the reactions was obtained by using singular value decomposition methods combined with a global analysis to obtain the rate constants of complex sequential reactions. The analysis also provided spectra for the derived absorbers, which allowed self-consistent calibration to the spectra of known myoglobin species. Using this method, the determined rate for trapping of HNO by metmyoglobin, which produces NO-myoglobin, is found to be 2.7 × 10(5)M(-1)s(-1) at pH7.0 and 1.1 × 10(5)M(-1)s(-1) at pH9.4. The reaction of deoxymyoglobin with HNO generates the adduct HNO-myoglobin directly, but is followed by a secondary reaction of that product with HNO yielding NO-myoglobin; the determined bimolecular rate constants for these reactions are 3.7 × 10(5)M(-1)s(-1) and 1.67 × 10(4)M(-1)s(-1) respectively, and are independent of pH. The derived spectrum for HNO-myoglobin is characterized by a Soret absorbance maximum at 423 nm with an extinction coefficient of 1.66 × 10(5)M(-1)cm(-1). The rate constant for unimolecular loss of HNO from HNO-myoglobin was determined by competitive trapping with CO at 8.9 × 10(-5)s(-1), which gives the thermodynamic binding affinity of HNO to deoxymyoglobin as 4.2 × 10(9)M(-1). These results suggest that the formation of HNO-ferrous heme protein adducts represents an important consideration in the biological action of HNO-releasing drugs.

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

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

CO binding and ligand discrimination in human myeloperoxidase

Despite the fact that ferrous myeloperoxidase (MPO) can bind both O(2) and NO, its ability to bind CO has been questioned. UV/visible spectroscopy was used to confirm that CO induces small spectral shifts in ferrous MPO, and Fourier transform infrared difference spectroscopy showed definitively that these arose from formation of a heme ferrous-CO compound. Recombination rates after CO photolysis were monitored at 618 and 645 nm as a function of CO concentration and pH. At pH 6.3, k(on) and k(off) were 0.14 mM(-1) x s(-1) and 0.23 s(-1), respectively, yielding an unusually high K(D) of 1.6 mM. This affinity of MPO for CO is 10 times weaker than its affinity for O(2). The observed rate constant for CO binding increased with increasing pH and was governed by a single protonatable group with a pK(a) of 7.8. Fourier transform infrared spectroscopy revealed two different conformations of bound CO with frequencies at 1927 and 1942 cm(-1). Their recombination rate constants were identical, indicative of two forms of bound CO that are in rapid thermal equilibrium rather than two distinct protein populations with different binding sites. The ratio of bound states was pH-dependent (pK(a) approximately 7.4) with the 1927 cm(-1) form favored at high pH. Structural factors that account for the ligand-binding properties of MPO are identified by comparisons with published data on a range of other ligand-binding heme proteins, and support is given to the recent suggestion that the proximal His336 in MPO is in a true imidazolate state.

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

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

The structure and function of plant hemoglobins

Plants, like humans, contain hemoglobin. Three distinct types of hemoglobin exist in plants: symbiotic, non-symbiotic, and truncated hemoglobins. Crystal structures and other structural and biophysical techniques have revealed important knowledge about ligand binding and conformational stabilization in all three types. In symbiotic hemoglobins (leghemoglobins), ligand binding regulatory mechanisms have been shown to differ dramatically from myoglobin and red blood cell hemoglobin. In the non-symbiotic hemoglobins found in all plants, crystal structures and vibrational spectroscopy have revealed the nature of the structural transition between the hexacoordinate and ligand-bound states. In truncated hemoglobins, the abbreviated globin is porous, providing tunnels that may assist in ligand binding, and the bound ligand is stabilized by more than one distal pocket residue. Research has implicated these plant hemoglobins in a number of possible functions differing among hemoglobin types, and possibly between plant species.

Influence of the protein matrix on intramolecular histidine ligation in ferric and ferrous hexacoordinate hemoglobins

Present in most organisms, hexacoordinate hemoglobins (hxHbs) are proteins that have evolved the capacity for reversible bis-histidyl heme coordination. The heme prosthetic group enables diverse protein functionality, such as electron transfer, redox reactions, ligand transport, and enzymatic catalysis. The reactivity of heme is greatly effected by the coordination and noncovalent chemical environment imposed by its connate protein. Of considerable interest is how the hxHb globin fold achieves reversible intramolecular coordination while still enabling high-affinity binding of oxygen, nitric oxide, and other small ligands. Here we explore this question by examining the role of the protein matrix on coordination behavior in a group of hxHbs from animals, plants, and bacteria, including human neuroglobin and cytoglobin, a nonsymbiotic hemoglobin from rice, and a truncated hemoglobin from the cyanobacterium Synechocystis. This is done with a set of experiments measuring the reduction potentials of each wild-type hxHb and its corresponding mutant protein where the reversibly bound histidine (the distal His) has been replaced with a noncoordinating side chain. These reduction potentials, coupled with studies of the mutant proteins saturated with exogenous imidazole, enable us to assess the effects of the protein matrices on histidine coordination. Our results show significant variation among the hxHbs, demonstrating flexibility in the globin moiety's ability to regulate reversible coordination. This regulation is particularly evident in the plant nonsymbiotic hemoglobins, where ferric state histidine coordination affinity is substantially lowered by the protein matrix.

Proximal influences in two-on-two globins: effect of the Ala69Ser replacement on Synechocystis sp. PCC 6803 hemoglobin

The cyanobacterium Synechocystis sp. PCC 6803 (S6803) expresses a two-on-two globin in which His46 (distal side) and His70 (proximal) function as heme iron axial ligands. His46 can be displaced by O2, CO, and CN-, among others, whereas His70 is not labile under native conditions. The residue preceding the proximal histidine has been implicated in controlling globin axial ligand reactivity; the details of the mechanism, however, are not well understood, and little information exists for bis-histidyl hexacoordinate proteins. In many vertebrate hemoglobins and in the Synechocystis protein, the position is occupied by an alanine, whereas, in myoglobins, it is a serine involved in an intricate hydrogen-bond network. We examined the role of Ala69 in S6803 hemoglobin through the effects of an Ala --> Ser replacement. The substitution resulted in minor structural perturbations, but the response of the holoprotein to temperature-, urea-, and acid-induced denaturation was measurably affected. Enhanced three-state behavior was manifested in the decoupling of heme binding and secondary-structure formation. Urea-gradient gel experiments revealed that the stability of the apoprotein was unchanged by the replacement and that a slight alteration of the folding kinetics occurred in the holoproteins. Cyanide-binding experiments were performed to assess trans effects. The apparent rate constant for association decreased 2-fold upon Ala69Ser replacement. This deceleration was attributed to a change in the lifetime of a state containing a decoordinated His46. The results demonstrated that, as in vertebrate globins and leghemoglobin, proximal influences operate to determine fundamental dynamic and thermodynamic properties of the protein.

Heme Protein Oxygen Affinity Regulation Exerted by Proximal Effects

Heme proteins are found in all living organisms and are capable of performing a wide variety of tasks, requiring in many cases the binding of diatomic ligands, namely, O(2), CO, and/or NO. Therefore, subtle regulation of these diatomic ligands' affinity is one of the key issues for determining a heme protein's function. This regulation is achieved through direct H-bond interactions between the bound ligand and the protein, and by subtle tuning of the intrinsic heme group reactivity. In this work, we present an investigation of the proximal regulation of oxygen affinity in Fe(II) histidine coordinated heme proteins by means of computer simulation. Density functional theory calculations on heme model systems are used to analyze three proximal effects: charge donation, rotational position, and distance to the heme porphyrin plane of the proximal histidine. In addition, hybrid quantum-classical (QM-MM) calculations were performed in two representative proteins: myoglobin and leghemoglobin. Our results show that all three effects are capable of tuning the Fe-O(2) bond strength in a cooperative way, consistently with the experimental data on oxygen affinity. The proximal effects described herein could operate in a large variety of O(2)-binding heme proteins-in combination with distal effects-and are essential to understand the factors determining a heme protein's O(2) affinity.

Dioxygen affinity in heme proteins investigated by computer simulation

We present an investigation of the molecular basis of the modulation of oxygen affinity in heme proteins using computer simulation. QM-MM calculations are applied to explore distal and proximal effects on O(2) binding to the heme, while classical molecular dynamics simulations are employed to investigate ligand migration across the polypeptide to the active site. Trends in binding energies and in the kinetic constants are illustrated through a number of selected examples highlighting the virtues and the limitations of the applied methodologies. These examples cover a wide range of O(2)-affinities, and include: the truncated-N and truncated-O hemoglobins from Mycobacterium tuberculosis, the mammalian muscular O(2) storage protein: myoglobin, the hemoglobin from the parasitic nematode Ascaris lumbricoides, the oxygen transporter in the root of leguminous plants: leghemoglobin, the Cerebratulus lacteus nerve tissue hemoglobin, and the Alcaligenes xyloxidans cytochrome c'.

Multiple active site conformers in the carbon monoxide complexes of trematode hemoglobins

Sequence alignment of hemoglobins of the trematodes Paramphistomum epiclitum and Gastrothylax crumenifer with myoglobin suggests the presence of an unusual active site structure in which two tyrosine residues occupy the E7 and B10 helical positions. In the crystal structure of P. epiclitum hemoglobin, such an E7-B10 tyrosine pair at the putative helical positions has been observed, although the E7 Tyr is displaced toward CD region of the polypeptide. Resonance Raman data on both P. epiclitum and G. crumenifer hemoglobins show that interactions of heme-bound ligands with neighboring amino acid residues are unusual. Multiple conformers in the CO complex, termed the C, O, and N conformers, are observed. The conformers are separated by a large difference (approximately 60 cm(-1)) in the frequencies of their Fe-CO stretching modes. In the C conformer the Fe-CO stretching frequency is very high, 539 and 535 cm(-1), for the P. epiclitum and G. crumenifer hemoglobins, respectively. The Fe-CO stretching of the N conformer appears at an unusually low frequency, 479 and 476 cm(-1), respectively, for the two globins. A population of an O conformer is seen in both hemoglobins, at 496 and 492 cm(-1), respectively. The C conformer is stabilized by a strong polar interaction of the CO with the distal B10 tyrosine residue. The O conformer is similar to the ones typically seen in mutant myoglobins in which there are no strong interactions between the CO and residues in the distal pocket. The N conformer possesses an unusual configuration in which a negatively charged group, assigned as the oxygen atom of the B10 Tyr side chain, interacts with the CO. In this conformer, the B10 Tyr assumes an alternative conformation consistent with one of the conformers seen the crystal structure. Implications of the multiple configurations on the ligand kinetics are discussed.

Bis-histidyl hexacoordination in hemoglobins facilitates heme reduction kinetics

Hexacoordinate hemoglobins are a class of proteins that exhibit reversible bis-histidyl coordination of the heme iron while retaining the ability to bind exogenous ligands. One hypothesis for their physiological function is that they scavenge nitric oxide, a reaction that oxidizes the protein and requires reduction of the heme iron to continue. Reduction kinetics of hexacoordinate hemoglobins, including human neuroglobin and cytoglobin, and those from Synechocystis and rice, are compared to myoglobin, soybean leghemoglobin, and several relevant mutant proteins. In all cases, bis-histidyl coordination greatly increases the rate of reduction by sodium dithionite when compared to pentacoordinate hemoglobins. In myoglobin and leghemoglobin, reduction is limited by the rate constant for electron transfer, whereas in the hexacoordinate hemoglobins reduction is limited only by bimolecular binding of the reductant. These results can be explained by differences in the reorganization energy for reduction between hexacoordinate and pentacoordinate hemoglobins.

Ligand binding properties of myoglobin reconstituted with iron porphycene: unusual O2 binding selectivity against CO binding

Sperm whale myoglobin, an oxygen storage hemoprotein, was successfully reconstituted with the iron porphycene having two propionates, 2,7-diethyl-3,6,12,17-tetramethyl-13,16-bis(carboxyethyl)porphycenatoiron. The physicochemical properties and ligand bindings of the reconstituted myoglobin were investigated. The ferric reconstituted myoglobin shows the remarkable stability against acid denaturation and only a low-spin characteristic in its EPR spectrum. The Fe(III)/Fe(II) redox potential (-190 mV vs NHE) determined by the spectroelectrochemical measurements was much lower than that of the wild-type. These results can be attributed to the strong coordination of His93 to the porphycene iron, which is induced by the nature of the porphycene ring symmetry. The O2 affinity of the ferrous reconstituted myoglobin is 2600-fold higher than that of the wild-type, mainly due to the decrease in the O2 dissociation rate, whereas the CO affinity is not so significantly enhanced. As a result, the O2 affinity of the reconstituted myoglobin exceeds its CO affinity (M' = K(CO)/K(O2) < 1). The ligand binding studies on H64A mutants support the fact that the slow O2 dissociation of the reconstituted myoglobin is primarily caused by the stabilization of the Fe-O2 sigma-bonding. The IR spectra for the carbon monoxide (CO) complex of the reconstituted myoglobin suggest several structural and/or electrostatic conformations of the Fe-C-O bond, but this is not directly correlated with the CO dissociation rate. The high O2 affinity and the unique characteristics of the myoglobin with the iron porphycene indicate that reconstitution with a synthesized heme is a useful method not only to understand the physiological function of myoglobin but also to create a tailor-made function on the protein.

The complex of apomyoglobin with the fluorescent dye coumarin 153

Understanding a protein's dielectric response requires both a theoretical model and a well-defined experimental system. The former has already been proposed by Song (J. Chem. Phys. 116, 9359 [2002]). We suggest that the latter is provided by the complex of coumarin 153 (C153) with apomyoglobin (ApoMb). C153 has been exhaustively studied and has proven to be an excellent probe of the solvation dynamics of polar solvents. Myoglobin is one of the most thoroughly studied proteins. Myoglobins from a wide range of species have been subject to X-ray structural analysis and site-directed mutagenesis. Here, we demonstrate the existence of a robust C153-apomyglobin system by means of molecular dynamics simulations, equilibrium binding studies using a Job's plot and capillary electrophoresis, circular dichroism and time-resolved fluorescence. The reorganization energy of C153 bound to ApoMb is compared with that of C153 in bulk solvent using the method of Jordanides et al. (J. Phys. Chem. B 103, 7995 [1999]).

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

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

The Crystal Structure of Synechocystis Hemoglobin with a Covalent Heme Linkage

The x-ray crystal structure of Synechocystis hemoglobin has been solved to a resolution of 1.8 A. The conformation of this structure is surprisingly different from that of the previously reported solution structure, probably due in part to a covalent linkage between the heme 2-vinyl and His117 that is present in the crystal structure but not in the structure solved by NMR. Synechocystis hemoglobin is a hexacoordinate hemoglobin in which the heme iron is coordinated by both the proximal and distal histidines. It is also a member of the "truncated hemoglobin" family that is much shorter in primary structure than vertebrate and plant hemoglobins. In contrast to other truncated hemoglobins, the crystal structure of Synechocystis hemoglobin displays no "ligand tunnel" and shows that several important amino acid side chains extrude into the solvent instead of residing inside the heme pocket. The stereochemistry of hexacoordination is compared with other hexacoordinate hemoglobins and cytochromes in an effort to illuminate factors contributing to ligand affinity in hexacoordinate hemoglobins.

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.

Plants, humans and hemoglobins

New developments have forced a re-evaluation of our understanding of the structure and function of hemoglobins. Leghemoglobins regulate oxygen affinity through a mechanism different from that of myoglobin using a novel combination of heme pocket amino acids that lower the oxygen affinity. The hexacoordinate hemoglobins are characterized by intramolecular coordination of the ligand binding site at the heme iron, and were first identified in plants as the 'non-symbiotic plant hemoglobins'. They are now known to be present in animals and bacteria. Many of these proteins are upregulated in both plants and animals during hypoxia or similar stresses. Therefore, there might be a common physiological function for hexacoordinate hemoglobins in plants and animals.

Direct measurement of equilibrium constants for high-affinity hemoglobins

The biological functions of heme proteins are linked to their rate and affinity constants for ligand binding. Kinetic experiments are commonly used to measure equilibrium constants for traditional hemoglobins comprised of pentacoordinate ligand binding sites and simple bimolecular reaction schemes. However, kinetic methods do not always yield reliable equilibrium constants with more complex hemoglobins for which reaction mechanisms are not clearly understood. Furthermore, even where reaction mechanisms are clearly understood, it is very difficult to directly measure equilibrium constants for oxygen and carbon monoxide binding to high-affinity (K(D) << 1 micro M) hemoglobins. This work presents a method for direct measurement of equilibrium constants for high-affinity hemoglobins that utilizes a competition for ligands between the "target" protein and an array of "scavenger" hemoglobins with known affinities. This method is described for oxygen and carbon monoxide binding to two hexacoordinate hemoglobins: rice nonsymbiotic hemoglobin and Synechocystis hemoglobin. Our results demonstrate that although these proteins have different mechanisms for ligand binding, their affinities for oxygen and carbon monoxide are similar. Their large affinity constants for oxygen, 285 and approximately 100 micro M(-1) respectively, indicate that they are not capable of facilitating oxygen transport.

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.