Location of the disulfide bonds in human plasma prekallikrein: the presence of four novel apple domains in the amino-terminal portion of the molecule

The location of 16 of the 18 disulfide bonds in human plasma prekallikrein was determined by amino acid sequence analysis of cystinyl peptides produced by chemical and enzymatic digestions. A unique structure, named the apple domain, was established for each of the four tandem repeats in the amino-terminal portion of the molecule. The apple domains (90 or 91 amino acids) contain 3 highly conserved disulfide bonds linking the first and sixth, second and fifth, and third and fourth half-cystine residues present in each repeat. The fourth tandem repeat contains an extra disulfide bond that forms a second small loop within the apple domain. The carboxyl-terminal portion of plasma prekallikrein containing the catalytic region of the molecule was found to have disulfide bonds located in positions similar to those of other serine proteases.     

Necessity of conserved asparagine residues in the leucine-rich repeats of platelet glycoprotein Ib alpha for the proper conformation and function of the ligand-binding region

The polypeptides of the platelet von Willebrand factor (vWf) receptor, the GP Ib-IX-V complex, each contain tandem repeats of a sequence that assigns them to the leucine-rich repeat protein family. Here, we studied the role of conserved Asn residues in the leucine-rich repeats of GP Ib alpha, the ligand-binding subunit of the complex. We replaced the Asn residue in the sixth position of the first or sixth leucine-rich repeat (of seven) either with a bulky, charged Lys residue or with a Ser residue (sometimes found in the same position of other leucine-rich repeats) and studied the effect of the mutations on complex expression, modulator-dependent vWf binding, and interactions with immobilized vWf under fluid shear stress. As predicted, the Lys substitutions yielded more severe phenotypes, producing proteins that either were rapidly degraded within the cell (mutant N158K) or failed to bind vWf in the presence of ristocetin or roll on immobilized vWf under fluid shear stress (mutant N41K). The binding of function-blocking GP Ib alpha antibodies to the N41K mutant was either significantly reduced (AK2 and SZ2) or abolished (AN51 and CLB-MB45). Ser mutations were tolerated much better, although both mutants demonstrated subtle defects in vWf binding. These results suggest a vital role for the conserved asparagine residues in the leucine-rich repeats of GP Ib alpha for the structure and functions of this polypeptide. The finding that mutations in the first leucine-rich repeat had a much more profound effect on vWf binding indicates that the more N-terminal repeats may be directly involved in this interaction.     

Second-order repeats in Xenopus laevis finger proteins

The primary structure of 342 finger repeats encoded in 42 different cDNA clones isolated from Xenopus laevis oocyte and gastrula cDNA libraries has been determined. Comparative sequence analysis of the predicted protein sequences results in a consensus repeat sequence that has an extended conserved segment of 16 amino acid residues, including the evolutionary conserved H/C link element, connected to a highly variable segment that is located in the finger loop region. Groups of tandem finger repeats are found to be organized in distinct higher-order structural units, with a pair of mutually distinct fingers being the most frequently observed second-order repeat unit. Structural features observed are discussed in respect to existing models for Zn finger structure and function.     

Solution structure of the human oncogenic protein gankyrin containing seven ankyrin repeats and analysis of its structure--function relationship

