Solution (1)H NMR study of the influence of distal hydrogen bonding and N terminus acetylation on the active site electronic and molecular structure of Aplysia limacina cyanomet myoglobin

The sea hare Aplysia limacina possesses a myoglobin in which a distal H-bond is provided by Arg E10 rather than the common His E7. Solution (1)H NMR studies of the cyanomet complexes of true wild-type (WT), recombinant wild-type (rWT), and the V(E7)H/R(E10)T and V(E7)H mutants of Aplysia Mb designed to mimic the mammalian Mb heme pocket reveal that the distal His in the mutants is rotated out of the heme pocket and is unable to provide a stabilizing H-bond to bound ligand and that WT and rWT differ both in the thermodynamics of heme orientational disorder and in heme contact shift pattern. The mean of the four heme methyl shifts is shown to serve as a sensitive indicator of variations in distal H-bonding among a set of mutant cyanomet globins. The heme pocket perturbations in rWT relative to WT were traced to the absence of the N-terminal acetyl group in rWT that participates in an H-bond to the EF corner in WT. Analysis of dipolar contacts between heme and axial His and between heme and the protein matrix reveal a small approximately 2 degrees rotation of the axial His in rWT relative to true WT and a approximately 3 degrees rotation of the heme in the double mutant relative to rWT Mb. It is demonstrated that both the direction and magnitude of the rotation of the axial His relative to the heme can be determined from the change in the pattern of the contact-dominated heme methyl shift and from the dipolar-dominated heme meso-H shift. However, only NOE data can determine whether it is the His or heme that actually rotates in the protein matrix.     

Structural and dynamic repercussions of heme binding and heme-protein cross-linking in Synechococcus sp. PCC 7002 hemoglobin

The recombinant two-on-two hemoglobin from the cyanobacterium Synechoccocus sp. PCC 7002 (S7002 rHb) is a bishistidine hexacoordinate globin capable of forming a covalent cross-link between a heme vinyl and a histidine in the C-terminal helix (H helix). Of the two heme axial histidines, His46 (in the E helix, distal side) and His70 (in the F helix, proximal histidine), His46 is displaced by exogenous ligands. S7002 rHb can be readily prepared as an apoglobin (apo-rHb), a non-cross-linked hemichrome (ferric iron and histidine axial ligands, rHb-R), and a cross-linked hemichrome (rHb-A). To determine the effects of heme binding and subsequent cross-linking, apo-rHb, rHb-R, and rHb-A were subjected to thermal denaturation and 1H/2H exchange. Interpretation of the latter data was based on nuclear magnetic resonance assignments obtained with uniformly 15N- and 13C,15N-labeled proteins. Apo-rHb was found to contain a cooperative structural core, which was extended and stabilized by heme binding. Cross-linking resulted in further stabilization attributed mainly to an unfolded-state effect. Protection factors were higher at the cross-link site and near His70 in rHb-A than in rHb-R. In contrast, other regions became less resistant to exchange in rHb-A. These included portions of the B and E helices, which undergo large conformational changes upon exogenous ligand binding. Thus, the cross-link readjusted the dynamic properties of the heme pocket. 1H/2H exchange data also revealed that the B, G, and H helices formed a robust core regardless of the presence of the heme or cross-link. This motif likely encompasses the early folding nucleus of two-on-two globins.     

Structural basis of human cytoglobin for ligand binding

