Trematode Myoglobins, Functional Molecules with a Distal Tyrosine*

(Received for publication, July 3, 1996, and in revised form, October 8, 1996)

Aftab K. Rashid Dagger §, Marie-Louise Van Hauwaert Dagger , Masoodul Haque , Ather H. Siddiqi , Ignace Lasters par , Marc De Maeyer par , Nathalie Griffon **, Michael C. Marden **, Sylvia Dewilde Dagger , Julius Clauwaert Dagger , Serge N. Vinogradov Dagger Dagger and Luc Moens Dagger §§

From the Dagger  Department of Biochemistry, University of Antwerp (UIA), B-2610 Wilrijk, Belgium, the  Department of Zoology, A. M. University, Aligarh, 202002, India, the par  Center for Transgene Technology and Gene Therapy, Vlaams Interuniversitair Instituut voor Biotechnologie, K. U. Leuven Campus Gasthuisberg O & N, B-3000 Leuven, Belgium, ** INSERM U299, Hopital de Bicêtre, 94275 Le Kremlin Bicêtre, France, and the Dagger Dagger  Department of Biochemistry, Wayne State University, School of Medicine, Detroit, Michigan 48201

ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
FOOTNOTES
Acknowledgments
REFERENCES


ABSTRACT

The myoglobins of two trematodes, Paramphistomum epiclitum and Isoparorchis hypselobagri, were isolated to homogeneity. The native molecules are monomeric with Mr 16,000-17,000 and pI 6.5-7.5. In each species, at least four different globin isoforms occur. Primary structure was determined at the protein level. The globin chains contain 147 amino acid residues. Although major determinants of the globin fold are conserved, characteristic substitutions are present. A Tyr residue occurs at the helical positions B10 and E7 (distal position). This is confirmed by NMR measurements (Zhang, W., Rashid, K. A., Haque, M., Siddiqi, A. H., Vinogradov, S. N., Moens, L. & La Mar, G. N. (1997) J. Biol. Chem. 272, 3000-3006). A distal Tyr normally provokes oxidation of the iron atom and the inability to bind oxygen, whereas a Tyr-B10 is indicative for a high oxygen affinity. In contrast, trematode myoglobins are functional molecules with a high oxygen affinity. Molecular modeling predicts two possible positions for the aromatic ring of Tyr-E7: one being outside the heme pocket making it freely accessible to the ligand and one within the heme pocket potentially able to form a second hydrogen bond with the iron-bound oxygen. A hydrogen bond between Tyr-B10 and the bound oxygen as in the Ascaris hemoglobin is predicted as well. The predicted structure may explain the high oxygen affinity of the trematode myoglobins.


INTRODUCTION

The influence of distal amino acid residues on the modulation of ligand binding characteristics to Hb has been the subject of much research (2, 3). Two key residues are His-E7, which stabilizes the iron-bound oxygen by hydrogen bonding and inhibits carbon monoxide by sterical hindrance, and Val-E11, which inhibits sterically the binding of all ligands. This hypothesis is mainly proven by the study of naturally occurring human mutants or engineered sperm whale Mb and by the comparison of the structure-function relationship in other species (4-7). For example, substitution of His-E7 to less polar residues results in reduction of the oxygen affinity due to a large increase of the dissociation rate relative to the association rate (6). Substitution of His-E7 right-arrow Tyr in the human Hb M mutants results in a strong decrease in oxygen affinity and an increased autoxidation rate leading to nonfunctional ferri-derivatives (8, 9).

In nonvertebrate globins, His-E7 is often substituted by other residues such as Gln (17.8%), Val (3.4%), Leu (3.4%), and Ile (0.68%) (10). Nevertheless, molecules such as the perienteric fluid Hb and Mb of Ascaris (Gln-E7) and the Mbs of the Aplysiidae (Val-E7) display the highest affinities (6, 11, 12). Therefore, mechanisms for the stabilization of the Fe2+-O2 complex, leading to high oxygen affinity, different from that of the vertebrate Hbs and Mbs, must exist.

Site-directed mutagenesis and x-ray analysis of domain I of Ascaris Hb show that the Fe2+-O2 complex is stabilized; hence, a high oxygen affinity is created by the formation of two hydrogen bonds: the usual one from the residue at position E7 (Gln) and an additional bond from the phenolate ion of Tyr at position B10 (11, 13, 14). This results in a decreased dissociation rate (koff) (Table I).

Table I.

Summary of kinetic constants for the reaction with O2 and CO of some hemoglobins and myoglobins (natural and engineered) with different distal (E7) amino acid residues


kon (10-6) koff lon (10-6) loff Ref.

Mb (His-E7) wild type 19 10 0.44 0.02 (58, 59)
Soybean legHb (His-E7) 116 5.55 12.7 0.01 (60)
Mb (Gln-E7) 24 130 0.98 0.01 (6)
Ascaris domI (Gln-E7) 2.8 0.01 0.3 (11)
Ascaris Mb (Gln-E7) 1.2 0.23 (55)
Mb (Val-E7) 250 23,000 7 0.05 (6)
Aplysia Mb (Val-E7) 15 70 0.5 0.02 (50)
Mb (Leu-E7) 98 4,100 26 0.02 (6)
Glycera HbII (Leu-E7) 186 1,800 26.8 0.04 (6)
Mb (Tyr-E7) 0.41 0.06 (5)
Dicrocoelium (Tyr-E7?) pH = 5 300 10 120 0.27 (30)
Dicrocoelium (Tyr-E7?) pH = 7 300 30 110 0.65 (30)

X-ray and NMR analysis of the globin of Aplysia demonstrate that the Fe2+-O2 complex is stabilized by the formation of a hydrogen bond between the oxygen and a guanidino-NH proton of Arg-E10 (7, 12, 15). Consequently, the oxygen affinity of Aplysia Mb is very near to that of vertebrate Mb (Table I).