Human gankyrin (226 residues, 24.4 kDa) is a liver oncoprotein that plays an important role in the development of human hepatocellular carcinomas. In this paper, its solution structure is reported, which is the largest ankyrin protein ever determined by NMR. The highly degenerate primary sequences of the seven ankyrin repeats presented a major challenge, which was overcome by combined use of TROSY experiments, perdeuterated samples, isotope-filtered NMR experiments, and residual dipolar couplings. The final structure was of high quality, with atomic rmsds for the backbone (N, C', and C(alpha)) and all heavy atoms (residues 4-224) of 0.69 +/- 0.09 and 1.04 +/- 0.09 A, respectively. Detailed analyses of NMR data suggested that the conserved TPLH motifs play important structural roles in stabilizing the repeating ankyrin scaffold. Gankyrin is conformationally more stable than the tumor suppressor p16(INK4A), possibly due to the structural roles of conserved residues evidenced by slowly exchanging backbone amides as well as hydrogen bonding networks involving labile side chain protons. Structural comparison with p16(INK4A) identified several residues of gankyrin that are potentially important for CDK4 binding, whereas observation of the thiol proton of C180 indicated a well-structured Rb-binding site in the helical region of the sixth ankyrin repeat. Interestingly, the CDK4-binding site and Rb-binding site located in N- and C-terminal regions, respectively, are separated by comparatively more stable ankyrin repeats and highly condensed positive surface charge. These results and analyses will shed light on the structural basis of the function of human gankyrin.     

Identification and analysis of novel amino acid sequence repeats and domains in Pyrobaculum aerophilum using computational tools

We have identified four repeats and five domains that are novel in proteins encoded by the Pyrobaculum aerophilum str. IM2 proteome using automated in silico methods. A "repeat" corresponds to a region comprising less than 55 amino acid residues that occurs more than once in the protein sequence and sometimes present in tandem. A "domain" corresponds to a conserved region comprising greater than 55 amino acid residues and may be present as single or multiple copies in the protein sequence. These correspond to (1) 85 amino acid residues AAG domain, (2) 72 amino acid residues GFGN domain, (3) 43 amino acid residues KGG repeat, (4) 25 amino acid residues RWE repeat, (5) 25 amino acid residues RID repeat, (6) 108 amino acid residues NDFA domain, (7) 140 amino acid residues VxY domain, (8) 35 amino acid residues LLPN repeat and (9) 98 amino acid residues GxY domain. A repeat or domain is characterized by specific conserved sequence motifs. We discuss the presence of these repeats and domains in proteins from other genomes and their probable secondary structure.     

Crystal structure of human 53BP1 BRCT domains bound to p53 tumour suppressor

The BRCT (BRCA1 C-terminus) is an evolutionary conserved protein-protein interacting module found as single, tandem or multiple repeats in a diverse range of proteins known to play roles in the DNA-damage response. The BRCT domains of 53BP1 bind to the tumour suppressor p53. To investigate the nature of this interaction, we have determined the crystal structure of the 53BP1 BRCT tandem repeat in complex with the DNA-binding domain of p53. The structure of the 53BP1-p53 complex shows that the BRCT tandem repeats pack together through a conserved interface that also involves the inter-domain linker. A comparison of the structure of the BRCT region of 53BP1 with the BRCA1 BRCT tandem repeat reveals that the interdomain interface and linker regions are remarkably well conserved. 53BP1 binds to p53 through contacts with the N-terminal BRCT repeat and the inter-BRCT linker. The p53 residues involved in this binding are mutated in cancer and are also important for DNA binding. We propose that BRCT domains bind to cellular target proteins through a conserved structural element termed the 'BRCT recognition motif'.     

Asparaginyl β-Hydroxylation of Proteins Containing Ankyrin Repeat Domains Influences Their Stability and Function

Recent reports have provided evidence that the β-hydroxylation of conserved asparaginyl residues in ankyrin repeat domain (ARD) proteins is a common posttranslational modification in animal cells. Here, nuclear magnetic resonance (NMR) and other biophysical techniques are used to study the effect of asparaginyl β-hydroxylation on the structure and stability of ‘consensus’ ARD proteins. The NMR analyses support previous work suggesting that a single β-hydroxylation of asparagine can stabilize the stereotypical ARD fold. A second asparaginyl β-hydroxylation causes further stabilization. In combination with mutation studies, the biophysical analyses reveal that the stabilizing effect of β-hydroxylation is, in part, mediated by a hydrogen bond between the asparaginyl β-hydroxyl group and the side chain of a conserved aspartyl residue, two residues to the N-terminal side of the target asparagine. Removal of this hydrogen bond resulted in reduced stabilization by hydroxylation. Formation of the same hydrogen bond is also shown to be a factor in inhibiting binding of hydroxylated ARDs to factor-inhibiting hypoxia-inducible factor (FIH). The effects of hydroxylation appear to be predominantly localized to the target asparagine and proximal residues, at least in the consensus ARD protein. The results reveal that thermodynamic stability is a factor in determining whether a particular ARD protein is an FIH substrate; a consensus ARD protein with three ankyrin repeats is an FIH substrate, while more stable consensus ARD proteins, with four or five ankyrin repeats, are not. However, NMR studies reveal that the consensus protein with four ankyrin repeats is still able to bind to FIH, suggesting that FIH may interact in cells with natural ankyrin repeats without resulting hydroxylation. Overall, the work provides novel biophysical insights into the mechanism by which asparaginyl β-hydroxylation stabilizes the ARD proteins and reduces their binding to FIH.     

ST proteins, a new family of plant tandem repeat proteins with a DUF2775 domain mainly found in Fabaceae and Asteraceae

Background

Many proteins with tandem repeats in their sequence have been described and classified according to the length of the repeats: I) Repeats of short oligopeptides (from 2 to 20 amino acids), including structural cell wall proteins and arabinogalactan proteins. II) Repeats that range in length from 20 to 40 residues, including proteins with a well-established three-dimensional structure often involved in mediating protein-protein interactions. (III) Longer repeats in the order of 100 amino acids that constitute structurally and functionally independent units. Here we analyse ShooT specific (ST) proteins, a family of proteins with tandem repeats of unknown function that were first found in Leguminosae, and their possible similarities to other proteins with tandem repeats.

