Next Article in Journal
Spatial Variability of Aroma Profiles of Cocoa Trees Obtained through Computer Vision and Machine Learning Modelling: A Cover Photography and High Spatial Remote Sensing Application
Previous Article in Journal
Hierarchical Identifier: Application to User Privacy Eavesdropping on Mobile Payment App
Previous Article in Special Issue
An Adaptive Bluetooth/Wi-Fi Fingerprint Positioning Method based on Gaussian Process Regression and Relative Distance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Method for Quantitative Estimation of Glucose and Fructose Concentrations in Aqueous Solutions Based on Infrared Nanoantenna Optics

1
Institute for System Dynamics and Research Center SCoPE, University of Stuttgart, Waldburgstr. 17/19, 70563 Stuttgart, Germany
2
4th Physics Institute and Research Center SCoPE, University of Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(14), 3053; https://doi.org/10.3390/s19143053
Submission received: 28 February 2019 / Revised: 19 March 2019 / Accepted: 24 March 2019 / Published: 11 July 2019

Abstract

:
In life science and health research one observes a continuous need for new concepts and methods to detect and quantify the presence and concentration of certain biomolecules—preferably even in vivo or aqueous solutions. One prominent example, among many others, is the blood glucose level, which is highly important in the treatment of, e.g., diabetes mellitus. Detecting and, in particular, quantifying the amount of such molecular species in a complex sensing environment, such as human body fluids, constitutes a significant challenge. Surface-enhanced infrared absorption (SEIRA) spectroscopy has proven to be uniquely able to differentiate even very similar molecular species in very small concentrations. We are thus employing SEIRA to gather the vibrational response of aqueous glucose and fructose solutions in the mid-infrared spectral range with varying concentration levels down to 10 g/l. In contrast to previous work, we further demonstrate that it is possible to not only extract the presence of the analyte molecules but to determine the quantitative concentrations in a reliable and automated way. For this, a baseline correction method is applied to pre-process the measurement data in order to extract the characteristic vibrational information. Afterwards, a set of basis functions is fitted to capture the characteristic features of the two examined monosaccharides and a potential contribution of the solvent itself. The reconstruction of the actual concentration levels is then performed by superposition of the different basis functions to approximate the measured data. This software-based enhancement of the employed optical sensors leads to an accurate quantitative estimate of glucose and fructose concentrations in aqueous solutions.

1. Introduction

The reliable detection and identification of minute analyte concentrations in complex environments has been a long-standing problem, in particular in health care, life science, and point of care applications. In recent years the detection of glucose has attracted particular interest, as it is a highly important molecule in the human diet and metabolism. The, ideally, continuous monitoring of blood sugar levels is crucial for the treatment of, e.g., diabetes mellitus, with an ever growing number of patients [1,2,3,4]. This sensing task is, in fact, a very challenging one. Not only are the concentrations of glucose in the human blood or in the ocular fluid (as an easy-to-access system) very small, their detection is also hurdled by the extremely large number of additional molecular species in these body fluids [5,6,7]. A particular challenge lies with the other monosaccharides, which have extremely similar physical properties and might obstruct the real glucose concentration.
A few particular interesting sensor concepts are related to plasmonics-based systems. Plasmonics, or noble metal nanoparticle optics, utilize the large number of quasi-free conduction electrons in noble metals [8,9]. If a gold nanoparticle is illuminated by an external light field of appropriate wavelength, the quasi-free conduction electrons oscillate with respect to the fixed ionic background of the particle. As a consequence, this resonance is associated with so-called local near-fields in the close vicinity of the particle, which are significantly stronger than the far-field plane wave intensity. These local fields can be utilized for so-called refractive index sensors [10,11,12]. If an analyte molecular enters the near-fields, the far-field resonance will spectrally shift. This shift can be used as a highly sensitive sensor down to a few or even a single molecule [13,14]. However, without appropriate functionalization, such a sensor cannot distinguish different analyte molecules, as it only reports refractive index changes. While functionalization is possible, it is often highly complex, still sensitive to a certain cross-talk with similar molecules, and also tends to degrade rather quickly [15,16]. All of these points are, in particular, valid for complex and highly multi-entity environments, such as the human blood and any other body fluids.
Infrared spectroscopy, on the other hand, is known for its supreme specificity [17]. The infrared spectra of biomolecules show characteristic absorption features which allow to unambiguously identify molecular species. This spectral region is therefore often called the fingerprint region and the spectral signatures the molecular fingerprints. However, while specificity is excellent, there is little sensitivity. In general, large amounts of analyte are needed in order to acquire infrared absorption spectra. Here, again, the complex nature of the investigated samples comes as a challenge.
Fortunately, it has been shown that the best of both of these techniques can be combined in surface enhanced infrared absorption (SEIRA) spectroscopy [18]. SEIRA spectroscopy has been studied in great detail in order to uncover the basic working principle [19]. Generally speaking, the utilized plasmonic nanoparticles or nanoantennas are designed such that their plasmonic resonances are close to the molecular absorption bands. When the target molecular species enters the locally enhanced near-fields of the antennas, the vibrational modes and the plasmonic mode couple. As a consequence, the very strong and bright optical response of the plasmonic antenna is imprinted with the signature of the molecular species [19]. An enhancement of up to six orders of magnitude can be observed [18]. This allows to measure very small amounts of analyte in nanoscale volumes while simultaneously allowing to unambiguously identify the molecular species based on the unique vibrational signatures [20,21,22,23]. The basic working principle of SEIRA has been fully understood, making it a routine method for the investigation of other systems, species, or entities. Besides the observation of molecular species, for example during folding and unfolding of proteins [24,25], SEIRA has also been used for the investigation of the plasmonic resonances themselves, e.g., for the extend and intensity of the local near-fields [26,27].

