Glycera dibranchiata hemoglobin. Structure and refinement at 1.5 A resolution

Purification of diverse hemoglobins by metal salt precipitation

Although donated blood is the preferred material for transfusion, its limited availability and stringent storage requirements have motivated the development of blood substitutes. The giant extracellular hemoglobin (aka erythrocruorin) of the earthworm Lumbricus terrestris (LtEc) has shown promise as a blood substitute, but an efficient purification method for LtEc must be developed to meet the potential large demand for blood substitutes. In this work, an optimized purification process that uses divalent and trivalent metal salts to selectively precipitate human, earthworm, and bloodworm hemoglobin (HbA, LtEc, and GdHb, respectively) from crude solutions was developed. Although several metal ions were able to selectively precipitate LtEc, Zn(2+) and Ni(2+) provided the lowest heme oxidation and highest overall yield of LtEc. In contrast, Zn(2+) was the only metal ion that completely precipitated HbA and GdHb. Polyacrylamide gel electrophoresis (PAGE) analysis shows that metal precipitation removes several impurities to provide highly pure hemoglobin samples. Heme oxidation levels were relatively low for Zn(2+)-purified HbA and LtEc (2.4±1.3% and 5.3±2.1%, respectively), but slightly higher for Ni(2+)-purified LtEc (8.4±1.2%). The oxygen affinity and cooperativity of the precipitated samples are also identical to samples purified with tangential flow filtration (TFF) alone, indicating the metal precipitation does not significantly affect the function of the hemoglobins. Overall, these results show that hemoglobins from several different species can be highly purified using a combination of metal (Zn(2+)) precipitation and tangential flow filtration.

Orthogonal Multipolar Interactions in Structural Chemistry and Biology

The past few decades of molecular recognition studies have greatly enhanced our knowledge on apolar, ion-dipole, and hydrogen-bonding interactions. However, much less attention has been given to the role that multipolar interactions, in particular those with orthogonal dipolar alignment, play in organizing a crystal lattice or stabilizing complexes involving biological receptors. By using results from database mining, this review attempts to give an overview of types and structural features of these previously rather overlooked interactions. A number of illustrative examples of these interactions found in X-ray crystal structures of small molecules and protein-ligand complexes demonstrate their propensity and thus potential importance for both, chemical and biological molecular recognition processes.

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.

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.

Protein structure alignment using a genetic algorithm

We have developed a novel, fully automatic method for aligning the three-dimensional structures of two proteins. The basic approach is to first align the proteins' secondary structure elements and then extend the alignment to include any equivalent residues found in loops or turns. The initial secondary structure element alignment is determined by a genetic algorithm. After refinement of the secondary structure element alignment, the protein backbones are superposed and a search is performed to identify any additional equivalent residues in a convergent process. Alignments are evaluated using intramolecular distance matrices. Alignments can be performed with or without sequential connectivity constraints. We have applied the method to proteins from several well-studied families: globins, immunoglobulins, serine proteases, dihydrofolate reductases, and DNA methyltransferases. Agreement with manually curated alignments is excellent. A web-based server and additional supporting information are available at http://engpub1.bu.edu/-josephs.

Non-functional conserved residues in globins and their possible role as a folding nucleus

Structure-based sequence alignment of 728 sequences of different globin subfamilies shows that in each subfamily there are two clusters of consensually conserved residues. The first is the well-known "functional" cluster which includes six heme-binding conserved residues (Phe CD1, His F8; aliphatic E11, FG5; hydrophobic F4, G5) and seven other conserved residues (Pro C2; aliphatic H19; hydrophobic B10, B13, B14, CD4, E4) that do not bind the heme but belong to its immediate neighborhood. The second cluster revealed here (aliphatic A8, G16, G12; aromatic A12; hydrophobic H8 and possibly H12) is distant from the heme. It is entirely non-polar and includes one turn (i, i+4 positions) from each of helices A, G, and H. It is known that A, G, and H helices formed at the earliest stage of apomyoglobin folding remain relatively stable in the equilibrium molten globule state, and are likely to be tightly packed with each other in this state. We have shown the existence of two similar conserved clusters in c -type cytochromes, heme-binding and distal from the heme. The second cluster in c -cytochromes includes one turn from each of the N and C-terminal alpha-helices. These N and C-terminal helices in cytochrome c are formed at the earliest stage of protein folding, remain relatively stable in the molten globule state, and are tightly packed with each other in this state, similar to the observed behavior of the globins. At least these two large protein families (c -type cytochromes and globins) have a close similarity in the existence and mutual positions of non-functional conserved residues. We assume that non-functional conserved residues are requisite for the fast and correct folding of both of these protein families into their stable 3D structures.

Structure of a biological oxygen sensor: a new mechanism for heme-driven signal transduction