Results

ST protein sequences were only found in dicotyledonous plants, limited to several plant families, mainly the Fabaceae and the Asteraceae. ST mRNAs accumulate mainly in the roots and under biotic interactions. Most ST proteins have one or several Domain(s) of Unknown Function 2775 (DUF2775). All deduced ST proteins have a signal peptide, indicating that these proteins enter the secretory pathway, and the mature proteins have tandem repeat oligopeptides that share a hexapeptide (E/D)FEPRP followed by 4 partially conserved amino acids, which could determine a putative N-glycosylation signal, and a fully conserved tyrosine. In a phylogenetic tree, the sequences clade according to taxonomic group. A possible involvement in symbiosis and abiotic stress as well as in plant cell elongation is suggested, although different STs could play different roles in plant development.

Conclusions

We describe a new family of proteins called ST whose presence is limited to the plant kingdom, specifically to a few families of dicotyledonous plants. They present 20 to 40 amino acid tandem repeat sequences with different characteristics (signal peptide, DUF2775 domain, conservative repeat regions) from the described group of 20 to 40 amino acid tandem repeat proteins and also from known cell wall proteins with repeat sequences. Several putative roles in plant physiology can be inferred from the characteristics found.     

Thermal stability of chicken brain α-spectrin repeat 17: a spectroscopic study

Spectrin is a rod-like multi-modular protein that is mainly composed of triple-helical repeats. These repeats show very similar 3D-structures but variable conformational and thermodynamical stabilities, which may be of great importance for the flexibility and dynamic behaviour of spectrin in the cell. For instance, repeat 17 (R17) of the chicken brain spectrin α-chain is four times less stable than neighbouring repeat 16 (R16) in terms of ?G. The structure of spectrin repeats has mainly been investigated by X-ray crystallography, but the structures of a few repeats, e.g. R16, have also been determined by NMR spectroscopy. Here, we undertook a detailed characterization of the neighbouring R17 by NMR spectroscopy. We assigned most backbone resonances and observed NOE restraints, relaxation values and coupling constants that all indicated that the fold of R17 is highly similar to that of R16, in agreement with previous X-ray analysis of a tandem repeat of the two domains. However, (15)N heteronuclear NMR spectra measured at different temperatures revealed particular features of the R17 domain that might contribute to its lower stability. Conformational exchange appeared to alter the linker connecting R17 to R16 as well as the BC-loop in close proximity. In addition, heat-induced splitting was observed for backbone resonances of a few spatially related residues including V99 of helix C, which in R16 is replaced by the larger hydrophobic tryptophan residue that is relatively conserved among other spectrin repeats. These data support the view that the substitution of tryptophan by valine at this position may contribute to the lower stability of R17.     

Identification and analysis of novel tandem repeats in the cell surface proteins of archaeal and bacterial genomes using computational tools

We have identified four novel repeats and two domains in cell surface proteins encoded by the Methanosarcina acetivorans genome and in some archaeal and bacterial genomes. The repeats correspond to a certain number of amino acid residues present in tandem in a protein sequence and each repeat is characterized by conserved sequence motifs. These correspond to: (a) a 42 amino acid (aa) residue RIVW repeat; (b) a 45 aa residue LGxL repeat; (c) a 42 aa residue LVIVD repeat; and (d) a 54 aa residue LGFP repeat. The domains correspond to a certain number of aa residues in a protein sequence that do not comprise internal repeats. These correspond to: (a) a 200 aa residue DNRLRE domain; and (b) a 70 aa residue PEGA domain. We discuss the occurrence of these repeats and domains in the different proteins and genomes analysed in this work.     

