Skip Navigation


MBE Advance Access originally published online on August 16, 2007
Molecular Biology and Evolution 2007 24(11):2424-2432; doi:10.1093/molbev/msm171
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrowOA All Versions of this Article:
24/11/2424    most recent
msm171v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Google Scholar
Right arrow Articles by Sutherland, T. D.
Right arrow Articles by Haritos, V. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sutherland, T. D.
Right arrow Articles by Haritos, V. S.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2007 The Authors.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.


Research Articles

Conservation of Essential Design Features in Coiled Coil Silks

Tara D. Sutherland*, Sarah Weisman*, Holly E. Trueman*, Alagacone Sriskantha*, John W. H. Trueman{dagger} and Victoria S. Haritos*

* CSIRO Entomology, Canberra, ACT, Australia
{dagger} School of Botany and Zoology, Australian National University, Canberra, ACT, Australia

E-mail: tara.sutherland{at}csiro.au.


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Silks are strong protein fibers produced by a broad array of spiders and insects. The vast majority of known silks are large, repetitive proteins assembled into extended ß-sheet structures. Honeybees, however, have found a radically different evolutionary solution to the need for a building material. The 4 fibrous proteins of honeybee silk are small (~30 kDa each) and nonrepetitive and adopt a coiled coil structure. We examined silks from the 3 superfamilies of the Aculeata (Hymenoptera: Apocrita) by infrared spectroscopy and found coiled coil structure in bees (Apoidea) and in ants (Vespoidea) but not in parasitic wasps of the Chrysidoidea. We subsequently identified and sequenced the silk genes of bumblebees, bulldog ants, and weaver ants and compared these with honeybee silk genes. Each species produced orthologues of the 4 small fibroin proteins identified in honeybee silk. Each fibroin contained a continuous predicted coiled coil region of around 210 residues, flanked by 23–160 residue length N- and C-termini. The cores of the coiled coils were unusually rich in alanine. There was extensive sequence divergence among the bee and ant silk genes (<50% similarity between the alignable regions of bee and ant sequences), consistent with constant and equivalent divergence since the bee/ant split (estimated to be 155 Myr). Despite a high background level of sequence diversity, we have identified conserved design elements that we propose are essential to the assembly and function of coiled coil silks.

Key Words: silk • coiled coils • bees • ants • social insects


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Silks are integral to invertebrate biology, providing building material for structures ranging from prey capture webs to protective cocoons for vulnerable juvenile stages. Silks have also attracted much attention in recent years due to their excellent mechanical properties as new generation biomaterials (Vollrath 2000Go; Scheller et al. 2001Go; Vollrath and Knight 2001Go; Lazaris et al. 2002Go). The best studied silks are the cocoon fibers of silkworms and related moths and the webs spun by orb-weaving spiders. Despite the evolutionary distance between these groups, their silk proteins share a number of convergent design principles (Bini et al. 2004Go). The architecture of the large silk proteins (>150 kDa) is an alternation of long hydrophobic blocks and short hydrophilic spacers, flanked by hydrophilic termini. The proteins are highly soluble in the gland due to their hydrophilic moieties and then are processed by hydrophobic aggregation into insoluble fibers with ß-crystalline structure. The primary sequences of the silk proteins excluding their termini consist of arrays of highly conserved repeat motifs.

It has long been known that several species of Hymenoptera, including honeybees, hornets, and solitary wasps, produce silk with an entirely different molecular structure to the ß-sheet silks described above. The cocoon silk of these species consists of coiled coils that run parallel to the fiber axis (Rudall 1962Go; Atkins 1967Go). Coiled coils are a protein-structural arrangement where multiple {alpha}-helices wind around each other. Their protein sequences are characterized by a 7-residue periodicity where large, nonpolar residues generally occupy the first and fourth position of each heptad (Woolfson 2005Go). This coiled coil silk structure is apparently confined to aculeate species with the silks of the nonaculeate Hymenoptera (parasitic wasps and sawflies) generally composed of proteins in predominantly a ß-sheet conformation (Rudall 1962Go). The paraphyly of parasitic wasps and sawflies with respect to the aculeate Hymenoptera suggests that coiled coil silks represent the derived condition within this insect order with the ß-sheet silks being plesiomorphic.

