Molecular Phylogeny of the Chipmunks Inferred from Mitochondrial Cytochrome b and Cytochrome Oxidase II Gene Sequences

There are currently 25 recognized species of the chipmunk genus Tamias. In this study we sequenced the complete mitochondrial cytochrome b (cyt b) gene of 23 Tamias species. We analyzed the cyt b sequence and then analyzed a combined data set of cyt b along with a previous data set of cytochrome oxidase subunit II (COII) sequence. Maximum-likelihood was used to further test the fit of models of evolution to the cyt b data. Other sciurid cyt b sequence was added to examine the evolution of Tamias in the context of other sciurids. Relationships among Tamias species are discussed, particularly the possibility of a current sorting event among taxa of the southwestern United States and the extreme divergences among the three subgenera (Neotamias, Eutamias, and Tamias).     

Novel primers for complete mitochondrial cytochrome b gene sequencing in mammals

Sequence-based species identification relies on the extent and integrity of sequence data available in online databases such as GenBank. When identifying species from a sample of unknown origin, partial DNA sequences obtained from the sample are aligned against existing sequences in databases. When the sequence from the matching species is not present in the database, high-scoring alignments with closely related sequences might produce unreliable results on species identity. For species identification in mammals, the cytochrome b (cyt b) gene has been identified to be highly informative; thus, large amounts of reference sequence data from the cyt b gene are much needed. To enhance availability of cyt b gene sequence data on a large number of mammalian species in GenBank and other such publicly accessible online databases, we identified a primer pair for complete cyt b gene sequencing in mammals. Using this primer pair, we successfully PCR amplified and sequenced the complete cyt b gene from 40 of 44 mammalian species representing 10 orders of mammals. We submitted 40 complete, correctly annotated, cyt b protein coding sequences to GenBank. To our knowledge, this is the first single primer pair to amplify the complete cyt b gene in a broad range of mammalian species. This primer pair can be used for the addition of new cyt b gene sequences and to enhance data available on species represented in GenBank. The availability of novel and complete gene sequences as high-quality reference data can improve the reliability of sequence-based species identification.     

Reconstructing mammalian phylogenies: a detailed comparison of the cytochrome B and cytochrome oxidase subunit I mitochondrial genes

The phylogeny and taxonomy of mammalian species were originally based upon shared or derived morphological characteristics. However, genetic analyses have more recently played an increasingly important role in confirming existing or establishing often radically different mammalian groupings and phylogenies. The two most commonly used genetic loci in species identification are the cytochrome oxidase I gene (COI) and the cytochrome b gene (cyt b). For the first time this study provides a detailed comparison of the effectiveness of these two loci in reconstructing the phylogeny of mammals at different levels of the taxonomic hierarchy in order to provide a basis for standardizing methodologies in the future. Interspecific and intraspecific variation is assessed and for the first time, to our knowledge, statistical confidence is applied to sequence comparisons. Comparison of the DNA sequences of 217 mammalian species reveals that cyt b more accurately reconstructs their phylogeny and known relationships between species based on other molecular and morphological analyses at Super Order, Order, Family and generic levels. Cyt b correctly assigned 95.85% of mammal species to Super Order, 94.31% to Order and 98.16% to Family compared to 78.34%, 93.36% and 96.93% respectively for COI. Cyt b also gives better resolution when separating species based on sequence data. Using a Kimura 2-parameter p-distance (x100) threshold of 1.5-2.5, cyt b gives a better resolution for separating species with a lower false positive rate and higher positive predictive value than those of COI.     

Comparison of evolutionary rates in the mitochondrial DNA cytochrome b gene and control region and their implications for phylogeny of the Cobitoidea (Teleostei: Cypriniformes)

