Abstract

It has been predicted that the nanomaterials of graphene will be among the candidate materials for postsilicon electronics due to their astonishing properties such as high carrier mobility, thermal conductivity, and biocompatibility. Graphene is a semimetal zero gap nanomaterial with demonstrated ability to be employed as an excellent candidate for DNA sensing. Graphene-based DNA sensors have been used to detect the DNA adsorption to examine a DNA concentration in an analyte solution. In particular, there is an essential need for developing the cost-effective DNA sensors holding the fact that it is suitable for the diagnosis of genetic or pathogenic diseases. In this paper, particle swarm optimization technique is employed to optimize the analytical model of a graphene-based DNA sensor which is used for electrical detection of DNA molecules. The results are reported for 5 different concentrations, covering a range from 0.01 nM to 500 nM. The comparison of the optimized model with the experimental data shows an accuracy of more than 95% which verifies that the optimized model is reliable for being used in any application of the graphene-based DNA sensor.

1. Introduction

In the recent decade, DNA biosensors have been widely considered as a promising means of diagnostic prediction in genetic research such as cancer or hereditary disease due to their rapid and sensitive response, as well as the simple and convenient operation [1, 2]. The interface between nanomaterials and biosystems is emerging as one of the most various and dynamic areas of intense research [3]. In order to overcome the limit of planar semiconductor devices regarding sensitivity and potential for integration [4, 5], these nanomaterials are commonly proposed in sensor platform. Nanoscale channel confinement of planar FETs makes them extremely sensitive to electric perturbations, due to their one-dimensional structures. Therefore, obtaining a high sensitivity does not provide a worthy and reliable sensor. The implementation of nanomaterials of graphene in graphene field effect transistor (GFETs) constitutes an important step toward low-cost, easy handle, and highly sensitive molecular diagnostics [6, 7]. On the other hand, graphene with a two-dimensional thick monoatomic block of a carbon allotrope offers a large detection area (as seen in Figure 1) for DNA sensing compared with other known nanomaterials [811].

The research to date tends to focus on the impact of graphene-based materials in electronic and optoelectronic devices along with their potential applications [6].

2. Materials and Methods

2.1. Graphene-Based DNA Sensor Structure

Over the past few decades, a lot of information has become available on field-effect devices for DNA detection which is reported in the literature [1221]. DNA sensors usually presented in diverse sensor configurations such as electrolyte-silicon (ES) structures, depletion- and enhancement-mode FETs, floating-gate FET, and FET devices with or without a reference electrode. SiO2, silanized SiO2, SiO2-Si3N4, SiO-poly-L-lysine were employed as various gate-insulator materials with different thicknesses (from 2 to 100 nm) in various DNA immobilization methods (such as adsorption, covalent attachment, biotin-vidin complexation, and linker molecules) [22]. As reported in [16], the densities of the immobilized DNA varies from 2.4 × 108 to 5 × 1013 molecules/cm2 and hybridization buffer solutions with different electrolyte concentrations (from 10 nM to 1 nM) were examined to obtain the higher adsorption. In addition, the mostly reported flat-band or threshold voltage shift induced upon the DNA hybridization reach from several mV up to around 1.9 V. The focus of this theoretical study will be on developing the DNA-sensors-based graphene nanomaterials which have emerged as one of the most popular devices holding great promise to provide a straightforward platform of patience diagnosis [23, 24].

As shown in Figure 2, the schematic of the DNA sensor consists of a 300-nm SiO2 layer as a back-gate dielectric and a doped silicon substrate as the back gate which was proposed for modeling of DNA sensor in our previous studies [9, 25]. Based on that, graphene layer acts as a conducting channel connected to the source and drain electrodes which are made from Au. An Ag/AgCl reference electrode commonly acts as a gate for measuring the electrical characteristics in a liquid-gated or electrochemically gated configuration that controls the current along the graphene sheet between the two electrodes [2628]. Considering that the gate leakage effect is not considered in the current study, the DNA sensors were exposed to a phosphate buffer solution (PBS) containing the DNA molecules. In order to understand whether DNA molecules are adsorbed on the graphene surface or not, source-drain conductance were measured before and after exposing the DNA sensor to the different concentration of DNA molecules [3].

It is noteworthy to explain the DNA adsorption effect on nanomaterials of graphene surface, as well as the proposed model. The FET devices are very practical for measuring the surface charge and therefore the charge changes because a DNA adsorption event can be measured with a high accuracy. It has been demonstrated that the interaction between DNA molecules and graphene layer causes a significant change in the conductance of the graphene-based DNA sensor through the nonelectrostatic gating mechanism [29, 30]. In Figure 3, the effect of DNA concentration on the gate voltage at nanomolar concentrations is explored. It is generally suggested that the electronic-doping (n-doping) effect by the adsorbed charge species has a significant impact on conductance change of graphene-based transistors [31]. Moreover, the conductivity of the graphene-based DNA sensor is influenced by the increased number of carriers in the channel. As the number of DNA molecules increased from 0.01 to 500 nM, due to the negative charges of DNA molecules, it could be expected that the gate voltage would shift leftwards to the lower amounts [32, 33]. Based on the detection mechanism, we recently proposed an analytical model for detection of DNA molecules in which the DNA concentration was modeled by a gate voltage [9].

