Controlling the Molecular Structure and Physical Properties of Artificial Honeybee Silk by Heating or by Immersion in Solvents

Honeybee larvae produce silken cocoons that provide mechanical stability to the hive. The silk proteins are small and non-repetitive and therefore can be produced at large scale by fermentation in E. coli. The recombinant proteins can be fabricated into a range of forms; however the resultant material is soluble in water and requires a post production stabilizing treatment. In this study, we describe the structural and mechanical properties of sponges fabricated from artificial honeybee silk proteins that have been stabilized in aqueous methanol baths or by dry heating. Aqueous methanol treatment induces formation of ß-sheets, with the amount of ß-sheet dictated by methanol concentration. Formation of ß-sheets renders sponges insoluble in water and generates a reversibly compressible material. Dry heat treatments at 190°C produce a water insoluble material, that is stiffer than the methanol treated equivalent but without significant secondary structural changes. Honeybee silk proteins are particularly high in Lys, Ser, Thr, Glu and Asp. The properties of the heat treated material are attributed to generation of lysinoalanine, amide (isopeptide) and/or ester covalent cross-links. The unique ability to stabilize material by controlling secondary structure rearrangement and covalent cross-linking allows us to design recombinant silk materials with a wide range of properties.


Introduction
Silk-based materials are of interest for many applications, with over 300 patents around silk protein materials approved in the past five years. The silk produced by honeybees shares many desirable properties with the better known silkworm or spider silks, but has a more varied amino acid composition giving it a different surface chemistry and a niche for different applications [1]. Natural honeybee silk is composed of four related fibrous proteins [1] which assemble into a coiled coil structure [2]. Honeybee silk proteins can be produced recombinantly in E. coli [3,4], and at high yield [4]. Solutions of single or all four recombinant honeybee silk proteins can be manufactured into material forms including fibers and films [4,5], electrospun mats [6] and sponges (this work). All forms require post fabrication treatment of some description to stabilize the material to render them stable in aqueous solutions.
Stabilization of honeybee silk protein has been achieved after inducing a partial molecular structural transition from coiled coil to b-sheet structure by immersion in 90% methanol solution [4], immersion in 70% methanol solution [5] or by water annealing [6], presumably by generating b-sheet cross-links between the coiled coil proteins. b-sheet structure in coiled coil silk materials has also been induced by dry heating to 215uC, well above the protein's glass transition temperature [5], or by autoclaving at 121uC [7]. Similar approaches have previously been used for stabilization of regenerated silkworm fibroin material [8,9,10,11,12].
Methanol induces formation of b-sheets by increasing the hydrophobicity of the solvent and therefore weakening internal protein-protein hydrophobic interactions. At the same time it decreases the availability of water for hydrogen bonding thereby driving protein-protein hydrogen bonding [13]. In mixed wateralcohol systems the underlying mechanisms are much more complex with several distinct physicochemical effects impacting such transitions [14]. The structural transition of silkworm fibroin to b-sheet structure is rapid, and more complete in aqueous methanol solutions than in dry methanol solutions [12]. Methanol treatment of reconstituted silkworm silk increases the compressive strength of silkworm fibroin sponges, but the rapid dehydration induced also makes them brittle [11]. Changes to the mechanical properties of recombinant honeybee silk materials due to methanol treatment have not been previously investigated.
Collagen has been used extensively as a biomedical material. Various chemical and physical treatments have been used to stabilize collagen, giving insolubility and durability. For example, stabilization of collagen can be achieved by dry heat treatments that involve heating the material under vacuum to temperatures of 100-120uC for several days [15]. The stability of the treated material is due to the covalent cross-linking of the proteins rather than secondary structure protein changes. The major covalent cross-links found are amides formed by condensation reactions between amine and carboxylic acid side groups [15]. A small amount of lysinoalanine formed by dehydration of serine (or desulfation of cysteine) followed by reaction with lysine, is also present [16]. Heat treatment significantly increases the strength, extensibility and elasticity of regenerated collagen fibers [17].
This work systematically studies the effects of different stabilization treatments of artificial honeybee silk on the material's protein structure and mechanical properties. The two treatment types studied are submersion in aqueous methanol of varying concentrations and dry heat treatment using varying temperatures and times. The honeybee material form chosen for study is freezedried sponges, as these have a high surface-to-volume ratio that allows good penetration of solvents and can be manufactured at physical sizes large enough to allow simple and robust mechanical property measurements.

