Patent Application
20080033705
An SPR sensor system AP1 in accordance with the present invention provides for determining when a single run suffices for determining kinetic parameters such as association rate (ka) and dissociation rate (kd), as well as parameters calculable from kinetic measurements, e.g., affinity (KD). When a positive determination is made, the method increases analytical throughput by requiring only one run per binding pair as opposed to the more typical five. In the context of drug discovery, even when the determination is negative, throughput need not be adversely affected, as the negative determination typically indicates a poor drug candidate. Other embodiments of the invention provide other advantages instead of or in addition to higher throughput.
SPR sensor system AP1 comprises a flow channel 11, a laser 13, a prism 15, a substrate 17, a photodetector 19, a computer 21, and a display 23. Substrate 17 is glass (or other dielectric) with a coating 25, including gold (or other noble metal). SPR sensor system AP1 is designed to track binding and dissociation between a fixed interactant 31 and a mobile (while unbound) interactant 33, for example, between a fixed target and a mobile drug candidate “analyte”. The fixed interactant is immobilized on coating 25 in that it is attached to the gold or to an intermediate layer on the gold. The mobile interactant is in solution and flowed over coating 25. When bound, the fixed and mobile interactants form bound pairs 35.
In SPR sensor system AP1, most of the light incident substrate 17 is reflected. However, over a range of angles of incidence, some of the incident light is converted to an evanescent wave in the vicinity of the gold coating, resulting in a reduction in the intensity of the reflected light. The angle at which this reduction is maximal is called the “resonance” angle. The resonance angle is sensitive to the refractive index of material near the gold coating 25 opposite substrate 17. The refractive index changes with the population of bound pairs 35. Accordingly, the binding activity, as represented in the population of bound pairs can be tracked by monitoring the shift in resonance angle over time.
In
Computer 21 provides the hardware and software to implement a method ME1, shown flow-charted in
1) an association phase 47 represents binding activity from a time T1 when the fixed interactant is unpopulated to a time T2 when equilibrium (i.e., when the rate of dissociation equals or nearly equals the rate of association) is reached or nearly reached.
2) a dissociation phase 49 during which interactant-free buffer is introduced to the fixed interactant and any bound mobile interactant from time T2 to a time T3 when the amount of binding pairs has reached or nearly reached a minimum (absent some purging activity).
At method segment M2, computer 21 determines a best-fit time-profile curve for the binding interactions. The curve is selected according to a model of the interaction being tracked. In this case, a reversible interaction obeying the laws of mass transport is involved, so exponential curves are fitted: 1) an “association” exponential curve 51 for association phase 47; and 2) a “dissociation” exponential curve 53 for dissociation phase 49. Of course, the best-fit exponential curves largely overlap the plots they are intended to fit, so the curves and plots are most readily distinguished for times when the fit is less optimal.
In general, the type of curve that is fitted is dependent on an interaction model. Models for cooperative binding (which can involve interaction of analytes to form complexes that bind to a target), and for coupled kinetic reactions (in which multiple equilibria are involved in an interaction) yield their own respective curve types.
At method segment M3, computer 21 evaluates the curvature (or non-linearity in time) of the association and dissociation time-profile curves 51 and 53. The curvature of a time-profile curve is taken relative to a linear rate of change. For example, the linear rate of change can be represented by a chord (i.e., a straight-line segment that intersects a curve at both its endpoints) that intersects association curve 47 at times T1 and T2. For every response level between these times, association curve 47 reaches this level before it is reached by the chord. The curvature can be measured as a maximum lead time, an average lead time, or weighted or normalized variations of maximum and average lead time, or other related parameters. In the case of an exponential curve R=e−kx, the exponential constant k can serve as a measure of curvature. In the illustrated case, the exponential constants for time-profile curves are used to evaluate time profile curves 51 and 53 independently. In addition, the product of the exponential constants serves as a combined evaluation parameter for the pair of curves 51 and 53.
At step M4, low curvature profiles are flagged; in other words, computer 21 provides some indication that a curve cannot be relied on to evaluate a parameter of interest with the required precision. For example, in
The evaluations can be separately considered for each parameter of interest. For example, the association phase curvature can be used to determine whether ka can be reliably determined, the dissociation phase curvature can be used to determine whether kd can be reliably determined, and some mathematical combination of the curvatures for the two phases can be used to determine whether KD can be reliably determined.
The implications of a low-curvature flag can be that the curve and associated data can vary: 1) the associated data can be discarded and replaced by data from another run, e.g., using a higher concentration of mobile interactant; 2) the data can be supplemented using data from different runs of the same or different concentrations, but the same interactants; or 3) the associated data can be used to evaluate parameters of interest with low precision.
The threshold between insufficient and sufficient curvature can be determined for each parameter of interest empirically for a given sensor system and target during calibration. For example, in
Table I indicates that the calculated values ka and KD are relatively close below the double line, but far off the mark for data above the double line. Since the dissociation values kd are reasonably uniform, the fault appears to lie with the exponential constants for the fitted exponential constants of association. Values of ko above 0.13 appear to be reliable, while those below may not me. Hence 0.13 can serve as a cut-off curvature value for observed exponential curves of association. Note that further runs at intermediate concentrations (between 2.5 μM and 0.625 μM) might yield lower thresholds, while runs with different interactants can be used to confirm the general applicability of the threshold.
When the curvature is sufficient, as determined at method segment M4, the data from a single run suffices for calculating the parameter or parameters of interest and no further runs are required. This can result in a substantial improvement in throughput versus the typical five runs used to calculate kinetics and affinity values. Coupled with advances in parallelism for SPR sensors that permit a hundred or more runs to be performed using a single SPR machine, the invention provides for a several-hundred fold improvement in throughput relative to non-parallel, multi-run approaches.
This improvement has significant implications for target-based drug discovery and it becomes practical to perform kinetics measurements on more drug candidates. This in turn permits the kinetics measurements to be obtained earlier in the drug discovery process. Since kinetic measurements “include” endpoint measurements, the need for some endpoint tests is obviated. This results in a further increase in throughput for the drug discovery process.
For example, a flow chart of a target-based drug-discovery method ME2 is shown in
However, whereas, screening typically involves secondary screening using endpoint measurements, the present invention provides for kinetic measurements early in secondary screening. Thus, in screening M23, a primary screen is performed at method subsegment M31 using low-precision endpoint measurements on a large (˜106 compounds) drug library 60, yielding a much reduced number (e.g., ˜103) of hits. At method segment M32 all or most hits are kinetically screened using system AP1 and method ME1, with the target as the fixed interactant and drug candidates as the mobile interactants for respective runs. The kinetics for each target-candidate pair are measured once. The hits yielding insufficient curvature are rejected in the screening process, along with any hits that are reliable determined to have undesirable kinetic and affinity binding values. The resulting leads (e.g., about ˜102) are then subjected to counter screens at method subsegment M33, either using endpoint or kinetic measurements.
By screening hits immediately for kinetics, those with unfavorable on-rates are excluded from subsequent testing. Complementarily, some hits with moderate affinity but favorable on-rates (that might have been rejected using endpoint testing) can be carried forward. Thus, the process of selecting candidates becomes more robust. Any decrease in throughput due to kinetic versus endpoint testing of the hits is offset by the reduced number of screen stages required and the reduced number of poor candidates that are subjected to repeated testing.
“Monitoring” herein refers to generating an effectively continuous signal or signals representing a succession of states for the subject being monitored. “Effectively continuous” means continuous, as in an analog signal, or involving a sufficient sample rate to adequately capture all state changes of interest for the subject. (For example, in video recorders, a 30 frame-per-second sample rate is effectively continuous.) “Tracking” is “monitoring-plus-recording” the states determined by said monitoring.
“Binding reaction” is a reaction in which plural, typically two, bio-chemical entities become physically coupled. In the context of drug discovery, the binding partners or interactants are typically a drug target and a drug candidate. From the perspective of an SPR sensor, one of the interactants is fixed to a sensor surface, while the other is typically in a solution that is introduced to the fixed partner. In this context, the binding reactions are typically non-covalent and reversible. However, the invention has applicability to irreversible reactions and reactions involving covalent bonds as well.
“Kinetics” refers to reaction rates and how reaction rates change over time. “Kinetic measurement” implies that non-equilibrium conditions are monitored, although equilibrium conditions can also be monitored. Thus, kinetics involves characterizing a profile of binding as a function of time, e.g., as interactants increasing associate or increasing dissociate until equilibrium reached.
A “time-profile” curve is a continuous n-dimensional shape in n+m dimensions, each point of which is a well-defined function of time. In practice, a time-profile curve can be a continuous one-dimensional shape in two Cartesian dimensions, one of which represents time and the other of which represents a mass or other parameter corresponding to the extent of binding between two interactants.
“Fitting” involves selecting from a set of curves the curve that best meets some criterion, such resulting in the lowest sum of least squares of distances of data points from the fitted curve. The set from which the time-profile curve is selected corresponds to a theoretical model for the binding reaction. For example, most reversible binding reactions follow the law of mass transport; in this case, if the concentration of the mobile interactant remains effectively constant, the model calls for an exponential time-profile curve. In this case, the exponential that best fits the data is the time-profile curve for that data.
“Curvature” herein refers to deviations from linearity, e.g., of a time-profile curve. A straight-line (corresponding to a constant binding rate) can be a time-profile curve, but its curvature is zero. For other shapes, curvature values depend on the curvature-related parameter being measured. Where the time profiles are exponential functions of time, the exponential constant can serve as the curvature-related parameter. However, other parameters can include maximal or average time-deflection from a linear reference, area between a linear reference and the time profile, magnitudes of second-derivative peaks (where a rate changes most dramatically), as well as weighted combinations and variations of the foregoing. The linear reference can be a chord connecting two points on the time-profile curve, for example, the points corresponding to 10% and 90% levels in the span of the time-profile curve.
It should be noted that an observed exponential constant of dissociation determined from a fitted time-profile curve should be representative of the dissociation rate constant (kd) for the binding pair. On the other hand, an observed exponential constant of association (ko) determined from fitted time-profile curve is not generally representative of the association rate constant (ka). Instead, the latter is equal to the former less the dissociation rate constant (kd), with the difference divided by the concentration of the mobile interactant in solution in the vicinity of the fixed interactant, i.e., ka=(ko−kd)/[concentration].
“Validity” herein refers to the confidence level that can be assigned to a fitted time profile as a representative of the binding reaction being studied. When there is more curvature in the kinetic data, the range of viable exponential curve fits is smaller than when there is less curvature. Therefore, there is more range for error when the data is more linear and less basis for confidence in the selected exponential time-profile. “Evaluating the confidence level” here means assigning a value to indicate the confidence level for the fit. The result of the evaluation can be a one-bit value (pass versus fail), or a value on a continuum.
The value resulting from evaluating the confidence level can be a scalar or a vector. For example, one scalar value can be assigned to the association segment of a time profile and another scalar value can be assigned to the dissociation segment. In that case, a decision can be made based on the vector that includes both the association and dissociation scalar values. For example, for the purpose of calculating affinity, the decision to use a time-profile curve can be a function of a product of scalar values for association and dissociation.
