Electronic Olfactory Sensor Based on A. mellifera Odorant-Binding Protein 14 on a Reduced Graphene Oxide Field-Effect Transistor.

An olfactory biosensor based on a reduced graphene oxide (rGO) field-effect transistor (FET), functionalized by the odorant-binding protein 14 (OBP14) from the honey bee (Apis mellifera) has been designed for the in situ and real-time monitoring of a broad spectrum of odorants in aqueous solutions known to be attractants for bees. The electrical measurements of the binding of all tested odorants are shown to follow the Langmuir model for ligand-receptor interactions. The results demonstrate that OBP14 is able to bind odorants even after immobilization on rGO and can discriminate between ligands binding within a range of dissociation constants from K(d)=4 μM to K(d)=3.3 mM. The strongest ligands, such as homovanillic acid, eugenol, and methyl vanillate all contain a hydroxy group which is apparently important for the strong interaction with the protein.

Odorant‐binding proteins (OBPs) are small acidic proteins (ca. 13–16 kDa) highly concentrated in the lymph of the chemosensillae of insects or in the nasal mucus of vertebrates.1 They act as carriers for volatile organic compounds, VOCs, (air‐borne odorants) shuttling them from the air–water interface to the membrane‐integral odorant receptor. Although their full function has not been completely clarified yet, OBPs certainly play a major role in detecting and recognizing olfactory stimuli.2

A large number of OBPs has been expressed in bacterial systems and their ligand‐binding properties have been investigated in solution by a fluorescent ligand displacement assay.3 Generally, dissociation constants, K d, are in the upper nanomolar or lower micromolar range for strong odorants.1c,d An important characteristic of OBPs for technical applications is their stability to extreme temperatures, solvents, and proteolysis, making such proteins ideal elements for biosensors to be used in medical application, for example, in breath analysis for cancer diagnostics, for food quality control, for crop‐disease detection, or in general environmental monitoring.4a–f

Various biosensor devices mimicking the olfactory system (artificial noses) have been developed but only few studies have used OBPs as functional elements for the design of a “bio‐electronic nose”.5 Herein we present the fabrication and functional characterization of a label‐free biosensor for odorant detection based on a reduced graphene oxide field‐effect transistor functionalized with the odorant‐binding protein 14 (OBP14) from the honey bee Apis mellifera.

Reduced graphene oxide field‐effect transistor (rGO‐FET) devices were fabricated according to established methods, schematically given in Figure 1 A.6 A scanning electron microscopic image of rGO flakes assembled onto the gate substrate before any further surface functionalization is shown in Figure 1 B. IR spectra of the linker monolayer, 1‐pyrenebutanoic acid succinimidyl ester (PBSE), that is typically used for protein immobilization on graphene substrates,7 with partial covalent immobilization of OBP taken at different times during the assembly from solution to the gate are given in Figure 1 C, while Figure 1 D summarizes the quantitative analysis of the functionalization process by monitoring the time‐dependent increase of the amide I and II bands of the protein and the corresponding decrease of the band at 1738 cm−1, assigned to the cleavage of the active ester group of the linker molecules. The entire fabrication process of the olfactory biosensor device, including the successful reduction of GO to rGO, the linker binding, and more details of the protein attachment are described in the Supporting Information. The cloning, expression, and purification of OBPs, the preparation of odorant solutions and the performance of the electrical measurements are also described in the Supporting Information.

Figure 1

A) Schematic illustration of the individual fabrication steps of the graphene biosensor device. B) Scanning electron microscopic image of the rGO‐FET before the functionalization with 1‐pyrene‐butanoic acid succinimidyl ester (PBSE) linker. C) Infrared spectra of the PBSE linker attached to gate area of the rGO surface, and OBP14 immobilized for 5 and 35 min, respectively, as indicated (spectra measured in ATR configuration and water corrected). D) Time‐dependent increase of the amide I and II bands, upon binding of OBP14 to the linker molecules at the gate surface, and the corresponding decrease of the band at 1738 cm−1 upon cleavage of the active ester of the linker molecules during the protein immobilization (cf. also Figure S6). The dashed red curves are guides to the eye.

The mode of operation of the device as a field‐effect transistor is demonstrated in Figure S5. Recording the ISD versus VG scans under different bulk solution conditions, in particular, in aqueous solutions with different ligand concentrations results in a slightly modified slope of the cathodic branch. We attribute this change to a slight modification of the dipolar layer upon binding of the ligands to the free binding sites in the OBP which act as receptors. The OBP protein monolayer is immobilized on the graphene gate and ligand binding causes a the slight reorientation of its the alpha‐helical parts. At the selected gate voltage of V G=−0.6 V the concentration dependency of the slopes in the I SD versus V G curves could be used to measure the binding of odorants to OBP14 in real‐time, resulting in the quantitative determination of the kinetic rate constants for the association process, k on, and for the dissociation process, k off, as well as for the affinity constant K A and the dissociation constant K d. As an example, Figure 2 A shows a global analysis measurement, that is, the time dependence of the change of the source‐drain current, ΔI SD, of the FET as a function of time upon binding of methyl vanillate, a strong binder, from solution to the OBP14‐functionalized gate surface at increasing and decreasing odorant bulk concentrations. The error of the kinetic rate constants, k, obtained was estimated to be around 20 % (cf. also Figure S8 in the Supporting Information).