Solution structure of the SEA domain from the murine homologue of ovarian cancer antigen CA125 (MUC16)

Human CA125, encoded by the MUC16 gene, is an ovarian cancer antigen widely used for a serum assay. Its extracellular region consists of tandem repeats of SEA domains. In this study we determined the three-dimensional structure of the SEA domain from the murine MUC16 homologue using multidimensional NMR spectroscopy. The domain forms a unique alpha/beta sandwich fold composed of two alpha helices and four antiparallel beta strands and has a characteristic turn named the TY-turn between alpha1 and alpha2. The internal mobility of the main chain is low throughout the domain. The residues that form the hydrophobic core and the TY-turn are fully conserved in all SEA domain sequences, indicating that the fold is common in the family. Interestingly, no other residues are conserved throughout the family. Thus, the sequence alignment of the SEA domain family was refined on the basis of the three-dimensional structure, which allowed us to classify the SEA domains into several subfamilies. The residues on the surface differ between these subfamilies, suggesting that each subfamily has a different function. In the MUC16 SEA domains, the conserved surface residues, Asn-10, Thr-12, Arg-63, Asp-75, Asp-112, Ser-115, and Phe-117, are clustered on the beta sheet surface, which may be functionally important. The putative epitope (residues 58-77) for anti-MUC16 antibodies is located around the beta2 and beta3 strands. On the other hand the tissue tumor marker MUC1 has a SEA domain belonging to another subfamily, and its GSVVV motif for proteolytic cleavage is located in the short loop connecting beta2 and beta3.     

Hydrogen bonds in rubredoxins from mesophilic and hyperthermophilic organisms

The extent and strength of the hydrogen bond networks in rubredoxins from the hyperthermophile Pyrococcus furiosus (PfRd), and its mesophilic analogue Clostridium pasteurianum (CpRd), are examined and compared using NMR spectroscopy. NMR parameters examined in this study include through-hydrogen bond (h3)J(NC)(') scalar couplings and (1)H, (13)C, and (15)N chemical shifts, as well as covalent (1)J(NH) and (1)J(NC)(') scalar couplings. These parameters have allowed the characterization in solution of 12 hydrogen bonds in each protein. Despite a 83% sequence homology and a low RMSD for the backbone heavy atoms (0.648 A) in the crystalline state, subtle, but definite, changes have been identified in the detailed hydrogen-bonding patterns. CpRd shows an increased number of hydrogen bonds in the triple-stranded beta-sheet and an additional hydrogen bond in the multiple-turn segment including residues 14-32. On the other hand, PfRd exhibits an overall strengthening of N-H...O=C hydrogen bonds in the loops involved at the metal binding site as well as evidence for an additional NH...S(Cys) hydrogen bond involving the alanine residue 44. These data, as well as temperature dependence of the NMR parameters, suggest that the particular NMR hydrogen bond pattern found in the hyperthermophile rubredoxin leads to an increased stabilization at the metal binding pocket. It seems to result from a subtle redistribution of hydrogen-bonding interactions between the triple-stranded beta-sheet and the actual metal binding site.     

Sequential unfolding of ankyrin repeats in tumor suppressor p16

The ANK repeat is a ubiquitous 33-residue motif that adopts a beta hairpin helix-loop-helix fold. Multiple tandem repeats stack in a linear manner to produce an elongated structure that is stabilized predominantly by short-range interactions between residues close in sequence. The tumor suppressor p16(INK4) consists of four repeats and represents the minimal ANK folding unit. We found from Phi value analysis that p16 unfolded sequentially. The two N-terminal ANK repeats, which are distorted from the canonical ANK structure in all INK4 proteins and which are important for functional specificity, were mainly unstructured in the rate-limiting transition state for folding/unfolding, while the two C-terminal repeats were fully formed. A sequential unfolding mechanism could have implications for the cellular fate of wild-type and cancer-associated mutant p16 proteins.     

An undecided coiled coil: the leucine zipper of Nek2 kinase exhibits atypical conformational exchange dynamics

