3.1 Highly specific and efficient GPC3-binding DNA aptamers
Based on our aptamer screening procedure, specific GPC3 binding aptamers were carefully isolated in this study. Single-strand DNA aptamer candidates against the recombinant GPC3 protein were generated using the modified SELEX protocol [40]. The details of this protocol are available in the Supplementary Information section. The eluted ssDNA concentration in each selection round was monitored using a NanoDrop spectrophotometer, which reflected the binding capacity between the selected ssDNA pool and the GPC3 target. ssDNA concentrations in each round increased until the sixth round of selection (Figure S1). Changes in ssDNA concentration indicate that the ssDNA pool was enriched after each selection round. After a negative round, ssDNA aptamer concentrations of the eighth, ninth, and tenth selection rounds were 488.4, 491.9, and 405.9 ng µL−1, respectively. At the tenth round, the ssDNA aptamer concentration decreased, suggesting that the ninth round was the optimal state. The GPC3 bound aptamer pool was amplified, cloned, and sequenced. A total of 36 clones were selected for DNA sequencing and 12 different DNA sequences were identified (Table 1).
Table 1
Sequences and affinity value of GPC3 aptamers. Affinity value derived ka, kd, and KD values of 12 selected GPC3 aptamers by measuring the surface plasmon resonance (SPR; Biacore X100)
Clone | Selected Sequence | ka (1/Ms) | kd (1/s) | KD (M) |
GPC3_1 | CATGATCAATCCCGTACATATTTTATCCTCAATTTCACA | 2.09 × 104 | 6.45 × 10−5 | 3.09 × 10−9 |
GPC3_2 | ATTTCATAGTGTTCTGTTTTTTTCCAGTTTTCTTATGTCG | 2.28 × 105 | 9.84 × 10−4 | 4.31 × 10−9 |
GPC3_3 | CCTATTCCTTATTATATTTTCTTTTTTTGTAATTTGGTCG | 3.67 × 105 | 4.40 × 10−5 | 1.2 × 10−10 |
GPC3_4 | CCATAATACTGTTTTTTCCGCTGCTATTAGTTACCTCACG | 1.11 × 104 | 7.14 × 10−5 | 6.46 × 10−9 |
GPC3_5 | GCAATCATATGATTTTCCATTACCTCGTTTTTACTTTATT | 2.20 × 105 | 8.53 × 10−4 | 3.88 × 10−9 |
GPC3_6 | TAGTTCTTTTAGGGGTTTTTTGTTTGTCTATGTGTTGTGG | 2.00 × 105 | 1.35 × 10−4 | 6.71 × 10−10 |
GPC3_7 | ATCTTTTAGTTTGTTTGTATATGTTTTTCTCATTTTCTCG | 5.94 × 103 | 6.30 × 10−5 | 1.06 × 10−8 |
GPC3_8 | ATACATTACTATACTGTGCCATTTTTTATATCCGTTCG | 2.75 × 104 | 5.00 × 10−5 | 1.82 × 10−9 |
GPC3_9 | ACTTCCTTTTTTTCGTCTTGTTATCGTTTCCTTATTCTCG | 2.5 × 105 | 5.17 × 10−5 | 2.07 × 10−10 |
GPC3_10 | TTCTTGTCATAAATATTTCCTATTTCCATTCGTACTCCG | 1.72 × 104 | 6.12 × 10−5 | 3.57 × 10−9 |
GPC3_11 | TTTCTTTTGCGTATTATTTACCTTATTTTCTTAGTGTG | 1.80 × 106 | 9.80 × 10−4 | 5.43 × 10−10 |
GPC3_12 | CTTTTTTCGTGTTTTTAGCTCTTTCCGAGTCTTTTCGAAC | 1.09 × 106 | 8.91 × 10−4 | 8.16 × 10−10 |
To determine the characteristics of GPC3-binding ssDNA aptamer candidates, an SPR assay was used to analyze the specific binding of 12 aptamers with different sequences to GPC3 and determine the binding affinity (KD values) of the ssDNA aptamer (Supplementary materials). KD values of ssDNA aptamer candidates were measured using BIA evaluation software (GE Healthcare). The results indicated that the KD values of these ssDNA aptamer candidates could be classified into three distinct groups. Of these ssDNA aptamer candidates, GPC3_1–GPC3_12 had low KD values (Table 1). The KD values of ssDNA aptamers ranged from 1.2×10−10 to 1.06×10−8 M. In addition, the results obtained from SPR confirmed that GPC3_3 had the highest affinity for GPC3 based on KD values (1.2×10−10 M), which would aid the detection of GPC3 targets.
3.2 3D structure analysis of the molecular interaction between DNA aptamer and GPC3 protein
The GPC3 protein is a member of the glypican family, which is attached to the cell surface and has a glycosylphosphatidylinositol (GPI) anchor region [41]. The protein comprises two subunits, the ~40 kDa α-subunit (light blue) and the ~30 kDa β-subunit (deep blue), which is cleaved between Arg358 and Ser359 [42] (Figure S2A and S2B). The GPC3 protein structure was predicted using the Glypican-1 structure model, 4YWT, in SWISS-MODEL by piecewise homology (Figure S2C). The GPC3 protein structure was used to determine the structural interactions with the DNA aptamer. The important GPC3 protein residues are shown in Figure S2B and S2C. The predicted GPC3 model consists of 16 α-helices and three major lobes [43].
These aptamers have complex secondary structures, including protruding loops and stems. The interaction conformation between the GPC3 protein and each aptamer candidate is shown in Figure S3A. The GPC3 protein was fixed in the same position for analysis, and each aptamer binding site of the GPC3 protein is shown in a different color. We also analyzed the detailed GPC3 protein–aptamer binding interaction using COCOMAPS/CONSRank and LigPlot+, and the amino acid sequence of the GPC3 protein is shown in Table S1. We further analyzed the distance of contact sequence between the GPC3 and 12 aptamer candidates, and indicated them using the heatmap; dark purple color indicates > 4 Å distance, while the bright color indicates < 4 Å distance or closer [44] (Figure S3B). All aptamer candidates showed a strong bright color in distinct areas of the GPC3 protein sequences, showing that the interaction between the GPC3 amino acid residues and DNA aptamers is in the proximation position. We analyzed the heatmap to measure the binding distance between the GPC3 amino acid residues. Heatmaps representing intermolecular contacts for individual residues indicated that frequent contact of the GPC3/GPC3_3 aptamer (green) interaction was comparable to that of the GPC3/GPC3_9 aptamer (magenta) interaction. However, the interaction mode of the GPC3/GPC3_2 aptamer (light purple) was completely different from those of the GPC3/GPC3_3 and GPC3/GPC3_9 aptamers (Figure 1A). To characterize the GPC3–DNA aptamer interaction comprehensively, visual contact performances using the COCOMAPS web server were obtained (Figure 1B). Despite the nucleotide sequence differences, two DNA aptamers (GPC3_3 and GPC3_9) bind to GPC3 with similarly. In addition, we studied the relative contribution of hydrogen bonds and hydrophobic interactions of GPC3 at the DNA aptamer binding sites. All 12 DNA aptamer candidates were analyzed for hydrogen bonding at the GPC3 binding interface using LigPlot+ software (Table S1). The average distance of the hydrogen bonds is shown in Figure 1C. This result demonstrates that the GPC3_3 aptamer had the shortest distance, which may help stabilize the GPC3 protein–aptamer complex by mediating strong hydrophobic interactions. Hydrophobic interactions can influence binding affinity and molecular functional activity [45, 46]. These results correspond well with those described in the SPR affinity analysis (Table 1). Based on the GPC3_3 structural interactions, which the binding region is comparable to GPC3_3, and GPC3_2, which is totally different, referred to “Aptaprobes.” We then attempted to determine the functional binding of GPC3 aptaprobes (GPC3_3, GPC3_2, and GPC3_9) and confirmed their structure-based functional differences. As shown in Table S1, the results of the structural interaction of the aptaprobes binding with GPC3 protein involves 3D interaction between the binding site on GPC3 and the nucleotide paratopic sequences (nt-paratope), the binding site on the aptaprobes. The binding geometry of GPC3_3 aptaprobes indicated that five nt-paratopes (T3, C13, A25, A97, and C99) were responsible for binding to the GPC3 aptatopes: the specific amino acid residues, Lys453, Lys467, and Val532, have a significant interaction with the 5′-end nt-paratope and Asp193 and Cys197 bind with the 3′-end nt-paratope (Figure 1D). GPC3 consists of an N-terminal secretory signal peptide and a C-terminal core domain. The GPC3 core forms a complex with Wnts, activating downstream signal transduction and stimulating HCC proliferation [47]. GPC3_3 aptaprobes bind directly to the GPC3 core region, suggesting that the aptaprobe may alter the GPC3–Wnt signaling pathway. Since GPC3 expression increases in HCC, aptaprobes can be suggested as potential therapeutics for HCC.
Figure 1E illustrates that the GPC3_2 aptaprobe binds to GPC3 at two binding sites different from GPC3_3 binding sites (Figure 1D). The GPC3_2 aptaprobe comprises two physically connected paratopic moieties (six nt-paratopes of GPC3_3 aptaprobes: T72, T78, T80, A4, A87, and G88) that can simultaneously interact with two different aptatopes, His412, Asn127, Asn124, His121, and Gln106 (Table S1), respectively. On the contrary, GPC3_9 aptaprobe and GPC3 binding are mediated by the interaction of the nt-paratope and aptatope interface, which partially overlap with those of the GPC3_3 aptaprobe (Figure 1E and F). The three-region nt-paratope of GPC3_9 aptaprobe from nucleotides C62, C93, and A94 generates specific GPC3 binding, corresponding to the aptatopic residues Lys247 and Thr4. In addition, the nt-paratopes A4, C8, C13, and T16, and linear aptatopes, Cys197, Arg199, and Lys347, are linked to each other. In particular, aptatopesLys486 and Asn554 aptatopes in the GPI region of the GPC3 protein interact with A25 and A30 nt-paratopes.
Structural analysis of the molecular interaction was performed using MOE-docking-based 3D prediction to determine whether the two aptaprobes occupy distinct GPC3 sites during binding. If an identical binding site is involved, the aptaprobe–GPC3–aptaprobe sandwich configuration would be disturbed. The data on non-ligand residues involved in hydrophobic contact of GPC3 protein–aptaprobe analysis (Table S1) indicate that GPC3_2 and GPC3_3 aptaprobes have different binding sites, while GPC3_9 and GPC3_3 aptaprobes have similar binding sites. Figure S4 depicts the interaction between the GPC3 protein and different aptaprobes (GPC3_3, GPC3_2, and GPC3_9). The GPC3_2 aptamer binds to different setsamino acid residues on GPC3 (Figure S4A and S4B), but the GPC3_9 and matchedGPC3_3 aptaprobes bind to 17 common amino acids on GPC3 (Figure S4A and S4C; Met1, Ala2, Leu192, Asp193, Ile194, Glu196, Cys197, Phe445, Pro455, Ser460, Glu490, Leu523, Tyr528, Leu530, Asp531, Asp533, and Asp534). Taken together, GPC3_2 and GPC3_3 were finally chosen for ALISA as capture and reporter, respectively. The results of this study highlight that the molecular docking can successfully predict the interaction between nt-paratopes and aptaprobes, which is promising for high-throughput aptamer screening. We present three GPC3 detection methods using: GPC3_3, GPC3_2, and GPC3_9 aptaprobes.
3.3 Aptablotting assay for the GPC3-aptaprobe interaction
The GPC3 protein is gradually expressed in the serum and tissues of HCC patients [48]. Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS–PAGE) is the most popular biomarker detection method, which identifies a specific biomarker protein in a complex matrix, such as tissue lysate, blood, or other clinical samples. This motivated us to use aptaprobes for GPC3 detection. Aptablotting has previously showed that the aptamer can help overcome the obstacles of current immunoblot assays [28]. In this study, we designed an aptablot assay for GPC3 detection. GPC3_3, GPC3_2, and GPC3_9 aptaprobes were chosen for establishing an aptablotting assay and chemically synthesized as biotin-labeled affinity aptamers (Figure 2A). Subsequently, we tested whether the aptamers non-specifically react with various proteins, such as serum proteins, human serum albumin (HSA), and bovine serum albumin (BSA). The purified recombinant protein GPC3 was used as a control (Figure 2B). The GPC3 proteins and test samples were separated using 10% SDS–PAGE and then transferred onto a nitrocellulose (NC) membrane. Following a simple blocking step, the aptablotting assay required only a single aptaprobe affinity detection using each biotinylated aptaprobe. It is important to note that three GPC3 aptaprobes could recognize their target in both SDS-denatured environments and complex biological samples. In contrast, GPC3 aptaprobes did not react with the serum abundant proteins HSA or BSA.
The aptaprobes could recognize GPC3 in its native state. Figure 2C and Figure S5A show the specificity and sensitivity of the aptamer probe to GPC3 concentration in the dot blot assay. The GPC5 protein, a member of the glypican family protein, is not specific to antibodies and aptaprobes with specificity [49]. When GPC3 was included to represent the complex protein mixture (human serum: H. serum; bovine serum: B. serum), it was clearly demonstrated that the aptamer-bound GPC3 proteins selectively and the observed binding signals were not due to non-specific adsorption to either NC membrane or abundant serum proteins. However, anti-GPC3 antibody showed a strong signal in human serum without GPC3, which was generated by the reaction of human serum with second antibody for detection because the first GPC3 antibody was sourced from humans. Finally, all aptaprobes, GPC3_3, GPC3_2, and GPC3_9, exhibited similar binding behaviors toward native GPC3 proteins. In particular, GPC3_3 aptaprobe reacted with 1 pmol GPC3 in both H. serum and B. serum.
Realizing that the isolated aptaprobes tended to recognize their target structure well in SELEX conditions where they were originally isolated or not, the aptablotting performances in SDS-denatured environment were far from their appropriate condition. Some proteins are known to be resistant to SDS-induced denaturation [50, 51]. It has been shown that aptamer binding was observed in the SDS–PAGE and blotting assays, suggesting that insufficient SDS amount and intrinsic structural rigidity contribute to SDS resistance [52, 53]. We concluded that the high level of GPC3 structural conservation from SDS denaturation may be used as the basis for the aptablotting assay to identify both SDS-denatured and native GPC3 proteins.
The aptablotting assay can directly detect GPC3, whereas western blotting detection by antibody-based experiments is generally stepwise with incubation and labeling, such as primary and secondary antibody use. In addition, negatively charged nucleic acid aptamers bind poorly to the PVDF/NC membrane, indicating that the aptaprobe-involved aptablotting assay can be improved and contribute to their ability to detect and image targets with various signal read-out fluorescent or colorimetric tags. It is expected that the immunoblotting with primary and secondary antibodies might lead to several non-specific signals in addition to the GPC3 signals (Figure S5B, see red arrow). This means that the GPC3 aptaprobe should not cross-react with non-target molecules in the sample mixture, and our aptablotting protocol simplifies the immunoblotting workflow. To be an effective cellular target binding reagent, an aptaprobe should ideally be achieved following binding to its aptatopes expressed on the cell surface of various cell lines.
3.4 Quantitative-ALISA platform for detecting GPC3
Interest in the use of aptamers in ALISA is on the rise [8, 30, 35]. ALISA relies on the formation of a sandwich configuration between two layers of the aptaprobe pair (capture and reporter probes), which binds to the target to create a sandwich complex. This indicates that aptaprobes should bind to distinct aptatopes, and their target recognition should not be mutually exclusive. To identify the aptaprobe pair, individual aptatopes should be tested pairwise in a sandwich matrix so that each aptatope is occupied with other candidates, both as capture and reporter aptatopes. However, this conventional method is extremely time-consuming and labor-intensive to complete, and less time-consuming approaches are required. It should be noted here that molecular docking can be considered as a valuable strategy to provide high-quality and functional aptaprobes and contribute to understanding the binding configuration of aptaprobe-target ligands.
ALISA-quantitative detection was performed as shown in the sandwich assay platform, which is easy to use and does not require any coupling of stepwise secondary antibodies for signal readouts (Figure 3). To confirm the use of aptaprobe pairs for detecting GPC3 in ALISA, three blanks (C1: microwell plate, C2: non-specific reporter aptaprobe binding, and C3: non-specific target binding without capture aptaprobe), and two negative controls (C4; biotinylated GPC3_9 and fluorescein amidite (FAM)-labeled GPC3_3, and C5; biotinylated GPC3_2 and FAM-labeled GPC3_3) were used. The experimental test groups were set as [Test1] to [Test4] (Figure 3A).
The normal sandwich assay platform was subjected to quantitative analysis after treatment with 5, 10, 25, 50, 100, 250, and 500 ng in 100 µL GPC3 protein. Therefore, the fluorescence signal combined with the reporter aptaprobe was confirmed only in the GPC3_2–GPC3–GPC3_3 configuration condition [Test4]. This is because GPC3_9 and GPC3_3 have the same binding sites in GPC3, which interferes with GPC3 binding when they were used for the ALISA [Test3]. The blanks, control (C1–C5) and Test1–Test3 groups had no effect on signal response (Figure 3A). Increasing GPC3 concentration increased the signal response, and ALISA demonstrated identified <1 ng/mL GPC3 in clinically relevant serum levels. The results of the ALISA followed the linear equation Y=616.92Log10(X)−342.75 (R2=0.9958) (Figure 3B) and corresponded significantly to the binding interaction results (Figure 3C).
3.5 Aptaprobe-based fluorescence imaging: Aptaoptical imaging
GPC3 has three parts: signal peptide for membrane translocation (SP), cysteine-rich domain (CRD), and glycosylphosphatidylinositol (GPI). The Notum cleavage site breaks the bond between amino acids 562 and 563 [47], causing proteins other than GPI to be secreted into the serum (Figure S2). FAM-labeled GPC3_3 and GPC3_2 aptaprobes successfully achieved targeted imaging in HepG2 cells, which are known to overexpress GPC3 protein on their cytoplasmic membrane (Figure 4). However, HepG2 cells were incubated with the FAM-labeled GPC3_9 aptaprobe and confocal imaging results demonstrated that the fluorescence signal was reduced (Figure S6). This may be due to the non-binding of the GPC3_9 aptaprobes that did not bind to their potential binding site (GPI site, Table S1) by the cell membrane-bound region. However, FAM-labeled GPC3_3 and GPC3_2 aptamers successfully achieved targeted imaging in HepG2. FACS experiments supporting the results of confocal imaging, blocked GPC3_9 aptaprobe binding. This is consistent with the results of previous confocal imaging experiments (Figure 4). As mentioned earlier, GPC3 is a heparan sulfate proteoglycan, and the relationship between protein glycosylation alteration and liver diseases has been reported [54]. Our GPC3_3 aptaprobes may be used to monitor glycosylation alterations during HCC progression.
The aptaprobe-binding capacity can be maintained within the actual cellular environment. We used two pancreatic cancer cells, PANC-1 and MIA-PaCa, as negative controls for aptamer binding, as these cells do not express GPC3. To confirm GPC3 expression, the mRNA expression levels were determined preferentially. GPC3 gene expression was observed in only two HCC cell lines, HepG2 and Hep3B, but not in the pancreatic cancer cells (Figure S7). To verify extracellular GPC3 targeting, GPC3-positive and -negative cell lines were incubated at 37°C for 30 min with FAM-labeled GPC3_3 aptaprobes. An intense fluorescent signal was observed in GPC3-positive cells, indicating that GPC3_3 aptamers specifically targeted the surface of the liver cancer cell membrane. FACS experiment using FAM-modified aptamers confirmed the result of confocal imaging. A shift ratio of <1% was observed in the GPC3-negative cell lines (HeLa, PANC-1, and MIA-PaCa), while a shift ratio of approximately 99% was observed in the GPC3 positive cell lines (HepG2 and Hep3B) (Figure 5). These results indicate that the aptaprobe can distinguish the GPC-3 positive cancer cell lines from the GPC-3 negative-cell lines. Chemical modifications (Cy3 and FAM) are beneficial for applications in HCC imaging. In the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide salt (MTT) assay, HepG2 cells incubated with GPC3_3 aptaprobes retained their high survival properties, indicating the low cytotoxicity of GPC3_3 aptaprobes (Figure S8A). Corresponding to confocal imaging, the HepG2 cells did not show any nuclear or cytoplasmic morphological changes. Additionally, the intracellular metabolic potential and mitochondrial respiration of the aptaprobe were confirmed. Since at least 90% cellular energy is produced by mitochondria, we investigated whether the aptaprobe affects mitochondrial structure, which reflects a disruption in cellular bioenergetics. We utilized the Seahorse XF analyzer mito stress test, which allows the real-time measurement of both extracellular acidification rate (ECAR), a measure of glycolysis, and oxygen consumption rate (OCR), a measure of oxidative phosphorylation [55, 56]. A 1 h treatment with 0.15 to 10 µM GPC3_3 aptaprobe slightly increased baseline OCR than that in the control, and no difference in baseline ECAR (Figure S8B and S8C). The 10 µM aptaprobe concentration did not affect toxicity, cytoplasmic morphology, or mitochondrial respiration. Taken together, the aptaprobe has low (or no) cytotoxicity and can be used in various in vitro and in vivo bio-imaging experiments.
Finally, we tested the ability of the GPC3_3 aptaprobe to bind GPC3-expressing tumor in HCC xenograft models in vivo. We established HCC xenograft tumors grown subcutaneously in BALB/c nude mice (Figure S9). To avoid the overlay between the tumor and other tissues in fluorescence imaging, an approximately 169.76 mm3 tumor was successfully induced between the right flank and thigh (data not shown). To confirm that the GPC3_3 aptaprobe targets the target GPC3 protein in vivo, the GPC3_3 aptaprobe and the same size of the random sequence was labeled at the 5′ end with Cy5.5 fluorescence and intravenously injected to measure the fluorescence distribution. The retention of the Cy5.5 signal was assessed using the IVIS system. We examined GPC3_3 aptaprobe binding to induce tumor intravenous injection. Figure 6A shows the fluorescence image of Balb/c nude mice bearing subcutaneous HCC tumors after intravenous injection of ControlCy5.5 and GPC3_3Cy5.5. All injected formulations were imaged immediately after injection, indicating the whole body through blood circulation after intravenous injection. The GPC3_3Cy5.5 signal shifted to the right flank and thigh, where the induced tumor was localized. The fluorescence signal intensity of normalized tumor area with GPC3_3Cy5.5 was significantly higher than that of ControlCy5.5 signal (1.4×109 vs 2.89×109 counts; P=0.002; Figure 6B). Then, ex vivo fluorescence imaging was conducted to examine the biodistribution profile 1 h post injection. We found that the fluorescence signal was also observed in the tumor tissue of the GPC3_3Cy5.5 injected mice group, but no ControlCy5.5 signal was not observed in the control group (Figure 6C). This indicates that the GPC3_3 aptaprobe, which has high affinity and specificity to the GPC3 protein, is quickly distributed and targeted to HCC tumor tissue. This verifies that the aptaprobe can have high specificity and discriminate the target protein in complex biological samples. In addition to tumor sites, the GPC3_3 aptaprobe (as indicated by fluorescence) was mainly distributed in the liver and kidney (Figure S10), consistent with the biodistribution data of ControlCy5.5 [57–59]. It should be emphasized that the relatively weak fluorescence of GPC3_3Cy5.5 aptaprobe in the tumor fluorescence image might be due to fluorescence quenching by in vivo biomolecules because the tumors were isolated > 1 h after the injection of GPC3_3Cy5.5 aptaprobe. After intravenous injection of GPC3_3Cy5.5 aptaprobe, the signal in the tumor region dramatically increased, leading to the visualization of the tumor at 1 h post injection (Figure 6A). The fluorescence signal of the aptaprobe in the tumor region appears faster than that in the in vivo imaging experiment using other materials; such nanoparticles appeared within 72 h, and biological molecules (peptides and aptamers) were observed within 1 to 2 h (Table S2). This should be mainly attributed to the fact that fluorescence signal imaging usually relies on in vivo experimental conditions. These data highlight that aptaprobe-based fluorescence imaging allows fast detection of tumors with high specificity and targeting over other materials.