Alice R.
Hewson
a,
Henry O.
Lloyd-Laney
a,
Tessa
Keenan
a,
Sarah-Jane
Richards
bc,
Matthew I.
Gibson
bc,
Bruno
Linclau
de,
Nathalie
Signoret
f,
Martin A.
Fascione
*a and
Alison
Parkin
*a
aDepartment of Chemistry, University of York, YO10 5DD, York, UK. E-mail: alison.parkin@york.ac.uk; martin.fascione@york.ac.uk
bDepartment of Chemistry, The University of Manchester, M13 9PL, UK
cManchester Institute of Biotechnology, The University of Manchester, M1 7DN, UK
dDepartment of Organic and Macromolecular Chemistry, Ghent University, Krijgslaan 281-S4, 9000 Gent, Belgium
eSchool of Chemistry, University of Southampton, Highfield, Southampton, SO17 1BJ, UK
fHull York Medical School, University of York, YO10 5DD, York, UK
First published on 13th September 2024
Glycans play a major role in biological cell–cell recognition and signal transduction but have found limited application in biosensors due to glycan/lectin promiscuity; multiple proteins are capable of binding to the same native glycan. Here, site-specific fluorination is used to introduce protein–glycan selectivity, and this is coupled with an electrochemical detection method to generate a novel biosensor platform. 3F-lacto-N-biose glycofluoroform is installed onto polymer tethers, which are subsequently immobilised onto gold screen printed electrodes, providing a non-fouling surface. The impedance biosensing platform is shown to selectively bind cancer-associated galectin-3 compared to control glycans and proteins. To improve the analytical capability, Bayesian statistical analysis was deployed in the equivalent circuit fitting of electrochemical impedance spectroscopy data. It is shown that Markov Chain Monte Carlo (MCMC) analysis is a helpful method for visualising experimental irreproducibility, and we apply this as a quality control step.
Proteins, carbohydrates, nucleic acids and lipids are the four major classes of biological macromolecules. While proteins (often in the form of enzymes and antibodies) and nucleic acids (e.g. aptamers) are commonly utilised as bioreceptors in impedimetric biosensors, and lipid layers immobilised on surfaces can be used to capture receptors for sensing,2–4 carbohydrates are not widely harnessed as bioreceptors.5 This is despite the established major roles that carbohydrate–protein interactions play in a range of biological processes including signalling, cell adhesion, agglutination, protein folding and the immune response.6–8 The carbohydrate binding proteins integral to these interactions are known as lectins.9 These exhibit reversible and non-covalent carbohydrate binding, both in solution and on surfaces, with substrates of varying complexity.9,10 Lectins are sub-classified into large families primarily based on the structure of their carbohydrate recognition domain (CRD), and their carbohydrate binding specificity.11 Although lectins can bind oligosaccharides with high specificity, they often bind simpler mono- and disaccharides less discriminately. This innate non-selective small molecule binding,12 combined with the fact that multiple lectins can also bind the same carbohydrate, is a significant barrier in the development of carbohydrate receptors for biosensor applications, with only a handful of examples of biosensor platforms in the literature.13–22 Such constraints are exemplified by the galectin family – proteins that are prime targets for biosensors because they have been implicated in adhesion and inflammation in several diseases,23 with galectin-3 specifically identified as a biomarker for acute and chronic heart failure, liver fibrosis, renal failure, and cancer of the prostate, bladder, thyroid, lung, breast, colon and rectum.24–27 The construction of a carbohydrate-based biosensor specific for galectin-3 is challenging because all galectin proteins share a conserved binding site within their CRD.23
An emerging solution to the challenge of achieving selective lectin-binding to mono- and disaccharides is the use of glycofluoroforms (GFFs),28 unnatural glycomimetics where fluorine atoms typically replace one or more hydroxyl groups in the carbohydrate scaffold. Although the presence of an electronegative fluorine can modulate hydrogen bonding capability and lipophilicity, the conformational effects are often minimal.29 This attractive combination of electronic modulation and minimal structural perturbation (a ‘molecular editing’ approach) has previously been applied in the construction of targeted probes and tools for structural biology, rational drug design and medicinal therapeutics,30,31 and was recently harnessed in the development of a library of GFFs based on the same carbohydrate scaffold, capable of selective binding to different protein targets.28,32,33 Notably, Richards et al. showed that in screening experiments using 3FGal-β(1,3)-GlcNAc, a lacto-N-biose glycofluoroform, 8-fold selective binding to galectin-3 (KD = 6 nM) over galectin-7 could be achieved.32
Herein, we demonstrate how strong, selective binding between 3F-Gal-β(1,3)-GlcNAc and galectin-3 can be harnessed to develop the first GFF electrochemical biosensor (Fig. 1), in which we utilise a thiol terminated polymer–GFF conjugate immobilised on a gold screen-printed electrode (Au-SPE). We note that Fig. 1 depicts a multimeric structure for galectin-3 and a multivalent interaction with the surface in line with the proposed biological function.35 The successful detection of galectin-3 by the GFF-biosensor is subsequently demonstrated using EIS, which measures the current–time response to a small amplitude sinusoidal voltage–time oscillation of variable frequency.36,37 The data is visualised using both Bode and Nyquist plots to represent changes in impedance, Z, the ratio of voltage to current, as a function of frequency. The impedance has a real, Z′, and imaginary, Z′′, component and the magnitude is denoted |Z|. The EIS response is influenced by changes in resistive forces associated with a binding event on the surface of the working electrode. This can be quantitatively determined by equivalent circuit fitting, and such data analysis is significantly more rapid than solving differential equations, as required for the quantitative analysis of voltammetric data. This combination of high surface sensitivity and ease of quantitative analysis makes EIS an attractive detection method for point-of-care diagnostic devices. Indeed, prior studies by other authors have already established that EIS sensors can be directly applied to biofluids with high recovery rates of analytes within clinically relevant ranges.38–41 Statistical analytic methods, namely Bayesian inference and Monte Carlo Markov chain, enable us to visualise the confidence and relationship between the modelling parameters used, which highlights the potential of the platform for impedimetric biosensing, and enables critical analysis of the data quality.
Fig. 1 An overview of the GFF-biosensor for galectin-3 detection. (a) The chemical structures and schematic representations of the core building blocks of the sensor: glycofluoroform Gal3F-β(1,3)-GlcNAc-N3 and polymeric thiol linker DBCO-(PHEA)25-SH. (b) Schematic diagram of the assembled biosensor and the proposed galectin-3 binding interaction. (c) Binding of an analogous native glycoform, LacNAc, on the surface of galectin-3 (PDB: 1KJL).34 |
The control sugar, GlcNAc-N3, formally 3-azidopropyl 2-acetamido-2-deoxy-β-D-glucopyranoside, was synthesised as described previously.32
The polymer, DBCO-(PHEA)25-SH, was synthesised as described previously.32
Polymer–sugar conjugates were made by reacting the thiol polymer and either the glycofluoroform, Gal3F-β(1,3)-GlcNAc-N3, or the control sugar, GlcNAc-N3, via a strain-promoted azide–alkyne cycloaddition (SPAAC) reaction. Appropriate volumes of polymer (0.32 μmol, 1 equiv.) and sugar (0.64 μmol, 2 equiv.) solutions (both in Milli-Q H2O) were added to a 2 mL Eppendorf tube and made up to 450 μL with Milli-Q H2O and shaken at room temperature in an anaerobic glovebox environment for 18–24 hours.
Unless otherwise stated, following the SPAAC reaction 450 μL aliquots of the reaction mixture in 2 mL Eppendorf tubes were used to modify two Au-SPEs simultaneously by immersing the SPEs to a depth where the working, counter and reference electrodes were covered. The SPE modification was left to proceed at room temperature in an anaerobic glovebox environment for 18–24 hours. Following Au-SPE modification, the SPAAC reaction aliquots are stored at −20 °C and then thawed for subsequent Au-SPE modifications.
To assess the protein binding affinity of a single electrode, if necessary, the electrode was first modified with sugar-polymer. Next, a protein-free initial electrochemical impedance spectroscopy (EIS) measurement was carried out, as described below. The electrode was then rinsed thoroughly to remove all potassium ferricyanide solution before being incubated in the appropriate protein solution, starting at the lowest concentration to be tested within the dataset. The electrode was immersed in approximately 400 μL of protein (sufficient to cover all working electrode surface area) for 30 min (timed using a stopwatch) at room temperature, before being rinsed with Milli-Q H2O and measured in a solution of buffered potassium ferricyanide using EIS (i.e. electrochemical measurements conducted without protein in solution). Following EIS, the electrode was then rinsed thoroughly with Milli-Q H2O before being incubated in the next concentration of protein, followed by another EIS measurement. This two-step, protein incubation followed by EIS measurement, process is repeated for each protein sample to be tested within the dataset, with an incremental increase in concentration each time.
To generate posterior distributions for the parameter values, Markov chain Monte Carlo (MCMC) analysis was utilised, using the adaptive covariance MCMC algorithm as implemented in the PINTS repository,45 using the same log-likelihood function as maximised by CMA-ES, and an uninformative log-prior using the boundaries for CMA-ES. Three chains were run independently for 10000 iterations, initialised from the maximum likelihood point determined by CMA-ES. The pooled samples from these chains are plotted as histograms, with the first 3000 samples discarded as “burn-in”. Convergence information can be found in the ESI.†
To prove the formation of gold-thiol bonds between the electrode and GFF-polymer conjugate, we first used EIS to compare the response of ‘bare’ Au-SPEs (i.e., new, from the packet electrodes) to SPEs that have been modified with the SPAAC polymer–glycofluoroform (P–GFF) reaction mixture. Electrochemical measurements carried out in the presence of solution-phase potassium ferricyanide showed the expected response for a solution EIS experiment in the Nyquist plot (Fig. 2),48 namely a semi-circle with a linear region at approximately 45°. This indicated that analysis using a modified Randles circuit (Fig. 2, inset) was appropriate, with computational analysis generating a good fit to the experimental data. Notably, a distinct difference can be observed between the bare Au-SPE and the P–GFF modified Au-SPE. The large magnitude of the Nyquist plot for the P–GFF modified Au-SPE relative to a bare Au-SPE indicated that incubation of a Au-SPE in SPAAC P–GFF reaction mixture resulted in gold–thiol bond formation, with the presence of the coating on the gold surface hindering electron transfer between solution ferricyanide and the electrode surface. Numerically, this can be quantified as the change in the electron transfer rate, as discussed later.49 Further evidence that incubation of an Au-SPE in P–GFF solution generates a gold–thiol bond is provided by electrochemical ‘stripping’ experiments (Fig. S1†).
Complementary control experiments were also carried out to assess the surface coverage of polymer on the electrode by replacing the glycofluoroform with a ferrocene derivative. A surface-confined ferrocene DCV signal was observed (Fig. S2†) and analysed using a third-order polynomial baseline subtraction to isolate the faradaic current. As detailed in the ESI,† from a simple integration of the total charge passed and assuming a one-electron transfer for each ferrocene molecule, the surface coverage in this experiment was found to be 0.35 pmol which is approximately three-fold lower than the estimated monolayer coverage, 1.01 pmol (Fig. S3†). Since the polymer structure is consistent, and the gold–thiol bond formation chemistry will be unchanged, we also interpret the ferrocene measurements as providing a proxy measurement of the coverage of GFF on the electrode surface. Other control experiments also indicate non-ideal, sub-monolayer formation (vide infra).
EIS experiments comparing bare and P–GFF modified electrodes were repeated extensively (Fig. S4 and S5)† and the Nyquist plots showed substantive variation within the P–GFF dataset that is not present in the bare Au-SPE dataset. This variation could not be reduced by careful control of experimental conditions. Even construction of a specialised tank to enable the simultaneous immersion of eight SPEs in a single aliquot of SPAAC P–GFF reaction mixture (Fig. 3), with agitation to ensure solution homogeneity, identical temperature, concentration, and incubation time, still resulted in a large variation in Nyquist response (Fig. 3(a)). From this data, we concluded electrode modification was inconsistent, which is unsurprising considering extensive pre-modification mechanical and electrochemical polishing is necessary to achieve highly consistent self-assembled Au–thiol monolayer formation on gold disc electrodes.50 Such extensive cleaning of the gold surface is not possible when using Au-SPEs.
Fig. 3 (a) Nyquist plots of experimental data (dots) and the computational fitting to a modified Randles circuit (line) for EIS measurements of eight P–GFF Au-SPEs modified simultaneously. (b) MCMC analysis of the data in (a). The same colour code is used to distinguish between the separate electrodes. The uppermost plots (along the top diagonal) comprise 1-D histograms (y-axis is frequency) for each parameter used in equivalent circuit fitting (see x-axis labels), showing the distribution of the best-fit values for multiple analysis runs for each individual SPE. The 2-D scatter plots show correlations between parameters. (c) Pictures of the tank used for simultaneous Au-SPE modification. (d) A box and whisker plot showing the variation in the extracted Rct parameters for experiments on different ‘bare’, i.e. unmodified, Au-SPE (purple, n = 12) and different ‘P–GFF’ modified Au-SPEs (green, n = 63) following equivalent circuit fitting to the large datasets shown in Fig. S4 and S5.† The limits of the solid box represent the 25th and 75th percentile (20.8 and 44.9 kΩ, respectively), the central horizontal line in the box shows the position of the median (31.7 kΩ) and the upper and lower horizontal “whisker” lines represent the 5th and 95th percentiles (5.5 and 79.1 kΩ, respectively). EIS measurements were performed in the presence of 10 mM K3[Fe(CN)6] in pH 7 aqueous buffer (100 mM sodium phosphate, 233 mM sodium chloride) with the following parameters: t(equilibration) = 300 s, Edc = 0.230 V vs. Ref, Eac = 10 mV, f = 0.1 Hz–10 kHz. |
The Nyquist plots of the dataset comprised of eight simultaneously modified Au-SPEs all exhibit significant variation from one another (Fig. 3(a)) but could all be well modelled by the same modified Randles circuit. To explore if there were any underlying trends in the best fit parameters, we conducted a Markov Chain Monte Carlo (MCMC) analysis to generate posterior fitting distributions that represent the different equivalent circuit model parameters that can give a “best-fit” to the experimental data (Fig. 3(b)).45 Within this, we observed the expected correlation between Q and α, the coupled parameters which together constitute the constant phase element (CPE) used to model double-layer capacitance (diagonal lines in the scatter plot with Q on the x-axis and α on the y-axis).51 However, the plots showing the frequency distribution for each best fit parameter (the uppermost graphs within Fig. 3(b)) reveal no clear trends in the electrode-to-electrode variation; i.e. the ordering of the best fit parameter values with respect to the SPE number (essentially the colour bar code) is different from one parameter frequency distribution to another, meaning there is no clear underlying synergistic trend in the way one parameter changes versus another. Alternatively, this lack of pairwise parameter compensation can also be seen by the fact that there is no clear trend in how one parameter value changes relative to another in the 2-D scatter plots where one equivalent circuit model parameter is on the x-axis, and another is on the y-axis.45 Thus, the Rct value of a GFF-modified Au-SPE cannot be normalised with respect to another parameter to enable the trivial correction for variation in the dataset.
Since most EIS biosensors are assayed in a manner which correlates change in the charge transfer resistance (Rct) with target analyte concentration,2 the variation in the extracted Rct parameters from all the unmodified and P–GFF modified Au-SPEs datasets (Fig. S4 and S5†) was compared (Fig. 3(d)). The 12-fold increase in the average Rct for the bare electrodes (2.8 kΩ) versus the P–GFF dataset (34.2 kΩ) is equivalent to a 12-fold decrease in the apparent ferricyanide electron transfer rate constant, kapp, of 3.0 μm s−1versus 0.2 μm s−1.49 The observed variation between P–GFF modified electrodes was substantive, and this means it is non-trivial to compare repeats across multiple electrodes in experiments described later. However, it is notable that this variability does not change as a function of the age of the SPAAC P–GFF reaction mixture (Fig. 4(a)), with no correlation between Rct and the age (days post click reaction) of the SPAAC P–GFF reaction mixture. This dataset is comprised of extracted Rct values from different Au-SPEs modified with the same aliquot of SPAAC P–GFF reaction mixture and demonstrates the remarkable stability of the P–GFF conjugate, which yields self-assembled monolayers over an approximate 10 months time scale despite multiple freeze–thaw cycles. In addition, the formation and stability of the P–GFF conjugate was confirmed by comparison with the EIS response for Au-SPE modified with thiol-polymer only (Fig. S6†). In these ‘GFF-free’ control experiments a significantly larger Nyquist response is observed for the ‘polymer-only’ dataset (average Rct = 89.5 kΩ). This is attributed to differences in the interactions of the solution-based ferricyanide with polymer (P) versus polymer–glycofluoroform (P–GFF); the anionic charge of the Fe(CN)64−/3− redox couple makes it highly sensitive to changes in electrode surface electrostatics.52 The difference between the control experiment data and data from the modification of different Au-SPE with a single P–GFF aliquot (Fig. 4(a)) supports our assertion that the glycofluoroform modified polymer is stable over a prolonged timescale.
Fig. 4 (a) Extracted Rct parameters (black dots) for individual P–GFF modified Au-SPEs plotted against the age of the P–GFF conjugate solution showing no correlation in Rct over time. The mean Rct values (ct) for both bare (light purple dashed line) and polymer only, “P”, modified (orange dashed line) Au-SPEs are shown for comparison. (b) Distribution plots overlaying the datasets from Fig. 3(d) (purple and green), (a) (grey) and S6† (orange). Lines represent Kernel smooth distributions and a 50% gap is used to separate the vertical bars. The lower portion of the plot is a “rug” representing the point values of the individual measurements. |
To further assist in the visualisation of differences between datasets, Fig. 4(b) presents an overlay of the distributions of Rct values from repeat experiments analysing bare Au-SPEs, the full P–GFF dataset (both from Fig. 3(d)), the subset of P–GFF experiments carried out across many different days using the same polymer–glycofluoroform solution to modify the electrode (data from Fig. 4(a)), and polymer-only control experiments (Fig. S6†). Complementary t-test analysis (see ESI†) supports the conclusion that there is a statistically significant difference between the unmodified Au-SPE “bare” electrodes and the combined P–GFF dataset, but the inter-day variability is insignificant when this sub-dataset (N = 20) is compared to the mean of the collated P–GFF dataset from all experiments (N = 63). We did not make comparison to the polymer-only dataset because as shown in Fig. 4(b) we have insufficient data to resolve the skewed distribution in this relatively small (N = 5) set of control experiments, and polymer-only modifications are not the focus of this study.
Given the substantive variation shown between P–GFF modified Au-SPEs, we repeated the galectin-3 titration experiments on a total of 15 different P–GFF modified SPEs (Fig. S7†), noting a general trend of increasing Rct with increasing galectin-3 across a protein concentration range of 15–30 μg mL−1 (Fig. S8†). MCMC analysis was used to identify three electrodes which had sufficiently similar Rct values in the initial protein-free EIS experiments (12.9, 7.7, and 10.4 kΩ) to justify amalgamation (Fig. S8†). The Rctvs. galectin-3 concentration plot for this combined dataset is shown in Fig. S9(a),† and this confirms that a reproducible galectin-3 binding response is achieved across repeat electrodes where they have a similar initial Rct. Applying the analysis described by Bandyopadhyay et al.,53 assuming a Langmuir isotherm, enables the dataset to be linearised as shown in Fig. 6 with the gradient of the plot corresponding to 1/Kd as shown in eqn (1), where C is the galectin-3 concentration and Rct,0 and Rct,i are the resistance to charge transfer at a protein-free stage and after protein incubation, respectively. Using this analysis and setting a y-intercept of zero yields a Kd of 6.3 μg mL−1 which equates to 240 nM assuming a molecular weight of 26 kDa for galectin-3. The Kd value for the GFF EIS biosensor is notably higher than the Kd of 6.0 nM achieved using the same P–GFF binding system on Au nanoparticles,32 we attribute this to the differences in the grafting density and the increased steric repulsion in the interaction between solution-phase galectin-3 and a planar Au-SPE surface compared to a solution suspension of Au nanoparticles.
(1) |
Fig. 6 Combined datasets from three different P–GFF (green squares) and P-CS (pink dots) modified Au-SPEs titrated against galectin-3. The P–GFF dataset is fit to a straight line (black dashed line, eqn (1)) to enable determination of Kd, as described in the text. Inset: chemical structures of GFF (green) and CS (pink). The datapoints represent the average values from the combined datasets, while vertical error bars represent the standard deviation. |
An alternative method of assessing ligand-protein binding affinity is to use the Hill equation, as has been previously reported in the electrochemical biosensor literature.54–56 As shown in Fig. S9(b),† such an analysis can also be used to generate a good fit to the data. The best-fit Hill coefficient value, n = 5.7, is consistent with the notion that the proposed multimeric–multivalent binding mode of galectin-3 (illustrated in Fig. 1) supports cooperative binding between the protein and the P–GFF surface.35 The best-fit value for KA, the ligand concentration producing half occupation, was 18.5 μg mL−1 (equivalent to 713 nM).
Finally, an ideal diagnostic biosensor would show minimal off-target binding to other biomolecules present in biological samples; this can be very challenging to achieve because serum contains albumins and immunoglobins, and plasma additionally contains fibrinogen57 which are all proteins capable of non-specific surface-binding interactions. We therefore explored whether the incorporation of the P–GFF within the biosensor would act to advantageously reduce any non-specific binding, using albumin (BSA) as a model analyte. Importantly, although EIS experiments using increasing concentrations of BSA (Fig. 7) show a significant change in Rct on bare electrodes, this non-specific binding to the gold electrode is reduced when using electrodes modified with the P–GFF.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sc04409f |
This journal is © The Royal Society of Chemistry 2024 |