Quantification of protein interactions and solution transport using high-density GMR sensor arrays.

Monitoring the kinetics of protein interactions on a high-density sensor array is vital to drug development and proteomic analysis. Label-free kinetic assays based on surface plasmon resonance are the current gold standard, but they have poor detection limits, suffer from non-specific binding, and are not amenable to high-throughput analyses. Here, we show that magnetically responsive nanosensors that have been scaled to over 100,000 sensors per cm² can be used to measure the binding kinetics of various proteins with high spatial and temporal resolution. We present an analytical model that describes the binding of magnetically labelled antibodies to proteins that are immobilized on the sensor surface. This model is able to quantify the kinetics of antibody-antigen binding at sensitivities as low as 20 zeptomoles of solute.

Affinity-based sensing of DNA hybridization, antigen-antibody binding, and DNA-protein interactions play a vital role in basic science research, clinical diagnostics, biomolecular engineering, and drug design[1,2,3,4,5,6,7,8,9,10]. As the state of the art advances, demand for accurate, sensitive, specific, high-throughput, and rapid methods for the determination of molecular identities and reaction details places increasing pressure on evolving analytical methods[11,12,13,14,15,16,17]. To meet these pressing needs, researchers have turned to nano-scale labels in order to enhance the limit of detection (LOD) and specificity for detecting low abundance molecules. Such labels, however, can alter diffusion and steric phenomena. In addition, high-throughput, or speed requirements often prohibit the use of classical equilibrium methods, so a precise understanding of reaction kinetics, transport phenomena, and the implications of surface immobilization becomes critical for extracting meaningful molecular reaction parameters for nanoparticle labeling methodologies. This report addresses these issues and demonstrates that nanoparticle labeled proteins offer unique advantages over label-free methods, making this system very effective for modeling and extracting binding kinetics and analyte transport.

Extant modeling of molecular interactions has predominantly been restricted to label-free binding in solution. Early work by Berg and Stenberg proposed some of the first kinetic models of surface antigen-antibody interactions that explained the new restrictions that labeled reagents introduce on surface reaction kinetics by altering rotational and translational motion[18,19,20]. Furthermore, they argued that the use of targets immobilized on sensor surfaces implies that diffusion can become problematic due to the existence of long range concentration gradients, which can require ligands to traverse macroscopic distances (>100 μm) prior to reaction. Though many of these details are elaborated by Waite[21] and Sheehan[22], their emphasis on numerical methods precludes the derivation of semi-analytical expressions. While binding kinetics of quantum-dot-labeled macromolecules in liquid phase has been studied with fluorescence cross-correlation spectroscopy[23,24], we found no similar literature describing reactions on a sensor surface. Our investigation provides new quantitative insight into the binding kinetics of labeled macromolecules interacting with targets immobilized on a sensor surface, addressing this gap in the literature.

GMR nanosensor platform and magnetic nanoparticle tags

Our approach utilizes giant magnetoresistive (GMR) biosensors, an emerging tool for both basic science research and clinical diagnostics. Their superior LOD, multiplex capacity, broad linear dynamic range, and real-time readout capabilities make them ideal for kinetic analysis measurements[25,26,27]. GMR nanosensors, initially utilized as read head elements in computer hard drives, operate by changing their electrical resistance in response to changes in the local magnetic field[28,29,30,31]. Recent work has adapted GMR sensors for detection of biological species in solution by implementing a traditional sandwich assay directly on GMR nanosensors. If a magnetic particle is introduced to label the biomolecule of interest, GMR sensors are capable of highly sensitive DNA and protein detection[32,33,34,35,36]. This prior work[25,26] has involved quantifying the amount of protein, but has provided little information about the kinetics of the biomolecular reaction. In the current research, we pre-label the soluble ligand with a magnetic nanoparticle (MNP) in order to monitor the real-time binding kinetics of the ligand-MNP complex to antigens immobilized on sensor surfaces (Fig. 1a). As the antibody-MNP complexes are captured, their magnetic fields induce changes in electrical resistance in the underlying GMR sensor. Using the rapid, real-time readout of our GMR sensor array[25], we can monitor and quantify the kinetics of binding, thus determining the associated kinetic rate constants. Each GMR sensor in the array covers a total area of 100 μm × 100 μm and is comprised of twelve parallel GMR sensor stripes that are connected in series six times, producing a total of 72 stripes per sensor (Fig. 1b). Each stripe is 750 nm wide, approximately 20 nm thick, and spans 100 μm in length. Using scanning electron microscopy, it is possible to resolve nanoparticles bound over each sensor stripe (Fig. 1b insert).