Hbs/Mbs have been found in the phylum Platyhelminthes, in the representatives of the classes Turbellaria, Trematoda, and Monogenea but are absent from the class Cestoda (16). In Fasciola hepatica, the pigment is intracellular preferentially around the vitellaria and uterine coils. Consequently, it has to be considered as a "tissue Hb" or Mb. The function of the Mbs in these animals is completely speculative as it is assumed that they are facultative aerobes (16). A function in some way in oxygen transport and/or storage is likely. However, the Mb of F. hepatica is extremely antigenic and can be detected in the serum of the host a few days after infection (17). A similar high antigenic character of globins is observed for the Hbs/Mbs of the nematodes where the pigment is found in the excretory/secretory products of Trichostrongylus colubriformis and Nippostrongylus brasiliensis (18-20). The extracellular Hbs of Ascaris suum and Chironomus thummi thummi are also documented as potent allergens (21, 22). The purpose of the high antigenic character of the Hbs/Mbs of the parasite and their occurrence in the excretory/secretory products is unknown.

Spectral analysis of trematode Mbs and their different derivatives clearly proves the presence of a protoporphyrin IX group in a modified heme environment (23). The major Mb fraction in different trematode species displays a native Mr 15,000-17,000 (24, 25). Minor fractions with Mr 32,000-68,000 composed of monomeric globin chains can also be detected during purification (Table II).

Table II.

Summary of the physicochemical characteristics of trematode myoglobins


Native myoglobin
Globin chains
Ref.
Mr pI Mra NT

Paramphistomum epiclitum
  Monomer (major)  ---b ~17,000 7-7.3 ~15-17,000 >2 (24)
 ---c 16,906 ± 845 ~16,500 4 This paper
 ---d 16,643.2 ± 1.58 This paper
 ---e 16,643.0 (25)
  Tetramer (minor)  ---b ~65,000 This paper
Isoparorchis hypselobagri
  Monomer (major)  ---b ~17,000 4.2 ~15-17,000 >3 (25)
  Dimer (minor)  ---b ~35,000
  Monomer  ---c 16,559 ± 828 This paper
 ---e 16,499.6 This paper
Gastrothylax crumenifer
  Monomer (major)  ---b ~18,000 7.6 ~16,000 (24)
  Dimer (minor)  ---b ~30,000 (24)
  Tetramer (minor)  ---b ~71,000 (24)
Dicrocoelium dendriticum 22,000 4.51-4.53 15,500 2 (26, 27, 29)
Fasciolopsis buski 15-17,500 ~15-17,500 2 (45)
Fasciola hepatica >200,000 ? ? (17)

a  SDS-PAGE.
b  Gel filtration.
c  Equilibrium centrifugation.
d  Mass spectrometry.
e  Calculated from sequence.

The Mb of Dicrocoelium dendriticum, the small liver fluke, is the best documented trematode globin (Tables I and II). The molecule is monomeric, independent of pH and ligation state (26). Partial sequence data and NMR analysis strongly suggest that the distal His is replaced by Tyr (27, 28). Nevertheless, this monomeric Mb is fully functional, has a high oxygen affinity with a p50 of 0.095 mm Hg, due to a high kon value, a Hill coefficient of 1.0, and a pronounced acid Bohr effect (26, 29, 30). The latter was proposed to involve the formation and disruption of a salt bridge between a carboxylate side chain and a Lys or Arg residue (26). This is completely in contrast with the human Hb M variants and engineered molecules with a Tyr-E7 which all display a high autoxidation rate very rapidly producing ferri-derivatives unable to bind oxygen (8, 31). Consequently, the Dd Mb must have a specific unknown mechanism to stabilize the Fe2+-O2 complex. Unfortunately, the primary structure is only partially known, and conflicting data are available for the E-helix (27, 28).

Trematode Mbs are the most primitive animal globins known today; therefore, their Mb must resemble closely the common globin ancestor of the various kingdoms. A phylogenetic analysis based on their primary structure may reveal interesting features.

We here report the primary structure of the Mbs of the trematodes Paramphistomum epiclitum and Isoparorchis hypselobagri.


MATERIALS AND METHODS

Biological Material

P. epiclitum (Platyhelminthes; Trematoda, Paramphistomatidae) parasitic in the rumen of the common Indian water buffalo Bubalus bubalis and I. hypselobagri (Platyhelminthes; Trematoda: Isoparorchidae) parasitic in the swim bladder of a catfish Wallago attu were obtained, respectively, from a local slaughterhouse and fish market in Aligarh, India. The trematodes were washed thoroughly with normal saline and incubated for 1 h at 37 °C in 0.15 M NaCl, 8 mM glucose to make them shed their eggs and gut contents. Trematodes were stored at -80 °C until use.

Purification of Trematode Myoglobins and Globin Chains

Trematodes were homogenized in 3 volumes of 50 mM Tris-HCl, pH 7.0, at 4 °C with 1 mM phenylmethylsulfonyl fluoride. The cleared homogenate was differentially precipitated with ammonium sulfate (70% and 90% saturation), and the crude Mb fraction was collected, desalted, and passed through a Sephadex G-75 column equilibrated in 50 mM Tris-HCl, pH 7.0. Mb containing fractions were pooled and concentrated by ultrafiltration (25).

Final purification of the Mbs or globin chains was performed by semipreparative isoelectric focusing on immobilized pH gradient gel strips (pH 3-10) under native or denaturing conditions (32). Pe Mb was dissolved in 0.2% ampholytes (pH 4-10) while Ih Mb was dissolved and run in 9 M urea, 1% dithiothreitol, 2% CHAPS, 1 0.8% ampholytes (pH 3-10). Isoelectric focusing was carried out for a total of 10,875 V-h at 20 °C.

Mb bands were localized either by their brilliant red color (native) or as white bands after trichloroacetic acid precipitation (denaturing). Bands of interest were excised, washed with distilled water, and extracted repeatedly with 6 M urea, 0.1% trifluoroacetic acid. The combined extracts were dialyzed against 0.1% trifluoroacetic acid and vacuum-dried. The resulting polypeptide chain was denatured by heating in 6 M guanidine HCl, 1% 2-mercaptoethanol, 50 mM Tris-HCl, pH 7.0, at 100 °C for 3 min and subjected to RP-HPLC on a Vydac C4 column (4.1 mm) developed in a 0.1% trifluoroacetic acid/acetonitrile system. Purity of the obtained Mb and globin fractions was analyzed by one- (SDS-PAGE) and two-dimensional electrophoresis (32, 33).

