Sequence variability in the retinal-attachment domain of mammalian rhodopsins

Phosphorylation-induced conformational changes of photoactivated rhodopsin probed by fluorescent labeling at Cys140 and Cys316

In order to monitor conformational changes following photoactivation and phosphorylation of bovine rhodopsin, the two reactive sulfhydryl groups at Cys¹⁴⁰ and Cys³¹⁶ were specifically labeled with the monobromobimane (mBBr) fluorophore. Although alterations in conformation after light exposure of rhodopsin were not detected by fluorescence excitation scans (300-450 nm) of the mBBr-labeled protein, the fluorescence signal was reduced ∼ 90% in samples containing photoactivated phosphorhodopsin. Predominant labeling at either Cys¹⁴⁰ or Cys³¹⁶ in light-activated and phosphorylated rhodopsin merely generated a decrease of ∼38% and 28%, respectively, in the fluorescence excitation intensity. Thus, neither mBBr-modified Cys¹⁴⁰ nor mBBr-modified Cys³¹⁶ were involved single-handedly in the remarkable fall seen on the signal following phosphorylation of the protein; rather, the incorporation of phosphate groups on the mBBr-labeled light-activated rhodopsin appeared to affect its fluorescence signal in a cooperative or synergistic manner. These findings demonstrated that the phosphorylation of specific hydroxyl groups at the carboxyl terminal tail of rhodopsin causes definite conformational changes in the three-dimensional fold of the protein. Apparently, amino acid residues that are buried in the interior of the inactive protein become accessible following illumination and phosphorylation of rhodopsin, quenching in turn the fluorescence excitation signal of mBBr-modified rhodopsin.

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models. The process for generating the inactive and active models was validated by the hit rates, yields, and enrichment factors determined for the selection of antagonists in the inactive model and for the selection of agonists in the active model from a set of nonadrenergic GPCR drug-like ligands in a virtual screen using ligand docking software. The simulation results provide new insights into the relationships observed between selected biochemical data, the crystal structure of rhodopsin, and the structural rearrangements that occur during activation.

Anti-bovine rhodopsin monoclonal antibody recognizing light-dependent structural change

The antigenic structure of the bovine rhodopsin molecule was investigated by using a bovine rhodopsin-specific monoclonal antibody designated Rh 29. Competition assay with sealed intact disks and broken disks indicated that the antibody-binding region was localized in the intradiscal surface. An antigenic peptide obtained by a cyanogene bromide cleavage of rhodopsin was purified and determined as residues 2-39 in the amino acid sequence. Further analysis suggested that the antigenic determinant included at least residues 21-25. These results were consistent with the structural model for membrane topology of rhodopsin. The antigenicity of the rhodopsin was compared among several states. The antibody bound to both ammonyx LO-solubilized unbleached and bleached rhodopsin. In contrast, upon membrane-embedded rhodopsin, unbleached one was 100-times less antigenic than bleached one. The results suggested that the segment around the determinant of membrane-embedded rhodopsin should undergo a structural change upon absorption of light. Rh 29 detected a band corresponding to bovine, porcine and octopus opsins in immunoblotting. Protein blot of crayfish rhabdome did not show any reactive band. These bands except for crayfish reacted with concanavalin A as well. The N-terminal structure may, therefore, conserved between mammal and erthropoda and diverge between them and cepharopoda.

G protein-coupled receptor activation: analysis of a highly constrained, "straitjacketed" rhodopsin

G protein-coupled receptor (GPCR) activation is generally assumed to result in a significant structural rearrangement of the receptor, presumably involving the rigid body movement of transmembrane helices. We have investigated the activation of the GPCR rhodopsin by the construction and analysis of a mutant which contains a total of four disulfide bonds connecting the cytoplasmic ends of helices 1 and 7, and 3 and 5, and the extracellular ends of helices 3 and 4, and 5 and 6. Despite the constraints imposed by four disulfides, this "straitjacketed" receptor retains the ability to activate the G protein transducin and, therefore, provides insight into the molecular mechanism of the initial step in signal transduction of this important class of receptors.

Disulfide structure of the pheromone binding protein from the silkworm moth, Bombyx mori