Figure 1

GMR nanosensor and nanoparticle system for kinetic analysis

a, Schematic representation of antibody-antigen binding. On the left, antibody labeled with a magnetic nanoparticle tag in solution at concentration C approaches the GMR sensor surface. When not bound, most diffusing magnetically labeled antibodies are too far from the GMR sensor to be detected. Antigens are immobilized on the sensor surface at an initial surface concentration of n. Once the magnetically labeled antibody binds to the antigen, as depicted on the right, the magnetic field from the magnetic tag is detected by the underlying, proximity-based GMR nanosensor. The captured antibody-antigen complex surface concentration is n. b, Optical micrograph showing the GMR sensor architecture comprising 72 stripes connected in parallel and in series. Insert: scanning electron microscopy image of one stripe of the GMR sensor with several bound magnetic nanoparticle tags. c, Schematic representation of a magnetically labeled antibody drawn to scale. The magnetic tag is comprised of a dozen iron oxide cores embedded in a dextran polymer and then functionalized with antibody or receptor.

The MNPs that label the protein or antibody of interest are comprised of approximately twelve 10 nm iron oxide cores embedded in a dextran polymer (Fig. 1c), as determined by TEM analysis[37]. The entire nanoparticle averages 46 ± 13 nm in diameter (from number weighted Dynamic Light Scattering). Based on the Stokes-Einstein relation, these particles have a translational diffusion coefficient of approximately 9.3×10-12 m2 s-1. The MNPs have a reported zeta potential of -11 mV[38]. These particles are superparamagnetic and colloidally stable, so they do not aggregate or precipitate during the reaction. Importantly, the GMR sensors operate as proximity-based detectors of the dipole fields from the magnetic tags; thus, only tags within ∼150 nm of the sensor surface are detected[25]. Therefore, unbound MNP tags contribute negligible signal in the absence of binding, making this unique nanosensor-MNP system ideal for real-time analysis of kinetic binding process at the sensor surface.

Our GMR assay using MNP labels has several advantages over surface plasmon resonance (SPR) - the standard method of monitoring protein binding interactions, which operates by measuring changes in the refractive index of a thin film when unlabeled solute molecules bind to the surface[39,40]. While recent studies have attempted to design higher density SPR instrumentation[41,42,43], their methods typically monitor only a few reactions in parallel, and their LOD and dynamic range are limited to ∼ 25 ng mL-1 and ∼ 2 logs, respectively. In contrast, the GMR biosensor array can simultaneously monitor hundreds to thousands of sensors at sensitivities of ∼ 1 pg mL-1 or below and dynamic ranges of 6 logs or more[26]. Furthermore, we have fabricated GMR sensor arrays with 1,008 sensors on a 1 mm2 chip area[44]. The calculated feature density is now over 100,000 GMR sensors per cm2, the highest density reported to date for biosensor technologies. Each sensor within a sub-array is individually addressable by row and column decoders via a shared 6-bit control bus fabricated with state of the art VLSI technology. Such highly integrated GMR sensor arrays allow for exceptionally sensitive, massively parallel multiplex monitoring of protein binding kinetics.

Mathematical model of binding kinetics