Generation of protein sponges
A honeybee silk protein (AmelF3; NCBI accession no: NP_001129680) was expressed into Escherichia coli (E. Coli) inclusion bodies and purified as previously described [4]. Briefly, purified inclusion bodies were solubilised in 3% aqueous sodium dodecyl sulfate (SDS) solution then treated with 300 mM KCl. The potassium serves to precipitate the dodecyl sulfate of SDS which can then be removed by centrifugation, thereby reducing SDS levels to around 0.3% (weight/volume) and producing highpurity AmelF3 solution [4]. The AmelF3 solution was dialysed against 15% polyethylene glycol (PEG), 0.25% SDS to remove salt and generate a solution containing 3.6% protein, 0.3% SDS. The solution also contained 0.1% salt comprising residual Na + derived from the SDS after precipitation of the dodecyl sulfate and K + with Cl 2 counter ions. SDS concentration was determined according to Rusconi et al. [18], salt concentration was determined by calculation of dilution during dialysis, and total solute weight was determined by weighing aliquots of dried solution. The solution was poured into silicone rubber moulds (146566 mm; RL060, ProSciTech, QLD), frozen at 220uC overnight, and placed in a freeze-dryer (FD355DMP, FTS Systems) for 24 hours to generate sponges with typical dimensions of 12.664.565.4 mm. The silk sponges were stored in sealed plastic bags at room temperature until use.
For comparison, stabilized collagen sponges were made from ovine skin collagen. Minced ovine skin was digested with 1 mg/ mL pepsin at 4uC for 24 hours in 100 mM acetic acid adjusted to pH 2.5 with HCl. The hydrolysate was purified using two 0.7 M NaCl precipitation steps [19]. SDS-PAGE showed that the collagen was predominantly type I, but a small quantity of type III collagen, ,5%, was also present. PEG 4000 was used to precipitate collagen as fibre-like aggregates [20] as previously described [21]. Purified collagen was dissolved at 2 mg/mL in 20 mM acetic acid then adjusted to pH 7.0 with 40 mM sodium phosphate buffer. PEG (40% w/v) was added slowly with stirring to a final concentration of 4 dry wt% solution and the solution was left at 4uC for 16 hours. The precipitate was collected as a loose pellet by centrifugation at 10006g for 15 min then re-suspended in 20 mM sodium phosphate buffer, pH 7.0, and the precipitate re-collected as above after 1 hour. The precipitate was re-suspended in 10 mM sodium phosphate buffer, pH 7.0, and the precipitate re-collected after 1 hour at 10006g for 30 minutes. The collagen paste was adjusted to 4% protein, transferred into silicone rubber moulds (RL060, ProSciTech, QLD) and freezedried. Collagen sponges were suspended in dry ethanol and air in the pores was removed by vacuum prior to cross-linking with a 10:1 w/w excess of 1-ethyl-3-[3-dimethylaminopropyl]carbodiimide hydrochloride for 16 h at room temperature. The sponges were then washed extensively in ethanol and air dried.

Stabilization of honeybee silk sponges and water solubility testing
Honeybee silk sponges were treated to decrease water sensitivity by either heating in air to temperatures of 160-190uC for times ranging between 10-120 minutes, or by immersion for 48 hours in aqueous methanol solutions, with methanol concentrations ranging from 0-100%. In order to test water stability, sponges were cut into 4 pieces and each sample, under ambient laboratory conditions, was weighed on a microbalance prior to heat or methanol treatment (initial mass). After methanol treatment, sponge samples were transferred through a dilution series of 40% and 20% methanol into water. All samples were incubated individually in 1 mL water at room temperature on an orbital shaker for time periods ranging from 1 to 32 days, dried at 60uC for 2 hours and then weighed to determine mass loss. Mass loss is reported as the percentage of the initial mass. Initial mass loss and rate of mass loss were determined from the y-intercept and slope of the linear regression through the mass loss vs. time data points, respectively. Four samples per treatment were assessed and data collected at 0, 1, 2, 4, 8, 16 and 32 days. SDS levels in the methanol solutions or water after sponges were removed were calculated according to Rusconi et al. [18], and protein concentration was determined using the QuantiPro BCA assay kit (Sigma; St Louis, MO).

