The invention relates to a method and a microelectronic sensor device for the determination of the amount of target particles in a sample, wherein the amount of target particles in a sensitive region is measured. Moreover, it relates to a magnetic sensor device for detecting magnetized particles.
From the WO 2005/010543 A1 and WO 2005/010542 A2 (which are incorporated into the present application by reference) a microelectronic magnetic sensor device is known which may for example be used in a microfluidic biosensor for the detection of molecules, e.g. biological molecules, labeled with magnetic beads. The microsensor device is provided with an array of sensor units comprising wires for the generation of a magnetic field and Giant Magneto Resistances (GMR) for the detection of stray fields generated by magnetized beads. The signal of the GMRs is then indicative of the number of the beads near the sensor unit. A problem of these and similar biosensors is that the concentration of the target substance is typically very low and that the measurement signals are therefore severely corrupted by different sources of noise. Moreover, the measurement signals are very sensitive to variations in the parameters of the read-out electronics, for example the sensitivity of the sensor unit.
Based on this situation it was an object of the present invention to provide means for improving the accuracy, robustness and/or signal-to-noise ratio of microelectronic sensor devices of the kind described above, particularly of magnetic biosensors, wherein these means shall preferably work for different concentrations of target substance.
This objective is achieved by a microelectronic sensor device according to claim 1, a method according to claim 2, and a magnetic sensor device according to claim 20. Preferred embodiments are disclosed in the dependent claims.
A microelectronic sensor device according to the present invention is intended for the determination of the amount of target particles in a sample. The target particles may for instance be biological molecules like proteins or oligonucleotides, which are typically coupled to a label like a magnetic bead or a fluorescent molecule that can readily be detected. The “amount” of the target particles may be expressed by their concentration in the sample, and the sample is typically a fluid, i.e. a liquid or a gas. The microelectronic sensor device comprises the following components:
The invention further relates to a method for the determination of the amount of target particles in a sample provided in a sample chamber, wherein the method comprises the following steps:
The microelectronic sensor device and the method described above have the advantage that their determination of the amount of target particles in the sample is based on a plurality of measurement signals that were consecutively sampled during some observation period. The determination can thus exploit a redundancy to achieve a higher accuracy than the single measurements that are usual in the state of the art. Moreover, an estimation of the measurement error can be provided by a statistical analysis of the sampled measurements.
In the following, preferred embodiments that apply to both the microelectronic sensor device and the method defined above will be described.
In a first preferred embodiment of the invention, the sensitive region comprises specific binding sites for the target particles. The sensitive region may for example be a part of the walls of the sample chamber that is coated with hybridization probes which can specifically bind to complementary biological target molecules. Thus target particles of interest can selectively be enriched in the sensitive region, making the measurement specific to the target particles and increasing the amplitude of the measurement signals.
In a further development of the aforementioned approach, the measurement signals that are provided by the at least one sensor unit are indicative of the amount of target particles bound to the binding sites. This can for example be achieved by making the sensitive region small enough such that it substantially comprises only a volume in which target particles can only be if they are attached to a binding site. Alternatively, the amount of free (unbound) target particles within the sensitive region—which also contribute to the measurement signals—may be estimated and subtracted from the whole measurement signal to determine the amount of bound target particles. Finally, it is possible to remove unbound target particles from the sensitive region by some washing step (e.g. a fluid exchange or a magnetic repulsion of free target particles) to make the measurement signals only depend on the bound target particles.
In another variant of the aforementioned embodiment, a parametric binding curve is fitted to the sampled measurement signals, wherein preferably at least one of the fitted parameters is directly indicative of the amount of target particles in the sample. The binding curve can for example be provided by theoretical models of the binding process or simply be taken from general purpose functions for curve fitting (e.g. polynomials, sine curves, wavelets, splines, etc.). As the amount of target particles in the sample obviously has a critical influence on the binding kinetics, the binding curve will particularly reflect this value that is to be determined.
A particularly important realization of the aforementioned approach comprises the application of a Langmuir isotherm as a binding curve, which describes a large variety of different binding processes.
The fitting of the parametric binding curve, i.e. the adjustment of its parameters, can in general be achieved by any method known for this purpose from mathematics. Preferably, the fitting is achieved by a linear or a weighted least squares regression. In a weighted least squares regression, the weights may for example be determined by the expected or theoretical noise level which normally goes with the square-root of the number of particles.
A central aspect of the approach described above is that the amount of target particles in the sample is determined from a series of measurement signals, wherein the redundancy of these measurements is used to improve the accuracy of the final result and to provide an error estimation. According to a further development of the invention, the series of measurement signals is further exploited to adjust dynamically (i.e. during the ongoing sampling process) the configuration and parameter settings of the measurement device for improving the signal-to-noise ratio of the final results. One particularly important example of a parameter that can dynamically be adjusted is the sampling rate, i.e. the frequency with which measurement signals indicative of the amount of bound target particles are generated by the sensor unit. A further parameter of particular importance is the size of the sensitive region. As this size has opposite effects on different kinds of noise, there exists an optimal value for which the generated noise is minimal.
In a preferred embodiment of the invention, the sampling rate is adjusted such that it is of the same order as or larger than the binding rate of target particles to binding sites in the sensitive region (i.e. larger than about 5% of the binding rate). Said binding rate describes the net number of target particles that are bound to the sensitive region per unit of time. Making the sampling rate as large as the binding rate or a larger guarantees that in the mean each binding event will be captured by the measurement signals, thus providing complete information about the binding process.
In the aforementioned embodiment, the sampling rate can be adjusted once at the beginning of the sampling process. The determination results can however be improved if the binding rate is estimated during the sampling process from the momentarily available measurement signals and if the sampling rate is dynamically adjusted according to these estimations of the binding rate. Thus it is possible to start a sampling process without any previous knowledge about the sample and the amount of target substance therein and to improve in one or more steps the sampling rate as one crucial parameter of the process based on the most recently available information.
The size of the sensitive region may optionally be adjusted based on a given value of the sampling rate, wherein said adjustment is typically done such that the theoretically or empirically determined signal-to-noise ratio is optimized. The given value of the sampling rate may for example be determined before the sampling process starts or dynamically during the ongoing sampling process according to the principles described above. The size of the sensitive region may then accordingly be adjusted once at the beginning of the sampling process or dynamically during this process based on the most recent values of the sampling rate.
A preferred way to adjust the size of the sensitive region is by functionally coupling various numbers of sensor units to one “super-unit”.
As was already mentioned, the sensor unit may particularly be adapted to measure magnetic fields. In a preferred embodiment of this variant, the sensor unit comprises at least one magnetic sensor element for measuring magnetic fields, wherein said sensor element may particularly comprise a coil, a Hall sensor, a planar Hall sensor, a flux gate sensor, a SQUID (Superconducting Quantum Interference Device), a magnetic resonance sensor, a magneto-restrictive sensor, or a magneto-resistive element like a GMR (Giant Magneto Resistance), a TMR (Tunnel Magneto Resistance), or an AMR (Anisotropic Magneto Resistance) element.
The sensor unit may further comprise at least one magnetic field generator for generating a magnetic excitation field in the sensitive region. Thus magnetic entities (e.g. target particles comprising magnetic beads) may be magnetized in order to detect their presence by excited reaction fields.
In a further development of the microelectronic sensor device and/or the method of the present invention, the measurement signals that are provided by the sensor unit are indicative of “events” that are by definition related to the movement of (at least) a limited number of target particles into the sensitive region, out of the sensitive region and/or within the sensitive region. Preferably, the limited number is “one”, i.e. the measurement signals can resolve events related to the movement of single target particles. The detection of events caused by single or a few target particles provides insights into the microscopic behavior of the system under investigation that can favorably be exploited to determine the amount of target particles in the sample. Particular embodiments of this approach will be described in more detail in the following.
Thus the evaluation unit may for example be adapted to detect and count the events indicated by the measurement signals. Detection of an event in a (quasi-) continuous measurement signal may for example be achieved via matched filters that are sensitive to the specific signal shapes of the events. Counting the detected events, which can readily be realized by e.g. a digital microprocessor, will then provide data that are directly related to the amount of target particles in the sensitive region. If the counted events correspond for example to the entrance of single target particles into or their escape from the sensitive region, the total number of target particles inside the sensitive region can be determined by observing the process from the beginning on, starting with a sensitive region free of target particles. The great advantage of this counting approach is that the detection of events is very robust with respect to variations in e.g. the sensor electronics, because an event can reliably be recognized even if its particular shape varies in a broad range. This is comparable to the high robustness of digital data encoding and processing with respect to analog procedures.
The evaluation unit may preferably be adapted to determine the changing rate and/or the amplitude step in the measurement signals that are associated with an event. The amplitude step obviously comprises information about the number of target particles that enter or leave the sensitive region. The changing rate with which such an amplitude step takes place may provide valuable information, too, because it is related to the movement velocity of the target particles. The determination of the changing rate may thus for example allow to determine the average velocity of the target particles in the sample.
According to another embodiment, the evaluation unit may be adapted to discriminate between events that correspond to the movement of single target particles and the movement of clustered target particles, respectively. The clustering of target particles, particularly particles labeled with magnetic beads, is often an undesired but unavoidable process taking place in a sample. The clustered target particles usually deteriorate the measurement results. A cluster of e.g. four target particles that is bound to one binding site may for example wrongly be interpreted as four single target particles occupying four binding sites. The accuracy of the measurement results may therefore be improved if the effects caused by clusters can be discriminated from the effects of single particles. Such a discrimination between single and clustered target particles may in the described embodiment for example be achieved based on differences in their movement velocity, which is typically larger for the clusters.
The evaluation unit may further be adapted to determine the amount of unbound target particles in the sensitive region from events corresponding to target particles entering and/or leaving the sensitive region. The target particles that are free to move, i.e. not fixed to binding sites in the sensitive region, will usually follow a random walk due to their thermal motion. The rate with which such target particles cross the interface between the sensitive region and the residual sample chamber depends on the amounts of target particles on both sides of said interface (or, more specifically, their concentrations). Detecting events of interface crossings will thus allow to estimate said amounts.
The invention further comprises a magnetic sensor device with an electrically driven magnetic sensor component for detecting magnetized particles in an associated (one-, two-, or three-dimensional) sensitive region, wherein the size of said sensitive region can dynamically be adjusted. In this context, “dynamical adjustment” is to be understood as a change of the sensitive region that can be made (and reversed) at arbitrary times by external commands or inputs; the term shall particularly distinguish the adjustments meant here from changes of the design at the time of the production of the magnetic sensor device or from physical reconstructions of the device, which are of course always possible. Moreover, it should be noted that the magnetic sensor component by definition needs the electrical energy it is driven with to provide measurement signals indicative of the detected magnetized particles.
The dynamical adjustment of the sensitive region allows to tune a parameter that has turned out to have a crucial influence on the detection of magnetized particles. The positive effects of this approach will be described in more detail in the following with respect to specific embodiments of the magnetic sensor device.
In general, there are many possibilities to change the size of the sensitive region in a magnetic sensor device of the kind described above. In a preferred realization, the magnetic sensor component comprises a plurality of magnetic sensor elements that can selectively be coupled in parallel and/or in series. By coupling different numbers and/or configurations of individual magnetic sensor elements to one “super-unit”, the resulting sensitive region, which is composed of the individual sensitive regions of all coupled magnetic sensor elements, can stepwise be adapted as desired. A change of the sensitive region can thus be achieved by a reconfiguration of the network of coupled magnetic sensor elements, for instance by closing/opening appropriate switches.
According to a further development of the aforementioned embodiment, the magnetic sensor elements can selectively be coupled in such a way that a predetermined distribution of coupled magnetic sensor elements is achieved in a given investigation region, wherein said distribution is preferably homogenous. Thus a whole investigation region can be covered with effectively different sizes of sensitive regions.
In a further development of the invention, the magnetic sensor device comprises an electrically driven magnetic field generator for generating a magnetic (excitation) field in an associated excitation region, wherein the size of said excitation region can dynamically be adjusted. The magnetic field generator uses the supplied electrical energy to generate the magnetic excitation field, which is preferably used to magnetize particles which shall thereafter be detected by the magnetic sensor component.
While there are again many possibilities to realize the dynamically adjustable excitation region, it is preferred that the magnetic field generator comprises a plurality of individual magnetic excitation elements that can selectively be coupled in parallel and/or in series. Furthermore, these magnetic excitation elements can preferably be coupled such that a predetermined (preferably homogenous) distribution of coupled magnetic excitation elements is achieved in a given investigation region.
In general, the sensitive region associated to the magnetic sensor component and the excitation region associated to the magnetic field generator may be separate. Preferably, these regions will however partially or completely overlap.
The adjustment of the sensitive region or the excitation region may be exploited for different purposes. Preferably, the size of the sensitive region and/or the size of the excitation region is adjusted such that the signal-to-noise ratio of the magnetic sensor device is optimized, as analysis shows that this ratio is significantly influenced by the size of said regions.
Moreover, the size of the sensitive region and/or the size of the excitation region may be adjusted such that a predetermined ratio between thermal (i.e. temperature dependent) noise and statistical noise (i.e. noise caused by the magnetized particles) is achieved in the overall signal of the magnetic sensor component, wherein said ratio optionally may vary between 80% and 120% of its nominal value. As will be described in more detail with reference to the Figures, the ratio of the noises typically has a crucial influence on the signal-to-noise ratio.
The magnetic sensor component may particularly comprise a coil, a Hall sensor, a planar Hall sensor, a flux gate sensor, a SQUID (Superconducting Quantum Interference Device), a magnetic resonance sensor, a magneto-restrictive sensor, or a magneto-resistive element like a GMR (Giant Magneto Resistance), a TMR (Tunnel Magneto Resistance), or an AMR (Anisotropic Magneto Resistance) element.
In a particular embodiment of the invention, the magnetic sensor device comprises an alternating sequence of resistances functioning as magnetic excitation element and magnetic sensor component, respectively. It may for example consist of a sequence “wire-GMR-wire-GMR- . . . ”, wherein the wires are individually addressable magnetic field generators and the GMRs are individually addressable sensors.
These and other aspects of the invention will be apparent from and elucidated with reference to the embodiment(s) described hereinafter. These embodiments will be described by way of example with the help of the accompanying drawings in which:
Like reference numbers or numbers differing by integer multiples of 100 refer in the Figures to identical or similar components.
The concentration of the target particles 2 which has to be measured can be very low, depending on the biochemical application. To reach a detection limit which is as low as possible, the sensor geometry, electronics and detection algorithms have to be optimized. Furthermore, preferably the device should be able to detect different kinds of target particles which requires multiple sensors onto a single die.
In the following it will first be shown that the signal-to-noise ratio (SNR) of a magnetic biosensor can be optimized by optimizing the size of its sensitive region, i.e. the “sensor area”, as different noise sources scale differently with sensor area. In the presented analysis, the SNR will be the performance indicator for which the optimization is carried out, and constant power dissipation will be assumed during the optimization process, because typically the total power dissipation is limited by temperature and battery lifetime considerations. Moreover, the scaling of the sensor area is discussed by describing the effect of combining multiple sensor units (e.g. the sensor units 10a to 10d of
To determine how the SNR scales with sensor area, the scaling effects on the sensor signal and the main noise sources will be discussed first.
The complete circuit of
The sensor signal S provided by an individual sensor element can be expressed as in equation (3), where Isense is the current through the sensor element, ssense is the sensitivity of the sensor element (dR/dH)H-0/R, Rsense is the resistance of the sensor element, Iexc is the current through the associated excitation element, nbead is the number of beads on the associated area of the sensor element, and Xbead is the magnetic susceptibility of a single bead.
In the same way the signal change S′ of the series/parallel-connected network can be expressed by equation (4). The factor 1/m expresses the reduction in the excitation current due to the distribution of the current over the series/parallel network. By substituting equations (1) and (2), the signal S′ can be expressed in terms of the signal S.
The thermal noise power, Nth2, of an individual sensor element can be expressed as in equation (5), where k is the Boltzmann constant, T is absolute temperature, and B is the bandwidth. The thermal noise power scales directly with the total resistance of the magnetic sensor component; for a network consisting of series and parallel connected units the thermal noise power can therefore be expressed as in equation (6).
There are a few other noise sources that also lead to variations in the sensor signal:
Since these noise sources scale equally with sensor area, they will be treated here all at once. The statistical noise power of a single sensor element, Nstat2, translates to a statistical noise contribution of the series/parallel-connected network as expressed in equation (7). The uncorrelated variation of all n times m sensor units in the total network therefore works out as in equation (8).
These statistical noise sources scale with the sensor signal per network element, therefore the noise contribution needs to be multiplied by the scaling of the currents Isense and Iexc per element, cf. equation (9).
The overall signal-to-noise ratio SNR′ can then be expressed as in equation (10). From this expression two very important conclusions can be drawn:
Expression (10) shows that for the SNR it does not matter whether multiple elements are connected in series or in parallel. The choice between series and parallel can thus be made in accordance with the read-out electronics.
The statistical noise is a function of the sensor signal and therefore its value changes with the bead concentration on the surface of the sensor. The thermal noise is constant in time. Therefore, the optimal sensor area is a function of the concentration of bound target: For large concentrations the signal is much larger than the thermal noise. By increasing the area (increasing n×m), the signal is reduced in favor of a better statistics.
In summary, it has been derived that the optimal sensor area can be optimized for the bead concentration on the sensor surface. However, there are situations where this bead concentration is not always the same. Thus different target concentrations will lead to different concentrations of bound beads at the sensor surface. For each concentration there is an optimal sensor area. To get optimal performance one should use a differently sized sensor for each target concentration. This is not very practical. What makes it even harder is that typically the target concentration is not known beforehand.
As will be derived below, it is advantageous to continuously measure the sensor signal during the binding process of the beads. This means that the concentration of beads on the sensor surface continuously increases over time. To maintain the optimal SNR during the experiment, the sensor area needs to scale with time.
To perform optimal measurements with a magnetic biosensor under these circumstances a sensor is required of which the (active) sensor area can be adapted dynamically. This can be realized by splitting up the entire sensor area into multiple blocks. Depending on the concentration of beads on the surface, one or multiple sensor blocks can be read-out. When the target concentration on the sensor surface increases over time, the optimal SNR can be maintained by distributing the total power over ever more sensor blocks.
Based on the above observations, a signal analysis method will be described in the following which increases the signal-to-noise ratio of the sensor device such that lower concentrations of target particles can be detected, decreases the required area of the sensor device, allowing more sensors onto one die, thus allowing a larger variety of substances which can be measured simultaneously, and makes the sensor design independent of the concentration of target particles.
In the sensor units 10a, 10b, . . . of a magnetic sensor device like that of
In order to detect one particular kind of target particle—in the following without loss of generality assumed to be a protein 2—, the surface of the sensor device is prepared with species (anti-bodies) such that only one particular kind of protein can attach, i.e. the binding or adsorption sites 3 are specific for the protein 2 of interest. In an unused sensor device, no magnetic beads will be detected by the sensor units since no proteins are yet present. Once the sample solution to be analyzed is brought into the sample chamber 1 and comes into contact with the sensor surface, the proteins 2 with magnetic label start reacting with the prepared sensitive region 14. With increasing time, more proteins 2 will be bound to the surface 14 and the sensor signal will increase in time. The rate at which the signal increases in time is dependent on the concentration of the proteins 2 in the sample solution which is the actual parameter which needs to be determined. After a certain time an equilibrium state is reached in which the rate at which the proteins 2 are bound to the sensitive region 14 is equal to the rate at which the proteins are released again. This time-dependent adsorption mechanism is called “Langmuir adsorption”, and
To describe the time-dependent occupation of the sensor surface, several parameters are of importance such as the concentration of the proteins 2 in the solution (target concentration [T], measured e.g. in mol per unit volume), the number of possible adsorption sites 3 on the surface (antibody concentration [Ab], measured e.g. in sites per unit area), a parameter which describes the chance of attachment of a protein 2 to an antibody 3 (the “forward” reaction constant kon), and a parameter which describes the release of the protein 2 from the antibody 3 (the “reverse” reaction constant koff). Formula (2) of
In an end-point measurement which is typically used by a number of known techniques, one would inject the sample solution into the sample chamber 1, wait for some time tm and do a readout of the sensor signal. From the signal, the number of proteins 2 on the surface can be determined and thus the concentration in the solution. However, from the signal no estimation of the error in the signal can be obtained other than the theoretically expected error. In the following it will be shown that the described magnetic biosensor device allows a dynamic measurement of the slope of the binding curve which a) will allow a more accurate determination of the slope than a single end-point measurement, and b) will give an estimation of the error in the slope as well.
The slope of the Langmuir isotherm at t=0 is linearly dependent on the target concentration [T]. By determining this slope, the concentration can therefore be calculated. If the Langmuir isotherm consists of a discrete number of n measured points, the slope can be determined by using a linear (or a weighted linear) regression. Assuming a total measurement time tm, each individual measurement lasts a sampling time Δt=tm/n and the sampling rate at which the signal is sampled is equal to 1/Δt. The sensor signal is linearly dependent on the number of bound proteins, so the signal Si from an individual measurement i in time can be written as in equation (5) with a proportionality constant a′.
The slope of the signal versus i·Δt is equal to a′·[T]. As previously described, the noise in the signal consists of two different kinds of noise: a) the thermal noise in the sensor units and electronics, which is independent of the number of particles and averages out better for longer sampling times Δt, and b) statistical noise. The latter noise signal scales with ✓[T]. The variances in the individual data points are described by equation (6).
By using linear regression on the n data points, it can be shown that the signal-to-noise ratio, SNR, of the slope a′·[T] can be written as in equation (7). For a relatively large number of data points (i.e. n→∞), this SNR reduces to equation (8).
Given a maximum allowed measurement time tm, the SNR of the biosensor has to be optimized with respect to the number of data points n (and thus the sampling rate) and the number of sensor units N. However, since the target concentration [T] is still present in the expression (8), the sensor can only be optimized for one specific concentration, which is disadvantageous.
To overcome this limitation, it is proposed to adapt the number n of data points (and thus the sampling rate) to the concentration [T]. More specific, the sampling rate, n/tm, is chosen equal or faster than the adsorption rate of the proteins according to equation (9). In words this means that the sampling rate should be fast enough to catch all adsorption events since every adsorption event carries information. Taking a sampling rate (orders of magnitude) slower than the adsorption rate misses information, sampling faster does not add extra information but also does not harm the SN-ratio. By substituting n from equation (9) in equation (8), the optimal value Nopt for which the signal-to-noise ratio is maximal with respect to N becomes independent of the target concentration [T], cf. equation (10).
Since the adsorption rate rads is unknown at the beginning of the measurement, it is further proposed to split the measurement into two or more parts:
The advantage of using a regression technique with optimized sampling rate over an end-point method where only one data point is taken at t=tm is four-fold:
The table of
In summary, the central aspects of the proposed method are:
Instead of performing an end-point measurement with the biosensor in which first the target molecules of interest are collected on the sensor surface, followed by the actual measurement of the target concentration, it is proposed here to dynamically measure the collection process with the advantage that the concentration measurement can be done much more accurately while also an estimation of the statistical error can be obtained.
In the following, further embodiments of the present invention will be described that are based on the detection of events related to the movement of single target particles (or at least a small number of target particles).
It should be noted that the sensor device may optionally comprise any combination of the features described with respect to the previous Figures (and vice versa). Moreover, it should be noted that the detection principle which will be described in the following with respect to the magnetic sensor unit 110 are also applicable to other types of sensors, for example optical sensors that use the principle of frustrated total internal reflection of an incident light beam at the bottom surface 4.
Conventionally, biosensors are operated in the linear regime, i.e. the sensor response is proportional to the density of target particles 2 (e.g. super-paramagnetic beads linked to target molecules in the sensitive region). In order to relate the sensor response to the exact concentration of target particles in the sample volume, the sensor sensitivity has to be calibrated. During measurements the sensor sensitivity or the properties of the read-out apparatus may slightly change and an additional control system is required to check and correct these variations.
To address these problems, a non-linear read-out method for microelectronic sensor devices is proposed that is based on the movement of target particles explained above. This method distinguishes from conventional linear read-out by detecting signal events, i.e. short time occurrences or persistent signal changes resulting from movements of target particles in the sensitive region.
By detecting and counting events in the sensor signal that correspond to the entrance of target particles into the sensitive region, particularly to their binding, the number of immobilized target particles on the sensor surface can be determined without (re-) calibration. The method further enables discrimination of signal events corresponding to single target particle binding or to the binding of clustered particles, thereby making the detection method robust to clustering.
Due to the fact that target particles move into and out of the sensitive region by thermal motion, the number of free target particles above the sensor can be determined by detecting and counting events in the sensor signal that correspond to target particles entering and leaving the sensitivity volume.
In the following, a number of signal realizations corresponding to particular events in the sensitive region is analyzed, but the proposed method is not limited to these particular events or analysis. In addition, the proposed signal analysis techniques can be operated in place of or complementary to linear detection methods.
By detecting and counting events in the sensor signal S that correspond to label binding, the number of immobilized target particle labels on the sensor surface can be determined. To that end, the rate at which the sensor response is sampled must be sufficiently high so that individual binding events can be distinguished.
Curve “S_a” of
The exact amplitudes Δ of the binding event signals are of secondary importance, thus making this non-linear method independent of the sensor sensitivity, its calibration or variations therein, as well as non-uniformity of the particle labels.
Curve “S_aa” of
In practice, the sensor signal S will be perturbed with noise. Based on prior knowledge of the signal shapes corresponding to binding events, filters can be constructed to match these signals (cf. e.g. L. A. Wainstein and V. D. Zubakov, Extraction of signals from noise, Prentice-Hall, Englewood Cliffs, UK, 1962). Matched filters can be applied in a signal post-processing system for the purpose of increasing the signal-to-noise ratio and thus the ability to detect binding events. The present invention encloses the application of matched filters to binding event detection, but is not limited to this technique. Other methods for the detection of binding events in the sensor signal are also included.
As was already explained, the target particles 2 may attach to each other forming more or less large clusters 2c. A sensor signal that corresponds to the entrance of such a cluster 2c (with N=2 particles) into the sensitive region 114 and its binding to the sensor surface 4 is shown in curve “S_c” of
A more detailed analysis of the aforementioned situation starts with the observation that the target particle velocity determines the rise time of the signal step, being defined as the time the signal requires to increase from its initial value to its persistent value. In the sensitive region 114, the target particle velocity v is dominantly governed by the magnetic force exerted by the excitation wires 111, 113. As can be seen from the formula of
Loosely speaking, a cluster of N target particles can be regarded as a single target particle having an N times larger volume, or equivalently having a N1/3 larger diameter. The velocity of said cluster thus scales with N2/3, and consequently the rise time of the signal increases with this factor, as illustrated by curve S_c of
The sensor response is proportional to the magnetic moment of a bead, and thus to the target particle susceptibility and volume. As a result the persistent signal from a cluster bound to the sensor surface is substantially larger than that of a single bead. In a first approximation, the amplitude of the signal step that is induced by a cluster of N particles is N-times larger than the step induced by a single particle.
By simultaneous analysis of the rise time and amplitude of the step signal, cluster binding events can be discriminated from single target particle binding events. Moreover, consideration of the rise time helps to discriminate the binding of N-particle clusters from occasionally occurring simultaneous bindings of N single particles, as is illustrated by curves S_aa and S_c in
Based on prior knowledge of the signal shapes corresponding to cluster binding events, filters can be constructed to match these signals. The present invention encloses the application of matched filter banks to both single binding event detection and cluster detection, but is not limited to this technique.
According to another aspect of the described approach, the number of free target particles 2 above the sensor can be determined by detecting and counting pulses in the sensor signal S that correspond to target particles 2 entering and leaving the sensitive region 114. Due to thermal motion, target particles 2 constantly move into and out of the sensitive region 114. The number of particles in the sensitive region is characterized as a spatial Poisson process, with mean and variance equal to the average number of particles in the volume. The sensor response to a target particle 2 migrating into and out of the sensitive region 114 will result in a signal pulse. Clearly, such a pulse does not have a persistent value, since the target particle will leave the sensor sensitivity zone, and can thus be distinguished from binding events. By counting the number of pulse events during the diffusion time, an estimate of the number of target particles in the volume can be obtained.
The average number of free target particles 2 in the sensitive region 114 is linearly related to the total number of target particles in the sample volume. In particular if an inhibition assay is used to detect small molecules, the knowledge of the number of target particles in the sample volume is essential.
It was already discussed that the rise time of a signal event is proportional to the target particle velocity. By examining the rise time distributions of the various signal classes enclosed in the previously described embodiments, the average velocity of target particles 2 can be determined. If the average properties of the target particles (or their labels) such as susceptibility and volume are known, then the average magnetic force acting on the target particles can be determined. From this information and the average velocity measurements, the fluid viscosity η may be obtained according to the formula of
The main advantages of the embodiments shown in
Finally it is pointed out that in the present application the term “comprising” does not exclude other elements or steps, that “a” or “an” does not exclude a plurality, and that a single processor or other unit may fulfill the functions of several means. The invention resides in each and every novel characteristic feature and each and every combination of characteristic features. Moreover, reference signs in the claims shall not be construed as limiting their scope.
Number | Date | Country | Kind |
---|---|---|---|
06113707.1 | May 2006 | EP | regional |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB07/51505 | 4/24/2007 | WO | 00 | 11/5/2008 |