2. Materials and Methods

Glucose and fructose, as our model system, show, in fact, distinct vibrational modes in the infrared spectral range [28]. This fact allows us to distinguish these two species. Figure 1 illustrates our measurement concept. Resonant gold nanoantennas are incorporated into a flow cell which allows us to flush different aqueous analyte solutions over the antennas while simultaneously acquiring infrared reflection spectra. The reflection spectra contain the signature of the plasmonic resonance as well as the imprinted vibrational modes of the analyte molecules. As sketched, these differ for glucose and fructose. Apart from the mere presence of the analyte molecules, it should be possible to also extract the concentration of the molecules from the measurements as a higher concentration leads to a larger modulation depth of the vibrational features.
While the general concept is straightforward, there are in fact a number of challenges when aiming for a quantitative analysis. Firstly, the contribution of the plasmonic resonance has to be removed from the data. This is in general achieved by the so-called baseline correction. However, the exact modulation of the remaining features is found to depend significantly on the exact baseline. While this is of no further concern for identifying a molecular species, it is a significant hurdle in determining quantitative concentrations. We thus need a stable routine that allows us to “calibrate” our sensor scheme for concentration levels, that is, correlate the modulation depth of the molecular features with the concentration of molecules present. This calibration must not be disturbed by varying baseline corrections. We are going to show that we are indeed able, using a global baseline correction and an adaptive algorithm for basis function approximation, to extract these absolute values.
Within the scope of this work, two different monosaccharides are examined with the proposed method. Due to their importance in health science, aqueous solutions of monosaccharaides are investigated. Glucose is a highly important biological molecule and its concentration in the blood or other body fluids needs to be determined with high accuracy. The primary hurdle lies in the presence of other biomolecules in, e.g., the blood or tear liquid. In particular, the other monosaccharaides cause problems in identifying the correct concentration [29]. As a model system, we therefore utilize aqueous solutions of different glucose and fructose concentrations in order to establish a robust routine for quantitative concentration determination. The proposed procedure to achieve this goal is schematically depicted in Figure 2.
As mentioned above, our sensor data is obtained from SEIRA measurements. We collected a large set of reflective SEIRA spectra of different pure glucose and fructose as well as mixed solutions, as published elsewhere [28]. In brief: We utilize a reflection flow cell in inverse geometry, which is flushed via attached tube connectors that transport the desired solutions into and out of the flow cell. The key parts of our sensor are the different linear gold antenna arrays which were fabricated with electron beam lithography on top of IR transparent calcium fluoride substrates. The choice of the plasmonic element is motivated by the following aspects: Ease of structure fabrication, chemical stability, quality of the plasmonic resonance in terms of linewidth and amplitude, as well as resonance position. While there are many highly sophisticated structures and systems, which include, among the noble metals, also graphene [20,30,31,32,33,34], we choose the most straightforward realization via dipolar gold nanoantennas. These antennas have shown excellent quality factors, large dipole moments, are chemically inert, and can be easily fabricated via many bottom-up and top-down techniques. All of these aspects are also highly relevant for future applications. The geometrical parameters of the nanoantennas are chosen such as to exhibit a plasmon resonance at the spectral position of the targeted molecular vibrations of glucose and fructose, resulting in a length of 3500 nm, width, and thickness of 100 nm (over a 2 nm chromium adhesion layer), periodicity 4500 nm in x direction and 3000 nm in y direction [35,36]. For the spectral measurements, we use a commercial FTIR spectrometer (Bruker VERTEX 80, Bruker Optik GmbH, 76275 Ettlingen, Germany) and an optical microscope (Bruker Hyperion 2000, Schwarzschild objective with 15-fold magnification, NA = 0.4, Bruker Optik GmbH, 76275 Ettlingen, Germany). The spectra are measured with a nitrogen-cooled mercury cadmium telluride (MCT) detector and referenced to a gold mirror. The measurements spot is about 90 μm × 90 μm. Top-down electron beam lithography is a very well established technique leading to highly uniform structures which result in highly uniform optical properties of the elements. Additionally, the measured area contains a significant number of elements, leading to an effective ensemble averaging, minimizing the contribution of the individual, single object.
In general each SEIRA spectrum is a combination of the signature of the plasmonic resonance and the vibrational features which are characteristic for the molecular species. Since we are interested in the pure vibrational spectrum of the examined specimens, the plasmonic background has to be removed first. This can be achieved with, e.g., asymmetric least squares smoothing (ALSS) for baseline correction [37]. Thereby, it is possible to reconstruct the unperturbed plasmonic resonance or at least approximate it sufficiently. Afterwards, the pure vibrational data r B C is obtained by dividing the measured SEIRA spectrum by the reconstructed baseline. The so-called baseline corrected spectra, then, only contain the pure vibrational information showing up as peaks at their specific wavenumbers ν i . In principle, these peaks contain not only information about the presence of the respective molecule but also information about the absolute concentration. In order to extract this quantitative information, we developed an algorithm to evaluate the contribution of each specimen to the height of these characteristic peaks. The basic assumption for our algorithm is that a superposition principle can be applied. If this is valid, each specimen can be evaluated separately and the overall sum of all individual parts represents the measured height.
We expect the contribution of pure water to be constant ( k i ), whereas the contribution of glucose ( φ i ) or fructose ( ψ i ) depends on the concentrations ( c ) of the dissolved monosaccharides. φ i and ψ i are the so-called basis functions which can be, in the most general form, arbitrarily shaped (e.g., polynomial, exponential, or other). To model the correlation between the level of monosaccharide-concentration and the resulting height of the peaks at the characteristic wavenumbers we can choose a suitable basis function for each specimen. This can be written as
Δ s g ( ν i ) = φ ˜ i ( c g ) k i = φ i ( c g )
and
Δ s f ( ν i ) = ψ ˜ i ( c f ) k i = ψ i ( c f ) ,
resulting in
Δ s ( ν i ) = Δ s g ( ν i ) + Δ s f ( ν i ) = φ i ( c g ) + ψ i ( c f ) ,
where Δ s is the height of the peak at the characteristic wavenumber. As we want to determine the absolute concentrations from our measurements, one has to solve a system of nonlinear equations for arbitrary basis functions φ i and ψ i . Introducing x = [ c g     c f ] T and using (3) we obtain
[ φ 1 ( x 1 ) + ψ 1 ( x 2 ) Δ s m ( ν 1 ) φ n ( x 1 ) + ψ n ( x 2 ) Δ s m ( ν n ) ] : = F ( x ) = 0 ,   n 2
with Δ s m being the measured peak heights and n being the number of characteristic wavenumbers which shall be evaluated. If n > 2 holds, a least squares solution [38] is pursued given by
x * = max   ( arg min x F ( x ) T F ( x ) , 0 ) .
If the calculations should yield negative concentrations, they are automatically set to zero, since they are from a physical point of view infeasible.
The developed algorithm allows for a precise estimation of the concentrations for glucose and fructose. Importantly, the results of further measurements can be used as additional feedback to adapt the basis functions in order to improve the process quality and accuracy on the fly.