Compression testing
The mechanical properties of honeybee silk sponges and stabilized collagen sponges were measured using compression tests on an Instron 5500R (Instron, USA) fitted with a 2.5 N static load cell. Methanol treated samples were treated in 60-100% methanol in water for 48 hours and then transferred through a dilution series of 80, 60, 40 and 20% methanol solutions each for 10 minutes and then soaked for 24 hours in phosphate buffered saline (PBS) prior to testing. Heat treated samples and collagen sponges were soaked for 5 minutes in PBS with gentle squeezing to ensure saturation prior to testing. Sponge samples were compressed to 190-200 mN at 2 mm/min then decompressed at 2 mm/min. The compression tests were repeated three times on each sample with 1 minute relaxation time between tests. The cross-sectional area of the samples, measured on a light microscope to be around 57 mm 2 , was used to convert force values to stress values. Sample thickness was determined during the compression test, as the distance between the base plate and the point at which the compression plate contacts the sample (F = 1 mN).

Raman spectroscopy
Raman spectra were obtained from dried honeybee silk sponges at a resolution of 4 cm 21 using a Bruker RFS-100 FT-Raman spectrometer (Karlsruhe, Germany) equipped with an Adlas Nd:YAG laser operating at 1.064 mm and 500 mW and with a liquid nitrogen cooled Germanium diode detector. Spectra were collected using 180u backscatter geometry with the samples held in a compression cell described elsewhere [22]. Data acquisition over 512 scans was performed using Bruker OPUS software (version 3.1). Four spectra were obtained from different areas of each sample and averaged to produce a final spectrum for analysis. All spectra obtained from a given sample were found to be in excellent agreement. Data manipulation was carried out using Grams AI v8.0 software (Thermo Electron Corp., USA). All spectra were normalized on the C-H bending mode at 1449 cm 21 as this mode is insensitive to protein conformation. The structural changes of the honeybee silk sponges after the different treatments was also assessed by the quantitation of the overlap area between the second derivatives of the amide I bands [23]. Briefly, the method, which was implemented using Matlab R2010a (Math Works, USA), involves the normalization of the spectra between 1620 and 1730 cm 21 followed by the calculation of their 2 nd derivatives using the method of Savitzky and Golay [24]. The two 2 nd derivative spectra to be compared were passed through a logical filter which takes the lowest intensity value of the pair and creates an array that constitutes a spectrum representing the overlap area.

Amino acid analysis
Quantitative amino acid analysis via acid hydrolysis was carried out on untreated and heat treated honeybee silk sponge samples. Digestions were carried out in triplicate, analyzed in duplicate and the results reported as an average. Deficits of amino acids in the treated samples relative to untreated samples were taken to indicate residues that had been modified by heat-treatment.
Modifications such as dehydration reactions producing isopeptide bonds would be reversed by acid hydrolysis. In order to identify residues involved in such modifications, amino acid analyses were also performed following enzymatic digestions of silk sponge samples. Enzymatic digestions were performed in duplicate on samples of heat-treated and non-heat-treated honeybee silk sponges that were first washed with 70% methanol. Methanol would not be expected to induce any covalent modification but removes detergent that might interfere with enzymatic digestion or mass spectral analysis (below).
Approximately 3 mg of sample was accurately weighed and resuspended in 100 mM HEPES buffer pH 7.5. The samples underwent a 24 hour enzymatic hydrolysis at 37uC with Pronase E, Leucine Aminopeptidase M, and Prolidase in a final volume of 500 mL. Samples were heated at 100uC for 10 minutes before undergoing a further 24 hour enzymatic hydrolysis at 37uC with carboxypeptidase.
Levels of amino acids after both treatments were analyzed using the Waters AccQTag Ultra chemistry on a Waters ACQUITY UPLC system. The non-standard amino acid e-(c-glutamyl)-lysine was quantified by reference to a standard (Sigma) and eluted with a retention time similar to cysteine, an amino acid that is conveniently absent from AmelF3.