Figure 2

A) Real‐time biosensor measurement of the binding of methyl vanillate to OBP14: the current increases with the bulk concentration of methyl vanillate increasing (from 100 nm to 200 μm) and then saturates. Blue arrows indicate runs with pure buffer, red arrows indicate experiments with methyl vanillate solutions. Red curves are the fitting of the raw data by kinetic simulations of the association and dissociation processes based on the Langmuir model. (For the estimation of the error limits, see the Supporting Information). B) Analysis of the reaction rate constants, k, obtained from the fitted data in (A) as a function of methyl vanillate concentration; different symbols from three different devices; error bars are ±20 %. C) Langmuir adsorption isotherm, obtained for three different samples; the red fit curve gives K d=20 μm (plotted are also error limits for K d of 20 %).

Upon plotting the resulting reaction rate constants k, as they were derived from the fits to the kinetic traces as a function of the bulk ligand concentration, c 0=c Methyl vanillate, (Figure 2 B), gives a straight line, the slope of which, according to k=k on  c 0+k off, with k on being the association and k off the dissociation rate constant, respectively, yields k on=235 m −1 s−1, and the intersection with the y‐axis gives k off=0.01 s−1. According to the Langmuir model the ratio k on/k off gives the affinity constant K A=2.3×104  m −1, which can be compared to the value derived from the equilibrium titration experiment given in Figure 2 C: according to θ=K A  c 0/(1+K A  c 0) the fit to the data (full red curve plus the error limits of ±20 %) gives the affinity constant K A=5×104  m −1 which compares quite well with the value obtained from the kinetic experiments thus confirming the Langmuir model for this binding process.

To exclude false signals arising from, in particular, non‐specific binding, several control experiments were performed. Firstly, a sensor was prepared which was covered only with linker molecules, without the coupling of the odorant‐binding protein. Even high concentrations of the strong binder homovanillic acid (cf. Figure S7A) or other ligands (Figure S7B) did not result in a significant signal. Coating the gate with OBP14 but exposing it to a totally uncorrelated small molecule, biotin, gave no signal (Figure S7C). And finally, immobilizing a protein (OBP9A from the red flour beetle, Tribolium castaneum) that is structurally similar but is not a receptor for these ligand odorants also gave a very weak signal (Figure S7D).

Further evidence for the specificity of the ligand receptor binding originates from a direct comparison of the binding of eugenol and methyl eugenol to the same chip, functionalized with OBP14. Figure 3 A shows a real‐time current trace recorded during the addition of methyl eugenol at 50 μm, then rinsing with pure buffer, and then a 5 μm solution of eugenol being rinsed through the flow cell. Despite the minor chemical variation between the two ligands a large difference in the affinity of their binding reaction to OBP14 confirms the specificity of the sensor. Figure 3 B summarizes the Langmuir isotherms for both ligands, demonstrating the significantly different dissociation constants: K d=40 μm for eugenol and K d=1400 μm for methyl eugenol.

Figure 3

A) Real‐time sensor response to the injection of a 50 μm solution of methyl eugenol, subsequent buffer wash, and the injection of 5 μm eugenol solution. B) Langmuir isotherms of eugenol and methyl eugenol.

Table 1 summarizes the quantitative data measured for a series of ligands and gives the kinetic rate constants, k on and k off, obtained as well as the dissociation constant K d. It could be confirmed also for this label‐free sensing format that the odorant binding affinity to OBP14 gradually decreases from homovanillic acid to citral and to geraniol.8 Molecules which are structurally related to eugenol also showed strong affinities to OBP14.

Table 1

Dissociation constants, K d, association rate constants, k on, and dissociation rate constants, k off, for a variety of odorants binding to OBP14 as obtained by the global analysis.

Odorant	Homovanillic acid	Methyl vanillate	Eugenol	Citral	Methyl eugenol	Geraniol
structural formula	WILEY-VCH	WILEY-VCH	WILEY-VCH	WILEY-VCH	WILEY-VCH	WILEY-VCH
K d [×10−6  m]	4	20	40	800	1400	3300
k on [m −1 s−1]	1130	235	170	9	6	3
k off [s−1]	0.008	0.01	0.006	0.003	0.006	0.008

Interestingly, also other odorant‐binding proteins bind their odorants with affinities in the same order of magnitude.9 It is remarkable to note that for all the ligands investigated the dissociation rate constants k off, differ by less than a factor of 2.5. The strongly differing affinity constants, varying by nearly three orders of magnitude, can be almost exclusively attributed to the differences in the association rate constants k on, (cf. Table 1).

The Spinelli group used X‐ray diffraction analysis to determine the structure of the eugenol–OBP14 complex. They found that the hydroxy group of eugenol interacts with the cavity wall of the OBP14 binding pocket by forming two hydrogen bonds.10

It has been speculated that the hydroxy group, together with the substituted aromatic backbone, plays a key role for the strong binding which is in agreement with the high affinities also of homovanillic acid and methyl vanillate (cf. Table 1).10 This hypothesis is further supported by the lower affinity measured for methyl eugenol, in which the hydroxy group is replaced by a methoxy group.

As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.

miscellaneous_information

Click here for additional data file.