Given our ability to monitor real-time binding kinetics using highly sensitive and multiplexed nanosensors, we first tested the widely used Langmuir absorption isotherm model, which assumes that there is negligible depletion of reactants and that the concentration of reactants near the surface is the same as that in the bulk. However, this did not adequately describe the binding kinetics in our system (Supplementary Fig. 1). A new model that would accurately describe the binding rates and would give an analytical solution is needed to explicitly describe the microscopic processes and offer insights into the mechanism of the overall reaction. We devised an analytical model capable of fitting the data with a high degree of generality across a variety of reaction conditions. This new model is capable of fitting real-time binding kinetics data for nanoparticle labeled macromolecules so that one can calculate the molecular association and dissociation rate constants, k and k.

Traditionally, protein binding is envisioned as a two-step process in which the target protein or antibody in the bulk solution first enters a surface compartment (reaction zone) via diffusion and flow, and then binds to or escapes from the immobilized targets via the chemical processes of association and dissociation (Supplementary Fig. 2)[45]. If the soluble ligand diffuses slowly, as for a large MNP tag, the diffusion into the reaction zone becomes negligible (see materials and methods and Supplementary Figs. 2-3 for complete derivation). Therefore, replenishment of soluble antibody-MNP via diffusion from the bulk compartment (bulk zone) to the surface compartment can be assumed to be zero in our model. Accordingly, the rate equation for the traditional two-compartment model can be reduced to the following (see materials and methods for derivation):

where n is the surface concentration of bound MNP-antibody-antigen complexes that have formed over the sensor, k is the association rate constant, k is the dissociation rate constant, C is the bulk concentration of magnetically tagged antibody, n is the maximum moles of surface bound complexes per area, V is the volume of the surface compartment and A is the reaction area (in this case it refers to the surface area of each sensor). The quantity of nA/V is equivalent to the bulk concentration in the reaction zone which would be consumed to produce the surface concentration, n, according to mass conservation. Conceptually, the two terms within the parentheses represent the depletion of tagged antibody in the reaction zone and reduction of available surface sites due to binding, respectively.

As both volume and surface concentrations are discussed, we must be aware of the different dimensions: n and n are expressed in mol m-2 whereas C and C are expressed in mol m-3. Note that n is limited by the maximum concentration of close packed MNPs, not by the maximum concentration of analyte on the surface. Therefore, in order to exclude steric effects related to MNP crowding on the sensor surfaces, the highest analyte surface concentration tested was significantly lower than the surface concentration of the close packed antibody-MNP complexes. In addition, by restricting the amount of protein deposited on the sensor surface to be under this limit, the proteins will be spaced across the sensor at distances greater than one MNP diameter. Therefore, each binding event monitored will be from one antibody on a MNP binding to one antigen on the sensor, and not from multiple antibodies binding to multiple antigens, mitigating avidity issues that may complicate the model.

We further simplify Equation 1 by assuming that k is zero, since antibody-antigen dissociation is typically negligible on the 400-1,000 second time scales of our experiments (Supplemental Fig. 3). Equation 1 now has the following analytical solution (derivation in online materials and methods):

If C ≫ nmaxA, which implies a vast excess of solution molecules over available surface sites, the kinetics reduce to Langmuir absorption dynamics (derivation in online materials and methods), demonstrating that the Langmuir model is merely a special case of our more general analytical model.

Experimental and theoretical binding kinetics measurements

To test this solution, we examined antibody-antigen binding kinetics for epithelial cell adhesion molecule (EpCAM) using this analytical model, and compared our experimental results to literature values. Anti-EpCAM antibody was selected because it was formulated into a chemotherapeutic drug, edrecolomab[46].

