Laura
Senior
ab,
Matthew P.
Crump
c,
Christopher
Williams
c,
Paula J.
Booth‡
a,
Stephen
Mann
c,
Adam W.
Perriman
d and
Paul
Curnow
*a
aSchool of Biochemistry, University of Bristol, Medical Sciences Building, University Walk, BS8 1TD, UK. E-mail: p.curnow@bristol.ac.uk; Fax: +44(0) 17 331 2168; Tel: +44(0) 117 331 1837
bBristol Centre for Functional Nanomaterials, University of Bristol, BS8 1TD, UK
cSchool of Chemistry, University of Bristol, BS8 1TD, UK
dSchool of Cellular and Molecular Medicine, University of Bristol, BS8 1TD, UK
First published on 12th February 2015
The 19-mer synthetic peptide known as R5 has been used widely in studies of peptide-driven silica condensation. Despite this, the structure and function of R5 have not yet been fully characterized. Here, we present a systematic study of R5 silicification focusing on three key variables: the concentration of the peptide, the concentration of the silica precursor silicic acid, and the solution pH. Additionally, we present the first study of R5 secondary structure in the presence and absence of silicic acid and introduce one-dimensional and two-dimensional solution NMR to probe both structure and higher-order peptide aggregation. We find that R5-directed silicification is linear with regard to silicic acid and H+ but, unexpectedly, that silicification appears to be cooperative with respect to peptide concentration. We also find that R5 is a random coil ensemble at subsaturating silicic acid concentrations and does not spontaneously self-assemble to form discrete aggregates in solution. These data contradict a model that invokes the functional micellization of R5 and provide a framework for future studies with the R5 peptide.
Diatoms (Bacillariophyceae) are a group of abundant and diverse algae that are sheathed in a silicified external cell wall known as the frustule. The biomolecules responsible for diatom silicification8,9 include silicon-responsive12,13 silicifying polypeptides known as silaffins that were originally isolated as proteinaceous components of the frustule14 and have now been characterized in some detail.15–19 In Cylindrotheca fusiformis, a group of homologous silaffin isoforms known collectively as silaffin-1A1 are produced in vivo by endoproteolysis of the gene product Sil1p.12–14 Sequence analysis of silaffin-1A1 led to the identification of the repeat sequence motif SSKKSGSYSGSKGSKRR(I/N)L. Mature native silaffin-1A1, from which the C-terminal RRIL group is apparently removed, has unusual posttranslational modifications so that 5–11 repeats of N-methylpropylamine are introduced at the ε-amino group of K3 and K15, K4 is ε-N,N-dimethyllysine, K12 is a phosphorylated ε-N,N,N-trimethyl-δ-hydroxylysine, and all of the serine residues are phosphorylated. Both native silaffin-1A1 and a 19-mer synthetic peptide based upon the entire repeat motif, termed R5, were shown to induce rapid silica precipitation from silicic acid solutions in phosphate/citrate buffer. Although silaffin-1A1 could induce silicification down to pH 5 and was even partly active at pH 4, R5 was only active above pH 6.14 There is evidence that both silaffins and R5 become entrapped within the precipitated silica, giving rise to peptide–silica precipitates.14,20
There have since been numerous attempts to exploit this peptide sequence for biotechnology. For example, R5-associated silicification has been shown to be a generic strategy for bulk enzyme encapsulation,21 and covalent conjugation to R5 allows the encapsulation of small molecules22 and peptides.23 Additionally, R5 has been introduced to enzymes and other proteins as a silicifying motif by chemical and transcriptional fusion to enable the encapsulation of specific target proteins24–26 and for the formation of novel hybrid materials27,28 and biosensors.29,30 Modifications to the peptide sequence can also expand the utility of R5; for example to induce silicification at biomembranes.31 The dependence of silicification upon the R5 peptide sequence has been investigated and the C-terminal sequence RRIL appears to be essential for mineralization,20,32 with suggestions that this motif facilitates functionally-important, spontaneous peptide self-assembly either via micellization20,33 or through salt-bridging.32 These models are consistent with the electrostatically-driven aggregation that appears to be necessary for silicification by the zwitterionic native silaffins.16 Interestingly, scrambling the RRIL motif or moving it to the N-terminal or centre of R5 have a marginal effect on peptide activity but induce substantial changes to the resultant silica morphology.32
Significantly, the solution properties of R5 have not yet been fully characterized. To address this, we present here a study of the structure and function of R5 in physiologically-relevant buffers with a physiologically-relevant substrate, silicic acid (Si(OH)4). In particular, we examine whether R5 adopts any secondary structure in solution and study silicification as a function of peptide concentration, silicic acid concentration, and pH. We also determine whether the peptide spontaneously self-assembles in the absence and presence of silicic acid.
NMR spectra were recorded on a 600 MHz Varian VNMRS spectrometer equipped with a triple resonance cryogenically cooled probe-head. NMR samples contained 0.5 mM R5 in 50 mM sodium phosphate buffer and 90% H2O/10% D2O at pH 7.0 or 5.5. High-resolution one-dimensional proton spectra and two-dimensional nuclear Overhauser enhancement spectroscopy (NOESY), double-quantum filtered correlation spectroscopy (DQF-COSY) and total correlation spectroscopy (TOCSY) spectra were recorded. Mixing times were 150 ms and 65 ms for NOESY and TOCSY experiments respectively. Two-dimensional data were acquired with a spectral width of 8012 Hz, 4096 complex points in f2 and 800 complex points in f1 with sculpted excitation water suppression for both the NOESY and TOSCY spectra.35 Spectra were processed and plotted using VnmrJ 4.0 (Agilent).
For scanning electron microscopy, vacuum-dried R5–silica and protamine–silica precipitates were fixed onto aluminium stubs using adhesive tabs and were sputter-coated with silver. All samples were examined using a JSM 5600LV SEM (JEOL) fitted with the ISIS 300 X-ray analysis system (Oxford Instruments).
As anticipated, the presence of either R5 or protamine resulted in rapid silica precipitation (within seconds) once silicic acid was raised above the saturating concentration of 2 mM at pH 7 (Fig. 1). No such rapid precipitation was observed in the absence of either peptide, confirming the roles of both R5 and protamine in promoting silica formation. The dry weight of these peptide–silica precipitates was found to be linear with silicic acid concentration, with the gradients of these linear curves being 0.36 ± 0.02 and 0.44 ± 0.02 mg product per ml per mM silicic acid for R5 and protamine, respectively (R2 = 0.96). Thus, under equivalent conditions, protamine generates approximately 1.2-fold greater dry weight of precipitate versus R5. This modest difference in activity presumably arises from the higher charge density carried by protamine (21 charged amino acids of 32 residues vs. 6 of 19 residues in R5). It was previously shown that the number of charged residues within the R5 sequence was a critical factor in controlling the rate and yield of silicification, and that the specific position of lysine and arginine groups within the peptide sequence also influenced silica morphology.20,32 The difference in activity between R5 and protamine may also reflect the side chain chemistries of the cationic amino acids within each peptide. The cation in protamine is exclusively arginine, and the arginine side-chain comprises a guanidino group with a geometry and charge distribution that are particularly favourable for the binding of tetrahedral oxoanions.40 This may make arginine-rich peptides more effective at silicification in vitro. We did not determine the contribution of any sequestered peptide to the weight of the peptide–silica products, but previous reports have suggested that R5 comprises about 25% of the precipitate by weight.20,33
The R5 results are shown in Fig. 2A–C. Fig. 2A shows that silica was only precipitated above 1 mM R5 at pH 7 and 8. No rapid precipitation by R5 was observed under any peptide concentrations tested at pH 6. The reaction appeared to be strongly cooperative and was well fit to the Hill equation with Hill slopes of 4.2 ± 1.0 and 3.5 ± 0.8 at pH 7 and 8, respectively (R2 ≥ 0.9). This is in contrast to previous experiments with native silaffins, where silica precipitation has a linear dependence on peptide concentration.14,16 From the curve-fitting shown in Fig. 2A, the maximal weight of precipitate that could be formed from 50 mM silicic acid was 19.3 ± 1.0 mg ml−1 at pH 7 and 22.5 ± 1.5 mg ml−1 at pH 8, in good agreement with Fig. 1. Half-maximal precipitation occurred at 1.3 ± 0.1 mM R5 at both pH 7 and 8. This R5 data could also be readily fit to a Boltzmann sigmoidal function. Similar behaviour was previously observed in a study of the peptide concentration dependence of R5-induced precipitation from tetramethoxysilane,20 confirming that this is not related to our preparation of the silicic acid precursor, reaction conditions or use of precipitate weight measurements. Although this prior study actually applied a linear fit to such data, we use here a wider R5 concentration range and so are able to confidently assign sigmoidal behaviour. SEM microscopy (Fig. 2B) confirmed that the presence of R5 resulted in the formation of R5–silica nanospheres with an average particle size (dav) of 540 ± 60 nm (N = 25) and EDX analysis of these precipitates (Fig. 2C) showed the expected elemental composition.
Fig. 2D–F show the analogous experiments for protamine. Protamine was again found to be more active than R5 in silicification, generating substantial silica at 0.6 mM peptide (the lowest concentration tested). Protamine exhibits a less pronounced, though still discernible, pH dependence with the highest concentrations of protamine generating ∼1.5-fold more precipitate at pH 7 and 8 than the equivalent reaction at pH 6. In contrast to the results with R5, protamine reactions were not strongly cooperative with lower Hill coefficients of 1.5 ± 0.2 at all pH values. Fitting to a hyperbolic equation yielded maximal precipitation values of 22.5 ± 1.8 mg ml−1 at pH 6, 35.5 ± 1.7 mg ml−1 at pH 7 and 41.4 ± 2.0 mg ml−1 at pH 8, and gave half-maximal precipitation values of 0.8 ± 0.2 mM at pH 6, 1.2 ± 0.2 mM at pH 7 and 1.4 ± 0.2 mM at pH 8. SEM micrographs from the resulting precipitates (Fig. 2E) confirmed that the addition of protamine also resulted in the formation of silica nanospheres; EDX analysis (Fig. 2F) gave the anticipated elemental composition. The average diameter of the silica spheres formed (dav) was slightly smaller for protamine compared to R5, being 380 ± 80 nm (N = 25).
The reactive species for silicification is Si(OH)3O−, giving a linked equilibrium between the ionization of silicic acid and silica precipitation. It seems plausible that increased silicification at higher pH (lower [H+]) arises from a linear increase in Si(OH)3O−. Silicification only occurs above pH 6.4 and given that the pKa for the first ionization state of silicic acid is 9.8, this equates to a ratio of ∼2500:1 Si(OH)4:Si(OH)3O− being the critical threshold condition for rapid silicification. Silica formation by native silaffins and polyamines depends upon the presence of phosphate groups, occurring either as phosphorylated serines introduced via post-translational modifications or through phosphate ions being present in the reaction buffer.16,41,42 This requirement for phosphate does not appear to extend to R5, where introducing phosphoserine residues into the peptide has the opposite effect of inhibiting silicification at neutral pH;22 nonetheless, buffer phosphate may still be required to promote peptide self-assembly by bridging between cationic sidechains. Since the pH range explored in this study does encompass the equilibrium between H2PO4− and HPO42− (pKa = 7.2) it is possible that changes in the protonation state of buffer phosphate could influence silicification across this pH range.
It is also possible that the pH-dependence of silicification arises from changes in peptide charge at increasing pH. However, the pKa values of the lysine and arginine sidechains are expected to be close to the literature values of ∼10 and ∼12 in the absence of significant secondary structure (see below) and will remain fully protonated over the pH range tested here. The terminal amino and carboxyl groups are expected to have pKa values close to 8 and 3, respectively, but the N-terminal amino group alone is not sufficient to induce silica precipitation.32 Thus although we cannot exclude changes in peptide ionization playing a role in silicification, the overall peptide charge is unlikely to change significantly over the range pH 6.0–8.0.
A computational prediction of the secondary structure of R5 was made by submitting the peptide sequence to the PSIPRED server.48 The results, shown in Fig. 4, suggest that R5 is largely unstructured, with the surprising possibility of a short length of β-strand at the C-terminal residues KRR. The confidence in this prediction is low, and it seems unlikely to be correct since independent β-strands are unstable. We extended our computational analysis by using the PEP-FOLD software49 to build multiple de novo models of the three-dimensional structure of R5, and representative models are shown in Fig. 4. The models did not converge upon any stable secondary structure in any part of the peptide sequence and this analysis suggests that the peptide is a diverse random coil ensemble in solution in agreement with experimental circular dichroism (CD) and NMR data presented here. In order to compare these random coils with experimental dynamic light scattering measurements (Fig. 7) we measured the distance between the two furthest separated main chain carbons for each of the models as shown. This measurement gives an approximate cross-section for the random coil of 1.4–2 nm, in good agreement with the hydrodynamic diameter of 2 nm determined using DLS (Fig. 7).
Fig. 4 Bioinformatic predictions of R5 secondary structure. Sequence-based prediction by the PSIPRED server48 gives a low-confidence prediction of a short region of β-strand at the C-terminal. However, de novo structure predictions using PEP-FOLD49 suggest that R5 is a disparate random coil ensemble in solution. Three representative models from the PEP-FOLD output are shown with N-terminal serine and C-terminal leucine indicated as S1 and L19, respectively. Black dotted lines show distances measured in Pymol,52 which was used to construct the figures. |
CD was used to assess whether R5 contained any substantial secondary structure, or whether such structure was induced upon interacting with subsaturating concentrations of silicic acid. The results are shown in Fig. 5. R5 exhibited a classical random coil signature at all pH values tested between 4.5–8.0 (Fig. 5A). Although the signal intensities were slightly varied between samples the shape of the spectra were identical with a single deflection at 198 nm characteristic of random coil.50 The presence of 2 mM silicic acid at pH 7 did not induce any structural change (Fig. 5B). It was not possible to collect data at higher silicic acid concentrations because of turbidity from precipitated silica.
1H solution NMR spectra were also collected to ascertain the secondary structure and aggregation state of R5. The 1D spectrum (Fig. 6A) showed sharp peaks in the aliphatic region (0–5 ppm) of the spectrum as well as clear signals for the aromatic protons of the single tyrosine residue present. At pH 7.0 the amide signals for all but the C-terminal isoleucine (I18) and leucine (L19) were broadened out to the baseline by solvent exchange effects. By comparison, sixteen of the nineteen amide signals were observable at pH 5.5. The aliphatic region of spectra collected at pH 7.0 and pH 5.5 were very similar with the exception of several intensity changes for side-chain βH protons of serine residues that sharpened at low pH due to reduced exchange effects from the side-chain (Oγ) hydroxyl group proton. At pH 7.0, all of the expected aliphatic signals could be identified in the two-dimensional TOCSY spectra (Fig. 6B), confirming solvent exchange as the source of line broadening and ruling out association or dynamic effects. The TOCSY experiment showed CαH and sidechain proton correlations for each amino acid type (Ser, Gly, Tyr, Lys, Arg, Leu and Ile) that were clustered around the expected values for a random coil peptide and did not show the dispersion associated with structured sequences. This apparent lack of secondary structure was also suggested by NOESY experiments (Fig. 6C) that showed no significant short-range or long-range NOEs at pH 5.5 or pH 7.0. At pH 5.5 two short stretches from G6 to Y8 and I18 to L19 could be assigned using sequential assignment methods but severe overlap and lack of dispersion prevented any further assignment.
CαH chemical shifts are strong indicators of secondary structure and the chemical shift index51 predicted random coil structure along the length of the peptide; the lack of variation of CαH chemical shift for each residue type with sequence position meant a chemical shift index prediction could be made without complete assignments.
Identical NMR spectra were recorded in the presence of 50 μM silicic acid and Fig. 6D shows the one-dimensional spectra in the absence and presence of silicic acid for comparison. Because of the long time-scale of the NMR experiments, we empirically determined that 50 μM silicic acid was the highest concentration that could be used without samples gelling or precipitating. Collectively, the NMR and CD data both suggest that R5 is a monomeric random coil in solution and that soluble silicic acid does not induce the formation of secondary structure at the concentrations and pH values examined here. Our data thus confirm that primary amino acid sequence, rather than any peptide structure, is the critical factor in R5-driven silicification.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4tb01679c |
‡ Current address: Department of Chemistry, King's College, London, SE1 1DB, UK. |
This journal is © The Royal Society of Chemistry 2015 |