3. Results

Three different SEIRA measurement cycles are recorded to test and evaluate the proposed procedure. The pure glucose and fructose solutions are used to train our model, that is, determine the basis functions as well as the respective peak contributions to each of the characteristic spectral signatures as a function of the absolute concentration. Afterwards, we can analyze the mixed solutions which we assume as the solutions of “unknown” concentration in order to test the accuracy of our evaluation procedure.
Within each cycle, eleven sets of dilutions are examined. In the glucose cycle concentrations of 10, 25, and 50 g/l glucose are dissolved in water. Before, after, and between each aqueous solution, the flow cell of the sensor is flushed with pure water. On the one hand, this shall remove all residuals of the monosaccharide which could spoil subsequent measurements. On the other hand, a reference measurement for comparison to the aqueous solutions is recorded [28]. The order of execution is depicted in Figure 3 in panels (a) to (c). Within each of the eleven sets thirty individual spectra are taken. This not only allows for checking repeatability but also averaging over these thirty samples to give a more accurate and smoother representative measurement for further evaluation. A brief statistical analysis is provided in the Appendix A (Figure A1 and Figure A2). In a second measurement cycle the order of execution is repeated using fructose instead of glucose. During a last cycle mixed aqueous solutions containing both specimens with different varying concentrations are examined. All measurement data and all MATLAB related files for their evaluation are provided in the Supplementary Materials for the interested reader.
Exemplary SEIRA reflectance spectra for all three sets are shown in the Figure 3 panels (d) to (f). We can clearly identify vibrational signatures characteristic for glucose at ν = 1034   cm 1 and ν = 1078   cm 1 (light green bars), whereas for fructose the vibrational signatures occur around ν = 1063   cm 1 and ν = 1080   cm 1 (light blue bars). The identified characteristic wavenumbers are in accordance with the ones reported in literature [39]. These vibrations can be identified as stretching vibrations of C-C and C-O bonds of the glucose ring or the fructose ring respectively. In case of the mixed solutions, we clearly observe the vibrational signatures of both molecular specimens. Examining the spectra closely, it is also obvious that the modulation depth is a function of the concentration. A higher respective concentration leads to a more pronounced vibrational feature. From these raw spectra it is also obvious that the mere presence of the molecular species can be easily and straightforwardly determined whereas the absolute concentration can only be estimated at best. This shortcoming we want to address next.

3.1. Pre-Processing the SEIRA Measurement

In order to extract the vibrational features of the present molecular specimen and determine its concentration we apply a slightly modified version of the ALSS algorithm [37] to the SEIRA data. The original version requires two weighting factors to be chosen. The first parameter λ determines how smoothly the reconstructed baseline is shaped. The second parameter p is responsible for the allowed asymmetry of the line shape. The algorithm is modified such that it allows to exclude parts of the fitting range within a certain neighborhood of the identified or known characteristic wavenumbers of the molecular species in the SEIRA measurement as the line is strongly deformed here. We performed a sensitivity analysis for the ALSS parameters which revealed that the absolute peak values of the vibration features in the baseline corrected spectra strongly depend on the chosen parameters. We have made different choices for the parameters λ and p , leading to different baselines. Panel (a) of Figure 4 shows a SEIRA measurement for a 50 g/l glucose solution and panel (b) a zoom into the spectral region of interest. The green curves in both panels represent different possible base lines. All of these baselines could be identified as reasonable but clearly different baselines, as is particularly obvious from the full reflectance spectrum and baselines depicted in panel (a). Panels (c) and (d) of Figure 4 display the baseline corrected SEIRA spectrum for the different baselines. The spectral position of the vibrational features is well retained but the absolute values, which are related to the actual concentration, are lost. For all the spectra it is possible to determine the presence of glucose but the concentration cannot be determined. These results clearly show that the choice of the baseline is of utmost importance and that the baseline must be chosen for all measurements simultaneously and consistently in order to retain quantitative information. Hence, the parameters are kept constant and set to λ = 500 , p = 0.99 , and ν e x c l . [ 1000 ,   1100 ] for comparability between all data sets.
Using these baseline correction values, we obtained the pure vibrational spectra. Figure 5 displays exemplary spectra for all three cases, that is, the pure glucose solutions in panel (a), the pure fructose solutions in panel (b), and the mixed solutions in panel (c). These spectra clearly exhibit an increasing modulation depth for increasing concentrations, as expected. The vibrational spectra of the mixed glucose and fructose solutions are rather complicated to interpret as they contain all the vibrational modes of glucose and fructose. The characteristic wavenumbers of glucose are determined to be at 1034 cm−1 and 1078 cm−1, whereas we obtain 1063 cm−1 and 1080 cm−1 for fructose. Since the second characteristic wavenumber of glucose and fructose are spectrally close together, we choose to only evaluate the peaks at ν = [ 1034 ,   1063 ,   1078 ] T   cm 1 in the following.
The next question arising is related to the relative weight of each component in a certain peak, that is, the question whether or not each peak can be decomposed in its relative composition. Only if this is possible, the absolute concentrations, particularly in the mixed solutions, can be determined. For this purpose, we apply the proposed algorithm to the preprocessed data.