Disulfide bond formation is the only known posttranslational modification of insect pheromone binding proteins (PBPs). In the PBPs from moths (Lepidoptera), six cysteine residues are highly conserved at positions 19, 50, 54, 97, 108 and 117, but to date nothing is known about their respective linkage or redox status. We used a multiple approach of enzymatic digestion, chemical cleavage, partial reduction with Tris-(2-carboxyethyl)phosphine, followed by digestion with endoproteinase Lys-C to determine the disulfide connectivity in the PBP from Bombyx mori (BmPBP). Identification of the reaction products by on-line liquid chromatography-electrospray ionization mass spectrometry (LC/ESI-MS) and protein sequencing supported the assignment of disulfide bridges at Cys-19-Cys-54, Cys-50-Cys-108 and Cys-97-Cys-117. The disulfide linkages were identical in the protein obtained by periplasmic expression in Escherichia coli and in the native BmPBP.

Tertiary interactions between the fifth and sixth transmembrane segments of rhodopsin

We have used cysteine scanning mutagenesis and disulfide cross-linking in a split rhodopsin construct to investigate the secondary structure and tertiary contacts of the fifth (TM5) and sixth (TM6) transmembrane segments of rhodopsin. Using a simple increase in pH to promote disulfide bond formation, three cross-links between residues on the extracellular side of TM5 (at positions 198, 200, and 204) and TM6 (at position 276) have been identified and characterized. The helical pattern of cross-linking observed indicates that the fifth transmembrane helix extends through residue 200 but does not include residue 198. Rhodopsin mutants containing these disulfides demonstrate nativelike absorption spectra and light-dependent activation of transducin, suggesting that large movements on the extracellular side of TM5 with respect to TM6 are not required for receptor activation.

Carboxyl-terminal processing of the cytoplasmic NAD-reducing hydrogenase of Alcaligenes eutrophus requires the hoxW gene product

Two open reading frames (ORFs) were identified immediately downstream of the four structural genes for the soluble hydrogenase (SH) of Alcaligenes eutrophus H16. While a mutation in ORF2 had no obvious effect on hydrogen metabolism, an in-frame deletion in ORF1, subsequently designated hoxW, led to a complete loss of SH activity and hence a significant retardation of autotrophic growth on hydrogen. Hydrogen oxidation in the hoxW mutant was catalyzed by the second hydrogenase, a membrane-bound enzyme. Assembly of the four subunits of the SH was blocked in mutant cells, and HoxH, the hydrogen-activating subunit, accumulated as a precursor which was still capable of binding nickel. Protein sequencing revealed that HoxH isolated from the wild type terminates at His-464, whereas the C-terminal amino acid sequence of HoxH from the hoxW mutant is colinear with the deduced sequence. Processing of the HoxH precursor was restored in vitro by a cell extract containing HoxW. These results indicate that HoxW is a highly specific carboxyl-terminal protease which releases a 24-amino-acid peptide from HoxH prior to progression of subunit assembly.

Nucleotide sequence of the canine rod-opsin-encoding gene

The major parts of two canine rod-specific opsin (Ops) transcripts have been cloned by polymerase chain reaction from retinal mRNA. Both transcripts are derived from the same gene. The 5' leader sequence of the transcripts was cloned from canine peripheral blood DNA. The transcripts code for a protein of 348 amino acids (aa), M(r) 38,962 (prior to any protein modification). The aa sequence suggests that in common with other sequenced Ops, canine rod Ops contains seven transmembrane domains, and residues believed essential for retinal pigment binding and for palmitate binding are conserved in the canine protein. Northern blotting using the central part of the ops gene as probe suggested that mature transcripts of three different sizes (about 1900, 2600 and 5500 bases) were found in retina. Of these, the 2600-base transcript was the most abundant. RACE cloning of the 3' end of ops showed that at least two of these size classes originate from differential transcript termination.

An analysis of the periodicity of conserved residues in sequence alignments of G-protein coupled receptors. Implications for the three-dimensional structure

Twenty-three sequences from the family of G-protein coupled receptors have been aligned according to the 'historical alignment' procedure of Feng and Doolittle. Fourier transform analysis of this reveals that parts of five of the seven putative membrane-spanning regions exhibit a periodicity of conserved/nonconserved residues which is compatible with the periodicity of the alpha-helix. This would place the conserved residues on one side of the helix, which may face the inside of the proposed seven membered helical bundle.

