Structural and thermodynamic aspects of cooperativity in the homodimeric hemoglobin from Scapharca inaequivalvis

The homodimeric cooperative hemoglobin from the mollusk Scapharca inaequivalvis displays an unusual subunit assembly with respect to vertebrate hemoglobins. The intersubunit contact region is formed by the two heme-carrying E and F helices, which bring the two hemes in contact with each other. At variance with tetrameric vertebrate hemoglobins, the ligand binding is not accompanied by a significant quaternary transition. The major ligand-linked changes are tertiary and are limited to the heme pocket and subunit interface. These unique structural features of HbI are not easily reconciled with the classical thermodynamic models used to describe cooperative ligand binding in vertebrate hemoglobins. The lack of distinct quaternary states and the absence of allosteric effectors suggested that cooperativity in HbI is entirely homotropic in origin. Thereafter, high resolution X-ray crystallographic data displayed the preferential binding of water molecules at the intersubunit interface in the unliganded protein with respect to the liganded one. These ordered water molecules were thus proposed to act as heterotropic effectors in HbI. The contribution of specific water binding to the observed cooperativity in HbI is discussed in the framework of the enthalpy-entropy compensation effect emerging from previous accurate equilibrium oxygen binding measurements.     

Reverse Engineering the Cooperative Machinery of Human Hemoglobin

Hemoglobin transports molecular oxygen from the lungs to all human tissues for cellular respiration. Its α₂β₂ tetrameric assembly undergoes cooperative binding and releasing of oxygen for superior efficiency and responsiveness. Over past decades, hundreds of hemoglobin structures were determined under a wide range of conditions for investigation of molecular mechanism of cooperativity. Based on a joint analysis of hemoglobin structures in the Protein Data Bank (Ren, companion article), here I present a reverse engineering approach to elucidate how two subunits within each dimer reciprocate identical motions that achieves intradimer cooperativity, how ligand-induced structural signals from two subunits are integrated to drive quaternary rotation, and how the structural environment at the oxygen binding sites alter their binding affinity. This mechanical model reveals the intricate design that achieves the cooperative mechanism and has previously been masked by inconsistent structural fluctuations. A number of competing theories on hemoglobin cooperativity and broader protein allostery are reconciled and unified.     

Cyclostome Hemoglobins Are Possibly Paralogous to Gnathostome Hemoglobins

Dimeric hemoglobins in lampreys are thought to be orthologous to gnathostome hemoglobins, comprising a familiar tetrameric assembly, despite their different subunit interface. To elucidate this contradictory problem, a phylogenetic analysis of vertebrate globins was conducted. The inferred maximum-likelihood trees revealed that the cyclostome hemoglobins are closely related to STAPs, recently identified stellate cell activation-associated proteins, and are paralogous to gnathostome hemoglobins. The quaternary structural difference between cyclostome and gnathostome hemoglobins is well understandable in light of their paralogous relationship.     

Structures of Invertebrate Hemoglobins

Hemoglobin is widely distributed among the invertebrates. Intracellularhemoglobins consist of relatively small molecules with mol wtsof 15–17,000 or dimeric, tetrameric or octameric aggregatesof 15–17,000 mol wt subunits. Sequence homology is presentbut not extensive in those pigments which have been studiedand the characteristic myoglobin fold of vertebrate hemoglobinoccurs in at least two invertebrate hemoglobins. The wide arrayof aggregation states among invertebrate hemoglobins providessome simple models for understanding homotropic functional propertiesexhibited by many of these pigments. Polymeric extracellularhemoglobins are present in annelids molluscs crustacean arthropodsand nematodes. Annelid extracellular hemoglobins and chlorocruorinsconsist of 3 x 10⁶ mol wt two-tiered hexagonal arrays of submultipleswhich in turn are based on polypeptide chain subunits of molwt 14–16 000. Molluscan extracellular hemoglobins areconstructed from a different subunit arrangement. In the planorbidsnail and clam extracellular hemoglobins the subunits appearto be 175 000 and 300 000 mol wt linear series of 15–17000 dalton oxygen binding domains respectively. Planorbid snailnative hemoglobin presents circular structures 200 Å indiameter in the electron microscope with 10-fold symmetry inat least one view, and clam extracellular hemoglobins are huge345 by 1200 Å rodlike structures. Crustacean extracellularhemoglobins are also polymeric pigments and at least in a fewspecies appear to have subunits which are tandemly linked oxygenbinding domains. The polymeric hemoglobins of nematodes havemolecular weights of about 330 000. The subunit molecular weightand heme content suggest a value of 40,000 daltons which setthe nematode pigments apart from all other hemoglobins so farstudied. An overview of invertebrate hemoglobin structures andsome of the questions they pose are presented in this paper.     

