Abstract

This study draws attention towards the application of identification nanoparticle (NPs) sensor based on back propagation (BP) neural network optimized by genetic algorithm (GA) in the early diagnosis of cancer cells. In this study, the traditional and optimized BP neural networks are compared in terms of error between the actual value and the predictive value, and they are further applied to the NP sensor for early diagnosis of cancer cells. The results show that the root mean square (RMS) and mean absolute error (MAE) of the optimized BP neural network are comparatively much smaller than the traditional ones. The particle size of silicon-coated fluorescent NPs is about 105 nm, and the relative fluorescence intensity of silicon-coated fluorescent NPs decreases slightly, maintaining the accuracy value above 80%. In the fluorescence imaging, it is found that there is obvious green fluorescence on the surface of the cancer cells, and the cancer cells still emit bright green fluorescence under the dark-field conditions. In this study, a phenolic resin polymer CMK-2 with a large surface area is successfully combined with Au. NPs with good dielectric property and bioaffinity are selectively bonded to the modified electrode through a sulfur-gold bond to prepare NP sensor. The sensor shows good stability, selectivity, and anti-interference property, providing a new method for the detection of early cancer cells.

1. Introduction

In the past few years, people’s increasingly high requirements for microscopic detection technology have led to an evolution, in the applications of various nanomaterials in the biological field. NPs refer to particles with a particle size between 1 and100 nm, and it is a new material for research at the molecular level, cellular level, and subcellular level [1]. The NP-based nanotechnology can be applied to the diagnosis and treatment of tumors [2].

Fluorophores can be added to NPs by covalent bonding or embedding to form fluorescent NPs [3], which are used as fluorescent dyes for detection or labeling. Compared with ordinary fluorescent dyes, fluorescent NPs have high brightness, good stability, strong water solubility, and good biocompatibility [4]. Cancer cell imaging is of great significance in the diagnosis of tumors. As a carrier material, fluorescent NPs can enhance biological signals to facilitate detection [5]. Hyaluronic acid (HA) has a good application prospect in the diagnosis and treatment of tumor diseases [6]. A nanodrug-carrying system prepared by combining NPs and HA can specifically bind to receptor proteins on the cell surface to detect cancer cells. Shen et al. [7] prepared diamine-hyaluronic acid-iron-modified rhodamine fluorescent NPs for cancer imaging detection. The results proved that fluorescent NPs had high sensitivity and stability.

Sensors play an important role in various fields, and chemical sensors are generally used to detect chemical substances. Electrochemical biosensors use active substances in organisms as sensitive elements. It converts them into electrochemical signals to output characteristic signals such as potential, current, and capacitance [8]. The principle of electrochemical biosensor is as follows. When the active substance in the organism diffuses to the biosensitive layer, relevant chemical reactions occur. The biological signals are then converted into electrical signals that can be quantitatively analyzed. Finally, the concentration of the active substance can be identified by processing the electrical signals [9]. Nanomaterials have many functions in electrochemical sensors. They can modify electrodes, accelerate the rate of electron transfer, promote catalytic reactions, and immobilize biomolecules. They can further be used as carriers for markers, control switches for chemical reactions, as well as reactants [10]. The introduction of NPs into electrochemical biosensors can effectively improve its sensitivity and reproducibility.

The genetic algorithm (GA) is a probabilistic adaptive iterative algorithm that can process the discontinuous function through the global search [11], as it finds the global optimal solution by not being affected by gradient information. It can be used to optimize BP network [12]. When the optimized BP neural network based on GA is applied to the NP sensor, it endows the nanoparticle sensor with ideal stability by adjusting the weight and threshold [13].

In the past studies depicted in this section, it is revealed that alone neural network-based approaches do not provides satisfactory outcomes for finding optimal global solution for diagnosis of cancer cells and this is a great research challenge. In the study, the GA is used to optimize the BP neural network. The fluorescent NPs are applied in electrochemical sensors, with their basic characteristics expounded. The modified electrodes of Me-β-Cy/Au NPs/CMK-2/GCE are introduced in the biosensor which is applied in the early diagnosis of cancer cells, aimed at exploring the diagnostic value of fluorescent NP sensors for cancer cells.

2. Experimental Methodology

2.1. Me-β-Cy/Au NPs/CMK-2/GCE Electrodes
2.1.1. Preparation of CMK-2