Amino acid sequence of the retinal binding site of squid visual pigment

The retinylpeptides of visual pigments of two species of squid were identified in invertebrate visual pigments. Their primary structures were identical: H-Phe-Ala-Lys-Ala-Ser-Ala-Ile-His-Asn-pro-Hse(Met)-OH. The sequence was homologous to those of the corresponding region of other visual pigments, but the eighth amino acid, His, was found in squid visual pigments. In this experiment the retinylpeptides of eleven amino acid residues were isolated by monitoring the absorbance spectrum of the reduced retinal Schiff base without using radio-active [3H]retinal. This method is valid for the isolation and identification of retinylpeptides of other invertebrate visual pigments in which the chromophore is not exchangeable.

Ectopic expression of a minor Drosophila opsin in the major photoreceptor cell class: distinguishing the role of primary receptor and cellular context

We have used P-element-mediated transformation to introduce the cloned Rh1 rhodopsin gene into the germ line of Drosophila and fully rescue the visual phenotype of mutant ninaE flies. A transcriptional fusion between the ninaE promoter and the structural gene for a minor opsin (Rh2) that is not normally expressed in the R1-R6 photoreceptor cells was used to demonstrate that Rh2 rhodopsin can photoactivate the R1-R6 transduction cascade, but with different spectral sensitivity. In addition, we show that two mutants that specifically affect the R1-R6 cells, ninaA and rdgB, do not directly affect expression of the ninaE gene.

Purification and characterization of a 22-kDa microsomal protein from rat parotid gland which is phosphorylated following stimulation by agonists involving cAMP as second messenger

Stimulation of secretion in exocrine glands by agonists involving cAMP as second messenger leads to the phosphorylation of the ribosomal protein S6 (protein I) and two other particulate proteins with apparent molecular masses of 24 kDa (protein II) and 22 kDa (protein III) [Jahn, R., Unger, C. & Söling, H. D. (1980) Eur. J. Biochem. 112, 345-352]. This report describes the purification and characterization of protein III. Solubilization studies indicate that protein III is an intrinsic membrane protein. It could be extracted from the endoplasmic reticulum membrane only with Triton X-100, SDS or concentrated formic or acetic acid. The purification of this protein involved extraction of the microsomes with Triton X-100, removal of the detergent by acetone precipitation, extraction of water-soluble proteins, lipids and lipoproteins, and preparative SDS polyacrylamide gel electrophoresis. The protein has a basic pI (greater than 8.7). For determination of the amino acid composition of protein III and for sequencing of its amino-terminal portion, the protein was electroeluted out off the gel, the detergent removed and the protein finally purified by reversed-phase HPLC. Protein III could be phosphorylated in vitro by the catalytic subunit of the cAMP-dependent protein kinase to a degree of approximately 0.14 mol phosphate/mol protein. The only phosphopeptide obtained after in vitro phosphorylation and subsequent tryptic or chymotryptic digestion was identical with the phosphopeptide obtained after stimulation of intact rat parotid gland lobules with isoproterenol. The sequence of this peptide was Lys-Leu-Ser(P)-Glu-Ala-Asp-Asn-Arg. It was confirmed by an analysis of the synthetic peptide following in vitro phosphorylation with cAMP-dependent protein kinase. The first 41 N-terminal residues of protein III were sequenced. So far no sequence homology with other known peptides or proteins could be found.

On the disulphide bonds of rhodopsins

Carboxymethylation using 14C- or 3H-labelled iodoacetic acid has been used to identify the cysteine residues in bovine rhodopsin involved in the formation of the two intramolecular disulphide bridges. Iodo[2-14C]acetic acid was used to modify 5.8-5.9 residues of cysteine under non-reducing conditions. After dialysis and reduction of disulphide bridges by 2-mercaptoethanol, iodo[2-3H]acetic acid was employed to covalently modify 3.3-3.6 residues of cysteine. Peptide purification and sequencing has unambiguously shown that cysteine residues 322 and 323 are only carboxymethylated after reduction of disulphide bridges. Indirect evidence presented, now coupled with the earlier finding [Findlay & Pappin (1986) Biochem. J. 238, 625-642] suggests that the other disulphide bridge is formed between cysteine residues 110 and 187. A comparison is made of all the sequences of mammalian rhodopsins and colour pigments and attention is drawn to the fact that whereas Cys-322 and Cys-323 are conserved only in three rhodopsins (bovine, ovine and human), the residues corresponding to Cys-110 and Cys-187 are found in all the visual proteins (from rods as well as human cones).