Sequence Determination

Purified globins were modified with maleic anhydride, and peptides were generated by CNBr cleavage and digestion with trypsin and endoproteinase Asp-N (34). Resulting peptides were separated by RP-HPLC on a Vydac C4 column (2.1 mm) developed in a 0.1% trifluoroacetic acid/acetonitrile system. Sequencing was performed on a ABI 471-B Sequencer operated in the pulsed liquid mode as recommended by the manufacturer.

Mass Determination

Electrospray mass spectra were recorded on a VG Quattro II triple quadrupole mass spectrometer (VG, Manchester, UK) equipped with a Kontron HPLC system (Kontron Instruments, Milan, Italy) consisting of a 325 pump, a HPLC autosampler 46%, and a 332 HPLC detector. Tuning of the instrument was done by injecting a 20 pmol/ml solution of Mb. The electrospray carrier solvent used was CH3CN/H2O (50/50, v/v) containing 0.1% HCOOH. Samples were dissolved in the same solvent and injected via a Rheodyne loop injector of 100 µl. The flow rate of the carrier solvent was 40 µl/min. The capillary voltage was set at +3.81 kV.

The source temperature was 80 °C. The flow rate of the nebulizing gas and the drying gas were 200 liters/h and 350 liters/h, respectively. The cone voltage was 25 V. The mass spectrometer was scanned from mass to charge 500-1500. Spectra were recorded in the multichannel aquisition mode; results are the averaging of 6 scans. The relative molecular mass (Mr) was calculated using Masslynx software.

Equilibrium Centrifugation

Equilibrium centrifugation was performed in a Beckman Optima XL-A analytical ultracentrifuge according to Yphantis (35). Mr was calculated assuming a partial specific volume, n = 0.725 ml/g.

Absorption Spectra

Absorption spectra were measured using a SLM-Aminco (DW2000) or HP 8354 spectrophotometer. The oxidized form was obtained by addition of a 10-fold excess of ferricyanide to a 5 µM protein sample, followed by incubation at 37 °C for at least 4 h. The spectra of the oxidized form versus pH were made with phosphate buffer (pH <8) or in Tris buffer (pH >8). Kinetics of ligand rebinding were measured after photolysis with 10-ns pulses at 532 nm. Samples were 10 µM on a heme basis, equilibrated under 1 atm of CO or air for the oxygen binding studies.

The ligand replacement reaction (O2 by CO) is not appropriate in this case involving a very high oxygen affinity. The replacement by CO is very slow and involves an increasing oxygen concentration. The oxygen off-rate was therefore measured by injecting an oxyHb sample into a cuvette containing about 2 mM sodium dithionite equilibrated under nitrogen. The transition from the oxy to deoxy form was followed spectrophotometrically using the HP 8354 with diode array detection. Spectra were recorded every 5 s.

Alignment

The obtained sequences were aligned with relevant nonvertebrate globin sequences using the published templates as a guideline (10, 36, 37). Phylogenetic trees were constructed by the TREECON program (38).

Table III.

Penalty scores of some globin sequences against the vertebrate (V, template 2) and nonvertebrate (NV) globin templates (36, 37)

Abbreviations: Phys, Physeter catodon; Lupi, Lupinus luteus leghemoglobin; Asc D1, Ascaris hemoglobin domain 1; Asc Mb, Ascaris myoglobin; Pe, P. epiclitum globin 3; Ih, I. hypselobagri globin 2. 
A-motif
BC-motif
E-motif
Fg-motif
H-motif
Total
NV V NV V NV V NV V NV V NV V

Phys 0.0 0.0 0.0 0.0 0.2 0.5 0.0 0.0 0.0 0.0 0.2 0.5
Lupi 0.0 0.0 0.0 1.0 1.2 0.0 0.0 0.0 1.0 1.0 2.2 2.0
Asc D1 1.0 1.0 0.0 0.2 1.0 3.9 0.0 2.4 0.0 0.7 2.0 8.2
Asc Mb 0.0 0.0 0.0 1.2 0.0 2.9 0.7 4.1 1.0 2.4 1.7 10.6
Pe 1.0 2.0 1.0 1.7 1.2 1.0 1.4 1.7 1.5 0.0 6.1 6.4
Ih 0.0 2.0 0.0 0.7 2.0 3.0 1.4 1.4 0.7 0.0 5.1 7.1

Computer Modeling

A model of the Mb structure of Pe was built using A. suum Hb domain 1 as a template. Insertion/deletions were modeled with the loop modeling algorithm of Ref. 39 which is part of the Brugel modeling package (40). Side chains were modeled using the Dead-End Elimination algorithm (41-43). All modeling work was carried out on a SGI R4000 computer.


RESULTS

Purification and General Characteristics

The Mb of Pe and Ih were purified as described under "Materials and Methods." As an example, the purification progress for the Pe Mb is documented in Figs. 1 and 2. Similar results were obtained for Ih Mb. The physicochemical characteristics of both Mb molecules are summarized together with data available from literature in Table II.


Fig. 1. Purification of Paramphistomum epiclitum myoglobin. Purification was performed as described under "Materials and Methods." A, gel filtration chromatography on a Sephadex G-75 column of crude myoglobin fraction. B, SDS-PAGE of fractions representative for the purification. 1, markers; 2, total extract; 3, 0-70% ammonium sulfate precipitate; 4, 70-90% ammonium sulfate precipitate; 5, SG 75 monomeric fraction; 6, globin-3 after focusing (Fig. 2); 7, globin-3 after focusing and RP-HPLC (Fig. 2).
[View Larger Version of this Image (33K GIF file)]



Fig. 2. Purification of P. epiclitum globins. Purification was performed as described under "Materials and Methods." A, separation of globin isoforms by isoelectric focusing. Separation was performed on linear IPG strips with a pH gradient from 4-10. Homogeneity was verified by running on a second dimension SDS-PAGE. Four globin isoforms can be distinguished. B, RP-HPLC of globin-3 after semipreparative one-dimensional isoelectric focusing.
[View Larger Version of this Image (28K GIF file)]


Absorption Spectra