Leucine zippers are oligomerization domains used in a wide range of proteins. Their structure is based on a highly conserved heptad repeat sequence in which two key positions are occupied by leucines. The leucine zipper of the cell cycle-regulated Nek2 kinase is important for its dimerization and activation. However, the sequence of this leucine zipper is most unusual in that leucines occupy only one of the two hydrophobic positions. The other position, depending on the register of the heptad repeat, is occupied by either acidic or basic residues. Using NMR spectroscopy, we show that this leucine zipper exists in two conformations of almost equal population that exchange with a rate of 17 s(-1). We propose that the two conformations correspond to the two possible registers of the heptad repeat. This hypothesis is supported by a cysteine mutant that locks the protein in one of the two conformations. NMR spectra of this mutant showed the predicted 2-fold reduction of peaks in the (15)N HSQC spectrum and the complete removal of cross peaks in exchange spectra. It is possible that interconversion of these two conformations may be triggered by external signals in a manner similar to that proposed recently for the microtubule binding domain of dynein and the HAMP domain. As a result, the leucine zipper of Nek2 kinase is the first example where the frameshift of coiled-coil heptad repeats has been directly observed experimentally.     

Solution structure of a truncated anti-MUC1 DNA aptamer determined by mesoscale modeling and NMR

Mucin 1 is a well-established target for the early diagnosis of epithelial cancers. The nucleotides of the S1.3/S2.2 DNA aptamer involved in binding to variable number tandem repeat mucin 1 peptides have been identified using footprinting experiments. The majority of these binding nucleotides are located in the 25-nucleotide variable region of the total aptamer. Imino proton and 2D NMR spectra of truncated and total aptamers in supercooled water reveal common hydrogen-bonding networks and point to a similar secondary structure for this 25-mer sequence alone or embedded within the total aptamer. NMR titration experiments confirm that the TTT triloop structure is the primary binding site and show that the initial structure of the truncated aptamers is conserved upon interaction with variable number tandem repeat peptides. The thermal dependence of the NMR chemical shift data shows that the base-paired nucleotides melt cooperatively at 47 ±?4°C. The structure of the 25-mer oligonucleotide was determined using a new combined mesoscale molecular modeling, molecular dynamics and NMR spectroscopy investigation. It contains three Watson-Crick pairs, three consecutive mispairs and four Watson-Crick pairs capped by a TTT triloop motif. The 3D model structures (PDB 2L5K) and biopolymer chain elasticity molecular models are consistent with both NMR and long unconstrained molecular dynamics (10 ns) in explicit water, respectively. Database Structural data are available in the Protein Data Bank and BioMagResBank databases under the accession numbers 2L5K and 17129, respectively.     

The leucine-rich repeat domain of Internalin B folds along a polarized N-terminal pathway

The leucine-rich repeat domain of Internalin B is composed of seven tandem leucine-rich repeats, which each contain a short beta strand connected to a 3(10) helix by a short turn, and an N-terminal alpha-helical capping motif. To determine whether folding proceeds along a single, discrete pathway or multiple, parallel pathways, and to map the structure of the transition state ensemble, we examined the effects of destabilizing substitutions of conserved residues in each repeat. We find that, despite the structural redundancy among the repeats, folding proceeds through an N-terminal transition state ensemble in which the extent of structure formation is biased toward repeats one and two and includes both local and interrepeat interactions. Our results suggest that the N-terminal capping motif serves to polarize the folding pathway by acting as a fast-growing nucleus onto which consecutive repeats fold in the transition state ensemble, and highlight the importance of sequence-specific interactions in pathway selection.     

Crystallographic and molecular dynamics analysis of loop motions unmasking the peptidoglycan-binding site in stator protein MotB of flagellar motor

Background

The C-terminal domain of MotB (MotB-C) shows high sequence similarity to outer membrane protein A and related peptidoglycan (PG)-binding proteins. It is believed to anchor the power-generating MotA/MotB stator unit of the bacterial flagellar motor to the peptidoglycan layer of the cell wall. We previously reported the first crystal structure of this domain and made a puzzling observation that all conserved residues that are thought to be essential for PG recognition are buried and inaccessible in the crystal structure. In this study, we tested a hypothesis that peptidoglycan binding is preceded by, or accompanied by, some structural reorganization that exposes the key conserved residues.

Methodology/Principal Findings

