The present invention is related to applications and methods which utilize surface plasmonic (SP) absorption at or near its resonance for detecting the presence and absence of molecules, biological cells or organelles on or near a surface, or structural and electronic changes of the molecules through their interactions with an analyte, particularly to a method of characterizing interactions between the analyte in the fluid and the ligand immobilized onto multiple sensing surface areas.
Many analytical techniques are used to characterize interactions between biomolecules. For example, the protein-drug interaction is an important type that is of fundamental importance to molecular biology, immunology, drug discovery, and pharmacokinetics and pharmacology. In this context, the analytical techniques focus on immobilization of a ligand, such as an antibody or other biomolecules, to a sensing surface, followed by exposing the ligand with an analyte, such as an antigen, via a fluidic delivering system. With exposure of the ligand to analyte, some characteristic change of the sensing surface is measured, and the change is indicative of the interaction, such as the ability of the ligand to bind the analyte and the rate of binding.
The Surface Plasmon Resonance (SPR) sensor system is a particularly popular method to monitor biomolecular interactions in real time without the necessity of attaching a label to the analyte. Many commercial biosensor systems based on SPR have become powerful tools for such interaction analysis. A more detailed description of SPR and theory of detection may be referred to [Reference 1]. All SPR systems use a flow system for sample delivery and optoelectronic elements to detect the analyte-ligand interaction at the sensing surface. By measuring the changes of the SPR angle, a plot (sensorgram) of the response signal vs. time is generated and analyzed to obtain relevant information on binding interactions, such as affinity and association and dissociation rates of the ligand-analyte binding complex.
A representative sensorgram for a reversible interaction at the sensing surface is presented in
After the sample injection is completed, the sensing surface is immediately flushed with an analyte-free buffer to allow dissociation of the analyte from the binding complex to take place, this is referred to as the dissociation phase, which causes the response signal 5 to decay as a function of time (in the form of e−kd t wherein kd is the dissociation constant). Furthermore, the sensing surface can be “washed” as shown in step 6 with an agent to “regenerate” the original surface without any analyte or even the pre-immobilized ligand, thereby allowing the sensor surface 1 to be reused for subsequent measurements. In practice, most sensor chips can only be “reused” for a limited time (generally less than 10 cycles) and in many cases the sensing surface cannot be regenerated at all. This aspect of SPR measurements significantly affects both the sample throughput and increases the assay costs and has not been sufficiently addressed in the past.
The sensorgram provides the essential information from which the binding kinetic data could be readily derived using pseudo-first order binding kinetics theory. This protocol has been widely used in all SPR biosensing measurement for the last 20 years. However, despite the development of highly advanced fluidic systems and sophisticated modification of sensor surface, as well as numerous publications on various applications of SPR to a wide range of biomolecules, optimization of SPR measurement parameters remain difficult and requires good knowledge of surface chemistry, mass transfer limitation, kinetic theory, and etc. [Reference 2-3]. Mistakes are often made with inaccurate kinetic data, due to the lack of understanding on how these secondary effects (mentioned above) impact the observed binding behaviors. To improve the quality and accuracy of the binding kinetic measurement, a wide range of concentrations of the analyte with numerous repeats and carefully designed references must be used to identify and eliminate the secondary effects in the measurement. [Reference 4]
In a typical SPR experiment on the “1:1” reaction, the ligand B is immobilized onto a sensor surface and its interaction with an analyte A injected onto the sensing surface through a microfluidic channel can be expressed as:
During the association phase (rate constant ka), the SPR signal (Ron) is given by
where C is the analyte concentration, ka the association rate constant, Rmax the maximum analyte coverage at the sensor, and t the reaction time.
To initiate the dissociation reaction, the injected analyte is rapidly replenished with a running buffer. The SPR signal (Roff) decays exponentially with time from R0, the signal at the moment when buffer replaces the analyte in the fluidic channel:
R
off
=R
0
e
−k
t (3)
If sufficient time is allowed during the sample injection for eq. 1 to reach the equilibrium, the SPR signal reaches a steady-state value Req, as shown by eq. 4:
To derive kinetic quantities of ka and kd, and thermodynamic parameters KD=kd/ka and Rmax, either several Rmax or C values in eq. 2 are needed for reliable fitting of the association phase as eq. 2 has multiple variables. Traditionally, a series of analyte concentrations have to be injected to generate many curves of Ron(t), which are subsequently fitted with the kinetic theory. Notice in eq.4, Req vs C should follow the Langmuir isotherm behavior and KD=kd/ka can be derived without knowing the ka and kd values. For example, as shown in
So far, there is no attempt on varying Rmax at a fixed C value in eq. 2 to calculate kinetic quantities. To vary Rmax and fix C, it will require an SPR system which has at minimum of 3 channels to immobilize the ligand with different and known surface densities.
It is important to note that eq. 1 assumes that the reaction has a pseudo-first order binding kinetic behavior, ignoring all the secondary effect such as analyte aggregation, nonspecific absorption, bulk refractive index change, mass transport limitation and etc. When secondary effects are considered, the kinetic behavior of eq.2 and eq.4 become
respectively. In eqs. 5 and 6. kf and kr represent the forward association and reverse dissociation rate constants. They are functions of ka, kb, the amount of ligand immobilized [B], and the rate of analyte mass transport to the surface (km). In addition, bulk refractive index change will lead to large baseline variations, and secondary effects will cause complicated binding behavior to change. Therefore, the pseudo-first order binding kinetic theory may become invalid. In this case, Req in eq.6 at equilibrium state is no longer directly proportional with Rmax. It has been extremely difficult to calculate kinetic quantities of ka and kd accurately in the presence of the secondary effects.
There is only one exception in the secondary effects which the pseudo-first order binding kinetic theory can still be used with minor corrections, that is when the binding event is partially controlled or limited by the mass transport (MT). Assuming MT is the only secondary effect present, kf and kr can be represented as:
By substituting eq.7 and eq.8 into eq. 6, Req becomes the same as eq.4. This implies that in the presence of the MT effect, the pseudo-first order binding kinetic theory with modified kf and kr can still describe the binding kinetic behavior [Reference 5].
Various practices for minimizing the secondary effects have been exercised over the past 15 years, which include: using the ligand density as low as possible to avoid complete or partial MT limitation; and for high-affinity binding involving small analyte molecules, using a high immobilized ligand density as the steric hindrance is small for small molecules that diffuse much faster than bulky biomolecules. To optimize the measurement condition and to minimize the secondary effects, injections of various analyte concentrations with regeneration cycles in between have to be conducted to validate the measurement and fitted with the pseudo-first order kinetics. Often the Langmuir isotherm curve is used to fit with Req vs C in order to validate the data quality. This brutal force approach is necessary because there is not a simple way to validate the measured data suitable for accurate pseudo-first order binding kinetic fit when different analyte concentration C are used. [Reference 6-7]
This widely accepted approach (i.e. fixing the ligand density and varying the analyte concentrations) is tedious and costly in terms of sample consumption and resources (analysis time, labor, reagents, sensor chips, etc.). A major problem inherent in the conventional approach is that the sensing surfaces must be “refreshed” or “regenerated”, which also degrades sensing surface and leads to error in measurements, and in some cases, is not possible.
Another problem is that the concentration of the analyte directly affects the time scale for binding kinetic to reach equilibrium as shown in equation 2. When a lower analyte concentration is used, the time takes to reach equilibrium state increases as the exponential coefficient in eq. 2 becomes smaller. This prolongs the measurement time significantly.
It is the principal objective of the present invention to overcome the above mentioned problems in the prior arts by providing a novel method to validate the binding data suitable for pseudo-first order binding kinetic fit in measuring the said interaction of an analyte with ligands, therefore, kinetic parameters can be quickly and accurately determined without making injections of multiple solutions and surface regeneration.
In the present invention, by focusing on varying the immobilization amount of a ligand (Rmax) instead of the analyte concentration (C) in eq. 2, a novel method is provided for validating the kinetic data in measuring the interaction of an analyte in the flowing stream with a ligand immobilized on the sensing surfaces. First, using a multi-channel SPR instrument to prepare multiple sensing areas through an in-line process to controllably immobilize graduated amounts of a ligand onto each area or zone. Then for binding kinetic measurements, an injection of analyte of a given concentration is made to flow over the said areas to obtain its binding kinetics. To validate that the measurement obeys the pseudo-first order binding kinetics, a plot of equilibrium binding responses (Req) vs. the amounts of immobilized ligand (Rimmob, which is propotional to Rmax) is generated. Only when Req is directly proportional to Rimmob will the pseudo-first order condition be valid. If a portion of the plot is not directly proportional to Rimmob due to the secondary effects, it will be discarded. In other words, only the data in the region wherein a proportional relationship exists will be used for calculating the binding kinetics.
As a result, the new approach removes many uncertainties without injection of multiple solutions and surface regeneration. Consequently, the method saves time, samples, and sensor chips, and allows equilibria to be established at the same time, and avoids any damages/changes that may occur at the sensor surface during its regeneration. It is worth noting that even though in theory only three channels are needed (two channels necessary to have two signals for solving two unknown ka and Rmax, and one channel for background signal recording), in practice five or more channels will be essential to exclude the interference caused by the secondary effects). This is because any two data points (two signals from two channels after subtraction of the background from the third channel) can render a liner plot. Only with at least 3 data points can a linear regression be meaningful. With four or more data points (i.e. five or six fluidic channels), the fidelity of the linear regression is further enhanced.
In one aspect, the present method is carried out by controlling the exposure (immobilization) time of the ligand to the sensing surface in each channel, so various ligand amounts can be immobilized. Afterward, loosely bound ligand molecules are washed off and the active or unreacted sites remaining on the sensor chip are blocked. The responses of the binding kinetics of ligand with analyte are measured by a single injection of the analyte into the flowing stream over the multiple sensing surface areas.
In another aspect, the present method keeps a fixed concentration and varies immobilization amount Rimmob. according to equation 2, all exposed sensing areas or channels reach equilibria Req at the same time because (kaC+kd )t is the same for all sensing areas. This allows quick affinity measurement to be carried out by using a simple linear regression fit of the data to
or eq 4 to obtain KD.
In still another aspect, the present invention provides an analytical method for studying molecular interactions, which comprises a computer processing procedure including program code and algorithm for performing the steps of said method.
These and other aspects of the invention will be evident upon reference to the accompanying drawings and the following detailed description.
For simplicity and clarity of illustration, the drawing figures illustrate the general manner of construction, and descriptions and details of well-known features and techniques may be omitted to avoid unnecessarily obscuring and invention. Additionally, elements in the drawing figures are not necessarily drawn to scale. For example, the dimensions of some of the elements in the figures may be exaggerated relative to other elements to help improve understanding of embodiments of the present invention. The same reference numerals in different figures denote the same elements.
As mentioned above, the present invention relates to a method for quickly validating the binding kinetic measurement for first-order binding kinetic region in assays or studies involving the detection of binding events of analyte at a sensing surface covered with immobilized ligand.
All measurements for binding kinetics in prior arts have been performed by measuring many binding curves with multiple injections of different analyte concentrations on the sensing surface with a fixed ligand density. With the availability of three or more channels in SPR systems such as the BI-4500 five-channel SPR system, it becomes possible to controllably vary immobilized ligand amounts instead of analyte concentration to perform kinetic measurements, thus providing many distinct advantages over prior art to validate the binding kinetic data based on the pseudo-first order binding kinetic theory.
In the following description of present art, various aspects of the present invention are disclosed more specifically for the purposes of illustration and not limitation.
In the present invention, with a multi-channel SPR system and a single injection of only one ligand concentration, multiple sensing surface areas with different amounts (densities) of a ligand can be immobilized by simply exposing the ligand solution to each sensing area for different periods of time. As shown in
To perform binding kinetic measurements, an analyte of a given concentration is then injected into the flowing stream over the immobilized sensing areas, again in a single injection. The responses of binding kinetics are measured simultaneously from the sensing areas.
R
eq
∝R
immob
In
In a variant of the present invention, eq.4 can also be used to find the right range of immobilized ligand densities that can be used to fit with pseudo-first order binding kinetic theory. As shown in
In another variant of the present invention, it is evident that accurate binding kinetic data can be obtained and validated with a single injection of analyte into multiple channels with various immobilized ligand densities. This makes it possible to measure the interaction with certain type of ligand molecules or antigen that could not be regenerated.
An example of binding kinetic measurement that is severely affected by the regeneration step is the interaction between a monoclonal antibody to the green fluorescent protein (GFP). GFP antibody is immobilized onto a sensor surface for binding kinetic measurement, but it cannot be regenerated because GFP would be greatly denatured by regeneration chemicals such as HCI or NaOH. Thus kinetic measurements of the interaction with injections of various analyte concentrations with regenerations in between could not be used in this case.
Such a problem is mitigated with the present invention. With a graduated immobilization of GFP onto four sensing areas (channels) 31 to 34 as shown in
In still another variant of the present invention, using a fixed concentration and varying immobilization amounts provides a quick and simple way to obtain affinity KD by using a simple linear regression fit with the said data set to equation 4 or Req=C/C+KDRimmob. According to eq. 2, at fixed C, all exposed sensing areas (channels) reach quasi equilibrium Req at the same time as e−(kC+k
Although the invention has been described with reference to specific embodiments, it will be understood by those skilled in the art that various changes may be made without departing from the spirit or scope of the invention, accordingly, the disclosure of embodiments of the invention is intended to be illustrative of the scope of the invention and is not intended to be limiting. It is intended that the scope of the invention shall be limited only to the extent required by the appended claims.
This application claims the benefit of United States Provisional Application No. 62/465,954 filed on Mar. 2, 2017
Number | Date | Country | |
---|---|---|---|
62465954 | Mar 2017 | US |