Absorption spectra of both Mbs in their oxy, deoxy, and carbonmonoxy forms are presented in Fig. 3, A and B. In contrast to the Pe Mb, the Ih Mb is oxidized to its ferric form during purification. Only the deoxy ferrous form was obtained, by addition of dithionite; subsequent addition of oxygen led to the ferric form.


Fig. 3. Absorption spectra of P. epiclitum and I. hypselobagri myoglobin. A, spectra of the oxy, deoxy, and carbon monoxide forms of Pe. B, spectra of the deoxy and met forms of Ih (note that the oxy form was not stable). C and D, spectra of the met form as a function of pH for Pe and Ih.
[View Larger Version of this Image (33K GIF file)]


The spectra of both Mbs in function of pH are presented in Fig. 3, C and D. The kinetics after photodissociation of Pe showed a high on-rate, for both oxygen and CO, as for the leghemoglobins. The CO kinetics were monoexponential to within 5% of the total amplitude and were independent of the fraction dissociation, varied by changing the laser intensity. There was thus no evidence for cooperativity within the dimeric forms. Injection of an oxy Pe sample into buffer with dithionite, equilibrated under nitrogen, revealed a slow oxygen dissociation process, requiring 30 s for 2/3 dissociation.

Primary Structure

The primary structure of Pe and Ih globin was reconstructed from N-terminal sequencing of relevant peptides generated by trypsin and protease Asp-N cleavage of the maleylated proteins (Fig. 4).


Fig. 4. Reconstruction of the primary structure of P. epiclitum and I. hypselobagri myoglobin from relevant peptides. Abbreviations: N-term, obtained by amino-terminal sequencing of the intact globin; asp-N and mal-trp, peptides obtained by digestion of the globins, respectively, with protease Asp-N and trypsin after maleylation of the substrate.
[View Larger Version of this Image (15K GIF file)]


Alignment

Using the vertebrate and nonvertebrate globin templates, the trematode globin sequences were aligned with globins having known tertiary structures as well as with partial sequence data from D. dendriticum (10, 27, 28, 37) (Fig. 5 and Table III). The presence of Tyr at the distal position (E7) was confirmed by NMR measurements (1).


Fig. 5. Alignment of P. epiclitum and I. hypselobagri myoglobin with globin sequences with known tertiary structure and with the partial sequence of D. dendriticum myoglobin (10, 27). In the heading of the alignment, the major characteristics of the globin fold are summarized: helical notation (A to H), heme contacts (p = proximal; d = distal), internal (lowercase), and surface (capitals) positions. Abbreviations and data base accession numbers: Phys, Physeter catodon Mb; GlycMCS, Glycera dibranchiata globin MII (GLB1); UrecFI, Urechis caupo FI (A25537); ScaI, Scapharca inequivalvis I (GGNKID); Aplim, A. limacina (GGGAA); CTTIII, C. thummi thummi III (GGICE3); AscaDI, A. suum Hb domain I (GLB_ASCSU); Glymaxa, Glycine maxima a (GPSYS); LupII, Lupinus luteus II (GPYL2).
[View Larger Version of this Image (27K GIF file)]


Computer Modeling

A model of the Mb structure of the trematode Pe was built by homology modeling using A. suum Hb domain I (14) as a template. This template (PDB code: 1ash) was taken from the Brookhaven data base (44). The choice of this template is motivated as follows. (i) These globins are of nonvertebrate origin. (ii) Both sequences have Tyr-B10 which in the Ascaris globin structure is known to directly hydrogen-bond to the heme-bound di-oxygen. (iii) In the vicinity of the di-oxygen, no insertions or deletions have to be modeled suggesting that the Ascaris globin structure can be taken as a reasonable template to model the heme environment of Pe Mb. The sequence alignment of Pe Mb relative to the Ascaris globin template is shown in Fig. 5. The Pe Mb model was built using the methods described above. The resulting model in the neighborhood of the heme is shown in Fig. 6A. Following visual inspection of the model it is seen that the Tyr-E7 can also adopt an alternate orientation wherein the phenol group hydrogen bonds with the iron-bound oxygen atom. Such orientation entails only a minor shift of the E-helix. Indeed, upon moving this helix by only 0.49 Å and performing an exhaustive conformational search exploring all chi  angles of Tyr-E7, it is seen that only two types of low energy solutions are found. One, where the Tyr group is turned away from the di-oxygen (as in Fig. 6A) and the other where the phenol group is indeed making the above described hydrogen bond (Fig. 6B).


Fig. 6. Stereoplot of the heme environment in the model of P. epiclitum myoglobin. A, in this view, the proximal and distal residues are situated, respectively, above and below the heme. The small spheres represent the di-oxygen ligand. B, same stereoplot as in A but wherein Tyr-E7 is moved in the heme pocket forming a hydrogen bond with the iron-bound oxygen atom.
[View Larger Version of this Image (34K GIF file)]


Phylogenetic Tree Construction

The alignment presented in Fig. 5 was extended with representative nonvertebrate globin sequences (10). From this data set, a phylogenetic tree was constructed using the neighbor-joining method (38) (Fig. 7).


Fig. 7. Phylogenetic tree of trematode myoglobins with some nonvertebrate globin sequences. P. epiclitum and I. hypselobagri myoglobin were aligned with 145 nonvertebrate globins (10), and an unrooted phylogenetic tree was constructed using the neighbor-joining method (38). Only some representative sequences are included. Abbreviations and accession numbers are as in Fig. 5 including: Tokunagayusurika akamusi V, Tylorrhynchus heterochaetus I (GGWNS), Paraspongia andersonii (A25054), Tetrahymena pyriformis (A36270), Nostoc commune (M92437), Chlamydomonas rheinhardtii, Escherichia coli (S15992), Vitreoscilla sp. (GGZLB), Candida norvegensis (S26964), Saccharomyces cerevisiae (A45323), Caenorhabditis elegans (S32218).
[View Larger Version of this Image (33K GIF file)]



DISCUSSION

