The goal of data fitting is to draw a line through a data set. Before fitting the data, we need a model. A model is a formal mathematical representation of a chemical or physiological idea. To be useful, the model must be expressed as an equation that defines Y, the outcome of a measurement, as a function of X and one or more variables to be fitted. The term variable is used to refer to the terms in the equation to be fitted. With non-linear regression, the term variable does not refer to X and Y (1).
For instance, the law of mass action is the model for the simplest binding reaction between two proteins (1),(2). The reaction between immobilized ligand (L) and an analyte (A) can be assumed to follow a pseudo first order kinetics (3),(4),(5). The reaction equation looks like this:
During the association phase the complex [LA] increases as a function of time. The differential rate equation describes the relation between the parameters:
In the equation for complex formation the k_{a} is the association rate constant which describes the rate of complex formation, i.e. the number of LA complexes formed per second in a one molar solution of L and A. The units of k_{a} are M^{-1}s^{-1} and are typically between 1.10^{3} and 1.10^{7} in biological systems.
k_{d} is the dissociation rate constant which describes the stability of the complex, i.e. the fraction of complexes that decays per second. The unit of k_{d} is s^{-1} and is typically between 1.10^{-1} and 1.10^{-6} in biological systems. A k_{d} of 1.10^{-2} s^{-1} = 0.01 s^{-1}. This means that 1 percent of the complexes decay per second.
Although the differential rate equation describes the one to one Langmuir reaction accurately, it is not useful as a model to fit to a curve. By integrating the differential rate equation, sensorgrams can be analysed directly using non-linear methods (6). This is because integrated rate equations describe the whole curve as opposed to the differential rate equations, which describe the slope of a curve.
The classical Langmuir surface adsorption (7) can be described with the following integrated rate equation(6)
While simple rate equations are easily integrated, the more complex (e.g. with mass transport parameters) are not. The more complex rate equations are resolved by numeric calculations. Using numerical integration (a.k.a. non-linear regression) has also the advantage that global analysis can be linked in. This means that several curves are fitted simultaneously giving a single and more robust answer for the rate constants (8).
Non-linear regression is a more general approach to data fitting because all models that define Y as a function of X can be fitted. The process will find the values for the variables giving the closest fit, reducing the sum-of-squares (SS) to a minimum. The SS is calculated from the vertical distances between the curve and the measured points.
Non-linear regression can fit complex models by an iterative process, which starts with initial values provided by the user. The iterative process goes on until the convergence criteria are matched. For instance, when two iterations change the SS by less than 0.01%.
The initial values are estimates of the real values. During the iterative process, the program adjusts these values to improve the fit. Initial values matter more when the data has a lot of scatter or does not span the whole curve. Depending on the initial values, different results may follow. This phenomenon is the result of the search for a localminimum.
In the figure 'Local minimum', there is a sinus curve. Searching for the true minimum can start on either side of the curve. Only one side will give the true minimum. So, when a fitting does not work, try different initial values in steps of 100 times difference.
Often the Chi^{2}-value is given after the fitting is done. The Chi^{2}-value shows information about the goodness of fit. For a good fitting, the Chi^{2} is in the same magnitude as the noise of the used system (but check the validation section).
The results of non-linear regression are meaningful only if the following assumptions are true (or nearly true):
(1) | Graphpad Software Graphpad Prism. (1998). Goto reference |
(2) | Elwing, H. Protein absorption and ellipsometry in biomaterial research. Biomaterials 19: 397-406; (1998). |
(3) | Bjorquist, P. and Bostrom, S. Determination of the kinetic constants of tissue factor/factor VII/factor VIIA and antithrombin/heparin using surface plasmon resonance. Thrombosis Research 85: 225-236; (1997). |
(4) | O'Shannessy, D. J. et al Determination of rate and equilibrium binding constants for macromolecular interactions using surface plasmon resonance: use of nonlinear least squares analysis methods. Analytical Biochemistry 212: 457-468; (1993). |
(5) | Johnsson, B. et al Immobilization of proteins to a carboxymethyldextran-modified gold surface for biospecific interaction analysis in surface plasmon resonance sensors. Analytical Biochemistry 198: 268-277; (1991). |
(6) | BIACORE AB BIACORE Technology Handbook. (1998). |
(7) | BIACORE AB Bia Journal Article. Bia Journal 2: 18(1998). |
(8) | Karlsson, R. and Falt, A. Experimental design for kinetic analysis of protein-protein interactions with surface plasmon resonance biosensors. Journal of Immunological Methods 200: 121-133; (1997). Goto reference |