Alignment of 700 globin sequences: extent of amino acid substitution and its correlation with variation in volume

Improving Protein Expression, Stability, and Function with ProteinMPNN

Natural proteins are highly optimized for function but are often difficult to produce at a scale suitable for biotechnological applications due to poor expression in heterologous systems, limited solubility, and sensitivity to temperature. Thus, a general method that improves the physical properties of native proteins while maintaining function could have wide utility for protein-based technologies. Here, we show that the deep neural network ProteinMPNN, together with evolutionary and structural information, provides a route to increasing protein expression, stability, and function. For both myoglobin and tobacco etch virus (TEV) protease, we generated designs with improved expression, elevated melting temperatures, and improved function. For TEV protease, we identified multiple designs with improved catalytic activity as compared to the parent sequence and previously reported TEV variants. Our approach should be broadly useful for improving the expression, stability, and function of biotechnologically important proteins.

Bacterial nitric oxide metabolism: Recent insights in rhizobia

Nitric oxide (NO) is a reactive gaseous molecule that has several functions in biological systems depending on its concentration. At low concentrations, NO acts as a signaling molecule, while at high concentrations, it becomes very toxic due to its ability to react with multiple cellular targets. Soil bacteria, commonly known as rhizobia, have the capacity to establish a N₂-fixing symbiosis with legumes inducing the formation of nodules in their roots. Several reports have shown NO production in the nodules where this gas acts either as a signaling molecule which regulates gene expression, or as a potent inhibitor of nitrogenase and other plant and bacteria enzymes. A better understanding of the sinks and sources of NO in rhizobia is essential to protect symbiotic nitrogen fixation from nitrosative stress. In nodules, both the plant and the microsymbiont contribute to the production of NO. From the bacterial perspective, the main source of NO reported in rhizobia is the denitrification pathway that varies significantly depending on the species. In addition to denitrification, nitrate assimilation is emerging as a new source of NO in rhizobia. To control NO accumulation in the nodules, in addition to plant haemoglobins, bacteroids also contribute to NO detoxification through the expression of a NorBC-type nitric oxide reductase as well as rhizobial haemoglobins. In the present review, updated knowledge about the NO metabolism in legume-associated endosymbiotic bacteria is summarized.

Lessons from the post-genomic era: Globin diversity beyond oxygen binding and transport

Vertebrate hemoglobin (Hb) and myoglobin (Mb) were among the first proteins whose structures and sequences were determined over 50 years ago. In the subsequent pregenomic period, numerous related proteins came to light in plants, invertebrates and bacteria, that shared the myoglobin fold, a signature sequence motif characteristic of a 3-on-3 α-helical sandwich. Concomitantly, eukaryote and bacterial globins with a truncated 2-on-2 α-helical fold were discovered. Genomic information over the last 20 years has dramatically expanded the list of known globins, demonstrating their existence in a limited number of archaeal genomes, a majority of bacterial genomes and an overwhelming majority of eukaryote genomes. In vertebrates, 6 additional globin types were identified, namely neuroglobin (Ngb), cytoglobin (Cygb), globin E (GbE), globin X (GbX), globin Y (GbY) and androglobin (Adgb). Furthermore, functions beyond the familiar oxygen transport and storage have been discovered within the vertebrate globin family, including NO metabolism, peroxidase activity, scavenging of free radicals, and signaling functions. The extension of the knowledge on globin functions suggests that the original roles of bacterial globins must have been enzymatic, involved in defense against NO toxicity, and perhaps also as sensors of O₂, regulating taxis away or towards high O₂ concentrations. In this review, we aimed to discuss the evolution and remarkable functional diversity of vertebrate globins with particular focus on the variety of non-canonical expression sites of mammalian globins and their according impressive variability of atypical functions.

Understanding molecular mechanisms in cell signaling through natural and artificial sequence variation

The functionally tolerated sequence space of proteins can now be explored in an unprecedented way, owing to the expansion of genomic databases and the development of high-throughput methods to interrogate protein function. For signaling proteins, several recent studies have shown how the analysis of sequence variation leverages the available protein-structure information to provide new insights into specificity and allosteric regulation. In this Review, we discuss recent work that illustrates how this emerging approach is providing a deeper understanding of signaling proteins.

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.

Human Sex Determination at the Edge of Ambiguity: INHERITED XY SEX REVERSAL DUE TO ENHANCED UBIQUITINATION AND PROTEASOMAL DEGRADATION OF A MASTER TRANSCRIPTION FACTOR

A general problem is posed by analysis of transcriptional thresholds governing cell fate decisions in metazoan development. A model is provided by testis determination in therian mammals. Its key step, Sertoli cell differentiation in the embryonic gonadal ridge, is initiated by SRY, a Y-encoded architectural transcription factor. Mutations in human SRY cause gonadal dysgenesis leading to XY female development (Swyer syndrome). Here, we have characterized an inherited mutation compatible with either male or female somatic phenotypes as observed in an XY father and XY daughter, respectively. The mutation (a crevice-forming substitution at a conserved back surface of the SRY high mobility group box) markedly destabilizes the domain but preserves specific DNA affinity and induced DNA bend angle. On transient transfection of diverse human and rodent cell lines, the variant SRY exhibited accelerated proteasomal degradation (relative to wild type) associated with increased ubiquitination; in vitro susceptibility to ubiquitin-independent ("default") cleavage by the 20S core proteasome was unchanged. The variant's gene regulatory activity (as assessed in a cellular model of the rat embryonic XY gonadal ridge) was reduced by 2-fold relative to wild-type SRY at similar levels of mRNA expression. Chemical proteasome inhibition restored native-like SRY expression and transcriptional activity in association with restored occupancy of a sex-specific enhancer element in principal downstream gene Sox9, demonstrating that the variant SRY exhibits essentially native activity on a per molecule basis. Our findings define a novel mechanism of impaired organogenesis, accelerated ubiquitin-directed proteasomal degradation of a master transcription factor leading to a developmental decision poised at the edge of ambiguity.

A bioinformatics insight to rhizobial globins: gene identification and mapping, polypeptide sequence and phenetic analysis, and protein modeling

Globins (Glbs) are proteins widely distributed in organisms. Three evolutionary families have been identified in Glbs: the M, S and T Glb families. The M Glbs include flavohemoglobins (fHbs) and single-domain Glbs (SDgbs); the S Glbs include globin-coupled sensors (GCSs), protoglobins and sensor single domain globins, and the T Glbs include truncated Glbs (tHbs). Structurally, the M and S Glbs exhibit 3/3-folding whereas the T Glbs exhibit 2/2-folding. Glbs are widespread in bacteria, including several rhizobial genomes. However, only few rhizobial Glbs have been characterized. Hence, we characterized Glbs from 62 rhizobial genomes using bioinformatics methods such as data mining in databases, sequence alignment, phenogram construction and protein modeling. Also, we analyzed soluble extracts from Bradyrhizobiumjaponicum USDA38 and USDA58 by (reduced + carbon monoxide (CO) minus reduced) differential spectroscopy. Database searching showed that only fhb, sdgb, gcs and thb genes exist in the rhizobia analyzed in this work. Promoter analysis revealed that apparently several rhizobial glb genes are not regulated by a -10 promoter but might be regulated by -35 and Fnr (fumarate-nitrate reduction regulator)-like promoters. Mapping analysis revealed that rhizobial fhbs and thbs are flanked by a variety of genes whereas several rhizobial sdgbs and gcss are flanked by genes coding for proteins involved in the metabolism of nitrates and nitrites and chemotaxis, respectively. Phenetic analysis showed that rhizobial Glbs segregate into the M, S and T Glb families, while structural analysis showed that predicted rhizobial SDgbs and fHbs and GCSs globin domain and tHbs fold into the 3/3- and 2/2-folding, respectively. Spectra from B. japonicum USDA38 and USDA58 soluble extracts exhibited peaks and troughs characteristic of bacterial and vertebrate Glbs thus indicating that putative Glbs are synthesized in B. japonicum USDA38 and USDA58.