Molecular Evolution of Cytochrome c Oxidase Subunit IV: Evidence for Positive Selection in Simian Primates

Cytochrome c oxidase (COX) is a multi-subunit enzyme complex that catalyzes the final step of electron transfer through the respiratory chain on the mitochondrial inner membrane. Up to 13 subunits encoded by both the mitochondrial (subunits I, II, and III) and nuclear genomes occur in eukaryotic organisms ranging from yeast to human. Previously, we observed a high number of amino acid replacements in the human COX IV subunit compared to mouse, rat, and cow orthologues. Here we examined COX IV evolution in the two groups of anthropoid primates, the catarrhines (hominoids, cercopithecoids) and platyrrhines (ceboids), as well as one prosimian primate (lorisiform), by sequencing PCR-amplified portions of functional COX4 genes from genomic DNAs. Phylogenetic analysis of the COX4 sequence data revealed that accelerated nonsynonymous substitution rates were evident in the early evolution of both catarrhines and, to a lesser extent, platyrrhines. These accelerated rates were followed later by decelerated rates, suggesting that positive selection for adaptive amino acid replacement became purifying selection, preserving replacements that had occurred. The evidence for positive selection was especially pronounced along the catarrhine lineage to hominoids in which the nonsynonymous rate was first faster than the synonymous rate, then later much slower. The rates of three types of ``neutral DNA' nucleotide substitutions (synonymous substitutions, pseudogene nucleotide substitutions, and intron nucleotide substitutions) are similar and are consistent with previous observations of a slower rate of such substitutions in the nuclear genomes of hominoids than in the nuclear genomes of other primate and mammalian lineages. Received: 22 May 1996 / Accepted: 24 November 1996     

Heart and T-Lymphocyte Cell Surfaces Both Exhibit Positive Cooperativity in Binding a Membrane-Lytic Toxin

Cobra venom cytotoxins (CTX) have been shown to disrupt cells as different as immunocytes, skeletal myocytes, erythrocytes and tumor cells. Nevertheless, even subpopulations of tumor cells are differentially susceptible to CTX by an order of magnitude. In the present study, our objective was to compare CTX-specific binding with cytolytic potency for two disparate cell types in vitro. We investigated the lytic activity of cytotoxin-III from Naja naja atra (NNA, fraction D) using heart cells and human leukemic T-cells (CEM cells). For both cell types, 50% cytolysis, assessed by tetrazolium dye conversion, occurred with μm concentrations of toxin (EC₅₀= 2.2 μm). We examined the binding of radiolabeled CTX III to both heart cells and CEM cells and found the apparent dissociation constant (K _Dapp) to be 0.69 μm and 0.75 μm, for CEM and heart cells respectively. The B _max for the CEM cells was 1.0 fmoles/cell and that for heart cells was 5.2 fmoles/cell, both exhibiting positive cooperativity between the sites (Hill coefficients 1.4, T-cells; 1.6, heart). Relatively modest dissociation constants plus high numbers of binding sites per cell are consistent with a model of CTX binding to plasma membranes by interaction with phospholipids in the bilayer. Our results suggest that the lytic activity of this cytotoxin follows its binding to a population of sites on the cells in a cooperative fashion. Received: 8 May 1995/Revised: 17 November 1995     

Phylogenetic Analysis of Components of the Eukaryotic Vesicle Transport System Reveals a Common Origin of Adaptor Protein Complexes 1, 2, and 3 and the F Subcomplex of the Coatomer COPI

Eukaryotic vesicular transport requires the recognition of membranes through specific protein complexes. The heterotetrameric adaptor protein complexes 1, 2, and 3 (AP1/2/3) are composed of two large, one small, and one medium adaptin subunit. We isolated and characterized the cDNA for Arabidopsisγ-adaptin and performed a phylogenetic analysis of all adaptin subunits (proteins) in the context of all known homologous proteins. This analysis revealed (i) that the large subunits of AP1/2/3 are homologous and (ii) two subunits of the heptameric coatomer I (COPI) complex belong to this gene family. In addition, all small subunits and the aminoterminal domain of the medium subunits of the heterotetramers are homologous to each other; this also holds for two corresponding subunits of the COPI complex. AP1/2/3 and a substructure (heterotetrameric, F-COPI subcomplex) of the heptameric COPI had a common ancestral complex (called pre-F-COPI). Since all large and all small/medium subunits share sequence similarity, the ancestor of this complex is inferred to have been a heterodimer composed of one large and one small subunit. The situation encountered today is the result of successive rounds of coordinated gene duplications of both the large and the small/medium subunits, with F-COPI being the first that separated from the ancestral pre-F-COPI. Received: 1 October 1998 / Accepted: 4 January 1999     

