Thermoglobin, oxygen-avid hemoglobin in a bacterial hyperthermophile

The hemoglobin family of proteins, ubiquitous in all domains of life, evolved from an ancestral protein of primordial function to extant hemoglobins that perform a myriad of functions with diverged biochemical properties. Study of homologs in bacterial hyperthermophiles may shed light on both mechanisms of adaptation to extreme conditions and the nature of the ancestral protein. A hemoglobin was identified in Aquifex aeolicus, cloned, recombinantly expressed, purified, and characterized. This hemoglobin is monomeric, resistant to thermal and chemical denaturation, pentacoordinate in the ferrous deoxygenated state, and oxygen-avid. The oxygen equilibrium dissociation constant is approximately 1 nm at room temperature, due in part to a hydrogen bond between the bound ligand and a tyrosine residue in the distal pocket. These biochemical properties of A. aeolicus thermoglobin, AaTgb, may have been shared by the ancestral hemoglobin, thus suggesting possible primordial functions and providing a starting point for consequent evolution of the hemoglobin family.     

Proton nuclear magnetic resonance studies of hemoglobin Malm?: implications of mutations at homologous positions of the alpha and beta chains.

The abnormal human hemoglobin Malm? (beta97FG4 His leads to Gln) has been studied and its properties are compared with those of normal adult hemoglobin A. The data presented here show that the ring-current shifted proton resonances of both HbCO and HbO2 Malm? are very different from the corresponding forms of Hb A. The hyperfine shifted proton resonances of deoxy-Hb Malm? do not differ drastically from those of deoxy-Hb A. This result, together with the finding that the exchangeable proton resonances of the deoxy form of the two hemoglobins are similar, suggests that unliganded Hb Malm? can assume a deoxy-like quaternary structure both in the absence and presence of organic phosphates We have also compared the properties of Hb Malm? with those of Hb Chesapeake (alpha92FG4 Arg leads to Leu). This allows us to study the properties of two abnormal human hemoglobins with mutations at homologous positions of the alpha and beta chains in the three-dimenstional structure of the hemoglobin molecule. Our present results suggest that the mutaion at betaFG4 has its greatest effect on the teritiary structure of the heme pocket of the liganded forms of the hemoglobin while the mutation at alphaFG4 alters the deoxy structure of the hemoglogin molecule but does not alter the teriary structure of the heme pockets of the liganded form of the hemoglobin molecule. Both hemoglobins undergo a transition from the deoxy (T) to the oxy (R) quaternary structure upon ligation. The abnormally high oxygen affinities and low cooperativities of these two hemoglobins must therefore be due to either the structural differences which we have observed and/or to an altered transition between the T and R structures.     

A membrane-bound hemoglobin from gills of the green shore crab Carcinus maenas

Most hemoglobins serve for the transport or storage of O(2). Although hemoglobins are widespread in "entomostracan" Crustacea, malacostracans harbor the copper-containing hemocyanin in their hemolymph. Usually, only one type of respiratory protein occurs within a single species. Here, we report the identification of a hemoglobin of the shore crab Carcinus maenas (Malacostraca, Brachyura). In contrast to the dodecameric hemocyanin of this species, C. maenas hemoglobin does not reside in the hemolymph but is restricted to the gills. Immunofluorescence studies and cell fractioning showed that C. maenas hemoglobin resides in the membrane of the chief cells of the gill. To the best of our knowledge, this is the first time that a membrane-bound hemoglobin has been identified in eukaryotes. Bioinformatic evaluation suggests that C. maenas hemoglobin is anchored in the membrane by N-myristoylation. Recombinant C. maenas hemoglobin has a hexacoordinate binding scheme at the Fe(2+) and an oxygen affinity of P(50) = 0.5 Torr. A rapid autoxidation rate precludes a function as oxygen carrier. We rather speculate that, analogous to prokaryotic membrane-globins, C. maenas hemoglobin carries out enzymatic functions to protect the lipids in cell membrane from reactive oxygen species. Sequence comparisons and phylogenetic studies suggested that the ancestral arthropod hemoglobin was most likely an N-myristoylated protein that did not have an O(2) supply function. True respiratory hemoglobins of arthropods, however, evolved independently in chironomid midges and branchiopod crustaceans.     

