Faisal Hossainab and
Michael J. Serpe
*a
aDepartment of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada. E-mail: serpe@ualberta.ca
bDepartment of Chemistry, Faculty of Science, University of Chittagong, Chattogram 4331, Bangladesh
First published on 5th February 2025
A method was developed for quantifying single-stranded DNA (ssDNA) through enzymatic digestion and using commercially available glucose test strips. The process involves the initial digestion of ssDNA using a combination of exonuclease 1 and alkaline phosphatase enzymes, leading to the liberation of phosphates from the ssDNA backbone as free orthophosphate. Subsequently, the orthophosphates react with maltose and maltose phosphorylase, producing equivalent amounts of glucose to orthophosphate. The resulting glucose, which can be related to the ssDNA concentration, can be measured amperometrically with an off-the-shelf glucose test strip connected to a mini potentiostat. This method offers versatility, allowing the determination of ssDNA, regardless of nucleotide-count or sequence, with increased sensitivity as the number of nucleotides (NT) in the DNA increases. The method exhibits a limit of detection of 780 nM for 22-NT, 527 nM for 53-NT, 422 nM for 75-NT, and 329 nM for 87-NT ssDNA, and a linear range of 0–2 μM. To selectively quantify a specific ssDNA target, a magnetic microparticle-based isolation step was incorporated, demonstrating high selectivity for quantifying a particular ssDNA target from a mixture. The method holds potential for label-free quantification of ssDNA that can have an impact in myriad fields.
Numerous analytical methods are available for quantification of DNA in solution, and the choice depends on the purpose of the analysis, required sensitivity, and resource availability. The most basic method that can be used for DNA quantification is UV-visible spectroscopy which utilizes the absorbance of DNA at 260 and 280 nm.9 However, this method is susceptible to errors due to potential interfering species in the solution that absorb at the same wavelengths and suffer from low sensitivity due to the measurement method itself. Quantitative polymerase chain reaction (qPCR) is a powerful tool for DNA quantification, relying on enzymatic DNA amplification and subsequent fluorescence detection.10 While qPCR offers high sensitivity, it necessitates the use of expensive equipment and qualified personnel for experimental execution. To overcome these problems, researchers have developed less resource-intensive and quantitative biosensors to quantify DNA in solution.5,7,11,12 For example, Kirimli et al. developed a lead magnesium niobate–lead titanate piezoelectric sensor to quantify single-stranded DNA (ssDNA) by utilizing a complementary probe DNA (pDNA) immobilized on a sensor plate.13 Upon addition of a sample, the target DNA (tDNA) binds to the immobilized probe DNA that causes a frequency change of the sensor that can be monitored in real-time. In another example, Can et al. developed an electrochemical hybridization biosensor to quantify target tDNA found in E. coli K12 using immobilized pDNA on a gold electrode and finally measuring amperometric signal changes upon target addition.14 Although most of the recent methods provide very good sensitivity for DNA detection, they mostly require sophisticated sensor fabrication processes, which are sometimes expensive and difficult to use in a resource-limited environment.15,16 Conventional colorimetric assays need laboratory-based spectrophotometers to quantify analytes and some field applications are semi-quantitative in nature and any background solution color can interfere with the analysis.
In this study, we present a method for quantifying ssDNA that eliminates the need for complicated sensor design and fabrication processes while not requiring sophisticated laboratory equipment. Instead, an off-the-shelf glucose test strip can be employed for ssDNA quantification. Specifically, the ssDNA that is being detected is digested using a combination of exonuclease 1 (EXO 1) and alkaline phosphatase (AP) to release most of the phosphates from its backbone. EXO 1 degrades the phosphodiester bonds in the DNA backbone, while AP cleaves any remaining phosphomonoester bonds, resulting in the generation of free orthophosphate in the solution. Previously, we developed an orthophosphate sensor utilizing an enzymatic reaction that involves maltose phosphorylase and maltose and the reaction generates glucose and beta-D-glucose 1-phosphate as the final products.17 The amount of glucose has been found to be equivalent to the concentration of orthophosphate and the resultant glucose was quantified utilizing off-the-shelf glucose test strips. Based on the previously optimized orthophosphate sensor, in this study, we utilized the sensor to generate glucose from the orthophosphate obtained from DNA digestion. We found that the concentration of the produced glucose is directly proportional to the initial concentration of ssDNA in the solution and its NT count. Subsequently, the produced glucose can be measured amperometrically using a commercially available Accu-Chek Guide glucose test strip from Roche Diabetic Care, connected to a handheld mini-potentiostat (BDTminiSTAT100, from Biodevice Technology, Ltd, Ishikawa, Japan). Commercially available electrochemical glucose test strips use glucose oxidase or glucose dehydrogenase enzyme to selectively oxidize glucose from blood samples. This oxidation reaction generates 2 electrons for each glucose molecule that is carried by the redox mediators to the working electrode of the test strip.18 The voltage applied across the working and counter electrode (450 mV) is large enough so that the reaction at the electrode is instantaneous and the generated amperometric signal is limited by diffusion process. The amperometric current generated by the test strip follows the Cottrell equation:18
i = (nFAD1/2C)/(πt)1/2 | (1) |
Type | Nucleotide count | Nucleotide sequence | Melting temperature, Tm (°C) in 50 mM NaCl |
---|---|---|---|
— | 22 | 5′-GCC GTT TGG GCC CAA GTT CGG C-3′ | 66.2 |
Model target ssDNA | 53 | 5′-GCA GAA CTT ACG ACC CAG GGG GGT GGA CAG GCG GGG GTT AGG GGG GTC GTA AG-3′ | 75.1 |
— | 75 | 5′-ATA CGA GCT TGT TCA ATA CGA AGG GAT GCC GTT TGG GCC CAA GTT CGG CAT AGT GTG GTG ATA GTA AGA GCA ATC-3′ | 71.2 |
— | 87 | 5′-TAT CGC ATA GCG ATA CGA GCT TGT TCA ATA CGA AGG GAT GCC GTT TGG GCC CAA GTT CGG CAT AGT GTG GTG ATA GTA AGA GCA ATC-3′ | 71.6 |
Probe DNA | 14 | 5′-GTC GTA AGT TCT GC/3 Bio-TEG/-3′ | 42.3 |
The second step involved the determination of tDNA through enzyme treatment and subsequent detection using a glucose test strip. To accomplish this, 30 μL of the isolated tDNA solution in DI water (concentrations 0–4 μM) was reacted with an 8 μL digestion enzyme mix at 37 °C for 1 h. The digestion enzyme mix contains 4 μL of 20000 U mL−1 EXO 1 enzyme, 2 μL of 10× nuclease buffer (10 mM 1,3-bis(tris(hydroxymethyl)methylamino)propane-HCl (Bis-Tris-Propane-HCl)), 10 mM MgCl2, 1 mM dithiothreitol, pH 7.0, and 2 μL of 1000 U mL−1 AP enzyme. During this reaction, nearly all of the single-stranded DNA breaks down into its nucleotides and generates free orthophosphate. The reaction product was then treated with 11 μL of a second reagent mixture (4 μL of 200 U mL−1 MP, 2 μL of 100 mM maltose, and 4 μL of 1 M citrate buffer, pH 4.6) at 30 °C for 90 min. In this step, orthophosphate in the solution reacts with maltose in the presence of MP to produce glucose. The solution containing glucose was introduced into the Accu-Chek Guide test strip connected to a potentiostat, and the amperometric signal was measured at 450 mV applied potential. The obtained current from 3 to 20 s was averaged. The ratio (relative current) of the averaged current of the sample to the corresponding averaged current obtained from the reagent blank was calculated and plotted against the sample concentration to generate calibration curves. The experimental steps are illustrated in Scheme 1.
![]() | ||
Scheme 1 Experimental steps for determination of a target ssDNA using magnetic microparticle isolation and quantification using an off-the-shelf glucose test strip. |
To investigate the impact of each digestion enzyme on a ssDNA, EXO 1 and AP were individually reacted with 4 μM 87-NT ssDNA following the procedure outlined in the Experimental section. Later, a mixture of both enzymes was reacted with 4 μM 87-NT ssDNA and detected with the Accu-Chek Guide test strip. When the mixture of EXO 1 (final concentration ∼1053 U mL−1) and AP (final concentration ∼106 U mL−1) was reacted with the 87-NT ssDNA, and subsequently reacted with MP and maltose, a significant increase in the amperometric signal was observed, indicating substantial breakdown of the ssDNA and the generation of a notable amount of orthophosphate/glucose (Fig. 1A). In contrast, individually reacting 4 μM 87-NT ssDNA with either AP (final concentration ∼106 U mL−1) or EXO 1 (final concentration ∼1053 U mL−1) resulted in very low relative currents, close to 1 (Fig. 1B), indicating that the cumulative action of EXO 1 and AP is essential for optimal digestion and the generation of free orthophosphate from the ssDNA backbone.
In further experiments, we investigated how AP concentration impacted the performance of the technique, i.e., we compared the signals obtained from ∼53 U mL−1 (AP*) to that obtained from ∼106 U mL−1 using the same amount of EXO 1 in the mixture (final concentration ∼1053 U mL−1). The relative currents for both AP concentrations were comparable (Fig. 1B), leading to the selection of ∼53 U mL−1 AP for subsequent experiments. To determine the optimal duration for the digestion reaction, 4 μM 87-NT ssDNA was exposed to the digestion enzyme mix for 0 to 120 min, followed by a subsequent reaction with maltose and MP for 120 min, following the procedure outlined in the Experimental section. The relative current exhibited an increase as reaction time increased, with the rise becoming less pronounced after 40 min of the digestion reaction (Fig. 1C). To ensure effective digestion with minimal time, a 60 min digestion reaction was selected. A comparison between the Accu-Chek Guide test strip, a newer version from Roche Diabetic Care, and Accu-Chek Aviva test strip used in our earlier reported methods17,25 revealed that the Accu-Chek Guide test strip provided significantly higher relative current values (Fig. 1C). The observed difference in current between the two test strips likely arises from a combination of different factors, including variations in enzyme type, mediator properties, and potentially electrode design or area. Our focus was on selecting a sensor based on experimental current outputs while maintaining a constant applied potential and without altering the test strips themselves. When selecting glucose test strips for this method, the primary considerations were selectivity and reproducibility. Since the method uses maltose as a substrate to generate glucose, it was crucial that the test strips demonstrate high selectivity for glucose over maltose. Both the Accu-Chek Aviva and Accu-Chek Guide strips use enzymes highly specific to glucose and exhibit a minimal response to maltose, enabling sensitive glucose detection in our approach. Of the two, the Accu-Chek Guide, being a newer version currently available on the market, provided higher sensitivity and was therefore chosen. The test strips were used in their original, unaltered form to ensure a user-friendly application. This method was designed so that users can purchase commercially available test strips and, without any modifications, perform DNA detection. By simply connecting the strips to a handheld mini-potentiostat, the process becomes accessible and straightforward. To assess the effect of reaction time on glucose generation from orthophosphate and maltose in the presence of MP, 4 μM 87-NT ssDNA was subjected to the digestion enzyme mix for 1 h, followed by reactions with maltose and MP for durations ranging from 0 to 150 min. The relative currents increased with reaction time and reached a plateau after 90 min (Fig. 1D). Hence, a 90 min duration for the glucose generation reaction was selected for further experiments.
To evaluate the method's performance across ssDNAs with various NT counts, solutions of ssDNA with NTs ranging from 22 to 75-NT were prepared at concentrations varying from 0 to 2 μM. Each ssDNA solution underwent reactions with the digestion enzyme mix, followed by reactions with maltose and MP, utilizing the same procedure described in the Experimental section. The absolute current obtained was found to be dependent on the temperature, humidity, and potentially other factors and it was found to be different from day to day measurement. On the other hand, if the current is measured against the corresponding reagent blank and a ratio is taken it improves the day-to-day reproducibility of the signal and the obtained calibration curve is reproducible across experiments. The relative currents were calculated by taking the ratio of the averaged current (nA) from 3 s to 20 s for each ssDNA sample and the averaged current obtained from the reagent blank (0 μM ssDNA). The relative current increased proportionally with the concentration of the ssDNA samples. Notably, for the same concentration of ssDNAs, the signal increased as the length of the DNA (and the number of NTs present in each ssDNA) increased, as can be seen in Fig. 2B. This can be attributed to the increased availability of intrinsic orthophosphate in long/high NT count ssDNA (compared to short/low NT count ssDNA), leading to the generation of a higher amount of glucose and consequently yielding a larger signal in the amperometric measurements. The linear regression equations and limit of detection, calculated by taking 3 times the standard deviation of the blank divided by the slope of the calibration curve, for 0–2 μM linear range for all of the ssDNA analyzed are shown in Table 2. Blank data point (0, 1) was not included in the linear regression calculation since the method uses relative current (to the blank) for each concentration data for generating the calibration curves. Following the trend of improved LOD with longer DNA, it stands to reason that extremely long DNA (e.g., genomic DNA) would yield very low LODs. The LOD in ng μL−1 unit, which is independent of the length of the ssDNA, is also reported in Table 2 for comparison between other DNA detection methods.
Nucleotide count | Regression equation | Regression coefficient, R2 | Limit of detection, LOD (nM) | Limit of detection, LOD (ng μL−1) |
---|---|---|---|---|
22 | y = (0.23 ± 0.02)x + (1.05 ± 0.02) | 0.9599 | 780 | 5.26 |
53 | y = (0.34 ± 0.04)x + (1.08 ± 0.04) | 0.9966 | 527 | 8.79 |
75 | y = (0.43 ± 0.01)x + (1.08 ± 0.01) | 0.9911 | 422 | 9.84 |
87 | y = (0.55 ± 0.03)x + (1.11 ± 0.03) | 0.9899 | 329 | 8.89 |
We compared the performance and simplicity of the method with the recently published method in Table S1.† Although other DNA quantification method provides higher sensitivity, they require a complex sensor fabrication process which is somewhat difficult in a resource-limited scenario. On the other hand, this method provides a simple and label-free detection technique that does not require any sensor fabrication process and a user can simply buy off-the-shelf glucose test strips for sensitive measurement. Moreover, we anticipate getting significantly higher sensitivity when measuring genomic DNA/long DNA targets improving the performance of this method to several folds.
To determine the extent of the ssDNA digestion using this method, the relative currents obtained from ssDNA were compared with the relative current obtained from standard orthophosphate solutions. To facilitate this comparison, 0–350 μM standard orthophosphate solutions were subjected to all the steps that ssDNA would go through to generate the free NTs, and finally measured using the Accu-Chek Guide test strip. By employing the equation representing the best fit line in Fig. 2C, we can determine the theoretical relative current for the complete digestion of a ssDNA. That is, for a given concentration of ssDNA, with a known NT count, the theoretical relative current (“y”) can be calculated for 100% efficient digestion and phosphate generation (“x”). For example, after 100% efficient/complete digestion of 2 μM 87-NT ssDNA, the concentration of generated orthophosphate (“x”) would be 174 μM. By plugging this ideal orthophosphate concentration into the regression equation from Fig. 2C, the calculated relative current obtained is expected to be 2.0 ± 0.1. When the experiment was carried out, experimental relative current observed for 2 μM 87-NT was 2.18 ± 0.06. Using these values, we can calculate the % digestion using eqn (2), which is calculated as (109 ± 6)%. The calculated % digestion ranged from (97 ± 7) to (109 ± 6) across all ssDNA sequences investigated here, signifying this method can digest ssDNA with nearly 100% efficiency, allowing for maximal sensitivity.
% Digestion = (Experimental relative current/Calculated relative current) × 100 | (2) |
ssDNA nucleotide count | Concentration (1 μM) | Concentration (2 μM) | ||||
---|---|---|---|---|---|---|
Experimental relative current | Calculated relative current | % digestion | Experimental relative current | Calculated relative current | % digestion | |
22 | 1.30 ± 0.06 | 1.33 ± 0.08 | 98 ± 7 | 1.5 ± 0.1 | 1.44 ± 0.08 | 104 ± 9 |
53 | 1.48 ± 0.09 | 1.48 ± 0.08 | 100 ± 8 | 1.7 ± 0.1 | 1.73 ± 0.09 | 98 ± 8 |
75 | 1.53 ± 0.08 | 1.58 ± 0.09 | 97 ± 7 | 1.92 ± 0.09 | 1.9 ± 0.1 | 101 ± 7 |
87 | 1.686 ± 0.001 | 1.64 ± 0.09 | 103 ± 6 | 2.18 ± 0.06 | 2.0 ± 0.1 | 109 ± 6 |
To determine the efficiency of hybridization of pDNA and 53-NT tDNA, 4 μM of pDNA was reacted with 4 μM of tDNA for 1 h at room temperature with gentle shaking without any magnetic microparticles. The final solution was analyzed in 14% non-denaturing PAGE as described in the Experimental section with two control solutions: ‘only pDNA’ (4 μM) and ‘only tDNA’ (4 μM). As expected pDNA band was further down in the gel while tDNA band was slower due to its larger size and shape as shown in Fig. 3A. In a non-denaturing PAGE, a smaller size and shape DNA is able to move faster than a large DNA with more complex secondary structures.
The hybridized dsDNA band migrated slower than the tDNA band, signifying hybridization between the pDNA and tDNA. Importantly, the pDNA band was not observed for the mixture DNA solution, signifying almost no free pDNA was left in the solution. Hence, the pDNA chosen can efficiently bind to the tDNA and is able to capture and isolate the tDNA from the mixture for downstream analysis.
Next, solutions of 0–4 μM tDNA were prepared and individually exposed to pDNA-modified magnetic microparticles following the procedure outlined in the Experimental section. After isolation of the particles, they were redispersed in water and heated at 100 °C for 10 min to “melt” the hybridized double stranded DNA and release the tDNA back into solution. The tDNA solution was then reacted with the digestion enzyme mix and reagent to generate glucose. The relative current of the 53-NT tDNA increased linearly with increasing concentration, which allowed a calibration curve to be constructed, and the concentration of a tDNA to be determined. The calibration curve obtained is shown in Fig. 3B. To check for any non-specific binding on the magnetic microparticles, 75-NT ssDNA was exposed to the same pDNA-modified magnetic microparticles specific for the 53-NT tDNA. As expected, the relative current obtained from the 75-NT ssDNA was very low and did not increase with increasing concentration (Fig. 3B). This implies that the method is selective for specific tDNA sequences (in this case, 53-NT tDNA) with minimum non-specific binding. Then both 53 and 75-NT ssDNA were mixed together in a 1:
1 ratio from concentration 0.25 μM to 2 μM and measured using the procedure described in the Experimental section. To check whether a mismatched DNA can interfere with the analysis of tDNA, we checked for any significant differences in each data points for a certain concentrations of tDNA to the mixture (53 and 75-NT). Paired t-tests were done on relative currents of the isolated 53-NT ssDNA and the mixture of 53-NT tDNA and 75-NT ssDNA for concentrations from 0.25 to 2 μM. The p values were: 0.43 for concentration 0.25 μM, 0.25 for concentration 0.5 μM, 0.16 for concentration 1 μM, and 0.14 for concentration 2 μM. Hence, the signals from the mixture were not significantly different (p > 0.01) from the signal that was obtained from the only 53-NT tDNA (Fig. 3B). Hence, this method can selectively isolate a certain tDNA from a mixture and allows for accurate determination.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sd00318g |
This journal is © The Royal Society of Chemistry 2025 |