Isolation and characterization of a vinculin cDNA from chick-embryo fibroblasts

A chick-embryo fibroblast lambda gt11 cDNA library was screened with affinity-purified antibodies to chick gizzard vinculin. One recombinant was purified to homogeneity and the fusion protein expressed in Escherichia coli strain C600. The fusion protein was unstable, but polypeptides that reacted with vinculin antibodies, but not non-immune immunoglobulin, were detected by Western blotting. The recombinant contained a single EcoRI fragment of 2891 bp with a single open reading frame. The deduced protein sequence could be aligned with that of six CNBr-cleavage peptides and two tryptic peptides derived from chicken gizzard vinculin. AUG-247 has tentatively been identified as the initiation codon, as it is contained within the consensus sequence for initiation sites of higher eukaryotes. The cDNA lacks 3' sequence and encodes 74% of the vinculin sequence, presuming the molecular mass of vinculin to be 130,000 Da. Analysis of the deduced sequence showed no homologies with other protein sequences, but it does display a triple internal repeat of 112 amino acid residues covering residues 259-589. The sequences surrounding the seven tyrosine residues in the available sequence were aligned with the tyrosine autophosphorylation consensus sequence found in protein tyrosine kinases. Tyr-822 showed a good match to this consensus, and may represent one of the two major sites of tyrosine phosphorylation by pp60v-sre. Northern blots showed that the 2.89 kb vinculin cDNA hybridized to one size of mRNA (approx. 7 kb) in chick-embryo fibroblasts, chick smooth muscle and chick skeletal muscle. Southern blots revealed multiple hybridizing bands in genomic DNA.

Homology between the pyrazine-binding protein from nasal mucosa and major urinary proteins

Sequence analysis of the pyrazine-binding protein from bovine olfactory mucosa reveals marked homology with a family of proteins of unknown function found in the urine of the adult male mouse and rat. In view of the dramatic biological responses to odorants transmitted in male rodent urines, it is proposed that these proteins play important roles in some aspects of odor transmission and reception.

The opsin family of proteins

Images

An opsin gene expressed in only one photoreceptor cell type of the Drosophila eye

We have isolated an opsin gene from D. melanogaster that is expressed specifically in photoreceptor cell 8 of the Drosophila compound eye. This opsin is 381 amino acid residues long and is 67% homologous to the ninaE opsin, which is expressed in photoreceptor cells 1-6. The gene is divided into four exons; only one of the intron positions is conserved with that of the ninaE gene.

Modification of ovine opsin with the photosensitive hydrophobic probe 1-azido-4-[125I]iodobenzene. Labelling of the chromophore-attachment domain

The hydrophobic photosensitive probe 1-azido-4-[125I]iodobenzene (AIB) partitioned preferentially into photoreceptor disc membranes and, upon u.v. irradiation, became covalently bound to opsin and phospholipid. The labelling of both protein and phospholipid was linearly related to AIB concentration. The amount of probe incorporated into protein was not significantly different when membranes were irradiated at -100 degrees, 4 degrees or 25 degrees C, but irreversible aggregation of monomeric opsin was dramatically reduced by performing the photolysis at -100 degrees C. Labelling of opsin after irradiation at -100 degrees or 4 degrees was not significantly reduced by the presence of lysine in the aqueous buffer, indicating that significant amounts of reactive species did not enter the aqueous phase. The incorporation into phospholipid, unlike that into opsin, decreased as the temperature of irradiation increased. Some labelling of opsin occurred on incubation with pre-photoactivated AIB, indicating that reaction may also occur with reactive species of longer lifetimes than the nitrene. Proteolysis of labelled opsin with Staphylococcus aureus V8 proteinase yielded two radiolabelled membrane-bound fragments. The location of the modified sites (cysteine, tryptophan, tyrosine, lysine and histidine residues: all nucleophiles) in the smaller fragment was entirely consistent with putative models for the protein derived from other studies.