Cytoglobin (Cgb), a newly discovered member of the vertebrate globin family, binds O(2) reversibly via its heme, as is the case for other mammalian globins (hemoglobin (Hb), myoglobin (Mb) and neuroglobin (Ngb)). While Cgb is expressed in various tissues, its physiological role is not clearly understood. Here, the X-ray crystal structure of wild type human Cgb in the ferric state at 2.4A resolution is reported. In the crystal structure, ferric Cgb is dimerized through two intermolecular disulfide bonds between Cys38(B2) and Cys83(E9), and the dimerization interface is similar to that of lamprey Hb and Ngb. The overall backbone structure of the Cgb monomer exhibits a traditional globin fold with a three-over-three alpha-helical sandwich, in which the arrangement of helices is basically the same among all globins studied to date. A detailed comparison reveals that the backbone structure of the CD corner to D helix region, the N terminus of the E-helix and the F-helix of Cgb resembles more closely those of pentacoordinated globins (Mb, lamprey Hb), rather than hexacoordinated globins (Ngb, rice Hb). However, the His81(E7) imidazole group coordinates directly to the heme iron as a sixth axial ligand to form a hexcoordinated heme, like Ngb and rice Hb. The position and orientation of the highly conserved residues in the heme pocket (Phe(CD1), Val(E11), distal His(E7) and proximal His(F8)) are similar to those of other globin proteins. Two alternative conformations of the Arg84(E10) guanidium group were observed, suggesting that it participates in ligand binding to Cgb, as is the case for Arg(E10) of Aplysia Mb and Lys(E10) of Ngb. The structural diversities and similarities among globin proteins are discussed with relevance to molecular evolutionary relationships.     

Invariant features of globin primary structure and coding of their secondary structure

Invariant features of the primary structure of 67 globins are analysed. These features may be responsible for the formation of the secondary structure of these proteins at the first stage of self-organization (in the unfolded chain). It is shown that in primary structures of globins there are 11 sites or regions of one to four residues in which at least one of the residues Asn, Asp, His, Pro, Ser or Thr is located in every globin (haem-linking His residues are excluded from these sites). An unambiguous correlation exists between the position of these regions and the secondary structure of globins: all these regions (except one) are located near the ends of helices in globins whose three-dimensional structure is known and the ends of all helices (except for the helix F) are coded by such regions. A decrease in the set of residues listed above leads to a sharp drop in the number of regions invariantly occupied by the residues, while an addition of residues such as Tyr and Gly to this set does not eventually increase the number of invariant regions. Five residues (Asn, Asp, His, Ser and Thr) of the six that code the ends of helical regions have polar side groups with a small number of degrees of freedom capable of forming hydrogen bonds with atoms of the backbone with a relatively small loss of entropy. One residue (Pro) has no NH-group and, therefore, has less chance of participating in the formation of hydrogen bonds between atoms of the backbone. This corroborates the hypothesis that competition between hydrogen bonds of short polar side groups and hydrogen bonds in the backbone is essential for the formation of the secondary structure in unfolded protein chains. Amino acid replacements in hydrophobic cores of the 67 globins are considered in the Appendix.     

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.     

Solution 1H NMR study of the active site molecular structure and magnetic properties of the cyanomet complex of the isolated, tetrameric beta-chain from human adult hemoglobin

The solution molecular structure and the electronic and magnetic properties of the heme pocket of the cyanomet complex of the isolated beta-chain of human adult hemoglobin, HbA, have been investigated by homonuclear 2D (1)H NMR in order to assess the extent of assignments allowed by (1)H NMR of a homo-tetrameric 65-kDa protein, to guide the future assignments of the heterotetrameric complex of HbA, and to compare the structure of the beta-chain to the crystallographically characterized complexes that contains the beta-chain. The target residues are those that exhibit significant (>|0.2| ppm) dipolar shifts, as predicted by a "preliminary" set of magnetic axes determined from a small set of easily assigned active site residues. All 104 target residues ( approximately 70% of total) were assigned by taking advantage of the temperature dependence predicted by the "preliminary" magnetic axes for the polypeptide backbone; they include all residues proposed to play a significant role in modulating the ligand affinity in the tetramer HbA. Left unassigned are the A-helix, the end of the G-helix and the beginning of the H-helix where dipolar shifts are less than |0.2| ppm. These comprehensive assignments allow the determination of a robust set of orientation and anisotropies of the paramagnetic susceptibility tensor that leads to quantitative interpretation of the dipolar shifts of the beta-chain in terms of the crystal coordinates of the beta-subunit in ligated HbA which, in turn, confirms a largely conserved molecular structure of the isolated beta-chain relative to that in the intact R-state HbA. The major magnetic axis, which is correlated with the tilt of the Fe-CN unit, is tilted approximately 10 degrees from the heme normal so that the Fe-CN unit is tilted toward the beta-meso-H in a fashion remarkably similar to the Fe-CO tilt in the beta-subunit of HbCO. It is concluded that a set of "preliminary" magnetic axes and the use of variable temperature 2D NMR spectra are crucial to effective assignments in the tetrameric cyanomet beta-chain and that this approach should be similarly effective in HbA.     