Although there are lots of works presented on the experimental progress, the detection mechanism is not understood quantitatively. On the other hand, modeling and simulation using partial differential equations (PDE) play a critical role in determining the current-voltage characteristics, sensitivity, and the behaviour of the sensing devices exposed to DNA molecules. Our proposed model is capable of performing the electrical detection of DNA molecules by modeling the conductance of the graphene sheets. The conductance of the large channel in graphene materials is obtained in the ohmic scaling law using Landauer formula [34].

Consider where is the electron charge, denotes Planck’s constant, represents the energy band structure, is the transmission probability, is the number of modes, and denotes the Fermi-Dirac distribution function [34, 35]. In other words, is the average probability of electron transmission in the channel from one electrode to another. Since the channel is assumed ballistic, the probability is considered equal to one [3638]. In order to highlight the effect of DNA adsorption on graphene, DNA concentration as a function of gate voltage is assumed and sensing factor is defined. High carrier mobility reported from experiments in the graphene leads to assume a completely ballistic carrier transportation in the graphene [39]. In the nonsaturation region, the DNA concentration model is employed as a function of gate voltage and the ideal current-voltage relation for the n-channel FET from [34] is modified as where  Å is carbon-carbon (C–C) bond length,  (eV) is the nearest neighbor C–C tight binding overlap energy, is the length of conducting channel, refers to the threshold voltage, is DNA sensing factor, and is DNA concentration. Different concentrations of DNA molecules were presented in the form of parameter. Thus, the DNA molecules adsorbed on graphene surface by iteration method are modeled as From extracted data, () parameters are calculated and obtained as , , and . Eventually, according to the proposed model of DNA sensor using nanostructured graphene layer, the current-voltage characteristic is modified as All theoretical characteristics of graphene FET for DNA concentration are plotted in [30] changing from 0.01 to 500 nM. It is stated that the sensor model with the proposed parameters represents the same trends as those reported as experimental data [3, 9]. Since the values of the parameters , , and in (3) are calculated based on trial and error, there is necessity of a methodological approach for obtaining a viable and accurate model which is reliable for being used in different applications of the graphene-based DNA sensor. To this aim, a swarm-based evolutionary algorithm (EA) called particle swarm optimization (PSO) is used for optimizing the mathematical model shown in (3). In important applications of DNA sensors such as hybridization detection, the behaviour of characteristic around is very important and determinant. Therefore, a proper and reliable model should be capable of providing the same features with experimental data around the mentioned area. The better the modelled curve fits the experimental result, the more the reliable model is achieved. In the current study, the optimized model using PSO provides a fitter characteristic compared with the previous mathematic model, especially around the point which is very important for the main applications of the model.

The PSO technique is widely used in optimizing different sorts of problems including fine materials, medical science, control theory, and energy issues [4043]. The important facts which make PSO popular among the researchers are its fastness, its ability to avoid being trapped in the local optima, and its capability of being employed in any type of optimization problem [4448]. Swarm intelligence techniques have several advantages that make them useful for optimization applications. Some of these advantages are scalability, flexibility, robustness, and low production costs. Scalability means particles can be added to the swarm without making any change in the overall behaviour of the colony. In other words, individuals typically have a limited range of interaction with their neighborhood, so that the introduction of a new particle does not directly affect the whole swarm [49]. That helps to change the scale and size of the swarm easily, based on the size of the search space defined. Flexibility gives the ability to the swarm to be adapted to different problems under different circumstances. A well-designed swarm algorithm can exploit and explore an unknown territory, even encountering a huge gap in the search space. With regard to the robustness, if a few faulty individuals exist in the swarm, they may not affect the global optimum found by the swarm. The fourth main advantage is having low production costs. This is respective to the computational cost and memory size per unit of the individuals which gives the possibility of being used in small cheap components. The main advantage of the particle swarm optimization technique over other global optimization methods such as simulated annealing is having a large number of agents which makes the particle swarm technique resilient to the local optima points of the search space. In other words, it helps the algorithm in not being trapped in the local optima points.

2.2. Particle Swarm Optimization (PSO)

The PSO algorithm is mainly based on the movement of a population of particles as a “Swarm.” This movement is inspired from a school of fish or a bird flock [5052]. The evaluation of the search space is being performed step by step. There are some parameters which need to be saved at each step, since they will be needed in calculating the parameters of the next step. Equations (4) and (5) show the velocity and location of th particle at th step.

Consider (number of particles); (maximum iteration number), where, is the iteration number; is the particle number; refers to the inertia weight coefficient which is normally chosen between 0.5 and 1.2; and are acceleration constants; and are random values between 0 and 1; is the velocity of particle at iteration ; and indicates the position of particle at iteration .

Besides the above parameters which should be calculated for each particle, there are some social parameters that lead the swarm to the region which contains the global optima. “Personal best” and “Global best” are two parameters which are being used in all the PSO versions. They are normally shown as “Pbest” and “Gbest,” respectively. The Pbest is known as the best location that each particle experienced so far during the runtime, and the Gbest is the best global optimum found so far by the swarm. Figure 4 shows a simple diagram of movement of a typical particle.

3. Implementing Particle Swarm Optimization in DNA Sensor

3.1. Representation of the Search Space

Swarm is the particles which are moving and giving solutions for solving the problem. The particles move in the domain of the problem space and each of them represents a solution for the problem. The parameters to be optimized in the DNA sensor model are , , and coefficients as in (2). Figure 5 illustrates a three-dimensional search space of the problem. If , , and are the vectors which return , , and , respectively, then particle flies in the search space to find the best possible values for the respective parameters. Vector is the movement velocity vector of the particle which is obtained from (4).

In the optimization process, each location represented by the particles is a possible solution. The search space consists of arrays and matrices. For example, is a matrix of the location of the particles in the proposed algorithm, while indicates the particle and is the index of the optimizing parameter. In Figure 5, the dimension of the matrix will be 200 × 3, if we consider the number of particles to be equal to 200.

3.2. Fitness Function and Parameter Selection

To evaluate the solutions proposed by the particles, we need to define a fitness function. The fitness function has to be able to determine which solution is better and more efficient after considering all the solutions obtained by the particles at each iteration. Normally the fitness function is being set to have the lowest possible value at an optimum point. In the current study, we also need to have the lowest possible value of the error between the experimental waveform and the proposed model; hence, the fitness function is proposed as follows: where represents the modelled current waveform for particle , is the experimental current of the DNA sensor, and is the fitness value of the th particle. Based on the fitness function chosen, the best solution will propose the most fitted curve for the current waveform which is desired for a suitable DNA sensor model. The parameters of the PSO algorithm are chosen based on the values tabulated in Table 1, and the strategy of the optimization is shown in a flowchart as in Figure 6.

3.3. Optimization Results for DNA Sensor Model

After the optimization process, the best values obtained for the parameters , , and are tabulated in Table 2. These values are the best among the 20 runs of the optimization algorithm, and the lowest fitness value is the indicator of choosing the best parameters. The lowest obtained fitness value is also shown in the table which determines how good the modelled waveform fits the experimental curve.

The train for the optimizing parameters is performed using the experimental curve of the bare DNA sensor. Figure 7 illustrates the experimental waveform and the optimized model of the probe DNA sensor.

To show the convergence quality of the optimization process using PSO technique, the convergence profile of the best fitness value is shown in Figure 8 which is based on the fitness value versus iteration. The graph shows that the algorithm could converge to the optimized values after around 60 iterations which show a good convergence speed of the proposed method.

Using the values obtained for , , and parameters, different concentrations of the DNA sensor are modelled and compared with the experimental waveforms with the same concentrations. To evaluate the quality of the obtained waveforms with respect to the experimental results, mean absolute percentage error (MAPE) index is used as an error evaluation parameter as follows: The results based on the MAPE index for different concentrations of the DNA sensor are shown in Table 3. The accuracy based on the MAPE value can be obtained as the result of deducting MAPE from 100 percent. Based on the results in the table, the accuracy of the optimized model is more than 95% for all the concentrations which is an acceptable range of accuracy.

Also, the optimized waveforms for different concentrations are seen in Figure 9, where for each concentration the optimized model and experimental waveforms are plotted. It needs to be mentioned that the middle part of the graph is very important, since the voltage used in the experiment is normally between 0 and 1 volts. The optimized model fits even better around the mentioned voltage range.

It needs to be mentioned that the (minimum gate voltage) is a very sensitive indicator for detection of DNA molecules and the optimized model shows a better performance around the (0 and 1 volt) voltage range.

4. Conclusion

Due to the unique characteristics of the graphene such as special C–C bonding, 1D structure, large surface-to-volume ratio, high conductivity, and biocompatibility, it has remarkable potential for biosensing material applications. The graphene-based DNA sensor has a similar structure to traditional metal-oxide-semiconductor FET (MOSFET) in which the gate is exposed to the electrolyte solution containing DNA molecules. In these sensors, DNA adsorption occurring on a graphene surface changes the characteristics of the transistors due to the intrinsic negative charge of DNA. In this research, DNA sensing factor which was presented in our previous work is optimized using particle swarm optimization (PSO) technique. The optimization process is performed in order to find the optimized values of fitting parameters in a current-voltage characteristic of DNA sensor model to predict graphene doping behavior in GFET devices. Based on the results obtained, an accuracy of more than 95% is obtained for the optimized model compared with the available experimental data. It can be concluded that the optimized grapheme-based DNA sensor model with the accuracy obtained is now suitable and reliable for being employed in different applications such as diagnosis of genetic and pathogenic diseases.

Acknowledgments

The authors would like to acknowledge the financial support from the Research University Grant of the Ministry of Higher Education of Malaysia (MOHE) under Project Grant: FRGS-4F149. Also, thanks to the Research Management Center (RMC) of Universiti Teknologi Malaysia (UTM) for providing an excellent research environment in which to complete this work.