The FixL proteins are biological oxygen sensors that restrict the expression of specific genes to hypoxic conditions. FixL's oxygen-detecting domain is a heme binding region that controls the activity of an attached histidine kinase. The FixL switch is regulated by binding of oxygen and other strong-field ligands. In the absence of bound ligand, the heme domain permits kinase activity. In the presence of bound ligand, this domain turns off kinase activity. Comparison of the structures of two forms of the Bradyrhizobium japonicum FixL heme domain, one in the "on" state without bound ligand and one in the "off" state with bound cyanide, reveals a mechanism of regulation by a heme that is distinct from the classical hemoglobin models. The close structural resemblance of the FixL heme domain to the photoactive yellow protein confirms the existence of a PAS structural motif but reveals the presence of an alternative regulatory gateway.

Solution structure and backbone dynamics of component IV Glycera dibranchiata monomeric hemoglobin-CO

The solution structure and backbone dynamics of the recombinant, ferrous CO-ligated form of component IV monomeric hemoglobin from Glycera dibranchiata (GMH4CO) have been characterized by NMR spectroscopy. Distance geometry and simulated annealing calculations utilizing a total of 2550 distance and torsion angle constraints yielded an ensemble of 29 structures with an overall average backbone rmsd of 0.48 A from the average structure. Differences between the solution structure and a related crystal structure are confined to regions of lower precision in either the NMR or X-ray structure, or in regions where the amino acid sequences differ. 15N relaxation measurements at 76.0 and 60.8 MHz were analyzed with an extended model-free approach, and revealed low-frequency motions in the vicinity of the heme, concentrated in the F helix. Amide proton protection factors were obtained from H-D amide exchange measurements on 15N-labeled protein. Patterns in the backbone dynamics and protection factors were shown to correlate with regions of heterogeneity and disorder in the ensemble of NMR structures and with large crystallographic B-factors in the X-ray structures. Surprisingly, while the backbone atoms of the F helix have higher rmsds and larger measures of dynamics on the microsecond to millisecond time scale than the other helices, amide protection factors for residues in the F helix were observed to be similar to those of the other helices. This contrasts with H-D amide exchange measurements on sperm whale myoglobin which indicated low protection for the F helix (S. N. Loh and B. F. Volkman, unpublished results). These results for GMH4 suggest a model in which the F helix undergoes collective motions as a relatively rigid hydrogen-bonded unit, possibly pivoting about a central position near residue Val87.

The mini-hemoglobins in neural and body wall tissue of the nemertean worm, Cerebratulus lacteus

Hemoglobin (Hb) occurs in circulating red blood cells, neural tissue, and body wall muscle tissue of the nemertean worm, Cerebratulus lacteus. The neural and body wall tissue each express single major Hb components for which the amino acid sequences have been deduced from cDNA and genomic DNA. These 109-residue globins form the smallest stable Hbs known. The globin genes have three exons and two introns with splice sites in the highly conserved positions of most globin genes. Alignment of the sequences with those of other globins indicates that the A, B, and H helices are about one-half the typical length. Phylogenetic analysis indicates that shortening results in a small tendency of globins to group together regardless of their actual relationships. The neural and body wall Hbs in situ are half-saturated with O2 at 2.9 and 4.1 torr, respectively. The Hill coefficient for the neural Hb in situ, approximately 2.9, suggests that the neural Hb self-associates in the deoxy state at least to tetramers at the 2-3 mM (heme) concentration estimated in the cells. The Hb must dissociate upon oxygenation and dilution because the weight-average molecular mass of the HbO2 in vitro is only about 18 kDa at 2-3 microM heme concentration. Calculations suggest that the Hb can function as an O2 store capable of extending neuronal activity in an anoxic environment for 5-30 min.

Solution and crystal structures of a sperm whale myoglobin triple mutant that mimics the sulfide-binding hemoglobin from Lucina pectinata

The bivalve mollusc Lucina pectinata harbors sulfide-oxidizing chemoautotrophic bacteria and expresses a monomeric hemoglobin I, HbI, with normal O2, but extraordinarily high sulfide affinity. The crystal structure of aquomet Lucina HbI has revealed an active site with three residues not commonly found in vertebrate globins: Phe(B10), Gln(E7), and Phe(E11) (Rizzi, M., Wittenberg, J. B., Coda, A., Fasano, M., Ascenzi, P., and Bolognesi, M. (1994) J. Mol. Biol. 244, 86-89). Engineering these three residues into sperm whale myoglobin results in a triple mutant with approximately 700-fold higher sulfide affinity than for wild-type. The single crystal x-ray structure of the aquomet derivative of the myoglobin triple mutant and the solution 1H NMR active site structures of the cyanomet derivatives of both the myoglobin mutant and Lucina HbI have been determined to examine further the structural origin of their unusually high sulfide affinities. The major differences in the distal pocket is that in the aquomet form the carbonyl of Gln64(E7) serves as a H-bond acceptor, whereas in the cyanomet form the amido group acts as H-bond donor to the bound ligand. Phe68(E11) is rotated approximately 90 degrees about chi2 and located approximately 1-2 A closer to the iron atom in the myoglobin triple mutant relative to its conformation in Lucina HbI. The change in orientation potentially eliminates the stabilizing interaction with sulfide and, together with the decrease in size of the distal pocket, accounts for the 7-fold lower sulfide affinity of the myoglobin mutant compared with that of Lucina HbI.

