Evolution of the Integrin α and β Protein Families

A phylogenetic analysis of vertebrate and invertebrate α integrins supported the hypothesis that two major families of vertebrate α integrins originated prior to the divergence of deuterostomes and protostomes. These two families include, respectively, the αPS1 and αPS2 integrins of Drosophila melanogaster, and each family has duplicated repeatedly in vertebrates but not in Drosophila. In contrast, a third family (including αPS3) has duplicated in Drosophila but is absent from vertebrates. Vertebrate αPS1 and αPS2 family members are found on human chromosomes 2, 12, and 17. Linkage of these family members may have been conserved since prior to the origin of vertebrates, and the two genes duplicated simultaneously. A phylogenetic analysis of β integrins did not clearly resolve whether vertebrate β integrin genes duplicated prior to the origin of vertebrates, although it suggested that at least the gene encoding vertebrate β4 may have done so. In general, the phylogeny of neither α nor β integrins showed a close correspondence with patterns of α–β heterodimer formation or other functional characteristics. One major exception to this trend involved αL, αM, αX, and αD, a monophyletic group of immune cell-expressed α integrins, which share a number of common functional characteristics and have evolved in coordinated fashion with their β integrin partners. Received: 22 June 2000 / Accepted: 11 September 2000     

Effect of Oxidation of αA- and αB-Crystallins on their Structure, Oligomerization and Chaperone Function

The purpose of this study was to investigate the effect of metal-catalyzed oxidation by H₂O₂ on the structure, oligomerization, and chaperone function of αA- and αB-crystallins. Recombinant αA-and αB-crystallins were prepared by expressing them in E. coli and purifying by size-exclusion chromatography. They were incubated with 1.5 mM H₂O₂ and 0.1 mM FeCl₃ at 37 ^∘C for 24 hrs and the reaction was stopped by adding catalase. Structural changes due to oxidation were ascertained by circular dichroism (CD) measurements and chaperone activity was assayed with alcohol dehydrogenase (ADH) and insulin as target proteins. The oligomeric nature of the oxidized proteins was assessed by molecular sieve HPLC. The secondary structure of the oxidized αA- and αB-crystallins has been substantially altered due to significant increase in random coils, in addition to decrease in β-sheet or α-helix contents. The tertiary structure also showed significant changes indicative of different mode of folding of the secondary structural elements. Chaperone function was significantly compromised as supported by nearly 50% loss in chaperone activity. Oxidation also resulted in the formation of higher molecular weight (HMW) proteins as well as lower molecular weight (LMW) proteins. Thus, oxidation leads to disintegration of the oligomeric structure of αA- and αB-crystallins. Chaperone activity of the HMW fraction is normal whereas the LMW fraction lacks any chaperone activity. So, it appears that the formation of the LMW proteins is the primary cause of the chaperone activity loss due to oxidation.     

Analysis of Chameleon Sequences by Energy Decomposition on a Pairwise Per-residue Basis

The conversion from α-helix to β-strand that has been widely observed in so-called chameleon sequences has received considerable attention since such a structural change may induce many amyloidogenic proteins to self-assemble into fibrils thus causing fatal diseases. Here we report a large scale-analysis of the energetics of secondary structural conversions in a collection of chameleon sequences retrieved from the Protein Data Bank. Major energetic contributions to the secondary structural conversion were analyzed by carrying out energy decomposition on a pairwise per-residue basis, i.e., (i,i), (i,i ± 1), (i,i ± 2), (i,i ± 3), (i,i ± 4) and > (i,i ± 4) intra-/inter-residual interactions. While the overall potential energy differences were subtle, individual residue-based interacting energy differences were observed to vary significantly depending on the specific type of secondary structural conversion. The average energy difference between α-helix and β-strand, <ΔE _α→β>, in the chameleon sequences varied significantly in (i,i), (i,i ± 1) and > (i,i ± 4) interactions. The major energetic factors in secondary structure conversions were electrostatic interactions and the polar term for solvation energy. In addition, residue-based average energy differences in α-helix → β-strand conversions were well-correlated to those in α-helix → random coil → β-strand conversions (R ² = 0.92). Assuming that three secondary structural elements can transform in either direction, this strong correlation indicates that the present energy decomposition method using database structures of chameleon sequences provides a reliable tool for the characterization of secondary structure fluctuations in amino acid sequences.     