In the first set of binding experiments, presented in Figure 2a, we performed a binding assay of MNP-labeled anti-EpCAM antibody to surface-bound EpCAM protein. C, n, V, and A in the model were fixed values determined from dimensions and concentrations, and k was determined via best fit of the predicted binding curves to experimental data. We first performed a binding assay for varying concentrations of surface bound EpCAM protein at a fixed concentration of MNP-anti-EpCAM antibody complexes. Twofold dilutions were used to prepare a series of sensor surfaces, beginning at a loading mass (i.e., the amount of protein that is actually bound to the sensor surface and functional) of 5 attomoles of EpCAM and diluting sequentially down to 20 zeptomoles. Therefore, the only parameter that varied between binding curves was n; all other parameters were unchanged. When this one parameter variation is implemented in the model, each experimental binding curve was fitted accurately (Fig. 2a). The values of the parameters (the undiluted case) were n = 9.5 × 10-10 mol m-2, C0 = 6.8 × 10-7 M, A = 5.4 × 10-9 m2 and V = 5.5 × 10-12 m3. Accordingly, k = 2.5 × 104 M-1 s-1 fit the data best (R2 = 0.98). In addition, after the MNP-antibody solution was washed away and replaced by antigen-loaded buffer, dissociation constants were calculated by fitting the subsequent data to a basic exponential decay model, n(t) = n0e, where n is the surface concentration of bound MNPs at the time of washing (derivation in online materials and methods). Accordingly, the anti-EpCAM antibody-antigen dissociation constant, k, was determined to be 2.0 × 10-6 s-1, supporting our assumption that k is essentially negligible compared to k.

Figure 2

Comparison of experimentally generated binding curves to kinetic model predictions

a, Binding curves for anti-EpCAM antibody binding to EpCAM antigens immobilized on the surface. Dotted lines are predictions using the analytical model in equation 2 (•••) and full lines are experimental data (–) obtained for surface loading amounts that vary from 5 attomoles (n) to 20 zeptomoles (n /256) in serial dilutions of 2×. The fitting error of all the curves in this experiment to the curves predicted by the model is R2=0.98. b, Binding curves for MNP-anti-EpCAM antibody binding to 833 zeptomoles (n/6) of EpCAM antigen immobilized on the sensor surface. Dotted lines show prediction using the analytical model and full lines are experimental data obtained for MNP-anti-EpCAM that are undiluted, diluted 2× and diluted 8×. The fitting error of all the curves to the model is R2=0.96. The y-axis is presented as changes in magnetoresistance (MR) normalized to the initial MR in parts per million (ppm).

In a second set of binding experiments, presented in Figure 2b, each sensor was immobilized with a constant load mass of 833 zeptomoles of EpCAM protein (1/6 of the maximum amount used in Figure 2a). The concentration of antibody-MNP complex applied to the sensors was varied between undiluted, twofold diluted, and eightfold diluted solutions of antibody-MNP complexes (corresponding to C0, C0/2 and C0/8 in the model). Since the antibody-antigen interaction should remain the same whether C0 or n is altered, the rate constants describing the interaction above should remain the same across these diverse experiments, as was observed. In fitting the model, we obtained the same k of 2.5 × 10[4] M-1 s-1 (R2=0.96), supporting the validity of our analytical model. Furthermore, these results lie within the expected range reported in the literature, confirming both the validity of our kinetic model to predict binding and the accuracy of the results derived (Table 1).

Table 1

Comparison of k for binding of biotin to streptavidin, binding of EpCAM antibody to EpCAM antigen, binding of CEA antibody to CEA antigen and binding of VEGF antibody to VEGF antigen when using the GMR sensor array and SPR. For both SPR and GMR sensor experiments, the exact same antibody pairs were used. In addition, each method of kinetic analysis was consistent with the literature.

	GMR sensor(×104 M-1 s-1)	SPR(×104 M-1 s-1)	Literature(×104 M-1 s-1)
Biotinylated DNA and streptavidin	467	550	300 – 4500[48]
EpCAM antigen and antibody	2.5	N/A	3.2 – 40[49]
CEA antigen and antibody	5.0	5.2	3.7 – 11[50]
VEGF antigen and antibody	1.6	N/A	0.5 – 7[51]