Globin gene family evolution and functional diversification in annelids

Globins are the most common type of oxygen-binding protein in annelids. In this paper, we show that circulating intracellular globin (Alvinella pompejana and Glycera dibranchiata), noncirculating intracellular globin (Arenicola marina myoglobin) and extracellular globin from various annelids share a similar gene structure, with two conserved introns at canonical positions B12.2 and G7.0. Despite sequence divergence between intracellular and extracellular globins, these data strongly suggest that these three globin types are derived from a common ancestral globin-like gene and evolved by duplication events leading to diversification of globin types and derived functions. A phylogenetic analysis shows a distinct evolutionary history of annelid extracellular hemoglobins with respect to intracellular annelid hemoglobins and mollusc and arthropod extracellular hemoglobins. In addition, dehaloperoxidase (DHP) from the annelid, Amphitrite ornata, surprisingly exhibits close phylogenetic relationships to some annelid intracellular globins. We have characterized the gene structure of A. ornata DHP to confirm assumptions about its homology with globins. It appears that it has the same intron position as in globin genes, suggesting a common ancestry with globins. In A. ornata, DHP may be a derived globin with an unusual enzymatic function.     

Interfacial and distal-heme pocket mutations exhibit additive effects on the structure and function of hemoglobin

Protein engineering strategies seek to develop a hemoglobin-based oxygen carrier with optimized functional properties, including (i) an appropriate O 2 affinity, (ii) high cooperativity, (iii) limited NO reactivity, and (iv) a diminished rate of auto-oxidation. The mutations alphaL29F, alphaL29W, alphaV96W and betaN108K individually impart some of these traits and in combinations produce hemoglobin molecules with interesting ligand-binding and allosteric properties. Studies of the ligand-binding properties and solution structures of single and multiple mutants have been performed. The aromatic side chains placed in the distal-heme pocket environment affect the intrinsic ligand-binding properties of the mutated subunit itself, beyond what can be explained by allostery, and these changes are accompanied by local structural perturbations. In contrast, hemoglobins with mutations in the alpha 1beta 1 and alpha 1beta 2 interfaces display functional properties of both "R"- and "T"-state tetramers because the equilibrium between quaternary states is altered. These mutations are accompanied by global structural perturbations, suggesting an indirect, allostery-driven cause for their effects. Combinations of the distal-heme pocket and interfacial mutations exhibit additive effects in both structural and functional properties, contribute to our understanding of allostery, and advance protein-engineering methods for manipulating the O 2 binding properties of the hemoglobin molecule.     

Ligand migration in the truncated hemoglobin of Mycobacterium tuberculosis

The truncated hemoglobin of Mycobacterium tuberculosis (Mt-trHbO) is a small heme protein belonging to the hemoglobin superfamily. Truncated hemoglobins (trHbs) are believed to have functional roles such as terminal oxidases and oxygen sensors involved in the response to oxidative and nitrosative stress, nitric oxide (NO) detoxification, O?/NO chemistry, O? delivery under hypoxic conditions, and long-term ligand storage. Based on sequence similarities, they are classified into three groups. Experimental studies revealed that all trHbs display a 2-on-2 α-helical sandwich fold rather than the 3-on-3 α-helical sandwich fold of the classical hemoglobin fold. Using locally enhanced sampling (LESMD) molecular dynamics, the ligand-binding escape pathways from the distal heme binding cavity of Mt-trHbO were determined to better understand how this protein functions. The importance of specific residues, such as the group II and III invariant W(G8) residue, can be seen in terms of ligand diffusion pathways and ligand dynamics. LESMD simulations show that the wild-type Mt-trHbO has three diffusion pathways while the W(G8)F Mt-trHbO mutant has only two. The W(G8) residue plays a critical role in ligand binding and stabilization and helps regulate the rate of ligand escape from the distal heme pocket. Thus, this invariant residue is important in creating ligand diffusion pathways and possibly in the enzymatic functions of this protein.     

Protein fold and structure in the truncated (2/2) globin family

Analysis of amino acids sequences and protein folds has recently unraveled the structural bases and details of several proteins from the recently discovered "truncated hemoglobin" family. The analysis here presented, in agreement with previous surveys, shows that truncated hemoglobins can be classified in three main groups, based on their structural properties. Crystallographic analyses have shown that all three groups adopt a 2-on-2 alpha-helical sandwich fold, resulting from apparent editing of the classical 3-on-3 alpha-helical sandwich of vertebrate and invertebrate conventional globins. Specific structural features distinguish each of the three groups. Among these, a protein matrix tunnel system is typical of group I, a Trp residue at the G8 topological site is conserved in groups II and III, and TyrB10 is almost invariant through the three groups. A strongly intertwined network of hydrogen bonds stabilizes the heme bound ligand, despite variability of the heme distal residues observed in the different proteins considered. Details of ligand recognition in the three groups are discussed at the light of residue conservation and of differing ligand diffusion pathways to the heme. Based on structural analyses of the family-specific fold, we endorse a recent proposal of leaving the "truncated hemoglobins" term, that does not represent properly the observed 2-on-2 alpha-helical sandwich fold, and adopting the simple "2/2Hb" term to concisely address this protein family.     

Hemoglobins and hemocyanins: comparative aspects of structure and function.

Comparative studies of protein structure and function can be quite interesting by themselves, and even more interesting when interpreted with respect to an animal's physiology. In the case of fish hemoglobins, some success in the latter has been achieved but there are still many unsolved problems. It appears that comparative physiology and biochemistry have entered an era where results from comparative studies can shed a great deal of light on biochemical mechanisms in general. The trout hemoglobin system is an example. Distinctive hemoglobins in this system are presently being used as high resolution probes of the ligand-binding mechanism. Characterization of the multiple, structurally distinct subunits of the 60S Limulus hemocyanin molecule may similarly aid in understanding its function. Our studies suggest the possibility of using Limulus hemocyanin and other hemocyanins as structural homologs and analogs of more complex macromolecular arrays. The rapid development of molecular structural data from X-ray crystallographers combined with the vast data of comparative physiology and biochemistry makes this one of the most exciting areas in present day science.     

Evolution of ruminant hemoglobins. Thermodynamic divergence of ox and buffalo hemoglobins.

The ligand-binding properties of hemoglobins from two homozygote phenotypes (AA and BB) of water buffalo (Bubalus bubalis) have been characterized by equilibrium and kinetic techniques. In the case of the BB phenotype, the two constituent hemoglobins have been purified and separately analysed. Buffalo hemoglobins display the reduced sensitivity to organic phosphates characteristic of ruminant hemoglobins, their physiological effector probably being the chloride ion. In contrast to the other known hemoglobins from ruminants, all the hemoglobins from the water buffalo display a significant temperature sensitivity, the delta H for oxygen binding in the presence of physiological effectors approaching that of human hemoglobin (delta H = -30.5 kJ/mol O2). This discrepancy with the other ruminant hemoglobins (e.g. ox, delta H = -10.4 kJ/mol O2), whose primary structure is very similar to that of buffalo, hemoglobins might be correlated to the different habitat and phylogenetic history of the two subfamilies (Bos and Bubalus) of Bovidae.     

A novel thermostable hemoglobin from the actinobacterium Thermobifida fusca

The gene coding for a hemoglobin-like protein (Tf-trHb) has been identified in the thermophilic actinobacterium Thermobifida fusca and cloned in Escherichia coli for overexpression. The crystal structure of the ferric, acetate-bound derivative, was obtained at 2.48 A resolution. The three-dimensional structure of Tf-trHb is similar to structures reported for the truncated hemoglobins from Mycobacterium tuberculosis and Bacillus subtilis in its central domain. The complete lack of diffraction patterns relative to the N- and C-terminal segments indicates that these are unstructured polypeptides chains, consistent with their facile cleavage in solution. The absence of internal cavities and the presence of two water molecules between the bound acetate ion and the protein surface suggest that the mode of ligand entry is similar to that of typical hemoglobins. The protein is characterized by higher thermostability than the similar mesophilic truncated hemoglobin from B. subtilis, as demonstrated by far-UV CD melting experiments on the cyano-met derivatives. The ligand-binding properties of Tf-trHb, analyzed in stopped flow experiments, demonstrate that Tf-trHb is capable of efficient O2 binding and release between 55 and 60 degrees C, the optimal growth temperature for Thermobifida fusca.     

Crystallographic analysis of synechocystis cyanoglobin reveals the structural changes accompanying ligand binding in a hexacoordinate hemoglobin

The crystal structures of cyanide and azide-bound forms of the truncated hemoglobin from Synechocystis are presented at 1.8 angstroms resolution. A comparison with the structure of the endogenously liganded protein reveals a conformational shift unprecedented in hemoglobins, and provides the first picture of a hexacoordinate hemoglobin in both the bis-histidyl and the exogenously coordinated states. The structural changes between the different conformations are confined to two regions of the protein; the B helix, and the E helix, including the EF loop. A molecular "hinge" controlling movement of the E helix is observed in the EF loop, which is composed of three principal structural elements: Arg64, the heme-d-propionate, and a three-residue extension of the F helix. Additional features of the structural transition between the two protein conformations are discussed as they relate to the complex ligand-binding behavior observed in hexacoordinate hemoglobins, and the potential physiological function of this class of proteins.     

Heme-ligand tunneling in group I truncated hemoglobins

Truncated hemoglobins (trHbs) are small hemoproteins forming a separate cluster within the hemoglobin superfamily; their functional roles in bacteria, plants, and unicellular eukaryotes are marginally understood. Crystallographic investigations have shown that the trHb fold (a two-on-two alpha-helical sandwich related to the globin fold) hosts a protein matrix tunnel system offering a potential path for ligand diffusion to the heme distal site. The tunnel topology is conserved in group I trHbs, although with modulation of its size/structure. Here, we present a crystallographic investigation on trHbs from Mycobacterium tuberculosis, Chlamydomonas eugametos, and Paramecium caudatum, showing that treatment of trHb crystals under xenon pressure leads to binding of xenon atoms at specific (conserved) sites along the protein matrix tunnel. The crystallographic results are in keeping with data from molecular dynamics simulations, where a dioxygen molecule is left free to diffuse within the protein matrix. Modulation of xenon binding over four main sites is related to the structural properties of the tunnel system in the three trHbs and may be connected to their functional roles. In a parallel crystallographic investigation on M. tuberculosis trHbN, we show that butyl isocyanide also binds within the apolar tunnel, in excellent agreement with concepts derived from the xenon binding experiments. These results, together with recent data on atypical CO rebinding kinetics to group I trHbs, underline the potential role of the tunnel system in supporting diffusion, but also accumulation in multiple copies, of low polarity ligands/molecules within group I trHbs.     

Amino acid sequences of hemoglobins I and II from root nodules of the non-leguminous Parasponia rigida-rhizobium symbiosis, and a correction of the sequence of hemoglobin I from Parasponia andersonii

The amino acid sequence of hemoglobins I (pI 6.15 as oxyhemoglobin) and II (pI 5.64 as oxyhemoglobin) from the nitrogen-fixing root nodules of Parasponia rigida have been determined by protein sequencing. The sequence of hemoglobin I (pI 6.16, as oxyhemoglobin) from Parasponia andersonii was re-examined and the corrected primary structure, now in agreement with that predicted from the DNA sequence, is reported. The three Parasponia hemoglobins contain 161 amino acid residues (Mr approximately equal to 18,700 including the heme) with a single cysteine residue and five methionine residues. The N-terminal serine is blocked by an acetyl group. The primary structure of the Parasponia hemoglobins is highly conserved. Hemoglobins I from the two species of Parasponia are identical; both show microheterogeneity at position 30 (Asp/Glu substitution) and hemoglobin I from P. rigida shows microheterogeneity at position 150 (Ala/Val) while hemoglobin I from P. andersonii has only an Ala at 150. P. rigida hemoglobin II shows no microheterogeneity at these positions, having Asp and Val residues respectively, and it contains a single amino acid change of a Gln for an Arg at position 85, which accounts for the 0.5 unit difference in isoelectric point observed between hemoglobins I and II. The sequence data are consistent with allelic heterogeneity at a single locus rather than different genes.     