It is widely accepted that mitochondrial DNA (mtDNA) control region evolves faster than protein encoding genes with few exceptions. In the present study, we sequenced the mitochondrial cytochrome b gene (cyt b) and control region (CR) and compared their rates in 93 specimens representing 67 species of loaches and some related taxa in the Cobitoidea (Order Cypriniformes). The results showed that sequence divergences of the CR were broadly higher than those of the cyt b (about 1.83 times). However, in considering only closely related species, CR sequence evolution was slower than that of cyt b gene (ratio of CR/cyt b is 0.78), a pattern that is found to be very common in Cypriniformes. Combined data of the cyt b and CR were used to estimate the phylogenetic relationship of the Cobitoidea by maximum parsimony, neighbor-joining, and Bayesian methods. With Cyprinus carpio and Danio rerio as outgroups, three analyses identified the same four lineages representing four subfamilies of loaches, with Botiinae on the basal-most clade. The phylogenetic relationship of the Cobitoidea was ((Catostomidae+Gyrinocheilidae)+(Botiinae+(Balitorinae+(Cobitinae+Nemacheilinae)))), which indicated that Sawada's Cobitidae (including Cobitinae and Botiinae) was not monophyletic. Our molecular phylogenetic analyses are in very close agreement with the phylogenetic results based on the morphological data proposed by Nalbant and Bianco, wherein these four subfamilies were elevated to the family level as Botiidae, Balitoridae, Cobitidae, and Nemacheilidae.     

Phylogenetic relationships among the phrynosomatid sand lizards inferred from mitochondrial DNA sequences generated by heterogeneous evolutionary processes

Nucleotide sequences of the mitochondrial protein coding cytochrome b (cyt b; 650 bp) and small-subunit 12S ribosomal RNA (approximately 350 bp) genes were used in analyses of phylogenetic relationships among extant phrynosomatid sand lizards, including an examination of competing hypotheses regarding the evolution of "earlessness." Sequences were obtained from all currently recognized species of sand lizards as well as representatives of the first and second outgroups and analyzed using both parsimony and likelihood methods. The cyt b data offer strong support for relationships that correspond with relatively recent divergences and moderate to low support for relationships reflecting more ancient divergences within the clade. These data support monophyly of the "earless" taxa, the placement of Uma as the sister taxon to the other sand lizards, and monophyly of all four taxa traditionally ranked as genera. All well-supported relationships in the 12S phylogeny are completely congruent with well-supported relationships in the cyt b phylogeny; however, the 12S data alone provide very little support for deeper divergences. Phylogenetic relationships within species are concordant with geography and suggest patterns of phylogeographic differentiation, including the conclusion that at least one currently recognized species (Holbrookia maculata) actually consists of more than one species. By independently optimizing likelihood model parameters for various subsets of the data, we found that nucleotide substitution processes vary widely between genes and among the structural and functional regions or classes of sites within each gene. Therefore, we compared competing phylogenetic hypotheses, using parameter estimates specific to those subsets, analyzing the subsets separately and in various combinations. The hypothesis supported by the cyt b data was favored over rival hypotheses in all but one of the five comparisons made with the entire data set, including the set of partitions that best explained the data, although we were unable to confidently reject (P < 0.05) alternative hypotheses. Our results highlight the importance of optimizing models and parameter estimates for different genes or parts thereof--a strategy that takes advantages of the strengths of both combining and partitioning data.     

Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

Sequence evolution and phylogenetic signal in control-region and cytochrome b sequences of rainbow fishes (Melanotaeniidae)

The nucleotide sequences of segments of the cytochrome b gene (351 bp), the tRNA(Pro) gene (49 bp), and the control region (approximately 313 bp) of mitochondrial DNA were obtained from 26 fish representing different populations and species of Melanotaenia and one species of Glossolepis, freshwater rainbow fishes confined to Australia and New Guinea. The purpose was to investigate relative rates and patterns of sequence evolution. Overall levels of divergence were similar for the cytochrome b and tRNA control-region sequences, both ranging from < 1% within subspecies to 15%-19% between genera. However, the patterns of sequence evolution differed. For the cytochrome b gene, transitions consistently exceeded transversions, the bias ranging from 4.2:1 to 2:1, depending on the level of sequence divergence. However, in the control-region sequence, a bias toward transitions (2:1) was observed only in comparisons between very similar sequences, and transversions outnumbered transitions in comparisons of divergent sequences. Graphic comparisons suggested that the control region was saturated for transitions at relatively low levels of sequence divergence but accumulated transversions at a greater rate than did the cytochrome b sequence. These distinct patterns of base substitution are associated with differences in A+T content, which is 70% for the tRNA control- region segment versus 50% for cytochrome b. A test for skewness in the distribution of lengths of random trees indicated that both segments contained phylogenetic signal. Parsimony analyses of the data from the two regions, with or without weighting schemes appropriate to the respective patterns of sequence evolution, identified the same five groupings of sequences, but the relationships among the groups differed. However, in most cases the branches uniting different combinations of groups were poorly supported, and the differences among topologies were insignificant. Considering the observed patterns of base substitution and the results of the phylogenetic analyses, we deduce that both the control region and cytochrome b are appropriate for population genetic studies but that the control region is less effective than cytochrome b for resolving relationships among divergent lineages of rainbow fishes.      