Trematodes contain respiratory proteins clearly different from the Hbs/Mbs of their hosts (16). As trematodes are acoelomates, their Mbs represent the most primitive globin type among multicellular organisms. The major Mb fractions of Pe and Ih are monomeric (Mr 16,000-17,000) as are all other trematode globins characterized so far (Table II). Concentration-dependent dimeric forms are observed as well. The molecular mass obtained by different techniques clearly confirm each other. The monomeric nature of these Mbs agrees with their intracellular localization (45). Isoelectric focusing under denaturing conditions reveals different isoforms of closely related pI in the range of pH 5-7. Pe, as well as Ih, displays at least four major globin isoforms (Table II; Fig. 2). It is unknown if these isoforms are the products of different genes or if they result from post-translational modifications. Different globin isoforms, in one species, are common in other trematodes as well as in other nonvertebrates (Table II) (46). In the nematodes N. brasiliensis, Syngamus trachea, and Ostertachia ostertagi cuticular and body wall globin isoforms have been reported (19, 20). Therefore, it cannot be excluded that the trematode globin isoforms are also tissue-specific.

The absorption spectra for Pe were typical for most mammalian Hbs, with peaks for the ferrous oxy, deoxy, and CO forms at 415, 432, and 420 nm, respectively (Fig. 3A).

The oxygen affinity was so high that the deoxy and CO samples still show evidence of a small percentage (~3%) oxy form. This is most obvious for the deoxy sample near 550 nm, where the deoxy absorption peak is flattened, due to a relative minimum of the oxy form. The CO Soret band is lower than for human Hb, also indicating a perturbation by the oxy form.

Injection of an oxyMb sample into a cuvette with sodium dithionite equilibrated under nitrogen allowed a measurement of the oxygen dissociation kinetics; a time constant of 30 s was determined for the oxy to deoxy transition. This is orders of magnitude higher than the millisecond kinetics for most Mb and Hb and therefore indicates a very high oxygen affinity as for Hb Ascaris.

The stock samples were in the oxy form. The samples were incubated at 37 °C to study the transition from oxy to the oxidized form. The rate for this process was 0.04/h, a value typical for most Hbs and Mbs. This value is therefore not consistent with the correlation of a lower oxidation rate for a higher oxygen affinity (31). The same is true for Hb Ascaris, indicating that another mechanism may cause the oxidation in the case of these very high affinity proteins.

From the absorption spectra we have to conclude that the Pe Mb is fully functional with a very high oxygen affinity. This confirms previous observations on Dd Mb (Table II) but is in contrast with the human Hbs M which have a distal Tyr as well, but are unable to bind oxygen. However, Ih Mb shows the same rapid oxidation as Hbs M.

The spectra of the oxidized form were not typical (Fig. 3B). While human Hb shows an acid (Soret at 406 nm with a water molecule as ligand) to alkaline (Soret at 420 nm) transition with a pK value near pH 8, the Pe Mb showed only a partial transition. This could indicate a very low pK value (below 4), but the results were hampered by a time dependence in the spectra, especially at the lower pH values. This would indicate a difference in the residues forming the distal side of the heme pocket.

Recombination of CO and oxygen to Pe Mb showed high on-rates, as for leghemoglobins. This combined with the very low oxygen dissociation rate, as for Hb Ascaris, leads to an extremely high oxygen affinity. As for Hb Ascaris, the partition coefficient (ratio of CO to oxygen affinity) is less than 1.

In Pe Mb, the iron atom is clearly in the Fe2+ (ferrous) form, binding oxygen very strongly. This Mb is very resistant to autoxidation during the purification process, whereas the Ih Mb is in the oxy form within the trematode body but as soon as it is exposed to air, autoxidation starts. This suggests a different heme environment in both Mbs. This is confirmed by the difference of the spectra of both Mbs as a function of pH (Fig. 3B). In the hydroxy-met form, the usual 540 and 575 nm maxima are present in both Mbs although in different intensities. However, both spectra are characteristic for the unusual presence of a shoulder at 603 nm. A similar shoulder is observed in the hydroxy-metMb spectrum of Dd, and it is proposed, based on analogy with the spectra of human Hb M mutants and hemin in dimethyl sulfoxide, that it substantiates the occurrence of an axially bound tyrosinate (phenoxide) to the central iron atom. Although the sequence of the heme environment of Dd Mb is uncertain, a distal Tyr is suggested based on optical spectra and 1H NMR analysis (47-49). The similarity in the hydroxy-metMb spectra of Pe and Ih Mb with the Dd Mb spectra suggest that in the former also a E7 Tyr may occur.

The primary structure of Pe globin-3 was reconstructed from relevant peptides generated by digestion of the maleylated polypeptide with trypsin and endoprotease Asp-N as documented in Fig. 4. All necessary overlaps are presented, and the majority of residues are sequenced at least twice. Moreover, the globin primary structure (C helix to carboxyl terminus) was confirmed by sequencing the coding gene.2 Globin-3 has 147 residues and a mass of 16,643.00 calculated from amino acid sequence. This is in good agreement with the mass determined by mass spectrometry (16,643.2 ± 1.58; Table II) and argues again for the correctness of the sequence. The primary structure of Ih globin-2 is determined in a similar way. It has 147 residues and a mass of 16,499.67 calculated from amino acid sequence. Not all peptide overlaps are presented. However, as the alignment with Pe globin-3 and Dd globin is unambiguous, we are confident that no residues/peptides are missing (Fig. 4).

Using existing templates, the Pe and Ih globin sequences together with the partial sequence of Dd globin are aligned with sperm whale Mb and nonvertebrate globins with known tertiary structure (10, 27, 36, 37) (Fig. 5). The alignment is confirmed (i) by the exclusion of polar residues out of the 33 invariant nonpolar sites listed in (51, 52); (ii) by the alignment of Pro-C2, which determines the folding of the BC corner; (iii) by the heme-linked His-F8; (iv) by the Phe-CD1, which is a heme contact invariably found in globin sequences; and (v) by the alignment of Gly or Ala residues B6 and E8, which are essential for the normal close contact between helices B and E. This results in a low penalty score for the globins against the vertebrate template II (36) and the nonvertebrate template (37) suggesting that the overall globin fold is conserved.