Structural basis for the heterotropic and homotropic interactions of invertebrate giant hemoglobin

The oxygen binding properties of extracellular giant hemoglobins (Hbs) in some annelids exhibit features significantly different from those of vertebrate tetrameric Hbs. Annelid giant Hbs show cooperative oxygen binding properties in the presence of inorganic cations, while the cooperativities of vertebrate Hbs are enhanced by small organic anions or chloride ions. To elucidate the structural basis for the cation-mediated cooperative mechanisms of these giant Hbs, we determined the crystal structures of Ca2+- and Mg2+-bound Hbs from Oligobrachia mashikoi at 1.6 and 1.7 A resolution, respectively. Both of the metal-bound structures were determined in the oxygenated state. Four Ca2+-binding sites and one Mg2+-binding site were identified in each tetramer subassembly. These cations are considered to stabilize the oxygenated form and increase affinity and cooperativity for oxygen binding, as almost all of the Ca2+ and Mg2+ cations were bound at the interface regions, forming either direct or hydrogen bond-mediated interactions with the neighboring subunits. A comparison of the structures of the oxygenated form and the partially unliganded form provides structural insight into proton-coupled cooperativity (Bohr effect) and ligand-induced transitions. Two histidine residues are assumed to be primarily associated with the Bohr effect. With regard to the ligand-induced cooperativity, a novel quaternary rotation mechanism is proposed to exist at the interface region of the dimer subassembly. Interactions among conserved residues Arg E10, His F3, Gln F7, and Val E11, together with the bending motion of the heme molecules, appear to be essential for quaternary rearrangement.     

Tryptophan substitutions at lipid-exposed positions of the gamma M3 transmembrane domain increase the macroscopic ionic current response of the Torpedo californica nicotinic acetylcholine receptor

Our previous amino-acid substitutions at the postulated lipid-exposed transmembrane segment M4 of the Torpedo californica acetylcholine receptor (AChR) focused on the alpha subunit. In this study we have extended the mutagenesis analysis using single tryptophan replacements in seven positions (I288, M291, F292, S294, L296, M299 and N300) near the center of the third transmembrane domain of the gamma subunit (γM3). All the tryptophan substitution mutants were expressed in Xenopus laevis oocytes following mRNA injections at levels close to wild type. The functional response of these mutants was evaluated using macroscopic current analysis in voltage-clamped oocytes. For all the substitutions the concentration for half-maximal activation, EC ₅₀, is similar to wild type using acetylcholine. For F292W, L296W and M299W the normalized macroscopic responses are 2- to 3-fold higher than for wild type. Previous photolabeling studies demonstrated that these three positions were in contact with membrane lipids. Each of these M3 mutations was co-injected with the previously characterized αC418W mutant to examine possible synergistic effects of single lipid-exposed mutations on two different subunits. For the γM3/αM4 double mutants, the EC ₅₀s were similar to those measured for the αC418W mutant alone. Tryptophan substitutions at positions that presumably face the interior of the protein (S294 and M291) or neighboring helices (I288) did not cause significant inhibition of channel function or surface expression of AChRs. Received: 29 January 2001/Revised: 14 May 2001     

Accumulation of Species-Specific Amino Acid Replacements That Cause Loss of Particular Protein Functions in Buchnera, an Endocellular Bacterial Symbiont