Similar real-time experiments were performed in order to quantify binding kinetics for MNP-anti-carcinoembryonic antigen (CEA) antibody to CEA, MNP-anti-vascular endothelial growth factor (VEGF) antibody to VEGF, and MNP-streptavidin to biotin binding kinetics (Supplementary Fig. 4). CEA and VEGF were chosen because they are among the most well-known clinical tumor markers, and anti-VEGF antibody drugs, such as bevacizumab (Avastin; Genentech/Roche), are highly effective anti-cancer drugs[46,47]. For anti-CEA antibody, binding and dissociation constants were monitored with the same reagents on both GMR sensors and SPR instruments, and the results were compared (Supplementary Figs. 4c,d). The GMR sensor array and SPR measurements yielded kinetic constants that were consistent with one another and with literature values (Table 1).

Scanning electron microscopy (SEM) confirmed that our analytical system also has the ability to quantify the precise number of proteins captured on each sensor. Thus, by calibrating GMR signal to an absolute number of magnetic tags bound to the sensor surface, we can derive the mass of protein bound and the signal generated per MNP. To accomplish this, a combination of the prior two experiments was performed and accurately modeled where EpCAM protein was serially diluted in twofold increments, starting at 2.5 attomoles and ending at 78 zeptomoles, on three to eight replica sensors (Fig. 3 left panel). A 20-fold dilution of the MNP-antibody complexes was subsequently added to all the sensors and the binding kinetics were monitored. After 20 minutes of incubation time, the solution of magnetically labeled antibody was washed away to terminate the binding reaction, at which point the sensors were imaged by SEM (Fig. 3 right panel). By normalizing the real-time experimental data and fitting it to the model, we were able to convert the sensor signal, measured as a change in magnetoresistance (MR) normalized to the initial MR (ΔMR/MR0), into the number of magnetic tags bound to each sensor. For example, the sensor functionalized with 2.5 attomoles of EpCAM protein captured 190,000 MNP tags within the time of the experiment. SEM imaging reveals that experimental results match with predictions from the kinetic model, thus extending the validity of the model for precisely quantifying the number of tags bound per sensor and determining the number of proteins that bind during a given reaction.

Figure 3

The kinetic model can predict the number of protein binding events

By fitting real-time binding curves to the model, it is possible to convert the signal generated from the GMR sensor into an absolute number of magnetic tags bound to the sensor surface. Here, EpCAM protein was loaded onto the sensors at a mass of 2.5 attomoles (n) and at masses serially diluted in twofold increments down to 78 zeptomoles (n/64). At least three replica sensors were used for each dilution. After 20 mins incubation with the 20-fold diluted solution of MNP-anti-EpCAM antibody (C0/20), the solution was washed away to terminate the binding reaction. Subsequently, a small section of the sensors was imaged with SEM (right panel; colour-coded boxes represent the different loading masses) to compare the number of MNPs bound in the experiment with that predicted by the model (left panel). The number of MNPs presented above each binding curve in the left panel represents the number of MNPs predicted to have bound to each corresponding sensor. Dotted lines are predictions using the analytical model in equation 2 (•••) and full lines are experimental data (–). SEM image of nmax/64 was not shown due to the low surface coverage of MNPs.

Furthermore, using our models, we estimate that every 150 MNPs produce ∼1 ppm of the normalized signal. The LOD of our sensors is ∼20 ppm (defined by the average background signal of a non-complementary antibody coated sensor plus two standard deviations). Therefore, we can detect as few as 0.6 particles μm-2. In addition, the model is capable of explaining when saturation of the sensor surface will occur. For example, when 10 attomoles of protein are deposited on the sensor surface, experimental data shows that the sensor surface is approaching saturation as judged by comparing signals at 5 attomoles and lower (Supplementary Fig. 5). A loading mass of 10 attomoles is described in the model as n being equal to 1.9×10-9 mol m-2. The saturation value predicted by the model is thus in close agreement with the maximum surface concentration of the MNPs in a close packed single layer, 1.0 × 10-9 mol m-2.