Cooperative mechanism in the homodimeric myoglobin from Nassa mutabilis

Oxygen binding and spectroscopic properties of the homodimeric myoglobin (Mb) from the prosobranchia sea snail Nassa mutabilis have been investigated. Oxygen equilibrium curves are pH-independent and cooperative with P50 = 5 +/- 1 mmHg and n approximately 1.5. Circular dichroism spectra of the oxygenated and deoxygenated form of N. mutabilis Mb are superimposable between 190 and 250 nm, suggesting a mechanism for cooperative ligand binding that does not involve changes in the alpha-helical content of the whole protein. The oxygen dissociation process is biphasic and pH-dependent, with different pKa values (=6.7 +/- 0.2 and 8.5 +/- 0.3) for the two phases. Moreover, the activation energy is essentially the same for both oxygen dissociation processes (Ea = 56.4 +/- 2.1 kJ/mol for the fast phase, and Ea = 53.8 +/- 1.9 kJ/mol for the slow phase), indicating that the rate difference for O2 dissociation between the diliganded and the monoliganded species is mostly dependent on a variation of the activation entropy. Ferrous nitrosylated N. mutabilis Mb shows, at alkaline and neutral pH, axial and rhombic X-band EPR signals, respectively, which display below pH 6 a three-hyperfine pattern typical of five-coordination. The results presented here suggest that in N.mutabilis Mb the kinetic control of cooperativity operates through a mechanism never observed before in other hemoproteins, which requires a ligand-linked large enhancement for the value of the oxygen association process in a molecule not undergoing changes in quaternary structure.

Alignment algorithm for homology modeling and threading

A DNA/protein sequence comparison is a popular computational tool for molecular biologists. Finding a good alignment implies an evolutionary and/or functional relationship between proteins or genomic loci. Sequential similarity between two proteins indicates their structural resemblance, providing a practical approach for structural modeling, when structure of one of these proteins is known. The first step in the homology modeling is a construction of an accurate sequence alignment. The commonly used alignment algorithms do not provide an adequate treatment of the structurally mismatched residues in locally dissimilar regions. We propose a simple modification of the existing alignment algorithm which treats these regions properly and demonstrate how this modification improves sequence alignments in real proteins.

Structural features of Glycera dibranchiata monomer hemoglobins. Primary sequences of monomer hemoglobin components II and III

Primary sequences for the remaining two members (GMH2, GMH3) of the group of three major monomeric hemoglobins from the marine annelid Glycera dibranchiata have been obtained. Full sequences of each 147-amino acid globin were achieved with a high degree of confidence using standard Edman technology in combination with molecular mass determinations of the intact globins and of the cyanogen bromide cleavage fragments using electrospray ionization mass spectrometry. When minor assumptions concerning Q/E identities are made these new results indicate the likely correspondence of GMG2 with the protein represented by the first Glycera dibranchiata monomer hemoglobin complete sequence [Imamura et al., (1972), J. Biol. Chem. 247, 2785-2797]. When these new sequences are combined with the previously determined primary sequence for the third major monomer hemoglobin, GMH4 [Alam et al., J. Protein Chem. (1994), 13, 151-164], it becomes clear that these three (GMG2-4) are truly distinct proteins, contrary to previous suggestions. Surprisingly, our results show that none of these three primary sequences is identical to the published sequence of the refined monomer hemoglobin crystal structure protein; however, there is a strong correspondence to the GMG2 sequence. The present sequencing results, in combination with the published GMH4 sequence, confirm the presence of a distal Leu in place of the more commonly encountered distal His in all three of the major monomer hemoglobins isolated in this laboratory and indicate that the unusual B10 Phe occurs only in GMH4. Analysis of the sequences presented here, along with comparison of amino acid content for Glycera dibranchiata monomer hemoglobins isolated from three different laboratories, and comparison of NMR results from two laboratories suggest further correspondence which unify disparate published isolations.

Globin and globin gene structure of the nerve myoglobin of Aphrodite aculeata