Screening and purification for novel cytochrome b5 from uncultured environmental micro-organisms

AIMS: We describe a sequence-based PCR method suitable for the isolation of a novel soluble heme-binding domain of cytochrome b(5) (cyt b(5)) gene directly from metagenomic DNA is described. METHODS AND RESULTS: Using the degenerate primer set, a cyt b(5) gene was isolated directly from metagenomic DNA. Based on the sequence-based PCR method, the similar conserved motif of cyt b(5) from Rhodopseudomonas palustris strain makes the novel target gene. The gene encoding cyt b(5) was cloned and expressed in Escherichia coli BL21 (DE3) using pET expression system. The expressed recombinant enzyme was purified by Ni-nitrilotriacetic acid affinity chromatography and characterized. CONCLUSIONS: Sequence-based strategy is an effective method for application of the novel gene from metagenomic DNA. SIGNIFICANCE AND IMPACT OF THE STUDY: Investigation of novel genes from metagenome, most of the micro-organism species are largely untapped, could represent an interesting and useful reservoir for biological processes.     

Purification and characterization of bovine adrenal cytochrome b561 expressed in insect and yeast cell systems

Bovine adrenal chromaffin granule cytochrome (cyt) b561 is a transmembrane hemoprotein that plays a key role in transporting reducing equivalents from ascorbate to dopamine-beta-hydroxylase for catecholamine synthesis. We have developed procedures for expression and purification of functional bovine adrenal cyt b561 in insect and yeast cell systems. The bovine cyt b561 coding sequence, with or without a hexahistidine-tag sequence at the C-terminus, was cloned into the pVL1392 transfer vector under the control of the polyhedrin promoter to generate recombinant baculovirus for protein expression in Sf9 insect cells (approximately 0.5 mg detergent-solubilized cyt b561/L culture). For the yeast system, the cyt b561 cDNA was modified with a hexahistidine-tag sequence at the C-terminus, and inserted into the pPICZB vector under the control of the alcohol oxidase promoter. The recombinant plasmid was transformed into Pichia pastoris GS115 competent cells to give methanol-inducible cyt b561 expression (approximately 0.7 mg detergent-solubilized cyt b561/L culture). Recombinant His-tagged cyt b561 expressed in Sf9 or Pichia cells was readily solubilized from membrane fractions with dodecyl maltoside and purified to electrophoretic homogeneity by one-step chromatography on Ni-NTA affinity resin. The purified recombinant cytochrome from both systems had a heme to protein ratio close to two and was fully functional, as judged by comparison with the spectroscopic and kinetic parameters of the endogenous cytochrome from chromaffin granules. A novel procedure for isolation of chromaffin granule membranes was developed to utilize frozen adrenal glands instead of fresh tissue.     

Diversity and phylogeny of mitochondrial cytochrome B lineages from six morphospecies of avian Haemoproteus (Haemosporida: Haemoproteidae)