Our real-time binding assay and kinetic model has been extended to visualize protein binding events in both space (due to the high density of the array architecture) and time (due to the rapid and real-time readout). As an example, we first immobilized the same anti-CEA capture antibody across a sensor array. We then added MNP-anti CEA conjugates in solution to the sensor well. Upon delivery of the CEA antigen, we were able to monitor the movement of CEA protein across the sensor array (Fig. 4a). With this method, we can visualize the binding of soluble CEA protein to each sensor, which is spatially distributed in the array, by means of the MNP-anti-CEA antibody binding (Fig. 4b and Supplementary Fig. 6). In this manner, each sensor monitors the reaction zone above it, and thus we can investigate protein binding, protein transport, and protein dynamics with high spatial and temporal resolution--an important and unique development. Such techniques will be applied to monitor localized cell-cell communications via cellular protein secretome analysis in future works.

Figure 4

Visualization of spatiotemporal resolution of the sensor array

a, Schematic depicting the GMR sensor array functionalized with monoclonal anti-CEA capture antibody (not to scale). The solution above the sensor array is comprised of magnetically labeled anti-CEA detection antibodies. The schematic includes a pipette tip containing a solution of CEA protein prior to injection. b, Once the CEA antigen is introduced into the solution above the sensor array, the transport of CEA antigen from the near center of the array radially was monitored in real-time as magnetically labeled detection antibodies bound to captured CEA protein, forming detectable sandwich structures. c, Visualization of CEA protein surface concentration at different times using a high density GMR sensor array. The units of the y-axis are presented in changes in MR normalized to the initial MR in parts per million (ppm).

Discussion and conclusions

We have presented a high-density, highly sensitive, real-time binding assay for quantifying protein binding kinetics and analyte transport at the surface of a biosensor array. Notably, we have also developed a novel analytical kinetics model that provides a precise and physically intuitive description of the dominant processes involved in labeled protein-protein interactions. We proved that the reduced rate of diffusion, that is consistent with the addition of nano-scale labels, can be exploited to derive generalizable kinetic binding models with analytical solutions. The combination of the GMR sensor technology and analytical model has enabled us to measure protein binding constants and quantify the number of proteins bound to a given sensor with unprecedented LOD. Furthermore, we have demonstrated the unique benefits of combining our GMR sensor array, kinetic binding assay, and analytical model to visualize protein transport in a reaction well.

The tools we have developed can be used broadly in basic science research for understanding receptor-ligand binding interactions involved in signal transduction in cell biology or for profiling the affinity of specific compounds of interest against an entire proteome on a high density array. In addition, applications to pharmacodynamics via receptor occupancy assays for drug development would open an entirely new and broad segment for our technology. MHC-peptide and T-cell receptor interactions could also be studied with exceptional speed and accuracy.

Further, the potential clinical applications of our method are vast, ranging from in-vitro clinical assay development[52], to targeted molecular imaging for early cancer diagnostics[53], to investigating drug on-target and off-target cross-reaction[54] binding kinetics. In each case, the affinity and cross-reactivity of reagents and targeting probes are critical issues to clinical adoption and can be easily addressed with our platform. In short, our method represents a significant advance in this domain as it has the potential to provide fundamentally important solutions for both cutting-edge basic science research and clinical practice.

Materials and Methods

GMR nanosensor array architecture