Cloning of α-amylase gene from Schwanniomyces occidentalis and expression in Saccharomyces cerevisiae

The cloning of α-amylase gene ofS. occidentalis and the construction of starch digestible strain of yeast,S. cerevisiae AS. 2. 1364 with ethanol-tolerance and without auxotrophic markers used in fermentation industry were studied. The yeast/E.coli shuttle plasmid YCEp1 partial library ofS. occidentalis DNA was constructed and α-amylase gene was screened in S.cerevisiae by amylolytic activity. Several transformants with amylolysis were obtained and one of the fusion plasmids had an about 5.0 kb inserted DNA fragment, containing the upstream and downstream sequences of α-amylase gene fromS. occidentalis. It was further confirmed by PCR and sequence determination that this 5.0 kb DNA fragment contains the whole coding sequence of α-amylase. The amylolytic test showed that when this transformant was incubated on plate of YPDS medium containing 1 % glum and 1 % starch at 30°C for 48 h starch degradation zones could be visualized by staining with iodine vapour. α-amylase activity of the culture filtratate is 740–780 mU/mL and PAGE shows that the yeast harboring fusion plasmids efficiently secreted α-amylase into the medium, and the amount of the recombinant α-amylase is more than 12% of the total proteins in the culture filtrate. These results showed that α-amylase gene can be highly expressed and efficiently secreted inS. cerevisiae AS. 2.1364, and the promotor and the terminator of α-amylase gene fromS. occidentalis work well inS. cercvisiac AS. 2.1364.     

A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences—exons in vertebrates and genes in prokaryotes—are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content. Correspondence to: J. L. Oliver     

Pathway of Oxidative Folding of a 3-Disulfide α-Lactalbumin May Resemble Either BPTI Model or Hirudin Model

Pathways of oxidative folding of disulfide proteins display a high degree of diversity and vary among two extreme models. The BPTI model is defined by limited species of folding intermediates adopting mainly native disulfide bonds. The hirudin model is characterized by highly heterogeneous folding intermediates containing mostly non-native disulfide bonds. αLA-IIIA is a 3-disulfide variant of α-lactalbumin (αLA) with a 3-D conformation essentially identical to that of intact αLA. αLA-IIIA contains 3 native disulfide bonds of αLA, two of them are located at the calcium binding β-subdomain (Cys⁶¹–Cys⁷⁷ and Cys⁷³–Cys⁹¹) and the third bridge is located within the α-helical domain of the molecule (Cys²⁸–Cys¹¹¹). We investigate here the pathway of oxidative folding of fully reduced αLA-IIIA with and without stabilization of its β-subdomain by calcium binding. In the absence of calcium, the folding pathway of αLA-IIIA was shown to resemble that of hirudin model. Upon stabilization of β-sheet domain by calcium binding, the folding pathway of αLA-IIIA exhibits a striking similarity to that of BPTI model. Three predominant folding intermediates of αLA-IIIA containing exclusively native disulfide bonds were isolated and structurally characterized. Our results further demonstrate that stabilization of subdomains in a protein may dictate its folding pathway and represent a major cause for the existing diversity in the folding pathways of the disulfide-containing proteins.     

Identifying a Folding Nucleus for the Lysozyme/α-Lactalbumin Family from Sequence Conservation Clusters

Four subfamilies of c-type lysozyme and one subfamily of α-lactalbumin are defined from 78 sequences, and their folding nucleus is identified with a method based on conserved residues and native structural contacts between pairs of conserved residues. One large cluster of 19 conserved residues is found which is mostly nonpolar, buried, and nonfunctional. It can be subdivided into three subclusters: (1) conserved residues in four helices; (2) conserved residues that stabilize the connector between the α and the β domains; and (3) a β-turn, sitting in the middle of a bowl of α-helix residues. It is proposed that this folding nucleus initiates four helices, A, B, C, and D, three β sheets, and the connector, which corresponds closely to the nucleation of the so-called fast folding track pathway. As the secondary structures propagate, nonconserved residues and functionally conserved residues would form additional contacts. The conserved residues are selected with a phylogenetic scheme in which single members of subfamilies are selected. Subfamilies are then equally weighted to obtain the consensus conservation. Received: 11 June 2001 / Accepted: 28 August 2001     