The globin of the nerve cord of the polychaete annelid Aphrodite aculeata was isolated and purified to homogeneity. The native molecule has a pI of 6.3 and acts as a dimer of two identical Mr 15, 644.5 polypeptide chains as determined by electrospray mass spectrometry. It has an average affinity for oxygen (P50 = 1.24 torr) resulting from fast association (kon = 170 X 10(6) M-1 . s-1) and dissociation rates (koff = 360 s-1). The partial primary structure of this nerve globin was determined at the protein level and completed and confirmed by translation of the cDNA sequence. The globin chain has 150 amino acid residues and a calculated Mr of 15, 602.69 strongly suggesting that the amino terminus is acetylated. The absence of a leader sequence and the lack of Cys at the positions NA2 and H9 needed for the formation of the high Mr complexes found in extracellular annelid globins classify the Aphrodite globin with the cellular globin species. The Aphrodite nerve globin is unlikely to represent a separate globin family, as cDNA derived primers detect globin messenger RNA in muscle, gut, and pharynx tissue as well. The gene encoding this globin species is interrupted by a single intron, inserted at position G7.0. Comparison to other globin gene structures strongly suggest that introns can be lost independently, rather than simultaneously as a result of a single conversion event as suggested previously (Lewin, R. (1984) Science 226, 328).

Iron ligand recognition by monomeric hemoglobins

Binding affinities of monomeric Glycera dibranchiata hemoglobin for some anions and heterocyclic amines, including imidazoles, pyrazole, triazole and tetrazole have been evaluated and compared with those of sperm whale and horse heart myoglobin. The proteins' affinities for substituted heterocyclic amines are strongly influenced by the steric bulk and flexibility of the aromatic ring. The ligand coordination mode depends on the heme oxidation state, iron(III) amine adducts being more stable than the iron(II) adducts, the higher affinities of stronger Brønsted-Lowry bases reflecting their essentially sigma-donor character. The bifunctional molecule morpholinoethylisocyanide acts as a redox-state-dependent ambidentate ligand, binding as an N-donor to iron(III), but as a C-donor to iron(II). pH-Dependences of the ESR and optical spectra of the azole adducts reveal iron-linked ionisations and spin-equilibria in the heme pocket. Enthalpy and entropy changes for the binding process were estimated for several ligands, and mutually compensatory behaviour is observed globally for delta H degree and delta S degree. At the compensation temperature theta, the binding affinities of monomeric Glycera dibranchiata hemoglobin and sperm whale myoglobin are similar and associated with free energy changes delta G degree (theta) approximately -9 +/- 1 kJ mol-1 for the heterocyclic and anionic ligands.

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

Seven-hundred globin sequences, including 146 nonvertebrate sequences, were aligned on the basis of conservation of secondary structure and the avoidance of gap penalties. Of the 182 positions needed to accommodate all the globin sequences, only 84 are common to all, including the absolutely conserved PheCD1 and HisF8. The mean number of amino acid substitutions per position ranges from 8 to 13 for all globins and 5 to 9 for internal positions. Although the total sequence volumes have a variation approximately 2-3%, the variation in volume per position ranges from approximately 13% for the internal to approximately 21% for the surface positions. Plausible correlations exist between amino acid substitution and the variation in volume per position for the 84 common and the internal but not the surface positions. The amino acid substitution matrix derived from the 84 common positions was used to evaluate sequence similarity within the globins and between the globins and phycocyanins C and colicins A, via calculation of pairwise similarity scores. The scores for globin-globin comparisons over the 84 common positions overlap the globin-phycocyanin and globin-colicin scores, with the former being intermediate. For the subset of internal positions, overlap is minimal between the three groups of scores. These results imply a continuum of amino acid sequences able to assume the common three-on-three alpha-helical structure and suggest that the determinants of the latter include sites other than those inaccessible to solvent.

Quantification of tertiary structural conservation despite primary sequence drift in the globin fold

The globin family of protein structures was the first for which it was recognized that tertiary structure can be highly conserved even when primary sequences have diverged to a virtually undetectable level of similarity. This principle of structural inertia in molecular evolution is now evident for many other protein families. We have performed a systematic comparison of the sequences and structures of 6 representative hemoglobin subunits as diverse in origin as plants, clams, and humans. Our analysis is based on a 97-residue helical core in common to all 6 structures. Amino acid sequence identities range from 12.4% to 42.3% in pairwise comparisons, and, despite these variations, the maximal RMS deviation in alpha-carbon positions is 3.02 A. Overall, sequence similarity and structural deviation are significantly anticorrelated, with a correlation coefficient of -0.71, but for a set of structures having under 20% pairwise identity, this anticorrelation falls to -0.38, which emphasizes the weak connection between a specific sequence and the tertiary fold. There is substantial variability in structure outside the helical core, and functional characteristics of these globins also differ appreciably. Nevertheless, despite variations in detail that the sequence dissimilarities and functional differences imply, the core structures of these globins remain remarkably preserved.

