Payel Sena and
Manisha Gupta*b
aDepartment of Chemical and Materials Engineering, University of Alberta, Edmonton, Canada. E-mail: payel@ualberta.ca
bDepartment of Electrical and Computer Engineering, University of Alberta, Edmonton, Canada. E-mail: mgupta1@ualberta.ca
First published on 3rd February 2021
Single nucleotide detection is important for early detection of diseases and for DNA sequencing. Monolayer (ML) MoS2 nanopores have been used to identify and distinguish single nucleotides with good signal-to-noise ratio in the recent past. Here, we use a bilayer (BL) MoS2 nanopore (∼1.3 nm thick) to detect distinct single nucleotides with high spatial resolution and longer dwell time. In this study, the performance of similar sized (<3 nm) ML and BL MoS2 nanopores for detection of a single nucleotide has been compared. Both single nucleotide and single stranded DNA translocations through them are studied. For single nucleotide detection, we observe that BL MoS2 nanopores demonstrate twice the dwell time as compared to ML MoS2 nanopores with 95% confidence. Single nucleotide detection rate for BL MoS2 nanopores (50–60 nucleotides per s) is five-fold higher as compared to ML MoS2 nanopores (10–15 nucleotides per s) in 10 pM analyte concentration. For single stranded DNA, we observe 89% (for 60 DNA molecules detected) single nucleotide detection efficiency with BL MoS2 nanopores as compared to 85% for ML MoS2. The DNA sequencing efficiency through BL MoS2 nanopores is also found to be 8–10% better than through ML MoS2 nanopores, irrespective of DNA sequencing orientation. Thus, owing to improved analyte/nanopore charge interaction BL MoS2 nanopores can be used for single nucleotide detection with high resolution due to longer dwell time, detection rate and efficiency. This study demonstrates the improved ability of BL MoS2 nanopores in sequencing DNA with 8–10% higher efficiency, two-times temporally resolved single-nucleotide current signatures and five-times higher detection rate, compared to ML MoS2 nanopores.
van der Waals force separate the layers of 2D materials like MoS2. Thus, the number of layers of MoS2 separated can be controlled, which allows tunability of the membrane thickness. Single nucleotide and polynucleotide translocation through monolayer MoS2 nanopores has been previously demonstrated experimentally.13 However, MoS2 nanopores also suffer from high translocation velocity similar to other solid-state nanopores. Although different techniques have been employed to slow down DNA translocation through solid state nanopore,28–31 the control of translocation speed is still challenging. The thickness of monolayer (ML) MoS2 (0.65 nm) is less than a single nucleotide (1.6–1.8 nm). Thus, the charge interaction along DNA/nanopore interface for ML MoS2 gets compromised making DNA translocation fast. Molecular translocation through ML MoS2 nanopore has been slowed down by using transmembrane viscosity gradient for producing temporally resolved current blockades capable of distinguishing different types of DNA nucleotides.13
In this work, we experimentally investigated the role of 2D material membrane properties on the translocation speed of DNA. It is found that the behavior of monolayer (ML) and bilayer (BL) MoS2 is different under a vertical electrical field applied using the cis and trans electrodes.32 The strong interlayer coupling in the BL MoS2 leads to a gradation in potential throughout the thickness as each layer will experience a different electric potential.32 Due to this potential gradient the phosphate groups are trapped and detrapped by the Molybdenum atoms in the membrane. An additional pull is created because of the increased negative (analyte DNA)–positive (Mo atoms lining the nanopore surface) charge interaction. This can help immobilize molecules at the nanopore for longer time producing highly resolved signals. Molecular dynamics simulation study conducted previously concluded that BL MoS2 can slow down peptide translocation producing stepwise distinguishable current blockades for single nucleotide sensing.22 Thus, we have conducted detailed single and poly nucleotide translocation experimental studies through ML and BL MoS2 nanopores. The translocations were conducted using picomolar concentration and the dwell times were then extracted from these measurements for further analysis.
300 mM filtered KCl solution at pH = 7 was prepared by adding 3 mM Tris–HCl buffer. The solution was then degassed in vacuum for 90 minutes and used for the experiments after attaining room temperature. 0.5 mM solutions of single nucleotides (deoxyadenosine triphosphate: dA, deoxythymidine triphosphate: dT, deoxycytidine triphosphate: dC and deoxyguanosine triphosphate: dG) were procured from Integrated DNA Technologies for nanopore translocation and dwell-time analysis thereof. 5 mM solution of single stranded (ss) customized DNA oligos (3′-ATCGATCGATCGATCGATCGATCGATCG-5′) were procured from Integrated DNA Technologies to evaluate the efficiency of ML and BL MoS2 nanopores for detection of single nucleotides from a DNA strand. Analyte solutions at 10 pM concentration were then prepared by dilution of the purchased single nucleotide and ssDNA solutions.
Custom-made Teflon half-cells were used to contain the electrolyte (KCl). As Teflon or poly-tetra-fluoro-ethylene are both hydrophobic polymers, they help in maintaining a clean environment. They also have high chemical resistance which ensures a well-insulated sensing environment. The half cells were sealed with polydimethylsiloxane (PDMS) to ensure that the electrolytic exchange between the chambers is only via the nanopore. PDMS, also reduces capacitive noise which helps in obtaining improved signals. The half-cells were thoroughly cleaned by ultra-sonication before each experiment. The nanopores/membranes were dipped in acetone and kept in vacuum for 30 min. The nanopore was cleaned by gradually pulling the acetone through it. The cleaning steps include three acetone-cleaning repeats followed by four IPA washes in vacuum (20 min each) to remove acetone residue. Fresh solvent is used for every cycle of acetone and IPA treatment. The membrane with the embedded nanopore was then mounted between the half cells. Two Ag/AgCl electrodes were cleaned using ethanol and DI water and re-chlorinated by treating them with bleach for an hour. The electrodes were then dipped in each half-cell to apply a voltage across the nanopore bearing membrane. The entire cell-assembly was enclosed in a Faraday cage (Warner instruments) to prevent electrical interference. Filtered and buffered 300 mM KCl solution (previously prepared) was introduced in both chambers and the ionic current at varying voltage (−200 mV to 200 mV) was recorded to evaluate the channel conductance. The nanopore was electrically conditioned at a constant 100 mV bias to obtain a stable baseline. The analytes were then added to the trans chamber and pulled electrophoretically to the cis chamber through the nanopore under a constant bias of 200 mV. Single channel blockade recordings obtained from translocating nucleotides and DNA were then filtered (20 kHz 8-pole Bessel filter), amplified (Axon MultiClamp 700B) and finally digitized (Axon Digidata 1550B) for further analysis. ESI Fig. S1b† presents schematic of the sensing assembly.
A customized MATLAB algorithm was used to identify translocation peaks. The raw traces were first loaded as an .abf file. The peaks were identified as change points and the blockade current was calculated by subtracting the mean blockade current from the mean baseline current. The mean was calculated for the 100 samples before and after the peak to account for shift in baseline during the translocation. The dwell time was calculated from the difference of the initial and final changepoint indices multiplied by the sampling frequency. The mean and median of the dwell time were then compared to test whether the distribution is normal or not. We then assumed a null hypothesis that the dwell times for ML and BL MoS2 nanopore is the same and then tested the hypothesis by conducting T-test at 95% confidence on the dwell time data. All statistical analysis was performed for a population of 300 nucleotides. For the mixed DNA, a MATLAB algorithm was built to first truncate the data at each molecular entry and exit. Then the dwell time and blockade current at each individual current drop steps within a molecular translocation event were obtained. The individual nucleotides were then identified by comparing the blockades with that obtained from single nucleotide translocation performed earlier through the same nanopore.
The ionic conductance of the fabricated ML and BL MoS2 nanopores were measured following the procedure described in the methods section to ensure properly cleaned and conducting nanopores. Fig. 2a and b shows MoS2 nanopores with diameters of 2.6 nm for ML and 2.7 nm for BL respectively. Ionic current through each nanopore was measured by applying a trans bias in the range of −200 mV to 200 mV in 300 mM KCl. 4 consecutive ionic conductance measurements were conducted through the nanopore along with intermediate cleaning steps. Here, we have reported the average conductance for each pore. Fig. 2c and d shows the IV plot for both the nanopores. The mean of ionic current recorded at each bias is displayed as scatter plot, with the error bars showing the current variation between the 4 measurements at each voltage. The mean nanopore conductance (calculated as the average of the ratio of ionic current and voltage at each bias value) was found to be 10.5 ± 0.016 nS for 2.6 nm ML MoS2 nanopore and 17.2 ± 0.023 nS for 2.7 nm BL MoS2 nanopore.
To study the effect of ML and BL MoS2 on single nucleotide translocation and sensing, 10 pM concentration of dA, dT, dC and dG were electrophoretically pulled through the above-mentioned ML and BL MoS2 nanopores. For these nucleotide translocations 200 mV bias was applied and the translocation traces were recorded. The ESI Fig. S3a–h† shows the recorded data for single nucleotide (for each of dA, dT, dC and dG) translocation for 5 s and 1.5 s through ML MoS2 and BL MoS2 nanopores respectively. We also observe a detection rate of 10–15 nucleotides per s and 50–60 nucleotides per s for ML and BL MoS2 nanopores respectively (see Fig. S3a–h†). The five-fold higher detection rate obtained for BL nanopores may be due to faster capture of the translocating nucleotides at the nanopore due to stronger electrophoretic pull towards the layered BL MoS2.
Fig. 3a–h presents truncated single nucleotide peaks obtained for ML and BL MoS2 nanopores for a direct comparison of dwell times. It is observed that the dwell times are higher for BL as compared for the ML MoS2 nanopores for all the different nucleotides. Blockade current is plotted as a function of dwell time for 3000 single nucleotide transport events in Fig. 3i. We observe four distinct blockade current regions for the different nucleotides. Thus, we can conclude that both ML and BL MoS2 nanopores are capable of detecting single nucleotides. The blockade current for the nucleotide translocation is plotted as histograms to observe their distribution (ESI Fig. S4†). We observe normal distribution for all the nucleotides for both ML and BL nanopores. Thus, the mean blockade current values along with their standard deviations can be obtained.
Table 1 lists the dwell time and blockade current values for all single nucleotides sensed using ML and BL MoS2 nanopores. The blockade current magnitudes measured are found to be the dA, dT, dC and dC in decreasing order which agrees with previous study.13 The mean blockade current obtained through BL MoS2 is found to be slightly higher than or comparable to ML MoS2, which may be due to better charge interaction at the van der Waals separated nanopore interface. Based on these numbers we can conclude that the current resolution obtained for ML is preserved in BL MoS2 nanopores. On the other hand, it was observed that the dwell time values of the translocations for BL nanopore are spread over a larger range as compared to the ML nanopore. The events detected through ML nanopores show mean dwell times in range of 0.053–0.061 ms with a 17–20% of deviation. For BL nanopore detection the mean dwell time range is 0.141–0.143 ms with a deviation of 28–30% (see Table 1 and ESI Fig. S5†). The mean and median dwell time values were then evaluated for ML and BL nanopores (see Fig. 3j) and the values were found to be pretty close to each other, supporting an assumption of a normal distribution. Therefore, T-tests were done on the values for a statistical comparison. T-Tests thus performed, indicate that for all the nucleotides, the dwell times for ML and BL nanopores are different at a 95% confidence level (shown in ESI Table T1†). The sample mean of the dwell time for the BL nanopore is approximately twice that of the ML nanopore. A higher dwell time indicates an improvement in sensing.
Analyte | ML MoS2 nanopore | BL MoS2 nanopore | ||
---|---|---|---|---|
Blockade current | Dwell time | Blockade current | Dwell time | |
dG | 0.595 ± 0.057 nA | 0.0529 ± 0.0092 ms | 0.608 ± 0.036 nA | 0.1409 ± 0.0424 ms |
dC | 0.987 ± 0.066 nA | 0.0577 ± 0.0102 ms | 1.256 ± 0.048 nA | 0.1430 ± 0.0424 ms |
dT | 1.454 ± 0.076 nA | 0.0610 ± 0.0110 ms | 1.809 ± 0.07 nA | 0.1422 ± 0.0405 ms |
dA | 2.245 ± 0.079 nA | 0.0615 ± 0.0122 ms | 2.434 ± 0.06 nA | 0.1412 ± 0.0401 ms |
To demonstrate the single nucleotide detection capability of ML and BL MoS2 nanopores mixed nucleotide translocations through ML and BL MoS2 nanopores were conducted. Fig. S6a and b† show translocation traces for 5 s and 1.5 s through ML and BL MoS2 nanopores respectively. The nucleotide distinction is made based on the current blockades obtained for individual nucleotides sensed separately. The dwell time range for all nucleotides (as observed in the scatter plot in Fig. S6c†) is found to be the same as that when detected separately through ML and BL MoS2 nanopores.
To demonstrate that the BL MoS2 improves the dwell time of the translocation and can detect single nucleotides, we used 10 pM single stranded (ss) DNA oligos with customized mixed sequence. The translocations were conducted using ML and BL nanopores. Fig. 4a and b show the translocation trace for a single ssDNA strand through the ML and BL MoS2 nanopore respectively. For both the nanopores distinct molecular current blockades were observed. From careful inspection of each molecular blockade, it can be clearly understood that the discrete blockade steps correspond to the constituent nucleotides. All the translocations for the ML and BL were analysed using the code developed and dwell time as a function of blockade current for 3 s translocations are plotted in Fig. 4c and d. For BL MoS2 nanopore, we observe four different current levels similar to observed for each nucleotide with dwell times ranging from 0.02–0.6 ms with the highest density obtained around 0.08–0.3 ms which is similar to what we obtained for single nucleotide translocations. On the other hand, for ML MoS2 nanopore, four different nucleotides are also detected. However, the time resolution obtained for dC and dG are particularly low (0.004–0.04 ms for dG and 0.004–0.06 ms for dG) compared to dA and dT, the latter nucleotides showing dwell times in the range of 0.01–0.07 ms, with the mean ∼0.05 ms, which is close to that obtained for single nucleotide triphosphates.
Based on the different dwell times and blockade current four different nucleotides can be distinctly identified. The blockade levels are then matched with the single nucleotide data recorded before to obtain and relate the DNA sequence. Fig. S7a and b† shows the sequence of nucleotides detected and left undetected by ML and BL MoS2 nanopore sensing in a color-coded format, thus representing the single nucleotide efficiency of detection in each case. ML MoS2 nanopore is found to demonstrate a single nucleotide detection efficiency of 84.9% ((total − undetected)/total nucleotides) considering 60 DNA molecules translocation. Similarly, the experimental analysis shows that BL MoS2 nanopore are capable of identifying single nucleotides with about 89.1% efficiency. For both ML and BL nanopores, it is observed that mostly dC and dG (especially at the tailing end or entering end of the DNA) are undetected or detected with poor time resolution. This may be due to the low density of these two nucleotides along with their low blockade current.
In addition, we have extracted the efficiency of ssDNA sequencing for both 3′–5′ and 5′–3′ orientations. BL MoS2 nanopores demonstrate 83 ± 1.07% and 92 ± 0.14% for 5′–3′ and 3′–5′-ssDNA sequencing efficiency respectively. These are 10% better than ML MoS2 for 5′–3′-orientation (73.35 ± 0.55%) and 8.4% better for 3′–5′ sequencing (83.6 ± 0.16%). Fig. S8 and S9† shows sequencing events for 3′–5′ and 5′–3′ ssDNA translocations for ML and BL MoS2 nanopores respectively.
Fig. 4e and f shows the efficiency of detection of individual molecules from ssDNA and ssDNA sequencing as a whole. It is found that the efficiency of BL MoS2 nanopores is better in both cases. The role of BL MoS2 is observed to be most prominent in improving the temporal resolution of the signature for dC and dG. It was observed that the 47/60 ssDNA translocations were in 3′–5′ orientation for ML MoS2 nanopore and 51/60 ssDNA translocations were in 3′–5′ orientation for BL MoS2 nanopore.
In this study, we experimentally explore the unaided capacity of BL MoS2 nanopores in improving real-time sensing resolution, detection rate and efficiency for ssDNA sequencing. We first analysed the blockade current obtained for single nucleotide using ML and BL MoS2 nanopores. Our ML MoS2 nanopore current signatures are similar to that obtained in previous study.13 However, we observed geometry defying blockade signatures for guanine. This may be due to the increased hydrogen bonding potential of guanine or increased physical adsorption energy of guanine on MoS2 which may increase the pore conductance, thus producing a lower overall ionic blockade event.37 For ssDNA sequencing, we also observed lower time resolution of guanine base of the ssDNA compared to other bases, which may indicate that the effect of surrounding A and C over nanopore blockade are dominant over G. Few studies have shown exceptions to the geometry-determined blockade when larger A and G have produced lower blockades than C and T, for both single nucleotide and polynucleotide sensing.13,20,38
We also obtain higher blockade current levels for BL MoS2 nanopore as compared to ML MoS2 nanopore. In general, for 3D bulk materials like SiNx, the conductance decreases with increase in pore/membrane thickness.39 However, in our study we find that 2D MoS2 nanopores do not demonstrate this behaviour. Unlike 3D bulk material nanopores, in 2D material nanopores interlayer charge storage, electro-activity and field confinement between the van der Waals coupled monolayers increases the capacitance of the pore.40 COMSOL Multiphysics simulation study indicates that the ionic conductance of MoS2 nanopores demonstrates the following trend 2-layers > 1-layer > 4-layers > 3-layers > 6-layers > 5-layers and so on above 2 nm diameter nanopore.41 Few studies on 2D materials nanopores like graphene also demonstrate similar behaviour.18,42,43
We further observe nearly double the dwell time for BL MoS2 as compared to ML MoS2 nanopore, with 95% confidence. To observe the effect of the bias on the blockade current we measured the blockade current for three different voltages. Fig. S11† shows blockade current vs. dwell time scatter plot for single nucleotide translocations through ML and BL MoS2 nanopores for varying transmembrane bias values (100 mV, 150 mV and 200 mV). However, significant change in dwell times is not observed for either of ML or BL MoS2 nanopore sensing. This further proves that the events considered are indeed translocation events and strengthens the inference drawn about improved ability of BL MoS2 nanopores in providing better resolution for single nucleotide sensing.
The effect of the temporal resolution improvement offered by BL MoS2 nanopore on sequencing single-stranded DNA oligos was investigated next. From these experiments, we found that the capability of BL MoS2 nanopores can resolve the nucleotide blockade currents and hence can distinguish single distinct nucleotides. Thus, we observed an overall single nucleotide (for 1800 nucleotides in 60 ssDNA strands) detection efficiency of 89% using BL MoS2 nanopore, a 4% improvement compared to that for ML MoS2 nanopores. We also observed an 8–10% improved ssDNA sequencing efficiency for both 3′–5′ and 5′–3′ ssDNA orientations. The interlayer potential gradient for BL MoS2 offers improved capture of analytes inside the nanopore even at low (picomolar) concentration leading to 5× better detection rate. Previous studies have shown that improved nanopore/DNA interaction can result in faster capture of DNA within the nanopore and also cause longer residence of DNA at the nanopore.44 We believe that the improved capture rate and nanopore dwell time is a result of the enhanced DNA/nanopore interaction experienced in BL MoS2 nanopore as compared to ML MoS2 nanopore.
Moreover, the behaviour of even and odd layers of MoS2 is found to differ under an applied vertical electrical field.45 The electric field confinement within the layers is more for odd-numbered layers than even-numbers, which thus improves the ionic conduction for even layers. Moreover, the BL MoS2 thickness is close to nucleotide size which helps maintain a good spatial resolution. However, >2 layers can hamper the spatial resolution. Additionally, our simulation results also show a reduction in blockade signal amplitude above two layers. All things considered; we believe >2-layers MoS2 may not be a better choice for DNA sequencing. The results conclusively suggest that BL MoS2 nanopores can distinguish four type of nucleotides with significantly higher dwell time. Thus, BL MoS2 nanopores can be a suitable choice for more efficient and faster DNA sequencing.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ra10222a |
This journal is © The Royal Society of Chemistry 2021 |