Construction and characterization of the hetero-oligomer of the group II chaperonin from the hyperthermophilic archaeon, Thermococcus sp. strain KS-1

The hyperthermophilic archaeon Thermococcus sp. strain KS-1 (T. KS-1) expresses two different chaperonin subunits, α and β, for the folding of its proteins. The composition of the subunits in the hexadecameric double ring changes with temperature. The content of the β subunit significantly increases according to the increase in temperature. The homo-oligomer of the β subunit, Cpnβ, is more thermostable than that of the α subunit, Cpnα. Since Cpnα and Cpnβ also have different protein folding activities and interactions with prefoldin, the hetero-oligomer is thought to exhibit different characteristics according to the content of subunits. The hetero-oligomer of the T. KS-1 chaperonin has not been studied, however, because the α and β subunits form hetero-oligomers of varying compositions when they are expressed simultaneously. In this study, we characterized the T. KS-1 chaperonin hetero-oligomer, Cpnαβ, containing both α and β in the alternate order, which was constructed by the expression of α and β subunits in a coordinated fashion and protease digestion. Cpnαβ protected citrate synthase from thermal aggregation, promoted the folding of acid-denatured GFP in an ATP-dependent manner, and exhibited an ATP-dependent conformational change. The yield of refolded GFP generated by Cpnαβ was almost equivalent to that generated by Cpnβ but lower than that generated by Cpnα. In contrast, Cpnαβ exhibited almost the same level of thermal stability as Cpnα, which was lower than that of Cpnβ. The affinity of Cpnαβ to prefoldin was found to be between those of Cpnα and Cpnβ, as expected.     

Effects of mutations in the WD40 domain of α-COP on its interaction with the COPI coatomer in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis>

Replacement of glycine 227 in the fifth WD40 motif of α-COP/Ret1p/Soo1p by charged or aromatic amino acids is responsible for the temperature-dependent osmo-sensitivity of Saccharomyces cerevisiae, while truncations of WD40 motifs exerted a reduction in cell growth rate and impairment in assembly of cell-wall associated proteins such as enolase and Gas1p. Yeast two-hybrid analysis revealed that the ret1-1/soo1-1 mutation of α-COP abolished the interaction with β- and ɛ-COP, respectively, and that the interaction between α-COP and β-COP relied on the WD40 domain of α-COP. Furthermore, although the WD40 domain is dispensable for interaction of α-COP with ɛ-COP, structural alterations in the WD40 domain could impair the interaction.     

Non-native hydrophobic interactions detected in unfolded apoflavodoxin by paramagnetic relaxation enhancement

Transient structures in unfolded proteins are important in elucidating the molecular details of initiation of protein folding. Recently, native and non-native secondary structure have been discovered in unfolded A. vinelandii flavodoxin. These structured elements transiently interact and subsequently form the ordered core of an off-pathway folding intermediate, which is extensively formed during folding of this α–β parallel protein. Here, site-directed spin-labelling and paramagnetic relaxation enhancement are used to investigate long-range interactions in unfolded apoflavodoxin. For this purpose, glutamine-48, which resides in a non-native α-helix of unfolded apoflavodoxin, is replaced by cysteine. This replacement enables covalent attachment of nitroxide spin-labels MTSL and CMTSL. Substitution of Gln-48 by Cys-48 destabilises native apoflavodoxin and reduces flexibility of the ordered regions in unfolded apoflavodoxin in 3.4 M GuHCl, because of increased hydrophobic interactions in the unfolded protein. Here, we report that in the study of the conformational and dynamic properties of unfolded proteins interpretation of spin-label data can be complicated. The covalently attached spin-label to Cys-48 (or Cys-69 of wild-type apoflavodoxin) perturbs the unfolded protein, because hydrophobic interactions occur between the label and hydrophobic patches of unfolded apoflavodoxin. Concomitant hydrophobic free energy changes of the unfolded protein (and possibly of the off-pathway intermediate) reduce the stability of native spin-labelled protein against unfolding. In addition, attachment of MTSL or CMTSL to Cys-48 induces the presence of distinct states in unfolded apoflavodoxin. Despite these difficulties, the spin-label data obtained here show that non-native contacts exist between transiently ordered structured elements in unfolded apoflavodoxin.     