Expression of recombinant monomer hemoglobins (component IV) from the marine annelid Glycera dibranchiata: evidence for primary sequence positional regulation of heme rotational disorder

A description of the efficient high-level expression of the monomer hemoglobin (GMG4) from Glycera dibranchiata is presented. The cDNA described by Simons and Satterlee [Simons, P.C., & Satterlee, J.D. (1989) Biochemistry 28, 8525-8530] was subcloned into an expression system, and conditions were found that led to the production of large amounts of soluble apoprotein (rec-gmg). These conditions included lowering the temperature during the induction period and growth in a rich medium with a higher ionic strength. Characterization of this reconstituted recombinant protein showed that it was not identical to the native GMH4 protein. Both UV-visible and 1H NMR data indicated differences within the holoprotein (rec-gmh) heme pocket compared to the native protein, the major difference being that two nonidentical heme orientations are significantly populated in rec-gmh. This phenomenon has been seen previously in other heme proteins, where these heme orientational isomers are described by a 180-deg rotation about the heme alpha-gamma meso axis. This work prompted the production of a complete chemical sequence for the native GMH4 [Alam S.L., Satterlee, J. D., & Edmonds, C. G. (1994) J. Protein Chem. 13, 151-164], which showed that the expressed rec-gmg protein differed at three primary sequence positions (41, 95, and 123) from the native component IV globin (GMG4). Subsequently, we have produced the triple-revertant mutations required to express the recombinant wild-type protein (recGMG4). The physical characteristics of the active site in the holoprotein (recGMH4) are identical to those of the native protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Volume changes in protein evolution

We have determined the variations in volume that occur during evolution in the buried core of three different families of proteins. The variation of the whole core is very small (approximately 2.5%) compared to the variation at individual sites (approximately 13%). However, by comparing our results to those expected from random sequences with no correlations between sites, we show that the small variation observed may simply be a manifestation of the statistical "law of large numbers" and not reflect any compensating changes in, or global constraints upon, protein sequences. We have also analysed in detail the volume variations at individual sites, both in the core and on the surface, and compared these variations with those expected from random sequences. Individual sites on the surface have nearly the same variation as random sequences (24% versus 28% variation). However, individual sites in the core have about half the variation of random sequences (13% versus 30%). Roughly, half of these core sites strongly conserve their volume (0 to 10% variation); one quarter have moderate variation (10 to 20%); and the remaining quarter vary randomly (20 to 40%). Our results have clear implications for the relationship between protein sequence and structure. For our analysis, we have developed a new and simple method for weighting protein sequences to correct for unequal representation, which we describe in an Appendix.

Complete amino acid sequence of the Glycera dibranchiata monomer hemoglobin component IV: structural implications

The globin derived from the monomer Component IV hemoglobin of the marine amnelid, Glycera dibranchiata, has been completely sequenced, and the resulting information has been used to create a structural model of the protein. The most important result is that the consensus sequence of Component IV differs by 3 amino acids from a cDNA-predicted amino acid sequence thought earlier to encode the Component IV hemoglobin. This work reveals that the histidine (E7), typical of most heme-containing globins, is replaced by leucine in Component IV. Also significant is that this sequence is not identical to any of the previously reported Glycera dibranchiata monomer hemoglobin sequences, including the sequence from a previously reported crystal structure, but has high identity to all. A three-dimensional structural model for monomer Component IV hemoglobin was constructed using the published 1.5 A crystal structure of a monomer hemoglobin from Glycera dibranchiata as a template. The model shows several interesting features: (1) a Phe31 (B10) that is positioned in the active site; (2) a His39 occurs in an interhelical region occupied by Pro in 98.2% of reported globin sequences; and (3) a Met41 is found at a position that emerges from this work as a previously unrecognized heme contact.

Role of hydrogen bonding to bound dioxygen in soybean leghemoglobin

Electron spin echo envelope modulation (ESEEM) spectroscopy was applied to oxy cobaltous soybean leghemoglobin (oxyCoLb) in D2O at various pH values to investigate electron nuclear superhyperfine coupling to N epsilon of the proximal histidyl imidazole and to exchangeable deuterons. Two spectroscopically distinct forms of oxyCoLb, acid and neutral, were identified. In the acid form, a 0.82-MHz hyperfine coupling to 2H was found, indicating the presence of a hydrogen bond to bound O2. No hyperfine-coupled 2H was found in the neutral form. Nuclear hyperfine and nuclear quadrupole couplings to the proximal histidyl N epsilon in the acid form are smaller than those in the neutral form: Aiso = 2.22 MHz and e2qQ = 1.98 MHz for the acid form; Aiso = 2.90 MHz and e2qQ = 2.22 MHz for the neutral form. The differences are believed to result from the presence of a hydrogen bond to bound O2 in the acid form. A discussion of the contribution of this hydrogen bond to the pH-dependent O2 affinity of leghemoglobin is presented.

