Enhancement of Heterogeneous Microfluidic Immunosensors Using New Sensing Area Shape with Electrothermal Effect

Abstract

In heterogeneous microfluidic immunosensors, the diffusion boundary layer produced on the sensing area represents a critical factor that limits the biosensor performance. A three-dimensional simulation using the finite element method on the binding reaction kinetics of C-reactive protein (CRP) has been performed. We present a new microfluidic biosensor based on a novel reaction-surface design without and with electrothermal force. Two reaction surface configurations were studied. The kinetic reaction rate was calculated with coupled Navier−Stokes, mass diffusion, energy, and Laplace equations. The numerical results reveal that the characteristics of a microfluidic biosensor are more enhanced by using the circular ring design of the sensing area coupled with the electrothermal force. The rate of initial slope related to the association phase is multiplied by a factor 2 when the voltage is increased from 10 to 15 V. The results prove to be valuable in designing new microfluidic biosensors.

1. Introduction

During the last decade, the development of nanofluidic devices has grown considerably owing to their ability to control small amounts of fluid. Their unique characteristics have attracted researchers’ interest in several fields such as chemical, biological, and biomedical applications [1,2,3,4]. In the biomedical field, many biochemical processes such as mixing [5] and sensing [6,7] can be manufactured and integrated in a single chip. Microfluidic biosensors offer several advantages such as reduced consumption of the reagents, short analysis time, and simultaneous detection of multiple analytes [8]. Immunosensors represent a class of biosensors using the sensitivity of antibody-antigen interactions for the detection of appropriate analytes. This technique is used to quantify the proteins and biological small molecules. Two types of immunoassays have been found: homogeneous immunoassay, where the interaction between antibodies and antigens is made in solution, and heterogeneous immunoassays where the interaction is made between antibodies immobilized on a solid membrane and the antigen present at the boundary layer [9]. In heterogeneous immunoassays, the concentration of the binding complex of antigen–antibody on the reaction surface has a crucial role. The chemical reaction time scale and the transport phenomena rate by diffusion and convection occurring in a heterogeneous immunosensors are related by the Damköhler number (Da). When Da is higher, mass transport is limited, while association kinetic is lower for low Da values [10]. For the microfluidic biosensors, the ratio of the diffusion velocity of molecules to the reaction surface is comparatively small. A diffusion boundary layer is developed and it leads to limiting the efficiency of biosensor [11,12]. Several methods have been tested to improve the phenomenon of transport of the analytes to the reaction surface such as microfluidic confinement [13], and electrothermal effect [14,15,16,17].

In our previous studies, we proposed two methods to improve the reaction rate: one based on inserting a cylindrical obstacle within the microchannel near to the reaction surface and the other based on flow confinement [18,19,20]. These studies are compared to the electrothermal flow and show its efficiencies. However, these numerical simulations have been done for two-dimensional configurations. This is somewhat constraining for spatial factors and practical design, especially the sensing area.

In the present work, three-dimensional simulations are performed. The effect of the sensing area shape and the geometrical location of the electrodes on the microfluidic biosensor performance are investigated.

2. Biosensor Geometry and Theory

2.1. Biosensor Design

The primary purpose of this work is to enhance the biosensor performance by proposing a new design of reaction surface and electrodes. This system consists of mixing a small concentration of analyte with the water in a microchannel having a biosensor on which antibody ligands are fixed. A binding reaction between the antibody fixed at the reaction surface and the antigen is developed.

Figure 1 illustrates the three-dimensional schematic of microfluidic biosensors. The microchannel used in this study was 250 μm length, 50 μm wide and 40 μm high. The locations of reaction surface and electrode pair are shown in Figure 1. The geometric parameters of two biosensor models are given in

Table 1. Geometrical parameters for type 1 of the biosensor model.

Parameter	Value
External radius of electrode Rext (µm)	25
Internal radius of electrode Rint (µm)	20
Radius of reaction surface Rs (µm)	15.96
Bottom electrode angle α (°)	40→160
Top electrode angle θ (°)	40→160
Distance from the inlet X (µm)	100

