Rapid and sensitive detection of viral particles by coupling redox cycling and electrophoretic enrichment.

The COVID-19 pandemic has highlighted the need for rapid, low-cost, and sensitive virus detection platforms to monitor and mitigate widespread outbreaks. Electrochemical sensors are a viable choice to fill this role but still require improvements to the signal magnitude, especially for early detection and low viral loads. Herein, finite element analysis of a novel biosensor concept for single virion counting using a generator-collector microelectrode design is presented. The proposed design combines a redox-cycling amplified electrochemical current with electrophoresis-driven electrode-particle collision for rapid virus detection. The effects of experimental (e.g. scan rate, collector bias) and geometric factors are studied to optimize the sensor design. Two generator-collector configurations are explored: a ring-disk configuration to analyze sessile droplets and an interdigitated electrode (IDE) design housed in a microchannel. For the ring-disk configuration, we calculate an amplification factor of ∼5 and collector efficiency of ∼0.8 for a generator-collector spacing of 600 nm. For the IDE, the collector efficiency is even larger, approaching unity. The dual-electrode mode is critical for increasing the current and electric field strength. As a result, the current steps upon virus capture are more than an order of magnitude larger compared to single-mode. Additionally, single virus capture times are reduced from over 700 s down to ∼20 s. Overall, the frequency of virus capture and magnitude of the electrochemical current steps depend on the virus properties and electrode configuration, with the IDE capable of single virus detection within seconds owing to better particle confinement in the microchannel.

Introduction

Prompt and accurate detection of viruses is critical for upholding public safety and minimizing widespread outbreaks of infectious diseases.[(Caygill et al., 2010; Fox, 2007)] This is especially relevant in the backdrop of the current COVID-19 pandemic – caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2)[ (Baig, 2021)] – which is responsible for more than 428 million infections and 5.9 million deaths globally at the time of writing this article.[(Johns Hopkins University, 2022)] Although the COVID-19 pandemic has been the most widespread in recent history, outbreaks caused by other viruses, such as influenza, dengue, human immunodeficiency virus (HIV), and Ebola place additional stress on the global healthcare system.[(Kisely et al., 2020)] Widespread testing is a cornerstone of the effort to mitigate the devastating impact of present and future pandemics. Simple-to-use, sensitive, and rapid diagnostics that can be manufactured at relatively low cost and at a large scale are critical to better monitor the spread of infection and minimize the burden on healthcare systems.

Current strategies for detecting viral infections include polymerase chain reaction (PCR), reverse transcription polymerase chain reaction (RT-PCR), [(Bustin, 2000)] enzyme linked immunosorbent assays (ELISA), [(Voller et al., 1976)] isothermal amplification, [(Kievits et al., 1991; Ito et al., 2006)] and immunochromatography, [(Arai et al., 1999)] among others. While effective, many of these techniques require significant overhead and have limited operational capacity outside of clinical settings. Therefore, electrochemical sensors are becoming increasingly used for virus detection, [(Zhao et al., 2021)] especially in point-of-need (PON) applications that require prompt results, minimal sample preparation, and simple-to-use equipment.[(Caygill et al., 2010)] The operating principle typically relies on changes in the current or impedance at the electrode-electrolyte interface from a redox molecule as virus particles attach to the electrode surface.[(Caygill et al., 2010)] To increase the appeal of electrochemical sensors for PON applications, electrode miniaturization is actively being pursued to further reduce sample volume requirements, [(Vishnu et al., 2020)] improve spatial resolution, especially in sensor arrays, [(Wolfrum et al., 2016)] and explore previously inaccessible regions, such as the interior of individual cells.[(Guo et al., 2022)] However, due to the reduced capture area of miniaturized electrochemical sensors, the number of virions reaching the sensor is reduced. In an effort to improve the sensitivity, redox cycling can be employed.[(Wolfrum et al., 2016; Ma et al., 2013; Lee et al., 2018)]