A heme c-peptide model system for the resonance Raman study of c-type cytochromes: characterization of the solvent-dependence of peptide-histidine-heme interactions

The visible absorption and Soret-excited resonance Raman spectra of ferrous microperoxidase-8 [MP8(II)], an octapeptide containing a heme c, are reported. These spectroscopies indicate that MP8(II), dissolved in aqueous buffered solutions, forms low-spin six-coordinated complexes in the 7-14 pH range. Intermolecular bonding interactions of MP8(II) in water account for this behavior. On the contrary, when the hemopeptide is dispersed in aqueous solutions containing detergent or an alcohol, the spectroscopic data show that the iron atom of MP8(II) is essentially high-spin five-coordinated in accordance with a monomeric structure of MP8(II). In addition to a high-spin signature to the heme skeletal modes, the high-frequency regions of resonance Raman spectra characterize an electronic influence of the thioether bridges on the frequency of stretching modes of C beta-C beta bonds (nu 2, nu 11, and nu 29). On the other hand, the low-frequency Raman spectra of monomeric MP8(II) at pH 7.5 present significant differences in the 150-250-cm-1 regions depending upon the solvent composition (pH, presence or absence of detergent, alcohol). These effects are attributed to frequency variations of the Fe-N(His)-involving mode which indicate changes in the H-bonding interactions of the axial His and therefore solvent-dependent changes of the octapeptide conformation. Our resonance Raman data further show that the axial His of monomeric MP8(II) could be totally deprotonated in aqueous cetyltrimethylammonium bromide solution at very alkaline pH (pKa = 13.3). The vibrational data (100-1700 cm-1) obtained for the various monomeric forms of MP8(II) are expected to be useful for determining the heme structure and environment in reduced c'-type cytochromes. Comparisons of resonance Raman data with X-ray crystallographic data available for different hemoproteins allow us to evaluate the ionization and H-bonding states of the His bound to the high-spin five-coordinated hemes. These data are discussed in terms of proximal influence of protein-His-heme interactions on the determination and the regulation of a particular biological function.

Adventitious variability? The amino acid sequences of nonvertebrate globins

1. The more than 140 amino acid sequences of non-vertebrate hemoglobins (Hbs) and myoglobins (Mbs) that are known at present, can be divided into several distinct groups: (1) single-chain globins, containing one heme-binding domain; (2) truncated, single-chain, one-domain globins; (3) chimeric, one-domain globins; (4) chimeric, two-domain globins; and (5) chimeric multi-domain globins. 2. The crystal structures of eight nonvertebrate Hbs and Mbs are known, all of them monomeric, one-domain globin chains. Although these molecules represent plants, prokaryotes and several metazoan groups, and although the inter-subunit interactions in the dimeric and tetrameric molecules differ from the ones observed in vertebrate Hbs, the secondary structures of all seven one-domain globins retain the characteristic vertebrate "myoglobin fold". No crystal structures of globins representing the other four groups have been determined. 3. Furthermore, a number of the one-, two- and multi-domain globin chains participate in a broad variety of quaternary structures, ranging from homo- and heterodimers to highly complex, multisubunit aggregates with M(r) > 3000 kDa (S. N. Vinogradov, Comp. Biochem. Physiol. 82B, 1-15, 1985). 4. (1) The single-chain, single-domain globins are comparable in size to the vertebrate globins and exhibit the widest distribution. (A) Intracellular Hbs include: (i) the monomeric and polymeric Hbs of the polychaete Glycera; (ii) the tetrameric Hb of the echiuran Urechis; (iii) the dimeric Hbs of echinoderms such as Paracaudina and Caudina; and (iv) the dimeric and tetrameric Hbs of molluscs, the bivalves Scapharca, Anadara, Barbatia and Calyptogena. (B) Extracellular Hbs include: (i) the multiple monomeric and dimeric Hbs of the larva of the insect Chironomus; (ii) the Hbs of nematodes such as Trichostrongylus and Caenorhabditis; (iii) the globin chains forming tetramers and dodecamers and comprising approximately 2/3 of the giant (approximately 3600 kDa), hexagonal bilayer (HBL) Hbs of annelids, e.g. the oligochaete Lumbricus and the polychaete Tylorrhynchus and of the vestimentiferan Lamellibrachia; and (iv) the globin chains comprising the ca 400 kDa Hbs of Lamellibrachia and the pogonophoran Oligobrachia. (C) Cytoplasmic Hbs include: (i) the Mbs of molluscs, the gastropods Aplysia, Bursatella, Cerithedea, Nassa and Dolabella and the chiton Liolophura; (ii) the three Hb of the symbiont-harboring bivalve Lucina; (iii) the dimeric Hb of the bacterium Vitreoscilla; and (iv) plant Hbs, including the Hbs of symbiont-containing legumes (Lgbs), the Hbs of symbiont-containing non-leguminous plants and the Hbs in the roots of symbiont-free plants.(ABSTRACT TRUNCATED AT 400 WORDS)