and

Table 2. Geometrical parameters for the second biosensor model.

Parameter	Value
Positive electrode radius $R_{El,int}$ (µm)	10
Internal radius of reaction surface $R_s^{int}$ (µm)	15
External radius of reaction surface $R_s^{ext}$ (µm)	22
Internal radius of the negative electrode $R_{El}^{int}$ (µm)	25
External radius of the negative electrode $R_{El}^{ext}$ (µm)	30

, respectively.

Based on our previous studies [18,19,20,21,22] and literature studies [23], it appears that the formation of the diffusion boundary layer is closely related to the length of the reaction surface and the position of the electrodes. For this reason, the above configurations are proposed in order to improve the microfluidic biosensor response and reduce the diffusion boundary layer.

2.2. Governing Equations

The application of a non-uniform Alternating Current (AC) electric field on a fluid provokes inhomogeneities of the electrical permittivity ε and the conductivity σ of the solution due to the temperature gradient subsequent to the Joule effect. These inhomogeneities of ε and σ give rise to the electrothermal force given by this expression:

$${\overrightarrow{F}}_E=-\frac{1}{2}\left\{\left(\frac{\nabla \sigma }{\sigma }-\frac{\nabla \epsilon }{\epsilon }\right)·\overrightarrow{E}\frac{\epsilon \overrightarrow{E}}{1+{\omega }^2{\epsilon }^2/{\sigma }^2}+\frac{1}{2}{|\overrightarrow{E}|}^2\nabla \epsilon \right\}$$

(1)

where $\omega =2\pi f$ represents the angular frequency of the electric field $\overrightarrow{E}$, $\sigma$ is the electrical conductivity and $\epsilon$ is permittivity of the fluid. For aqueous media at 293 K, $\nabla \epsilon /\epsilon =-0.04\nabla T$ and $\nabla \sigma /\sigma =0.02\nabla T$ [24] where T is temperature of the fluid.

The electrothermal flow inside the microfluidic chips is assumed to be steady and laminar. The fluid is Newtonian. The flow is governed by the continuity, Navier-Stokes and the energy balance equations:

$$\{\begin{array}{l}\nabla ·\overrightarrow{U}=0\\ \rho \overrightarrow{U}·\nabla \overrightarrow{U}=-\nabla p+\mu {\nabla }^2\overrightarrow{U}+{\overrightarrow{F}}_E\\ \rho C_p\overrightarrow{U}·\nabla T=k{\nabla }^{2}T+\sigma {|\overrightarrow{E}|}^2 \end{array}$$

(2)

where $\overrightarrow{U}$, p, ρ, μ, k, and C_p are the velocity vector, pressure, density, dynamic viscosity, thermal conductivity, and specific heat of the fluid, respectively.

The transport of the antigen in bulk liquid phase from the left inlet is expressed by the Fick second law involving both conduction and convection:

$$\frac{\partial C}{\partial t}+\overrightarrow{U}·\nabla C=D{\nabla }^{2}C$$

(3)

Here, C is the bulk concentration (mol/m³), D is the diffusion coefficient of the antigen, and t is the time (s).

Equation (3) reveals clearly that the analyte concentration in the bulk is affected by the velocity field which is itself influenced by the electrothermal force described in Equation (1). Therefore, the concentration at the biosensor surface will be affected by this force.

The reaction between the antibody fixed at the reaction surface and the suspended target analyte (antigen) is modeled by the so-called first-order Langmuir adsorption equation. The binding reaction at the sensitive surface is symbolized according to:

$$C+\theta \rightleftharpoons C_s$$

where θ and C_s are the surface concentration (mol/m²) of the active sites and the antigen-antibody complexes, respectively. The rate constants related to adsorption and desorption are k₁ and k₂, respectively.

The rate equation governing the time evolution of antigen-antibody complexes C_s at the sensor surface is given by the following expression:

$$\frac{\partial C_s}{\partial t}=k_{1}C\left({\theta }_0-C_s\right)-k_2{C}_s$$

(4)

where θ₀ denotes the total number of active sites accessible on the sensor surface.

In order to solve the above governing equations (Equations (1)–(4)), we used the boundary conditions presented in

Table 3. Boundary conditions used in this simulation.

	Electric Field	Temperature	Velocity	Analyte
Inlet	Insulation	$\overrightarrow{n}·\left(k\nabla T\right)=q_{01}$	Parabolic profile	$C=C_0$
Outlet	Insulation	$\overrightarrow{n}·\left(k\nabla T\right)=q_{02}$	No viscous stress	$\overrightarrow{n}·\left(D\nabla C\right)=0$
Walls (except reaction surface and electrodes)	Insulation	Insulation	No-slip	Impermeable
Reaction surface	Insulation	Insulation	No-slip	$D\nabla C=\frac{\partial C_s}{\partial t}$
Electrodes	$\pm V_{rsm}$	$T=T_0$	No-slip	Impermeable

$\overrightarrow{n}$ is the unit normal vector to the surface, $q_{01}=k\left(T_{in}-T_0\right)$ and $q_{02}=k\left(T_{out}-T_0\right)$ are the convective heat fluxes at the inlet and outlet. Initially, the concentrations of analyte C and complex at the C_s reaction surface were set to zero.

.

In the above conditions (Table 3), different meanings may be assigned to the word insulation owing to the considered equation. In the matter of the boundary conditions related to the energy equation, this word means thermal insulation, that is, the conduction heat flux through the considered surface is equal to zero. Concerning the Poisson equation, electrical insulation is used at all the boundaries except the electrodes. In other words, at these boundaries, the electrical field is equal to zero to insure its continuity.

Concerning the analyte transport equation, the word impermeable means that there is no diffusion through the walls (except the reaction surface).

The numerical method used to discretize the above equations and perform the simulations is the finite element method.

Table 4. Numerical values of the parameters used in the current simulations [12,23,25].

Parameter	Value	Unit
Adsorption rate constant $k_1$	${10}^4$	m3/(mol·s)
Desorption rate constant $k_2$	$2.6\times {10}^{-2}$	s−1
Ligand concentration ${\theta }_0$	$1.4\times {10}^{-8}$	mol/m2
Diffusion coefficient $D$	$2.175\times {10}^{-11}$	m2/s
Inlet concentration $C_0$	6.4	μmol/m3
Thermal conductivity $k$	0.6	W/(K·m)
Fluid density $\rho$	1000	kg/m3
Dynamic viscosity $\mu$	$1.08\times {10}^{-3}$	Pa·s
Specific heat $C_p$	4.184	kJ/(kg·K)
Electrical conductivity $\sigma$	$5.75\times {10}^{-2}$	S/m
Relative permittivity ${\epsilon }_r$	80.2
Frequency $f$	100	kHz

illustrates the values of the physicochemical parameters used for this numerical study. These values are taken from the literature [12,23,25].

3. Results and Discussion

3.1. Model Validation

The chemical binding kinetics of the analyte (Anti-rabbit IgG) with the ligand (Rabbit IgG) at the reaction surface has been validated with experimental results of Hofmann et al. [13] without confinement flow and for the case of no AC electric field. Figure 2a presents the transient evolution of the normalized concentration of the complex. We notice good agreement between the simulation results with the experimental data of Hofmann et al. [13].

Figure 2b compares the results obtained by the present studies with those obtained by Huang et al. [23] in the cases without and with the electrothermal force. A good agreement has been observed. It should be noted that several components of the computation code has also been validated in our previous papers [18,19,20,21,22].

In order to optimize microfluidic biosensor performance, several parameters were tested in this work. Two biosensor geometries with two forms of reaction surface were studied.

3.2. Effect of Surface Reaction Shape

Figure 1 exhibits the geometry of the two biosensors having different shapes of reaction surface and electrode arrangement. Figure 3 and Figure 4 show the transient evolution of the CRP complex of these two biosensors for u_ave = 100 µm/s without applied voltage (V_rms = 0 V) and with applied voltage (V_rms = 15 V), respectively. The geometric parameters of the two biosensor models are given in Table 1 and Table 2.

For the first model of the biosensor, the reaction surface was a disk of radius R_s = 15.96 μm while for the second model it was a circular ring of outer and inner radius $R_s^{ext}$ = 25 µm and $R_s^{int}$ = 20 µm, respectively. In both models of biosensors, the area of the reaction surface was taken equal to S = 800 µm². The electrodes angle θ and α for the first biosensor model were equal to 40°.

Figure 3 and Figure 4 present the evolution of the normalized complex concentration C_s during the association and dissociation phases. During the first phase (association), the complex concentration increases almost linearly until reaching saturation. At this stage (saturation), all the ligands are occupied by the analytes injected at the inlet of channel. After a given time ($t$ = 800 s), the injection of analyte is stopped and therefore the kinetic bond dissociation will take place.

In the case of V_rms = 0 V (Figure 3), the binding reaction rate related to the circular reacting surface without the application of the electrothermal force was smaller than that related to the reaction surface having the form of a circular ring.

In the case of V_rms = 15 V (Figure 4), a significant improvement in the binding reaction was observed for the second model compared to the first model. The electrothermal force induced by the electrodes for the second model contributed effectively to the increase of the reaction rate. The shape of the reaction surface and the electrode arrangement have an important role in improving the biosensor response.

3.3. Effect of the Applied Voltage

Figure 5 illustrates the effect of electrodes angles (θ, α) on the transient evolution of CRP complex concentration for V_rms = 15 V and u_ave = 100 µm/s. Four values of the angle θ = α (namely, 40, 80, 120, and 160°) were considered in order to study the effect of these angles on the biosensor response. The results are also compared to those obtained in the case without the electrothermal effect (V_rms = 0 V). As expected, this figure reveals that when the electrodes angle increases, the electrothermal effect becomes more important and subsequently there is an improvement of the reaction rate.

Table 5. Initial slopes and enhancement factors of the binding reactions under 15 V for different values of the electrode aperture angles.

	Initial Slope (Association) × 10−11	Initial Slope (Dissociation) × 10−11	Enhancement Factor (Association)	Enhancement Factor (Dissociation)
$\alpha =\theta =40°$	4.37	3.37	1.22	1.14
$\alpha =\theta =80°$	4.85	3.66	1.36	1.24
$\alpha =\theta =120°$	4.43	3.98	1.52	1.35
$\alpha =\theta =160°$	5.90	4.25	1.65	1.44

illustrates the initial slopes for association and dissociation binding reaction corresponding to four values of the electrode aperture angles θ and α. The initial slope was calculated by choosing few points near the initial time of each phase (i.e., near $t$ = 0 s for the association phase and near $t$ = 800 s for the dissociation phase). Then, a linear interpolation based on the least square method was used to compute the slope. The slope of the curves gives the binding rate and is related to the response time of the biosensor.

The enhancement factors of these initial slopes are also presented. It is clear that the improvement of the binding reaction is significant during both the association phase and the dissociation phase. The most significant enhancement factor for the adsorption and desorption of CRP protein pair under 15 V was obtained for α = θ = 160°.

The initial slopes depend not only on the electrode aperture angles but also on the applied voltage. Figure 6 presents the three-dimensional surface plot of the slope at the origin of time related to the association and dissociation stages versus the applied voltage and the electrode aperture angle for a fixed value of the average inlet velocity u_ave = 100 µm/s. The large value of the initial slope was achieved for the largest applied voltage and for the largest electrode aperture angle.

Figure 7 illustrates the temperature rise versus the applied voltage and electrode aperture angle for u_ave = 100 µm/s. From this figure, it is noted that the temperature rise is strongly influenced by the applied voltage. However, the electrode aperture angle has a small and linear effect on the temperature rise.

Furthermore, it clearly appears that the temperature rise due to the Joule heating has a quadratic variation versus the applied voltage. This result is expected due to the source term in the heat equation. By increasing the angle θ, the electrode surface increases causing an increase in the electric field and subsequently the temperature. For example, the electrical field for V_rms = 15 V increases from 4 MV/m for θ = α = 40° to 6.2 MV/m for θ = α = 160°. Thus, there is an increase in Joule heating.

For microfluidic biosensors under AC electric field, the Joule heating could provoke an important temperature rise. The exceeding of the safe limit temperature can damage the thermally labile biological samples (say about 3 K) [23]. For this reason, it is necessary to work with a voltage lower than 15 V.

To avoid exceeding the limit temperature of the biological samples, we work at V_rms = 15 V. For this type of biosensor, the optimum efficiency is achieved for a larger electrode angle and high applied voltage.

Figure 8 depicts the temporal evolution of average surface CRP complex concentration without and with electrothermal effect for u_ave = 100 µm/s. Three different AC voltages, 5, 10, and 15 V, were tested in this study. We noticed a significant improvement in reaction rate, especially for the larger applied voltage. For example, the increase rate of initial slope for the association phase was about 5% when the voltage V_rms increased from 0 to 5 V and 105% when the voltage V_rms increased from 10 to 15 V (Figure 9).

In

Table 6. Enhancement factor for association and dissociation phases and temperature rise.

Applied Voltage (V)	0	5	10	15
Temperature rise (K)	0	0.3	1.3	3.0
Enhancement factor for association	-	1.05	1.88	3.84
Enhancement factor for dissociation	-	1.01	1.58	2.25

, we calculated the temperature rise for different applied voltages. For V_rms = 15 V, ΔT = 3.0 K represents the limit value of the biological samples. Indeed, the maximum temperature rise obtained for the first model of the biosensor for α = θ = 160° is ΔT = 2.7 K under the same conditions of applied voltage and inlet velocity.

3.4. Effect of Inlet Velocity

To study the impact of the average inlet velocity on the response of both biosensor models, Figure 10 shows the evolution of the slope at the origin of time related to the association phase versus the average inlet velocity with and without electrothermal field. In absence of the electrothermal effect (i.e., V_rms = 0 V), increasing the flow velocity improves the convective mass transport of the analyte towards the sensitive surface and subsequently enhances the reaction rate. In Figure 10, the initial slope related to the association stage with electrothermal effect (V_rms = 15 V) for the first model of the biosensor is presented for α = θ = 160°.

This figure exhibits two regions: a linear region for low inlet velocities and an almost saturation region for larger inlet velocities. In the first region, the reaction rate is greatly accelerated for any value of the applied voltage because the analyte in the bulk has time to diffuse to the reaction surface. In the second region, the enhancement factor of the binding reaction decreases because the amount of analyte carried by the flow is large [23]. We can conclude that the electrothermal effect is dominant for low inlet velocities and will be negligible for the large flow rates.

For comparison, the initial slope of the association phase for both biosensor models without and with electrothermal field versus the average inlet velocity reveals that the second model of biosensor is more efficient than the first one without and with an applied voltage.

4. Conclusions

The present paper performs three-dimensional simulations based on the finite element method and investigating the binding reaction kinetics of C-reactive protein CRP. A new microfluidic biosensor based on a novel reaction surface design without and with electrothermal force has been analyzed. Two reaction surface configurations have been studied: one in the form of a disc surrounded by two electrodes and the other in the form of a circular ring, one electrode of which is located inside the ring and the other outside. The effects of several important parameters were discussed, namely, the electrode angle, the applied voltage, the inlet flow velocity, and the reaction surface shape. The numerical results prove that the characteristics of the microfluidic biosensor can be further upgraded by using the second design of the sensing area (circular ring) coupled with the electrothermal force. In summary, the sensing area shape and the electrode location are able to improve the kinetic efficiency. The presented numerical simulation should be useful in designing microfluidic biosensors.