A comparison of the heme binding pocket in globins and cytochrome b5.

Of the 85 three-dimensionally characterized residues of cytochrome b5, 51 are found to be structurally and topologically equivalent to the globin fold. When these proteins have been superimposed, the heme irons are found to be less than 1.4 A separated and the heme normals are inclined by less than 9.5 degrees. The proximal histidine of the globins and two adjacent helices are equivalent to the sixth iron ligand and adjacent helices of cytochrome b5. Larger differences in structure are observed on the distal side of the heme, coincident with the most changeable part of the globin structures. The heme itself is rotated by 53 degrees about its normal but such a change is energetically minimal and conservative as the heme side groups are not directly involved in the function of the molecules. The beta-sheet of cytochrome b5 is inserted into a corresponding cavity of the globins forming an additional lining to the heme pocket. The roughly 50 residues missing at the carboxy end of the known cytochrome b5 fragment could correspond in part to the H helix in the globins. While it would seem probable that these similarities represent divergent evolution from a primordial heme-binding protein, the possibility of structural convergence to a functionally satisfactory protein cannot be excluded.     

Solution 1H NMR study of the active site molecular structure and magnetic properties of the cyanomet complex of the isolated alpha-chain from human hemoglobin A

The solution electronic and molecular structure for the heme pocket of the cyanomet complex of the isolated alpha-chain of human adult hemoglobin (HbA) has been investigated by homonuclear two-dimensional 1H NMR in order to establish an assignment protocol for the dimeric chain that will guide similar assignments in the intact, heterotetrameric HbA complex, and to compare the structures of the alpha-chain with its subunit in HbA. The target residues are those that exhibit significant (>0.2 ppm) dipolar shifts, as predicted by a "preliminary" set of magnetic axes determined from a small set of easily assigned active site residues. All 97 target residues (approximately 70% of total) were assigned by taking advantage of the temperature dependence predicted by the "preliminary" magnetic axes for the polypeptide backbone; they include all residues proposed to play a significant role in modulating the ligand affinity in the tetramer HbA. Left unassigned are the A-helix, the end of the G-helix and the beginning of the H-helix where dipolar shifts are less than 0.2 ppm. The complete assignments allow the determination of a robust set of orientation and anisotropies of the paramagnetic susceptibility tensor that leads to quantitative interpretation of the dipolar shifts of the alpha-chain in terms of the crystal coordinates of the alpha-subunit in ligated HbA which, in turn, confirms a largely conserved molecular structure of the isolated alpha-chain relative to that in the intact HbA. The major magnetic axis, which is correlated with the tilt of the Fe-CN unit, is tilted approximately 10 degrees from the heme normal so that the Fe-CN unit is tilted toward the beta-meso-H in a fashion remarkably similar to the Fe-CO tilt in HbACO. It is concluded that a set of "preliminary" magnetic axes and the use of variable temperature two-dimensional NMR spectra are crucial to effective assignments in the cyanomet alpha-chain and that this approach should be similarly effective in HbA.     

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.     

The polymeric hemoglobin molecule of Artemia. Interpretation of translated cDNA sequence of nine domains.