Species of Haemoproteus (Haemosporida: Haemoproteidae), avian haemosporidians, have traditionally been described based on morphology of their gametocytes and on limited experimental information on their vertebrate host specificity. We investigated to what extent the morphological species are represented by monophyletic groups based on DNA sequence data using 2 different fragment lengths of the cytochrome b (cyt. b) gene. Phylogenetic reconstructions of obtained cyt. b lineages from 6 morphospecies of Haemoproteus showed that all lineages formed monophyletic clusters matching the morphospecies. Comparing our data with a recently published study showed that this is not always the case; the morphospecies H. belopolskyi consists of 2 distinct clusters of lineages that apparently have converged in morphology. However, the overall broad congruence between the molecular and morphological clustering of lineages will facilitate the integration of the knowledge obtained by traditional and molecular parasitology. Mean between morphospecies variation was 10-fold higher than the within species variation (5.5% vs. 0.54%), suggesting that Haemoproteus lineages with a genetic differentiation >5% are expected to be morphologically differentiated in most cases. When investigate the utility of 2 different fragment sizes of the cyt. b gene, the partial, 479-bp, cyt. b protocol picked up all mitochondrial (mt)DNA lineages that are found when using the full cyt. b gene, 1073 bp, suggesting that this protocol is sufficient for identification of most mtDNA lineages. All of the mtDNA lineages were associated with unique alleles when amplification was possible at a nuclear locus, strengthening the hypothesis that the designation of lineages based on mtDNA is largely genome-wide representative. We, therefore, propose the use of a cyt. b fragment of this length as a standard gene fragment for a DNA bar-coding system for avian Haemoproteus species.     

Stability and folding kinetics of structurally characterized cytochrome c-b562

The four-helix-bundle protein fold can be constructed from a wide variety of primary amino acid sequences. Proteins with this structure are excellent candidates for investigations of the relationship between folding mechanism and topology. The folding of cytochrome b(562), a four-helix-bundle heme protein, is hampered by heme dissociation. To overcome this complication, we have engineered a variant of cytochrome b(562) (cyt c-b(562)) featuring a c-type linkage between the heme and the polypeptide chain. The replacement of the native cyt b(562) leader sequence in this protein with that of a c-type cytochrome (cyt c(556)) led to high yields of fully matured and correctly folded cyt c-b(562). We have determined the X-ray crystal structure of cyt c-b(562) at 2.25 A and characterized its physical, chemical, and folding properties. These measurements reveal that the c-type linkage does not perturb the protein fold or reduction potential of the heme group. The covalent attachment of the porphyrin to the polypeptide does, however, produce a substantial change in protein stability and folding kinetics.     

Low mitochondrial DNA variation among American alligators and a novel non-coding region in crocodilians

We analyzed 1317-1823 base pairs (bp) of mitochondrial DNA sequence beginning in the 5' end of cytochrome b (cyt b) and ending in the central domain of the control region for 25 American alligators (Alligator mississippiensis) and compared these to a homologous sequence from a Chinese alligator (A. sinensis). Both species share a non-coding spacer between cyt b and tRNA(Thr). Chinese alligator cyt b differs from that of the American alligator by 17.5% at the nucleotide level and 13.8% for inferred amino acids, which is consistent with their presumed ancient divergence. Only two cyt b haplotypes were detected among the 25 American alligators (693-1199 bp surveyed), with one haplotype shared among 24 individuals. One alligator from Mississippi differed from all other alligators by a single silent substitution. The control region contained only slightly more variation among the 25 American alligators, with two variable positions (624 bp surveyed), yielding three haplotypes with 22, two, and one individuals in each of these groups. Previous genetic studies examining allozymes and the proportion of variable microsatellite DNA loci also found low levels of genetic diversity in American alligators. However, in contrast with allozymes, microsatellites, and morphology, the mtDNA data shows no evidence of differentiation among populations from the extremes of the species range. These results suggest that American alligators underwent a severe population bottleneck in the late Pleistocene, resulting in nearly homogenous mtDNA among all American alligators today.     

Formation of plastocyanin and cytochrome c-553 in different species of blue-green algae