3.2. Basis Function Approximation

When evaluating the extracted peaks of pure water measurements we obtain good agreement with the expected constant values for the glucose and fructose cycle which are depicted in panels (a) to (c) of Figure 6. The peaks obtained from the mixed cycle measurements deviate more around their average. A possible explanation could be that not all residues from previously measured aqueous solutions have been flushed out of the flow cell or are still sticking to the sensor surface. These deviations behave quite similar for all three examined wavenumbers. However, we also observe variations of the calculated mean values k i when comparing the different measurement cycles. In particular, this is visible at 1034 cm−1 in panel (a) and at 1063 cm−1 in panel (b). At this point we are not sure about the origin of these deviations. In order to include this influence, we decided to design the analysis routine adapted to the sensor data. This means that for each measurement cycle a calibration is carried out to identify the constants k i of Equations (1) and (2), characterizing the influence of pure water.
Subsequently, we analyze the SEIRA data for aqueous solutions containing only one monosaccharide. In a first step, we subtract the influence from pure water following Equations (1) and (2). Then, we have to choose suitable basis functions which are able to fit the data points best. Polynomials provided the optimum results. Starting from a simple linear function, the degree of the polynomials was tested up to an order of three. The polynomial basis functions φ i and ψ i are parametrized to match the data points best in a least squares sense. An exemplary combination of fitted basis functions is depicted in Figure 7. The shape differs clearly at each of the three characteristic wavenumbers. Panel (a) at ν 1 = 1034   cm 1 demonstrates the dominant influence of the glucose level as expected. However, interestingly, the influence of fructose is not negligible. At ν 2 = 1063   cm 1 it is just the opposite as displayed in panel (b). In the vicinity of ν 3 = 1078   cm 1 we expect resonances for both specimens visualized in panel (c). Overall, the basis functions of glucose and fructose can be matched quite well with quadratic polynomials.

3.3. Validation of Algorithm for Quantitative Concentration Estimate

Having identified suitable basis functions φ i and ψ i , the system of nonlinear Equation (4) can be solved in a least squares sense. For each of the thirty three measurement sets the concentrations are calculated using Equation (5) and compared to the expected value. For different combinations of polynomials as basis functions a comparison of the resulting RMS errors and mean errors are listed in Table 1. The best results are achieved with φ i and ψ i being quadratic polynomials with regards to the error mean value. The estimated concentrations are depicted in Figure 8 and exhibit very good agreement with the set values (Figure 3). Estimation of the concentrations of the glucose and fructose cycle, shown in panel (a) and (b), are expected to match well, since they have been used as training data for the algorithm. However, the algorithm also proves to be successful using the validation data (mixed cycle) depicted in panel (c). Overall, analyzing the concentration estimate for all three cycles a mean deviation of 0.71 g/l is achieved with the proposed procedure. With regards to the highest employed concentration of 60 g/l this mean error is below 1.2%.

4. Discussion

The results demonstrate that a concentration of monosaccharides such as glucose and fructose can be reliably and quantitatively sensed down to at least 10 g/l. The method not only works for aqueous solutions with pure glucose or pure fructose but also for mixed solutions containing both specimens. The non-invasive sensor principle relies on surface enhanced infrared absorption (SEIRA) spectroscopy of glucose and fructose in the mid-infrared spectral range. After preprocessing the SEIRA data with a baseline correction method, a new and adaptive algorithm is employed to reconstruct the concentration levels of glucose and fructose within aqueous solutions. Being able to choose arbitrarily shaped basis functions for any number of characteristic wavenumbers renders this algorithm highly flexible. The extension to incorporate additional monosaccharides is straightforward. Moreover, the structure of the algorithm facilitates its application to machine learning. Adapting or learning the basis functions and their parametrization could be achieved in an “online fashion”. This could be especially interesting for a sensor which senses glucose concentrations in the human ocular fluid or interstitial fluid that can be directly correlated to the blood glucose concentrations [7]. However, the measurement setup has to be improved in order to increase the sensitivity further to achieve sensitivity below 10 g/l, eventually down to physiological concentrations of mg/l.
In order to explore the limitations of the proposed methods regarding robustness and generality further experiments are required. Other monosaccharides, as well as maltose as a biologically relevant double sugar, should be added as analytes, not only to identify their fingerprint resonances but also to validate whether the method is robust enough to accurately predict the concentrations of glucose and fructose even in more complex environments. Furthermore, measurements with aqueous solutions containing glucose or fructose with concentrations below 10 g/l have to be taken in future research.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-8220/19/14/3053/s1.