Redox cycling leverages two independently biased working electrodes spaced in close proximity to one another. One electrode is biased above the formal potential (E ) of the redox probe while the other is biased below, allowing the oxidized and reduced species to cycle between the electrodes.[ (Bard and Faulkner, 2001)] This recycling amplifies the current and can lead to increased signal-to-noise ratios. As a result, redox cycling has been applied for the detection of a variety of analytes, including dopamine, [(Kätelhön et al., 2010)] catechol, [(Wolfrum et al., 2008)] ferri/ferrocyanide ([Fe(CN)6]3-/4-), [(Ma et al., 2013; Straver et al., 2012)] paracetamol, [(Goluch et al., 2009; Paixão et al., 2006)] nitrate/nitrite, [(Gross et al., 2015)] and clozapine [(Ben-Yoav et al., 2014; Winkler et al., 2017)] with amplification factors ranging from 1 to 2 to well over 6000.[(Wolfrum et al., 2016; Straver et al., 2012)] For example, Lee et al. recently demonstrated an amperometric redox cycling sensor based on an interdigitated electrode (IDE) design for detecting the hepatitis B virus surface antigen.[(Lee et al., 2018)] The sensor was combined with a commercial ELISA kit and demonstrated more than 10-fold higher sensitivity in dual-electrode mode compared to single-electrode mode, highlighting the potential of this technique for improving virus detection.

The charged surface of virions makes them prime candidates for electrophoretic manipulation. The surface charge is highly dependent on the medium properties, especially pH. Typically, viruses have an isoelectric point below pH 7, making them negatively charged at physiological pH values.[(Michen and Graule, 2009)] As a result, the virus particles are preferentially attracted to the anode where they can attach to the electrode surface. Such particle collisions can alter the electrochemical landscape of the electrode and cause detectable changes in the local capacitance [(Hellberg et al., 2002, 2005)] and charge transfer.[(Rees et al., 2004; Clegg et al., 2006; Xiao and Bard, 2007; Xiao et al., 2008)] The particle collision method has been applied in the detection of bacteria, such as Escherichia coli, [(Lee et al., 2016; Frkonja-Kuczin et al., 2018)] mammalian cells, [(Dick, 2016)] and viruses, [(Dick et al., 2016)] with the ultimate end goal being single entity detection in as little time as possible.[(Goines and Dick, 2020)]

In this work, we present the finite element analysis (FEA) of a novel biosensor concept that combines the advantages of redox cycling and electrophoretic enrichment for rapid detection and counting of single virus particles. The detection scheme is based on a current quenching strategy whereby the electrochemical current from the cyclic reduction and oxidation of a redox couple ([Fe(CN)6]3-/4- as model) is suppressed as virus particles are captured at the electrode surface. Similar “current blocking” strategies have been reported previously using individual microelectrodes.[(Quinn et al., 2004; Boika et al., 2012; Lebègue et al., 2015)] Here, two mushroom-like, generator-collector microelectrode configurations, a ring-disk configuration and IDE design, are demonstrated. The ring-disk configuration could be suitable for application in digital microfluidic platforms, whereas the IDE configuration is more often utilized within microfluidic channels. The mushroom-like design can be achieved based on a scalable, templated electrodeposition strategy and yields tunable geometries that improve the redox cycling characteristics through reduced generator-collector spacing. Such designs have been experimentally reported recently [(Chen et al., 2015; Chandramohan et al., 2016; Cerquido et al., 2018)] and show potential for biosensing applications. We explore various design parameters, including electrode geometry, testing geometry (e.g. microchannel, droplet), and experimental conditions (e.g. voltammetric scan rate, electrolyte concentration), to better understand their effect on the amplification factor (Γ) and collector efficiency (η) and optimize the proposed platform. Depending on the experimental parameters and electrode geometry, single virus collisions can be detected on the order of seconds. The computational results herein provide design rules and assist in guiding the experimental development of biosensors suitable for single particle counting, with potential applications including virus detection, immunoassays, and nanoparticle-based molecular diagnostics.

Methods

Model geometry

The generator-collector ring-disk configuration consists of two concentric 3D microelectrodes, an inner “disk” and outer “ring” for the droplet-based design (Fig. 1 ). For this configuration, a 1 μl hemispherical droplet (unless otherwise specified) containing 1 mm [Fe(CN)6]4- in a supporting KCl electrolyte is used. For the interdigitated electrode (IDE) design (Fig. 1), two electrodes with comb-like fingers are placed in an alternating pattern and contained within a microfluidic channel (not shown). At the top of all electrodes is a mushroom-like cap that could be formed by a template-driven over-electrodeposition.[(Chen et al., 2015)] The width of the electrode base (w ) is 2 μm which is within the resolution of standard photolithography methods. The electrode cap distance (d ) can be controlled through the electrodeposition parameters (e.g. deposition time, temperature, applied potential).[(Cerquido et al., 2018)]

Fig. 1