Translated cDNA for Artemia hemoglobin provided sequence data for almost nine domains, from the fourth residue of the A helix of one domain through 1405 residues to a stop codon after the ninth domain. The domain sequences were all different (homology between pairs 17-38%) but aligned well with each other and with conventional globins, satisfying the requirements for Phe at CD1, His at F8 and most other highly conserved features of globins including His at E7. Features found to be characteristic of Artemia globin and present in all nine domains were Phe at B10, Tyr at C4, Gly at F5, Phe at G5 and Gly at H22. Approximately 14 residues including a consensus -Val-Asp-Pro-Val-Thr-Gly-Leu- were available to form the linker between each pair of domains. The Artemia sequence data were compared with the crystal structures of Chironomus thummi thummi erythrocruorin III and sperm whale myoglobin in order to identify features of structural similarity and to examine the consequences of the differences. The Artemia sequences were compatible with the main helices and critical features of the globin fold. Possible modifications to the C helix, FG turn, and GH turn were studied in terms of molecular coordinates.     

The solution conformation of the antibacterial peptide cecropin A: a nuclear magnetic resonance and dynamical simulated annealing study

The solution conformation of the antibacterial polypeptide cecropin A from the Cecropia moth is investigated by nuclear magnetic resonance (NMR) spectroscopy under conditions where it adopts a fully ordered structure, as judged by previous circular dichroism studies [Steiner, H. (1982) FEBS Lett. 137, 283-287], namely, 15% (v/v) hexafluoroisopropyl alcohol. By use of a combination of two-dimensional NMR techniques the 1H NMR spectrum of cecropin A is completely assigned. A set of 243 approximate interproton distance restraints is derived from nuclear Overhauser enhancement (NOE) measurements. These, together with 32 distance restraints for the 16 intrahelical hydrogen bonds identified on the basis of the pattern of short-range NOEs, form the basis of a three-dimensional structure determination by dynamical simulated annealing [Nilges, M., Clore, G.M., & Gronenborn, A.M. (1988) FEBS Lett. 229, 317-324]. The calculations are carried out starting from three initial structures, an alpha-helix, an extended beta-strand, and a mixed alpha/beta structure. Seven independent structures are computed from each starting structure by using different random number seeds for the assignments of the initial velocities. All 21 calculated structures satisfy the experimental restraints, display very small deviations from idealized covalent geometry, and possess good nonbonded contacts. Analysis of the 21 converged structure indicates that there are two helical regions extending from residues 5 to 21 and from residues 24 to 37 which are very well defined in terms of both atomic root mean square differences and backbone torsion angles. For the two helical regions individually the average backbone rms difference between all pairs of structures is approximately 1 A. The long axes of the two helices lie in two planes, which are at an angle of 70-100 degrees to each other. The orientation of the helices within these planes, however, cannot be determined due to the paucity of NOEs between the two helices.     

Proximal cavity, distal histidine, and substrate hydrogen-bonding mutations modulate the activity of Amphitrite ornata dehaloperoxidase

Dehaloperoxidase (DHP) from Amphitrite ornata is the first globin that has peroxidase activity that approaches that of heme peroxidases. The substrates 2,4,6-tribromophenol (TBP) and 2,4,6-trichlorophenol are oxidatively dehalogenated by DHP to form 2,6-dibromo-1,4-benzoquinone and 2,6-dichloro-1,4-benzoquinone, respectively. There is a well-defined internal substrate-binding site above the heme, a feature not observed in other globins or peroxidases. Given that other known heme peroxidases act on the substrate at the heme edge there is great interest in understanding the possible modes of substrate binding in DHP. Stopped-flow studies (Belyea, J., Gilvey, L. B., Davis, M. F., Godek, M., Sit, T. L., Lommel, S. A., and Franzen, S. (2005) Biochemistry 44, 15637-15644) show that substrate binding must precede the addition of H2O2. This observation suggests that the mechanism of DHP relies on H2O2 activation steps unlike those of other known peroxidases. In this study, the roles of the distal histidine (H55) and proximal histidine (H89) were probed by the creation of site-specific mutations H55R, H55V, H55V/V59H, and H89G. Of these mutants, only H55R shows significant enzymatic activity. H55R is 1 order of magnitude less active than wild-type DHP and has comparable activity to sperm whale myoglobin. The role of tyrosine 38 (Y38), which hydrogen bonds to the hydroxyl group of the substrate, was probed by the mutation Y38F. Surprisingly, abolishing this hydrogen bond increases the activity of the enzyme for the substrate TBP. However, it may open a pathway for the escape of the one-electron product, the phenoxy radical leading to polymeric products.     