Display of Candida antarctica lipase B on Pichia pastoris and its application to flavor ester synthesis

Two alternative cell-surface display systems were developed in Pichia pastoris using the α-agglutinin and Flo1p (FS) anchor systems, respectively. Both the anchor cell wall proteins were obtained originally from Saccharomyces cerevisiae. Candida antarctica lipase B (CALB) was displayed functionally on the cell surface of P. pastoris using the anchor proteins α-agglutinin and FS. The activity of CALB displayed on P. pastoris was tenfold higher than that of S. cerevisiae. The hydrolytic and synthetic activities of CALB fused with α-agglutinin and FS anchored on P. pastoris were investigated. The hydrolytic activities of both lipases displayed on yeast cells surface were more than 200 U/g dry cell after 120 h of culture (200 and 270 U/g dry cell, respectively). However, the synthetic activity of CALB fused with α-agglutinin on P. pastoris was threefold higher than that of the FS fusion protein when applied to the synthesis of ethyl caproate. Similarly, the CALB displayed on P. pastoris using α-agglutinin had a higher catalytic efficiency with respect to the synthesis of other short-chain flavor esters than that displayed using the FS anchor. Interestingly, for some short-chain esters, the synthetic activity of displaying CALB fused with α-agglutinin on P. pastoris was even higher than that of the commercial CALB Novozyme 435.     

New proteins orthologous to cerato-platanin in various <Emphasis Type="Italic">Ceratocystis</Emphasis> species and the purification and characterization of cerato-populin from <Emphasis Type="Italic">Ceratocystis populicola</Emphasis>

Natural variants of cerato-platanin (CP), a pathogen associated molecular pattern (PAMP) protein produced by Ceratocystis platani (the causal agent of the plane canker stain), have been found to be produced by other four species of the genus Ceratocystis, including five clones of Ceratocystis fimbriata isolated from different hosts. All these fungal strains were known to be pathogenic to plants with considerable importance in agriculture, forestry, and as ornamental plants. The putative premature proteins were deduced on the basis of the nucleotide sequence of genes orthologous to the cp gene of C. platani; the deduced premature proteins of Ceratocystis populicola and Ceratocystis variospora reduced the total identity of all the others from 87.3% to 60.3%. Cerato-populin (Pop1), the CP-orthologous protein produced by C. populicola, was purified and characterized. Pop1 was a well-structured α/β protein with a different percentage of the α-helix than CP, and it self-assembled in vitro in ordered aggregates. Moreover, Pop1 behaved as PAMP, since it stimulated poplar leaf tissues to activate defence responses able to reduce consistently the C. populicola growth. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Molecular dynamics simulations of the folding of poly(alanine) peptides

The secondary structures and the shapes of long-chain polyalanine (PA) molecules were investigated by all-atom molecular dynamics simulations using a modified Amber force field. Homopolymers of polyaminoacids such as PA are convenient models to study the mechanism of protein folding. It was found that the conformational structures of PA peptides are highly sensitive to the chain length. In the absence of solvent, straight α-helices dominate in short (n ∼ 20) peptides at room temperature. A shape transition occurs at a chain length n of 40–45; the compact helix-turn-helix structure (the double-leg hairpin) becomes favored over a straight α-helix. For n = 60, double-leg and the triple-leg hairpins are the only structures present in PA molecules. An exploration of a chain organization in a cubic cavity revealed a clear predisposition of PA molecules for additional breaks in α-helices and the formation of multifolded hairpins. Furthermore, under confinement the hairpin structure becomes much looser, the antiparallel positions of helical stems are disturbed, and a sizeable proportion of the helical stems are transformed from α-helices into 3₁₀-helices.     

Evolution of cis-Acting Elements in 5′ Flanking Regions of Vertebrate Actin Genes