A schematic representation of the ring disk (cross-section) and IDE sensor geometries along with the current quenching sensing mechanism before and after virus attachment demonstrated in this work. When a virus particle attaches, it blocks the mass transfer that drives the redox cycling. As a result, the magnitude of the electrochemical current is reduced. An example oxidation reaction is included in the inset. The reduced species (R) travels from the bulk solution to the electrode surface where it transfers an electron (e−) to form the oxidized species (O), which then travels back to the bulk solution.

Due to the azimuthal symmetry of the ring-disk configuration, the electrochemical characterization in Section 3.1 is performed using a 2D axisymmetric geometry. For the simulations with virus particles presented in Section 3.2, a 2D cross-sectional geometry is employed. This is because of discrepancies in the electric field at the axis of symmetry when using an axisymmetric geometry, which causes incorrect particle behavior. Note that the current for the ring-disk configuration in Section 3.2 is still computed by integrating over a revolved geometry to ensure correct current values. For the IDE, a 2D cross-sectional geometry is used with an out-of-plane thickness of 0.5 mm. The IDE is 1 mm in length and is comprised of nearly 250 digits. Each electrode digit has a width w  = 2 μm with 2 μm spacing between adjacent digits (a). The height of the microfluidic channel above the IDE is set to 15 μm and 1.07 mm in length in the simulations. Additional details pertaining to the model geometry can be found in Section S.1 of the Supplementary Information.

Electrochemical simulations

Finite element analysis is performed using COMSOL Multiphysics versions 5.5 and 5.6 with the Electrochemistry module. The Tertiary Current Distribution physics node is used. An electroneutrality condition is maintained throughout the electrolyte (i.e. , where z and c are the charge number and concentration of species i). As a result, the concentration of Cl− in the KCl electrolyte is determined by the solver such that this condition is upheld. The K+ concentration in the supporting electrolyte is made up of the KCl concentration, along with contributions from the [Fe(CN)6]4-. Electrochemical analysis is carried out in a generator-collector configuration whereby the two working electrodes are biased independently. For all measurements in Section 3.1, the collector is held at a bias E  = −0.15 V while the generator (E ) is swept across a range of potentials, except when studying the effect of collector bias where other values of E are used (see Fig. S4). The amplification factor (Γ) and collector efficiency (η) are calculated at the end of the voltammetric sweep, where E  = 0.6 V. For the chronoamperometric measurements of Section 3.2, E  = −0.15 V and E  = 0.6 V are used. All simulations are referenced to an Ag/AgCl reference electrode (E  = 0.197 V). The electrolyte potential (φ ) is set equal to E at points at the top-center of the droplet (x = 0), or at the top of the microchannel in the IDE configuration (x = 0). See Fig. S1 for more details.

The mass transport of species i is governed by the Nernst-Planck equation for diffusive and migrative transport,where J is the flux, D is the diffusion coefficient, u is the mobility, F is Faraday's constant, and φ is the electrolyte potential. The mobility u of species i is related to the diffusion coefficient through the Nernst-Einstein relationship,where R is the molar gas constant, and T is the temperature. The one-electron electrochemical reaction at the electrode-electrolyte interface can be written as follows,which has a formal potential E  = 0.173 V vs. Ag/AgCl.[(Rock, 1966)] The resulting electrochemical current from this reaction is modeled using Butler-Volmer kinetics,where i is the exchange current (dependent upon k , the heterogeneous rate constant), α is the anodic transfer coefficient, n is the number of electrons transferred in the reaction, and η is the overpotential. The overpotential is expressed as,where φ is the potential on the electrode boundary, φ is the potential in the electrolyte domain, and E is the equilibrium potential of the redox reaction as determined by the Nernst equation,where c and c represent the concentration of oxidized and reduced species, respectively. The values of electrochemical parameters used are presented in Table 1 .

Table 1

The electrochemical parameters used in simulations [(Ma et al., 2013; Konopka and McDuffie, 1970)].

Parameter	Value
cred*	1 mm
k0	0.1 cm s−1
Dox	7.2 × 10−6 cm2 s−1
Dred	6.4 × 10−6 cm2 s−1
αa	0.50
n	1

Joule heating and heat transfer effects are coupled to the other physics modules through the Heat Transfer in Fluids physics module. The heat generated by current flow through the electrolyte is dissipated through the electrode and bottom substrate boundaries. A heat transfer coefficient (h) of 10 W m−2 K−1 is used.[(Erickson et al., 2003)] An external temperature (T ) of 298 K is used. The remaining boundaries are assumed to be insulating.

