Diabetes mellitus is a chronic metabolic disorder that affects over 347 million people worldwide . Although there are many causes for diabetes, the primary symptom of this disease is excessive glucose levels in the blood. The long-term complications associated with diabetes include cardiovascular disease, chronic renal failure, and diabetic retinopathy, all of which can be reduced significantly by controlling blood glucose levels using a blood glucose meter (BGM). Traditional BGMs employ bio-electrochemical schemes that use glucose oxidase (GOx) or glucose dehydrogenase (GDH) immobilized on the surface of a disposable electrode . These monitoring schemes are non-toxic and have high selectivity; however, the detection limit of commercial BGMs is ∼2.8 mm (50 mg/dL) , which prevents their use in fluids such as tears or saliva. These fluids, while not requiring invasive extraction are characterized by glucose concentrations that are typically two orders of magnitude lower than what is found in blood.
Over the past decade, extensive research has focused on optical, rather than electrochemical methods, with the goal of reducing the limit of detection and achieving non-invasive sensing [4–12]. Optical non-invasive detection schemes assume a direct correlation between the glucose levels found in blood and other bodily fluids; this correlation is well established for interstitial fluids  and saliva [14, 15]. Saliva is a particularly attractive candidate for non-invasive glucose sensing since it is almost always present and easy to acquire (unlike tears). These optical biosensors can be further improved by employing metal nanoparticles or planar metal films with nano-corrugations in order to favorably manipulate the incident light beam [16–21]. Example techniques include surface plasmon resonance (SPR) [22, 23], grating couplers , and plasmonic interferometers [25, 26]. In all these techniques, incident light is coupled to surface plasmon polaritons (SPPs) that are electromagnetic waves supported by electron density fluctuations in noble metals (e.g., silver and gold) . SPPs are generated from localized sources (e.g., grooves) and propagate away from the source and along the metal/dielectric interface to interfere with each other and produce local field intensity minima and maxima. Since SPPs are confined at the metal surface, they are highly sensitive to changes in the refractive index of the dielectric material, a property useful in developing ultra-sensitive, compact and real-time biosensors.
Herein we demonstrate a sensing scheme denominated “plasmonic cuvette” (Figure 1) obtained by coupling the spectroscopic fingerprinting capabilities of a groove-slit-groove (GSG) plasmonic interferometer (Figure 1G) to an Amplex–red/Glucose–oxidase/Glucose (AR/GOx/Glucose) enzyme-driven dye assay (Figure 1A–C) that adds selectivity for glucose detection, and further enhances the device sensitivity by one order of magnitude. The resulting sensor offers real-time sensitivity toward glucose in extremely small sensing volumes (i.e., ≤12 pL), and exhibits glucose selectivity in complex mixtures such as a 50 mm sodium phosphate buffer solution and “artificial” saliva, over the physiological range of glucose concentration in saliva (20–240 μm) . These results demonstrate the viability of measuring the concentration of glucose in saliva, which is a complex mixture of proteins, salts, and urea . More details about the experimental methods can be found in Appendix A.
The plasmonic cuvette (shown in Figure 1G) offers several advantages over a traditional UV–Visible (UV–Vis) cuvette (Figure 1E). First, the sensed volume in a plasmonic cuvette is 12 pL (compared to ≥240 μL for the UV–Vis cuvette), thus reducing the amount of sample and reagents required by seven orders of magnitude. Second, the plasmonic cuvette can be integrated with a microfluidic channel to achieve more uniform mixing than what is normally possible in a UV–Vis cuvette. Third, the plasmonic cuvette has a very short effective path length (<20 μm) compared to a UV-Vis cuvette (∼1 cm). The shorter path length of a plasmonic cuvette allows the response of the system to remain linear over a wider range of analyte concentrations and avoids the necessity to either dilute the solution or reduce the path length of the UV-Vis cuvette at higher analyte concentrations in order to have a measurable transmittance.
Figure 1A provides an overview of the assay reactions that, starting from one molecule of glucose, end up generating one molecule of optically active resorufin, which can be identified and quantified more easily using optical methods, such as plasmonic interferometry. Typically, the fluorescence intensity from resorufin is used to quantify the initial glucose concentration. However, the emission efficiency strongly depends on the pH of the solution, as well as on the chemical components dissolved in solution, thus affecting the quantum yield of the emitter, and making it difficult to directly correlate the measured fluorescence intensity with the initial glucose concentration. Here, we propose the use of the absorption (rather than emission) properties of resorufin to reliably determine the glucose concentration. We characterize the absorption cross section of the emitter in solution, and calibrate the assay for accurate glucose detection. Plasmonic interferometers are used to amplify and monitor the absorption changes in solution as well as quantify the glucose concentration in very small sample volumes, with high sensitivity and selectivity.
Details of each enzymatic reaction occurring in the assay are shown in Figure 1B–D. Reaction 1 (Figure 1B) describes the oxidation of β-D-glucose to D-gluconolactone by oxygen catalyzed by GOx to yield H2 O2 as a product . The H2 O2 generated in Reaction 1 is utilized by horseradish peroxidase (HRP) to oxidize Amplex red (AR), a colorless, non-fluorescent compound, to resorufin, a red, fluorescent compound (Reaction 2, Figure 1C) [31–33]. Reaction 3 (Figure 1D) is a much slower, H2 O2-limited side reaction that further oxidizes resorufin into an optically inactive product (OIP), thought to be a complex polymer [30, 33, 34]. Reactions 1 and 2 are commonly used to detect glucose in a wide variety of applications, such as determining glucose consumption by cancerous cells to better understand tumorigenesis , or quantifying starch and glucose concentrations in food . In these applications, the fluorescence properties of resorufin are preferred because of the high signal-to-noise ratio (Figure 1E). By contrast, this work employs the absorption properties of resorufin (Figure 1F) in order to determine a quantifiable change in the optical transmission response of the plasmonic cuvette, which can be correlated to the glucose concentration initially present in solution. The foundation of the plasmonic cuvette, a groove-slit-groove (GSG) plasmonic interferometer, is shown in Figure 1G . When light of wavelength λ is incident upon a groove, an SPP wave is generated and propagates through the sample solution toward the center slit, which is a distance p away from the groove (p is also known as the arm length of the interferometer). The counter propagating SPP waves originating from the two grooves interfere constructively or destructively at the slit location depending on p, λ, and the refractive index n of the dielectric resting on top of the metal. The output transmission spectrum will be a convolution of the incident light and the net SPP fields at the slit location. Changes in the intensity or shifts in the wavelength of the output spectrum provide information about the concentration and types of analytes present (Figure 1H). We also report a detailed study of the enzyme kinetics involved in the AR/GOx/Glucose assay: the effective rate constants are found and then used to tune the reaction time and range of applicability of the assay to cover the physiological range of glucose in saliva, in real-time.
2 Results and discussion
2.1 Enhancing the sensitivity of plasmonic interferometers
Both the plasmonic and UV-Vis cuvettes rely on the same figure of merit to quantify sensitivity, which itself depends on the ability of the solution to absorb incident light. The relevant figure of merit is the relative intensity change (ΔI/I0), defined as
where Iglucose is the light intensity transmitted through the slit of the plasmonic interferometer at a non-zero glucose concentration, while I0 corresponds to the absence of glucose. The light intensity transmitted through the slit of a GSG interferometer, normalized to the light intensity transmitted through a single slit (SS) is :
where p1 and p2 are the lengths of the left and right arms of the GSG plasmonic interferometer, respectively, λ is the free space incident wavelength, kSPP=2πnSPP/λ and k=2πn/λ, where nSPP is the complex refractive index of the SPP given by
Figure 2A shows the measured relative intensity change defined in Eq. (1) as a function of incident wavelength, when the assay is present (red solid line) or absent (blue dashed line); the initial glucose concentration is [G]t=0=250±6 μm in both cases. The AR/GOx/Glucose assay produces 1 mol of resorufin for every mol of glucose consumed. Consequently, any change in light intensity corresponds directly to a change in glucose concentration – this figure clearly shows that the assay is effective in increasing the sensitivity of the device. In this device, two distinct optical processes are occurring: direct absorption of the incident light (characterized by the absorption cross-section curve in Figure 1F) and absorption of the groove-generated surface plasmon polaritons (SPPs). In this regard, the spectrum can be divided into three distinct regions. At shorter wavelengths (450 nm<λ<510 nm), the direct absorption of light by resorufin is low enough that the two processes can be resolved. The pronounced peaks in this region of the spectrum can be attributed to interference between counter-propagating SPP waves; the shift in peak position as well as the intensity change can be correlated to a refractive index change, as initially reported in Ref. . At wavelengths between 510 nm and 590 nm, resorufin is strongly absorbing; thus, bulk absorption dominates over plasmonic interference and absorption. The presence of this absorption band, which closely resembles the shape of the absorption spectrum of resorufin, ensures the device is selective towards glucose and can be attributed to the high specificity of GOx to glucose in the assay. At λmax=571 nm, the absolute value of the relative intensity change is at a maximum and goes from 2.8% (without assay) to 41.7% (with assay). At wavelengths above 590 nm, resorufin does not absorb light significantly, so any relative intensity change is due solely to changes in the SPPs propagation constant and interference behavior, which depends on the refractive index of the solution. This observation suggests that differences in the two spectra are a result of the enzymes and resorufin altering the refractive index of the solution (both the real and imaginary parts). In theory, any wavelength in this spectrum could be used for detection, but the best sensitivity is achieved at λmax=571 nm. Note that Figure 1H displays the same data, “with assay”, but over a range of glucose concentrations between 0 and 250 μm, demonstrating the expanded range of applicability provided by the coupling of dye chemistry with plasmonic interferometry.
Figure 2B shows the relative intensity change as a function of initial glucose concentration with and without the dye assay present in the microfluidic channel. When the assay is absent (blue circles), the calibration curve yields a sensitivity of 0.2×105%/m. When the assay is present (red squares), a linear fit of the data yields a sensitivity of 1.7×105%/m, which corresponds to an 8.5× increase (N.B.: the sensitivity is defined as the absolute value of the slope of the linear fit). The measurements were performed on the same plasmonic interferometer, but the detection wavelengths were chosen as λ=571 nm (with assay) and λ=628 nm (without assay) to maximize the sensitivity of the device in each case. For comparison, the sensitivity of our previously reported GSG plasmonic interferometer – with p1=5.70 μm and p2=9.75 μm, λ=590 nm, no assay present – yielded a sensitivity of 0.4×105%/m . The data presented here demonstrates that not only can the assay be used with different plasmonic interferometers (i.e., different arm lengths, p) to detect glucose, but that the sensitivity and selectivity of the device is significantly enhanced by the presence of the assay.
2.2 Development of an optical ruler for resorufin
Since the absorption properties of resorufin affect the output signal of the plasmonic interferometer, an optical “ruler” is first determined based on the optical properties of resorufin and the buffer composition in order to optimize the performance of the plasmonic cuvette. In order to exploit the 1:1 stoichiometry of the dye assay, the optical parameters are measured at steady-state concentration of resorufin by optical means to retrieve the initial glucose concentration in solution.
One important parameter is the absorption cross-section (σ) of resorufin, in units of cm2, defined as the effective cross-sectional area of a resorufin molecule when light is incident upon it. The absorption cross-section can be defined at all wavelengths and it can slightly depend on the nature and composition of the buffer solution. Knowledge of σ is required to determine the absolute concentration of resorufin, and therefore glucose initially present in solution. Related to σ is the extinction coefficient (ε) of resorufin, in units of m-1 cm-1, which is a measure of how strongly the optically-active molecules absorb incident light. Shown in Figure 1F is a plot of σ as a function of wavelength on the left vertical axis. The extinction coefficient (ε) of resorufin is plotted on the right axis and is related to σ by
where c and ρ are the molar concentration and the density of resorufin, respectively (see Appendix B, for a full derivation, and Appendix C for the experimental determination of σ from absorption measurements in solution). The incident wavelength chosen for this work was λmax=571 nm where σ571nm=(2.05±0.08)×10-16 cm2 and ε571nm=(5.4±0.2)×104m-1 cm-1 (compared to the literature value of ε571nm=(5.4±0.4)×104m-1 cm-1) [31, 37]. The accuracy of the ε571nm value reported here confirms the validity of this method to obtain σ and ε at other wavelengths. Figure 1F expands the range of applicable wavelengths for resorufin detection to more than just its peak absorption wavelength. This method can be generalized for deriving the absorption properties for any optically active molecule – once σ or ε has been found, the optimal wavelength can be selected to use with the sensor. Additionally, because every molecule has its own unique σ and ε wavelength dependence, the plasmonic cuvette is capable of multiplex sensing, i.e., two or more optically active biochemical analytes that are present on the chip could be identified and quantified simultaneously.
2.3 Mathematical model of AR/GOx/Glucose assay to determine the reaction constants K1, K2, and K3
The AR/GOx/Glucose assay typically requires an incubation time of >30 min for the reactions to proceed to completion . The incubation time needs to be shortened significantly if the plasmonic cuvette is to be used as a real-time sensor for glucose. Efforts have been made to characterize the factors that affect the incubation time of the coupled GOx- and HRP-dependent reactions . A study by Nakajima et al.reduced the incubation time to 4.8 s by photochemically arranging the HRP and GOx molecules in a specific pattern inside a microfluidic PDMS channel .
Because all of the reagents are simultaneously present in solution in a plasmonic cuvette sample, the best approach is to increase the rates of the three reactions (see Figure 1B–D, corresponding to rate constants K1, K2 and K3, respectively). The values for these rate constants can be obtained by fitting the change in the concentration of resorufin as the reaction progresses. Using first-order kinetics, a model for the time-dependent concentrations can be established with the following rate equations:
where G is glucose, AR is Amplex red, R is resorufin, and OIP is the optically inactive product generated in Reaction 3. This non-linear system of differential equations was solved numerically to fit the experimental data, with the initial concentrations of the six reagents as initial parameters and the rate constants as fitting parameters. The concentration of resorufin at time t, [R]t was found by measuring the transmittance (T) as a function of time t, and using the following formula (see Appendix B for a derivation):
where ε is the extinction coefficient of resorufin and L is the path length of the cuvette.
2.4 Determination of K1
To obtain the experimental concentration of resorufin as a function of time, the transmittance recorder at each time t (Figure 3A) was monitored for reaction mixtures in a 50 mm sodium phosphate buffer (pH 7.4) with 280±6 μm AR, 5.5±0.1 nm HRP and 82.5±0.7 nm GOx, and then converted into resorufin concentration vs. time using Eq. (5), as reported in Figure 3B. The concentration of glucose was varied between 0 and 210 μm. The initial concentration of Amplex red was set to 280 μm for two reasons: (i) to ensure a 1:1 stoichiometric ratio between AR and glucose while allowing the initial concentration of glucose to vary across the entire physiological range of glucose in human saliva, (ii) to prevent precipitation of AR (i.e., the solubility limit of AR is 300 μm in the aqueous buffer solutions used here) .
Although it is possible to solve for K1, K2, and K3 simultaneously, using too many fitting parameters at once (including [G]t=0) could lead to large errors. To simplify the calculations, the system was solved using a bottom-up approach: Reaction 3 was studied by itself to obtain a value for K3; using this value, Reactions 2 and 3 were run to obtain K2 (see Appendix D–F for the kinetic analysis; see Table 1 for the final values). However, when solving the full system of reactions for K1, K3 was found to be effectively zero, i.e., Reaction 3 does not occur. This assumption is supported by the fact that Reaction 3 is slower than Reaction 2 (K3<K2) and that Reaction 3 will occur only after there is a buildup of H2 O2 in the reaction mixture. However, H2 O2 is readily used up by Reaction 2, Eqs. (4C) and (4D), so that it does not contribute significantly to Reaction 3, i.e., Eqs. (4E) and (4F). Zhou et al. demonstrated that a decrease in resorufin concentration (i.e., the initiation of Reaction 3 to generate OIP) is observed only when the ratio of solution concentrations of H2 O2 to AR is 2:1 . In the full assay, the initial concentrations of resorufin and H2 O2 always start at zero. Using the model to simulate the concentration of all six reagents as a function of time, the concentration of H2 O2 never exceeds 10 μm (see Appendix G, Figure 10) when the highest initial concentration of glucose (210 μm) was used. This concentration of H2 O2 is significantly lower than that of AR in solution.
Rate constants for the reactions in the AR/GOx/Glucose assay.
|Rate constant||Enzyme concentrations||Mean value (m-1 s-1)|
|K1||HRP: 5.5 nm; GOx: 82.5 nm||7.3±0.3|
|K2||HRP: 5.5 nm||135±5|
|K3||HRP: 5.5 nm||56±7|
Data taken at 25°C in 50 mm sodium phosphate buffer (pH 7.4), [AR]t=0=280±6 μm, [HRP]=5.5±0.1 nm, [GOx]=82.5±0.7 nm, and [G]t=0=0–200 μm.
To confirm that resorufin is the only optically active molecule in the buffer solution, the spectral absorption (defined as A=1-T in the case of negligible reflection) was measured at different times and normalized to the absorption curve after the system reached saturation. The absorption spectrum does not change over the course of the experiment, thus confirming that other optically active products were not generated. This observation gives further credence to setting K3=0 when solving for K1. If K3 was not zero, then resorufin would undergo further oxidation to OIP, which would change the absorption profile of the reaction over time, which is not observed (see Appendix H, Figure 11).
Shown in Figure 3A are transmittance curves plotted as a function of time for several different initial concentrations of glucose. Figure 3B plots the converted transmittance data to resorufin concentration over time using Eq. (5). Although [G]t=0 was known a priori, this quantity was used as a variable parameter in the fitting code to test the ability of our plasmonic cuvette as an effective sensing scheme to retrieve the glucose concentration initially present in solution; the numbers reported in Figure 3B refer to the fits (reported as solid lines). In all cases, the fitted value matches the concentration of resorufin at saturation, confirming the accuracy of the model. For each value of [G]t=0, a fitted value for K1 and [G]t=0 was found. The weighted average of all K1 data yields K1=7.3±0.3 m-1 s-1 (see Appendix F, Figure 9).
The fitted curves show an excellent correlation to the experimental data, with a small discrepancy occurring when [G]t=0>114 μm at long reaction times. There are two possible explanations for this slight discrepancy between the model and experiment. First, there may be local regions of depleted oxygen, which can be prevented by bubbling oxygen into the buffer solution prior to experimentation. Second, Amplex red may be photo-oxidized to resorufin [40–42]. In order to test this effect, three neutral density filters with different attenuation factors were used to reduce incident intensity. However, no significant change was observed in the measured concentration of resorufin (see Appendix I, Figure 12), thus eliminating this possibility.
2.5 Tuning the reaction time of the AR/GOx/Glucose assay
When applying the plasmonic cuvette to real-time glucose sensing, there are three time-dependent steps that must be considered: (i) the reaction time of the AR/GOx/Glucose assay (typically between 30 and 50 min for standard assays); (ii) the flow time required to deliver the reacted solution to the surface of the sensor (∼10 s); and (iii) the data acquisition time (∼60 s). Since the reaction time of the assay is the longest, this property was investigated further and reduced by increasing the GOx concentration of Reaction 1; this reaction was found to have the smallest rate constant and therefore, was recognized as the rate-limiting step in the assay (Table 1).
Figure 4A reports the concentration of resorufin measured as a function of time for various concentrations of GOx; each curve is normalized by its own saturation value. As expected, increasing [GOx] makes the reaction run faster, but the effect saturates as [GOx] increases. To quantify this effect, Figure 4B plots the reaction time (τ, defined as the time it takes for [R]t to reach 63.2% of its saturation value) as a function of [GOx]. For example, increasing [GOx] by a factor of 10 (from 82.5±0.7 nm to 825±7 nm) shortens the reaction time from 583.5 s to 136.5 s, i.e., a factor of 4.3×. A non-linear dependence between K1 and [GOx] is observed when the data from each concentration of GOx was fitted for a new K1 value, with K2 and K3 unchanged. The non-linear behavior may be caused by aggregation of enzyme at higher concentrations, thus reducing the amount of active sites available to react with glucose (see Appendix J, Figure 13A). At the other extreme, when Reaction 1 occurs very quickly, Reaction 2 takes over as the rate-limiting step in the overall array. A similar experiment was performed for Reaction 3 only. The value for K3 was found as a function of [HRP] by fitting the experimental data to Eqs. (4E) and (4F) with [AR]=0 at all times. The relationship between K3 and the concentration of HRP was found to be linear, as expected from Eq. (4E), with a slope of
2.6 Selectivity test of glucose in artificial saliva
Although the plasmonic cuvette demonstrated high selectivity for glucose in a 50 mm sodium phosphate buffer solution, clinical samples of saliva are more complex and include a mixture of glycoproteins, enzymes, urea and electrolytes . The following experiments were performed in modified Fusayama artificial saliva (MFAS), which differs from actual saliva by the absence of large (>12 kDa) proteins and enzymes (e.g., immunoglobin and α-Amylase, respectively). MFAS contains 16.7 mm of urea; in contrast, the concentration of urea in actual saliva is typically 2–3 mm. Figure 5 shows the production of resorufin over time using the AR/GOx/Glucose assay in MFAS to demonstrate the selectivity of glucose in a complex mixture and in 50 mm sodium phosphate buffer solution as a control. The black squares correspond to the experiment in MFAS (diluted 7× to increase the solubility of the urea and salts in the mixture) without glucose; as expected, in the absence of glucose molecules, the reaction to produce resorufin does not occur. When glucose is included in the reaction mixture, the assay proceeds normally (red circles). Since the composition of MFAS is different than that of 50 mm sodium phosphate buffer solution, the value of ε571nm for MFAS was found to be (4.6±0.2)×104m-1 cm-1. To quantify the difference between MFAS and the buffer solution, the experiment was repeated with the same [G]t=0 value in buffer solution (blue triangles). The two curves are quite similar: at the saturation time (t=1855 s) the difference between the two curves is ∼7.4%. Nothing in the MFAS solution should react with resorufin; thus, the discrepancy between the two curves can be attributed to a change in the overall kinetics of the assay, causing resorufin to be generated at a slightly faster rate in MFAS. A kinetic analysis of the AR/GOx/Glucose assay in complex mixtures such as human saliva at different pH and temperatures will be reported elsewhere.
In conclusion, this work introduces the concept of a plasmonic cuvette: a device that couples an enzyme-driven dye assay, specifically tailored for glucose sensing, to plasmonic interferometry to produce a highly sensitive, selective, real-time, compact optical glucose meter. The sensitivity was found to be 1.7×105%/m within the physiological range of glucose in saliva (20–240 μm). This high sensitivity can be attributed to the strongly absorbing properties of resorufin, which is characterized by unique spectral fingerprints. Our plasmonic cuvette is able to detect glucose concentrations with a resolution (limit of detection, LOD) of ∼0.1 μm. By varying the arm length of the plasmonic interferometer as well as the enzyme concentrations in the dye assay, the plasmonic cuvette studied in this report has been specifically optimized for detection of glucose concentrations in the 20–240 μm range, typically found in human saliva. Our method enables detection limits that are four orders of magnitude smaller than what is currently achievable using conventional electrochemical blood glucose meters (LOD=2.8 mm), employing ∼100 times lower sample volumes. The assay time (∼2 min) was reduced by altering the enzyme concentration of the rate-limiting reaction (i.e., Reaction 1), after determining specific reaction rate constants that underlie the assay and fitting a first-order non-linear kinetic model to the experimental data. Finally, the plasmonic cuvette was tested with artificial saliva to demonstrate the selective detection of glucose in a complex mixture of urea and salts typically found in human saliva. This report demonstrates the first step towards a sensitive and selective sensor for glucose in saliva. However, the methodology reported here is just a proof of concept of a more general approach involving plasmonic interferometry coupled to other dye assays and sensing schemes to screen for various analytes in bodily fluids such as serum, saliva and tears.
The authors gratefully acknowledge funding from the National Science Foundation (NSF Grant No. CBET-1159255, DMR-1203186 and HRD-0548311) and the Juvenile Diabetes Research Foundation (JDRF Grant No. 17-2013-483). This work was performed at the Brown Microelectronics Facility, a member of the Materials Research Facilities Network, which is supported by the National Science Foundation (NSF Grant No. DMR-0520651).
10-Acetyl-3,7-dihydroxyphenoxazine (Amplex red, or AR) and 7-Hydroxy-3H-phenoxazin-3-one sodium salt (resorufin sodium salt) were purchased from Invitrogen (Carlsbad, CA). A stock solution of 19.5 mm AR and 4.5 mm resorufin sodium salt was prepared in anhydrous dimethylsulfoxide (DMSO), aliquoted and stored at -20°C. Horseradish peroxidase (HRP) from Armoracia rusticana (253 purpurogallin U/mg, P-8375, EC184.108.40.206), glucose oxidase (GOx) from Aspergillus niger (100 kU/mg, G7141, EC220.127.116.11), D-(+)-Glucose and urea were purchased from Sigma (St. Louis, MO). Stock concentrations of HRP (10 U/mL) and GOx (1 U/mL) were determined spectrophotometrically using molar extinction coefficients ε403nm=102,000 m-1 cm-1 , and ε280nm=267,200 m-1 cm-1 , respectively. A stock solution of 0.5 mm D-(+)-Glucose was prepared and subsequently allowed to stand for 1 h to complete the mutarotation reaction. The concentration of hydrogen peroxide (35% w/w H2 O2 from Acros Organics) was determined spectrophotometrically at λ=240 nm using a molar extinction coefficient of ε240nm=43.6 m-1 cm-1 . Initial reaction mixtures of the AR/GOx/Glucose assay were prepared at ambient temperature and consisted of 280±6 μm AR, 5.5±0.1 nm HRP and 82.5±0.7 nm GOx dissolved in a 50 mm sodium phosphate buffer solution (pH 7.4) unless otherwise stated. The enzyme stock solutions were thawed and diluted immediately prior to measurement.
Plasmonic interferometers: Fabrication and optical characterization
Arrays of groove-slit-groove (GSG) plasmonic interferometers with arm lengths p1=7.00 μm and p2=7–10 μm in steps of 25 nm (as well as single slit devices for normalization purposes) were fabricated on a 300-nm-thick Ag film deposited on top of a 1 inch by 1 inch quartz slide using focused ion beam (FIB) milling. Each slit (groove) was 100 nm (200 nm) wide and 10 μm long. A 5-nm-thick, optically transparent layer of Al2 O3 was deposited on top of the milled films using atomic layer deposition (ALD) to prevent the silver from undergoing chemical modification in the presence of the buffer solution. A collimated beam of white light (λ=400–800 nm) illuminates the plasmonic interferometer at normal incidence. A 40× objective lens (numerical aperture NA=0.6) was used to collect the transmitted light, which is then dispersed by a single-grating spectrograph and detected by a charge-coupled device (CCD) camera. The spectral resolution of the optical setup was ∼0.4 nm; the number of counts at each wavelength and number of acquired spectra per experiment were adjusted to ensure a statistical error <0.1% in the measured transmitted intensity. An initial reaction mixture of the AR/GOx/Glucose assay was reacted with various concentrations of glucose (0–250 μm) in the dark for 50 min to allow the reaction to go to full completion. The reacted solution was then delivered to the plasmonic interferometer via a polydimethylsiloxane (PDMS) microfluidic channel (2-cm long×1-cm wide×70-μm high) using a syringe pump at a flow rate of 150 μL/min.
UV-Visible absorption spectroscopy
Spectroscopic measurements were performed with a dual-beam UV-Visible Cary 500 Spectrophotometer (Agilent Technologies, Chesterfield, MO) using a quartz cuvette with a path length of L=0.2 cm. The same initial reaction mixture that was used with the plasmonic interferometers was used for both the optical characterization of resorufin and kinetic experiments. For the optical characterization of resorufin, the reaction mixture was initiated with 100 μL of glucose at a concentration of 10.0±0.4 μm, 25±1 μm, 50±2 μm or 100±4 μm, reacted for 50 min in the dark; the transmittance was measured for incident wavelengths between 450 and 700 nm. For the kinetic experiments, the reaction was initiated with the addition of various concentrations of H2 O2 or glucose (both in 100 μL solutions). The reaction was monitored at the wavelength where resorufin exhibits maximum absorption (λmax=571 nm) in 0.1- or 1-s intervals. All experiments were performed in triplicate at standard temperature (25°C) and pressure (760 mm Hg); the concentration of dissolved oxygen ([O2]) was assumed to be constant and equal to 261 μm (calculated using Henry’s law) . All samples were illuminated with a 50 W tungsten light source; the total intensity at the sample surface was measured with a calibrated silicon photodiode and found to be 75 μW cm-2.
Selectivity test for glucose in artificial saliva
Modified Fusayama artificial saliva (MFAS) was prepared to determine the selectivity for glucose in the AR/GOx/Glucose assay. This mixture contains NaCl (6.8 mm), KCl (5.4 mm), CaCl2.2H2 O (5.4 mm), NaH2 PO4.H2 O (5.0 mm), KSCN (3.1 mm), Na2 S.9H2 O (6.9 μm) and urea (16.7 mm) at a pH of 7.4 . A 600-μL solution for measuring glucose concentration was prepared by mixing an initial reaction mixture of the AR/GOx/Glucose assay with 100 μL of MFAS containing 100 μm glucose. Diluting the MFAS solution by 7× increases the solubility of urea and the salts found in the mixture, while still allowing the assay reactions to occur. Control experiments used the same reaction mixture initiated with either 100 μL of the modified Fusayama artificial saliva without glucose or 100 μL of 50 mm sodium phosphate buffer with 100 μm glucose. The transmitted light intensity was measured every 1 s for 2000 s for each reaction and each condition was studied in triplicate.
The plasmonic cuvette merges two disparate fields into one device: dye chemistry and plasmonic interferometry. Because of the various conventions and standards that exist in the realms of chemistry and physics, it is important to be clear with the notation in this manuscript.
There are two equivalent conventions for discussing the amount of molecules in a solution. Concentration c, in units of moles/volume (called a molar, m), is typically used by chemists. Density ρ is typically used by physicists, and is expressed in (length)-3. The conversion factor from one to the other is:
where NA is Avogadro’s constant (NA=6.022×1023 mol-1).
For both cuvettes mentioned in this report (plasmonic and UV-Vis), the quantity that the system measures directly is the spectral transmittance, defined as:
where IT is intensity transmitted through the device and Iref is the reference intensity (i.e., the intensity of the reference beam for the standard UV-Vis cuvette, and the light intensity transmitted through an isolated slit for a plasmonic interferometer). IT and Iref refer to the background-subtracted signals. For a UV-Vis cuvette, under the assumption that no significant reflection occurs, the law of conservation of energy states that:
where A is the absorbance. The Beer-Lambert law defines a closely related quantity, the chemical absorbance (sometimes called absorptance) Ab:
A is typically used by physicists, whereas Ab is typically used by chemists. Although they both share the property of 0<A, Ab<1, they are mathematically distinct.
Next, the absorption coefficient α is defined. This quantity, in units of (length)-1, depends on the intrinsic properties of the absorbing molecule and its concentration in the solution. The first principles definition is:
where L is the path length of the cuvette (in cm). This expression appears in both physics and chemistry. However, it is often useful to define a quantity that can help quantify absorption and absorbance, but does not depend on the concentration of a species with a solution (that is, it is purely an intrinsic property of the absorbing molecule). Once again, physics and chemistry diverge on this point. Chemists prefer to use the extinction coefficient ε (also known as the molar absorptivity), which is defined as:
and is expressed in (molar)-1 (length)-1. Physicists prefer to use the absorption cross-section σ (in units of (length)2), which is defined as:
There is a linear relationship between ε and σ; either of these can be used to relate the time- and wavelength transmission to the measured concentration or density of the optically-active resorufin species. This relationship can be obtained by plugging Eq. (11) into Eq. (9):
This is Eq. (5) in the main text; the notation t refers to the time-dependent concentration of the species within the bracket. Next, plug in the definition of T from Eq. (10):
Using the definition of α from Eq. (12):
This is Eq. (3) in the main text. Finally, we plug in Eq. (6) to obtain the relationship between ε (extinction coefficient) and σ (absorption cross section):
where C1 and C2 are fitting parameters. It was found that C1=48.5±4.2 m-1 s-1 and C2=(4.5±0.9)×10-7m. Although this curve resembles the functional form of the Michaelis-Menten curve, it needs to be emphasized that Michaelis-Menten model is used to find the relationship between the rates of product generation (i.e., velocity v) as a function of substrate concentration ([S]). Instead, Eq. (18) is being used to determine the reaction rate (K1) as a function of enzyme concentration ([GOx]) . The Michaelis-Menten model is not directly applicable to this data, because one of the foundational assumptions used to simplify the model – that the initial substrate concentration is much greater than the concentration of enzyme – is not applicable here. In this system, the initial concentration of substrate (i.e., glucose) is at most only two orders of magnitude higher than the enzyme concentration [GOx] with [G] rapidly declining as it is converted into product (B). Extracted K3 values as a function of HRP concentration indicates a linear dependence with a slope of
World Health Organization. “Diabetes,” Fact Sheet 312 (2013). Retrieved from http://www.who.int/mediacentre/factsheets/fs312/en/.
Wang J. Electrochemical glucose biosensors. Chem Rev 2008;108:814–25.
Wang J. Electrochemical glucose biosensors. Chem Rev 2008;108:814–25.)| false 10.1021/cr068123a
Cornelius T, Pohlmeier H, Behnke T, Schmid V, Grenningloh M, Forst T, Pfützner A. Accuracy evaluation of five blood glucose monitoring systems obtained from the pharmacy: a European multicenter study with 453 subjects. Diabetes Technol Therap 2012;14:330–7.
Cornelius T, Pohlmeier H, Behnke T, Schmid V, Grenningloh M, Forst T, Pfützner A. Accuracy evaluation of five blood glucose monitoring systems obtained from the pharmacy: a European multicenter study with 453 subjects. Diabetes Technol Therap 2012;14:330–7.)| false 10.1089/dia.2011.0170
Lipson H, Bernhardt J, Block U, Freeman WR, Hofmeister R, Hristakeva M, Lenosky T, McNamara R, Petrasek D, Veltkamp D, Waydo S. Requirements for calibration in noninvasive glucose monitoring by Raman spectroscopy. J Diabetes Sci Technol 2009;3:233–41.
Lipson H, Bernhardt J, Block U, Freeman WR, Hofmeister R, Hristakeva M, Lenosky T, McNamara R, Petrasek D, Veltkamp D, Waydo S. Requirements for calibration in noninvasive glucose monitoring by Raman spectroscopy. J Diabetes Sci Technol 2009;3:233–41.)| false 10.1177/193229680900300203
Lee S, Lin R, Chen H, Liu K. US Patent No. 6043492, 2000.
Amir O, Weinstein D, Zilberman S, Less M, Perl-Treves D, Harel P, Weinstein A, Gabis E, Fikhte B, Karasik A. Continuous noninvasive glucose monitoring technology based on ‘occlusion spectroscopy’. J Diabetes Sci Technol 2007;1:463–9.
Amir O, Weinstein D, Zilberman S, Less M, Perl-Treves D, Harel P, Weinstein A, Gabis E, Fikhte B, Karasik A. Continuous noninvasive glucose monitoring technology based on ‘occlusion spectroscopy’. J Diabetes Sci Technol 2007;1:463–9.)| false 10.1177/193229680700100403
Harman-Boehm I, Gal A, Raykhman AM, Zahn JD, Naidis E, Mayzel Y. Noninvasive glucose monitoring: a novel approach. J Diabetes Sci Technol 2009;3:253–60.
Harman-Boehm I, Gal A, Raykhman AM, Zahn JD, Naidis E, Mayzel Y. Noninvasive glucose monitoring: a novel approach. J Diabetes Sci Technol 2009;3:253–60.)| false 10.1177/193229680900300205
Liu R, Chen W, Gu X, Wang RK, Xu K. Chance correlation in non-invasive glucose measurement using near-infrared spectroscopy. J Phys D Appl Phys 2005;38:2675–81.
Liu R, Chen W, Gu X, Wang RK, Xu K. Chance correlation in non-invasive glucose measurement using near-infrared spectroscopy. J Phys D Appl Phys 2005;38:2675–81.)| false 10.1088/0022-3727/38/15/020
Weiss R, Yegorchikov Y, Shusterman A, Raz I. Noninvasive continuous glucose monitoring using photoacoustic technology – results from the first 62 subjects. Diabetes Technol Therap 2007;9:68–74.
Weiss R, Yegorchikov Y, Shusterman A, Raz I. Noninvasive continuous glucose monitoring using photoacoustic technology – results from the first 62 subjects. Diabetes Technol Therap 2007;9:68–74.)| false 10.1089/dia.2006.0059
Lyandres O, Yuen JM, Shah NC, Van Duyne RP, Walsh JT, Glucksberg MR. Progress toward an in vivo surface-enhanced Raman spectroscopy glucose sensor. Diabetes Technol Therap 2008;10:257–65.
Purvinis G, Cameron BD, Altrogge DM. Noninvasive polarimetric-based glucose monitoring: an in vivo study. J Diabetes Sci Technol 2011;5:380–7.
Purvinis G, Cameron BD, Altrogge DM. Noninvasive polarimetric-based glucose monitoring: an)| false in vivostudy. J Diabetes Sci Technol 2011;5:380–7. 10.1177/193229681100500227
Gabbay RA, Sivarajah S. Optical coherence tomography-based continuous noninvasive glucose monitoring in patients with diabetes. Diabetes Technol Therap 2008;10:188–93.
Thennadil SN, Rennert JL, Wenzel BJ, Hazen KH, Ruchti TL, Block MB. Comparison of glucose concentration in interstitial fluid, and capillary and venous blood during rapid changes in blood glucose levels. Diabetes Technol Therap 2001;3: 357–65.
Thennadil SN, Rennert JL, Wenzel BJ, Hazen KH, Ruchti TL, Block MB. Comparison of glucose concentration in interstitial fluid, and capillary and venous blood during rapid changes in blood glucose levels. Diabetes Technol Therap 2001;3: 357–65.)| false 10.1089/15209150152607132
Yamaguchi M, Mitsumori M, Kano Y. Noninvasively measuring blood glucose using saliva. IEEE Eng Med Biol 1998;17:59–63.
Yamaguchi M, Mitsumori M, Kano Y. Noninvasively measuring blood glucose using saliva. IEEE Eng Med Biol 1998;17:59–63.)| false 10.1109/51.677170
Panchbhai AS. Correlation of salivary glucose level with blood glucose level in diabetes mellitus. J Oral Maxil Res 2012;3:1–7.
Panchbhai AS. Correlation of salivary glucose level with blood glucose level in diabetes mellitus. J Oral Maxil Res 2012;3:1–7.)| false 10.5037/jomr.2012.3303
Agrawal A, Susut C, Stafford G, Bertocci U, McMorran B, Lezec HJ, Talin AA. An integrated electrochromic nanoplasmonic optical switch. Nano Lett 2011;11:2774–8.
Agrawal A, Susut C, Stafford G, Bertocci U, McMorran B, Lezec HJ, Talin AA. An integrated electrochromic nanoplasmonic optical switch. Nano Lett 2011;11:2774–8.)| false 21612218
Homola J. Surface plasmon resonance sensors for detection of chemical and biological species. Chem Rev 2008;108:462–93.
Homola J. Surface plasmon resonance sensors for detection of chemical and biological species. Chem Rev 2008;108:462–93.)| false 10.1021/cr068107d
Haes AJ, Van Duyne RP. A nanoscale optical biosensor: sensitivity and selectivity of an approach based on localized surface plasmon resonance spectroscopy of triangular silver nanoparticles. J Am Chem Soc 2002;124:10596–604.
Haes AJ, Van Duyne RP. A nanoscale optical biosensor: sensitivity and selectivity of an approach based on localized surface plasmon resonance spectroscopy of triangular silver nanoparticles. J Am Chem Soc 2002;124:10596–604.)| false 10.1021/ja020393x
Lal S, Link S, Halas NJ. Nano-optics from sensing to waveguiding. Nat Photon 2007;1:641–8.
Yonzon CR, Jeung E, Zou S, Schatz GC, Mrksich M, Van Duyne RP. A comparative analysis of localized and propagating surface plasmon resonance sensors: the binding of concanavalin A to a monosaccharide functionalized self-assembled monolayer. J Am Chem Soc 2004;126:12669–76.
Yonzon CR, Jeung E, Zou S, Schatz GC, Mrksich M, Van Duyne RP. A comparative analysis of localized and propagating surface plasmon resonance sensors: the binding of concanavalin A to a monosaccharide functionalized self-assembled monolayer. J Am Chem Soc 2004;126:12669–76.)| false 10.1021/ja047118q
Brolo AG. Plasmonics for future biosensors. Nat Photon 2012;6:709–13.
Brolo AG. Plasmonics for future biosensors. Nat Photon 2012;6:709–13.)| false 10.1038/nphoton.2012.266
Hsieh HV, Pfeiffer ZA, Amiss TJ, Sherman DB, Pitner JB. Direct detection of glucose by surface plasmon resonance with bacterial glucose/galactose-binding protein. Biosens Bioelectron 2004;19:653–60.
Lam WW, Chu LH, Wong CL, Zhang YT. A surface plasmon resonance system for the measurement of glucose in aqueous solution. Sensor Actuat B-Chem 2005;105:138–43.
Lam WW, Chu LH, Wong CL, Zhang YT. A surface plasmon resonance system for the measurement of glucose in aqueous solution. Sensor Actuat B-Chem 2005;105:138–43.)| false 10.1016/j.snb.2004.04.088
Ozbay E. Plasmonics: merging photonics and electronics at nanoscale dimensions. Science 2006;311:189–93.
Ozbay E. Plasmonics: merging photonics and electronics at nanoscale dimensions. Science 2006;311:189–93.)| false 10.1126/science.1114849
Pacifici D, Lezec HJ, Atwater HA. All-optical modulation by plasmonic excitation of CdSe quantum dots. Nat Photon 2007;1:402–6.
Pacifici D, Lezec HJ, Atwater HA. All-optical modulation by plasmonic excitation of CdSe quantum dots. Nat Photon 2007;1:402–6.)| false 10.1038/nphoton.2007.95
Feng J, Siu VS, Roelke A, Mehta V, Rhieu SY, Palmore GTR, Pacifici D. Nanoscale plasmonic interferometers for multispectral, high-throughput biochemical sensing. Nano Lett 2012;12:602–9.
Raether HR. Surface plasmons on smooth and rough surfaces and on gratings. New York: Springer, 1988.
Raether HR. Surface plasmons on smooth and rough surfaces and on gratings. New York: Springer, 1988.)| false 10.1007/BFb0048317
Jurysta C, Bulur N, Oguzhan B, Satman I, Yilmaz TM, Malaisse WJ, Sener A. Salivary glucose concentrations and excretion in normal and diabetic subjects. J Biomed Biotechnol 2009;2009:430426.
Jurysta C, Bulur N, Oguzhan B, Satman I, Yilmaz TM, Malaisse WJ, Sener A. Salivary glucose concentrations and excretion in normal and diabetic subjects. J Biomed Biotechnol 2009;2009:430426.)| false 10.1155/2009/430426
Edgar WM. Saliva: its secretion, composition and functions. Brit Dent J 1992;172:305–12.
Edgar WM. Saliva: its secretion, composition and functions. Brit Dent J 1992;172:305–12.)| false 10.1038/sj.bdj.4807861
Golub E, Freeman R, Niazov A, Willner I. Hemin/G-quadruplexes as DNAzymes for the fluorescent detection of DNA, aptamer-thrombin complexes, and probing the activity of glucose oxidase. Analyst 2011;136:4397–401.
Golub E, Freeman R, Niazov A, Willner I. Hemin/G-quadruplexes as DNAzymes for the fluorescent detection of DNA, aptamer-thrombin complexes, and probing the activity of glucose oxidase. Analyst 2011;136:4397–401.)| false 10.1039/c1an15596b
Zhou M, Diwu Z, Panchuk-Voloshina N, Haugland RP. A stable non-fluorescent derivative of resorufin for the fluorometric determination of trace hydrogen peroxide: applications in detecting the activity of phagocyte NADPH oxidase and other oxidases. Anal Biochem 1997;253:162–8.
Gorris HH, Walt DR. Mechanistic aspects of horseradish peroxidase elucidated through single-molecule studies. J Am Chem Soc 2009;131:6277–82.
Gorris HH, Walt DR. Mechanistic aspects of horseradish peroxidase elucidated through single-molecule studies. J Am Chem Soc 2009;131:6277–82.)| false 10.1021/ja9008858
Towne V, Will M, Oswald B, Zha Q. Complexities in horseradish peroxidase-catalyzed oxidation of dihydroxyphenoxazine derivatives: appropriate ranges for pH values and hydrogen peroxide concentrations in quantitative analysis. Anal Biochem 2004;334:290–6.
Towne V, Will M, Oswald B, Zha Q. Complexities in horseradish peroxidase-catalyzed oxidation of dihydroxyphenoxazine derivatives: appropriate ranges for pH values and hydrogen peroxide concentrations in quantitative analysis. Anal Biochem 2004;334:290–6.)| false 10.1016/j.ab.2004.07.037
Brotea GP, Thibert RJ. Fluorometric determination of hydrogen peroxide using resorufin and peroxidase. Microchem J 1988;37:368–76.
Brotea GP, Thibert RJ. Fluorometric determination of hydrogen peroxide using resorufin and peroxidase. Microchem J 1988;37:368–76.)| false 10.1016/0026-265X(88)90150-6
Ramanathan A, Wang C, Schreiber SL. Perturbational profiling of a cell-line model of tumorigenesis by using metabolic measurements. Proc Natl Acad Sci USA 2005;102:5992–7.
Ramanathan A, Wang C, Schreiber SL. Perturbational profiling of a cell-line model of tumorigenesis by using metabolic measurements. Proc Natl Acad Sci USA 2005;102:5992–7.)| false 10.1073/pnas.0502267102
Mazura P, Fohlerová R, Brzobohatý B, Kiran NS, Janda L. A new, sensitive method for enzyme kinetic studies of scarce glucosides. J Biochem Biophys Meth 2006;68:55–63.
Mazura P, Fohlerová R, Brzobohatý B, Kiran NS, Janda L. A new, sensitive method for enzyme kinetic studies of scarce glucosides. J Biochem Biophys Meth 2006;68:55–63.)| false 10.1016/j.jbbm.2006.03.018
Molecular probes, Inc. Amplex® red glucose/glucose oxidase assay kit instructions. Eugene, OR: Molecular probes, Inc., 2006.
Hall MB, Keuler NS. Factors affecting accuracy and time requirements of a glucose oxidase-peroxidase assay for determination of glucose. J AOAC Int 2009;92:50–60.
Nakajima H, Ishino S, Masuda H, Nakagama T, Shimosaka T, Uchiyama K. Photochemical immobilization of protein on the inner wall of a microchannel and its application in a glucose sensor. Anal Chim Acta 2006;562:103–9.
Nakajima H, Ishino S, Masuda H, Nakagama T, Shimosaka T, Uchiyama K. Photochemical immobilization of protein on the inner wall of a microchannel and its application in a glucose sensor. Anal Chim Acta 2006;562:103–9.)| false 10.1016/j.aca.2005.12.051
Bueno C, Villegas ML, Bertolotti SG, Previtali CM, Neumann MG, Encinas MV. The excited-state interaction of resazurin and resorufin with amines in aqueous solutions. Photophysics and photochemical reaction. Photochem Photobiol 2002;76:385–90.
Bueno C, Villegas ML, Bertolotti SG, Previtali CM, Neumann MG, Encinas MV. The excited-state interaction of resazurin and resorufin with amines in aqueous solutions. Photophysics and photochemical reaction. Photochem Photobiol 2002;76:385–90.)| false 10.1562/0031-8655(2002)0760385TESIOR2.0.CO2
Zhao B, Summers FA, Mason RP. Photooxidation of Amplex red to resorufin: implications of exposing the Amplex red assay to light. Free Radical Bio Med 2012;53:1080–7.
Zhao B, Summers FA, Mason RP. Photooxidation of Amplex red to resorufin: implications of exposing the Amplex red assay to light. Free Radical Bio Med 2012;53:1080–7.)| false 10.1016/j.freeradbiomed.2012.06.034
Wang Z, King TL, Branagan SP, Bohn PW. Enzymatic activity of surface-immobilized horseradish peroxidase confined to micrometer- to nanometer-scale structures in nanocapillary array membranes. Analyst 2009;134:851–9.
Wang Z, King TL, Branagan SP, Bohn PW. Enzymatic activity of surface-immobilized horseradish peroxidase confined to micrometer- to nanometer-scale structures in nanocapillary array membranes. Analyst 2009;134:851–9.)| false 10.1039/b815590a
George P. Chemical nature of the secondary hydrogen peroxide compound formed by cytochrome-C peroxidase and horseradish peroxidase. Nature 1952;169:612–3.
George P. Chemical nature of the secondary hydrogen peroxide compound formed by cytochrome-C peroxidase and horseradish peroxidase. Nature 1952;169:612–3.)| false 10.1038/169612a0
Fasman GD. CRC handbook of biochemistry and molecular biology. Cleveland, OH: CRC Press, 1990.
Noble RW, Gibson QH. The reaction of ferrous horseradish peroxidase with hydrogen peroxide. J Biol Chem 1970;245:2409–13.
Dean JA. Lange’s handbook of chemistry. New York: McGraw-Hill, 1990.
Huang H-H, Chiu Y-H, Lee T-H, Wu S-C, Yang H-W, Su K-H, Hsu C-C. Ion release from NiTi orthodontic wires in artificial saliva with various acidities. Biomaterials 2003;24:3585–92.
Lehninger AL, Nelson DL, Cox MM. Lehninger principles of biochemistry. 4th ed. New York: W.H. Freeman, 2005.