Regulation of the vertebrate actin multigene family involves the recognition of various regulatory sequences (cis-acting elements) that specify the distinct tissue type and developmental program of expression for each actin paralogue, which implies that the distribution of cis-acting elements may be unique for each paralogue gene. To elucidate the evolution of these unique distribution patterns, we improved a method to scan for cis-acting elements in the 5′ flanking regulatory region of genes and used it to analyze five cis-acting elements (SRE, MyoD binding site, Elk-1 binding site, positive and negative YY1 binding sites) of six actin paralogue genes (β and γ cytoplasmic actins, α and γ smooth muscle actins, and α skeletal and α cardiac actins) among various vertebrates. It was shown that although an element(s) may exist in all paralogue genes of the same species, its numbers, compositions, and distribution patterns or even sequences vary remarkably among paralogues, which contributes to their different tissue- and developmental-specific expression. However, each pair of coexpressed paralogues has some certain similarity in distribution patterns. Furthermore, among various orthologues of actin genes derived from diverse vertebrates, the sequences, numbers, and distribution patterns of these cis-acting elements are highly conserved or even identical in the long run of phylogeny of vertebrates. Taken together, the results described above strongly indicate that not only the structures of actins but also their expression patterns are essential in both the phylogeny and the physiology of vertebrates. The distribution patterns of cis-acting elements of various actin genes can be regarded as indicators of both horizontal (paralogous) and vertical (orthologous) evolution of actins. Received: 1 March 1999 / Accepted: 6 August 1999     

Evidence for a Diverse Cys-Loop Ligand-Gated Ion Channel Superfamily in Early Bilateria

The genome sequences of Caenorhabditis elegans and Drosophila melanogaster reveal a diversity of cysteine-loop ligand-gated ion channels (Cys-loop LGICs) not found in vertebrates. To better understand the evolution of this gene superfamily, I compared all Cys-loop LGICs from rat, the primitive chordate Ciona intestinalis, Drosophila, and C. elegans. There are two clades of GABA receptor subunits that include both verterbate and invertebrate orthologues. In addition, I identified nine clades of anion channel subunits found only in invertebrates, including three that are specific to C. elegans and two found only in Drosophila. One well-defined clade of vertebrate cation channel subunits, the α7 nicotinic acetylcholine receptor subunits (nAChR), includes invertebrate orthologues. There are two clades of invertebrate nAChRs, one of α-type subunits and one of non-α subunits, that are most similar to the two clades of vertebrate neuronal and muscle α and non-α subunits. There is a large group of divergent C. elegans nAChR-like subunits partially resolved into clades but no orthologues of 5HT3-type serotonin receptors in the invertebrates. The topology of the trees suggests that most of the invertebrate-specific Cys-loop LGIC clades were present in the common ancestor of chordates and ecdysozoa. Many of these disappeared from the chordates. Subsequently, selected subunit genes expanded to form large subfamilies. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Rafael Zardoya]     

Four hydrophobic amino acid residues in the C-terminal effector domain of the yeast Mig1p repressor are important for its in vivo activity

The Mig1 repressor is a zinc finger protein that mediates glucose repression in yeast. Previous work in Saccharomyces cerevisiae has shown that two domains in Mig1p are required for repression: the N-terminal zinc finger region and a C-terminal effector domain. Both domains are also conserved in Mig1p homologs from the distantly related yeasts Kluyveromyces lactis and K. marxianus, and these Mig1 proteins can fully replace the endogenous Mig1p in S. cerevisiae. We have now made a detailed analysis of the conserved C-terminal effector domain in Mig1p from K. marxianus, using expression in S. cerevisiae to monitor its function. First, a series of small deletions were made within the effector domain. Second, an alanine scan mutagenesis was carried out across the effector domain. Third, double, triple and quadruple mutants were made that affect certain residues within the effector domain. Our results show that four conserved residues within the effector domain, three leucines and one isoleucine, are particularly important for its function in vivo. The analysis further revealed that while the C-terminal effector domain of KmMig1p mediates a seven- to nine-fold repression of the reporter gene, a five- to sixfold residual effect also exists that is independent of the C-terminal effector domain. Similar results were obtained when the corresponding mutations were made in ScMig1p. Moreover, we found that mutations in these residues affect the interaction between Mig1p and the general corepressor subunit Cyc8p (Ssn6p). Modeling of the C-terminal effector domain using a protein of known structure suggests that it may be folded into an α-helix. Received: 30 March 1998 / Accepted: 18 August 1998     