Modeling virus particles in an electric field

The virus particles are subject to a random release at t = 0 as the simulation begins. Each virus particle experiences drag (F ), electrophoretic (F ), and dielectrophoretic (F ) forces. The Stokes drag force equation can be expressed as,where μ is the fluid viscosity, r is the particle radius, is the fluid velocity, and is the particle velocity. For the electrolyte solution in the droplet configuration, a 30% wt. water-glycerol mixture is used to limit evaporation during practical experiments. In this case, a viscosity of 2.13 10−3 Pa s is calculated based on the formulation presented in ref.[(Cheng, 2008)] The density of the solution is 1.07 g cm−3, [(Takamura et al., 2012)] and the dielectric constant used is 68.[(Åkerlöf, 1932)] The virus particle density is 1.3 g cm−3, which is similar to previous reports.[(Ermolina et al., 2006)] For IDE design, which is embedded in a microfluidic channel, simulations are done in an aqueous KCl electrolyte without added glycerol. The enclosed nature of the microfluidic channel eliminates the need for the use of glycerol.

The electrophoretic force on a particle is given as,where e is the elementary charge, q is the charge number of the particle, and φ is again the electrolyte potential, which is obtained by coupling with the Electrochemistry module. The dielectrophoretic force can be written as,where ε is the permittivity of free space, ε is the dielectric constant of the medium, K is the Clausius-Mossotti function ( where ε is the dielectric constant of the particle), and is the electric field. The asterisk denotes the complex dielectric constant given as where j, σ, and ω represent the imaginary unit, electrical conductivity, and angular frequency, respectively. The dielectric constant of the droplet solution (glycerol-water mixture) is 68, whereas the dielectric constant for the aqueous electrolyte in the IDE microchannel is obtained from the COMSOL material library.

The virus particles are modeled as an inner sphere with an encapsulating outer layer with dielectric constants of ε = 75 and ε  = 7.5, respectively, unless otherwise noted.[(Hughes et al., 2002)] The effective dielectric constant of a sphere with an outer shell can be calculated with Eq. (10),where r is the total particle radius from center to outer shell, r is the radius of just the inner sphere, and ε and ε are the dielectric constants of the shell and core, respectively.

Results & discussion

Electrochemical analysis of the ring-disk (droplet) configuration

Before introducing the virus particles into the models, the baseline electrochemical properties are evaluated in the absence of virus particles. The disk and ring microelectrodes are located at the center of a 1 μl hemispherical droplet containing 1 mm [Fe(CN)6]4- (E  = 0.173 V vs. Ag/AgCl) in a supporting KCl electrolyte (1 mm) (see Fig. S2 for the effect of [Fe(CN)6]4- concentration). At the top, both electrodes have a mushroom-like cap that can be created via a template-driven overgrowth during electrodeposition, similar to reports from Garimella et al.[(Chen et al., 2015; Chandramohan et al., 2016)] An advantage of this overgrowth is that the spacing between generator and collector electrode caps (d ) is tunable and can be decreased beyond the predefined pillar spacing by controlling the electrodeposition current and duration.[(Cerquido et al., 2018)] Decreasing the generator-collector spacing is crucial for the diffusion-driven redox cycling process, [(Zafarani et al., 2016)] while the high aspect ratio of the three-dimensional electrodes enables better confinement of redox species [(Morita et al., 2014; Odijk et al., 2008)] and results in an increased amplification factor (), [(Zafarani et al., 2016)] where I and I represent the steady-state current with and without redox cycling, respectively.

Fig. 2a demonstrates the enhancement of the dual-mode Faradaic current, resulting from the recycling of redox-active species, compared to single electrode-mode. The linear sweep voltammograms (LSV) in Fig. 2a are carried out at a scan rate of 50 mV s−1 (see Fig. S3 for effect of scan rate) with E ranging from −0.2 to 0.6 V while E is either held constant at −0.15 V in dual-mode (see Fig. S4 for effect of different generator biases) or left floating in single-mode. The current shows the typical steady-state behavior encountered with microelectrodes. The subsequent concentration profiles for the oxidized species (c ) are shown in Fig. 2b and c at four points during the potential scan. As expected in the dual-mode configuration, a steep concentration gradient forming at E  = 0.2 V is seen between generator and collector electrodes that is not seen in single-mode. Because the flux of a species at an electrode is dependent on the concentration gradient, redox cycling can result in substantial amplification factors and explains the increased generator current seen in Fig. 2a.

Fig. 2