Recently, honeybee silk has been shown to consist of 4 small (~30 kDa) proteins that contain extensive coiled coil motifs (Sutherland et al. 2006Go). These proteins are paralogues but their sequence has diverged extensively so that they retain only around 30% identity. Within the proteins there are no primary sequence repeats. Clearly, the honeybee silks are based on completely different design principles to moth and spider silks. We analyzed silks produced by a range of members of the Aculeata to determine the distribution of coiled coil silks. Among the Aculeata producing coiled coil silk, we selected exemplar species separated by the greatest evolutionary distance. We obtained the silk genes of these species and compared the sequences of the encoded silk in order to identify conserved features that may be necessary for the assembly and function of coiled coil silk.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Insects and Tissues
Honeybee (Apis mellifera) larvae and silk were obtained from a commercial hive. Bumblebee (Bombus terrestris) larvae and silk were obtained from Biobees (New Zealand). Weaver ant (Oecophylla smargdina) larvae and silk were collected from Darwin, Northern Territory, Australia. Bulldog ant (Myrmecia forficata) larvae and silk were collected from the Tinderry Nature Reserve outside Canberra, Australia. Goniozus sp. wasp silk was kindly provided by Santhosh S. Nair, University of Calicut, and Prosierola sp. wasp silk was kindly provided by Robert Wharton, Texas A&M University. The labial gland from larvae in the final instar was dissected under phosphate-buffered saline, immediately transferred to RNAlater (Ambion, Austin, TX), and subsequently stored at 4 °C.

Fourier Transform Infrared Spectroscopy Methods
The bee silks were washed extensively in chloroform to remove wax. All silks were washed extensively in warm water containing 2% sodium dodecyl sulfate to remove water-soluble contaminants and then rinsed in distilled water. Cleaned and dried silk samples were examined using a Bruker Tensor 37 Fourier transform infrared spectrometer with a Pike Miracle diamond attenuated total reflection accessory. The amide I region of the spectra was analyzed to identify the protein secondary structure. Spectra from samples that were degummed (boiled for 30 min in 0.05% sodium carbonate solution, a standard procedure for degumming silkworm silk) were unchanged.

cDNA Library Construction
Total RNA was isolated from 4 bumblebee (2 µg), 4 bulldog ant (3 µg), and ~100 weaver ant (0.4 µg) late-larval labial glands using the RNAqueous4PCR kit from Ambion. mRNA was isolated from the total RNA using the Micro-FastTrack 2.0 mRNA Isolation kit from Invitrogen (Carlsbad, CA) into a final volume of 10 µl water. cDNA libraries were constructed from the mRNA using the CloneMiner cDNA kit of Invitrogen with amounts of reagents reduced from that described in the standard protocol to accommodate the small amounts of starting material (as described in Sutherland et al. 2006Go). The cDNA libraries comprised ~2 x 107 (bumblebee), 5 x 107 (bulldog ant), and 5,000 (weaver ant) colony-forming units with less than 1% original vector for the bulldog ant and bumblebee libraries and greater than 80% original vector in the weaver ant library. The average insert size within the libraries was 1.3 kbp.

Quantitative Reverse Transcriptase–Polymerase Chain Reaction Analysis
The cDNA library ligation mix was analyzed by quantitative reverse transcriptase–polymerase chain reaction (RT-PCR) in the Applied Biosystems 7,000 sequencing detection machine using iTaq plus ROX master mix (BioRad, Hercules, CA) according to manufacturer's instructions using the following primer pairs (shown in 5'–3' direction)—honeybee: AmelF1 TCGCCGAGATTGATAAAAAG, GCGGCACGCGCATTA; AmelF2 TAACGCCGATTCAATCATCAC, TTGCCGCCTTAACTTGAGCTT; AmelF3 GGGATGGAGCGATTATAGGACTT, CGCCTGTGCAAGCAACTG; AmelF4 TGGCAGGTGAGATGAAACCA, TCTTCCTTCGTTGACGTTTTCC; bumblebee: BBF1 CGCATCCGCGATGTCAT, TGATTCTTGGCCGTCTTTTCA; BBF2 CGGCCGCCGGTTTC, GTTTAGTGGCCTCGGATTCG; BBF3 CAACGCAGGCACAGTTGCT, CAGTTGTTGCTCAGATTGAGTTT; BBF4 ACCGCAGTGAACGCCATT, TGCTACACCAGCGGTGTTCTC; bulldog ant: BAF1 CTCGTGCAGGCCAATGC, GCTCTAGACTTCGAATCCAATTCG; BAF2 CAGGCCGCAAGTCAGATACA, GCGGTGTTGGTCTCTGATTGT; BAF3 AATCTTCACGTGTTGCCATTACC, TTGATGTCGAACACCGACTTG; BAF4 AACGCCGATCTGTCCAAGAG, CTGATTTTCTTGACTTCCTGTACTTGA; and weaver ant: WAF1 CAATCAGGAGATCGTTGAAATCG, TGTGCCGGTCGAGATCACT; WAF2 CGACCAGGGCGTTGCT, TTTCTGGGTAAGTTCGGCTTGA; WAF3 CGTTAGCGTCCACCGTCAA, CAGAAGAACTGGCGGTATTTCC; WAF4 AGCTACCGAAGCCGCTTCTAC, GGATTGGTGATCTCGCCTTTC. Primers were designed using the Primer Express software of Applied Biosystems and chosen based on software scores and relative position on the gene. Relative levels of cDNA were determined by comparing a threshold cycle for amplification for the 4 primer pairs at different library concentrations (1/1,000–1/100,000 dilutions) within the same experiment.