<Emphasis Type="Italic">Pichia pastoris</Emphasis> is a valuable host for the expression of genes encoding membrane proteins from the hyperthermophilic Archeon <Emphasis Type="Italic">Pyrococcus abyssi</Emphasis>

We present here the experimental strategies, first results and identified bottlenecks of a structural genomics initiative on membrane proteins of the hyperthermophilic archaea Pyrococcus abyssi. Five ORFs coding for putative membrane proteins have been cloned and expressed in the methylotrophic Pichia pastoris expression system, using two different constructs, with or without the signal sequence α-mating factor of Saccharomyces cerevisiae. A c-myc epitope and 6 His codons were added at the 3′-end of the targeted genes to allow immunodetection of the recombinant proteins and to facilitate their further purification. We have selected at least one producer clone for each protein of interest and for almost every construction. All the membrane proteins were produced in Erlenmeyer flasks culture and in fed-batch cultivation for large-scale preparation. The proteins were detected in the membrane fractions of P. pastoris. Production efficiencies were relatively low in both production conditions but the quantities of biomass obtained during fed-batch cultivation have allowed us to collect sufficient amount of material for further purification. The proteins were extracted, solubilized and partially purified. Large-scale purification will be necessary for further structural work.     

Structural comparison between wild-type and P25S human cystatin A by NMR spectroscopy. Does this mutation affect the α-helix conformation ?

The effect of substituting Pro²⁵, located in the α-helical region of the cystatin A structure, with Ser has been studied. The structures of wild type and P25S cystatin A were determined by multidimensional NMR spectroscopy under comparable conditions. These two structures were virtually identical, and the α-helix between Glu¹⁵-Lys³⁰ exists with uninterrupted continuity, with a slight bend at residue 25. In order to characterize the possible substitution effects of Pro²⁵ with Ser on the α-helix, the chemical shifts of the amide nitrogens and protons, the generalized order parameters obtained by the analyses of the ¹⁵N-¹H relaxation data, the amide proton exchange rates, and the NOE networks among the α-helical and surrounding residues were carefully compared. None of these parameters indicated any significant static or dynamic structural differences between the α-helical regions of the wild-type and P25S cystatin A proteins. We therefore conclude that our previous structure of the wild-type cystatin A, in which the α-helix exhibited a sharp kink at Pro²⁵, must be revised. The asymmetric distribution of hydrophobic interactions between the side-chain residues of the α-helix and the rolled β-sheet surface, as revealed by NOEs, may be responsible for the slight bend of the α-helix in both variants and for the destabilized hydrogen bonding of the α-helical residues that follow Pro²⁵/Ser²⁵, as evidenced by increased amide exchange rates. This revised version was published online in July 2006 with corrections to the Cover Date.     

Specific Selection Pressure at the Third Codon Positions: Contribution to 10- to 11-Base Periodicity in Prokaryotic Genomes

Prokaryotic sequences are responsible for more than just protein coding. There are two 10- to 11-base periodical patterns superimposed on the protein coding message within the same sequence. Positional auto- and cross-correlation analysis of the sequences shows that these two patterns are a short-range counter-phase oscillation of AA and TT dinucleotides and a medium-range in-phase oscillation of the same dinucleotides, spanning distances of up to ∼30 and ∼100 bases, respectively. The short-range oscillation is encoded by the amino acid sequences themselves, apparently, due to the presence of amphipathic α-helices in the proteins. The medium-range oscillation, related to DNA folding in the cell, is created largely by a special choice of the bases in the third positions of the codons. Interestingly, the amino acid sequences do contribute to that signal as well. That is, the very amino acid sequences are, to some extent, degenerate to serve the same oscillating pattern that is associated with the degenerate third codon positions. [Reviewing Editor: Dr. Richard Kliman]