a) A representative voltammogram of the ring-disk microelectrode configuration in single- and dual-mode (1 mm [Fe(CN)6]4-, 1 mm KCl, scan rate = 50 mV s−1, E = −0.15 V or floating). The resulting concentration profiles of [Fe(CN)6]3- are shown for b) single-mode, and c) dual-mode over the course of the voltammetric scan. Note that the formal potential of the redox probe is around 0.17 V.

Typically, in electrochemical experiments, a sufficiently conductive electrolyte is used to minimize the ohmic drop in solution. Fig. S5 demonstrates the effect of KCl electrolyte concentration on the voltammetric current. As the KCl concentration increases, the steady-state current magnitude slightly decreases. This may be due to less migrative transport of the charged redox probe as electrolyte concentration increases. As a result, at lower KCl concentration, a larger portion of the current carrying ions will be [Fe(CN)6]3-/4- rather than K+ and Cl−. The lower electrolyte concentration also supports a larger electric field, which will be advantageous for electrophoretic transport of the virus particles discussed later. Moreover, because Joule heating of the electrolyte scales with conductivity σ (, where Q represents the heat generated), [(Cornejo and Linz, 2021; Yan et al., 2008)] it is advantageous to use a lower conductivity electrolyte to minimize these effects. Fig. S6 shows the temperature in the electrolyte after a 1200 s amperometric scan. Only minor temperature differences of a few mK are seen in the immediate region surrounding the electrode, which are not likely to cause any significant degradation in operational performance.

Fig. 3 highlights the effect of electrode geometry, which plays an important role in determining the redox cycling capabilities.[(Wolfrum et al., 2016; Ma et al., 2013; Zafarani et al., 2016; Xiong et al., 2015)] In particular, the impact of cap spacing (d , left column of Fig. 3), electrode height (h , middle column of Fig. 3), and electrode spacing (d , right column of Fig. 3) are evaluated. Regarding the cap spacing, the amplification factor Γ ranges from ∼4.2 at d  = 0.05 μm down to ∼2.05 for d  = 1.8 μm, the largest spacing studied. Likewise, the collector efficiency () follows a similar trend, ranging from 0.78 for d  = 0.05 μm down to 0.55 for d  = 1.8 μm. Although a smaller d leads to improved redox cycling behavior, decreasing d beyond a few 100 nm could be difficult to control during the overgrowth electrodeposition. To improve process robustness and consistency, spacings of 0.1 μm and above may be more suitable for experimental realization of the proposed system from a technical standpoint.

Fig. 3

The effect of generator-collector geometry on the redox cycling characteristics. a-c) The linear sweep voltammograms (Parameters: 1 mm [Fe(CN)6]4-, 1 mm KCl, scan rate = 50 mV s−1, E = −0.15 V) demonstrating the effect of cap spacing (d), electrode height (h), and electrode spacing (d), respectively. d-f) The subsequent amplification factors (Γ) and g-i) collector efficiencies (η) as functions of electrode geometry.

More substantial differences in η and Γ are seen when varying the electrode spacing d . For d  = 20 μm, Γ is less than 1.2 and η ∼0.13, whereas decreasing d to 0.6 μm increases Γ nearly four-fold to ∼4.8 and η to ∼0.8. Interestingly, around d  = 8 μm and larger, the cycling between collector and generator is minimized due to decreased overlap in diffusion layers (see Fig. S7). This is evident by the onset of an oxidation peak in the LSV data seen at 8 μm, and especially at 15 and 20 μm, as opposed to the steady-state sigmoidal behavior typically seen with microelectrodes.[(Xiong et al., 2015)] This points to mass transfer limitations similar to those seen with macroelectrodes. The onset of an oxidation peak can also be seen in the single-mode configurations presented in Fig. S8, thus signifying the importance of redox cycling in maintaining steady-state conditions and minimizing depletion of the redox species. To ensure adequate amplification and collection is upheld, an electrode spacing of 2 μm or less is recommended.

The electrode height (h ) does have implications on η and Γ, albeit much less of an impact than d . Across the range of h from 1 to 10 μm studied in this work, Γ only varies from ∼2.1 to 2.6 with the collector efficiency following a similar trend. However, this may be advantageous as increasing the height beyond 10 μm could present fabrication challenges associated with the photoresist thickness (i.e. electrodeposition mold layer). With the above factors taken into consideration, in addition to potential challenges that could occur during the device fabrication process, a generator-collector electrode spacing of d  = 2 μm, electrode height of h  = 8.5 μm, and cap spacing of d  = 0.5 μm is used for subsequent simulations.