We determined the structure of a new crystalline form (Form B) of Helicobacter pylori MotB-C. Comparisons with the existing Form A revealed conformational variations in the petal-like loops around the carbohydrate binding site near one end of the β-sheet. These variations are thought to reflect natural flexibility at this site required for insertion into the peptidoglycan mesh. In order to understand the nature of this flexibility we have performed molecular dynamics simulations of the MotB-C dimer. The results are consistent with the crystallographic data and provide evidence that the three loops move in a concerted fashion, exposing conserved MotB residues that have previously been implicated in binding of the peptide moiety of peptidoglycan.

Conclusion/Significance

Our structural analysis provides a new insight into the mechanism by which MotB inserts into the peptidoglycan mesh, thus anchoring the power-generating complex to the cell wall.     

Mapping of three unique Ca(2+)-binding sites in human annexin II.

Site-directed mutagenesis was employed to map and characterize Ca(2+)-binding sites in annexin II, a member of the annexin family of Ca(2+)- and phospholipid-binding proteins which serves as a major cellular substrate for the tyrosine kinase encoded by the src oncogene. Several single amino acid substitutions were introduced in the human annexin II and the various mutant proteins were scored for their affinity towards Ca2+ in different assays. The data support our previous finding [Thiel, C., Weber, K. and Gerke V. (1991) J. Biol. Chem. 266, 14,732-14,739] that a Ca(2+)-binding site is present in the third of the four repeat segments which comprise the 33-kDa protein core of annexin II. In addition to Gly206 and Thr207, which are localized in the highly conserved endonexin fold of the third repeat, Glu246 is involved in the formation of this site. Thus the architecture of this Ca(2+)-binding site in solution is very similar, if not identical, to that of Ca2+ sites identified recently in annexin V crystals [Huber, R., Schneider, M., Mayr, I., R?misch, J. and Paques, E.-P. (1990) FEBS Lett. 275, 15-21]. In addition to the site in repeat 3, we have mapped sites of presumably similar architecture in repeats 2 and 4 of annexin II. Again, an acidic amino acid which is located 40 residues C-terminal to the conserved glycine at position 4 of the endonexin fold is indispensable for high-affinity Ca2+ binding: Asp161 in the second and Asp321 in the fourth repeat. In contrast, repeat 1 does not contain an acidic amino acid at a corresponding position and also shows deviations from the other repeats in the sequence surrounding the conserved glycine. These results on annexin II together with the crystallographic information on annexin V reveal that annexins can differ in the position of the Ca2+ sites. Ca(2+)-binding sites of similar structure are present in repeats 2, 3, and 4 of annexin II while in annexin V they occur in repeats 1, 2, and 4. We also synthesized an annexin II derivative with mutations in all three Ca2+ sites. This molecule shows a greatly reduced affinity for the divalent cation. However, it is still able to bind Ca2+, indicating the presence of (an) additional Ca2+ site(s) of presumably different architecture.     

Patterns of sequence variation in the mitochondrial D-loop region of shrews

Direct sequencing of the mitochondrial displacement loop (D-loop) of shrews (genus Sorex) for the region between the tRNA(Pro) and the conserved sequence block-F revealed variable numbers of 79-bp tandem repeats. These repeats were found in all 19 individuals sequenced, representing three subspecies and one closely related species of the masked shrew group (Sorex cinereus cinereus, S. c. miscix, S. c. acadicus, and S. haydeni) and an outgroup, the pygmy shrew (S. hoyi). Each specimen also possessed an adjacent 76-bp imperfect copy of the tandem repeats. One individual was heteroplasmic for length variants consisting of five and seven copies of the 79-bp tandem repeat. The sequence of the repeats is conducive to the formation of secondary structure. A termination-associated sequence is present in each of the repeats and in a unique sequence region 5' to the tandem array as well. Mean genetic distance between the masked shrew taxa and the pygmy shrew was calculated separately for the unique sequence region, one of the tandem repeats, the imperfect repeat, and these three regions combined. The unique sequence region evolved more rapidly than the tandem repeats or the imperfect repeat. The small genetic distance between pairs of tandem repeats within an individual is consistent with a model of concerted evolution. Repeats are apparently duplicated and lost at a high rate, which tends to homogenize the tandem array. The rate of D- loop sequence divergence between the masked and pygmy shrews is estimated to be 15%-20%/Myr, the highest rate observed in D-loops of mammals. Rapid sequence evolution in shrews may be due either to their high metabolic rate and short generation time or to the presence of variable numbers of tandem repeats.