All three trematode Mb sequences show a high degree of homology (~32%) especially in the first half of the molecule. The presence of all helical segments, including a short D-helix, can be accepted. The physicochemical nature of residues at key positions as NA2, B6, B14, C2, CD1, F8, and G5, as well as the residues in the proximity of the heme group, is conserved (Fig. 5). In all three trematode globin sequences, many side chains in the proximity of the heme pocket, especially at the distal side, are aromatic. Deviations of the classical globin pattern also occurs; Leu-A12, Tyr-B10, Tyr-C4, Tyr-E7. Tyr-F4, and Leu-H8 are rather exceptional.

At A12 and H8, the usual Trp is substituted in all three trematode sequences to a Leu, as in nematode globins, which might suggest that the interaction between A and E/H helix is more hydrophobic. The Tyr at position C4 is only precedent in the globin domains of the extracellular Hb of Artemia (53). The occurrence at this position of Ile in the Chironimid Hbs proves that there is space to harbor a bulky side chain (54).

Although an aromatic side chain at F4 is rather common, a Tyr at this position is only shared with T. heterochaetus globin I. A Pro residue is present in the middle of the G and H helix (G10, H13) so that a bend of the helices can be expected.

The three trematode globins share the unusual Tyr-B10. This also occurs in legHbs and in the globins of nematodes, protists, and yeasts. All these Hbs share a relative high oxygen affinity (10, 19, 20, 37). This is due to either a high kon (legHbs) or a low koff (Ascaris) value (6, 11). The low koff value, hence the high oxygen affinity, of the Ascaris perienteric fluid Hb is due to the formation of a strong hydrogen bond between Tyr-B10 and the distal oxygen of the ligand, combined with a weak hydrogen bond between Gln-E7 and the proximal oxygen whereas the total structure is stabilized by a third hydrogen bond between both residues involved (11, 13, 14). Preliminary results indicate a high oxygen affinity for the Pe Mb, and this confirms the measurements on the Dd Mb (30). In contrast, legHb3 shows a normal oxygen koff rate and a moderately higher oxygen affinity due to a higher kon rate.

Most remarkable, however, is the presence of a Tyr-E7 in all three trematode Mb sequences. This is already suggested by the hydroxy-met spectra described above (Fig. 3B) and confirms previous work on Dd Mb (47-49, 56). As a His is present at the E6 position, confirmation of the nature of the distal residue is essential. 1H NMR measurements on Pe Mb clearly confirm the presence of Tyr-E7, Leu-E11, and Tyr-B10 as proposed by the alignment (1). Therefore, this is the first functional monomeric Mb with confirmation of a Tyr at the distal position. As human Hbs with a distal Tyr always occur as their ferri derivatives (Hb M), unable to bind oxygen, the molecular architecture of the Pe Mb must be so that high oxygen affinity is created and the iron atom protected from oxidation.

Computer modeling of Pe Mb suggests two possible orientations for the Tyr-E7, each of which is at about the same distance (<5 Å) from the nearest atom in Phe-CD1. In one orientation, the Tyr is turned away from the heme (Fig. 6A), thereby making the heme pocket accessible to oxygen. In the other orientation (Fig. 6B), the hydroxyl group of Tyr-E7 forms a hydrogen bond with the iron-bound oxygen atom (distance between oxygen atoms equals 2.5 Å) thereby possibly stabilizing the O2 molecule in the heme pocket. The other oxygen atom hydrogen binds to Tyr-B10 (distance between oxygen atoms equals 2.7 Å) in the same way as observed in the Ascaris structure (14). Leu-E11 and Phe-G12 are in close contact with each other (2.25 Å between Leu-E11.HD21 and Phe-G12.HE2) and are buried deep into the heme pocket with the closest distance of 4.6 Å between Leu-E11.CD1 and heme-CMB. As Tyr-E7 is close to Leu-E11, this strongly suggests the localization of Tyr-E7 within the heme pocket. This hypothesis is confirmed by NMR data presented in the accompanying paper (1).

Similar reorientation of the distal residue (E7), inside or outside the heme pocket, is observed in the difference Fourier maps of the deoxy- and cyanomet form of Hb III of C. thummi (54). The same is suggested to occur, though to a lesser extent, in the leg Hbs (57). Both globin types display a high oxygen affinity due to a high kon rate created by an accessibility of the heme pocket approaching diffusion limits (6).

Another interesting feature in our model of Pe Mb is that the aromatic Tyr-F4 is not oriented parallel to His-F8. This cannot be due to the presence of an aromatic side chain (Phe) at the nearby G5 position. Indeed, in Aplysia limacina (PDB code 1mba), a Phe is observed at both G5 and F4 while Phe-F4 adopts a parallel orientation to His-F8. We believe that Tyr-F4 is turned away from His-F8 because otherwise its hydroxyl group would be devoided from hydrogen bonding interactions. In the orientation resulting from the present modeling (Fig. 6), the hydroxyl is positioned in the solvent. Interestingly, the small nearby Ala residue at E14 may favor this well solvated position of the phenol group even more. It is satisfying to observe that, independently from the modeling work, the NMR analysis detects a NOE contact between Tyr-F4 and Ala-E14 which provides additional evidence that Tyr-F4 is turned away from His-F8 (1).

The phylogenetic tree presented in Fig. 7 clearly demonstrates the primitive nature of the trematode Mbs. Both trematode globins have 43.84% identity with each other and share no more than 15-20% with all other globin sequences. Interpretation of the tree, however, must be done with caution as the oldest branches are not supported by bootstrapping.


FOOTNOTES

*   This work was supported in part by Belgian Fund for Joint Basic Research Grant 2.0023.94 (to M. L.), INSERM (France), and the Direction de la Recherche et de la Technologie (France). The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

The nucleotide sequence(s) reported in this paper has been submitted to the GenBankTM/EMBL Data Bank with accession number(s) A58345[GenBank], B58345[GenBank], P80721[GenBank], and P80722[GenBank].