Liquid chromatography-tandem mass spectrometry
Targeted liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to test for the presence of e-(c-glutamyl)lysine, a dipeptide that would remain following enzymatic digestion of honeybee silk sponges containing Glu-Lys iso-peptide crosslinks. The method is based on a published method [25], showing tandem mass spectra that distinguish the isolated isopeptide crosslink from isomeric dipeptides of lysine and glutamate. Heat-treated and non-heat-treated material was digested as above, acidified with formic acid (1% v/v final) and 1 ml was separated on an Agilent 1200 capillary LC system coupled by electrospray ionization to an Agilent XCT ion trap mass spectrometer. The column (Agilent Zorbax 300SB-C18, 3.5 mm, 5060.3 mm) was eluted isocratically with 3% acetonitrile/0.1% formic acid at 5 mL/min for 5 minutes before washing with 90% acetonitrile/0.1% formic acid and re-equilibration with 3% acetonitrile/0.1% formic acid. The ion trap was tuned to retain only ions of 276 m/z, corresponding to the MH + ion from e-(c-glutamyl)-lysine. Fragmentation energy was adjusted to produce two dominant ions at 131 and 148 m/z from a standard solution of e-(c-glutamyl)-lysine (Sigma) to obtain spectra consistent with that previously reported [25]. The standard eluted at 1-2 minutes. Mass spectra from the same retention time were averaged for digested honeybee silk sponges for comparison with the spectrum from the iso-peptide standard.

Efficiency of stabilization treatments
The efficiency of stabilization treatments on recombinant honeybee silk materials could be assessed by the performance of the material in water, with untreated sponges swelling and rapidly dissolving. Sponges that were treated either by heating dry material for various times or by immersion in aqueous methanol solutions were incubated in water for up to 32 days and their water solubility monitored (Table 1). Treatment with aqueous solutions containing .50% methanol but ,100% methanol, or heat treatment at 190uC, imparted a high degree of water stability to the sponges.
When a sponge was subjected to its initial wet exposure, either as part of its treatment and/or during the first day of water stability testing, a mass loss of 18-23% occurred (Table 1). This loss was primarily due to leaching of SDS from the sponges. The detergent was detected in the methanol solutions or water at levels equivalent to what was originally in the sponge (data not shown). The remaining mass lost is likely to be made up of leached salts and low molecular weight peptides present as contaminants or degradation products, and bound water loss during the dry heat treatment process (4-6%) or resulting from molecular rearrangement.
Samples treated with aqueous methanol solutions containing ,50% methanol dissolved in the treatment solutions. Sponges stabilized for up to 48 hours in solutions containing very low levels of water (100% methanol treatment), swelled extensively when transferred to water and continually lost protein at a rate of around 1% per day ( Table 1). The swollen sponges were too weak to be lifted out of the solution intact and hence could not undergo mechanical testing. The water stability of the swollen sponges, as determined after the sponge fragments were collected by centrifugation and weighing, was significantly greater than the untreated material (Table 1) despite their mechanical weakness.
Samples that had been heat treated at temperatures ,190uC (160uC, 180uC) were significantly less water stable than those treated for equivalent times at 190uC, the samples treated for 10 minutes dissolving completely in the first 24 hours of incubation (Table 1). Samples heated for 10 minutes at 190uC continually lost protein, at a rate of around 0.3% per day. By contrast the samples heated for longer time periods at 190uC were effectively impervious to water after their initial mass loss ( Table 1).

Effect of aqueous methanol stabilization on protein molecular structure
Raman spectroscopy is a well established technique for assessing changes in protein conformation [26,27,28], with the semiquantitative estimation of protein conformation by Raman spectroscopy in good agreement with those determined by X-ray diffraction and infrared spectroscopy [28]. Raman spectra of aqueous methanol treated sponges obtained immediately after treatment are shown in Figure 1. Spectra obtained from sponges that had been soaked in water for 24 hours after aqueous methanol then dried and analyzed were effectively identical to the spectra collected from dry sponges immediately after treatment, indicating that water immersion did not affect the protein secondary structure after stabilization (data not shown). Examination of the amide I region of the spectra indicated a distinct change in protein conformation as a consequence of the aqueous methanol treatment (Fig. 1). There is an increase in the peak at 1666 cm 21 , attributed to ß-sheet, and a corresponding decrease in the peak at 1655 cm 21 , attributed to coiled coil structure [5,27,29]. As it is not possible to resolve ß-sheet and disordered protein structure in the amide I region [29] there is a possibility that both are actually present. In the amide III region (not shown) where the ß-sheet and disordered proteins are better resolved, there is a slight increase in a feature at 1258 cm 21 confirming that a small level of disordered protein is being formed [29], particularly at the lower methanol concentrations. There is a trend towards decreasing structural change with increasing methanol concentration (Fig. 1, inset), so that the sample treated with 60% methanol has the highest ß-sheet content (corresponding with the lowest level of water swelling in Table 1), and the sample treated with 100% methanol has the lowest ß-sheet content (corresponding with highest level water swelling and some water solubility in Table 1).
The structural change upon aqueous methanol treatment was further investigated using the band overlap approach developed by Kendrick et al. [23] for the analysis of amide I bands ( Table 2, Fig. 2). The band overlap analysis provides a single measure of difference and is ideal for assessing minor conformational changes in a series of samples, and thus was used in preference to the more traditional full band deconvolution method [26,27,28]. Band overlap analysis confirmed that the greatest conformational change (26%) was observed for the 60% methanol treated sponge while the least (11%) was found for the 100% methanol treated sponge. The percentage conformational change was found to be a linear function of methanol concentration with an r 2 of 0.9838. The changes consisted of a major loss in coiled coil structure (1653 cm 21 ) off-set by a major gain in ß-sheet structure (1669 cm 21 ). There was also a minor loss in ß-turn structure observed at 1688 cm 21 [27].
In aqueous methanol treatments, water acts as a plasticizer, lowering the glass transition temperature of the protein and thus increasing the mobility of the protein and allowing structural rearrangement to occur more readily. Since water is also a solvent for the protein, stabilization of the material requires a compromise to be reached; too much water and the dissolution process dominates, too little or no water and the rate of structural rearrangement becomes unacceptably slow. Thus sponges treated in less than 50% methanol dissolved in the high aqueous environment. Sponges treated in 50-90% methanol were in an environment with sufficient water to allow protein structural rearrangement but also sufficient methanol to keep the protein precipitated and in the solid form. In the absence of water (100% methanol treatment), the sponges did not undergo sufficient protein rearrangement during the treatment period (48 hours) to generate a stable material and the material swelled extensively when transferred to water ( Table 1).
Effect of dry heat stabilization on protein molecular structure Previously, it had been shown that dry heating of artificial honeybee silk protein films to 215uC led to substantial ß-sheet formation [5]. In this study, the Raman spectra of sponges heated to 190uC remained similar to that of untreated sponges with only a slight structural transition in the direction of ß-sheet indicated in the amide I region (Fig. 3) and no significant differences observed in the amide III region. In the time series of samples treated at 190uC, the effect is smallest in the sample treated for 10 minutes then approximately constant for treatments of 30 minutes or longer (Fig. 3, insert). In all heat treated sponges the magnitude of the structural transition is less than that seen in samples treated with 100% methanol. However, the heat treated samples have high water stability, significantly better than sponges treated with 100% methanol ( Table 1), suggesting that the heat treated samples are undergoing stabilizing changes other than the dramatic protein structural transitions observed in methanol treated samples.
Amide I band overlap analysis results are presented in Table 2 and Figure 4. The conformational similarity of the sponges treated at 160 and 180uC for 60 minutes and 190uC for 10 minutes showed all underwent structural rearrangements of the order of 6% upon treatment. These changes can be associated with the loss of coiled coil structure coupled with gains in ß-sheet structure. Sponges treated at 190uC for more than 30 minutes also have very Table 1. Stability of honeybee silk sponges in water over 32 days. Standard deviations are given in brackets, ns = not significant at the 95% confidence limit.

Mechanical properties of treated sponges
Honeybee silk sponges treated by heating or by immersion in aqueous methanol solutions were subjected to compression tests in phosphate buffered saline (PBS) solution to assess their wet mechanical properties. For comparison, sponges of collagen, an established biomedical material, chemically cross-linked to a level that allows the sponge to return to its original size, were also tested. Representative compression curves are shown in Figure 5A and summarized results are shown in Table 3. Since the sponges have been compressed to a constant force, the compressive strain at this load is inversely related to the stiffness of the samples. The best measure of stiffness, however, is the stress at a given strain. Clearly, the honeybee silk sponges are softer than the collagen sponges requiring less stress to compress the sponge to 25% of its initial thickness (i.e. at 25% strain, Table 3). There is no significant difference between the honeybee samples treated in 60% or 80% methanol, both of which are softer than the samples heated at 190uC for greater than 10 minutes (compare stress at 25% and 50% strain). The energy required to compress the sponges is influenced by both the stiffness and shape of the curve and so is best used only as an indicator of the recovery of the samples between successive compressions (Table 3 and Fig. 5D). The samples that were heat treated show an initial J-shaped compression curve after 10 minutes heating (Fig. 5B). At longer heat treatment times a small yield point is evident (Fig. 5B) and the sample stiffness (stress at 50% strain) and energy required for compression increase logarithmically with time of treatment ( Fig. 5C and Table 3).
Resilience is measured as the ratio of the energy required to deform a sample, to the energy recovered when the load is removed. In other words, it is a measure of a material's ability to recover from deformation. Resilience values of the stabilized honeybee silks vary between 40 and 50%, substantially lower than the 80% shown by collagen. It should be noted, however, that the collagen sample, being much stiffer, was only compressed to about 40% compared to 70-90% for the honeybee silk samples. As a general rule, resilience decreases with increasing energy of deformation [30], therefore comparing samples which require vastly different energies of deformation can be misleading. It is important to note too, that the values in Table 3 are a measure of the instant recovery during the cyclic test and that given more time viscoelastic materials will recover further. Thus, although the   honeybee silk sponges only recover about 50% of the deformation energy during the cyclic test, if they are allowed one minute relaxation time then they show only a 10% decrease in energy to compress (Fig. 5D). If a sponge is left in solution for 24 hours after a series of compressions its mechanical properties recover fully to those seen in the initial compression. Table 3 includes data for the seventh compression of a 120 minute heat treated sample compressed seven times over a period of 32 days in PBS. Whilst the sample has maintained its integrity and the energy required to compress it remains unchanged, its properties have changed. The sponge has become softer and more resilient after multiple compressions and a lengthy immersion in PBS.

Mechanism of stabilization in heat treated materials
As described above, the secondary structural changes seen in the heat treated materials are insufficient to explain their physical and mechanical properties. The observed properties strongly suggest covalent cross-linking as the stabilization mechanism. Amino acid analysis was performed on acid hydrolyzed control and heat treated (60 minutes at 190uC) sponges and results are shown in Figure 6. Significantly, heat treatment resulted in losses in lysine and serine. Lysine and serine residues are known to form lysinoalanine crosslinks in collagen [16] and other proteins [31] subjected to heat treatment, through dehydration of the serine to form a dehydroa-lanine electrophile which can then undergo nucleophilic attack by the lysine's primary amine [31]. The analogous reaction with threonine, which also exhibits a significant reduction in heat treated sponges, is known to occur, forming methyl-lysinoalanine. Both lysinoalanine and methyl lysinoalanine are difficult to detect by Raman spectroscopy as the strongest bands, those associated with the symmetric C-N-C stretching vibrations, are expected in the 850 to 900 cm 21 region [32] which is already rich [29] in protein skeletal vibrations. No significant spectral changes were observed in this region as a function of thermal treatment.
Heat treatment will also generate isopeptide (amide) [33] or ester bonds between acidic residues such as Asp and Glu in close proximity to either amine-functional (Lys, Arg, His, Asn, Gln), or hydroxyl-functional (Ser, Thr, Tyr) residues, respectively [34]. Isopeptide and ester bonds cannot be inferred following acid hydrolysis as this treatment will hydrolyse all amide and ester bonds, returning the original amine, hydroxyl and carboxylate functions, except in the case where Arg could be thermally converted to a primary amine and form isopeptide bonds. In this event, hydrolysis would, of course, only return the primary amine. Indeed, there does appear to be a slight reduction in Arg content upon heat treatment at 190uC for 60 minutes which may be indicative of this manner of cross-linking having occurred. Isopeptide bonds are difficult to observe directly in the Raman spectra owing to their relatively low abundance compared to the backbone amides which will absorb in the same region [29]. While the carbonyl stretching vibrations of ester groups are expected on the high frequency side of the amide I band, its weak intensity in the Raman compared to that in the infrared [32] makes it hard to detect at low concentrations.
Therefore, to measure the contribution of isopeptide and ester bonds to the heat treated material's properties, amino acid analysis was conducted on enzymatically hydrolysed control and heat treated (60 minutes at 190uC) sponges (Fig. 7). The enzymes used to hydrolyse the protein's peptide bonds will not cleave isopeptide bonds and comparative reduction in amino acid levels will indicate that the amino acid is involved in cross-linking reactions. The enzymatically hydrolysed heat treated material showed significant losses in lysine, asparagine, aspartic acid and glutamic acid in comparison to the control sponge. The concentrations of histidine and tyrosine were not found to decrease. These differences as well as that observed for serine are represented as mole percentages in Figure 8. These results are consistent with the formation of isopeptide linkages between lysine/asparagine and aspartic acid/ glutamic acid as well as the formation of lysinoalanine linkages between lysine and serine and ester linkages between aspartic acid/glutamic acid and serine Enzymatic hydrolysis of material containing Lys-Glu isopeptide bonds will liberate e-(c-glutamyl)-lysine. Targeted LC-MS/MS directly demonstrated the presence of e-(c-glutamyl)-lysine in enzymatic digests of heat-treated honeybee silk sponges but not in non-heat-treated sponges (Fig. 9). Following calibration with a e-(c-glutamyl)-lysine standard, the amino acid analysis traces were re-examined. A peak corresponding to e-(c-glutamyl)-lysine was present in the heat-treated sponge at about 4.4 mg/g sponge corresponding to approximately one e-(c-glutamyl)-lysine moiety per protein molecule or two Lys-Glu isopeptide links on average between protein molecules. Heat treatment is also expected to result in analogous isopeptide bonds between lysine and aspartate residues.

Conclusions
Recombinant honeybee silk proteins can be fabricated into a range of forms that could be used as biomedical materials. The performance of the material can be tailored by stabilization treatments. Aqueous methanol treatment renders the material water insoluble and generates a reversibly compressible sponge. The methanol treatment generates ß-sheets in the material, presumably leading to ß-sheet cross-linking between proteins. The amount of ß-sheet formation can be controlled by varying  (2) All samples were compressed to a maximum load of 0.2 N. All values are for the first cycle except the sample marked by the asterisk (*) which is for the first cycle of the 3 rd compression experiment after 32 days in PBS. Values in brackets are standard deviations. doi:10.1371/journal.pone.0052308.t003 Figure 6. Amino acid analysis results obtained by acid hydrolysis, comparing a control untreated honeybee silk sponge sample to a sponge sample heat treated at 1906C for 60 minutes. Error bars represent one standard deviation. doi:10.1371/journal.pone.0052308.g006 methanol concentration. Dry heat treatments also render the material water stable with heated samples being stiffer in aqueous solution than methanol treated counterparts. In contrast to the methanol treatment, the changes in physical properties do not correlate with secondary structure molecular changes, with Raman spectra remaining broadly the same. Amino acid analysis on acid digested samples revealed a substantial loss of lysine and marginal losses of serine and threonine in dry heated sponges indicating the possible formation of lysinoalanine and/or methyllysinoalanine cross-links, respectively. Amino acid analysis on enzyme digested samples showed an additional loss of asparagine, aspartic acid and glutamic acid, consistent with the formation of isopeptide and ester cross-links as a consequence of heat treatment. The presence of the isopeptide e-(c-glutamyl)-lysine was confirmed directly by LC-MS/MS. The ability to readily control the formation of two distinct crosslinking mechanisms that occur intrinsic to the material is unusual. Whilst silkworm or tussah silk fibroins readily form ß-sheet crosslinks, heat treatment is unlikely to result in significant cross-linking because these proteins are poor in the charged residues that can participate in dehydration reactions. Intrinsic covalent cross-links occur in collagen upon heat treatment. However, structural rearrangements are constrained due to the existing native crosslinks that stabilize its coiled coil structure. Recombinant honeybee silk proteins are rich in charged residues and contain regions that readily form ß-sheet, so stabilization treatments can cause covalent cross-linking, structural transitions or a combination of the two. In  contrast to collagen-based biomedical materials, these honeybee silk proteins can be effectively stabilized without the need for any chemical crosslinking treatments that may leave cytotoxic chemical residues in the sample.