Structure-function relationships in human testis-determining factor SRY: an aromatic buttress underlies the specific DNA-bending surface of a high mobility group (HMG) box

Human testis determination is initiated by SRY, a Y-encoded architectural transcription factor. Mutations in SRY cause 46 XY gonadal dysgenesis with female somatic phenotype (Swyer syndrome) and confer a high risk of malignancy (gonadoblastoma). Such mutations cluster in the SRY high mobility group (HMG) box, a conserved motif of specific DNA binding and bending. To explore structure-function relationships, we constructed all possible substitutions at a site of clinical mutation (W70L). Our studies thus focused on a core aromatic residue (position 15 of the consensus HMG box) that is invariant among SRY-related HMG box transcription factors (the SOX family) and conserved as aromatic (Phe or Tyr) among other sequence-specific boxes. In a yeast one-hybrid system sensitive to specific SRY-DNA binding, the variant domains exhibited reduced (Phe and Tyr) or absent activity (the remaining 17 substitutions). Representative nonpolar variants with partial or absent activity (Tyr, Phe, Leu, and Ala in order of decreasing side-chain volume) were chosen for study in vitro and in mammalian cell culture. The clinical mutation (Leu) was found to markedly impair multiple biochemical and cellular activities as respectively probed through the following: (i) in vitro assays of specific DNA binding and protein stability, and (ii) cell culture-based assays of proteosomal degradation, nuclear import, enhancer DNA occupancy, and SRY-dependent transcriptional activation. Surprisingly, however, DNA bending is robust to this or the related Ala substitution that profoundly impairs box stability. Together, our findings demonstrate that the folding, trafficking, and gene-regulatory function of SRY requires an invariant aromatic "buttress" beneath its specific DNA-bending surface.

The diversity of 2/2 (truncated) globins

Small size globins that have been defined as 'truncated haemoglobins' or as '2/2 haemoglobins' have increasingly been discovered in microorganisms since the early 1990s. Analysis of amino acid sequences allowed to distinguish three groups that collect proteins with specific and common structural properties. All three groups display 3D structures that are based on four main α-helices, which are a subset of the conventional eight-helices globin fold. Specific features, such as the presence of protein matrix tunnels that are held to promote diffusion of functional ligands to/from the haem, distinguish members of the three groups. Haem distal sites vary for their accessibility, local structures, polarity, and ligand stabilization mechanisms, suggesting functional roles that are related to O2/NO chemistry. In a few cases, such activities have been proven in vitro and in vivo through deletion mutants. The issue of 2/2 haemoglobin varied biological functions throughout the three groups remains however fully open.

The evolution of land plant hemoglobins

This review discusses the evolution of land plant hemoglobins within the broader context of eukaryote hemoglobins and the three families of bacterial globins. Most eukaryote hemoglobins, including metazoan globins and the symbiotic and non-symbiotic plant hemoglobins, are homologous to the bacterial 3/3-fold flavohemoglobins. The remaining plant hemoglobins are homologous to the bacterial 2/2-fold group 2 hemoglobins. We have proposed that all eukaryote globins were acquired via horizontal gene transfer concomitant with the endosymbiotic events responsible for the origin of mitochondria and chloroplasts. Although the 3/3 hemoglobins originated in the ancestor of green algae and plants prior to the emergence of embryophytes at about 450 mya, the 2/2 hemoglobins appear to have originated via horizontal gene transfer from a bacterium ancestral to present day Chloroflexi. Unlike the 2/2 hemoglobins, the evolution of the 3/3 hemoglobins was accompanied by duplication, diversification, and functional adaptations. Duplication of the ancestral plant nshb gene into the nshb-1 and nshb-2 lineages occurred prior to the monocot-dicot divergence at ca. 140 mya. It was followed by the emergence of symbiotic hemoglobins from a non-symbiotic hemoglobin precursor and further specialization, leading to leghemoglobins in N₂-fixing legume nodules concomitant with the origin of nodulation at ca. 60 mya. The transition of non-symbiotic to symbiotic hemoglobins (including to leghemoglobins) was accompanied by the alteration of heme-Fe coordination from hexa- to penta-coordination. Additional genomic information about Charophyte algae, the sister group to land plants, is required for the further clarification of plant globin phylogeny.

Phylogenetic relationships of 3/3 and 2/2 hemoglobins in Archaeplastida genomes to bacterial and other eukaryote hemoglobins

Land plants and algae form a supergroup, the Archaeplastida, believed to be monophyletic. We report the results of an analysis of the phylogeny of putative globins in the currently available genomes to bacterial and other eukaryote hemoglobins (Hbs). Archaeplastida genomes have 3/3 and 2/2 Hbs, with the land plant genomes having group 2 2/2 Hbs, except for the unexpected occurrence of two group 1 2/2 Hbs in Ricinus communis. Bayesian analysis shows that plant 3/3 Hbs are related to vertebrate neuroglobins and bacterial flavohemoglobins (FHbs). We sought to define the bacterial groups, whose ancestors shared the precursors of Archaeplastida Hbs, via Bayesian and neighbor-joining analyses based on COBALT alignment of representative sets of bacterial 3/3 FHb-like globins and group 1 and 2 2/2 Hbs with the corresponding Archaeplastida Hbs. The results suggest that the Archaeplastida 3/3 and group 1 2/2 Hbs could have originated from the horizontal gene transfers (HGTs) that accompanied the two generally accepted endosymbioses of a proteobacterium and a cyanobacterium with a eukaryote ancestor. In contrast, the origin of the group 2 2/2 Hbs unexpectedly appears to involve HGT from a bacterium ancestral to Chloroflexi, Deinococcales, Bacilli, and Actinomycetes. Furthermore, although intron positions and phases are mostly conserved among the land plant 3/3 and 2/2 globin genes, introns are absent in the algal 3/3 genes and intron positions and phases are highly variable in their 2/2 genes. Thus, introns are irrelevant to globin evolution in Archaeplastida.

Nitrosyl hydride (HNO) as an O2 analogue: long-lived HNO adducts of ferrous globins

Nitrosyl hydride, HNO or nitroxyl, is the one-electron reduced and protonated form of nitric oxide. HNO is isoelectronic to singlet O(2), and we have previously reported that deoxymyoglobin traps free HNO to form a stable adduct. In this report, we demonstrate that oxygen-binding hemoglobins from human, soy, and clam also trap HNO to form adducts which are stable over a period of weeks. The same species can be formed in higher yields by careful reduction of the ferrous nitrosyl adducts of the proteins. Like the analogous O(2)-Fe(II) adducts, the HNO adducts are diamagnetic, but with a characteristic HNO resonance in (1)H NMR at ca. 15 ppm that splits into doublets for H(15)NO adducts. The (1)H and (15)N NMR resonances, obtained by HSQC experiments, are shown to differentiate subunits and isoforms of proteins within mixtures. An apparent difference in the reduction rates of the NO adducts of the two subunits of human hemoglobin allows assignment of two distinct nitrosyl hydride peaks by a combination of UV-vis, NMR, and EPR analysis. The two peaks of the HNO-hHb adduct have a persistent 3:1 ratio during trapping reactions, demonstrating a kinetic difference between HNO binding at the two subunits. These results show NMR characterization of ferrous HNO adducts as a unique tool sensitive to structural changes within the oxygen-binding cavity, which may be of use in defining modes of oxygen binding in other heme proteins and enzymes.

Rapid comparison of properties on protein surface