§   Recipient of a postdoctoral grant from the University of Antwerp.
§§   To whom correspondence should be addressed: Dept. of Biochemistry, University of Antwerp (UIA), Universiteitsplein 1, B-2610 Antwerp, Belgium. Tel.: 32-03-820-2323; Fax: 32-03-820-2248; E-mail: lmoens{at}uia.ua.ac.be.
1    The abbreviations used are: CHAPS, 3[(3-cholamidopropyl)dimethylammonia]-1-propanesulfonate; PAGE, polyacrylamide gel electrophoresis; RP-HPLC, reversed-phase high performance liquid chromatography; Pe, Paramphistomum epiclitum; Ih, Isoparorchis hypselobagri; Dd, Dicrocoelium dendriticum.
2    S. Dewilde, K. A. Rashid, and L. Moens, manuscript in preparation.
3    Q. H. Gibson and L. Parkhurst, unpublished data.

Acknowledgments

Dr. E. Esmans and Ir. W. Van Dongen are kindly acknowledged for mass determination.


REFERENCES

  1. Zhang, W., Rashid, K. A., Haque, M., Siddiqi, A. N., Vinogradov, S. N., Moens, L., and La Mar, G. N. (1997) J. Biol. Chem. 272, 3000-3006 [Abstract/Free Full Text]
  2. Olson, J., Mathews, A. J., Rohlfs, R. J., Springer, B. A., Egeberg, K. D., Sligar, S. G., Tame, J., Renaud, J., and Nagai, K. (1988) Nature 336, 265-266 [CrossRef][Medline] [Order article via Infotrieve]
  3. Johnson, K. A., Olson, J. S., and Philips, G. N. (1989) J. Mol. Biol. 207, 459-463 [Medline] [Order article via Infotrieve]
  4. Mathews, A. J., Rohlfs, R. J., Olson, J. S., Tame, J., Renaud, J.-P., and Nagai, K. (1989) J. Biol. Chem. 264, 16573-16583 [Abstract/Free Full Text]
  5. Springer, B. A., Egeberg, K. D., Sligar, S. G., Rohlfs, R. J., Mathews, A. J., and Olson, J. (1989) J. Biol. Chem. 264, 3057-3060 [Abstract/Free Full Text]
  6. Rohlfs, R. J., Mathews, A. J., Carver, T. E., Olson, J. S., Springer, B. A., Egeberg, K. D., and Sligar, S. G. (1990) J. Biol. Chem. 265, 3168-3176 [Abstract/Free Full Text]
  7. Yamamoto, Y., Iwafune, K., Chujo, R., Inoue, Y., Imai, K., and Suzuki, T. (1992) J. Mol. Biol. 228, 343-346 [Medline] [Order article via Infotrieve]
  8. Pulsinelli, P. D., Perutz, M. F., and Nagel, R. L. (1972) Proc. Natl. Acad. Sci. U. S. A. 70, 3870-3874
  9. Dickerson, R. E., and Geis, E. (1983) Hemoglobin: Structure, Function, Evolution, and Pathology, Benjamin/Cummings Publishing Co. Inc., Menlo Park, CA
  10. Kapp, O., Moens, L., Vanfleteren, J., Trotman, C., Suzuki, T., and Vinogradov, S. (1995) Protein Sci. 4, 2179-2190 [Abstract/Free Full Text]
  11. De Baere, I., Perutz, M. F., Kiger, L., Marden, M., and Poyart, C. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 1594-1597 [Abstract]
  12. Bolognesi, M., Onesti, S., Gatti, G., Coda, A., Ascenzi, P., and Brunori, M. (1989) J. Mol. Biol. 205, 529-544 [CrossRef][Medline] [Order article via Infotrieve]
  13. De Baere, I., Liu, L., Moens, L., Van Beeumen, J., Gielens, C., Richelle, J., Trotman, C., Finch, J., Gerstein, M., and Perutz, M. (1992) Proc. Natl. Acad. Sci. U. S. A. 89, 4638-4642 [Abstract]
  14. Yang, J., Kloek, A., Goldberg, D., and Scott Mathews, F. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 4224-4228 [Abstract]
  15. Qin, J., La Mar, G., Ascoli, F., Bolognesi, M., and Brunori, M. (1992) J. Mol. Biol. 224, 891-897 [Medline] [Order article via Infotrieve]
  16. Barrett, J. (1981) Biochemistry of Parasitic Helminths, MacMillan Publishers Ltd., London
  17. McGonigle, K., and Dalton, J. P. (1995) Parasitology 111, 209-215 [Medline] [Order article via Infotrieve]
  18. Frenkel, M. J., Dopheide, T. A. A., Wagland, B. M., and Ward, C. W. (1992) Mol. Biochem. Parasitol. 50, 27-36 [Medline] [Order article via Infotrieve]
  19. Blaxter, M. L., Ingram, L., and Tweedie, S. (1994) Mol. Biochem. Parasitol. 68, 1-14 [Medline] [Order article via Infotrieve]
  20. Blaxter, M. L., Vanfleteren, J. R., Xia, J., and Moens, L. (1994) J. Biol. Chem. 269, 30181-30186 [Abstract/Free Full Text]
  21. Baur, X., Aschauer, H., Mazur, G., Dewair, M., Prelicz, H., and Steigeman, W. (1986) Science 233, 351-354 [Medline] [Order article via Infotrieve]
  22. Wittenberg, J. B., and Wittenberg, B. A. (1981) Methods Enzymol. 76, 29-42 [Medline] [Order article via Infotrieve]
  23. Haider, S. A., and Siddiqi, A. H. (1976) J. Helminthol. Soc. Wash. 50, 259-266
  24. Haque, M., Rashid, K. A., Stern, M. S., Sharma, P. K., Siddiqi, A. H., Vinogradov, S. N., and Walz, D. A. (1992) Comp. Biochem. Physiol. 101B, 673-676
  25. Rashid, K. A., Haque, M., Siddiqi, A. H., Stern, M. S., Sharma, P. K., Vinogradov, S. N., and Wals, D. A. (1993) Comp. Biochem. Physiol. 106B, 993-998
  26. Smit, J. D., Sick, H., Peterhans, A., and Gersonde, K. (1986) Eur. J. Biochem. 133, 231-237
  27. Kunz, P. A. (1975) Isolation, Characterization, and Partial Sequence Determination of a Hemeprotein from the Small Liver Fluke Dicrocoelium dendriticum. Ph.D. thesis, University of Zurich
  28. Smit, J. D., Lecomte, J., La Mar, G., and Winterhalter, K. H. (1990) in Invertebrate Dioxygen Carriers (Preaux, G., and Lontie, R., eds), pp. 61-66, Leuven University Press, Leuven, Belgium
  29. Tuchschmid, P. E., Kunz, P. A., and Wilson, K. J. (1978) Eur. J. Biochem. 88, 387-394 [Abstract]
  30. Di Iorio, E. E., Meier, U. T., Smit, J. D. G., and Winterhalter, K. H. (1985) J. Biol. Chem. 260, 2160-2164 [Abstract]
  31. Brantley, R. E., Jr., Smerdon, S. J., Wilkinson, A. J., Singleton, E. W., and Olson, J. S. (1993) J. Biol. Chem. 268, 6995-7010 [Abstract/Free Full Text]
  32. Görg, A., Boguth, G., Obermaier, C., Posch, A., and Weiss, W. (1995) Electrophoresis 16, 1079-1086 [Medline] [Order article via Infotrieve]
  33. Laemmli, U. K. (1970) Nature 227, 680-685 [Medline] [Order article via Infotrieve]
  34. Allen, G. (1989) in Laboratory Techniques in Biochemistry and Molecular Biology (Burdon, R. H., and Van Knippenberg, P. H., eds), Vol. 9, Elsevier North Holland, Amsterdam
  35. Yphantis, D. A. (1964) Biochemistry 3, 297-317
  36. Bashford, D., Chothia, C., and Lesk, A. (1987) J. Mol. Biol. 196, 199-216 [Medline] [Order article via Infotrieve]
  37. Moens, L., Vanfleteren, J., Van de Peer, Y., Peeters, K., Kapp, O., Czelusniak, J., Goodman, M., Blaxter, M., and Vinogradov, S. (1996) Mol. Biochem. Evol. 13, 324-333
  38. Van de Peer, Y., and De Wachter, R. (1994) Comput. Appl. Biosci. 10, 569-570 [Medline] [Order article via Infotrieve]
  39. Claessens, M., Van Cutsem, E., Lasters, I., and Wodak, S. (1989) Protein Eng. 2, 335-345 [Abstract]
  40. Delhaise, P., Bardiaux, M., and Wodak, S. J. (1984) J. Mol. Graphics 2, 103-106 [CrossRef]
  41. Desmet, J., De Mayer, M., Hazes, B., and Lasters, I. (1992) Nature 356, 539-542 [CrossRef]
  42. Lasters, I., and Desmet, J. (1993) Protein Eng. 6, 717-722 [Abstract]
  43. Lasters, I., De Mayer, M., and Desmet, J. (1995) Protein Eng. 8, 815-822 [Abstract]
  44. Berenstein, F. C., Koetzle, T. F., Williams, G. J. B., Meyer, E. F., Jr., Brice, M. D., Rodgers, J. R., Kennard, O., Shimanouchi, T., and Tasumi, M. (1977) J. Mol. Biol. 112, 535-542 [Medline] [Order article via Infotrieve]
  45. Cain, G. D. (1969) J. Parasitol. 55, 301-306 [Medline] [Order article via Infotrieve]
  46. Vinogradov, S., Walz, D. A., Pohajdak, B., Moens, L., Kapp, O. H., Suzuki, T., and Trotman, C. (1993) Comp. Biochem. Physiol. 106, 1-26 [CrossRef]
  47. Ascenzi, P., Brunori, M., Giacometti, G. M., Winterhalter, K. H., and Smit, J. D. (1983) J. Mol. Biol. 168, 181-191 [Medline] [Order article via Infotrieve]
  48. Lecomte, J., La Mar, G., Winterhalter, K. H., and Smit, J. D. (1984) J. Mol. Biol. 180, 357-370 [Medline] [Order article via Infotrieve]
  49. Lecomte, J., Smit, J. D., La Mar, G., and Winterhalter, K. H. (1989) J. Mol. Biol. 209, 235-247 [Medline] [Order article via Infotrieve]
  50. Wittenberg, B. A., Okazaki, T., and Wittenberg, J. B. (1965) Biochim. Biophys. Acta 111, 485-495 [Medline] [Order article via Infotrieve]
  51. Perutz, M. F., Kendrew, J. C., and Watson, H. C. (1965) J. Mol. Biol. 13, 661-678
  52. Lesk, A. M., and Chothia, C. (1980) J. Mol. Biol. 136, 225-270 [Medline] [Order article via Infotrieve]
  53. Trotman, C. N. A., Manning, A. M., Moens, L., and Tate, W. P. (1991) J. Biol. Chem. 266, 13789-13795 [Abstract/Free Full Text]
  54. Steigeman, W., and Weber, E. (1979) J. Mol. Biol. 127, 309-338 [Medline] [Order article via Infotrieve]
  55. Gibson, Q. H., and Smith, M. H. (1965) Proc. R. Soc. Lond. B Biol. Sci. 163, 206-214 [Medline] [Order article via Infotrieve]
  56. Lecomte, J. T., La Mar, G. N., Smit, J. D., Winterhalter, K. H., Smith, K. M., Langry, K., and Leung, H. K. (1987) J. Mol. Biol. 197, 101-110 [Medline] [Order article via Infotrieve]
  57. Arutyunyan, E. H., Tovlis, A. B., Crelenko, A. I., Voronova, A. A., Nekreasov, Y. V., and Vainstain, B. K. (1980) Kristallografiya 25, 526-531
  58. Kraus, D. W., and Wittenberg, J. B. (1990) J. Biol. Chem. 265, 16043-16053 [Abstract/Free Full Text]
  59. Kraus, D. W., Wittenberg, J. B., Jing-Fen, L., and Peisach, J. (1990) J. Biol. Chem. 265, 16054-16059 [Abstract/Free Full Text]
  60. Appleby, C. A., Bradbury, J. H., Morris, R. J., Wittenberg, B. A., Wittenberg, J. B., and Wright, P. E. (1983) J. Biol. Chem. 258, 2254-2259 [Abstract/Free Full Text]

©1997 by The American Society for Biochemistry and Molecular Biology, Inc.