Electrochemical analysis in the presence of virus particles (droplet configuration)

After evaluating the electrochemical properties of the microelectrodes, we introduce virus particles into the simulation to better understand the detection capability of the proposed system. The virus particles are modeled as spherical core-shell structure, each layer having unique dielectric properties.[(Hughes et al., 2002; MacCuspie et al., 2008)] Typically, the outer shell has a dielectric constant on the order of ε  = 5–10 while the ε is around 70–80.[(Hughes et al., 1998, Morgan et al., 1999, Yan et al., 2008)] Of particular interest for electrophoretic transport are the virus charge (q), virus count (N), and virus size (d ). Some reports have estimated the net surface charge of different viruses, such as Cowpea Mosaic Virus (CPMV), which is estimated to have a charge of (negative) 200e (where e is the elementary unit of charge) per virion (or 0.013 C m−2).[(Ermolina et al., 2006)] Similarly, for Cowpea Chlorotic Mottle Virus (CCMV), a surface charge of (negative) 150e - 200e is estimated based on electrophoretic mobility (μ ) data at neutral pH.[(Vega-Acosta et al., 2014)] For Tobacco Mosaic Virus (TMV) the surface charge density is even higher, at 0.043 C m−2[(Ermolina et al., 2003)] The charged nature of the virus particles serves as the basis for electrophoretic enrichment as demonstrated here.

We simulate the chronoamperometric response (see Figs. S9–S11 for unprocessed amperometric data) with four different particle charges – 0, −100, −200, and −300 (in units of e) – and three particle sizes – 30, 120, and 200 nm – to better understand how the charge and size of a particle impacts the sensor response, specifically particle capture time and number of capture events. The chronoamperometric simulations are performed in dual-mode configuration with E  = −0.15 V and E  = 0.6 V, which significantly reduces the virus capture time (from ∼740 s down to ∼20 s), due to a stronger electric field, and results in larger current steps (∼0.265 pA vs. 0.013 pA), due to the larger amplified current compared to single-mode (see Fig. S12 for results of simulations done in single-mode). As shown in Fig. 4 a, in the absence of any surface charge (q = 0), no virus particles attach to the electrode. As the charge increases, the electrophoretic force (F ) increases proportionally, resulting in more capture events (2 for q = −300 opposed to 1 for other charge values) and reduced capture time. The current steps shown in Fig. 4 are proportional to the surface area of the virus particle ()[ (Bonezzi and Boika, 2017)] and result in a ∼1 pA current step for a 70 nm particle, which is similar in magnitude to previously reported values using murine cytomegalovirus.[ (Dick et al., 2016)] More pronounced steps in current are seen for larger particles, shown in Fig. 4b. It is worth noting that viruses such as SARS-CoV and SARS-CoV-2 have been reported to be on the order of 70–100 nm, [(Goldsmith et al., 2004; Menter et al., 2020)] with reports ranging from around 50 nm for SARS-CoV [(Shieh et al., 2005)] up to 140 nm for SARS-CoV-2.[(Zhu et al., 2020)] Other viruses, such as herpes simplex virus, can be 250 nm in size [(Hughes et al., 2002)] while hepatitis A virus particles have been reported to be just 27 nm in size.[ (Grom et al., 2006)] Sepunaru et al. demonstrated single-particle capture events of influenza virus (measured size ∼ 125 nm) conjugated with silver nanoparticles.[ (Sepunaru et al., 2016)] They showed current spikes on the order of 10 pA due to oxidation of the nanoparticles upon capture at the electrode. In our case, for a particle size of 200 nm, we see step sizes of ∼5 pA, which is comparable to previous reports using silica nanoparticles.[(Boika et al., 2012)]

Fig. 4

The effect of various virus properties on the amperometric detection of capture events in a 1 μl droplet. a) The background subtracted current (ΔI) obtained from chronoamperometric measurements demonstrating the effect of the virus surface charge (q), b) The effect of virus size (d), and c) The effect of virus number (N). The solid and dashed lines represent the generator and collector current, respectively. Parameters: 1 mm [Fe(CN)6]4-, 1 mm KCl, E = −0.15 V, E = 0.6 V.