Author Contributions

B.S. performed modeling and mathematical description of the data. L.K. und M.H. fabricated the samples and took the experimental data. H.G. and C.T. conceived the project idea and supervised the project. All authors contributed to the discussion and interpretation of the data as well as the writing of the manuscript.

Funding

L.K., M.H., and H.G. gratefully acknowledge financial support by the ERC Advanced Grant (COMPLEXPLAS), Deutsche Forschungsgemeinschaft (SPP1839), Baden-Württemberg Stiftung (PROTEINSENS), MWK Baden-Württemberg (IQST, ZAQuant), and BMBF (03IS2101E).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

A brief statistical analysis is done to proof the repeatability of the measurements. Within each measurement set thirty individual spectra are taken. For an exemplary set, namely the one with 50 g/l glucose and 50 g/l fructose dissolved in water, the individual spectra are depicted in panel (a) of Figure A1. Within the interested fingerprint region all spectra exhibit a very similar qualitative behavior (panel (b)). For the shown set, which is a worst case example, only one spectrum differs significantly from the others. Since we normalize the reflectance spectra with the baseline correction method afterwards in our data analysis procedure this offset becomes negligible anyways. Within the fingerprint region around 1000 to 1100 cm−1 the mean value of the variance (panel (d)) and the standard deviation (panel (f)) are below 10−5 and 4 × 10−3 respectively. Only for wavenumbers below 1000 cm−1 or above 6000 cm−1 the variance and standard deviation increase which is shown in panels (c) and (e). These inaccuracies in the border regions can be explained with the sensor design. It is tuned to provide best measurement results for the fingerprint region.
Repeating this statistical analysis for all measurement sets provides a base to judge the quality of the measurements. The results are depicted in Figure A2. In panels (a) to (c) the concentrations of glucose and fructose are visualized. The calculated mean values of the variances and the standard deviations within the fingerprint region are illustrated in panels (d) to (f) and panels (g) to (i) respectively. Only the first water measurement of the glucose and fructose cycle exhibit a slightly higher variance compared to the other sets of the corresponding cycle. Measurement sets six and eleven of the mixed cycle display the highest variance. For both sets, an individual spectrum had an almost constant offset in the reflectance as visualized in panel (b) of Figure A1. However, this offset becomes negligible due to the averaging over the thirty spectra and due to the normalization through the baseline correction method.
Figure A1. Statistical analysis of exemplary measurement data. The shown reflectance (a,b) results from a mixed solution with 50 g/l glucose and 50 g/l fructose. The measurement data consist of thirty individual spectra. The variance and the standard deviation between those spectra are calculated for each wavenumber and depicted in panels (c,d) and (e,f) respectively. The chosen measurement set reveals one of the largest variances because the reflectance of one spectra exhibits an almost constant offset. However, the qualitative features are preserved.
Figure A1. Statistical analysis of exemplary measurement data. The shown reflectance (a,b) results from a mixed solution with 50 g/l glucose and 50 g/l fructose. The measurement data consist of thirty individual spectra. The variance and the standard deviation between those spectra are calculated for each wavenumber and depicted in panels (c,d) and (e,f) respectively. The chosen measurement set reveals one of the largest variances because the reflectance of one spectra exhibits an almost constant offset. However, the qualitative features are preserved.
Sensors 19 03053 g0a1
Figure A2. Statistical analysis of the measurement sets for each of the three cycles. The glucose/fructose concentration for each set is depicted in panels (ac). The corresponding mean value of the variance within the fingerprint region ranging from 1000–1100 cm−1 (compare to panel (d) of Figure A1) is visualized in panels (df). The corresponding mean value of the standard deviation is depicted in panels (gi).
Figure A2. Statistical analysis of the measurement sets for each of the three cycles. The glucose/fructose concentration for each set is depicted in panels (ac). The corresponding mean value of the variance within the fingerprint region ranging from 1000–1100 cm−1 (compare to panel (d) of Figure A1) is visualized in panels (df). The corresponding mean value of the standard deviation is depicted in panels (gi).
Sensors 19 03053 g0a2