The mapping of physicochemical characteristics onto the surface of a protein provides crucial insights into its function and evolution. This information can be further used in the characterization and identification of similarities within protein surface regions. We propose a novel method which quantitatively compares global and local properties on the protein surface. We have tested the method on comparison of electrostatic potentials and hydrophobicity. The method is based on 3D Zernike descriptors, which provides a compact representation of a given property defined on a protein surface. Compactness and rotational invariance of this descriptor enable fast comparison suitable for database searches. The usefulness of this method is exemplified by studying several protein families including globins, thermophilic and mesophilic proteins, and active sites of TIM beta/alpha barrel proteins. In all the cases studied, the descriptor is able to cluster proteins into functionally relevant groups. The proposed approach can also be easily extended to other surface properties. This protein surface-based approach will add a new way of viewing and comparing proteins to conventional methods, which compare proteins in terms of their primary sequence or tertiary structure.

Modulation of the enzymatic efficiency of ferredoxin-NADP(H) reductase by the amino acid volume around the catalytic site

Ferredoxin (flavodoxin)-NADP(H) reductases (FNRs) are ubiquitous flavoenzymes that deliver NADPH or low-potential one-electron donors (ferredoxin, flavodoxin, adrenodoxin) to redox-based metabolic reactions in plastids, mitochondria and bacteria. Plastidic FNRs are quite efficient reductases. In contrast, FNRs from organisms possessing a heterotrophic metabolism or anoxygenic photosynthesis display turnover numbers 20- to 100-fold lower than those of their plastidic and cyanobacterial counterparts. Several structural features of these enzymes have yet to be explained. The residue Y308 in pea FNR is stacked nearly parallel to the re-face of the flavin and is highly conserved amongst members of the family. By computing the relative free energy for the lumiflavin-phenol pair at different angles with the relative position found for Y308 in pea FNR, it can be concluded that this amino acid is constrained against the isoalloxazine. This effect is probably caused by amino acids C266 and L268, which face the other side of this tyrosine. Simple and double FNR mutants of these amino acids were obtained and characterized. It was observed that a decrease or increase in the amino acid volume resulted in a decrease in the catalytic efficiency of the enzyme without altering the protein structure. Our results provide experimental evidence that the volume of these amino acids participates in the fine-tuning of the catalytic efficiency of the enzyme.

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

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

Wide diversity in structure and expression profiles among members of the Caenorhabditis elegans globin protein family

BACKGROUND

The emergence of high throughput genome sequencing facilities and powerful high performance bioinformatic tools has highlighted hitherto unexpected wide occurrence of globins in the three kingdoms of life. In silico analysis of the genome of C. elegans identified 33 putative globin genes. It remains a mystery why this tiny animal might need so many globins. As an inroad to understanding this complexity we initiated a structural and functional analysis of the globin family in C. elegans.

RESULTS

All 33 C. elegans putative globin genes are transcribed. The translated sequences have the essential signatures of single domain bona fide globins, or they contain a distinct globin domain that is part of a larger protein. All globin domains can be aligned so as to fit the globin fold, but internal interhelical and N- and C-terminal extensions and a variety of amino acid substitutions generate much structural diversity among the globins of C. elegans. Likewise, the encoding genes lack a conserved pattern of intron insertion positioning. We analyze the expression profiles of the globins during the progression of the life cycle, and we find that distinct subsets of globins are induced, or repressed, in wild-type dauers and in daf-2(e1370)/insulin-receptor mutant adults, although these animals share several physiological features including resistance to elevated temperature, oxidative stress and hypoxic death. Several globin genes are upregulated following oxygen deprivation and we find that HIF-1 and DAF-2 each are required for this response. Our data indicate that the DAF-2 regulated transcription factor DAF-16/FOXO positively modulates hif-1 transcription under anoxia but opposes expression of the HIF-1 responsive globin genes itself. In contrast, the canonical globin of C. elegans, ZK637.13, is not responsive to anoxia. Reduced DAF-2 signaling leads to enhanced transcription of this globin and DAF-16 is required for this effect.

CONCLUSION

We found that all 33 putative globins are expressed, albeit at low or very low levels, perhaps indicating cell-specific expression. They show wide diversity in gene structure and amino acid sequence, suggesting a long evolutionary history. Ten globins are responsive to oxygen deprivation in an interacting HIF-1 and DAF-16 dependent manner. Globin ZK637.13 is not responsive to oxygen deprivation and regulated by the Ins/IGF pathway only suggesting that this globin may contribute to the life maintenance program.

Exploiting a list of protein sequences

We describe a software program to help exploit a database of aligned protein sequences. In addition to the classical lists of sequences, a graphical representation is used to get a better overview of the information. As natural parameters, the type of amino acid and sequence position are used. Various plots or 3D representations are then updated. Examples are shown based on globin sequences from various species and on the abnormal human hemoglobins. The software should be of interest to protein engineers who need to know what variants are already known.

A model of globin evolution

Putative globins have been identified in 426 bacterial, 32 Archaeal and 67 eukaryote genomes. Among these sequences are the hitherto unsuspected presence of single domain sensor globins within Bacteria, Fungi, and a Euryarchaeote. Bayesian phylogenetic trees suggest that their occurrence in the latter two groups could be the result of lateral gene transfer from Bacteria. Iterated psiblast searches based on groups of globin sequences indicate that bacterial flavohemoglobins are closer to metazoan globins than to the other two lineages, the 2-over-2 globins and the globin-coupled sensors. Since Bacteria is the only kingdom to have all the subgroups of the three globin lineages, we propose a working model of globin evolution based on the assumption that all three lineages originated and evolved only in Bacteria. Although the 2-over-2 globins and the globin-coupled sensors recognize flavohemoglobins, there is little recognition between them. Thus, in the first stage of globin evolution, we favor a flavohemoglobin-like single domain protein as the ancestral globin. The next stage comprised the splitting off to single domain 2-over-2 and sensor-like globins, followed by the covalent addition of C-terminal domains resulting in the chimeric flavohemoglobins and globin-coupled sensors. The last stage encompassed the lateral gene transfers of some members of the three globin lineages to specific groups of Archaea and Eukaryotes.

Molecular characterization of two myoglobins of Paragonimus westermani

Myoglobins (Mbs), globin proteins, are present in high concentrations in trematodes. In Paragonimus westermani, 2 cDNAs were found to encode Mbs. The first clone, Pwmyo1, codes a total of 149 amino acids with a calculated mass of 16.6 kDa. The second, Pwmyo2, encodes a 146-amino acid protein with a calculated mass of 16.2 kDa. The predicted secondary structures showed the presence of 8 helices, which is the basic characteristic of Mbs. Sequence alignment revealed a high homology with the other trematode Mbs. The 2 clones contained the characteristic tyrosyl residues at helical positions B10 and distal E7, which are substitutions that have been previously shown to contribute to the high oxygen affinity of Mbs. Polyclonal antibodies against the recombinant Mbs were raised with no cross-reactivity observed. Immunolocalization revealed the proteins to be distributed generally throughout the parenchymal tissues, but absent from the tegument and reproductive organs. The cell mass of the eggs of the worm stained positive to Pwmyo2 but not Pwmyo1, suggesting the stage-specific expression of these Mbs.

A phylogenomic profile of globins

Background

Globins occur in all three kingdoms of life: they can be classified into single-domain globins and chimeric globins. The latter comprise the flavohemoglobins with a C-terminal FAD-binding domain and the gene-regulating globin coupled sensors, with variable C-terminal domains. The single-domain globins encompass sequences related to chimeric globins and «truncated» hemoglobins with a 2-over-2 instead of the canonical 3-over-3 α-helical fold.

Results

A census of globins in 26 archaeal, 245 bacterial and 49 eukaryote genomes was carried out. Only ~25% of archaea have globins, including globin coupled sensors, related single domain globins and 2-over-2 globins. From one to seven globins per genome were found in ~65% of the bacterial genomes: the presence and number of globins are positively correlated with genome size. Globins appear to be mostly absent in Bacteroidetes/Chlorobi, Chlamydia, Lactobacillales, Mollicutes, Rickettsiales, Pastorellales and Spirochaetes. Single domain globins occur in metazoans and flavohemoglobins are found in fungi, diplomonads and mycetozoans. Although red algae have single domain globins, including 2-over-2 globins, the green algae and ciliates have only 2-over-2 globins. Plants have symbiotic and nonsymbiotic single domain hemoglobins and 2-over-2 hemoglobins. Over 90% of eukaryotes have globins: the nematode Caenorhabditis has the most putative globins, ~33. No globins occur in the parasitic, unicellular eukaryotes such as Encephalitozoon, Entamoeba, Plasmodium and Trypanosoma.

Conclusion

Although Bacteria have all three types of globins, Archaeado not have flavohemoglobins and Eukaryotes lack globin coupled sensors. Since the hemoglobins in organisms other than animals are enzymes or sensors, it is likely that the evolution of an oxygen transport function accompanied the emergence of multicellular animals.

Three globin lineages belonging to two structural classes in genomes from the three kingdoms of life

Although most globins, including the N-terminal domains within chimeric proteins such as flavohemoglobins and globin-coupled sensors, exhibit a 3/3 helical sandwich structure, many bacterial, plant, and ciliate globins have a 2/2 helical sandwich structure. We carried out a comprehensive survey of globins in the genomes from the three kingdoms of life. Bayesian phylogenetic trees based on manually aligned sequences indicate the possibility of past horizontal globin gene transfers from bacteria to eukaryotes. blastp searches revealed the presence of 3/3 single-domain globins related to the globin domains of the bacterial and fungal flavohemoglobins in many bacteria, a red alga, and a diatom. Iterated psi-blast searches based on groups of globin sequences found that only the single-domain globins and flavohemoglobins recognize the eukaryote 3/3 globins, including vertebrate neuroglobins, alpha- and beta-globins, and cytoglobins. The 2/2 globins recognize the flavohemoglobins, as do the globin coupled sensors and the closely related single-domain protoglobins. However, the 2/2 globins and the globin-coupled sensors do not recognize each other. Thus, all globins appear to be distributed among three lineages: (i) the 3/3 plant and metazoan globins, single-domain globins, and flavohemoglobins; (ii) the bacterial 3/3 globin-coupled sensors and protoglobins; and (iii) the bacterial, plant, and ciliate 2/2 globins. The three lineages may have evolved from an ancestral 3/3 or 2/2 globin. Furthermore, it appears likely that the predominant functions of globins are enzymatic and that oxygen transport is a specialized development that accompanied the evolution of metazoans.

A new progressive-iterative algorithm for multiple structure alignment

MOTIVATION

Multiple structure alignments are becoming important tools in many aspects of structural bioinformatics. The current explosion in the number of available protein structures demands multiple structural alignment algorithms with an adequate balance of accuracy and speed, for large scale applications in structural genomics, protein structure prediction and protein classification.

RESULTS

A new multiple structural alignment program, MAMMOTH-mult, is described. It is demonstrated that the alignments obtained with the new method are an improvement over previous manual or automatic alignments available in several widely used databases at all structural levels. Detailed analysis of the structural alignments for a few representative cases indicates that MAMMOTH-mult delivers biologically meaningful trees and conservation at the sequence and structural levels of functional motifs in the alignments. An important improvement over previous methods is the reduction in computational cost. Typical alignments take only a median time of 5 CPU seconds in a single R12000 processor. MAMMOTH-mult is particularly useful for large scale applications.

AVAILABILITY

http://ub.cbm.uam.es/mammoth/mult.

Neuroglobin and cytoglobin: Two new entries in the hemoglobin superfamily*

Neuroglobin (Ngb) and cytoglobin (Cygb) are two newly discovered intracellular members of the vertebrate hemoglobin (Hb) family. Ngb, predominantly expressed in nerve cells, is of ancient evolutionary origin and is homologous to nerve-globins of invertebrates. Cygb, present in many different tissues, shares common ancestry with myoglobin (Mb) and can be traced to early vertebrate evolution. Ngb and Cygb display the classical three-on-three α-helical globin fold and are endowed with a hexa-coordinate heme Fe atom, in both their ferrous and ferric forms, having the heme distal HisE7 residue as the endogenous sixth ligand. Reversible intramolecular hexa- to penta-coordination of the heme Fe atom modulates Ngb and Cygb ligand-binding properties. In Ngb and Cygb, ligand migration to/from the heme distal site may be assisted by protein/matrix tunnel cavity systems. The physiological roles of Ngb and Cygb are poorly understood. Ngb may protect neuronal cells from hypoxic-ischemic insults, may act as oxidative stress-responsive sensor protein, and may sustain NO/O(2) scavenging and/or reactive oxygen species (ROS) detoxification. Cygb, located in the cytoplasm of fibroblasts, chondroblasts, osteoblasts, and hepatic stellate cells, has been hypothesized to be involved in collagen synthesis. In neurons, Cygb, located in both cytoplasm and nucleus, may provide O(2) for enzymatic reactions, and may be involved in a ROS (NO)-signaling pathway(s). Here, we review current knowledge on Ngb and Cygb in terms of their structure, function, and evolutionary links to the well-known human HbA and Mb.

Ancestral hemoglobins in Archaea

Hemoglobins are ubiquitous in Eukarya and Bacteria but, until now, have not been found in Archaea. A phylogenetic analysis of the recently revealed microbial family of globin-coupled heme-based sensors suggests that these sensors descended from an ancient globin-only progenitor, or a protoglobin (Pgb). Here, we report the discovery and characterization of two Pgbs from the Archaea: ApPgb from the obligately aerobic hyperthermophile Aeropyrum pernix, and MaPgb from the strictly anaerobic methanogen Methanosarcina acetivorans. Both ApPgb and MaPgb bind molecular oxygen, nitric oxide, and carbon monoxide by means of a heme moiety that is coordinated to the protein through the F8 histidine (histidine 120). We postulate that these archaeal globins are the ancestors of contemporary hemoglobins.

An electrospray ionization mass spectrometric study of the subunit structure of the giant hemoglobin from the leech Nephelopsis oscura

The subunit structure of the giant, extracellular hexagonal bilayer (HBL) hemoglobin (Hb) from the leech Nephelopsis oscura was investigated by electrospray ionization mass spectrometry (ESI-MS) employing a maximum entropy deconvolution of its complex, multiply charged ESI spectra. The denatured unreduced Hb consisted of three monomer globin chains (M), a1 = 16535 Da, a2 = 17171 Da and a3 = 17315 Da, five nonglobin linker chains, L1 = 24512 Da, L2 = 24586 Da, L3 = 24979 Da, L4 = 25006 Da, and L5 = 25566 Da and two subunits of 32950 Da and 33125 Da. ESI-MS of the denatured, reduced Hb showed that the latter were disulfide-bonded heterodimers (D) of globin chains b1 = 16322 Da and b2 = 16499 Da with chain c = 16632 Da. Time-of-flight ESI-MS of the Hb at pH 3.8, 4.5, 5.0, 5.8 and 7.0 revealed a distribution of charge states from 32(+) to 37(+) with masses decreasing from 211 to 208.5 kDa with increase in cone voltage from 60 to 160 V, indicating the presence of a subassembly comprising 12 globin chains. The subunit composition 6M + 3D + 12h, where M = 16993 Da and D = 33004 Da are the weighted masses and h = 616.5 Da, provides a calculated mass, 208.37 kDa that is closest to 208.5 kDa. Our experimental findings are consistent with the bracelet model of HBL Hbs, verified by the recent low-resolution crystal structure of Lumbricus Hb, wherein an HBL arrangement of 12 globin dodecamer subassemblies is tethered to a central complex of 36 linker chains for a total mass of 208.37 x 12 + 24.94 x 36 = 3398 kDa.

The redox state of the cell regulates the ligand binding affinity of human neuroglobin and cytoglobin

Neuroglobin and cytoglobin reversibly bind oxygen in competition with the distal histidine, and the observed oxygen affinity therefore depends on the properties of both ligands. In the absence of an external ligand, the iron atom of these globins is hexacoordinated. There are three cysteine residues in human neuroglobin; those at positions CD7 and D5 are sufficiently close to form an internal disulfide bond. Both cysteine residues in cytoglobin, although localized in other positions than in human neuroglobin, may form a disulfide bond as well. The existence and position of these disulfide bonds was demonstrated by mass spectrometry and thiol accessibility studies. Mutation of the cysteines involved, or the use of reducing agents to break the S-S bond, led to a decrease in the observed oxygen affinity of human neuroglobin by an order of magnitude. The critical parameter is the histidine dissociation rate, which changes by about a factor of 10. The same effect is observed with human cytoglobin, although to a much lesser extent (less than a factor of 2). These results suggest a novel mechanism for the regulation of oxygen binding; contact with an appropriate electron donor would provoke the release of oxygen. Hence the oxygen affinity would be directly linked to the redox state of the cell.

1H NMR study of the molecular structure and magnetic properties of the active site for the cyanomet complex of O2-avid hemoglobin from the trematode Paramphistomum epiclitum

The solution molecular and electronic structures of the active site in the extremely O2-avid hemoglobin from the trematode Paramphistomum epiclitum have been investigated by 1H NMR on the cyanomet form in order to elucidate the distal hydrogen-bonding to a ligated H-bond acceptor ligand. Comparison of the strengths of dipolar interactions in solution with the alternate crystal structures of methemoglobin establish that the solution structure of wild-type Hb more closely resembles the crystal structure of the recombinant wild-type than the true wild-type met-hemoglobin. The distal Tyr66(E7) is found oriented out of the heme pocket in solution as found in both crystal structures. Analysis of dipolar contacts, dipolar shift and paramagnetic relaxation establishes that the Tyr32(B10) hydrogen proton adopts an orientation that allows it to make a strong H-bond to the bound cyanide. The observation of a significant isotope effect on the heme methyl contact shifts confirms a strong contact between the Tyr32(B10) OH and the ligated cyanide. The quantitative determination of the orientation and anisotropies of the paramagnetic susceptibility tensor reveal that the cyanide is tilted approximately 10 degrees from the heme normal so as to avoid van der Waals overlap with the Tyr32(B10) Oeta. The pattern of heme contact shifts with large low-field shifts for 7-CH3 and 18-CH3 is shown to arise not from the 180 degrees rotation about the alpha-gamma-meso axis, but due to the approximately 45 degrees rotation of the axial His imidazole ring, relative to that in mammalian globins.

A TyrCD1/TrpG8 hydrogen bond network and a TyrB10TyrCD1 covalent link shape the heme distal site of Mycobacterium tuberculosis hemoglobin O

Truncated hemoglobins (Hbs) are small hemoproteins, identified in microorganisms and in some plants, forming a separate cluster within the Hb superfamily. Two distantly related truncated Hbs, trHbN and trHbO, are expressed at different developmental stages in Mycobacterium tuberculosis. Sequence analysis shows that the two proteins share 18% amino acid identities and belong to different groups within the truncated Hb cluster. Although a specific defense role against nitrosative stress has been ascribed to trHbN (expressed during the Mycobacterium stationary phase), no clear functions have been recognized for trHbO, which is expressed throughout the Mycobacterium growth phase. The 2.1-A crystal structure of M. tuberculosis cyano-met trHbO shows that the protein assembles in a compact dodecamer. Six of the dodecamer subunits are characterized by a double conformation for their CD regions and, most notably, by a covalent bond linking the phenolic O atom of TyrB10 to the aromatic ring of TyrCD1, in the heme distal cavity. All 12 subunits display a cyanide ion bound to the heme Fe atom, stabilized by a tight hydrogen-bonded network based on the (globin very rare) TyrCD1 and TrpG8 residues. The small apolar AlaE7 residue leaves room for ligand access to the heme distal site through the conventional "E7 path," as proposed for myoglobin. Different from trHbN, where a 20-A protein matrix tunnel is held to sustain ligand diffusion to an otherwise inaccessible heme distal site, the topologically related region in trHbO hosts two protein matrix cavities.

Globin and linker sequences of the giant extracellular hemoglobin from the leech Macrobdella decora

A detailed electrospray ionization mass spectrometric study of the approximately 3.5-MDa hexagonal bilayer hemoglobin (HBL Hb) from the pond leech Macrobdella decora has shown it to consist of at least six approximately 17-kDa globin chains, of which two are monomeric and the remaining four occur as disulfide-bonded heterodimers, and three approximately 24-kDa nonglobin linker chains (Weber et al., J. Mol. Biol. 251: 703-720, 1995). The cDNA sequences of the five major constituent chains, globin chains IIA, IIB, B, and C and linker chain L1, are reported here. The globins and linkers share 30%-50% and 20%-30% identity, respectively, with other annelid sequences. Furthermore, IIB and C align with strain A of annelid sequences, whereas IIA and B align with the strain B sequences. Although chains B and C are monomeric, chains IIA and IIB form the main disulfide-bonded dimer. They also have some unusual features: the distal His (E7) is replaced by Phe in IIA, and the highly conserved CD1Phe is replaced by Leu in IIB. In spite of these unusual features, the functional properties of Macrobdella Hb are comparable to those of other HBL Hbs. A phylogenetic analysis of the globin sequences from Macrobdella, the polychaete Tylorrhynchus, the oligochaete Lumbricus, and the vestimentiferan Lamellibrachia, indicates that the two strains originated by gene duplication followed by additional duplication of each of the two strains. The mutation rate of the linkers appeared to be faster than that of the globin chains. The phylogenetic trees constructed using the Maximum Likelihood, Neighbor-Joining and Fitch methods showed the Macrobdella globin sequences to be closest to Lumbricus, in agreement with a view of annelid evolution in which the divergence of the polychaetes occurred before the divergence of the leeches from oligochaetes.

Structure and function of the globin and globin gene from the Antarctic mollusc Yoldia eightsi

The mechanism of adaptation of haemoglobin from the Antarctic mollusc Yoldia eightsi to its low-temperature environment is a decrease in the oxygen affinity via an increased ligand-dissociation rate. At 2 degrees C this haemoglobin has an oxygen affinity similar to other haemoglobins at 25 degrees C. At 25 degrees C, Yoldia haemoglobin shows a low oxygen affinity, resembling that of human deoxyhaemoglobin. The mechanism involves a lower binding energy to oxygen, suggesting a loss or weakening of the usual hydrogen bond, leading to a higher oxygen-dissociation rate. However, Yoldia haemoglobin has the usual distal and proximal histidines, so the primary structure alone does not provide an obvious explanation for the low affinity. The CO-binding kinetics are biphasic, with the fraction of slow phase increasing at higher protein concentrations, indicating the formation of dimers or a higher level of polymerization. The protein-protein interaction appears to be of hydrophobic nature, since it can be partially reversed by addition of ethylene glycol as co-solvent. While the CO-association rates differ by a factor of 10, the oxygen equilibrium data could be simulated with a single affinity. The Yoldia haemoglobin gene contains three introns, interrupting the coding region at position NA1.2, B12.2 and G7.0. The conservation of the B12.2 and G7.0 introns is in contrast with the unprecedented NA1.2 intron. Phylogenetic analyses reveal a gene tree where the Yoldia haemoglobin gene is separated from other mollusc globin genes, confirming the specific adaptation of the Yoldia haemoglobin.

Neuroglobin and cytoglobin. Fresh blood for the vertebrate globin family

Neuroglobin and cytoglobin are two recently discovered members of the vertebrate globin family. Both are intracellular proteins endowed with hexacoordinated heme-Fe atoms, in their ferrous and ferric forms, and display O2 affinities comparable with that of myoglobin. Neuroglobin, which is predominantly expressed in nerve cells, is thought to protect neurons from hypoxic-ischemic injury. It is of ancient evolutionary origin, and is homologous to nerve globins of invertebrates. Cytoglobin is expressed in many different tissues, although at varying levels. It shares common ancestry with myoglobin, and can be traced to early vertebrate evolution. The physiological roles of neuroglobin and cytoglobin are not completely understood. Although supplying cells with O2 is the likely function, it is also possible that both globins act as O2-consuming enzymes or as O2 sensors. Here, we review what is currently known about neuroglobin and cytoglobin in terms of their function, tissue distribution and relatedness to the well-known hemoglobin and myoglobin. Strikingly, the data reveal that O2 metabolism in cells is more complicated than was thought before, requiring unexpected O2-binding proteins with potentially novel functional features.

Crystal structures of unligated and CN-ligated Glycera dibranchiata monomer ferric hemoglobin components III and IV

Erythrocytes of the marine annelid, Glycera dibranchiata, contain a mixture of monomeric and polymeric hemoglobins. There are three major monomer hemoglobin components, II, III, IV (also called GMH2, 3, and 4), that have been highly purified and well characterized. We have now crystallized GMH3 and GMH4 and determined their structures to 1.4-1.8 A resolution. The structures were determined for these two monomer hemoglobins in the oxidized (Fe3+, ferric, or met-) forms in both the unligated and cyanide-ligated states. This work differs from two published, refined structures of a Glycera dibranchiata monomer hemoglobin, which has a sequence that is substantially different from any bona fide major monomer hemoglobins (GMH2, 3, or 4). The high-resolution crystal structures (presented here) and the previous NMR structure of CO-ligated GMH4, provide a basis for interpreting structure/function details of the monomer hemoglobins. These details include: (1) the strong correlation between temperature factor and NMR dynamics for respective protein forms; (2) the unique nature of the HisE7Leu primary sequence substitutions in GMH3 and GMH4 and their impact on cyanide ion binding kinetics; (3) the LeuB10Phe difference between GMH3 and GMH4 and its impact on ligand binding; and (4) elucidation of changes in the structural details of the distal and proximal heme pockets upon cyanide binding.

Evolution of the sulfide-binding function within the globin multigenic family of the deep-sea hydrothermal vent tubeworm Riftia pachyptila

The giant extracellular hexagonal bilayer hemoglobin (HBL-Hb) of the deep-sea hydrothermal vent tube worm Riftia pachyptila is able to transport simultaneously O(2) and H(2)S in the blood from the gills to a specific organ: the trophosome that harbors sulfide-oxidizing endosymbionts. This vascular HBL-Hb is made of 144 globins from which four globin types (A1, A2, B1, and B2) coevolve. The H(2)S is bound at a specific location (not on the heme site) onto two of these globin types. In order to understand how such a function emerged and evolved in vestimentiferans and other related annelids, six partial cDNAs corresponding to the six globins known to compose the multigenic family of R. pachyptila have been identified and sequenced. These partial sequences (ca. 120 amino acids, i.e., 80% of the entire protein) were used to reconstruct molecular phylogenies in order to trace duplication events that have led to the family organization of these globins and to locate the position of the free cysteine residues known to bind H(2)S. From these sequences, only two free cysteine residues have been found to occur, at positions Cys + 1 (i.e., 1 a.a. from the well-conserved distal histidine) and Cys + 11 (i.e., 11 a.a. from the same histidine) in globins B2 and A2, respectively. These two positions are well conserved in annelids, vestimentiferans, and pogonophorans, which live in sulfidic environments. The structural comparison of the hydrophobic environment that surrounds these cysteine residues (the sulfide-binding domain) using hydrophobic cluster analysis plots, together with the cysteine positions in paralogous strains, suggests that the sulfide-binding function might have emerged before the annelid radiation in order to detoxify this toxic compound. Moreover, globin evolutionary rates are highly different between paralogous strains. This suggests that either the two globin subfamilies involved in the sulfide-binding function (A2 and B2) have evolved under strong directional selective constraints (negative selection) and that the two other globins (A1 and B1) have accumulated more substitutions through positive selection or have evolved neutrally after a relaxation of selection pressures. A likely scenario on the evolution of this multigenic family is proposed and discussed from this data set.

The 109 residue nerve tissue minihemoglobin from Cerebratulus lacteus highlights striking structural plasticity of the alpha-helical globin fold

A very short hemoglobin (CerHb; 109 amino acids) binds O(2) cooperatively in the nerve tissue of the nemertean worm Cerebratulus lacteus to sustain neural activity during anoxia. Sequence analysis suggests that CerHb tertiary structure may be unique among the known globin fold evolutionary variants. The X-ray structure of oxygenated CerHb (R factor 15.3%, at 1.5 A resolution) displays deletion of the globin N-terminal A helix, an extended GH region, a very short H helix, and heme solvent shielding based on specific aromatic residues. The heme-bound O(2) is stabilized by hydrogen bonds to the distal TyrB10-GlnE7 pair. Ligand access to heme may take place through a wide protein matrix tunnel connecting the distal site to a surface cleft located between the E and H helices.

Globin-coupled sensors: a class of heme-containing sensors in Archaea and Bacteria

The recently discovered prokaryotic signal transducer HemAT, which has been described in both Archaea and Bacteria, mediates aerotactic responses. The N-terminal regions of HemAT from the archaeon Halobacterium salinarum (HemAT-Hs) and from the Gram-positive bacterium Bacillus subtilis (HemAT-Bs) contain a myoglobin-like motif, display characteristic heme-protein absorption spectra, and bind oxygen reversibly. Recombinant HemAT-Hs and HemAT-Bs shorter than 195 and 176 residues, respectively, do not bind heme effectively. Sequence homology comparisons and three-dimensional modeling predict that His-123 is the proximal heme-binding residue in HemAT from both species. The work described here used site-specific mutagenesis and spectroscopy to confirm this prediction, thereby providing direct evidence for a functional domain of prokaryotic signal transducers that bind heme in a globin fold. We postulate that this domain is part of a globin-coupled sensor (GCS) motif that exists as a two-domain transducer having no similarity to the PER-ARNT-SIM (PAS)-domain superfamily transducers. Using the GCS motif, we have identified several two-domain sensors in a variety of prokaryotes. We have cloned, expressed, and purified two potential globin-coupled sensors and performed spectral analysis on them. Both bind heme and show myoglobin-like spectra. This observation suggests that the general function of GCS-type transducers is to bind diatomic oxygen and perhaps other gaseous ligands, and to transmit a conformational signal through a linked signaling domain.

Two-domain hemoglobin gene of the water flea Moina macrocopa: duplication in the ancestral Cladocera, diversification, and loss of a bridge intron