Occurrence of two architectural types of hexagonal bilayer hemoglobin in annelids: comparison of 3D reconstruction volumes of Arenicola marina and Lumbricus terrestris hemoglobins

A 3D reconstruction at 25 A resolution of native hemoglobin of the polychaete worm Arenicola marina was carried out from frozen-hydrated specimens examined in the electron microscope. The reconstruction volume of this large extracellular multimeric respiratory pigment appears as a hexagonal bilayer structure with eclipsed vertices in its upper and lower hexagonal layers. Conversely, in hemoglobins of oligochaetes, achaetes, and vestimentiferans and in chlorocruorins of the Sabellidae (polychaete) family, the vertices of the upper layer are 16 degrees clockwise rotated with respect to those of the lower layer. The fact that two other polychaete hemoglobins (Alvinella pompejana and Tylorrhynchus heterochaetus) have the same architecture as Arenicola led us to define two types of hexagonal bilayer hemoglobins/chlorocruorins: (i) type-I present in oligochaete, achaete, and vestimentiferan hemoglobins and in Sabellidae chlorocruorins; and (ii) type-II present in polychaete hemoglobins. A comparative study of the hemoglobins of Lumbricus terrestris (type-I) and Arenicola marina (type-II) showed that only two small differences located in the c4 and c5 linking units are responsible of the important architectural difference present in oligomers. A likely scheme proposed to explain the phylogenic distribution of the two types suggests that Clitellata, Sabellida (polychaete), and vestimentiferan hemoglobins and chlorocruorins derive from a type-I ancestral molecule, while Terebellida (Alvinella), Phyllodocida (Tylorrhynchus), and Scolecida (Arenicola) and possibly other polychaetes derive from an ancestor molecule with type-II hemoglobin. The architectures of the hollow globular substructures are highly similar in Arenicola and Lumbricus hemoglobins, with 12 globin chains and three linking units (c3a, c3b, and c4). The central piece of Arenicola hemoglobin is an ellipsoid while that of Lumbricus is a toroid. No phylogenic correlation could be found between the structure of the central pieces and the architecture type.     

Structure of the unliganded cAMP-dependent protein kinase catalytic subunit from Saccharomyces cerevisiae

The structure of TPK1delta, a truncated variant of the cAMP-dependent protein kinase catalytic subunit from Saccharomyces cerevisiae, was determined in an unliganded state at 2.8 A resolution and refined to a crystallographic R-factor of 19.4%. Comparison of this structure to that of its fully liganded mammalian homolog revealed a highly conserved protein fold comprised of two globular lobes. Within each lobe, root mean square deviations in Calpha positions averaged approximately equals 0.9 A. In addition, a phosphothreonine residue was found in the C-terminal domain of each enzyme. Further comparison of the two structures suggests that a trio of conformational changes accompanies ligand-binding. The first consists of a 14.7 degrees rigid-body rotation of one lobe relative to the other and results in closure of the active site cleft. The second affects only the glycine-rich nucleotide binding loop, which moves approximately equals 3 A to further close the active site and traps the nucleotide substrate. The third is localized to a C-terminal segment that makes direct contact with ligands and the ligand-binding cleft. In addition to resolving the conformation of unliganded enzyme, the model shows that the salient features of the cAMP-dependent protein kinase are conserved over long evolutionary distances.     

Evolutionary Conserved Positions Define Protein Conformational Diversity