Endosymbiotic bacteria live in animal cells and are transmitted vertically at the time of the host's reproduction. In view of their small and asexual populations with infrequent chances of recombination, these endocellular bacteria are expected to accumulate mildly deleterious mutations. Previous studies showed that the DNA sequences of these bacteria evolved faster than those of free-living bacteria. In this study, we compared all the ORFs of Buchnera, an endocellular bacterial symbiont of aphids, with those of 34 other prokaryotic organisms and estimated the effect of the accelerated evolution of Buchnera on the functions of its proteins. It was revealed that Buchnera proteins contain many mutations at the sites where sequences are conserved in their orthologues in many other organisms. In addition, amino acid replacements at the conserved sites are mostly changes to physicochemically different amino acids. These results suggest that functions and conformations of Buchnera proteins have been seriously impaired or strongly modified. Indeed, extensive loss of functional motifs was observed in some Buchnera proteins. In many Buchnera proteins mutations were not detected evenly throughout each molecule but tended to accumulate in some functional units, possibly leading to loss of specific functions. As Buchnera has an unusual and limited gene repertory, it is conceivable that the manner of interactions among its proteins has been changed, and thus, functional constraints over their amino acid residues have also been changed during evolution. This may account for the loss of some functional units only in the Buchnera proteins. We obtained evidence that amino acid replacements in Buchnera were not always deleterious, but neutral or, in some cases, even positively selected. Received: 14 December 2000 / Accepted: 12 March 2001     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

Enzymic and Structural Studies on Drosophila Alcohol Dehydrogenase and Other Short-Chain Dehydrogenases/Reductases

Enzymic and structural studies on Drosophila alcohol dehydrogenases and other short-chain dehydrogenases/reductases (SDRs) are presented. Like alcohol dehydrogenases from other Drosophila species, the enzyme from D. simulans is more active on secondary than on primary alcohols, although ethanol is its only known physiological substrate. Several secondary alcohols were used to determine the kinetic parameters k_cat and K_m. The results of these experiments indicate that the substrate-binding region of the enzyme allows optimal binding of a short ethyl side-chain in a small binding pocket, and of a propyl or butyl side-chain in large binding pocket, with stereospecificity for R(−) alcohols. At a high concentration of R(−) alcohols substrate activation occurs. The k_cat and K_m values determined under these conditions are about two-fold, and two orders of magnitude, respectively, higher than those at low substrate concentrations. Sequence alignment of several SDRs of known, and unknown three-dimensional structures, indicate the presence of several conserved residues in addition to those involved in the catalyzed reactions. Structural roles of these conserved residues could be derived from observations made on superpositioned structures of several SDRs with known structures. Several residues are conserved in tetrameric SDRs, but not in dimeric ones. Two halohydrin-halide-lyases show significant homology with SDRs in the catalytic domains of these enzymes, but they do not have the structural features required for binding NAD⁺. Probably these lyases descend from an SDR, which has lost the capability to bind NAD⁺, but the enzyme reaction mechanisms may still be similar. Received: 23 May 2000 / Accepted: 4 January 2001     

Scapharca inæquivalvis Tetrameric Hemoglobin A and B Genes: Evidence for a Minigene

Vertebrate and many invertebrate globin genes have a three-exon/two-intron organization, with introns in highly conserved positions. According to the ``intron early' hypothesis, introns are the vestigial segments which flank previously independent coding sequences, thus providing evidence for the assembly of the ancient proteins by ``exon shuffling.' In this paper, we report the analysis of the genes of the bivalve mollusk Scapharca inaequivalvis tetrameric hemoglobin (HbII), which support this hypothesis, at least for the hemoglobin genes. We show the existence of ``minigenes' in the IIA and IIB globin genes, spanning part of the first and second introns, ``in frame' with the heme-binding domain coded by the second exon. Further support for the exon shuffling hypothesis can be found in the degree of identity of the ``new' translated sequences with those flanking the central protein domain of some invertebrate hemoglobins. Received: 31 July 1997 / Accepted: 12 December 1997     

Yeast holoenzyme of protein kinase CK2 requires both β and β′ regulatory subunits for its activity

Protein kinase CK2 is a highly conserved Ser/Thr protein kinase that is ubiquitous among eucaryotic organisms and appears to play an important role in many cellular functions. This enzyme in yeast has a tetrameric structure composed of two catalytic (α and/or α′) subunits and two regulatory β and β′ subunits. Previously, we have reported isolation from yeast cells four active forms of CK2, composed of αα′ββ′, α₂ββ′, α′₂ββ′ and a free α′-catalytic subunit. Now, we report that in Saccharomyces cerevisiae CK2 holoenzyme regulatory β subunit cannot substitute other β′ subunit and only both of them can form fully active enzymatic unit. We have examined the subunit composition of tetrameric complexes of yeast CK2 by transformation of yeast strains containing single deletion of the β or β′ regulatory subunits with vectors carrying lacking CKB1 or CKB2 genes. CK2 holoenzyme activity was restored only in cases when both of them were present in the cell. Additional, co-immunoprecypitation experiments show that polyadenylation factor Fip1 interacts with catalytic α subunits of CK2 and interaction with beta subunits in the holoenzyme decreases CK2 activity towards this protein substrate. These data may help to elucidate the role of yeast protein kinase CK2β/β′ subunits in the regulation of holoenzyme assembly and phosphotransferase activity.     