Crystal structure of human heme oxygenase-1.

Heme oxygenase catalyzes the first step in the oxidative degradation of heme. The crystal structure of heme oxygenase-1 (HO-1) reported here reveals a novel helical fold with the heme sandwiched between two helices. The proximal helix provides a heme iron ligand, His 25. Conserved glycines in the distal helix near the oxygen binding site allow close contact between the helix backbone and heme in addition to providing flexibility for substrate binding and product release. Regioselective oxygenation of the alpha-meso heme carbon is due primarily to steric influence of the distal helix.     

Structure of beef liver catalase

The three-dimensional structure of beef liver catalase has been determined to 2.5 å resolution by a combination of isomorphous and molecular replacement techniques. Heavy-atom positions were found using vector search and difference Fourier methods. The tetrameric catalase molecule has 222 symmetry with one of its dyads coincident with a crystallographic 2-fold axis. The known polypeptide sequence has been unambiguously fitted to the electron density map. The heme is well buried in a hydrophobic pocket, 20 Å below the surface of the molecule, and accessible through a hydrophobic channel. Residues that line the heme pocket belong to two different subunits. Tyr357 is the proximal heme ligand and the catalytically important residues on the distal side are residues His74 and Asnl47. The tertiary structure consists of four domains: an extended non-globular amino-terminal arm, which stabilizes the quaternary structure; an anti-parallel, eight-stranded β-barrel providing the residues on the distal side of the heme; a rather random “wrapping domain” around the subunit exterior including the proximal heme ligand; and a final λ-helical structure resembling the E, F, G and H helices of the globins.     

The solution structure of the oxidized bovine microsomal cytochrome b(5) mutant V61H

Using 1488 NOE constraints, 19 stereo-specific assignments, 13 pairs of H-bond constraints, and 140 pseudo-contact shift constraints, a family of 35 structures of bovine microsomal cytochrome b(5) mutant V61H has been obtained through the program PSEUDYANA. The family has been further refined by restrained energy minimization to give a family of final structures. The RMSD values of final structures with respect to the average structure are 0.45+/-0.11 and 0.96+/-0.10A for backbone and heavy atoms, respectively. The final Deltachi(ax) and Deltachi(rh) values are 2.34 x 10(-32) and -0.67 x 10(-32)m(3), respectively. The comparisons between the solution structures of mutant V61H and WT cytochrome b(5), and X-ray structure of the mutant V61H show that the global folding of the molecule in solution is unchanged and the side-chain of His61 deviates from the heme pocket and extends into the solvent like in its crystal structure. However, the helices around the heme pocket undergo outward global displacement while their local conformations are well maintained. Meanwhile, the heme ring shows a little off the heme pocket, which accounts for the lower stability of the mutant. Additionally, the axial ligand rings counterclockwise rotate around His39 N-Fe axis due to the mutation, which is confirmed by variation of the hyperfine shifts of the heme protons of V61H compared to those of WT cytochrome b(5).     

Novel subunit structure observed for noncooperative hemoglobin from Urechis caupo