Two cDNAs encoding the two-domain hemoglobin (Hb) chains of a crustacean Cladocera, Moina macrocopa, were cloned and their nucleotide (nt) sequences were determined. The amino acid (aa) sequences of both the gene products deduced from the nt sequences consisted of 348 residues and showed 98% identity with each other. These sequences together with the NH(2)-terminal aa sequences of the Hb chains determined after separation by two-dimensional gel electrophoresis showed that the Hb chains are synthesized as a secretory precursor with a signal peptide of 17 aa residues. The aa sequences of M. macrocopa Hb chains shared the following features with those of Daphnia Hb chains. Firstly, the signal peptide is followed by an NH(2)-terminal extension containing a threonine-rich sequence that might play a role in the multimerization of subunit chains. Secondly, the identity between the aa sequences of the first and second domains is exceptionally low. These facts suggest that duplication of the cladoceran Hb gene occurred before the divergence of families Moinidae and Daphniidae. Analysis of genomic DNA showed that the M. macrocopa Hb genes consist of two large repeated regions, encoding the first and second domains of Hb chains, respectively. The intron-exon organization of the first region of the M. macrocopa Hb genes was similar to that found in the Daphnia Hb genes, having the three-exon, two-intron structure characteristic of animal Hb genes. However, the intron bridging the two regions and the most downstream intron in the second region were missing in the Moina genes, providing a new example of intron loss. The following elements in the 5' flanking region were conserved in the Moina and Daphnia genes: (1) TATAAA, a typical TATA box sequence accompanied by a downstream sequence, GAAXAGCATCAGTT (the fourth residue X was G or A in Daphnia and absent in Moina); (2) CCAAT boxes, located upstream of the TATA box; (3) the binding sites for HIF-1 and GATA-1, also located upstream of the TATA box, that may be responsible for up-regulation of the cladoceran Hb genes under hypoxia.

The primary structure of globin and linker chains from the chlorocruorin of the polychaete Sabella spallanzanii

Annelid hemoglobins are organized in a very complex supramolecular network of interacting polypeptides, the structure of which is still not wholly resolved. We have separated by two-dimensional electrophoresis the 4-MDa chlorocruorin of Sabella spallanzanii and identified its components by amino-terminal sequencing. This work reveals a high rate of heterogeneity of constituent chains in a single animal as well as in the Sabella population. Using a cDNA library prepared from the hematopoietic tissue of this worm, we have isolated and fully sequenced most globin and linker cDNAs. The primary structure features of these polypeptides have been characterized by comparison with model globin and linker sequences.

The evolution of extracellular hemoglobins of annelids, vestimentiferans, and pogonophorans

The evolution of extracellular hemoglobins of annelids, vestimentiferans, and pogonophorans was investigated by applying cladistic and distance-based approaches to reconstruct the phylogenetic relationships of this group of respiratory pigments. We performed this study using the aligned sequences of globin and linker chains that are the constituents of these complex molecules. Three novel globin and two novel linker chains of Sabella spallanzanii described in an accompanying paper (Pallavicini, A., Negrisolo, E., Barbato, R., Dewilde, S., Ghiretti-Magaldi, A., Moens, L., and Lanfranchi, G. (2001) J. Biol. Chem. 276, 26384--26390) were also included. Our results allowed us to test previous hypotheses on the evolutionary pathways of these proteins and to formulate a new most parsimonious model of molecular evolution. According to this novel model, the genes coding for the polypeptides forming these composite molecules were already present in the common ancestor of annelids, vestimentiferans, and pogonophorans.

Very high resolution structure of a trematode hemoglobin displaying a TyrB10-TyrE7 heme distal residue pair and high oxygen affinity

Monomeric hemoglobin from the trematode Paramphistomum epiclitum displays very high oxygen affinity (P(50)<0.001 mm Hg) and an unusual heme distal site containing tyrosyl residues at the B10 and E7 positions. The crystal structure of aquo-met P. epiclitum hemoglobin, solved at 1.17 A resolution via multiwavelength anomalous dispersion techniques (R-factor=0.121), shows that the heme distal site pocket residue TyrB10 is engaged in hydrogen bonding to the iron-bound ligand. By contrast, residue TyrE7 is unexpectedly locked next to the CD globin region, in a conformation unsuitable for heme-bound ligand stabilisation. Such structural organization of the E7 distal residue differs strikingly from that observed in the nematode Ascaris suum hemoglobin (bearing TyrB10 and GlnE7 residues), which also displays very high oxygen affinity. The oxygenation and carbonylation parameters of wild-type P. epiclitum Hb as well as of single- and double-site mutants, with residue substitutions at positions B10, E7 and E11, have been determined and are discussed here in the light of the protein atomic resolution crystal structure.

Multimeric hemoglobin of the Australian brine shrimp Parartemia

The hemoglobin molecule of the commercially important brine shrimp Artemia sp. has been used extensively as a model for the study of molecular evolution. It consists of nine globin domains joined by short linker sequences, and these domains are believed to have originated through a series of duplications from an original globin gene. In addition, in Artemia, two different polymers of hemoglobin, called C and T, are found which differ by 11.7% at the amino acid level and are believed to have diverged about 60 MYA. This provides a set of data of 18 globin domain sequences that have evolved in the same organism. The pattern of amino acid substitution between these two polymers is unusual, with pairs of equivalent domains displaying differences of up to 2.7-fold in total amino acid substitution. Such differences would reflect a similar range of molecular-clock rates in what appear to be duplicate, structurally equivalent domains. In order to provide a reference outgroup, we sequenced the cDNA for a nine-domain hemoglobin (P) from another genus of brine shrimp, Parartemia zietziana, which differs morphologically and ecologically from Artemia and is endemic to Australia. Parartemia produces only one hundredth the amount of hemoglobin that Artemia produces and does not upregulate production in response to low oxygen partial pressure. Comparison of the globin domains at the amino acid and DNA levels suggests that the Artemia globin T gene has accumulated substitutions differently from the Parartemia P and Artemia C globin genes. We discuss the questions of accelerated evolution after duplication and possible functions for the Parartemia globin.

Nonvertebrate hemoglobins: functions and molecular adaptations

Hemoglobin (Hb) occurs in all the kingdoms of living organisms. Its distribution is episodic among the nonvertebrate groups in contrast to vertebrates. Nonvertebrate Hbs range from single-chain globins found in bacteria, algae, protozoa, and plants to large, multisubunit, multidomain Hbs found in nematodes, molluscs and crustaceans, and the giant annelid and vestimentiferan Hbs comprised of globin and nonglobin subunits. Chimeric hemoglobins have been found recently in bacteria and fungi. Hb occurs intracellularly in specific tissues and in circulating red blood cells (RBCs) and freely dissolved in various body fluids. In addition to transporting and storing O(2) and facilitating its diffusion, several novel Hb functions have emerged, including control of nitric oxide (NO) levels in microorganisms, use of NO to control the level of O(2) in nematodes, binding and transport of sulfide in endosymbiont-harboring species and protection against sulfide, scavenging of O(2 )in symbiotic leguminous plants, O(2 )sensing in bacteria and archaebacteria, and dehaloperoxidase activity useful in detoxification of chlorinated materials. This review focuses on the extensive variation in the functional properties of nonvertebrate Hbs, their O(2 )binding affinities, their homotropic interactions (cooperativity), and the sensitivities of these parameters to temperature and heterotropic effectors such as protons and cations. Whenever possible, it attempts to relate the ligand binding properties to the known molecular structures. The divergent and convergent evolutionary trends evident in the structures and functions of nonvertebrate Hbs appear to be adaptive in extending the inhabitable environment available to Hb-containing organisms.

A novel two-over-two alpha-helical sandwich fold is characteristic of the truncated hemoglobin family

Small hemoproteins displaying amino acid sequences 20-40 residues shorter than (non-)vertebrate hemoglobins (Hbs) have recently been identified in several pathogenic and non-pathogenic unicellular organisms, and named 'truncated hemoglobins' (trHbs). They have been proposed to be involved not only in oxygen transport but also in other biological functions, such as protection against reactive nitrogen species, photosynthesis or to act as terminal oxidases. Crystal structures of trHbs from the ciliated protozoan Paramecium caudatum and the green unicellular alga Chlamydomonas eugametos show that the tertiary structure of both proteins is based on a 'two-over-two' alpha-helical sandwich, reflecting an unprecedented editing of the classical 'three-over-three' alpha-helical globin fold. Based on specific Gly-Gly motifs the tertiary structure accommodates the deletion of the N-terminal A-helix and replacement of the crucial heme-binding F-helix with an extended polypeptide loop. Additionally, concerted structural modifications allow burying of the heme group and define the distal site, which hosts a TyrB10, GlnE7 residue pair. A set of structural and amino acid sequence consensus rules for stabilizing the fold and the bound heme in the trHbs homology subfamily is deduced.

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.