Silk Preparation, Amino Acid Analysis and Liquid Chromatography/Mass Spectrometry Methods
The amino acid compositions of clean silks were determined after 24 h gas phase hydrolysis at 110 °C using the Waters AccQTag chemistry by the Australian Proteome Analysis Facility Ltd. (Macquarie University, Sydney). The proteins present in the clean silks were identified by liquid chromatography followed by tandem mass spectrometry (MS/MS) as described previously (Sutherland et al. 2006Go). Briefly, silk was digested with trypsin and the resultant peptides were analyzed by MS/MS . The experimentally observed sets of peptide mass fragments were compared with predicted proteins encoded by cDNAs from the silk gland libraries.

Cladistic analysis and domain predictions
Signal peptides were predicted by SignalP 3.0 (Bendtsen et al. 2004Go), and coiled coil regions were predicted by MARCOIL (Delorenzi and Speed 2002Go; http://www.isrec.isb-sib.ch/webmarcoil/webmarcoilC1.html). The intrinsically unstructured nature of the flanking regions was determined using the protein prediction program DisEMBLE (Linding et al. 2003Go). Protein alignments of the 210 residue coiled coil region were performed with MUSCLE (Edgar 2004Go).

Plots of the pairwise distances against transition/transversion ratios in the nucleotides encoding the coiled coil regions showed no decline in transition/transversion ratio with pairwise distances other than for ant–ant and bee–bee comparisons within each fibroin, indicating that the third codon positions are saturated except between the species pairs Apis and Bombus, and Oecophylla and Myrmecia. Likewise, the first and second codon positions are saturated in transitions and in synonymous changes (unpublished data). Thus, phylogenetic information is restricted to the pattern of nonsynonymous mutations that was analyzed using equally weighted protein parsimony, and support for the tree was analyzed by doing 1,000 bootstrap runs using PAUP* version 4.0b10 (Swofford 2002Go).

Distances within heptad positions among the 4 taxa were calculated for the 4 fibroin groups using the Tajima-Nei nucleotide model (DNA distance) excluding the third codon position. The analysis assumed uniform rates of changes among sequences. These molecular evolutionary analyses were conducted using MEGA version 3.1 (Kumar et al. 2004Go).


    Results
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Coiled Coil Silks Are Found in Bees and Ants but Not in Aculeate Parasitic Wasps
We used Fourier transform infrared spectroscopy (FTIR) to identify the silk structure for species within each of the 3 superfamilies of the infrasuborder Aculeata (Hymenoptera) including bumblebees (B. terrestris; Apoidea), weaver ants (O. smargdina; Vespoidea), Australian bulldog ants (M. forficata; Vespoidea), and 2 species of aculeate parasitic wasps (Goniozus and Prosierola sp.; Chrysidoidea). A schematic depicting the relationship between these species is shown in figure 1. The silks of bumblebees, weaver ants, and bulldog ants all have FTIR spectra very similar to that of honeybee (A. mellifera; Apoidea) silk (fig. 2). We, therefore, expect these silks to have coiled coil structure, as does honeybee silk (Atkins 1967Go). Consistent with this, the spectral maxima in these 4 FTIR spectra (black lines, fig. 2) are at 1645–1646 cm–1, shifted ~10 cm–1 lower than a classical {alpha}-helical signal and broadened, as is typical of coiled coil proteins (Heimburg et al. 1999Go). In contrast, the silks of Goniozus and Prosierola wasps had spectra (gray lines, fig. 2) with maxima at ~1635 cm–1, suggesting predominantly ß-sheet structure.


Figure 1
View larger version (8K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 1.— Cladogram of the Hymenopteran infrasuborder Aculeata. The monophyly of the Aculeata and the radiation of the 3 superfamilies Chrysidoidea, Vespoidea, and Apoidea are well supported (Grimaldi and Engel 2005Go).

 

Figure 2
View larger version (17K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 2.— The Fourier transform infrared spectra of Aculeata silks. Spectra of honeybee, bumblebee, weaver ant, and bulldog ant silks are shown by black lines, and spectra of Goniozus and Prosierola wasp silks are shown by gray lines.

 
Bee and Ant Silks Are Composed of 4 Fibroins
Sequence data were obtained from 117, 131, and 23 random clones from cDNA libraries constructed from the silk glands of bumblebee, bulldog ant, and weaver ant larvae, respectively (table 1). As we obtained low levels of mRNA from the silk glands of the tiny weaver ants, the library we constructed from this species contained a high proportion of vector-only sequence and only 23 nonvector sequences were present in the approximately 500 clones examined. We identified the silk genes from the bumblebee, bulldog ant, and weaver ant by matching mass spectral data obtained from trypsin-digested silk to the predicted proteins encoded by the silk gland cDNAs. The Spectrum Mill software used (Agilent) generated scores for the quality of each match, and the sequence matches reported in this paper scored between 88 and 228 (table 1), considerably over the default cutoff for confident acceptance of valid matches (>20).


View this table:
[in this window]
[in a new window]

 
Table 1 Properties of Proteins Identified in the Silks of Honeybee, Bumblebee, Bulldog Ant, and Weaver Ant

 
We found that bumblebee silk comprised 5 proteins: 4 fibroins (BBF1–4) and a glue-like protein (BBSA1), homologous to those found in the honeybee (Sutherland et al. 2006Go). In the bulldog ant and weaver ant silks, we found 4 proteins (BAF1–4, WAF1–4; table 1) homologous to the fibroin proteins of the honeybee, but there were no homologs of the bee glue-like proteins. The measured amino acid compositions of the ant silks are consistent with the silk comprising only the 4 fibroin proteins in equimolar ratios (fig. 3). With the exception of the honeybee fibroins, the bumblebee and ant fibroins are unlike any other published sequences.


Figure 3
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 3.— Measured and predicted amino acid composition of the silk proteins of (A) bulldog ant cocoons and (B) weaver ant nests. The mean measured amino acid content is shown by filled bars, and predicted amino acid composition assuming equimolar amounts of each of the 4 silk fibroins is shown by open bars.

 
In all species, the silk proteins corresponded to the most abundant cDNAs in the silk gland libraries. Individual silk cDNAs accounted for 4–35% of the message in the libraries, and combined, the silk proteins accounted for 73% (bumblebee), 87% (weaver ant), and 99% (bulldog ant) of the silk gland cDNAs (table 1). We used quantitative PCR to show that in each species, including honeybees, the 4 silk genes were expressed in approximately equimolar amounts within the silk gland (table 1).

The Silk Fibroins Have Similar Architecture but Low Primary Sequence Identity
The architecture of the bee and ant fibroin proteins is very similar (fig. 4). The silk fibroins are all small (309–443 amino acids) and do not contain primary sequence repeats. The mature proteins each contain an approximately 210-amino acid coiled coil region (predicted by MARCOIL; Delorenzi and Speed 2002Go). The properties of the coiled coil protein sequences are described in detail below. The coiled coil regions are flanked by regions of variable length (N-terminal region, 33–105 amino acids; C-terminal region, 23–125 amino acids) that contain high levels of charged amino acids (average 20.1%). Generally, the flanking regions are longer in the ant fibroins (average 90 residues) than in the bee fibroins (average 50 residues) in part due to terminally located serine-rich regions of 15–23 residues that are intrinsically unstructured (predicted by DisEMBLE; Linding et al. 2003Go). Peptides from protease-digested mature silk matched to the N- and C-terminal diverged regions of all bee and ant fibroins, confirming that these regions are retained in the processed silk (unpublished data).


Figure 4
View larger version (5K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 4.— Common architecture of the 16 silk proteins. The proteins contain a signal peptide (shaded bar) and a coiled coil region containing around 30 heptad repeats flanked by regions of variable length (length range indicated by the dotted line).

 
The bee and ant silk genes have diverged extensively as we were unable to align the nucleotide sequences. We were also unsuccessful in aligning the protein sequences of the N- and C-flanking regions. Although sequence identity was detected between some regions in species in the same superfamily, primary sequence identity in these regions between ant and bee fibroins, or within a species, was unconvincing. We were able to align the 210-amino acid section of predicted coiled coil protein sequence from the 16 fibroins (fig. 5) using the protein alignment algorithm MUSCLE (Edgar 2004Go). As this alignment did not fragment the coiled coil heptad motifs predicted by MARCOIL (Delorenzi and Speed 2002Go) and agreed with manual alignments of the honeybee silk genes described in Sutherland et al. (2006)Go, it was used for all subsequent analysis. Pairwise protein identity was low, ranging from 21% to 68% identity between the ant species, 25% to 71% between the bee species, and 20% to 45% between fibroins from the different superfamilies.


Figure 5
View larger version (42K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 5.— Amino acid sequence alignment of the coiled coil regions of the silk proteins. HB, honeybee; BB, bumblebee; BA, bulldog ant; WA, weaver ant; F1–4, silk fibroins 1–4. Residues conserved between species within a fibroin set are indicted in bold typeface and residues conserved in all fibroins within a species have shaded boxes, excluding alanine conservation in any position, which is likely to arise through chance.

 
Unusual Composition of Coiled Coil Regions in Fibroins
Coiled coil protein sequences are characterized by a 7-residue [(abcdefg)n] repeating motif usually containing hydrophobic residues in the a and d positions and hydrophilic residues in the remaining positions (Woolfson 2005Go). The MARCOIL algorithm predicted that the bee and ant silk fibroins contain approximately 30 contiguous heptad repeats. The amino acid composition at the predicted heptad positions of the silks is very high in alanine in comparison to classic coiled coils especially in the a and d core positions (fig. 6; average composition—a: 65%, b: 35%, c: 37%, d: 51%, e: 39%, f: 33%, g: 30%). The core positions were also surprisingly abundant in small polar residues (serine and threonine) and low in large hydrophobic residues (leucine, isoleucine, and methionine) (fig. 6).


Figure 6
View larger version (19K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 6.— Summary of the properties of amino acids occupying each of the predicted heptad positions (a–g) in the coiled coil fibroins that comprise the silk of ants and bees. The graph shows the distribution of amino acids populating each heptad position for each fibroin (F1–F4). Beneath the graph is the average hydrophobicity (kJ/mol) and average side-chain length (atoms) of the residues occupying each position.

 
We calculated the average side-chain hydrophobicity of residues at each of the heptad positions for each fibroin (F1–F4) using the free energy of transfer between octanol and water for each amino acid relative to glycine (Eisenberg and McLachlan 1986Go). As expected, in all fibroins, the most hydrophobic positions were the classic core positions a and d (fig. 6). On average, the e and b positions were more hydrophobic than the remaining positions. We also compared the average side-chain length of residues at each of the fibroin heptad positions, calculated by counting the maximum number of linearly bonded atoms, excluding hydrogens. Average residue size was smallest in the a and d positions in all fibroin groups (fig. 6). Overall, residues in the e and b positions were also smaller than in the remaining positions.

Bees and Ants Have Retained Single Copies of 4 Ancestral Silk Genes
The pairwise distances between the nucleotides encoding the coiled coil regions of the 16 fibroin proteins were compared with the transition/transversion ratios in the same region. The result showed that the third codon positions of each fibroin gene were saturated except between the species pairs Apis and Bombus, and Oecophylla and Myrmecia. Likewise, the first and second codon positions were saturated in transitions and in synonymous changes (unpublished data). Phylogenetic information was therefore restricted to the pattern of nonsynonymous mutations in the coiled coil region and was analyzed using equally weighted protein parsimony (fig. 7).


Figure 7
View larger version (11K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 7.— Equally weighted protein parsimony analysis showing the relationship of the coiled coil regions of the honeybee (HB), bumblebee (BB), bulldog ant (BA), and weaver ant (WA) silk proteins (F1–F4). The area bound by the dotted line indicates divergence that occurred before the split of the ants and wasps (Vespoidea) from the bees (Apoidea) in the Late Jurassic (155 Myr; Grimaldi and Engel 2005Go).

 
Cladistic analysis of the coiled coil regions of the silk fibroins revealed that each species contained a single copy of each of the 4 fibroin proteins found in the honeybee (fig. 7). The presence of orthologs in each species suggests that after initial duplication events more than 155 MYA (Grimaldi and Engel 2005Go) in the common ancestor of the Apoidea and Vespoidea, there have been no further gene duplications or deletions in any of the species. At least since the split of the ants and bees, the regions of the fibroin genes encoding the coiled coil motif have been diverging in an approximately equivalent manner (fig. 7).


    Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
We have established that coiled coil silks are produced by members of the Apoidea and the Vespoidea, but not by insects of the Chrysidoidea, the oldest lineage of the Aculeata. Thus, this radically different building material replaced the plesiomorphic ß-sheet silk over 155 MYA in the period after the split of the Chrysidoidea from the social Hymenopteran lineages and before the split of the ants and bees.

In order to understand the important elements governing the assembly and function of coiled coil silks, we identified and sequenced the silk genes from bumblebees (Apoidae) and 2 species from Vespoidea. Along with the previously described honeybee silk sequences (Sutherland et al. 2006Go), this allows us to compare 2 representative sets of sequences that are separated by the greatest possible evolutionary distance. In addition, these silks have diverse applications: honeybees build silk/wax composite hives, bulldog ant larvae produce cocoons for protection during pupation, bumblebee larvae spin cocoons within waxy hives and then reuse the cocoons to store pollen and honey, and silk from weaver ant larvae is used to fasten fresh plant leaves together to form large communal nests.

We found that the silk from each species was made of 4 coiled coil proteins (fibroins), each an ortholog of a silk protein identified in honeybee silk. The bee silks also contained an additional protein but, as this protein could be washed out of the silk while the silk retained its fibrous nature, it is not considered to be an integral component of the silken fibers. The fibroin proteins contained extensive coiled coil regions of conserved length flanked by largely unstructured termini. The hydrophilicity of the termini is similar to the terminal regions of ß-sheet silks (Bini et al. 2004Go) and probably serves the same function of improving solubility in the gland.

The fibroin proteins had no primary sequence repeats and had diverged considerably between species. Only the protein sequence from the predicted coiled coil regions had sufficient identity to be aligned and even within this region pairwise similarity was low, at less than 50% between ant and bee fibroins. This degree of sequence divergence contrasts with the extreme conservation of primary sequence motifs seen between, and within, silk fibroins with ß-sheet structure, such as among some orb-weaving spider fibroins, which have been separated for 125 Myr of evolution (Gatesy et al. 2001Go). We attribute the sequence diversity within the coiled coil silks to the nature of their structure. For coiled coils, the character of the residues in particular heptad positions is important rather than the retention of identical sequence repeats. Hence, there is great latitude for amino acid substitution in the sequence while retaining the coiled coil structure.

The amino acid composition of the coiled coil regions of the silks is unusual, with high levels of alanine, a small hydrophobic amino acid, in the a and d core positions (fig. 6). Most known coiled coils contain predominantly large hydrophobic residues at these positions to maximize the hydrophobic forces stabilizing the core (Woolfson 2005Go). We ascribe the atypical composition of the coiled coils in the bee and ant silks to metabolic constraints. Silks must be produced in large quantities over a short period and are not recycled, so insects may be constrained to using nonessential amino acids for silk production. Of residues able to form alpha-helices, alanine is the most hydrophobic of the nonessential amino acids for bees (DeGroot 1953Go).

In physiological conditions, alanine is insufficiently hydrophobic to hold together a classic dimeric or trimeric coiled coil with a 2-residue (a and d positions) core (Liu and Lu 2002Go). X-ray diffraction data indicate that honeybee silk comprises tetrameric coiled coils (Atkins 1967Go), which are stabilized by a wider hydrophobic interface. Fibroin sequence analysis shows that residues at heptad positions a, d, and e are generally smaller and more hydrophobic than are those at the other positions (fig. 6). We propose that all 3 of these positions are included within the core of the coiled coil silks, as shown in figure 8. In fibroins 3 and 4, the b heptad position is also small and hydrophobic, which is consistent with its proposed partially submerged position in the core. Three-residue cores are expected to dictate a highly stable, antiparallel, tetrameric variant of the coiled coil structure (Deng et al. 2006Go).


Figure 8
View larger version (28K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 8.— A structural model for coiled coil silks. The {alpha}-helical strands corresponding to each of the fibroins are arranged in an antiparallel tetrameric configuration (direction indicated by arrows). Three residues (a, d, e) from each heptad repeat are buried in the core.

 
Other design elements of the fibroins are expected to contribute to the stability of the coiled coil silks. Alanine is a strongly {alpha}-helix–favoring residue (Pace and Scholtz 1998Go), so the overall high alanine content of the proteins would stabilize helices and facilitate coiled coil formation. The predicted coiled coil regions of the fibroins are long and uninterrupted, features known to increase stability (Su et al. 1994Go; Brown et al. 1996Go). In addition, the uniform small size of the residues in the core (fig. 6) allows tight packing of the core and strong van der Waals interactions.

The conservation of 4 orthologous fibroin genes in each species for 155 Myr suggests that each gene product performs a nonredundant function in the silk. In addition to the absolute number of fibroin genes being retained, there is strict conservation of approximately 10% of amino acid residues in the coiled coil region between orthologous proteins from the different species, despite a background of extreme sequence divergence (fig. 5). This conservation suggests that the function of each orthologue has been retained in each species. It is most likely that each protein is exclusively located in 1 of the 4 strands of the coiled coil, forming a hetero-tetrameric structure. Our quantitative PCR data indicate that the fibroin genes are expressed at equivalent levels in the silk gland suggesting that the gene products are required in equimolar amounts, as expected for a hetero-tetrameric structure. Within orthologous fibroin groups (i.e., F1 from each species), there is complete conservation of between 14 and 21 residues (fig. 5), not including alanine residues (which are highly abundant and likely to occur in the same position by chance). This is significant given the high levels of sequence diversity observed in the remaining protein sequence, and we propose that these conserved residues are involved in specific polar interactions with neighboring strands, as are often observed in coiled coil proteins (Woolfson 2005Go). The conserved residues may dictate either partner selection within the coiled coil or intercoil interactions that facilitate fiber synthesis.

It is unclear how coiled coil silk proteins could arise from ß-sheet structured silk proteins, which are large and repetitive with completely different functional constraints. Homology between the coiled coil fibroin genes suggests that they evolved from a single gene that emerged in the common ancestor of the Apoidea and Vespoidea, after the split from the Chrysidoidea (fig. 7). As coiled coils are a common and versatile motif, present in an estimated 10% of all eukaryote proteins (Liu and Rost 2001Go) and in many viral coat proteins, there are many templates that this original protein may have evolved from. The simplicity of the structure also opens the possibility that they evolved de novo. Why coiled coil silks arose in Apoidea and Vespoidea in place of ß-sheet silks is not known, but honeybee silk has different functional properties to ß-sheet silks that have been examined. Honeybee silk is tougher, primarily as a function of its greater extensibility, and unlike silkworm silk, it retains its properties when wet (Hepburn et al. 1979Go).

The social Aculeata are a spectacularly successful group; they dominate the terrestrial ecosystem in terms of animal biomass and energy consumption. Availability of an efficient building material is a prerequisite for the ability to build nests and protect juveniles, which is intrinsically linked to sociality (Wilson 1971Go). We propose that the development of a tough coiled coil silk has contributed to the huge evolutionary success of these insects.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
This work was financially supported by the Australian Grains Research and Development Corporation. We thank D. Anderson, A. Zwick, L. Malon, J. Todd, R. Newcomb, and J. LaSalle for assistance with insects and R. Wharton, Texas A&M University, and S. S. Nair, University of Calicut, for the kind loan of Chrysidoidea larval silk cocoons.


    Footnotes
 
Barbara Holland, Associate Editor


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 

    Atkins EDT. A four-strand coiled coil model for some insect fibrous proteins. J Mol Biol. (1967) 24:139–141.[CrossRef][Web of Science]

    Bendtsen JD, Nielsen H, von Heijne G, Brunak S. Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. (2004) 340:783–795.[CrossRef][Web of Science][Medline]

    Bini E, Knight DP, Kaplan DL. Mapping domain structures in silks from insects and spiders related to protein assembly. J Mol Biol. (2004) 335:27–40.[CrossRef][Web of Science][Medline]

    Brown JH, Cohen C, Parry DAD. Heptad breaks in {alpha}-helical coiled coils: stutters and stammers. Proteins (1996) 26:134–145.[CrossRef][Web of Science][Medline]

    DeGroot AP. Protein and amino acid requirements of the honeybee (Apis mellifera). Physiologia Comparata et D'Ecologia (1953) 3:197–285.

    Delorenzi M, Speed T. An HMM model for coiled coil domains and a comparison with PSSM-based predictions. Bioinformatics (2002) 18:617–625.[Abstract/Free Full Text]

    Deng Y, Liu J, Zheng Q, Eliezer D, Kallenbach NR, Lu M. Antiparallel four-stranded coiled coil specified by a 3-3-1 hydrophobic heptad repeat. Structure (2006) 14:247–255.[Medline]

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. (2004) 32:1792–1797.[Abstract/Free Full Text]

    Eisenberg D, McLachlan AD. Solvation energy in protein folding and binding. Nature (1986) 319:199–203.[CrossRef][Medline]

    Gatesy J, Hayashi C, Motriuk D, Woods J, Lewis R. Extreme diversity, conservation, and convergence of spider silk fibroin sequences. Science (2001) 291:2603–2605.[Abstract/Free Full Text]

    Grimaldi D, Engel MS. Evolution of the insects (2005) Cambridge: Cambridge University Press.

    Heimburg T, Schünemann J, Weber K, Geisler N. FTIR-spectroscopy of multistranded coiled coil proteins. Biochemistry (1999) 38:12727–12734.[CrossRef][Medline]

    Hepburn HR, Chandler HD, Davidoff MR. Extensometric properties of insect fibroins: the green lacewing cross-ß, honeybee {alpha}-helical and greater waxmoth parallel-ß conformations. Insect Biochem. (1979) 9:69–77.[CrossRef]

    Kumar S, Tamura K, Nei M. MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform. (2004) 5:150–163.[Abstract/Free Full Text]

    Lazaris A, Arcidiacono S, Huang Y, Zhou JF, Duguay F, Chretien N, Welsh EA, Soares JW, Karatzas CN. Spider silk fibers spun from soluble recombinant silk produced in mammalian cells. Science (2002) 295:472–476.[Abstract/Free Full Text]

    Linding R, Jensen LJ, Diella F, Bork P, Gibson TJ, Russell RB. Protein disorder prediction: implications for structural proteomics. Structure (2003) 11:1453–1459.[Medline]

    Liu J, Lu M. An alanine-zipper structure determined by long range intermolecular interactions. J Biol Chem. (2002) 277:48708–48713.[Abstract/Free Full Text]

    Liu J, Rost B. Comparing function and structure between entire proteomes. Protein Sci. (2001) 10:1970–1979.[CrossRef][Web of Science][Medline]

    Pace CN, Scholtz JM. A helix propensity scale based on experimental studies of peptides and proteins. Biophys J. (1998) 75:422–427.[Web of Science][Medline]

    Rudall KM. Silk and other cocoon proteins. In: Comparative biochemistry Vol. IV—Florkin M, Mason HS, eds. (1962) New York: Academic Press. 397–433.

    Scheller J, Gührs KH, Grosse F, Conrad U. Production of spider silk proteins in tobacco and potato. Nat Biotechnol. (2001) 19:573–577.[CrossRef][Web of Science][Medline]

    Su JY, Hodges RS, Kay CM. Effect of chain length on the formation and stability of synthetic {alpha}-helical coiled coils. Biochemistry (1994) 33:15501–15510.[CrossRef][Medline]

    Sutherland TD, Campbell PM, Weisman S, Trueman HE, Sriskantha A, Wanjura WJ, Haritos VS. A highly divergent gene cluster in honeybees encodes a novel silk family. Genome Res. (2006) 16:1414–1421.[Abstract/Free Full Text]

    Swofford DL. PAUP*. Phylogenetic analysis using parsimony (*and Other Methods) (2002) Sunderland(MA): Sinauer Associates.

    Vollrath F. Strength and structure of spider silks. Rev Mol Biotechnol. (2000) 74:67–83.[CrossRef]

    Vollrath F, Knight DP. Liquid crystalline spinning of spider silk. Nature (2001) 410:541–548.[CrossRef][Medline]

    Wilson EO. The insect societies (1971) Cambridge(MA): Harvard University Press.

    Woolfson DN. The design of coiled coil structures and assemblies. In: Fibrous proteins: coiled-coils, collagen and elastomers—Parry DAD, Squire JM, eds. (2005) San Diego (CA): Elsevier. 79–112.

Accepted for publication August 14, 2007.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrowOA All Versions of this Article:
24/11/2424    most recent
msm171v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Google Scholar
Right arrow Articles by Sutherland, T. D.
Right arrow Articles by Haritos, V. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sutherland, T. D.
Right arrow Articles by Haritos, V. S.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?