References

  1. Alberti, K.G.; Zimmet, P.Z. Definition. Diagnosis and Classification of Diabetes Mellitus and Its Complications. Part 1: Diagnosis and Classification of Diabetes Mellitus. Provisional Report of a WHO Consultation. Diabet. Med. 1998, 15, 539–553. [Google Scholar] [CrossRef]
  2. American Diabetes Association. Diagnosis and Classification of Diabetes Mellitus. Diabetes Care 2014, 37, 81–90. [Google Scholar] [CrossRef] [PubMed]
  3. Ritz, E.; Rychlík, I.; Locatelli, F.; Halimi, S. End-Stage Renal Failure in Type 2 Diabetes: A Medical Catastrophe of Worldwide Dimensions. Am. J. Kidney Dis. 1999, 34, 795–808. [Google Scholar] [CrossRef]
  4. Eberle, C.; Amet, C. Real-Time State Estimation and Long-Term Model Adaptation: A Two-Sided Approach toward Personalized Diagnosis of Glucose and Insulin Levels. J. Diabetes Sci. Technol. 2012, 6, 1148–1158. [Google Scholar] [CrossRef] [Green Version]
  5. Alexeev, V.L.; Das, S.; Finegold, D.N.; Asher, S.A. Photonic Crystal Glucose-Sensing Material for Noninvasive Monitoring of Glucose in Tear Fluid. Clin. Chem. 2004, 50, 2353–2360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pleitez, M.A.; Hertzberg, O.; Bauer, A.; Seeger, M.; Lieblein, T.; Lilienfeld-Toal, H.V.; Mäntele, W. Photothermal Deflectometry Enhanced by Total Internal Reflection Enables Non-Invasive Glucose Monitoring in Human Epidermis. Analyst 2015, 140, 483–488. [Google Scholar] [CrossRef]
  7. Bruen, D.; Delaney, C.; Florea, L.; Diamond, D. Glucose Sensing for Diabetes Monitoring: Recent Developments. Sensors 2017, 17, 1866. [Google Scholar] [CrossRef] [PubMed]
  8. Halas, N.J.; Lal, S.; Chang, W.-S.; Link, S.; Nordlander, P. Plasmons in strongly coupled metallic nanostructures. Chem. Rev. 2011, 111, 3913–3961. [Google Scholar] [CrossRef] [PubMed]
  9. Maier, S.A. Plasmonics: Fundamentals and Applications; Springer: Berlin, Germany, 2007. [Google Scholar]
  10. Shih, W.C.; Santos, G.M.; Zhao, F.; Zenasni, O.; Arnob, M.M.P. Simultaneous Chemical and Refractive Index Sensing in the 1–2.5 Μm near-Infrared Wavelength Range on Nanoporous Gold Disks. Nano Lett. 2016, 16, 4641–4647. [Google Scholar] [CrossRef] [PubMed]
  11. Liu, N.; Mesch, M.; Weiss, T.; Hentschel, M.; Giessen, H. Infrared perfect absorber and its application as plasmonic sensor. Nano Lett. 2010, 10, 2342–2348. [Google Scholar] [CrossRef]
  12. La Spada, L. Metasurfaces for Advanced Sensing and Diagnostics. Sensors 2019, 19, 355. [Google Scholar] [CrossRef]
  13. Liu, N.; Tang, M.L.; Hentschel, M.; Giessen, H.; Alivisatos, A.P. Nanoantenna-enhanced gas sensing in a single tailored nanofocus. Nat. Mater. 2011, 10, 631–636. [Google Scholar] [CrossRef]
  14. Brolo, A.G. Plasmonics for future biosensors. Nat. Photonics 2012, 6, 709–713. [Google Scholar] [CrossRef]
  15. Mesch, M.; Zhang, C.; Braun, P.V.; Giessen, H. Functionalized hydrogel on plasmonic nanoantennas for noninvasive glucose sensing. ACS Photonics 2015, 2, 475–480. [Google Scholar] [CrossRef]
  16. Lee, Y.-J.; Pruzinsky, S.A.; Braun, P. V Glucose-Sensitive Inverse Opal Hydrogels: Analysis of Optical Diffraction Response. Langmuir 2004, 20, 3096–3106. [Google Scholar] [CrossRef]
  17. Griffiths, P.R.; de Haseth, J.A. Fourier Transform Infrared Spectrometry; John Wiley & Sons: Hoboken, NJ, USA, 2007. [Google Scholar]
  18. Neubrech, F.; Pucci, A.; Cornelius, T.; Karim, S.; García-Etxarri, A.; Aizpurua, J. Resonant plasmonic and vibrational coupling in a tailored nanoantenna for infrared detection. Phys. Rev. Lett. 2008, 101, 157403. [Google Scholar] [CrossRef] [PubMed]
  19. Neubrech, F.; Huck, C.; Weber, K.; Pucci, A.; Giessen, H. Surface-enhanced infrared spectroscopy using resonant nanoantennas. Chem. Rev. 2017, 117, 5110–5145. [Google Scholar] [CrossRef]
  20. Rodrigo, D.; Limaj, O.; Janner, D.; Etezadi, D.; Garcia de Abajo, F.J.; Pruneri, V.; Altug, H. Mid-infrared plasmonic biosensing with graphene. Science 2015, 349, 165–168. [Google Scholar] [CrossRef] [Green Version]
  21. Adato, R.; Altug, H. In-situ ultra-sensitive infrared absorption spectroscopy of biomolecule interactions in real time with plasmonic nanoantennas. Nat. Commun. 2013, 4, 2154. [Google Scholar] [CrossRef]
  22. Le, F.; Brandl, D.W.; Urzhumov, Y.A.; Wang, H.; Kundu, J.; Halas, N.J.; Aizpurua, J.; Nordlander, P. Metallic nanoparticle arrays: A common substrate for both surface-enhanced Raman scattering and surface-enhanced infrared absorption. ACS Nano 2008, 2, 707–718. [Google Scholar] [CrossRef]
  23. Wang, T.; Nguyen, V.H.; Buchenauer, A.; Schnakenberg, U.; Taubner, T. Enhanced Infrared Spectroscopy with Gold Strip Gratings. Opt. Express 2013, 21, 9005. [Google Scholar] [CrossRef] [PubMed]
  24. Semenyshyn, R.; Hentschel, M.; Stanglmair, C.; Teutsch, T.; Tarin, C.; Pacholski, C.; Giessen, H.; Neubrech, F. In Vitro Monitoring Conformational Changes of Polypeptide Monolayers Using Infrared Plasmonic Nanoantennas. Nano Lett. 2018, 19, 1–7. [Google Scholar] [CrossRef]
  25. Etezadi, D.; Warner, J.B.; Lashuel, H.A.; Altug, H. Real-time in situ secondary structure analysis of protein monolayer with mid-infrared plasmonic nanoantennas. ACS Sens. 2018, 3, 1109–1117. [Google Scholar] [CrossRef]
  26. Neubrech, F.; Beck, S.; Glaser, T.; Hentschel, M.; Giessen, H.; Pucci, A. Spatial extent of plasmonic enhancement of vibrational signals in the infrared. ACS Nano 2014, 8, 6250–6258. [Google Scholar] [CrossRef] [PubMed]
  27. Dregely, D.; Neubrech, F.; Duan, H.; Vogelgesang, R.; Giessen, H. Vibrational near-field mapping of planar and buried three-dimensional plasmonic nanostructures. Nat. Commun. 2013, 4, 2237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kühner, L.; Semenyshyn, R.; Hentschel, M.; Neubrech, F.; Tarin, C.; Giessen, H. Vibrational Sensing Using Infrared Nanoantennas: Toward the Noninvasive Quantitation of Physiological Levels of Glucose and Fructose. ACS Sens. 2019. [Google Scholar] [CrossRef]
  29. Isensee, K.; Müller, N.; Pucci, A.; Petrich, W. Towards a Quantum Cascade Laser-Based Implant for the Continuous Monitoring of Glucose. Analyst 2018, 143, 6025–6036. [Google Scholar] [CrossRef] [PubMed]
  30. La Spada, L.; Vegni, L. Electromagnetic Nanoparticles for Sensing and Medical Diagnostic Applications. Materials 2018, 11, 603. [Google Scholar] [CrossRef] [PubMed]
  31. Liberal, I.; Engheta, N. Near-zero refractive index photonics. Nat. Photonics 2017, 11, 149–158. [Google Scholar] [CrossRef]
  32. La Spada, L.; Spooner, C.; Haq, S.; Hao, Y. Curvilinear MetaSurfaces for Surface Wave Manipulation. Sci. Rep. 2019, 9, 3107. [Google Scholar] [CrossRef]
  33. Soavi, G.; Wang, G.; Rostami, H.; Purdie, D.G.; De Fazio, D.; Ma, T.; Luo, B.; Wang, J.; Ott, A.K.; Yoon, D.; et al. Broadband, electrically tunable third-harmonic generation in graphene. Nat. Nanotechnol. 2018, 13, 583–588. [Google Scholar] [CrossRef] [Green Version]
  34. La Spada, L.; Vegni, L. Metamaterial-based wideband electromagnetic wave absorber. Opt. Express 2016, 24, 5763. [Google Scholar] [CrossRef] [PubMed]
  35. Adato, R.; Yanik, A.A.; Wu, C.-H.; Shvets, G.; Altug, H. Radiative engineering of plasmon lifetimes in embedded nanoantenna arrays. Opt. Express 2010, 18, 4526. [Google Scholar] [CrossRef] [PubMed]
  36. Bagheri, S.; Weber, K.; Gissibl, T.; Weiss, T.; Neubrech, F.; Giessen, H. Fabrication of Square-Centimeter Plasmonic Nanoantenna Arrays by Femtosecond Direct Laser Writing Lithography: Effects of Collective Excitations on SEIRA Enhancement. ACS Photonics 2015, 2, 779–786. [Google Scholar] [CrossRef]
  37. Eilers, P.H.C.; Boelens, H.F.M. Baseline Correction with Asymmetric Least Squares Smoothing. Leiden Univ. Med. Cent. Rep. 2005, 1, 5. [Google Scholar]
  38. Shardt, Y.A.W. Statistics for Chemical and Process Engineers; Springer International Publishing: Berlin, Germany, 2015. [Google Scholar]
  39. Ibrahim, M.; Alaam, M.; El-Haes, H.; Jalbout, A.F.; de Leon, A. Analysis of the structure and vibrational spectra of glucose and fructose. Eclét. Quim. 2006, 31, 15–21. [Google Scholar] [CrossRef]