Interestingly, the time required for the first capture event increases with particle size, likely a result of the larger drag force (F ). Although particle size has implications in both the drag force () and dielectrophoretic force () acting on the particle, for DC measurements such as those here, we find F does not contribute significantly to the transport of the virus particles (see Fig. S13 for a comparison of average forces acting on the virus particles). Increasing the number of particles results in more capture events as shown in Fig. 4c, although no response is seen for N = 1, there is only one capture event for N = 5.

In an effort to improve the number of capture events and achieve the ultimate goal of single-entity detection, we confine the droplet volume from l μl down to 2 nl. Droplets of this scale and even smaller can be created and manipulated using a microcapillary pressure injection system.[(Salm et al., 2013)] Alternatively, controlled evaporation of the droplet has been utilized as a pre-concentration strategy as well.[ (Ebrahimi et al., 2013)] Nanoliter volumes are also commonly achieved in microfluidic channels which will be discussed in regard to the IDE structure in Section 3.3. Fig. 5 a shows the amperometric current for N = 2 for a particle size of d  = 70 nm in a 2 nl droplet. With 2 virus particles present, one particle is sufficiently far away such that it is not attracted to the electrode over the time frame tested. The other particle can be seen attaching after ∼40 min. The corresponding time evolution of the virus capture is shown in Fig. 5b with the electric field.

Fig. 5

a) The background subtracted current (ΔI) obtained from chronoamperometric measurements within a 2 nl droplet demonstrating the current steps caused by virus particle attachment to the electrode surface for N = 2. Solid and dashed lines represent the generator and collector currents, respectively. Inset shows a close up image of the capture event near 2375 s. b) The electric field at a sampling of time points showing the evolution of the particle attachment to the electrode surface. Parameters: 1 mm [Fe(CN)6]4-, 1 mm KCl, E = −0.15 V, E = 0.6 V, d = 70 nm, q = −200.

One possible drawback of performing electrochemical measurements in such confined geometries on the timescales shown here is the extent of the diffusion layer. Estimating the diffusion layer as , one obtains a layer of roughly 1 mm, well beyond the confines of the droplet. In the absence of redox cycling, the analyte supply is depleted after about 60 s, and the current diminishes to ∼0 (see Fig. S14). However, our configuration takes advantage of redox cycling to constantly resupply [Fe(CN)6]3-/4-, leading to a finite steady-state current and allowing for virus detection to continue within the droplet. Due to these constraints associated with the droplet geometry, the amperometric current of the 2 nl droplet is slightly less than the 1 μl droplet (see Fig. S15).

Electrochemical analysis in the presence of virus particles (IDE configuration)

After demonstrating the virus capture ability of an isolated generator-collector ring-disk pair in a droplet configuration, we investigate an interdigitated electrode (IDE) design, which is more suitable for integration with microfluidic channels rather than individual droplets. As discussed below, the IDE also improves detection at lower particle number, which may be an issue with the droplet configuration. The IDE structure is a common generator-collector configuration whereby sub-micrometer spacing (a) is traditionally achieved using more costly fabrication methods, such as electron beam lithography.[(Wolfrum et al., 2016)] The proposed template-driven overgrowth method is also applicable here and can be used to further reduce the lateral spacing between generator and collector. In our analysis, the IDE with overgrown cap is placed within a microfluidic channel of 15 μm in height, which inherently confines the virus particles close to the electrodes. Fig. 6 a shows the electric field distribution of an exemplary cross section of the IDE with electrode spacing d  = 2 μm, h  = 2 μm, and d  = 0.5 μm. Again, E  = −0.15 V and c  = 1 mm (see Figs. S16 and S17 for effect of collector bias and KCl concentration). Electric fields above 1 kV m−1 in magnitude are seen in the immediate vicinity of the electrode, which is sufficient for electrophoretic manipulation. Fig. 6b highlights the concentration profile of [Fe(CN)6]4- at the generator and collector. Significant overlap in the diffusion layers of adjacent electrodes is apparent, which allows the IDE to reach a steady-state current during the LSV scan that is not seen in single-mode. In fact, without redox cycling, the [Fe(CN)6]4- is depleted very rapidly during LSV, resulting in the current dropping to near 0 before the scan is complete (see Figs. S18 and S19).

Fig. 6

a) The electric field distribution of a cross section sample of the interdigitated electrode (IDE) and microchannel. b) The concentration profile of c showing the overlapping diffusion layers of adjacent generator and collector electrodes. c) The background subtracted current (ΔI) obtained from chronoamperometric measurements demonstrating the current steps caused by virus particle attachment to the electrode surface. Note that the attachment occurs within 1 s for both N = 1 and N = 2. Solid and dashed lines represent the generator and collector currents, respectively. Parameters: 1 mm [Fe(CN)6]4-, 1 mm KCl, E = −0.15 V, E = 0.6 V, d = 70 nm, q = −200.