The giant magnetoresistive (GMR) sensor used in our experiment has a bottom spin valve structure of the type: Si/Ta(5)/seed layer/IrMn(8)/CoFe(2)/Ru/(0.8)/CoFe(2)/Cu(2.3)/ CoFe(1.5)/Ta(3), all numbers in parenthesis are in nanometers. Each chip contains an array of GMR sensors, which are connected to peripheral bonding pads by a 300 nm thick Ta/Au/Ta lead. To protect the sensors and leads from corrosion, two passivation layers were deposited by ion beam sputtering: first, a thin passivation layer of SiO2(10 nm)/Si3N4(20 nm)/ SiO2(10nm) was deposited above all sensors and leads, excluding only the bonding pad area; then a thick passivation layer of SiO2(100 nm)/Si3N4(150 nm)/SiO2(100 nm) was deposited on top of the reference sensors and leads, but excluding the active sensors and bonding pad area. The magnetoresistive ratio is approximately 12% after patterning. The pinning direction of the spin valve is in-plane and perpendicular to the sensor strip. The easy axis of the free layer is set by the shape anisotropy to be parallel with the sensor strip. This configuration allows the GMR sensors to work at the most sensitive region of their MR transfer curves.

Due to the GMR effect, the resistance of the sensor changes with the orientation of the magnetization of the two magnetic layers, which are separated by a copper spacer layer:

Here R0 is the resistance under zero magnetic field, δRmax is the maximum resistance change and θ is the angle between the magnetization of the two magnetic layers. In the bottom spin valve structure, the magnetization of the bottom magnetic layer (pinned layer) is pinned to a fixed direction, while the magnetic orientation of the top magnetic layer (free layer) can freely rotate with the external magnetic field. As a result, the stray field from magnetic labels can change the orientation of the free layer magnetization and therefore change the resistance of the sensor.

Magnetic labels

The magnetic labels we used are commercial magnetic nanoparticles from Miltenyi Biotech Inc., referred to as “MACS”. Each MACS particle is a cluster of 10nm Fe2O3 nanoparticles held together by a matrix of dextran. Due to the small size of the Fe2O3 nanoparticles, the MACS particle is superparamagnetic, with a hydrodynamic diameter of ∼46 nm, and contains ∼10% magnetic material (wt/wt). MACS particles can be functionalized with specific affinity molecules corresponding to the analyte being studied. For EpCAM (CEA, VEGF) experiments, the MACS particles are functionalized with anti-EpCAM (CEA, VEGF) antibodies; for biotin-streptavidin experiments, the MACS particles are functionalized with streptavidin.

Surface chemistry

The sensor surface was first cleaned using acetone, methanol and isopropanol and subsequently exposed to oxygen plasma for three minutes. A 2% (w/v) polyallylamine solution in Milli-Q water was applied to the sensor for 5 minutes. The chips were then rinsed with Milli-Q water and baked at 150°C for 45 minutes. A 10% (w/v) solution of EDC and 10% (w/v) solution of NHS were then added to the sensor surface at room temperature for 1 hour, after which the sensor was rinsed and dried. Capture protein EpCAM (960-EP-050 from R&D Systems), CEA (4128-CM-050 from RD Systems) or VEGF (293-VE165 from R&D Systems) or capture antibody (antibody to EpCAM (ab20160 from Abcam or 960 from R&D), CEA (5910 from BiosPacific) or VEGF (ab69479 from Abcam) was robotically deposited (Scienion sciFlexarrayer from BioDot) over each sensor in 360 picoliter droplets 3 times (total volume of ∼1 nanoliter). In order to monitor reproducibility, from 3 to 8 sensors, randomly distributed across the GMR sensor array, were incubated with the same capture protein. The control sensors were immobilized in a similar fashion with either BSA at 1 mg/mL or a non-complementary antibody (typically anti-survivin antibody (H00000332-P01 from Novus Biologicals, LLC) at 500 μg/mL. Finally, the entire surface of the sensor array was blocked with 1 mg/mL of BSA in phosphate buffered saline for 30 minutes.

Kinetic Assay

After the sensor surface is functionalized with the appropriate capture protein, the GMR sensor array was placed in the test station. The BSA blocking buffer was washed away and a 50 μL solution of the magnetically labeled detection antibody (described above) was added to the reaction well. The GMR sensor array was monitored in real-time as the magnetically labeled detection antibodies bound to the corresponding surface-immobilized proteins. The binding curves, unique to each protein, were then plotted and the binding constants could be determined. An assay is typically run for approximately 5 minutes.