The amino acid sequence of hemoglobin III from the symbiont-harboring clam Lucina pectinata

The cytoplasmic hemoglobin III from the gill of the symbiont-harboring clam Lucina pectinata consists of 152 amino acid residues, has a calculated Mm of 18,068, including heme, and has N-acetyl-serine as the N-terminal residue. Based on the alignment of its sequence with other vertebrate and nonvertebrate globins, it retains the invariant residues Phe45 at position CD1 and His98 at the proximal position F8, as well as the highly conserved Trp16 and Pro39 at positions A12 and C2, respectively. The most likely candidate for the distal residue at position E7 is Gln66. Lucina hemoglobin III shares 95 identical residues with hemoglobin II (J. D. Hockenhull-Johnson et al., J. Prot. Chem. 10, 609-622, 1991), including Tyr at position B10, which has been shown to be capable of entering the distal heme cavity and placing its hydroxyl group within a 2.8 A of the water molecule occupying the distal ligand position, by modeling the hemoglobin II sequence using the crystal structure of sperm whale metmyoglobin. The amino acid sequences of the two Lucina globins are compared in detail with the known sequences of mollusc globins, including seven cytoplasmic and 11 intracellular globins. Relative to 75% homology between the two Lucina globins (counting identical and conserved residues), both sequences have percent homology scores ranging from 36-49% when compared to the two groups of mollusc globins. The highest homology appears to exist between the Lucina globins and the cytoplasmic hemoglobin of Busycon canaliculatum.

Yeast flavohemoglobin is an ancient protein related to globins and a reductase family

The hemoglobin of yeast is a two-domain protein with both heme and flavin prosthetic groups. The nucleotide sequences of the cDNA and genomic DNA encoding the protein from Saccharomyces cerevisiae show that introns are absent and that both domains are homologous with a flavoheme protein from Escherichia coli. The heme domains are also homologous with those of O2-binding heme proteins from several other distantly related bacteria, plants, and animals; all appear to be members of the same globin superfamily. Although the homologous hemoglobin of the bacterium Vitreoscilla sp. is a single-domain protein, several bacteria have related O2-binding heme proteins whose second domains have different structures and enzymatic activities: dihydropteridine reductase (E. coli), cytochrome c reductase (Alcaligenes eutrophus), and kinase in the O2 sensor of Rhizobium meliloti. This indicates that one evolutionary pathway of hemoglobin is that of a multipurpose domain attached to a variety of unrelated proteins to form molecules with different functions. The flavin domain of yeast hemoglobin is homologous with members of a flavoprotein family that includes ferredoxin reductase, nitric oxide synthase, and cytochrome P-450 reductase. The correspondence of yeast and E. coli flavohemoglobins indicates that the two-domain protein has been conserved intact for as long as 1.8 billion years, the estimated time of divergence of prokaryotes and eukaryotes provided that cross-species gene transfer has not occurred.

Binding mode of azide to ferric Aplysia limacina myoglobin. Crystallographic analysis at 1.9 A resolution

The binding mode of azide to the ferric form of Aplysia limacina myoglobin has been studied by X-ray crystallography. The three-dimensional structure of the complex has been refined at 1.9 A resolution to a crystallographic R-factor of 13.9%, including 126 ordered solvent molecules. Azide binds to the heme iron, at the sixth co-ordination position, and is oriented towards the outer part of the distal site crevice. This orientation is stabilized by an ionic interaction with the side-chain of Arg66 (E10) which, from an outer orientation in the 'aquo-met' ligand-free myoglobin, folds back towards the distal site in the presence of the anionic ligand. In the absence of a hydrogen bond donor residue at the distal E7 position in Aplysia limacina myoglobin, a different polar residue, Arg66 at the E10 topological position, has been selected by molecular evolution in order to grant ligand stabilization.

Amino acid sequence of the monomer subunit of the giant extracellular hemoglobin of the aquatic oligochaete, Tubifex tubifex