Fifteen species from different genera of blue-green algae have been examined for their formation of plastocyanin (PC) and cytochrome c-553 (cyt c-553) in high or low Cu media. In addition to species which contain only cyt c-553 and those which completely exchange their cyt c-553 by PC, a new regulatory type was detected in which this exchange was incomplete. By comparing different species, it could be shown that either this incomplete exchange of cyt c-553 by PC as well as lack of PC in some other blue-green algae is not caused by restricted Cu uptake but is due to different biosynthetic and regulatory properties. Occurrence of PC and cyt c-553 cannot be used as a taxonomic criterium to classify blue-green algae. However, formation of either one or both of these redox components fits well into a line of evolution of the photosynthetic apparatus from the blue-green algae via green algae to higher plants.Abbreviations PC plastocyanin; cyt c-553, cytochrome c-553     

The role of the length and sequence of the linker domain of cytochrome b5 in stimulating cytochrome P450 2B4 catalysis

Cytochrome b(5) (cyt b(5)) is a 15-kDa amphipathic protein with a cytosolic amino-terminal catalytic heme domain, which is anchored to the microsomal membrane by a hydrophobic transmembrane alpha-helix at its carboxyl terminus. These two domains are connected by an approximately 15-amino acid linker domain, Ser(90)-Asp(104), which has been modified by site-directed mutagenesis to investigate whether the length or sequence of the linker influences the ability of cyt b(5) to bind ferric cytochrome P450 2B4 and donate an electron to oxyferrous (cyt P450 2B4), thereby stimulating catalysis. Because shortening the linker by 8 or more amino acids markedly inhibited the ability of cyt b(5) to bind cyt P450 2B4 and stimulate catalysis by this isozyme, it is postulated 7 amino acids are sufficient to allow a productive interaction. All mutant cyts b(5) except the protein lacking the entire 15-amino acid linker inserted normally into the microsomal membrane. Alternatively, lengthening the linker by 16 amino acids, reversing the sequence of the amino acids in the linker, and mutating conserved linker residues did not significantly alter the ability of cyt b(5) to interact with cyt P450 2B4. A model for the membrane-bound cyt b(5)-cyt P450 complex is presented.     

Comparative evolution of the mitochondrial cytochrome b gene and nuclear beta-fibrinogen intron 7 in woodpeckers

Most molecular phylogenetic studies of vertebrates have been based on DNA sequences of mitochondrial-encoded genes. MtDNA evolves rapidly and is thus particularly useful for resolving relationships among recently evolved groups. However, it has the disadvantage that all of the mitochondrial genes are inherited as a single linkage group so that only one independent gene tree can be inferred regardless of the number of genes sequenced. Introns of nuclear genes are attractive candidates for independent sources of rapidly evolving DNA: they are pervasive, most of their nucleotides appear to be unconstrained by selection, and PCR primers can be designed for sequences in adjacent exons where nucleotide sequences are conserved. We sequenced intron 7 of the beta-fibrinogen gene (beta-fibint7) for a diversity of woodpeckers and compared the phylogenetic signal and nucleotide substitution properties of this DNA sequence with that of mitochondrial-encoded cytochrome b (cyt b) from a previous study. A few indels (insertions and deletions) were found in the beta-fibint7 sequences, but alignment was not difficult, and the indels were phylogentically informative. The beta-fibint7 and cyt b gene trees were nearly identical to each other but differed in significant ways from the traditional woodpecker classification. Cyt b evolves 2.8 times as fast as beta-fibint7 (14. 0 times as fast at third codon positions). Despite its relatively slow substitution rate, the phylogenetic signal in beta-fibint7 is comparable to that in cyt b for woodpeckers, because beta-fibint7 has less base composition bias and more uniform nucleotide substitution probabilities. As a consequence, compared with cyt b, beta-fibint7 nucleotide sites are expected to enter more distinct character states over the course of evolution and have fewer multiple substitutions and lower levels of homoplasy. Moreover, in contrast to cyt b, in which nearly two thirds of nucleotide sites rarely vary among closely related taxa, virtually all beta-fibint7 nucleotide sites appear free of selective constraints, which increases informative sites per unit sequenced. However, the estimated gamma distribution used to model rate variation among sites suggests constraints on some beta-fibint7 sites. This study suggests that introns will be useful for phylogenetic studies of recently evolved groups.     

中国地鼠线粒体Cyt b基因测序及其分子进化