(1) Preparation of the Phenolic Resin Prepolymer. 3.25 g of phenol is dissolved in 0.5 g of 20 wt% NaOH solution. 15 minutes later, the uniform mixture is slowly dripped into 5.5 g of 37 wt% formaldehyde solution, with the temperature adjusted to 75°C, followed by reflux for 1 h. After the mixture cools to 25°C, 0.4 mol/L HCl is used to adjust , followed by evaporation treatment. The mixture is then diluted with absolute ethanol. After the white particles of NaCl are removed by filtering, the remaining light-yellow solution is the phenolic resin prepolymer.

(2) Preparation of CMK-2. 1.5 g of polyether F127 is added to the solution prepared above, with the temperature adjusted to 45°C, followed by stirring for 30 min to form a light-yellow transparent solution. The mixture is then transferred to a petri dish at 45°C and dried for 8 hours to remove ethanol. The temperature is adjusted again to 100°C for thermal polymerization for 24 hours to form the transparent film, which is pulverized and fired at 350°C for 5 hours under nitrogen conditions to finally obtain dark brown powder. The powder is then dispersed in 98% concentrated sulfuric acid solution, followed by reflux for 24 hours. After the template polyether F127 is removed, CMK-2 is finally obtained, as shown in Figure 1.

2.1.2. Preparation of Modified Electrode

The glassy carbon electrode (GCE) is polished with 0.2 μm and 0.05 μm alumina wet powder, respectively, and then undergoes ultrasonic cleaning in absolute ethanol and ultrapure water, respectively. The control volume (CV) method is used for scanning activation in a 5 mmol/L sulfuric acid solution [14], with potential controlled at -0.5-1.6 V, and the scanning rate set as 100 mV/s. The electrode is then cleaned with ultrapure water and dried at room temperature.

1.8 mg of CMK-2 powder prepared above is dissolve in 1 mL of dimethylformamide (DMF), followed by ultrasound treatment for 1 h to form a uniformly dispersed black suspension. Subsequently, 2.5 μL of the suspension is applied to bare GCE, dried at room temperature, and recorded as a CMK-2/GCE electrode [15]. The above electrode is placed in 0.6 mmol/L chloroauric acid (HAuCl4) (Sigma-Aldrich, USA) solution for electrodepositing gold (Au NPs), with potential controlled at -0.5-1.6 V. The CV method is used to scan 10 times, with the scanning speed set as 50 mV/s, followed by rinse step using ultrapure water. It is then dried at room temperature and recorded as Au NPs/CMK-2/GCE. The above electrode is immersed in a phosphate buffer solution containing 1 mmol/L Me-β-Cy (Binzhou Zhiyuan Biotechnology, China) for assembling for 12 h, followed by rinse with ultrapure water. It is then dried at room temperature [16], denoted as Me-β-Cy/Au NPs/CMK-2/GCE, as shown in Figure 2.

2.2. The BP Neural Network Improved by GA

The BP neural network is a kind of feedforward neural network [17]. It is generally composed of an input layer, a hidden layer, and an output layer (Figure 3). The BP neural network has good self-learning and self-adaptation capabilities and strong fault tolerance [18].

In the study, the GA is used to optimize the initial weight threshold of the BP neural network, with the general initial weight search replaced by global search, so that the BP algorithm can accurately debug the neural network in the solution space, and search for the global optimal solution [19]. The optimization includes population initialization, fitness function selection, selection operation, crossover operation, and mutation operation [20], as shown in Figure 4.

After the training of the BP neural network, the sum of the residual absolute value between the predicted output and the expected output is set as the fitness function . where represents the coefficient, and it is 1 in the study. represents the value of the output node, represents the expected output value of the point, and represents the predicted output value of the point. The probability of each individual is expressed as follows: where represents the fitness function and represents the probability of an individual . In this study, the real number crossover method is applied, and the crossover operation equation of the chromosome and the chromosome at the point is as follows: where is within . If the gene of the individual is selected for mutation processing, the operating equation is as follows: where is the maximum of , is the minimum, is the number of iterations, is the maximum number of evolutions, and and are within .

The best optimized weight values are obtained by using GA methodology, and the weight and threshold values are optimized. It is further followed by the calculation of absolute values of residual, and then, the network weights are further modified. The accuracy requirements are then checked, and the entire optimization process of the BP neural network is then terminated.

2.3. The Fluorescent NPs Modified by HA
2.3.1. The Silicon-Coated Fluorescent NPs

1.5 mg of fluorescein isothiocyanate (FITC) (Sigma-Aldrich, USA) is dissolved in 15 mL of double-distilled absolute ethanol, followed by ultrasound treatment. The solution is then dripped into anhydrous ethanol without distillation treatment, followed by magnetic stirring for 24 hours. Subsequently, the mixture is centrifuged and washed with anhydrous ethanol and n-hexanol alternately until there is no fluorescence in the supernatant. It is then dried in vacuum and stored in darkness [21]. 7.5 mL of cyclohexane, 2.5 mL of n-hexanol, and 1.88 mL of Triton x-100 are mixed uniformly, followed by magnetic stirring until it is clear. Subsequently, 250 μL of ultrapure water is added, followed by stirring for 15 min to form a relatively stable microemulsion system [22].

Figure 5 shows the preparation process of FITC-APTES. 0.25 mg of FITC and APTES (Sigma-Aldrich, USA) is dissolved in a mixture of 0.6 mL cyclohexane and 0.6 mL n-hexanol, followed by ultrasound treatment to fully dissolve them. The mixture is then dripped into the microemulsion system prepared above. After stirring for 30 min, nucleation occurs. 125 μL mixture of ethyl orthosilicate and 3-aminopropyltriethoxysilane is added [23], together with 50 μL ammonia water. After magnetic stirring for 24 h, a few drops of acetone are added to demulsify. Subsequently, the sample is washed alternately with absolute ethanol and ultrapure water [24]. Finally, H2N-SiO2/FITC-APTES-NPs are obtained and stored in refrigerator at 4°C for later use.

2.3.2. The HA Covalent Modification

2.5 mg of sodium HA (NaA) is dissolved in 20 mL of ultrapure water by ultrasound treatment, with 0.02 mg of carbodiimide (EDC) and N-hydroxysuccinimide (NHS) (Ballingweier Reagent Company, China) added, respectively, followed by magnetic stirring for 2 h to activate the carboxyl group on HA [25]. 6 mL of the fluorescent NPs prepared above is then added, followed by stirring for 24 h under darkness. The ultrapure water is used for multiple washing. The unreacted HA and other small molecule impurities are removed through centrifugation [26]. Finally, HA covalently modified silicon-coated fluorescent NPs are obtained and stored at 4°C for later use.

2.3.3. Characterization and Determination of Silicon-Coated Fluorescent NPs

The electron microscope images of HN-SiO2/FITC-APTES-NPs are measured using JEOL-2100 electron microscope (JEOL Company, Japan). A Nicolet Avatar 370 DTGS infrared spectrometer (Nicolet Avatar, USA) is used to obtain an infrared spectrum. The fluorescence spectrum is measured with Hitachi F-4600 fluorescence spectrophotometer (Hitachi, Japan) instrument. LSM 5 Pascal Laser Module laser scanning confocal microscope (Carl Zeiss, Germany) is used to characterize the results of cell incubation. A Bio-Rad Model 680 spectrophotometer (Bio-Rad Company, USA) is used to test the MTT results.

2.4. Cultivation and Survival Rate of Cancer Cells
2.4.1. Cultivation of Cancer Cells

Human cervical cancer Hela cells were cultured in DMEM containing 12% fetal bovine serum, penicillin, and streptomycin at 37°C [27]. Hela cells in the proliferation phase were treated with trypsinization and seeded (104-105 cells/well). 60 mg/L H2N-SiO2/FITC-APTES-NPs were then added to incubate the cells for 4 h. After the supernatant was removed, the sample was rinsed with phosphate buffer and visualized using LSM 5 Pascal Laser Module laser scanning microscope.

2.4.2. The Survival Rate Detection

The H2N-SiO2/FITC-APTES-NPs were diluted using a phosphate buffer to configure solutions of 50, 100, 150, 200, 250, and 300 mg/mL, respectively. Hela cells were seeded at 104-105 cells/well in a 96-well plate, and the abovementioned different concentrations of phosphate buffer is then transferred to the 96-well plate, followed by incubation at 37°C overnight [28]. After the medium is removed, 10 μL of 1 mg/mL MTT (Sigma-Aldrich, USA) solution is added to each well to culture the cells for 4 hours, with the supernatant removed. 100 μL sodium dodecyl sulfate (SDS) was then added to each well to culture the cells for 15 minutes, so that the MTT can fully react with the cells. The Bio-Rad Model 680 microplate reader is used to measure the absorbance at 580 nm.

2.5. Statistics

In this study, the data was processed by SPSS19.0, and all experimental data were expressed as (). The root mean square (RMS) was applied to analyze the deviation between the predicted value and the expected value along with the evaluation of average absolute error (MAE). The -test was used for comparison between groups. was the threshold for significance.

3. Experimental Results

3.1. The BP Neural Network Optimized by GA
3.1.1. The Displacement Measurement

In the experiment of displacement measurement, the temperature in the control room is constant. The end surface of the probe and cellophane is wiped using alcohol to prevent dust and stains from affecting the measurement results. The angle between the probe and the target substance is 90°. Before the measurement, the laser is preheated. The probe fits closely to the aluminum workpiece surface, and the movement distance is controlled at 0.05 mm each time. The upper position machine is used to read the power value of the channels CH1 and CH2, and a curve is drawn to show the relationship between the displacement and the received power (Figure 6).

3.1.2. Comparative Analysis of Optimized GA-BP and Traditional BP Approach in Terms of Accuracy and Error Value

In the parameters of the BP neural network, 2 means input and 1 means output. 7 hidden layer nodes are set, and the target error is 0.0001. The number of training is 1500, and the training function is Trainingfu. The population size is set as 100 by GA, with a crossover probability of 0.2 and a mutation probability of 0.05. After setting of parameters, the power value and displacement value are used as the input and output of the neural network for repeated training, respectively. At the same time, the relevant parameters of the network are adjusted to make the training result reach the predicted value. To test the data, the electric displacement stage moves at 0.1 mm/step. The power meter is used to obtain the predicted output value and the predicted output value of the BP neural network optimized by GA.

Further, for evaluating the accuracy rates of traditional and optimized GA-BP-based approach, the accuracy parameters have been evaluated for the testing displacements.

From the graphical representation in Figure 7, it is depicted that the accuracy value increases with a small change in testing displacement for both the traditional and optimized cases. However, it should be noted that the value of accuracy for the optimized GA-BP neural network is much better compared to the traditional BP neural network-based approach. The final accuracy value obtained for the GA-BP method is 94.54% which is superior to the final accuracy value of 86.78% obtained for the traditional BP neural network method.

For the comparison of the traditional BP neural network and the GA optimized BP neural network, the RMS and MAE error values are introduced. It is found that the RMS and MAE of the optimized BP neural network are smaller, as shown in Figure 8.

Figure 8 depicts the error value with respect to testing displacement in (mm) for the normal BP predicted value and the true value. It also depicts the error rate between the GA-BP predicted values and the true values for the corresponding testing displacement. It can be seen that the error value is minimized near 0 for the case of the GA-BP-based approach which is not minimized much while utilizing the BP predictions alone. Also, it can be noted that the variations in the error values are minimized while utilizing the GA-based optimized BP neural network approach, as compared to the tradition BP neural network-based method.

3.2. Characterization of Silicon-Coated Fluorescent NPs
3.2.1. Morphology and Particle Size

It is evident from Figure 9(a) that the silicon-coated fluorescent NPs have good morphology, good dispersion, and relatively uniform particle size, which are sphere in shape. It is evident from Figure 9(b) that the particle size of NPs is about 105 nm.

3.2.2. The Infrared Spectrum

The infrared spectrum is obtained for measurement of FITC, FITC-APTES precursors, H2N-SiO2/FITC-APTES-NPs, and HA-SiO2/FITC-APTES-NPs. It is evident form Figure 10(a) that the FITC has the characteristic peak of the isothiocyano group at 2186 cm-1 where the FITC-APTES precursor does not possess the characteristic peak. It further indicates that a coupling reaction occurs between FITC and APTES, the isothiocyano group is combined with the guanidine group, and the FITC-APTES precursor is successfully prepared.

It is evident from Figure 10(b) that H2N-SiO2/FITC-APTES-NPs have the characteristic peak of amino group appears at 3433 cm-1. The characteristic peak of silica at 1100 cm-1 shows that the prepared NPs contain amino groups. The HA-SiO2/FITC-APTES-NPs have the stretching vibration peak of NH at 3058 cm-1. The stretching vibration peak of CH is at 2990 cm-1, and characteristic peaks of carboxyl groups are at 1711 cm-1 and 1736 cm-1. These peaks are all caused by the group’s stretching vibration, indicating that HA successfully binds to the surface of fluorescent NPs through amide bonds.

3.2.3. The Fluorescence Property

It is evident from Figure 11(a) that after FITC reacts with APTES to form the FITC-APTES precursor, the overall fluorescence intensity decreases. Though the fluorescence intensity of H2N-SiO2/FITC-APTES-NPs and HA-SiO2/FITC-APTES-NPs decreases slightly, the peak value is still comparatively high. This indicates that the H2N-SiO2/FITC-APTES-NPs and HA-SiO2/FITC-APTES-NPs have relatively good fluorescence properties. The fluorescence confocal microscopic image of HA-SiO2/FITC-APTES-NPs shows that the fluorescence NPs are uniformly distributed (Figure 11(b)).

3.2.4. The Chemical Stability

Whether the formation of FITC-APTES will cause the leakage of fluorescent dyes is explored. The HA-SiO2/FITC-APTES-NPs are immersed in ultrapure water for 80 h and compared with the HA-SiO2/FITC-NPs for fluorescence intensity. It is found that the fluorescence intensity of HA-SiO2/FITC-NPs decreases obviously in the first 10 hours, indicating that the dye leakage is serious. The curve gradually becomes stable after 10 hours. The fluorescence intensity of HA-SiO2/FITC-APTES-NPs decreases slowly. The relative fluorescence intensity of HA-SiO2/FITC-APTES-NPs in each time period higher than that of HA-SiO2/FITC-NPs is maintained over 80% (Figure 12). It shows that the HA-SiO2/FITC-APTES-NPs have good fluorescence stability.

3.2.5. Detection of Hela Tumor Cells

During diagnosis of tumor diseases, specific markers are often used as auxiliary means. HA can specifically recognize the unique HA adhesion factor protein receptor on the surface of cancer cells and binds to the surface of cancer cells. In this study, the HA-modified fluorescent NPs were visualized using a laser confocal microscope (Figure 13). The left was a bright-field optical image, and the right was a dark-field fluorescent image. The cells in the two corresponded to each other. In the fluorescence bright-field image, there was obvious green fluorescence on the surface of cancer cells, which indicated that HA-modified fluorescent NPs covered the surface of cancer cells (Figure 13(a)). In the dark-field image, the cells emitted bright green fluorescence (Figure 13(b)). Therefore, the HA-modified fluorescent NPs are a new method for early diagnosis of cancer cells.

3.2.6. Cell Survival Rate

In this study, the MTT method was used to determine the survival rate of cells. After MTT treatment of living cells, the mitochondrial-specific dehydrogenase in the cell would cause MTT to form water-insoluble blue-purple crystal formazan, which would not be formed in dead cells. To measure the absorbance of different concentrations of HA-modified fluorescent NPs at 580 nm, cells cultured in DMEM were set as blank group, and the cells cultured with fluorescent NPs were set as the control group. A lower survival rate indicates greater toxicity of HA-modified fluorescent NPs to cancer cells. It was found that as the concentration increased, the survival rate of cancer cells was gradually decreasing (Figure 14). To maintain good fluorescence intensity and cell survival rate, 50 mg/mL fluorescent NPs were selected for cell incubation experiments.

3.3. The Characterization of the Modified Electrode
3.3.1. Characterization of CMK-2

The Rigaku Corporation SA-HF3 (Rigaku Corporation, Japan) is used to obtain the X-ray diffraction image of CMK-2. It is found that there is an obvious diffraction peak that appears at 2θ, indicating that CMK-2 has an obvious two-dimensional hexagonal structure (Figure 15(a)). The JEOL-2100 instrument (JEOL Company, Japan) is used to obtain the transmission EMS image of CMK-2. It is noted that it has a clear stripe structure and uniform size, indicating that the CMK-2 has a relatively good stripe structure (Figure 15(b)).

3.3.2. Electrochemical Properties of Modified Electrodes

Electrochemical impedance spectroscopy (EIS) can be used for the characterization of different electrodes. In this study, the impedance change during the modification of the electrode is explored, and further, the electrochemical impedance spectroscopy is observed. The electrochemical impedance of bare glass carbon electrode (GCE) CMK-2/GCE, Au NPs/CMK-2/GCE, and Me-β-Cy/Au NPs/CMK-2/GCE was measured in 1 mmol/L Fe(CN)63-/4- phosphate buffer (Figure 16(a)). The frequency change range is 0.01-10 kHz.

The electron transfer resistance of bare GCE is very small, about 500 Ω. After it is modified by CMK-2, the electron transfer resistance rapidly increases to 3800 Ω. The reason may be that the mesoporous silicon material itself does not conduct electricity. After electrodeposition Au, the impedance value is reduced, which shows that the Au NPs has good conductivity and promotes the transfer of Fe(CN)63-/4-. After it is combined with Me-β-Cy, the electron transfer resistance continues to increase to 2100 Ω. It may be that the organic matter Me-β-Cy inhibits the transfer of Fe(CN)63-/4-. Taken together, it can be confirmed that Me-β-Cy/Au NPs/CMK-2/GCE is successfully prepared in the study. Figure 16(b) shows the CV spectra of bare GCE, CMK-2/GCE, Au NPs/CMK-2/GCE, and Me-β-Cy/Au NPs/CMK-2/GCE in 1 mmol/L Fe(CN)63-/4- phosphate buffer solution. The impedance changes in Figure 16(a) corresponded to that in Figure 16(b).

3.3.3. The Peak Current and Scanning Speed

Figure 17 shows the CA spectra of Me-β-Cy/Au NPs/CMK-2/GCE in a phosphate buffer at different scanning speeds. The scanning speed is controlled at 100-400 mv/s, and the results show that the scanning speed continues to increase providing the current peak value highest at 0.35 V.

3.3.4. The Optimal pH of the Solution

In this study, the influence of the pH value of the solution on the modified electrode is determined along with the relationship curve between the peak current value and pH drawn. The base solution is 15 μmol/L phosphate buffer with different pH levels. It is noted that as the pH value in the solution gradually increases from 1.0 to 4.0, the peak current signal also increases. When the pH value exceeds 4.0, the response signal decreases rapidly (Figure 18(a)). It may be that when the acidity of the solution increases, the number of protons is gradually decreasing, affecting the reaction on the electrode. When the pH is controlled in the range of 1.0-6.0, the reduction peak potential has a linear relationship with the pH value of the test solution. The slope of the fitting line is -60 mV, which is very close to the theoretical value -59 mV. It suggests that during the electrochemical reaction, the number of protons is equal to that of electrons (Figure 18(b)). When , the peak current is the highest. Hence, was selected as the best pH value.

3.3.5. The Time for Modification

Me-β-Cy bonds to the gold NPs through S-Au bond, but this kind of covalent bond takes a certain time to achieve a stable state. Therefore, the time needed to reach a stable state is explored during Me-β-Cy modification. The reduction peak current in 15 μmol/L phosphate buffer is recorded with time, as shown in Figure 19.

The results show that when the modification time of Me-β-Cy is less than 15 h, the reduction peak current gradually increases with the time. After 15 h, the peak current tends to be stable. In this study, the modification time is set as 15 h.

4. Conclusion

In the study, the GA is used to optimize the BP neural network, and the identification NP sensor is used to diagnose cancer cells, aimed at realizing the early diagnosis of cancer cells by fluorescence imaging. The RMS and MAE of the optimized BP neural network are analyzed which were observed to be smaller than the traditional one. Furthermore, the modified electrode Me-β-Cy/Au NPs/CMK-2/GCE is successfully prepared and applied to the sensor. This sensor shows good stability, selectivity, and anti-interference property, providing a new method for the detection of early cancer cells. The prepared HA-modified fluorescent NPs have good fluorescence stability, and HA covalently modified fluorescent NPs can specifically bind to cancer cells. This method can be used in fluorescence imaging detection. This study provides a good theoretical basis for clinical diagnosis and prevention of cancer. However, the practical implementation of results using cervical cancer cells are universally applicable to other tumor cells, and its further verification will be the future research directions of this research work.

Data Availability

All data has been provided in the main manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.