Creating dynamic nanocavity by optical plasmonic trapping for efficient single-molecule SERS detections in solutions
To develop an optical plasmonic tweezer-coupled SERS platform, we constructed a simple plasmonic junction between two AgNP-coated silica microbeads and placed it at the focus of two laser beams to trap an AgNP, which could form a SERS-active nanocavity among the trapped AgNP and the coated AgNP as illustrated in Fig. 1a. Supplementary Fig. 1 shows the experimental set-up combing a 532 nm (2.7 mW) excitation laser and a 1064 nm trapping laser (8.0 mW) via stereo double-layer-pathways in a microscope for simultaneous optical manipulations and spectroscopic measurements. First, silica microbeads were functionalized with amino groups to chemically coat AgNP in sparse distribution (Supplementary Fig. 2), which were further bonded with the 5’-thiol-modified oligonucleotide L1 and the 5’-thiol-modified oligonucleotide L2, respectively. Next, the oligonucleotide L3 was added to tether L1 and L2 by complementary base pairing to generate AgNP-coated bead assemblies, giving a nanogap of around 20 nm. These bead assemblies were then deposited on the bottom of a microfluidic chamber, as visualized in the microscopic images of Fig. 1b (zoom-in) and Supplementary Fig. 3 (zoom-out). Before the spectroscopic measurements, the solutions of analyte molecules and AgNP were injected into the microfluidic chamber in adjacent channels under a consistent flow rate. Upon laser irradiation, a freely diffusing AgNP could be trapped at the plasmonic junction of an AgNP-coated bead dimer to form a dynamic nanocavity among the trapped AgNP and the coated AgNP, sandwiching the analyte molecules in the confined detection volume with a further enhanced local field for sensitive SERS characterizations. Moreover, the AgNP could be trapped with the 1064 nm trapping laser (on state) and released without it (off state), which enables dynamic switching of the trapped AgNP for high-throughput sampling. (Fig. 1)
As a proof of concept, we switched the 1064 nm trapping laser between on and off states and utilized the 532 nm excitation laser to conduct continuous SERS measurements at the plasmonic junction of an AgNP-coated bead dimer in 500 nM Nile Blue A (NBA) solution, where the formation of dynamic nanocavity could be controlled by the optical plasmonic trapping of an AgNP in close contact of the coated AgNP at the junction. The trapping laser was switched 12 times continuously and 60 real-time SERS spectra with 1 s integration time per spectrum were recorded accordingly in Supplementary Fig. 4. By taking the time-averaged SERS spectra at each switching state, Fig. 1c demonstrates a regular pattern that the characteristic peak of NBA at ~ 1645 cm-1 (ring stretching)33 was emerged at on states and vanished at off states. Moreover, the histogram analysis of the peak intensities at each state gives the relative standard deviations (RSD) within 20%. These results verify the reproducibility of creating the dynamic nanocavity on the optical plasmonic tweezer-coupled SERS platform. To further evaluate the effectiveness and robustness of the platform, numerical simulation using FDTD method was carried out (Supplementary Note.1), which shows the electric field distribution of the plasmonic junction from the AgNP-coated bead dimer upon excitation in Supplementary Fig. 6. Since the trapping force is proportional to the square of electric field in Rayleigh regime24,34, a freely diffusing AgNP is attracted towards the highest-intensity in the plasmonic junction. However, the size of the nanogap between AgNP-coated beads is determined by the oligonucleotide linkers (20 nm), which is smaller than the diameter of AgNP (70 nm). Consequently, the trapped AgNP is confined against the coated AgNP in the plasmonic junction with the action force of optical plasmonic trapping and the reaction force at the contact surface as shown in Fig. 1d. The numerical simulation gives a time-averaged total force of 8.8 pN upon the excitation of 532 nm and 1064 nm laser beams at the plasmonic junction. Moreover, the trapping potential distribution with respect to the Z positions in Fig. 1d indicates the depth of potential well at the equilibrium position as 3.5×10–18 J, which is sufficient to overcome the thermal energy of Brownian motion (kBT = 4×10–21 J at room temperature) for the stable AgNP trapping35,36. Besides, optical plasmonic trapping requires less laser power to exert stronger force in comparison with optical trapping, which could avoid potential thermal damages to protect biological samples37. More importantly, this platform enables the mechanical confinement of the trapped AgNP at a nanoscale precision since it is no longer subject to the optical diffraction limit. As shown in Fig. 1e, the dynamic nanocavity among the trapped AgNP and the coated AgNP in the plasmonic junction significantly enhances the localized electric field and provides the enhancement factors up to 109, which is adequate to detect a single molecule38.
To verify the single-molecule sensitivity and identify a single molecule in crowded environments, we exploited bi-analytes SERS approach (BiASERS)39 in the solution phase on the optical plasmonic tweezer-coupled SERS platform. Methylene Blue (MB) and NBA were chosen as the bi-analytes partner because of their comparable Raman cross-section and fair detection probability under 532 nm excitation. Equimolar solutions of each dye (10− 8 M) were mixed and injected into the sample channel in adjacent of the AgNP channel inside the microfluidic chamber. Similar to the above procedures, we switched the 1064 nm trapping laser between on and off states and utilized the 532 nm excitation laser to conduct the SERS measurements under on states at the plasmonic junction of an AgNP-coated bead dimer. Figure 2a displays the SERS mapping out of 3600 spectra continuously collected from dynamic nanocavity in time series, showing the distinct spectral features in the regions of 590–610 cm-1 and 1620–1660 cm-1. Statistical analysis of the data indicated that ~ 5% spectra are considered to be above the cutoff of noise and classified into the catalog of ‘single-MB event’, ‘single-NBA event’, or ‘dual-MB and NBA event’, based on the characteristic Raman peaks of MB or NBA at a high concentration of 10− 5 M in Supplementary Fig. 7. Specifically, the appearance of the Raman peak at 1650 cm-1 accompanied by the emergence of a peak at 600 cm-1 represents the detection of ‘single-NBA’ event (Fig. 2b blue), while the appearance of a peak at 1630 cm-1 and the vanish of the peak at 600 cm-1 represents the detection of ‘single-MB event’ (Fig. 2b red)33. ‘dual-MB and NBA event’ refers to the presence of the dual peaks at around 1640 cm-1 (Fig. 2b black). Representative full spectra of each event are shown in Supplementary Fig. 8, and the fluctuation of peaks position was analyzed in Supplementary Fig. 9 to facilitate accurate identification among different events. The statistical histogram in Fig. 2c is predominated by ‘single-dye-event’ (NBA or MB), suggesting that the signals were highly likely originated from a single molecule at the dynamic nanocavity. Furthermore, we repeated the BiASERS experiments by increasing the concentration of MB by 10 times (10− 7 M) while keeping the concentration of NBA (10− 8 M) unchanged. Most of the spectra contain the characteristic peaks of MB at 1630 cm-1 with a large number of MB as the overflowing background molecules. Nevertheless, we still captured a few spectra in which the subtle peaks of NBA at 1650 cm-1 emerged from the predominant peaks of MB at 1630 cm-1 (Fig. 2d), implying the presence of an NBA molecule in the MB-dominated nanocavity40. As the low-populated species, NBA was differentiated from the dominating MB at the single-molecule level. Hence, the optical plasmonic tweezer-coupled SERS platform could achieve the single-molecule SERS detection and differentiation among overflowing background molecules in aqueous environments.
Monitoring the dynamic pH-dependent structures of tyrosine in solution
To further verify the aqueous stability and compatibility, we employed the optical plasmonic tweezer-coupled SERS platform to characterize tyrosine (Tyr) in different pH environments, considering that the side chains of amino acids dictate the protein structure–function relationships while their hydrogen bonding and electrostatic interactions are affected by environmental factors, such as pH41. The real-time SERS spectra of Tyr were monitored by the dynamic nanocavity during the change from acidic to basic conditions, as presented in Fig. 3a. At the beginning, the microfluidic chamber was filled with 50 µM Tyr solution at pH 1, generating the characteristic peak at 1620 cm-1 as the in-plane ring stretching (v8a mode) of Tyr at + 1 charged state42,43. Then 2 M NaOH was added through a side channel to adjust the environments from pH 1 to pH 13, meanwhile SERS spectra were continuously collected. There is a downshift of v8a mode to 1602 cm-1 along with the gradual conversion of Tyr to -2 charged state42,43 and a change in the intensity ratio of Tyr doublet located at ~ 854 cm-1 and ~ 830 cm-1, which is assigned to the Fermi resonance between the symmetric ring-breathing (v1) and the overtone of the out-of-plane ring deformation (2v16a)44. The intensity of 830 cm-1 increases and the ratio of I830/I854 changes from approximately 1:1 to 2:143,44. These differences are attributed to the increase of negative charge on the phenolic hydroxyl group (pKa = 10.5) upon the adjustment from acidic to basic conditions44. Other changes, such as the increasing intensity of the peak at 1071 cm-1 (C-N stretching) and the peak shifting from 1153 cm-1 to 1161 cm-1 (C-C-N asymmetric stretching), are also observed and consistent with literatures45. On the other hand, there is no obvious signal at 930 cm-1 (C-COO- stretching), 1390 cm-1 (COO- symmetric stretching), and 1090 cm-1 (amino group vibration) related to Ag-Tyr interactions via the carboxylate and amino groups46,47, owing to the attachment of the nucleotides on the AgNP surface as a protected layer for Tyr. Although the oligonucleotide linkers might be liable to depurination and cleavage at extreme conditions, the geometry and interparticle distance between the AgNP-coated beads were fixed upon deposition in the microfluidic chamber. Moreover, at the midpoint of the pH adjustment process, the spectra measured with a long integration time (10 s) show the spectral broadening of v8a mode at 1610 cm-1 in Fig. 3a (green), while the spectra acquired with a short integration time (1 s) exhibit sharp and fluctuated peaks in the range of 1620 to 1602 cm-1 in Fig. 3b, indicating the dynamic snapshots of single or a few Tyr molecules with distinct charge states at the dynamic nanocavity analogues to the previous single-molecule SERS studies47.
To understand the chemical environment of Tyr in aqueous solutions, we compared the SERS spectra of Tyr with the calculated Raman spectra of various Tyr-H2O clusters in Fig. 3c. Using density functional theory (DFT) simulations, the construction of Tyr-H2O clusters was initially based on B. Hernández’s work by involving seven water molecules in the vicinity to coordinate with the amino group (pKa = 9.2), the carboxyl group (pKa = 2.2), and the phenol hydroxyl group (pKa = 10.5) of Tyr (Supplementary Note. 2)48, then their structures were further optimized at + 1 and − 2 charged states to match the experimental results at pH 1 and pH 13, respectively. Supplementary Fig. 10 presents possible hydrogen-bonds between Tyr and water molecules and the charges of Tyr nitrogen and oxygen accordingly. The electrostatic potential (ESP) mapping49 in Figs. 3d and e were generated from the interaction energy between a unit charge on 0.01 a.u. iso-surface and the Tyr molecule, demonstrating the charge distribution of Tyr with the electrophilic moieties in red and the nucleophilic moieties in blue. These results indicate that both the phenol hydroxyl group and the carboxyl group of Tyr become negatively charged after the deprotonation at pH 13, which depolarizes the whole molecule and increased the Raman activity of in-plane ring stretching to account for the downshift and enhancement of v8a mode43,45. Moreover, the hydroxyl group serves as both hydrogen-bond donor and acceptor in acidic and neutral environments, while its oxygen is exposed as a sole hydrogen-bond acceptor under basic conditions. The loss of hydrogen and the increase in negative charges on phenolic oxygen lead to the enhanced intensity at 836 cm-1 and the unequal intensity ratio of the Tyr doublet44. With rich pH-dependent vibrational characters, Tyr could be an indicator to reveal the structural detail and micro-environment of biomolecules in solution by the optical plasmonic tweezer-coupled SERS platform.
Characterizing transient species of hIAPP under different physiological pH
With efficiently reduced detection volume and integration time, we utilized the optical plasmonic tweezer-coupled SERS platform to probe the SERS signal of hIAPP at two physiological pH conditions (pH 5.5 and pH 7.4) in the manner of small-size sampling to overcome the ensemble and temporal averaging of heterogeneous species23. As an amyloidogenic IDP, free hIAPP monomers possess the propensity to assemble into heterogeneous aggregates and amyloid fibrils10. Its ensemble oligomers in dilute solutions and amyloid fibrils in solid forms are well characterized by SERS studies50,51. Yet, it is unclear how the environmental factors affect hIAPP intermediate species and aggregation pathways52,53. Specifically, low pH is documented as the inhibitor against amyloid formation6. The CD spectra of 100 µM hIAPP in the PBS buffer of pH 5.5 at different incubation times (0 h, 2 h, 24 h) present a persistent negative peak at around 203 nm in Fig. 4a (black, green, blue), indicating that the ensemble conformation is dominated by random coils6. Whereas, the CD spectrum of hIAPP after the 24 h incubation under pH 7.4 displays a broad peak with negative ellipticity centered at 218 nm in Fig. 4b (blue), suggesting the formation of β-sheet structure6. Furthermore, we incubated hIAPP at pH 7.4 for 2 h, then adjusted the pH back to 5.5 for continuing incubation to obtain the CD spectra in Fig. 4c. Interestingly, the β-sheet structure was still detected at t = 24 h after the post-adjustment of pH from 7.4 to 5.5, implying that the critical intermediate species associated with amyloid developments already formed before the pH adjustment at the timepoint of t = 2 h. However, no considerable change was observed in the CD spectra of hIAPP at pH 7.4 from t = 0 h to t = 2 h, regarding as the macroscopic lag phase of hIAPP fibrillation54. Thus, to distinguish the low-populated species buried under the ensemble averaging in quantity and time, we utilized the dynamic nanocavity to characterize the structural features of hIAPP in a small sample size at different incubation time-points under pH 5.5, pH 7.4, and post-adjustment of pH from 7.4 to 5.5 conditions.
Following the same experimental procedure in the previous Tyr section, 10 µM hIAPP solution was injected into the microfluidic chamber with AgNP-coated bead dimers deposited at the bottom, which is adjacent to the AgNP channel to allow the formation of the dynamic nanocavity and sandwich hIAPP species for SERS characterizations. This concentration mimics a near-physiological milieu involving moderate intra- and inter-molecular interactions that are crucial to initiate the aggregation process and stabilize the transient intermediates of hIAPP55–57. Although it is below the detection threshold of spontaneous Raman spectroscopy (Supplementary Fig. 11), we have obtained the clear SERS signals of hIAPP on the optical plasmonic tweezer-coupled SERS platform. Figure 4d shows the representative SERS spectra of hIAPP under pH 5.5 incubation at t = 0 h (black) and t = 2 h (green). The characteristic peaks are assigned to protein backbones including amide I band (1656 cm-1), CH2 deformation (1450 cm-1) and amide III band (1250 cm-1 and 1287 cm-1), as well as specific residues such as Phe (1006 cm-1, 1585 cm-1) and Tyr (830 cm-1 and 850 cm-1, 1605 cm-1)58,59. The details of peak assignments are listed in Supplementary Table. 1. To better analyze protein secondary structures, we conducted the second derivative analysis59 for the amide I region from 1550 cm-1 to 1750 cm-1 in Fig. 4e to plot Fig. 4f. The spectra at t = 0 and t = 2 h exhibit the identical amide I band centered at 1656 cm-1, which are attributed to the persistence of the helical and disordered hIAPP50,51. It is further supported by the amide III bands at 1250 cm-1 and 1287 cm-1, assigned to random coil and α-helix structures58, respectively. Since hIAPP has four positive charged sites (α-amino group, K1, R11, H18) at pH 5.5, as illustrated in Fig. 4g, the electrostatic repulsions would give rise to the swollen coil conformation and plausible inhibitory effect on fibril development6,60,61. The co-existing α-helical structural component could be attributed to the moderately amphipathic segment 8–16 (Fig. 4g, dashed box), based on previous studies60,62. In addition, the Tyr doublets with the value of I830/I854 ≈ 1 in the SERS spectra of hIAPP under pH 5.5 (Fig. 4d) imply that the C-terminal Tyr residue acted as both hydrogen-bond donor and acceptor in exposure to water molecules44,59, in consistent with the extended and coil-rich conformations of hIAPP in the physiological acidic environment63,64.
The representative SERS spectra of hIAPP under pH 7.4 incubation at t = 0 h and t = 2 h are shown in Fig. 4h, followed by the plots of their amide I region from 1550 cm-1 to 1750 cm-1 in Fig. 4i and the corresponding second derivative spectra in Fig. 4j, respectively. At the beginning, the spectrum of hIAPP at pH 7.4 (black) demonstrates the amide I band at 1656 cm-1, which is attributed to the helix-coil structure and similar to the initial conformation at pH 5.5. However, after 2-hours incubation at pH 7.4, we observed three types of SERS spectra of hIAPP from the small-size sampling measurements at the dynamic nanocavity with 1 s accumulation time: the predominant type and two rare types. The most frequently appeared spectrum (Fig. 4h, green) exhibits a typical amide I band at 1656 cm-1, representing the predominant helix-coil structure of hIAPP at 7.4. Intriguingly, two distinct types of spectra were obtained occasionally and reproducibly, as shown in Fig. 4h, purple and red. The type I of rare spectra (purple) features the broaden amide I band at 1668 cm-1, attributed to an emerging partial β-sheet conformation51. Moreover, the intensity ratio of the Tyr doublets I830/I854 ≈ 1 in the type I of rare spectra indicates the exposure of the C-terminal Tyr-37 in an extended conformation44,59. While the type II of rare spectra (red) shows the characteristic amide I band at 1655 and 1674 cm-1, implying the co-existence of helix-coil and β-sheet structures50,51. In addition, the Tyr doublets with the intensity ratio I830/I854 > 1 suggest that the C-terminal Tyr-37 slightly shifted toward hydrogen donating as interrelated to its surrounding chemical environment (i. e., buried hydrophobic core, as shown in Fig. 4k)59,65,66. As evident from these spectral features, the type I and type II of rare spectra might be putatively assigned to the intermediate species of hIAPP containing the partial β-sheet with the C-terminal in extended or bent conformations, respectively. Since His-18 (pKa3 = 6.0) deprotonated at pH 7.4, it reduces the electrostatic repulsion of the whole hIAPP molecule to promote intra- and inter-molecular interactions61. As a result, the hIAPP transient species with a bent C-terminal could be stabilized by the special hydrogen bonding environment near His-18 and the ring stacking between aromatic residues, which is indicated by the distinct change of Tyr doublets in the type II of rare spectra (red)43,67. This assignment is supported by the previous studies that the interaction between His-18 and Tyr-37 would minimize the entropic cost of forming a transient β-sheet in the amyloidogenic region (roughly S20-S29, Fig. 4k, dashed box)60, 68–70. Although existing in a low population, the two β-sheet-containing transient species of hIAPP were directly differentiated from the dominating helix-coil species of hIAPP at the early lag phase of the amyloid formation at pH 7.4, which reveals the influence of pH on the structural conversions of hIAPP and demonstrates the effectiveness of our SERS platform in characterizing a few heterogeneous proteins involving in molecular interactions under physiological pH environment.
After the 24 h-incubation at pH 5.5, the SERS spectrum of hIAPP in Supplementary Fig. 12 demonstrates similar spectral characteristics to the spectra at t = 0 and t = 2 h, indicating that most hIAPP maintained the helix-coil structure in consistent with the ensemble CD spectroscopy analysis. In contrast, hIAPP fibrils were observed and confirmed by AFM image (Supplementary Fig. 13) at t = 24 h of the incubation at pH 7.4. No obvious SERS signal of hIAPP was detected at the nanocavity, due to the size limit50. We then performed the bulk characterization of hIAPP fibrils using AgNP colloids to serve as a reference. Supplementary Fig. 14 clearly demonstrates the amide I band at 1672 cm-1 and amide III band at 1226 cm-1, which are attributed to the ordered β-sheet structure in hIAPP fibrils51,71. Furthermore, identical amide I band at 1674 cm-1 and amide III band at 1226 cm-1 were observed from the SERS spectra of hIAPP incubated at pH 7.4 for the first 2 h then at pH 5.5 for the remaining 22 h. The SERS spectrum and AFM image (Supplementary Fig. 14) imply the formation of hIAPP fibrils with the ordered β-sheet structure. More importantly, the spectral features of hIAPP further verifies the assignments of the type I and type II of rare spectra of hIAPP under pH 7.4 incubation at t = 2 h to its β-sheet containing transient species. Since amyloidogenesis is driven by nucleation60,72, the pH-regulated formation of these intermediates might slightly shift the dynamic conversions of hIAPP species in amyloid-competent and alternative conformations among the predominant helix-coil structures to facilitate the subsequent aggregation even after the post-adjustment of pH from 7.4 to 5.5. Figure 4l illustrates the influence of pH on the amyloid developments of hIAPP, based on the small-size sampling results on our SERS platform. At t = 0 h, hIAPP possesses helix-coil conformations, which is preserved without apparent aggregation during the incubation under pH 5.5 for 24 h. The pH 7.4 environment maintains the predominant helix-coil hIAPP in the macroscopic lag phase, but induces the formation of the low-populated transient species of hIAPP containing the partial β-sheet with the C-terminal in extended or bent conformations during 2 h-incubation, which eventually converts into amyloid fibrils. The slight shift in the equilibrium between different hIAPP species at the early stage of incubation under pH 7.4 subsequently results in irreversible amyloid formations even after the post-adjustment of pH from 7.4 to 5.5. It is worth noticing that the environmental pH plays a vital role in the structural conversions of hIAPP, where His-18 (highlighted in yellow in Fig. 4g and 4k) is deprotonated in the change from pH 5.5 to pH 7.4 to reduce its electrostatic repulsive interactions. However, the formation of low-populated hIAPP transient species drives the subsequent amyloid development irreversibly and independently upon pH post-adjustments, which is highly likely mediated by the C-terminal residues (Tyr-37, highlighted in cyan in Fig. 4g and 4k) in addition to His-1868,73. Thus, the previous strategies that focused on a single site of hIAPP might not be optimal in blocking hIAPP amyloidosis, whereas targeting its on-pathway transient species could be a promising strategy74,75. Furthermore, the direct structural characterizations of the early-stage transient species of hIAPP involving in intra- and inter- molecular interactions at physiological concentrations would provide profound mechanistic insights to understanding the protein folding and assembly in amyloid formation, which also showcases the power of the optical plasmonic tweezer-coupled SERS platform in detecting low-populated species in heterogeneous mixtures.