The extracellular hemoglobin of the aquatic oligochaete Tubifex tubifex consists of four subunits: a monomer of 16.5 kDa, a disulfide-bonded trimer of about 50 kDa and at least two subunits of about 30 kDa. The complete amino acid sequence of the monomeric subunit was determined: it consists of 141 amino acid residues and has a molecular mass of 16,286 Da including a heme group. 39 residues (28%) were found to be identical with those in the corresponding positions in the monomeric globin chains from Lumbricus terrestris, Pheretima sieboldi, and Tylorrhynchus heterochaetus. Tubifex and Lumbricus are most similar, with 75 amino acid identities (53%). There are eight invariant residues amongst these monomeric globins and the intracellular monomeric globin of Glycera and the human beta-globin. The monomeric globin from Tubifex aligns best with those of group A, globins which have a Cys in their second position and an invariant Lys-Val-Lys at positions 9-11 [Gotoh et al. (1987) Biochem. J. 241, 441-445]. The two cysteine residues, at positions 2 and 131, appear to be disulfide-bonded.

The heterogeneity of the polymeric intracellular hemoglobin of Glycera dibranchiata and the cDNA-derived amino acid sequence of one component

The erythrocytes of the marine polychaete Glycera dibranchiata contain a number of different, single-chain hemoglobins, some of which self-associate into a 'polymeric' fraction. An oligodeoxynucleotide probe was synthesized based on partial amino acid sequences determined by chemical methods, and used to screen a cDNA library constructed from the poly(A+)mRNA of Glycera erythrocytes (Simons, P.C. and Satterlee, J.D. (1989) Biochemistry 28, 8525-8530). The longest positive inserts found were sequenced using the dideoxy nucleotide chain termination method. One complete clone was obtained: clone 5A, 816 bases long, contained 59 bases of 5'-untranslated RNA, an open reading frame of 441 bases coding for 147 amino acids and a 3'-untranslated region of 316 bases. The derived amino acid sequence of Glycera globin P1 was in agreement with the partial amino acid sequences obtained by chemical methods. Three additional inserts obtained in the screening were also sequenced: the inferred amino acid sequences proved to be partial globin sequences which were different from each other and from the sequence of P1. Thus, the 'polymeric' fraction of the intracellular hemoglobin of Glycera probably consists of at least four different globin chains much like the 'monomeric' fraction. Comparison of the 'polymeric' sequence with the two known 'monomeric' sequences, M-II and M-IV, shows that they share 54 identical residues. At 74 positions, the identical residues in M-II and M-IV differ from the corresponding residue in P1, including at E-7, where P1 has a distal His, in contrast to Leu in M-II and M-IV. The alignment of Bashford et al. ((1987) J. Mol. Biol. 196, 199-216) and their templates were used to examine the principal differences between the two types of Glycera globin sequences. They appear to consist of uncommon surface amino acid residues at positions C6 (Phe vs. Ala), E10 (Val vs. Lys), E17 (Lys vs. Val), G1 (Arg vs. Lys), G10 (Met vs. Ala) and H5 (Arg vs. Lys). One or more of these residues could be responsible for the self-association exhibited by the 'polymeric' Glycera globins.

Horse heart myoglobin reconstituted with a symmetrical heme. A circular dichroism study

Proton NMR studies on myoglobins and hemoglobins reconstituted with non-natural hemes, possessing different side chains in the pyrrolic rings, have provided interesting information for the understanding of the mechanism governing heme reorientation in the globin pocket, during synthesis of the native protein in vivo or in the reconstitution process in vitro. More recently, circular dichroism (CD) studies have been reported as a qualitative, alternative tool, with respect to 1H-NMR for detecting heme disorder in a reconstituted myoglobin or hemoglobin. In this paper, a CD study is reported on the reconstitution of horse heart myoglobin with protoheme XIII, a heme possessing true rotational symmetry about its alpha, gamma-meso axis. The results obtained show that the reconstitution product with this heme, which binds to the apoprotein with high affinity, not dissimilar from that of the natural heme, is characterized by a CD spectrum with bands possessing rotational strengths much lower than in the native protein. Furthermore, the CD changes detected as a function of time, during heme reorientation, in the case of natural heme, are absent when the apoprotein is reconstituted with protoheme XIII. These data provide independent evidence for reorientation of the natural heme, which follows its insertion into the protein matrix.

Comparison of the structures of globins and phycocyanins: evidence for evolutionary relationship

Globins and phycocyanins are two classes of proteins with different function, different ligands, and no substantial sequence similarity, yet the conformations of their polypeptide chains show very similar folding patterns. Does this arise from a genuine, albeit very distant, evolutionary relationship, or does it represent a common solution of a structural problem? We address this question by a very detailed comparison of the structures of the two protein families. An analysis of the helices and their interactions shows many features common to globins and phycocyanins, including some exceptional features of the globins such as a 3-10 C helix and the unusual "crossed-ridge" packing pattern at the B/E helix interfaces. We conclude that the evidence supports the hypothesis of distant evolutionary relationship between globins and phycocyanins.