Open Access Article
Paolo
Passaretti
*a,
Yiwei
Sun
b,
Inam
Khan
c,
Kieran
Chan
a,
Rania
Sabo
a,
Henry
White
d,
Timothy R.
Dafforn
e and
Pola Goldberg
Oppenheimer
*a
aSchool of Chemical Engineering, University of Birmingham, Birmingham, B15 2TT, UK. E-mail: PXP561@student.bham.ac.uk; P.GoldbergOppenheimer@bham.ac.uk; Tel: +44 (0) 778 49 68 145 Tel: +44 (0) 121 41 45 287
bSchool of Physics and Astronomy, Queen Mary University of London, London, E1 4NS, UK
cSchool of Metallurgy and Materials, University of Birmingham, Birmingham, B15 2TT, UK
dBAE-Systems-Air Sector, Buckingham House, FPC 267, Filton, Bristol, BS34 7QW, UK
eSchool of Biosciences, University of Birmingham, Birmingham, B15 2TT, UK
First published on 28th June 2019
Graphene, since its successful exfoliation and characterisation has been continuously drawing extensive research interests due to its potential for a broad range of applications ranging from energy, microelectronics, through polymer fillers and sensors to environmental and biomedical devices. Exploitation of its unique chemical and physical properties for the manufacturing of functional materials, requires careful structural control and scaling-up into three-dimensional morphologies. Here, a facile method is established to create and control the bottom-up self-assembly of graphene oxide nano-sheets via unprecedented integration with a highly versatile bio-ingredient, the filamentous bacteriophage M13, into hierarchical, three-dimensional, porous sponges of GraPhage13. This study explores the interplay of the GraPhage13 structure formation and studies the mechanisms that give rise to the controllable self-assembly. The straightforward fabrication of robust hierarchical micro-nano-architectures further lays a platform for applications in energy storage and conversion, catalysis and sensing.
In this work, we demonstrate a novel route for bottom-up self-assembly of GO into a 3D porous structure spanning the micro- and nanoscales, yielding a unique hierarchical GO sponge based on a synthetic biology approach, by exploiting the bacteriophage M13, acting as a cross-linker between the graphene nanosheets.
M13 is a filamentous bacteriophage that infects Escherichia coli (E. coli) strains showing the conjugative pilus. This virus establishes sustainable infections without killing the host, releasing new copies of itself to propagate the infection to the nearby E. coli cells.26,27 M13 phage has typical dimensions of approximately 880 nm in length and 6.6 nm in width, and it is composed of a circular, single-stranded DNA molecule (ssDNA), encoding ten genes encapsulated in a flexible rod-like shape capsid, consisting of 2700 copies of its major coat protein PVIII. At the ends, there are five copies of each minor coat protein, PIII and PVI on one side and PVII and PIX on the other (Fig. 1b).28–30
Our facile method to fabricate GO-M13 sponges (GraPhage13), is performed without the need of extreme conditions of temperature, pressure or pH rendering it an environmentally friendly and scalable material, further characterised by ultra-low-density and high-surface-area. The fabrication of GraPhage13 is based on four main steps: (i) assembly reaction in aqueous buffers, (ii) hydrogel separation through centrifugation, (iii) drop casting onto a silicon wafer and finally (iv) drying under vacuum (Fig. 1c and d). Furthermore, M13 can be functionalised through the many specific chemical groups on its surface with a broad range of functional molecules such as peptides, antibodies, proteins and enzymes.31–33 These hybrid micro-nano sponges, therefore have potential applications for the development of electrochemical energy storage and conversion devices, absorbers, sensors or composite scaffolds.1
Under these conditions, it was not possible to distinguish between the GO precipitation due to the self-assembly processes or its clustering, therefore, hampering the spectrophotometric analysis of the hydrogel formation. To improve its stability, GO of 5 mg mL−1 from the stock solution, was sonicated for 5 min (GO-5′), obtaining a new stock. Moreover, the concentrations of the samples prepared by diluting both stocks (GO and GO-5′ of 5 mg mL−1), were estimated using calibration curves calculated by measuring the absorbance of serial dilutions (ESI Fig. S1†). The sonication allowed the separation of the layered flakes as well as obtaining smaller GO sheets, achieving an even dispersion (Fig. 2b).34
AFM height images (Fig. 2d–f) and the corresponding cross-sections (Fig. 2g–i) of the dispersed GO, GO-5′ and M13 bacteriophage clearly show that while the GO samples consist of large flakes and aggregates of heterogeneous thickness (Fig. 2d), the GO-5′ is predominantly composed of small, well-distributed nanosheets with a typical thickness of ∼1 nm (Fig. 2e).24,25 The M13 bacteriophages measure, on average, 880 μm in length, 60 μm in width and 3 μm in height (Fig. 2f). The apparent non-circular cross-section of the filamentous bacteriophage is due to a variety of conditions and in particular, the relative humidity of the sample, the applied force on the surface during the scan and tip convolution effects, caused by the curvature radius of the tip itself.35
UV-Vis spectra of GO exhibit a π → π* plasmon peak at 230 nm and a shoulder corresponding to the n → π* plasmon peak around 310 nm (Fig. 2a).36 The spectra of M13, show a high absorbance in the region between 200–230 nm, mostly due to the π → π* transitions in the peptide bonds37 and the typical peak at 269 nm of M13 bacteriophages which is due to the combination of the viral DNA and its major coat protein PVIII.38
The comparison of the absorbance intensities before and after the centrifugation enabled the stability of GO dispersion in DIW to be studied (Fig. 2a–c). It was found that only 7.05 ± 0.80% of the GO remained in the supernatant, while GO-5′, on the other hand, exhibited considerably better stability with a higher spectral absorbance intensity, corresponding to 94.51 ± 1.54% of the GO-5′ remaining within the supernatant. Also, 96.89 ± 6.60% of M13 remained dispersed in the supernatant after centrifugation (ESI Table S1 and Fig. S2†). Therefore, GO-5′ has been identified as the optimal GO preparation for the fabrication of GraPhage13.
To establish the optimal weight ratio between the two components for the formation of GraPhage13 hydrogel (GPH) self-assembly, GO was mixed with M13 in DIW at different weight ratios of 2
:
1, 1
:
1 and 1
:
2. GO concentration was kept constant at 0.3 mg mL−1 in all the samples, while M13 was of 0.15 mg mL−1 (GO-M1303015-5′), 0.3 mg mL−1 (GO-M130303-5′) and 0.6 mg mL−1 (GO-M130306-5′), respectively (Fig. 3).
GO-M1303015-5′ exhibited a minor aggregation and the characteristic spectrum of pure GO in solution, although, due to the presence of M13 in solution, with a slightly increased absorbance at 200 nm. Following the centrifugation, a dark pellet and a transparent light brown supernatant were formed indicating the presence of GO (Fig. 3a). GO-M130303-5′ exhibited an evident formation of the hydrogel (Fig. 3b), which is clearly separated from the supernatant with the corresponding spectrum showing a continuous increase of absorbance due to the macroaggregates of the hydrogel, responsible for the light scattering. The supernatant itself appeared transparent, and the lack of characteristic GO and M13 peaks in the UV-Vis spectrum, suggests the nearly complete precipitation of both components in the pellet. GO-M130306-5′ on the other hand, appeared clear with no apparent hydrogel formation and the corresponding spectrum exhibiting the fingerprint features of GO and M13 individually (Fig. 3c). Furthermore, following the centrifugation, only a small pellet was formed, and the supernatant spectrum was found to completely overlap with the spectrum of the GO-M130306-5′, suggesting that the two components did not interact, and they were still fully dispersed in the solution.
Given that the absorbance of two or more components in solution is derived by the sum of the individual absorbance values,39 the spectra of each sample were compared to the one obtained by the sum of their respective components, measured individually (Fig. 3d–f), to preliminary investigate the interaction between the two components. Interestingly, there was a reduction in the total absorbance in each case and in particular, for GO-M130303-5′ (Fig. 3e). Furthermore, this composition does not exhibit the characteristic peaks of absorbance of the individual GO and M13 components. In contrast, the GO-M130306-5′ spectrum almost entirely overlaps with the corresponding one obtained by the sum of the respective components (Fig. 3f), suggesting that GO and M13 were absorbing the light in the same manner as they would if they were freely dispersed in DIW. These results indicate that GO and M13 in the form of a hydrogel scatter the light, not exhibiting their typical absorption peaks.
The absence of the hydrogel formation in the case of GO-M130306-5′, whereas a pellet and an excess of M13 was expected in the supernatant (Fig. 3c), suggested that pH might play an important role in the GPH self-assembly. Subsequently, the effect of pH was systematically verified by adding set volumes of 1 M NaOH or HCl to the aqueous GO-M130303-5′ solution. While the excess of HCl led to hydrogel assembly, NaOH yielded a reverse outcome of disassembly. Moreover, the pH values of GO-M1303015-5′, GO-M130303-5′ and GO-M130306-5′ in DIW were measured and found to be 3.77 ± 0.02, 4.85 ± 0.02 and 5.7 ± 0.02, respectively (ESI Table S2†).
GO in DIW at a concentration of 0.3 mg mL−1 exhibited a pH of 3.49 ± 0.02, due to its acidic properties given by the presence of phenolic (pKa = 9.8) and carboxylic (pKa1 = 6.6 and pKa2 = 4.3) functional groups on its surface.40 These groups are responsible for its negative charge at pHs ranging from 2 to 12, and a stable dispersion can be typically achieved at pH ≥ 3.40,41 On the other hand, M13 in the same conditions showed a pH of 6.97 ± 0.06 and stability at pH ≥ 4.8. M13 is almost entirely comprised of its major coat proteins PVIII, which is made of 50 amino acids, and 13 of these, including the N- and the C-terminus contribute to the total protein charge (Fig. 4a). The C-terminus and the residues K40, K43, K44 and K48 face the internal cavity of the phage, stabilising the negative charge of the ssDNA during the capsid self-assembly mechanism.42 Moreover, Y21 and Y24 residues, mostly embedded in the protein structure, and exhibiting high pKa values (ESI Table S3†), do not contribute to the total charge of the virus. Therefore, just the N-terminus and the residues E2, D4, D5, K8 and E20 could be responsible for the M13's overall charge distributed on a surface area of ∼18
300 nm2, calculated based on the assumption that the bacteriophage is of a shape of a regular cylinder (Fig. 4a).43,44
![]() | ||
| Fig. 4 M13 PVIII structure analysis and GPHs characterisation. (a) M13 PVIII major coat protein sequence and 3D-structure (top left-hand side). On the right-hand side, a section of the phage structure made of 35 copies of the PVIII, assembled in the quaternary structure (PDB: 2MJZ).45 The hydrophilic surface (green) corresponds to the amino acids Arg, Lys, His, Glu, Asp, Asn, Gln, Thr, Ser and Cys, while the hydrophobic surface (orange) corresponds to the amino acids Ala, Gly, Val, Ile, Leu, Phe and Met. The cavity in two different orientations shows the location of the nucleic acid encapsulated within the M13. (b) The ζ-potential of M13 (red) at pH values ranging from 2 to 11 is overlapping the theoretically calculated values of M13's charge spectra (blue). (c) The amount of GO and M13 left in the supernatant after centrifugation at different buffer conditions, including the DIW and citrate buffer at pH 3.5, 4.9 and 6.9. (d) The comparison between four different GO-M13 formulations at pH 4.9 and 6.9 at which GO and M13 show self-assembly (green) as well as the respective control samples (blue), where there was no observed interaction taking place. (e) The amount of GO (brown) and M13 (grey) left in the supernatant after centrifugation of the four GO-M13 formulations in citrate buffer pH 4.9. (f) Representative AFM height image of the M13 bacteriophages interacting with GO. The unbound phage was washed with citrate buffer (10 mM, pH 4.9). Schematic representation of two different GraPhage13 hydrogel/aerogel compositions: (g) GO- and (h) M13-based. | ||
This hypothesis was further corroborated by measuring the ζ-potential of M13 and comparing it to the theoretically predicted values (Fig. 4b). The charge of each amino acid was assigned at each corresponding pH value, based on its pKa and following the Henderson–Hasselbalch equation.46 The theoretical values of the pKa of each amino acid were calculated by taking into consideration its chemical surrounding (ESI Table S3†). The isoelectric point (pI) of M13 was established as pH 4.1 from both the ζ-potential and the theoretical curves. Similar values were previously measured for fd bacteriophage, which has an identical, mature PVIII protein, excluding the substitution N35D.30,45,47
To establish the pH at which the assembly of GPH occurs, an acid titration, from pH 12 to 2, was performed. At pH 5.4, a minor aggregation started occurring, which further has been dissolved upon mixing followed by a more substantial non-soluble aggregation at pH ≤ 4.8. Notably, the self-assembly occurred at pH values where both components are stably dispersed in solution. The calculated ζ-potential and pH values indicate that the interaction between the two components is most probably due to the reduction of repulsion between the deprotonated carboxylic groups of the GO and the negatively charged amino acids present on M13, at the pH range from 5.4 to 4.8. Moreover, after the protonation of E20, the interaction between the positively charged N-terminus and K8 with the carboxylic groups on the GO can trigger the self-assembly process, further enabling the self-assembly of the GPH. Although this suggests that the electrostatic interactions dominate the self-assembly process, other noncovalent intramolecular interactions involving π-systems are not excluded.
The significant effect of pH on the interaction between the two components was further used to determine the optimal reaction conditions to form the GPH. Therefore, the stability of GO and M13 was tested in 10 mM citrate buffer (CB) at three specific values of pH 3.5, 4.9 and 6.9. UV-Vis analysis demonstrates that both components are stable at all the above-analysed pH conditions, except at pH 3.5, at which M13 has precipitated considerably due to its pI (Fig. 4c). Further UV-Vis analyses were performed to establish and quantify the assembly-ratio between the two components at optimal assembly conditions. GO was transferred into CB and mixed with a 30 times weight excess of M13 (GO-M1300103-5′) and reciprocally, M13 was mixed with 30 times excess of GO inside the buffer (M13-GO00103-5′). Moreover, two other samples were prepared by inverting the order of addition of the two components to evaluate if the mixing was affecting the reaction ratio (M13-GO03001-5′, GO-M1303001-5′). Due to the instability of M13 at pH ≤ 4.8, the self-assembly of GPH was carried-out at pH 4.9, which is slightly above its pI, and compared to the control solutions at pH 6.9, at which no interaction was expected to occur (ESI Fig. S3†).
The UV-Vis analyses further corroborate the significance of the pH on the interaction between the GO and M13 (Fig. 4d). Furthermore, the mixing appeared of significant importance (Fig. 4e) for the self-assembly process.
At optimal pH values for the self-assembly, it is possible to tune the reaction by varying the ratios of each component, to yield the two different assembly conditions called, GO- and M13-based (Fig. 4g and h), depending on which of the two is in excess. GO
:
M13 weight ratio ranges from 4
:
1 to 1
:
5, nevertheless, inverting the order of addition of the two components the M13-based hydrogel ratio can increase up to 1
:
12 (ESI Fig. S3 and Table S4†). These extreme assembly conditions are highly-appealing since they can potentially enable tailoring the properties and functional applications of GraPhage13 aerogel (GPA) such as porosity, surface area and stiffness.
The drop casting (iii) and the drying (iv) steps included the deposition of 50 μL of the hydrogel onto a silicon wafer and the subsequent placing in a vacuum chamber. As the drying process progressed, caused by the vacuum evaporating the water, the hydrogel expanded into an extremely light 3D porous aerogel, showing ultra-low-density of 8.82 ± 0.03 mg cm−3 and high-surface area of 325.147 ± 34.307 m2 g−1 (ESI Fig. S5, Tables S5 and S6†).
Scanning electron microscopy (SEM) coupled with energy dispersive X-ray (EDX) spectroscopy was used to characterise the 3D morphology and composition of the formed GPA (Fig. 5), which exhibited a layered and porous 3D sponge-like structure with a random distribution of micro and nanopores (Fig. 5a–e). GO, on the other hand, exhibited a very fine and fragile structure (Fig. 5f), while the M13 appeared to be of an amorphous morphology with tangled fibres (Fig. 5g). The unique structure of the sponge-like GPA arises from the intertwined M13 nanowires, acting as threads in-between the thin GO nanosheets, forming the hierarchical 3D architectures with embedded microporous morphology. While GO was predominantly comprised of carbon and oxygen, M13 and GPA exhibited the presence of N, S, P, Na and Cl from the amino acids, DNA and traces of salts (Fig. 5i–k). The Si peak is due to the underlying supporting substrate on which the free-standing 3D sponge with an average thickness of ∼2 mm was generated. EDX analysis at three random locations across the macro-sponge yielded identical elemental composition throughout, indicating a homogeneous distribution of GO and M13 inside the 3D sponge (ESI Tables S7–9†). Raman spectroscopy was further carried out on the GO, M13 and GPA to quantify the effects of the interaction between GO and M13 in the sponge on the original bonds of each component. While Raman spectra of GO and GPA exhibited two clear peaks at 1350 cm−1 and 1600 cm−1 at the excitation wavelengths of 633 and 785 nm, accordingly (Fig. 6a), M13 exhibited three strong peaks at 1310, 1650 and 1590 cm−1 (Fig. 6a), agreeing with the literature.48 For the control sample comprised of stacked M13 on top of GO, with no sponge formation but of merely physical contact (PC) between the two components, Raman spectra exhibited highly similar fingerprints to those of GO and GPA (Fig. 6a).
From Raman analysis, it is evident that the spectra of both GPA and GO-M13PC contain contribution mainly from GO rather than the M13. It was anticipated that due to the interaction between GO and M13 during the assembly of the sponge, two main peaks of GPA would be shifted, while those of GO-M13PC will not, indicating the modification of the GO bonds. Considering that all the spectra exhibited broad peaks with a large background, and by employing spectra fitting procedure, a subjective choice of the number of peaks and the background could affect the position of fitted peaks. Baseline-subtraction was not applied to the acquired spectra in Fig. 6a. Maximum likelihood estimation was employed to obtain the peak positions in the model (the peak shape, number of peaks, parameters of each peak and background) optimal for the acquired data ESI Tables S10–13.† Each spectrum was initially subtracted using a straight line for flattening, and subsequently, a corresponding number of Lorentzians and Fourier series of a range of terms were applied to give a fit for each spectrum (Fig. 6b and c).
Finally, by comparing the Bayesian information criteria (BIC) of the fits, the optimal fit was established for each spectrum. To note, no subjective input was provided to the BIC fitting, except setting two rational constraints: amplitudes of Lorentzians are positive, and coefficients of the Fourier series are comparatively small and, therefore, that they do not fit the peaks. The optimal fitting parameters of all the spectra are shown in ESI Tables S10–13.† The BIC fitting shows that the spectrum of the GO-M13PC consists of two peaks, which correspond to the unmodified sp2 (G mode) and sp3 (D mode) bonds of the GO.49 In the spectrum of the GPA, Raman G mode of GO splits, with approximately 30% upshifted and the rest remains unchanged. While the effect of charge transfer usually shifts the entire Raman peak, this observation can be interpreted as approximately 30% of the sp2 bond of GO stiffened by ∼2%, due to the interaction between M13 and GO. These results confirm that the interaction between GO and M13 in the fabricated sponge is strong enough to modify the original sp2 bonds.
295g) and then the pellet containing E. coli was discarded. The supernatant containing M13 was subsequently incubated with PEG 6000 (Sigma-Aldrich®) and stirred for 90 min at 4 °C. Next, the supernatant/PEG solution was centrifuged at 10
000 rpm (Beckman Coulter JA 10®, RCF = 17
649g) and this time the supernatant was discarded. The obtained white pellet mainly consisting of M13, was resuspended in DIW. Finally, the concentration was determined with a UV-Vis spectrophotometer (Thermo Fisher Scientific® Evolution 600) measuring the absorbance at 269 nm and using an extinction coefficient (ε) of 3.84 ± 0.06 cm2 mg−1.29 The stock solution used for all the experiments was of a concentration of 5 mg mL−1.
000 rpm (SciSpin MICRO Centrifuge, RCF = 15
100g) for 1 min.
:
1 and 10 mM APTES.50 Subsequently, two different silicon-APTES wafers were soaked for 3 s in a solution of GO and GO-5′ respectively. Both solutions had a concentration of 1 μg mL−1. 20 μL of M13 25 μg mL−1 in DIW were drop-cast onto a ∼25 mm2 silicon wafer and dried at room temperature. The images were acquired in tapping mode with BRUKER® RTESPA-300 probe (T: 3.4 μm; L: 125 μm; W: 40 μm; f0: 300 kHz; k: 40 N m−1).
000 rpm (SciSpin MICRO Centrifuge, RCF = 15
100g) for 1 min. A P-type silicon wafer (Sigma-Aldrich®) was cleaved (∼25 mm2) and cleaned with a CO2 snow jet on a hot plate at a temperature of 250 °C. After the centrifugation, 900 μL of the supernatant was removed, and the pellet was well resuspended. 50 μL of the final hydrogel were deposited onto the Si substrate and finally, dried for 1 h in a vacuum chamber (Edwards® Rotary vacuum pump E2M2) until reaching a pressure of ∼90 Pa.
437 (cone angle: 20°0′14′′, cone diameter 20 mm and truncation 55 μm) varying the shear rate from 0.1 to 1000 s−1. Moreover, the hydrogel density was predicted (1.003 mg cm−3) and calculated by weighing three samples of 1.0 mL on an analytical balance (OHAUS® Pioneer™ PA114C).
Footnote |
| † Electronic supplementary information (ESI) available: Additional figures and discussion. See DOI: 10.1039/c9nr03670a |
| This journal is © The Royal Society of Chemistry 2019 |