Figure 1. Measurement principle and sensor design. Aqueous pure glucose and fructose solutions as well as mixed solutions are flushed over plasmonic nanoatennas resonant with the characteristic mid-infrared vibrational modes of these monosaccharaides. The tailored flow cell allows to simultaneously acquire reflectance spectra, as sketched on the right hand side. The reflectance peaks of the antennas are modified by the vibrational signatures of the molecular species present in the analyte. From their spectral position the molecular species can be determined, from their modulation depth, in principle, the concentration, given appropriate evaluation and calibration, which is the key result reported here.
Figure 1. Measurement principle and sensor design. Aqueous pure glucose and fructose solutions as well as mixed solutions are flushed over plasmonic nanoatennas resonant with the characteristic mid-infrared vibrational modes of these monosaccharaides. The tailored flow cell allows to simultaneously acquire reflectance spectra, as sketched on the right hand side. The reflectance peaks of the antennas are modified by the vibrational signatures of the molecular species present in the analyte. From their spectral position the molecular species can be determined, from their modulation depth, in principle, the concentration, given appropriate evaluation and calibration, which is the key result reported here.
Sensors 19 03053 g001
Figure 2. Schematic description of the developed adaptive process for evaluation of surface enhanced infrared absorption (SEIRA) measurement data to quantitatively estimate glucose and fructose concentrations in aqueous solutions.
Figure 2. Schematic description of the developed adaptive process for evaluation of surface enhanced infrared absorption (SEIRA) measurement data to quantitatively estimate glucose and fructose concentrations in aqueous solutions.
Sensors 19 03053 g002
Figure 3. Employed measurement cycles with different concentrations of glucose (a), fructose (b), or both (c). Panels (d) to (f) show the corresponding relative reflectance of each of the eleven SEIRA measurement sets with the expected fingerprint region for molecular vibrations of glucose and fructose highlighted in light green or light blue, respectively.
Figure 3. Employed measurement cycles with different concentrations of glucose (a), fructose (b), or both (c). Panels (d) to (f) show the corresponding relative reflectance of each of the eleven SEIRA measurement sets with the expected fingerprint region for molecular vibrations of glucose and fructose highlighted in light green or light blue, respectively.
Sensors 19 03053 g003
Figure 4. (a) Relative reflectance measured for 50 g/l aqueous glucose together with a set of baselines for varying parameters λ ; (b) Zoom in on the relevant range of wavenumbers at which vibrational resonances are expected to appear for glucose and fructose; (c,d) Baseline corrected signals and the influence of λ on the vibrational features.
Figure 4. (a) Relative reflectance measured for 50 g/l aqueous glucose together with a set of baselines for varying parameters λ ; (b) Zoom in on the relevant range of wavenumbers at which vibrational resonances are expected to appear for glucose and fructose; (c,d) Baseline corrected signals and the influence of λ on the vibrational features.
Sensors 19 03053 g004
Figure 5. Relative reflectance of baseline corrected signals with λ = 500 and p = 0.99 for glucose cycle (a), fructose cycle (b), and mixed cycle (c).
Figure 5. Relative reflectance of baseline corrected signals with λ = 500 and p = 0.99 for glucose cycle (a), fructose cycle (b), and mixed cycle (c).
Sensors 19 03053 g005
Figure 6. Signal peaks Δ s evaluated at the three characteristic wavenumbers ν 1 = 1034   cm 1 (a), ν 2 = 1063   cm 1 (b), and ν 3 = 1078   cm 1 (c) originating from pure water for each of the three measurement cycles (gc = glucose cycle, fc = fructose cycle, and mc = mixed cycle).
Figure 6. Signal peaks Δ s evaluated at the three characteristic wavenumbers ν 1 = 1034   cm 1 (a), ν 2 = 1063   cm 1 (b), and ν 3 = 1078   cm 1 (c) originating from pure water for each of the three measurement cycles (gc = glucose cycle, fc = fructose cycle, and mc = mixed cycle).
Sensors 19 03053 g006
Figure 7. Signal peaks Δ s evaluated at the three characteristic wavenumbers ν 1 = 1034   cm 1 (a), ν 2 = 1063   cm 1 . (b), and ν 3 = 1078   cm 1 (c) originating from pure glucose or pure fructose, respectively. The content of pure water is already subtracted using the corresponding average value. The displayed basis functions φ i and ψ i are quadratic polynomials.
Figure 7. Signal peaks Δ s evaluated at the three characteristic wavenumbers ν 1 = 1034   cm 1 (a), ν 2 = 1063   cm 1 . (b), and ν 3 = 1078   cm 1 (c) originating from pure glucose or pure fructose, respectively. The content of pure water is already subtracted using the corresponding average value. The displayed basis functions φ i and ψ i are quadratic polynomials.
Sensors 19 03053 g007
Figure 8. Estimated levels of concentration for the glucose cycle (a), fructose cycle (b), and mixed cycle (c) using cycle-specific constants k i , quadratic polynomials φ i and linear basis functions ψ i .
Figure 8. Estimated levels of concentration for the glucose cycle (a), fructose cycle (b), and mixed cycle (c) using cycle-specific constants k i , quadratic polynomials φ i and linear basis functions ψ i .
Sensors 19 03053 g008
Table 1. Comparison of RMS errors, mean errors, and maximum errors for the estimated concentrations resulting from different combinations of polynomial basis functions φ i and ψ i with order n g l u c o s e or n f r u c t o s e respectively.
Table 1. Comparison of RMS errors, mean errors, and maximum errors for the estimated concentrations resulting from different combinations of polynomial basis functions φ i and ψ i with order n g l u c o s e or n f r u c t o s e respectively.
n g l u c o s e n f r u c t o s e e r m s   [ g / l ] e m e a n   [ g / l ] e m a x   [ g / l ]
111.270.813.86
121.270.773.98
131.890.9510.02
211.200.734.13
221.230.714.29
231.670.888.36
311.510.808.03
321.530.798.24
331.690.906.50