Coronavirus IBV: partial amino terminal sequencing of spike polypeptide S2 identifies the sequence Arg-Arg-Phe-Arg-Arg at the cleavage site of the spike precursor propolypeptide of IBV strains Beaudette and M41

The spike protein of avian infectious bronchitis coronavirus comprises two glycopolypeptides S1 and S2 derived by cleavage of a proglycopolypeptide So, the nucleotide sequence of which has recently been determined for the Beaudette strain (Binns, M.M. et al., 1985, J. Gen. Virol. 66, 719-726). The order of the two glycopolypeptides within So is aminoterminus(N)-S1-S2-carboxyterminus(C). To locate the N-terminus of S2 we have performed partial amino acid sequencing on S2 from IBV-Beaudette labelled with [3H]serine and from the related strain labelled with [3H]valine, leucine and isoleucine. The residues identified and their positions relative to the N-terminus of S2 were: serine, 13; valine, 6, 12; leucine, none in the first 20 residues; isoleucine, 2, 19. These results identified the N-terminus of S2 of IBV-Beaudette as serine, 520 residues from the N-terminus of S1, excluding the signal sequence. Immediately to the N-terminal side of residue 520 So has the sequence Arg-Arg-Phe-Arg-Arg; similar basic connecting peptides are a feature of several other virus spike glycoproteins. It was deduced that for IBV-Beaudette S1 comprises 519 residues (Mr 57.0K) or 514 residues (56.2K) if the connecting peptide was to be removed by carboxypeptidase-like activity in vivo while S2 has 625 residues (69.2K). Nucleotide sequencing of the cleavage region of the So gene of IBV-M41 revealed the same connecting peptide as IBV-Beaudette and that the first 20 N-terminal residues of S2 of IBV-M41 were identical to those of the Beaudette strain. IBV-Beaudette grown in Vero cells had some uncleaved So; this was cleavable by 10 micrograms/ml of trypsin and of chymotrypsin. Partial N-terminal analysis of S1 from IBV-M41 identified leucine and valine residues at positions 2 and 9 respectively from the N-terminus. This confirms the identification, made by Binns et al. (1985), of the N-terminus of S1 and the end of the signal sequence of the IBV-Beaudette spike propolypeptide. N-terminal sequencing of [3H]leucine-labelled IBV-Beaudette membrane (M) polypeptide showed leucine residues at positions 8, 16 and 22 from the N-terminus; these results confirm the open reading frame identified by M.E.G. Boursnell et al. (1984, Virus Res. 1, 303-313) in the nucleotide sequence of M. The N-terminus of the nucleocapsid (N) polypeptide appeared to be blocked.

Molecular biology of the visual pigments

With the identification and structural characterization of several visual pigments has come a new era of investigation. The above comparisons of amino acids sequences predict specific functional domains that may be tested to tell us how visual pigments function to absorb light and transform this "signal" to trigger a neural response. The details of how rod and cone pigments differ are now known for human pigments. The striking similarities between vertebrate and invertebrate pigments are remarkable for pigments that have been subject to divergence for over 500 million years. There are yet challenges ahead of us. The true tertiary structure of visual pigments must be obtained from a 3-dimensional crystal structure. The predictions for functional domains of interaction with the GTP binding protein must be confirmed or redefined. A rigorous definition of the chromophore environment and the properties that control the wavelength of absorption of 11-cis retinal chromophore are certainly still on the drawing boards. Specific genetic alteration through in vitro mutagenesis promises much insight, but the technology for expressing these membrane proteins in functional form has yet to be achieved. We may expect, however, these problems will be addressed and in the next few years facts should replace what are now speculations. Finally, it is a delightful observation that nature has capitalized on a general biochemical mechanism for control of second messengers in the cytoplasm of cells. Protein structural data deduced from genetic information now document the hypothesis that the structure and function of receptors for the catecholamines and that of visual pigments are similar. The receptors for serotonin, leukotrienes, prostaglandins, histamine and acetylcholine (muscarinic) are expected to belong to this same family. The lessons learned about visual pigments can be applied broadly to a general set of membrane receptors.

The Drosophila ninaE gene encodes an opsin

The Drosophila ninaE gene was isolated by a multistep protocol on the basis of its homology to bovine opsin cDNA. The gene encodes the major visual pigment protein (opsin) contained in Drosophila photoreceptor cells R1-R6. The coding sequence is interrupted by four short introns. The positions of three introns are conserved with respect to positions in mammalian opsin genes. The nucleotide sequence has intermittent regions of homology to bovine opsin coding sequences. The deduced amino acid sequence reveals significant homology to vertebrate opsins; there is strong conservation of the retinal binding site and two other regions. The predicted protein secondary structure strikingly resembles that of mammalian opsins. We conclude the Drosophila and vertebrate opsin genes are derived from a common ancestor.

Phosphorylation of ovine rhodopsin. Identification of the phosphorylated sites

Light-dependent phosphorylation of sheep opsin was obtained in purified discs to which was added a partially purified preparation of rhodopsin kinase. A maximum ratio of 1.8 mol of phosphate/mol of rhodopsin bleached was obtained. Perturbing the lipid bilayer did not alter the phosphorylation ratio. Dephosphorylation in both segments and discs was only achieved when the supernatant fraction from a retina homogenate was added. Complete dephosphorylation required the presence of the detergent dodecyltrimethylammonium bromide in the incubation medium. Treatment of phosphorylated disc membranes with Staphylococcal aureus V8 proteinase generated two membrane-bound fragments, only one of which (V8-S, Mr 12 000) was labelled, together with a soluble seven-residue peptide that contained [32P]phosphoserine. Peptide sequencing, together with subdigestion procedures, localized the phosphorylation sites to serine residues at positions 334, 338 and 343 in the whole sequence and threonine residues at positions 335 and 336.

Labelling of the cytoplasmic domains of ovine rhodopsin with hydrophilic chemical probes

The disposition of polypeptide chain of ovine rhodopsin in the photoreceptor disc membrane was investigated by using two hydrophilic reagents, 3,5-di-[125I]iodo-4-diazobenzenesulphonate [( 125I]DDISA) and [14C]succinic anhydride. Both reagents were used to modify rhodopsin in intact disc membranes under conditions where no loss of A500 occurred. Reaction of [125I]DDISA with rhodopsin approached completion after 30 min. Binding was saturated at a 75-fold molar excess of reagent, which gave binding ratios of up to 2 mol/mol of rhodopsin. Proteolysis of rhodopsin, using Staphylococcus aureus V8 proteinase, yielded two membrane-bound fragments, both of which contained bound radioactive probe. Subsequent CNBr cleavage of these fragments produced five radiolabelled peptides which corresponded to the C-terminal region and cytoplasmic loops of rhodopsin. Similar studies with [14C]-succinic anhydride also gave binding ratios of up to 2 mol/mol of rhodopsin. Sequencing of the [14C]succinylated peptides identified the location of the reactive sites as lysine residues 66, 67, 141, 245, 248, 311, 325 and 339 in the polypeptide chain. Non-permeability of both probes was demonstrated by the absence of any radioactivity associated with the intradiscal N-terminal glycopeptide. Sonication of membranes in the presence of [125I]DDISA led to the incorporation of label in this peptide.

The structure of mammalian rod opsins

Ovine rhodopsin is organised in disc membranes as a monomer. The determination of its amino acid sequence has permitted the utilisation of structure prediction programmes which indicate the probable disposition of the polypeptide chain in the bilayer. This putative model is consistent with labelling data using the chemical probes, [14C]succinic anhydride, [125I]diazodiido sulphanilic acid and [125I]iodophenyl azide, and with the cleavage points for several proteases. More surprisingly the predicted structure points to the occurrence of breaks/distortions in the transmembrane helical segments. These distorted regions may be of primary functional importance to the protein and at least one is associated with the attachment point of the chromophore. This particular part of the structure is also identified as a "mutational hot spot", for bovine, equine, ovine and porcine opsins exhibit different sequences (but conserved molecular volumes) in the four residues following the retinyllysine. In an otherwise highly conserved protein with no obvious functional differences between the four species, the high substitution rate in this region is unexplained.