The N-terminal region of the heme-regulated eIF2alpha kinase is an autonomous heme binding domain

The N-terminal domain (NTD) of the heme-regulated eukaryotic initiation factor (eIF)2alpha kinase (HRI) was aligned to sequences in the NCBI data base using ENTREZ and a PAM250 matrix. Significant similarity was found between amino acids 11-118 in the NTD of rabbit HRI and amino acids 16-120 in mammalian alpha-globins. Several conserved amino acid residues present in globins are conserved in the NTD of HRI. His83 of HRI was predicted to be equivalent to the proximal heme ligand (HisF8) that is conserved in all globins. Molecular modeling of the NTD indicated that its amino acid sequence was compatible with the globin fold. Recombinant NTD (residues 1-159) was expressed in Escherichia coli. Spectral analysis of affinity purified recombinant NTD indicated that the NTD contained stably bound hemin. Mutational analysis indicated that His83 played a critical structural role in the stable binding of heme to the NTD, and was required to stabilize full length HRI synthesized de novo in the rabbit reticulocyte lysate. These results indicate that the NTD of HRI is an autonomous heme-binding domain, with His83 possibly serving as the proximal heme binding ligand.

1H NMR investigation of the distal hydrogen bonding network and ligand tilt in the cyanomet complex of oxygen-avid Ascaris suum hemoglobin

The O(2)-avid hemoglobin from the parasitic nematode Ascaris suum exhibits one of the slowest known O(2) off rates. Solution (1)H NMR has been used to investigate the electronic and molecular structural properties of the active site for the cyano-met derivative of the recombinant first domain of this protein. Assignment of the heme, axial His, and majority of the residues in contact with the heme reveals a molecular structure that is the same as reported in the A. suum HbO(2) crystal structure (Yang, J., Kloek, A., Goldberg, D. E., and Mathews, F. S. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 4224-4228) with the exception that the heme in solution is rotated by 180 degrees about the alpha,gamma-meso axis relative to that in the crystal. The observed dipolar shifts, together with the crystal coordinates of HbO(2), provide the orientation of the magnetic axes in the molecular framework. The major magnetic axis, which correlates with the Fe-CN vector, is found oriented approximately 30 degrees away from the heme normal and indicates significant steric tilt because of interaction with Tyr(30)(B10). The three side chain labile protons for the distal residues Tyr(30)(B10) and Gln(64)(E7) were identified, and their relaxation, dipolar shifts, and nuclear Overhauser effects to adjacent residues used to place them in the distal pocket. It is shown that these two distal residues exhibit the same orientations ideal for H bonding to the ligand and to each other, as found in the A. suum HbO(2) crystal. It is concluded that the ligated cyanide participates in the same distal H bonding network as ligated O(2). The combination of the strong steric tilt of the bound cyanide and slow ring reorientation of the Tyr(30)(B10) side chain supports a crowded and constrained distal pocket.

Iterated sequence databank search methods

Iterated sequence databank search methods were assessed from the viewpoint of someone with the sequence of a novel gene product wishing to find distant relatives to their protein and, with the specific searches against the PDB, also hoping to find a relative of known structure. We examined three methods in detail, spanning a range from simple pattern-matching to sophisticated weighted profiles. Rather than apply these methods 'blindly' (with default parameters) to a large number of test queries, we have concentrated on the globins, so allowing a more detailed investigation of each method on different data subsets with different parameter settings. Despite their widespread use, regular-expression matching proved to be very limited-seldom extending beyond the sub-family from which the pattern was derived. To attain any generality, the patterns had to be 'stripped-down' to include only the most highly conserved parts. The QUEST program avoided these problems by introducing a more flexible (weighted) matching. On the PDB sequences this was highly effective, missing only a few globins with probes based on each sub-family or even a single representative from each sub-family. In addition, very few false-positives were encountered, and those that did match, often only did so for a few cycles before being lost again. On the larger sequence collection, however, QUEST encountered problems with maintaining (or achieving) the alignment of the full globin family. psi-BLAST also recognised almost all the globins when matching against the PDB sequences, typically, missing three or four of the most distantly related sequences while picking-up a few false-positives. In contrast to QUEST, psi-BLAST performed very well on the larger databank, getting almost a full collection of globins although still retaining the same proportion of false-positives. SAM applied to the PDB sequences performed reasonably well with the myoglobin and hemoglobin families as probes, missing, typically several of the more difficult proteins but performed poorly with the leghemoglobin probe. Only with the full family range as a probe did it produce results comparable to psi-BLAST and QUEST. With the larger databank, SAM produced a good result but, again, this was only achieved using the full range of sequence variation with the default regulariser and use of Dirichlet mixtures completely failed in this situation.

Heterogeneity and differential expression under hypoxia of two-domain hemoglobin chains in the water flea, Daphnia magna

Hemoglobin (Hb) purified from the water flea, Daphnia magna, reared under hypoxia was analyzed by two-dimensional gel electrophoresis. The Hb was shown to be composed of six major subunit chain species (designated as DHbA to DHbF). The NH2-terminal amino acid sequences of DHbA, DHbB, DHbC, and DHbF are different from one another, indicating that at least four Hb genes are present in D. magna. The NH2-terminal amino acid sequences of DHbD and DHbE are the same as those of DHbA and DHbB, respectively. The six Hb chains were also found in the animal reared under normoxia in small amounts and with altered composition; the extent of decrease under normoxia was higher in the amounts of DHbC, DHbD, and DHbF than those of others. These results indicate that the Hb genes are differentially regulated by the ambient oxygen concentration. Four Hb genes constituting a cluster in the order, dhb4, dhb3, dhb1, and dhb2, were found on the chromosome of D. magna. The complete nucleotide sequences of the dhb1, dhb2, and dhb3 genes and their cDNAs showed that the genes have a seven-exon, six-intron structure. The structure consists of an intron separating an exon encoding a secretory signal sequence, two large repeated regions of a three-exon, two-intron structure that encode each a domain containing a heme-binding site, and an intron bridging the two repeated regions. The deduced amino acid sequences of the gene products showed higher than 79% identity to one another and showed unique features conserved in D. magna Hb chains. The analysis also suggested that DHbB (or DHbE), DHbF, and DHbC are encoded by the dhb1, dhb2, and dhb3 genes, respectively.

Structural, functional, and genetic characterization of Gastrophilus hemoglobin

Hemoglobin of Gastrophilus intestinalis (Insecta, Diptera), was purified and characterized. At least two isoforms have been identified by isoelectrofocusing, mass spectrometry, and genomic Southern blotting. Functional studies show a high oxygen affinity due to a low ligand dissociation rate (koff = 2.4 s-1) and a relatively high autoxidation rate (t1/2 = 1.6/h). The globins were separated under denaturing conditions, and the sequence of Hb1 (Mr = 17,965 +/- 2) was determined at the protein and DNA level. The open reading frame codes for a polypeptide of 150 amino acids. Although the globin is distantly related to globins from other species, it has a low penalty score against globin templates. Freshly isolated hemoglobin was crystallized from polyethylene glycol. Crystals contain two hemoglobin molecules per asymmetric unit. Solution of the three-dimensional structure by molecular replacement could not be achieved, possibly due to the presence of three protein isoforms in the crystals. In order to determine its three-dimensional structure, G. intestinalis Hb1 was overexpressed in Escherichia coli, resulting in a fully functional molecule as confirmed by ligand binding affinity. The globin gene contains two introns at positions D7.0 and G7.0. The D7.0 intron is unprecedented, suggesting that globin gene evolution is much more complex than originally thought.