Tetrameric hemoglobin from the "fat innkeeper" worm Urechis caupo possesses a novel subunit arrangement having an "inside out" quaternary structure in that the G/H helices are located on the outer surface of the tetramer. A 5-A resolution crystal structure reveals that although the individual subunits are beta-like, having a distinct D helix and the general myoglobin fold, the subunit contacts are very different from those previously observed for hemoglobins. Furthermore, the hemoglobin from U. caupo is also quite different from the unusual hemoglobin tetramer from clam which also has its G/H helices on the outer surface but with the hemes in close proximity through E-F helical contacts (Royer, W. E., Jr., Love, W. E., and Fenderson, F. F. (1985) Nature 316, 277-280).     

Truncated hemoglobins: trimming the classical ‘three-over-three’ globin fold to a minimal size

Truncated hemoglobins (trHbs) host the heme in a “two-over-two’ α-helical sandwich which results from extensive editing of the classical ‘three-over-three’ globin fold. The three-dimensional structure of trHbs is based on four main α-helices, arranged in a sort of α-helical bundle composed of two antiparallel helix pairs (B/E and G/H). Most notably, trHbs deviate from the conventional globin fold in that they display an extended loop substituting for the heme proximal F-helix observed in globins. Moreover, since efficient adaptation of a 110–130 amino acid trHb chain to host the porphyrin ring firstly requires specific chain flexibility, trHbs contain three invariant Gly-based motifs. Inspection of the trHb three-dimensional trHb structures shows that an apparent protein cavity or tunnel would connect the protein surface to an inner region very close to the heme distal site. Such a structural feature, never observed before in (non) vertebrate globins, may have substantial implications for ligand diffusion and binding properties in trHbs.     

The heme pocket of the globin domain of the globin-coupled sensor of Geobacter sulfurreducens — An EPR study

The globin-coupled sensor (GCS) of Geobacter sulfurreducens is unique amongst GCSs in that its signalling domain is a transmembrane domain with yet unknown function. In the present work we use X-band continuous-wave and pulsed electron paramagnetic resonance (EPR) to investigate the ferric form of the globin domain of the G. sulfurreducens GCS (GsGCS¹⁶²) at pH 8.5. This form shows a unique bis-histidine coordination of the heme with the F8His and E11His. In contrast with previous crystal structure data, where three conformers of the heme structure were identified, ferric GsGCS¹⁶² assumes only one conformation in frozen solution. The EPR data of ferric GsGCS162 are compared in detail with those of other bis-histidine coordinated globins, including other GCS systems.     

Solution 1H NMR study of the active site structure for the double mutant H64Q/V68F cyanide complex from mouse neuroglobin

Solution (1)H NMR spectroscopy has been carried out to investigate the molecular and electronic structures of the active site in H64Q/V68F double mutant mouse neuroglobin in the cyanomet form. Two heme orientations resulting from a 180 degrees rotation about the alpha-gamma-meso axis were observed with a population ratio about 1:1, and the clearly distinguished B isomer was used to perform the study. Based on the analysis of the dipolar shifts and paramagnetic relaxation constants, the distal Gln(64)(E7) side chain is obtained to adopt an orientation that may produce hydrogen bond between the N(epsilon)H(1) and the Fe-bound cyanide. The side chain of Phe(68)(E11) is oriented out of the heme pocket just like that in triple mutant of cyanide complex of sperm whale myoglobin. A 15 degrees rotation of the imidazole ring in axial His(96) is observed, which is close to the varphi angle determined from the crystal structure of NgbCO. The quantitative determinations of the orientation and anisotropies of the paramagnetic susceptibility tensor reveal that cyanide is tilted by 8 degrees from the heme normal which allows for contact to the Gln(64)(E7) N(epsilon)H(1). The E7 and E11 residues appear to control the direction and the extent of tilt of the bound ligand. Furthermore, the tilt of the ligand has no obvious influence on the heme heterogeneity of cyanide ligation for isomer A/B of the wild type and mutant protein, indicating that factors other than steric effects, such as polarity of heme pocket, impacts on ligand binding affinity.     

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.