N-type Inactivation in the Mammalian Shaker K+ Channel Kv1.4

Mammalian voltage-gated K⁺ channels are oligomeric proteins, some of which may be composed in vivo of subunits derived from several similar genes. We have studied N-type inactivation in the rapidly inactivating Kv1.4 channel and, in specific, heteromultimers of this gene product with Kv1.5 noninactivating subunits. Heteromultimeric channels were analyzed for the stoichiometry of Kv1.4:Kv1.5 subunits by observing shifts in the midpoints of steady-state availability from that of homomultimeric channels. This analysis was employed to examine inactivation of heteromultimeric channels expressed in Xenopus oocytes using two model systems: by expression of a Kv1.4–Kv1.5 tandem fusion construct and by coexpression of native Kv1.4 and Kv1.5 channels across a wide relative concentration range of microinjected mRNA. Additionally, inactivation was examined in coexpression experiments of N-terminal deletion mutants of Kv1.4. We found that (i) a single inactivating subunit conferred inactivation in all hetero-multimers studied; (ii) the rate of inactivation could not be distinguished in channels containing two inactivating subunits from those containing one inactivating subunit; and (iii) large deletions in the linker region between the N-terminal inactivation region and the first membrane-spanning domain had no effect on the rate of inactivation. These data confirm the importance of the proximal N-terminal region in the inactivation of mammalian Kv1.4 channels, and suggest that the inactivation particle remains in close proximity to the permeation pathway even when the channel is in the open state. Received: 24 August 1995/Revised: 7 February 1996     

Conformational dynamics of the tetrameric hemoglobin molecule as revealed by hydrogen exchange: 2. Effect of intersubunit contact cleavage

IR spectroscopy was used to study the rate of hydrogen-deuterium (H-D) exchange of peptide NH atoms in isolated α and β subunits of human hemoglobin (Hb) at pH 5.5–9.0 and 20°C. The H-D exchange occurs by the EX2 mechanism. The retardation factor of subunit exchange rate (P) is in a range of approximately 10²–10⁷. Compared to tetrameric Hb, the probability of local fluctuations (1/P) increases to a slightly greater extent in monomeric α subunits than in tetrameric β subunits. Unlike in the whole Hb molecule, oxygenation of its subunits has no effect on the probability of local fluctuations, and the subunits show no pH-dependent changes in 1/P values (observed for liganded Hb). Probable mechanisms accounting for overall intensification of local fluctuations upon the cleavage of contacts between subunits of the tetrameric Hb molecule are discussed with regard to structural crystallographic data.     

Erratic Evolution of Glycerol-3-Phosphate Dehydrogenase in Drosophila,Chymomyza, and Ceratitis

We have studied the evolution of Gpdh in 18 fruitfly species by sequencing 1,077 nucleotides per species on average. The region sequenced includes four exons coding for 277 amino acids and three variable-length introns. Phylogenies derived by a variety of methods confirm that the nominal genus Zaprionus belongs within the genus Drosophila, whereas Scaptodrosophila and Chymomyza are outside. The rate of GPDH evolution is erratic. The rate of amino acid replacements in a lineage appears to be 1.0 × 10⁻¹⁰/site/year when Drosophila species are considered (diverged up to 55 million years ago), but becomes 2.3 × 10⁻¹⁰ when they are compared to Chymomyza species (divergence around 60 My ago), and 4.6 × 10⁻¹⁰ when species of those two genera are compared with the medfly Ceratitis capitata (divergence around 100 My ago). In order to account for these observations, the rate of amino acid replacement must have been 15 or more times greater in some lineages and at some times than in others. At the nucleotide level, however, Gpdh evolves in a fairly clockwise fashion. Received: 13 June 1996 / Accepted: 16 August 1996     

Identification of Major Phylogenetic Branches of Inhibitory Ligand-Gated Channel Receptors