目的测定中国地鼠线粒体DNA细胞色素b基因部分序列,分析其分子系统进化关系。方法提取中国地鼠肝脏的总基因组DNA。设计合成特异引物进行PCR扩增,经检测进行测序。用Blast与GenBank中啮齿类其他常用实验动物的物种细胞色素b基因进行同源序列比较,分析其碱基组成及变异情况,并用邻接法、最大简约法、最小进化法构建了分子系统树,在分子水平上探讨中国地鼠和常用啮齿类实验动物的进化关系。结果获得了中国地鼠线粒体Cytb基因的部分序列,共936bp。结论中国地鼠和金黄地鼠的亲缘关系最近,与小鼠、大鼠存在的差异相对大,与豚鼠的亲缘关系最远,与传统的分类地位基本吻合。     

Molecular phylogeny of avian genus Syrmaticus based on the mitochondrial cytochrome B gene and control region

Mitochondrial DNA cytochrome b (cyt b) and control region (CR) nucleotide sequences were used to study the molecular phylogeny of the genus Syrmaticus. We found that the substitution rates among the three codon positions of cyt b were heterogeneous and the transition-transversion ratio was highly biased. As to CR sequences of the genus, most variable sites were in the peripheral domains. All molecular phylogenetic trees based on the two genes showed that: 1) the Syrmaticus was monophyletic and included five species with the following cladistic relationship: (S. reevesii, (S. soemmerringii, (S. mikado, (S. humiae and S. ellioti)))). Using the TN genetic distance of cyt b, we inferred the divergence time of the five species according to putative molecular clock and found that values were largely in agreement with the geological scenarios. The origin and speciation processes of the studied group were inferred by combining molecular and biogeographical evidences.     

A Brownian dynamics study of the interactions of the luminal domains of the cytochrome b6f complex with plastocyanin and cytochrome c6: the effects of the Rieske FeS protein on the interactions

The availability of the structures of the cytochrome b6f complex (cyt b6f), plastocyanin (PC), and cytochrome c6 (cyt c6) from Chlamydomonas reinhardtii allowed us, for the first time, to model electron transfer interactions between the luminal domains of this complex (including cyt f and the Rieske FeS protein) and its redox partners in the same species. We also generated a model structure in which the FeS center of the Rieske protein was positioned closer to the heme of cyt f than observed in the crystal structure and studied its interactions with both PC and cyt c6. Our data showed that the Rieske protein in both the original crystal structure and in our modeled structure of the cyt b6f complex did not physically interfere with binding position or orientation of PC or cyt c6 on cyt f. PC docked on cyt f with the same orientation in the presence or the absence of the Rieske protein, which matched well with the previously reported NMR structures of complexes between cyt f and PC. When the FeS center of the Rieske protein was moved close to the heme of cyt f, it even enhanced the interaction rates. Studies using a cyt f modified in the 184-191 loop showed that the cyt f structure is a more important factor in determining the rate of complex formations than is the presence or the absence of the Rieske protein or its position with respect to cyt f.     

Mapping the binding interface of the cytochrome b5-cytochrome c complex by nuclear magnetic resonance

The interaction between bovine cytochrome b(5) (cyt b(5)) and horse heart cytochrome c (cyt c) is investigated by NMR spectroscopy. Chemical shifts of cyt b(5) backbone resonances and side chain methyl resonances were monitored as a function of cyt c concentration. The shifts are small but saturatable and indicate that the binding of cyt b(5) with cyt c is in fast exchange. An equilibrium association constant of (6 +/- 3) x 10(4) M(-1) was obtained with a lower limit of 180 s(-1) for the dissociation rate of the complex. To resolve considerable ambiguities in the interpretation of the chemical shift mapping, (15)N relaxation experiments and cross-saturation experiments were used as alternative methods to map the cyt b(5)-cyt c binding interface. Results from the three experiments combined demonstrate that the conserved negatively charged region of cyt b(5) surrounding the solvent-exposed heme edge is involved in the interaction with cyt c. These data support the models proposed by Salemme and Mauk [(1976) J. Mol. Biol. 102, 563-568; (1993) Biochemistry 32, 6613-6623].