The IDE structure is evaluated with N = 1 and 2 virus particles of 70 nm size and q = −200, similar to our previous studies with the droplet configuration (Section 3.2). The background subtracted current for both N = 1 and 2 is presented in Fig. 6c (the unprocessed chronoamperometry data are shown in Fig. S20). The particle attaches very quickly in the IDE configuration – after just 1 s. This is a significant improvement over the isolated generator-collector ring-disk geometry, which only recorded one capture event with N = 2 and no capture events with N = 1 in the 2 nl droplet. The main drawback of using such a large area electrode is the reduced blocking effect caused by the attachment of the virus particle compared with the much smaller ring-disk geometry. As a result, the relative reduction in current per virus particle attachment is smaller.

The proposed biosensor scheme has the potential to enable rapid and sensitive electrochemical “particle counting” with low fabrication costs. While the focus of this study was to demonstrate detection time and sensitivity, it can be complemented by other methods to enable selective detection as well. This can be achieved by incorporating biorecognition elements specific to the virus of interest. One example from Dick et al. used glucose oxidase (GOx)-modified antibodies that attached to the surface of the target virus.[(Dick et al., 2016)] As the virus-antibody-GOx unit adsorbed on the electrode surface, the GOx acted as a reducing agent for the FcMeOH+ being generated at the electrode surface, which allowed for re-oxidation of the FcMeOH and significant amplification of the signal. An alternative approach could be to extract the target virus from the sample matrix using magnetic nanoparticles (MNPs) conjugated with the corresponding antibody or aptamer. The virus-MNP conjugate can then be selectively resuspended in the test solution for downstream counting. If needed, by using a thermally reversible conjugation mechanism, such as the Diels-Alder approach (Gregoritza and Brandl, 2015), the virus particles could be removed from the MNPs via heating.

Conclusion

In this work, we propose a novel biosensor concept for single virus counting using two tunable generator-electrode configurations. The proposed design combines redox cycling with electrophoresis-driven electrode-particle collision for rapid and sensitive virus detection. The impact of various experimental and geometric parameters on the electrochemical signal is investigated using detailed finite element analysis. These designs (either droplet-based or channel-based designs) have the potential for integration with microfluidic and other lab-on-a-chip platforms for portable, point-of-care virus detection. Amplification factors and collector efficiencies are highly dependent on the generator-collector spacing and range from beyond 4 to as little as ∼1.1 for the amplification factor and from 0.8 down to less than 0.2 for the collector efficiency. While the redox amplification reported here is modest compared to previous reports which often use nanofabrication methods, we should emphasize that the dimensions of the proposed device is well within the capabilities of standard photolithography techniques. Moreover, the overgrowth strategy is compatible with scalable electrodeposition methods. Strategies to improve the amplification, such as nanocavity-based sensors, could be employed, but may limit the reach of the electric field into the bulk solution and make the electrophoretic enrichment of viral particles challenging.

The detection scheme detailed in this report relies on the virus particle blocking the mass transfer of the electrochemical redox probe upon colliding with the electrode, which causes a detectable step in the amperometric current. The nature of the collisions is found to depend on the virus surface charge, size, and quantity. Although the ring-disk design (droplet configuration) is more sensitive to virus attachment due to smaller surface area, the IDE design (microchannel configuration) takes advantage of a large capture area and particle confinement within the microchannel to detect virus particles on the order of seconds. The present FEA results provide valuable guidance for the practical fabrication of biosensors suitable for single particle counting, which has potentially broad applications ranging from electrochemical virus detection to immunoassays and molecular diagnostics.

Finally, while this work focused on establishing general design rules for sensing based on redox cycling and electrophoretic particle enrichment, future computational investigations into how the location of particle attachment impacts the measured current step size would be beneficial. Such studies would shed light on how particles of the same size can result in different current step magnitudes due to inhomogeneous flux of electroactive species across the electrode surface, therefore better quantifying the range of step sizes to be expected experimentally.

CRediT authorship contribution statement

Derrick Butler: Conceptualization, Methodology, Investigation, Visualization, Writing – original draft, preparation, Writing – review & editing. Aida Ebrahimi: Conceptualization, Methodology, Supervision, Visualization, Funding acquisition, Writing – review & editing.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.