Conformational diversity of the native state plays a central role in modulating protein function. The selection paradigm sustains that different ligands shift the conformational equilibrium through their binding to highest-affinity conformers. Intramolecular vibrational dynamics associated to each conformation should guarantee conformational transitions, which due to its importance, could possibly be associated with evolutionary conserved traits. Normal mode analysis, based on a coarse-grained model of the protein, can provide the required information to explore these features. Herein, we present a novel procedure to identify key positions sustaining the conformational diversity associated to ligand binding. The method is applied to an adequate refined dataset of 188 paired protein structures in their bound and unbound forms. Firstly, normal modes most involved in the conformational change are selected according to their corresponding overlap with structural distortions introduced by ligand binding. The subspace defined by these modes is used to analyze the effect of simulated point mutations on preserving the conformational diversity of the protein. We find a negative correlation between the effects of mutations on these normal mode subspaces associated to ligand-binding and position-specific evolutionary conservations obtained from multiple sequence-structure alignments. Positions whose mutations are found to alter the most these subspaces are defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. These positions are shown to be evolutionary conserved, mostly buried aliphatic residues localized in regular structural regions of the protein like β-sheets and α-helix.     

A hemoglobin with an optical function

Hemoglobins are best known as oxygen transport proteins. Here we describe a hemoglobin from the parasitic nematode Mermis nigrescens (Mn-GLB-E) that has an optical, light shadowing function. The protein accumulates to high concentration as intracellular crystals in the ocellus of mature phototactic adult females while also being expressed at low concentration in other tissues. It differs in sequence and expression pattern from Mn-GLB-B, a second Mermis globin. It retains the structure and oxygen-binding and light-absorbing properties typical of nematode hemoglobins. As such, recruitment to a shadowing role in the eye appears to have occurred by changes in expression without modification of biochemistry. Both globins are coded by genes interrupted by two introns at the conserved positions B12.2 and G7.0, which is in agreement with the 3exon/2intron pattern model of globin gene evolution.     

Persistently conserved positions in structurally similar, sequence dissimilar proteins: Roles in preserving protein fold and function

Many protein pairs that share the same fold do not have any detectable sequence similarity, providing a valuable source of information for studying sequence-structure relationship. In this study, we use a stringent data set of structurally similar, sequence-dissimilar protein pairs to characterize residues that may play a role in the determination of protein structure and/or function. For each protein in the database, we identify amino-acid positions that show residue conservation within both close and distant family members. These positions are termed "persistently conserved". We then proceed to determine the "mutually" persistently conserved (MPC) positions: those structurally aligned positions in a protein pair that are persistently conserved in both pair mates. Because of their intra- and interfamily conservation, these positions are good candidates for determining protein fold and function. We find that 45% of the persistently conserved positions are mutually conserved. A significant fraction of them are located in critical positions for secondary structure determination, they are mostly buried, and many of them form spatial clusters within their protein structures. A substitution matrix based on the subset of MPC positions shows two distinct characteristics: (i) it is different from other available matrices, even those that are derived from structural alignments; (ii) its relative entropy is high, emphasizing the special residue restrictions imposed on these positions. Such a substitution matrix should be valuable for protein design experiments.     

Epistatic interactions between ancestral genotype and beneficial mutations shape evolvability in Pseudomonas aeruginosa

The idea that interactions between mutations influence adaptation by driving populations to low and high fitness peaks on adaptive landscapes is deeply ingrained in evolutionary theory. Here, we investigate the impact of epistasis on evolvability by challenging populations of two Pseudomonas aeruginosa clones bearing different initial mutations (in rpoB conferring rifampicin resistance, and the type IV pili gene network) to adaptation to a medium containing l ‐serine as the sole carbon source. Despite being initially indistinguishable in fitness, populations founded by the two ancestral genotypes reached different fitness following 300 generations of evolution. Genome sequencing revealed that the difference could not be explained by acquiring mutations in different targets of selection; the majority of clones from both ancestors converged on one of the following two strategies: (1) acquiring mutations in either PA2449 (gcsR, an l ‐serine‐metabolism RpoN enhancer binding protein) or (2) protease genes. Additionally, populations from both ancestors converged on loss‐of‐function mutations in the type IV pili gene network, either due to ancestral or acquired mutations. No compensatory or reversion mutations were observed in RNA polymerase (RNAP) genes, in spite of the large fitness costs typically associated with mutations in rpoB. Although current theory points to sign epistasis as the dominant constraint on evolvability, these results suggest that the role of magnitude epistasis in constraining evolvability may be underappreciated. The contribution of magnitude epistasis is likely to be greatest under the biologically relevant mutation supply rates that make back mutations probabilistically unlikely.