The gene superfamily of ligand-gated ion channel (LGIC) receptors is composed of members of excitatory LGIC receptors (ELGIC) and inhibitory LGIC receptors (ILGIC), all using amino acids as ligands. The ILGICs, including GABA_A, Gly, and GluCl receptors, conduct Cl⁻ when the ligand is bound. To evaluate the phylogenetic relationships among ILGIC members, 90 protein sequences were analyzed by both maximum-parsimony and distance matrix-based methods. The strength of the resulting phylogenetic trees was evaluated by means of bootstrap. Four major phylogenetic branches are recognized. Branch I, called BZ, for the majority of the members are known to be related to benzodiazepine binding, is subdivided into IA, composed of all GABA_A receptor α subunits, and IB, composed of the γ and ε subunits, which are shown to be tightly linked. Branch II, named NB for non–benzodiazepine binding, and consisting of GABA_A receptor β, δ, π, and ρ subunits, is further subdivided into IIA, containing β subunits; IIB, containing δ, and π subunits; and IIC, containing ρ subunits. Branch IIIA, composed of vertebrate Gly receptors, is loosely clustered with Branch IIIB, composed of invertebrate GluCl receptors, to form Branch III, which is designated NA for being non–GABA responsive. Branch IV is called UD for being undefined in specificity. The existence of primitive forms of GABA_A receptor non-β subunits in invertebrates is first suggested by the present analysis, and the identities of sequences p25123 from Drosophila melanogaster, s34469 from Lymnaea stagnalis, and u14635 and p41849 from C. aenorhabditis elegans are determined to be different from their previously given annotations. The proposed branching classification of ILGICs provides a phylogenetic map, based on protein sequences, for tracing the evolutionary pathways of ILGIC receptor subunits and determining the identities of newly discovered subunits on the basis of their protein sequences. Received: 15 April 1997 / Accepted: 11 March 1998     

The Prokaryote-to-Eukaryote Transition Reflected in the Evolution of the V/F/A-ATPase Catalytic and Proteolipid Subunits

Changes in the primary and quarternary structure of vacuolar and archaeal type ATPases that accompany the prokaryote-to-eukaryote transition are analyzed. The gene encoding the vacuolar-type proteolipid of the V-ATPase from Giardia lamblia is reported. Giardia has a typical vacuolar ATPase as observed from the common motifs shared between its proteolipid subunit and other eukaryotic vacuolar ATPases, suggesting that the former enzyme works as a hydrolase in this primitive eukaryote. The phylogenetic analyses of the V-ATPase catalytic subunit and the front and back halves of the proteolipid subunit placed Giardia as the deepest branch within the eukaryotes. Our phylogenetic analysis indicated that at least two independent duplication and fusion events gave rise to the larger proteolipid type found in eukaryotes and in Methanococcus. The spatial distribution of the conserved residues among the vacuolar-type proteolipids suggest a zipper-type interaction among the transmembrane helices and surrounding subunits of the V-ATPase complex. Important residues involved in the function of the F-ATP synthase proteolipid have been replaced during evolution in the V-proteolipid, but in some cases retained in the archaeal A-ATPase. Their possible implication in the evolution of V/F/A-ATPases is discussed. Received: 27 August 1997 / Accepted: 14 January 1998     

Selection and Methionine Accumulation in the Fat Body Protein 2 Gene (FBP2), a Duplicate of the Drosophila Alcohol Dehydrogenase (ADH) Gene

The Drosophila fat body protein 2 gene (Fbp2) is an ancient duplication of the alcohol dehydrogenase gene (Adh) which encodes a protein that differs substantially from ADH in its methionine content. In D. melanogaster, there is one methionine in ADH, while there are 51 (20% of all amino acids) in FBP2. Methionine is involved in 46% of amino acid replacements when Fbp2 DNA sequences are compared between D. melanogaster and D. pseudoobscura. Methionine accumulation does not affect conserved residues of the ADH-ADH^r-FBP2 multigene family. The multigene family has evolved by replacement of mildly hydrophobic amino acids by methionine with no apparent reversion. Its short-term evolution was compared between two Drosophila species, while its long-term evolution was compared between two genera belonging respectively to acalyptrate and calyptrate Diptera, Drosophila and Sarcophaga. The pattern of nucleotide substitution was consistent with an independent accumulation of methionines at the Fbp2 locus in each lineage. Under a steady-state model, the rate of methionine accumulation was constant in the lineage leading to Drosophila, and was twice as fast as that in the calyptrate lineage. Substitution rates were consistent with a slight positive selective advantage for each methionine change in about one-half of amino acid sites in Drosophila. This shows that selection can potentially account for a large proportion of amino acid replacements in the molecular evolution of proteins. Received: 12 December 1994 / Accepted: 15 April 1996