Share and Cite

MDPI and ACS Style

Schuler, B.; Kühner, L.; Hentschel, M.; Giessen, H.; Tarín, C. Adaptive Method for Quantitative Estimation of Glucose and Fructose Concentrations in Aqueous Solutions Based on Infrared Nanoantenna Optics. Sensors 2019, 19, 3053. https://doi.org/10.3390/s19143053

AMA Style

Schuler B, Kühner L, Hentschel M, Giessen H, Tarín C. Adaptive Method for Quantitative Estimation of Glucose and Fructose Concentrations in Aqueous Solutions Based on Infrared Nanoantenna Optics. Sensors. 2019; 19(14):3053. https://doi.org/10.3390/s19143053

Chicago/Turabian Style

Schuler, Benjamin, Lucca Kühner, Mario Hentschel, Harald Giessen, and Cristina Tarín. 2019. "Adaptive Method for Quantitative Estimation of Glucose and Fructose Concentrations in Aqueous Solutions Based on Infrared Nanoantenna Optics" Sensors 19, no. 14: 3053. https://doi.org/10.3390/s19143053

APA Style

Schuler, B., Kühner, L., Hentschel, M., Giessen, H., & Tarín, C. (2019). Adaptive Method for Quantitative Estimation of Glucose and Fructose Concentrations in Aqueous Solutions Based on Infrared Nanoantenna Optics. Sensors, 19(14), 3053. https://doi.org/10.3390/s19143053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop