GB1601436A

GB1601436A - Quantitative protein analysis by immunodiffusion

Info

Publication number: GB1601436A
Application number: GB14576/78A
Authority: GB
Original assignee: Aladjem F J; Ayengar P K
Current assignee: Aladjem F J; Ayengar P K
Priority date: 1978-04-13
Filing date: 1978-04-13
Publication date: 1981-10-28

Description

(54) QUANTITATIVE PROTEIN ANALYSIS BY IMMUNODIFFUSION (71) We, FREDERICK J. ALADJEM & PADMASINE KRISHNA AYENGAR, United States citizens both of 845 Las Palmas Road, Pasadene, California 91102 U.S.A., do hereby declare the invention, for which we pray that a patent may be granted to us, and the method by which it is to be performed. to be particularly described in and by the following statement: The invention relates generally to the quantitative measurement of proteins in mixtures, especially when only small amounts of material are available.

In view of the current rapid expansion of knowledge concerning the role of proteins in health and disease, there is an increasing need for a general, rapid and relatively economical quantitative method for protein analysis of such fluids as serum, spinal fluid, tissue extracts and the like.

The proteins occurring in such fluids are frequently identified by immunochemical procedures which depend upon precipitation of each protein by an antibody specific to the particular protein. The production of such specific antibodies is stimulated when foreign proteins (antigens) are introduced into a living body. Antisera can be prepared, containing known mixtures of such antibodies. By reacting a protein sample in virtro with such an antiserum and observing the resulting precipitation, or lack of precipitation, useful information may be obtained as to the types of proteins in the sample.

A primary object of the present invention is to obtain quantitative values for the concentration of one or more proteins in an antigen solution, utilizing immunodiffusion procedures which have previously been regarded as only qualitative techniques.

According to the invention a quantitative measure of the concentration of a protein in an antigen sample is obtained after the sample has been subjected to immunodiffusion with an antibody source containing an antibody specific to said protein, the protein and antibody diffusing into reactive contact generally parallel to an axis of and through an area of a supporting medium initially free of both reactants to produce at least one elongated precipitation zone of limited length, the method of obtaining such quatitative measure comprising scanning the zone by optical means to develop electrical signals responsive to precipitin concentration at a two-dimensional array of positions distributed partly within the zone and partly outside the zone. deriving electronically from said electrical signals the corresponding value of a selected parameter of the zone which varies with said protein concentration, and comparing the resulting parameter value with a set of reference parameter values derived correspondingly from reference zones produced by equivalent immunodiffusion of a plurality of reference antigen solutions containing respective known concentrations of said protein.

The described collection of experimental data and the computations made with them can be carried out manually, if desired. They are also well adapted for semi-automatic data collection by optical and electronic scanning devices. Also, the necessary data processing can be made fully automatic with the aid of a general purpose computer of moderate capacity.

In preferred form of the invention, the proteins in the initial antigen mixture are first partially fractionated by causing them to migrate in one dimension at rates that differ characteristically among the various proteins. Such selective migration may, for example, utilize simple diffusion, electrophoresis, or more complex techniques such as chromatography without fundamentally altering the nature of the precipitin zones produced by the subsequent step of immunodiffusion. In the illustrative case of electrophoresis, differences of electrohoretic mobility between different proteins cause the proteins to become distributed in the direction of the electrical field in accordance with their mobility.

Following such initial fractionation, the resulting essentially linear distribution of proteins is brought into contact with the antiserum by relative movement in another dimension, typically by mutual diffusion in a suitable agar or agarose support medium. The precipitation zones of the respective proteins are then typically entirely separate and can be clearly distinguished.

Useful separation of the precipitation zones of a plurality of distinct proteins is also attainable without an initial step of fractionation, if the antigen and antibody are allowed to diffuse toward each other in a manner to form elongated precipitation zones of limited length extending transversely of the primary direction of diffusion. The mobilities for diffusion of different proteins, and/or of their antibodies, are ordinarily sufficiently different that such precipitation zones are clearly distinguishable. Any overlapping that may occur is usually limited to portions of the zones, the zone end points being usually clearly separated. The quantitative procedures of the present invention are usefully applicable to immunodiffusion of that type.

Similarly, if antigen and antibody are placed into wells for immunodiffusion, and an electric field is applied to accelerate the movement of the antigen and antibody to ward each other, as in the procedure known as electroimmunodiffusion, the resulting precipitin zones retain the same basic forms as in absence of an electric field. The same is true if the immunodiffusion step of immunoelectrophoresis is aided by a suitably directed electric field. Accordingly, the term "immunodiffusion" in the present specification and claims refers to diffusion with or without an accelerating electric field.

In the following description of certain illustrative members of carrying the invention into practice, reference will be made to the accompanying drawings in which: Figure 1 is a schematic plan of an immunoelectrophoresis plate, Figure 1A illustrating a typical distribution of proteins following electrophoresis, and Figures 1B, 1C and ID illustrating successive stages of the subsequent diffusion and immunoprecipitation reaction; Figure 2 is a schematic axial section representing typical apparatus for measuring a slide; Figure 3 is a schematic drawing illustrating properties of a precipitation zone relating to the invention; Figure 4 is a graph representing typical dependence of zone end points upon time; Figure 5 is a graph representing typical dependence of zone length upon time;; Figure 6 is a graph representing typical dependence of time of initial zone appearance upon protein concentration; Figures 7 and 8 are graphs representing typical dependence of zone length upon time and upon protein concentration; Figure 9 is a graph showing actual intensity scans across a precipitation zone at two incubation times; Figure 10 is a graph representing intensity scans across precipitation zones formed by respective protein concentraions; Figure 11 is a graph representing illustrative dependence of the parameter Iz upon time; Figure 12 is a graph representing illustrative dependence of the parameter Iz and of the parameter dIz/dt upon protein concentration;; Figure 13 is a schematic graph illustrating derivation of protein concentration values by known additions of the protein to the unknown sample; and Figure 14 is a schematic block diagram representing electronic scanning apparatus for carrying out the invention.

Many aspects of the invention are well illustrated by its embodiment in the preferred process of immunoelectrophoresis, and the following description will emphasize that form of the invention, but without implying any limitation of scope.

Immunoelectrophoresis is well known as a qualitative procedure, and many forms of apparatus have been decribed for carrying it out, differing in detail rather than in principle.

The electrophoresis and subsequent diffusion and immunoreaction are typically performed in a single layer of gel from a fraction of a millimeter to several millimeters thick carried on an optically transparent sheet. A currently preferred supporting medium is agarose saturated with barbital buffer of approximately pH 8.6 and ionic strength 0.1. The active materials are typically inserted into wells, which may be cut from the gel coating or formed when the gel is molded on the carrier.

Figure 1 shows schematically a plate 20 with a typical arrangement of wells, comprising the two circular antigen wells 21 and 22, spaced equally on opposite sides of the elongated antibody well or trough 24 with axis 25, which extends parallel to the direction of electrophoresis. With that geometric arrangement of wells two distinct antigen solutions, or two specimens of the same solution, can be run simultaneously against the same antibody solution on each plate. If preferred, the capacity of each plate can be multiplied, as by providing additional antibody wells outward of the two illustrated antigen wells, with additional antigen wells outward of them. Similarly, additional antigen wells may be added, spaced far enough from wells 21 and 22 in the direction of electrophoresis to prevent overlapping of the patterns.

The shaded areas 26 and 28 in Figure 1A represent typical approximate distributions of four varieties of protein a, b, c and d from duplicate samples in the respective wells 21 and 22 following a period of electrophoresis in the direction of the arrow 23. Although all proteins usually migrate in the same direction through the liquid medium, the solvent itself tends to carry a net charge and to have a resultant flow or electroosmosis relative to the gel.

Hence the proteins may have a net movement in either direction relative to the well.

Following termination of electrophoresis and addition of antibody solution to well 24, precipitation zones are produced by mutual diffusion of the proteins a, b, c and d of Figure 1A and the corresponding antibodies. Figures 1B, 1C and 1D show typical zones at respective stages of development.

The precipitin arcs of unrelated antigens, reacting with the antibodies against them, form independently, and may cross when they are sufficiently close together, as shown in Figure 1D; whereas the arcs of antigens that are immunochemically related typically join in a continuous reaction. The zone of precipitation d' in Figures 1C and 1D indicates that the area d of Figure 1A actually contained two distinct proteins, illustrating the fact that even antigens which have identical electrophoretic mobility may form distinct zones of precipitation. The zones d and d' may alternatively be viewed as representing two distinct proteins that were initially places in well 21 and were subjected to immunodiffusion without initial electrophoresis. The clear separation of the respective zone ends would then be due to different rates of diffusion of the proteins or of their antibodies.In either case, each such zone can be analyzed independently, using any or all of the methods of analysis to be described.

In accordance with the present invention, the zones of precipitation resulting from immunoelectrophoresis are subjected to direct quantitative measurements. Such measurements may determine only the physical position on the plate of certain selected features of each precipitation zone of interest, or may include quantitative photooptical measurements of the light intensity. For both types of measurement, the time is noted and may be used as an integral part of the observational data. If the incubation process is allowed to reach equilibrium, producing a nearly static zone configuration, the exact time of measurements becomes immaterial.

Position measurements on slide 20 can be made, for example, with the aid of a low power microscope. As represented schematically in Figure 2, plate 20 is placed over an adjustable opening 32 in the top of a light-tight box 30 with the lamps 36 and the dark backing 34 of light absorbing material such as black velvet. The microscope 40 has the objective lens 41, the eye piece 42 and a set of cross hairs or other reference reticle at 43 in the focal plane.

The microscope is mounted above light box 30 on a double slide mechanism 45 with screw drive 46, and with accurate scales, not explicitly shown, for reading the microscope position in two coordinates. For clarity, only one coordinate of the movement is directly indicated in the drawing. It is usually convenient to select coordinates having the x-axis, say, parallel to the direction of electrophoresis and to the length of the antibody well, and having the origin of coordinates at or near the well axis 25, as indicated in Figure 1A.

For measurements of the light intensity the microscope typically includes an oblique beam-splitting mirror 48 which sends part of the light to the eye piece while another part forms a real image in the plane of a diaphragm 52. Diaphragm 52 then transmits to the photosensitive transducer 50 only radiation from the elemental area of plate 20 that coincides optically with the cross hair image. Transducer 50 is electrically connected to the amplifying circuit 54 and the meter 56. The latter may be observed visually and the value recorded manually, or the meter may embody means such as analogue-to-digital circuitry and printout mechanism for automatically recording the light intensity in response to a command signal.Many changes can be made in the illustrative apparatus of Figure 2, including, for example, replacing the dark field illumination by direct lighting or by top lighting such that light reflected from the zone is observed or measured. Illustrative apparatus for making completely automated measurements is described below.

Figure 3 illustrates schematically certain preferred precipitation zone features which are selected by the present invention for position measurement as the zone develops. The horizontal line 61 at y = YAb represents the adjacent edge of the antibody well in the immunoelectrophoresis slide. The points E and F at the coordinates (xe, ye) and (xf,yf) represent the left and right end points of the elongated zone 60.

In addition to the zone end points E and F, we have found it useful to measure coordinates of a number of intermediate points of each zone. The single point G at the coordinates (xg,y ) in Figure 3 is illustrative of such points.

Since the zone tias some width in the y direction, the y coordinate of each intermediate point such as G can be conveniently placed either at 63 at the leading edge of the zone closest to the antibody well, at 65 at the trailing edge of the zone, or at one or more points such as 64 within the zone, typically including the point of maximum intensity.

As the zone develops with increasing time of incubation, the points E, F and G progressively change both their relative and their absolute positions. Illustrative values of xe and Xf as functions of time are plotted in Figure 4 for a typical protein concentration. Such position values may be used directly for determining protein concentration, as by comparing the values of xe and Xf obtained with unknown solutions to corresponding values obtained with known concentrations of proteins at measured times. However, more reliable and accurate results are ordinarily obtainable by deriving from such initial measurements one or more functions, which will be referred to for convenience as parameters.

An important parameter of the precipitation zone employed by the present invention is the coordinate difference Xf - xe, which is a measure of the length L of the precipitation zone at the particular time t of measurement. The variation of that parameter with time of incubation is plotted in Figure 5 for the typical data of Figure 4.

Position parameters other than the zone length L may be computed from measurements of the zone as it develops. For example, the zone curvature and its variation along the length of the zone are useful parameters, as well as providing information as to presence of protein abnormalities (see below). A rough measure of curvature for the zone arc as a whole, or for a selected zone segment, can be obtained by relatively simple comparisons of the x and y coordinates for three mutually spaced points on the zone axis. To obtain more accurate values of zone curvature, the zone axis is typically fitted approximately by a curve of the type y = f(x), where f(x) may represent any suitable function of x. The radius of curvature R is then given by the general formula

where y' and y" represent the first and second derivatives of y with respect to x.

A suitable illustrative function for curve fitting represents a parabola, typically with axis parallel to the y axis. Such a parabola may be expressed in either of the equivalent forms y = aO + a1x + a2x2 (2a) and y = A(x - B)2 + C (2b) where A = a2 B = -a1/2a2 C = a0 - al2/4a2 The values of the constants in (2a) or (2b) can be found directly from the coordinates of any three points of the zone axis, or may be fitted by least squares or other known procedure to any desired number of such points. The axis of symmetry of the parabola is at x = B, and the curve at that axis is a distance C from the x axis. Using the formula (1), the radius of curvature may be expressed as

That radius has its maximum value R0 at the axis of symmetry, where (3) reduces to R0 = 2A (4) Any of the above quantities, which may be derived as indicated from as few as three points of the zone axis, may be employed as parameters in accordance with the invention.

Another parameter that is useful for determining protein concentration in accordance with the invention is the time To of first appearance of the zone, That time is difficult to determine by direct observation. One aspect of the invention provides a practical way of obtaining a reliable and reproducible value for the time of first appearamce.

In Figures 4 and 5 the solid lines represent typical plots of direct experimental values. The figures also include extrapolations of the solid curves toward earlier times. The extrapolations are shown as dashed lines. The point 66 at which the extrapolated curves of Figure 4 meet represents a time at which the zone must have had zero length. The extrapolated curve of Figure 5 intersects the time axis at the point 67, giving an equivalent procedure for finding T,. Such extrapolation represents a reasonable and highly useful definition of the time of first appearance of the precipitation zone. That method of determining To has the advantage that continuous observation of the plate is not required.

Turning now to the quantitative determinations of the zone light intensity, it has been discovered that a single light intensity reading does not usually provide a useful measure of protein concentration. That is primarily due to the variability of zone shape and speed of formation, and the tendency of the intensity to change as the zone expands.

On the other hand, we have found that the variability of such factors can be largely compensated by taking a series of intensity readings at suitably selected locations and treating them collectively to evaluate an intensity parameter. A preferred procedure is to take such readings at uniform intervals along a line extending linearly across the zone of precipitation, typically in the y direction at a particular value of x. Such a series preferably includes several intensity values outside of the zone at each side. Those offset values are then averaged to provide a measure of the background intensity, which is subtracted from each of the intensity readings within the zone. The resulting adjusted intensity values are effectively summed, yielding essentially a linear integral of the intensity along a line crossing the zone at a selected value of x.We have found that such a linear intensity sum Ix tends to increase with incubation time in a regular and reproducible manner, the value at any given time increasing with the concentration of the reacting protein over a wide range of experimental conditions.

Typical plots of the relative intensity observed during cross-zone scans are shown in Figure 9. The two curves were plotted by semi-automatic apparatus of the type described below, scanning in the y direction at xg, the point of closest approach of the precipitation zones to the antibody well. The peaks at 74 and 77 in Figure 9 are due to the zones formed on opposite sides of an antibody well by identical protein samples in the two antigen wells of a plate similar to that of Figure 1. The two small peaks at 75 and 76 are due to the respective edges of the antibody well, providing a convenient reference from which to measure distances such as D from selected zone points to that well. The parameter Ix defined above corresponds essentially to the area under a peak such as those at 74 or 77 of Figure 9.

Peaks 74 and 77 of Figure 9 were made with identical samples of human serum albumin.

They illustrate typical development of the precipitation zones between 2 hours of incubation (solid line curve 70) and 4 hours (dashed line curve 72). Although the zone position remains remarkably stationary during the time between those two sets of measurements, the area of each peak grows appreciably. The near identity of the peaks at 74 and 77 is noteworthy. The two curves are offset vertically by an arbitrary distance for clarity of illustration.

Figure 10 is a schematic plot illustrating typical scans in the y direction on plates made with different concentrations of protein, all measured at substantially the same time of incubation. The graph brings out clearly the increasing area of the individual peaks and the shift of the entire zone toward the antibody well with increasing protein concentration progressing from peaks A to C.

Whereas the parameter Ix is highly useful for determining protein concentrations, still better results are obtained from a multiple intensity parameter, obtained by summing or averaging such linear intensity sums at several different x values, A typical procedure is to compute linear sums at xg and at values spaced on each side of xg by a selected interval.

Averaging or summing a uniform predetermined number of such linear sums reduces experimental error and improves the overall accuracy of the determination of protein concentration.

In accordance with another, generally preferred procedure, each time the plate is scanned the number of linear sums included in the computation is increased as the length of the precipitation zone increases. As illustrative procedure of that type is to determine the linear sum of the intensity at xg and to continue to compute such sums on each side of xg until the vaue of the sum falls below a selected threshold. Addition of all the linear sums provides a parameter Iz which is essentially the integral of the intensity of the zone of precipitation at the time of the scan. That approximation can be obtained as precisely as is desired, within the limits of resolution of the instrumentation, by reducing the x and y increments at which measurements are made.The value of Iz varies especially steeply as a function of protein concentration over a wide range of experimental conditions. That is because, as the concentration is increased, both the x and y dimensions of the zone increase. and the average intensity of the zone also tends to increase. That steeper dependence upon protein concentration makes the total intensity parameter especially effective as a criterion for determining the concentration.

After obtaining experimental values for one or more parameters for an antigen sample to be analyzed, those values are compared with suitable sets of standard values for the respective parameters, obtaind under closelv similar experimental conditions but with a series of protein solutions containing respective known concentrations of the protein of interest. To obtain such an array of standard values, standard runs are carried out with such standard protein solutions, and measurements are made at corresponding points of the respective plates at successive times as their incubation proceeds. The standard runs are preferably carried out with all conditions as nearly identical as possible to those of the experimental runs for which thev are to provide reference values.In fact, a distinct set of reference values is preferably obtained for every group of experimental runs. However. for routine measurements for which the reference curve slopes are known from previous experience, even a single standard ruI1 may sometimes suffice.

Standard values of the selected parameter are derived from the results of those standard measurementss typically for each concentration and at several times. Each measured standard value of the parameter is therefor considered as a function of both concentration and time. When individual linear intensity sums lx are to be considered separately, a full identification also requires specification of the value of x.

An advantage of using the time of first appearance of the zone as a parameter is that standard values of To do not involve time as a variable. That is, although the evaluation of -rt, by the methods described requires measurements at a series of definite times, once T, has been derived from those measurements the times of the respective measurements become immaterial. Thus. the values of 'lt, can be plotted as a function of protein concentration C, yielding a single standard curve. Such a curve is represented schematically in Figure 6, typically based on values obtained for respective concentrations by the extrapolation method described in connection with Figures 4 ind 5.Concentration values are directly readable from the curve of Figure 6 for an unknown simple once its I ?',, has been measured.

In the case of parameters such as L. found l, preparation of standard curves may be less direct. Since values depend upon the times nt which the measurements are made, the series of reference standards must be prepared in such form that they cover a range of times. It is not always feasible to measure data for all concentrations at the same moment. Each measured value is therefore associated with its time of measurement, and the resulting standard values for the respective protein concentrations are then plotted on separate curves as functions of time.

Figure 7 shows a typical family of such curves in which standard values of the parameter for three concentrations are plotted against time. It is noted in passing that the indicated extrapolation of those curves to L = () can provide standard values ofT,, for plotting Figure 6. or can provide experimental values of I, for comparison with Figure 6. Vertical lines are drawn on Figure 7 at a series of arbitrary times, shown at t1, t2 and t3. Their intersections with the curves then provide a set of L values for different concentrations, all corresponding to the same time.Each such set of l values is replotted as a separate curve as a function of concentration. The result is an array of curves, each showing 1, as a function of protein concentrations for a particular time. Such nn array is shown schematically in Figure 8, and is found more convenient for comparison with an experimental parameter value than the plot against time of Figure 7.

Standard values for comparison with experimental values of other parameters are typically obtained and treated in a manner analogous to that described for the parameter 1..

We have discovered, however that the total intensity parameter I,, when measured for given protein concentration at successive times of zone development, increases linearly with the time. That linear relation is illustrated in Figure l l, which is a plot of l, against time for four solutions containing the indicated known concentrations of the protein immunoglobu lin. The indicated values of lZ were derived by a suitably programmed general purpose computer from intensity values measured automatically at a two-dimensional array of zone positions in the general manner to be described. The points shown in Figure Il were originally plotted automatically, and the straight-line curves were fitted to each set of points, by the computer. The figure has been redrawn manually and is repn)duced at greatly reduced scale.

The linear relation shown in Figure 11 aids the preparation of standard or reference plots from which to read the protein concentration corresponding to an experimentil value ()f the parameter 17. One such curve, derived from the data of Figure 11 and showing I, as a function of concentration for the time 2.4 hours, is shown in Figure 12 at 7().

The linear dependence of I, upon time also means that the time derivative dl,/dt, or rate of change of Iz with time, is constant. It thus provides a parameter having the practical advantage that it does not depend upon time. Curve 72 in Figure 12 illustrates typical behaviour of that parameter as a function of protein concentration, each point representing the slope of one of the curves of Figure 11. As already explained in connection with the parameter To and Figure 6, a single curve such as 72 of Figure 12 can serve effectively as reference curve for the parameter dlz/dt.

A further aspect of the invention provides improved accuracy, especially when the sample under study contains the protein or proteins of interest at relatively low concentrations. The sample solutions are then preferably supplemented by adding known amounts of such proteins directly to aliquot portions of the sample itself, and running the resulting solutions in parallel with the original solution, and with standards containing definitely known concentrations of the protein. Values are obtained for the desired parameter for the respective solutions all at a uniform time, using the above described technique for interpolating with respect to time if necessary. The parameter values P obtained for the regular standard solutions are plotted against protein concentration in the usual way, as indicated schematically by curve 80 in Figure 13.Also plotted are the values of P for the original antigen sample, and for the portions of that sample to which protein was added in known amount, indicated typically as the respective points h, i, j and k. Those points are plotted relative to the horizontal concentration axis as if the solutions contained only the protein that was added. Thus the value h for the original sample is plotted at C = 0.

Curve 81 is drawn through those points. With that illustrative arrangement of the data, the protein concentration in the original sample can be evaluated in several ways, which should give essentially the same value. Thus, any desired number of those procedures may be used and the resulting values averaged.

First, horizontal projection from point h to intersect curve 80 at h' yields the concentration value indicated at Ch, which corresponds to the general comparison procedure previously described. Further, similar projection of each value i, j and k to intersect curve 80 provides a measure of the concentration in terms of the lengths of the lines from i to i', j to j' and k to k', all of which lengths are theoretically equal. If the value of P at he is somewhat uncertain, for example because the observed zone is faint, an average of all four intervals gives a more reliable value.

Another procedure has the advantage that it can not only improve the accuracy of an uncertain determination of h, but can also provide an answer even if the original sample contains less than the threshold value of protein needed to produce a measurable precipitation zone. In the latter case, the upper portion of curve 81 is drawn through the available points i, j and k, and is extrapolated to the P axis, thus providing a determination of P at h. In thus projecting curve 81, standard curve 80 is a helpful guide. For example, curve 80 is shifted bodily to the left in Figure 13 till it represents as well as possible the points i, j and k. The intersection of the shifted curve 80 with the P axis then gives a good measure of point h, from which the concentration Ch is found.Also, further extrapolation of curve 81 to the negative C axis at --Ch gives a direct reading for the concentration, which should agree with the value C11.

The described reliance upon the supplementary standard curve 81 has the great advantage of using the identical medium for the experimental and reference solutions.

Especially when that medium is the human serum of a particular patient, that advantage more than overcomes any possibility of slight error in extrapolating the standard curve. By increasing the number of different amounts of added protein, the three shown illustratively in Figure 13 being merely illustrative and by using several runs for each amount, the reliability of the extrapolation can be increased almost without limit, and a good indication can be obtained of the experimental error that it may involve.

The invention further provides methods for the automated detection of abnormalities of certain proteins in an antigen sample. Normal and abnormal gamma globulin, for example, typically have slightly different ranges of electrophoretic mobility, but react with the same antibody, leading to precipitation zones with irregular and typically unsymmetrical distribution of precipitin along the length of the zone. Such irregularity can be detected by comparing experimental values of the intensity parameter Ix, for example, for different values of x.Observation of unsymmetrical or otherwise irregular variation of such values indicates that an abnormality has been encountered; and that indication is greatly strengthened if several independent determinations are carried out at each x value and if the computed experimental errors show the variation to be statistically significant.

That method may be considered as a special case of the general procedure of comparing the experimentally obtained values of different parameters which tend to respond differently to the abnormality of interest. For some parameters with such behaviour it may be more effective to compare the values of protein concentration that correspond to the respective observed parameter values, rather than directly comparing the parameter values themselves. For example, the desired irregular intensity distribution along the zone in presence of abnormal gamma globulin tends to cause the zone to appear sooner than normal, leading to higher values of computed concentration; whereas the total intensity parameter Iz tends to yield approximately equal concentration values for the normal and abnormal proteins.Thus significantly poor agreement between concentration values obtained with To and with Iz indicates presence of abnormal protein.

A further parameter which responds sensitively to presence of abnormal protein is the curvature of the longitudinal axis of the zone, In presence of abnormality the curvature tends to vary irregularly and unsymmetrically along the length of the zone, while the overall curvature tends to be less than normal.

Further advantages are nearly always obtainable by employing multiple parameters which comprise specific functions of two or more parameters, such as the described position or intensity parameters. For example, the sum of the zone length L and an intensity parameter, each appropriately weighted according to its experimental error of measurement, provides a new parameter which tends to give more reliable results than either one of its components when used alone.

Another advantageous combination comprises a differential function of two parameters, one of which increases with increasing protein concentration, while the other decreases.

Thus, a quotient or a difference of such parameters provides a multiple parameter with a steeper dependence upon concentration than either of its components. The time To of first appearance of the zone is an example of a parameter having inverse relation to protein concentration, and is especially useful for constructing such quotients. The distance from the antibody well to the zone at any desired value of the x coordinate, for example the distance of closest approach at xg, also depends inversely upon the concentration and may be employed as a parameter in such quotients. It is generally preferred to divide the parameter with direct dependence upon concentration by that with inverse dependence, so that the resulting multiple parameter has a direct rather than inverse net dependence.

Suitable reference values for comparison to experimentally determined multiple parameters of unknown samples are typically derived from reference values obtained as already described for the respective component parameters.

Whereas the position and intensity measurements and the parameter derivations that have been described can be carried out by direct visual and manual operations, as has been indicated, an advantage of the invention is that those operations are particularly well adapted for partial or substantially complete automation.

An especially convenient and effective procedure for making the measurements required by the invention is to scan the slide by optical means, such as a television camera, a charge coupled scanning device or equivalent means, for producing a video signal representing the apparent brightness of the scanned image at a two-dimensional array of elemental areas.

That signal for each element is typically converted to digital form and electronically stored in association with signals representing the x and y coordinates of the corresponding position on the plate and the time of observation. The entire array of such position signals, or the signal for a designated position, can then be recovered from memory as needed for further processing. Also, any desired portion of the slide image can be displayed either during the scanning operation or at a later time by means of a cathode ray tube or equivalent display device. Systems for scanning, storing and reproducing an image, and for extracting a video signal in digital form for a selected image point are well known in the electronic art. and are available commercially in forms well adapted for the present purpose.

Figure 14 represents schematically such apparatus in illustrative form, comprising the television camera 90 with lens 91 for imaging the plate 20 on the photosensitive surface of the camera, the monitor cathode ray tube 92, the scan control apparatus 94 and the general purpose computer 100. Plate 20 can be imaged at any desired magnification by conventional adjustment of lens 91. A desired portion of the plate can be centered in the field by shifting its position on its illuminated support, such as light box 30 of Figure 2, or by conventional electronic bias controls in the camera circuits.

More particularly, the scanning movements in the scanning device and in the synchronized monitor cathode ray tube, if used, are typically driven stepwise in both coordinates, with steps so small that the movement appears virtually continuous on the screen. The image is thereby divided into elemental areas, which can be identified, for example, by specifying their x and y coordinates. For manual or semi-automatic operation, the monitor display typically includes a cursor comprising an electronically produced bright spot which designates on the screen the elemental area from which the video sample is being extracted once each complete scan. The operator is provided with manually controlled electronic switching means by which he can move the cursor about on the screen under visual control.Means may also be provided for directly placing the cursor at a desired point, typically in response to x and y signals in digital form constituting a command address. Such signals may be derived under manual control, or supplied from a computer in response to conventional program instructions.

As illustratively shown in Figure 14, the counter 102 counts pulses received from the clock 103, producing on the multiple lines 104 digital signals which represent successive x coordinate values, say. Circuit 106 develops on the line 107 an analogue step voltage in response to those signals, with a step corresponding to each elemental x value. The dividing circuit 108 effectively counts the beam sweeps in the x direction, producing on the multiple lines 109 digital signals representing the y coordinate for each sweep. Circuit 110 develops on the line 111 an analogue step voltage in response to that count, with a step corresponding to each elemental y value. The step voltages on lines 107 and 111 control the x and y sweep movements in camera 90 and tube 92, insuring their synchronization.The video signal from camera 90 is supplied via the line 112 to monitor tube 92, reproducing on the tube screen the brightness variations of slide 20.

The selection circuit 118 typically comprises x and y counters which are selectively shiftable up or down either by a single count or in a countinuing series of counts in response to manual movement of the joy stick indicated schematically at 119. The resulting digital signals, representing x and y coordinates of a point selected for exploration, are stored in register 117. The comparison circuit 114 continuously compares the coordinates of that selected exploration address from register 117 with those of the scanning beams from lines 104 and 109. When the scanning spot reaches the stored address, typically once during each complete scan, comparison circuit 114 supplies an enabling pulse to the switching circuit 120. The video signal from camera 90 is thereby transmitted to analogue-to-digital circuit 122.The register 125 thus receives on the lines 123 a digital representation of the brightness of the scanned slide 20 at the selected exploration point. The enabling pulse from circuit 114 is supplied also to the cursor control circuit 126, superimposing upon the video signal a beam intensifying pulse which identifies that selected point on the monitor screen, as indicated schematically at 127.

Register 125 also receives the x and y coordinate signals from lines 115 and 116 and a digital signal continuously representing the time, developed under clock control by the time circuit 128. All of those signals are typically store in register 125 for delivery to computer 100 in response to unload signals which may be supplied either under manual control via the line 129 from selector 118 or via the line 130 under control of the computer.

Computer 100 is preferably provided with means for independently designating an exploration address and thus obtaining input information for the spot of plate 20 that corresponds to specific requirements of any program under which it is operating. Such exploration means may comprise circuitry basically similar to selector 118, register 117 and comparison circuit 114, the selector operating in response to electronic signals which selectively indicate the required direction of movement of the exploration spot or the digital address of its required position. Rather than duplicating such apparatus, Figure 14 indicates computer control of selector 118 by electronic signals supplied via the multiple lines 132, supplanting joy stick 119 when the switch 133 is shifted from manual position to automatic.

With such control available the computer is readily programmed to carry out position and intensity measurements and parameter derivations of the general types that have been described.

For automatic scanning of zones, the plates of a set may be transferred in sequence from the incubation chamber to an accurately defined scanning position with suitable illumination. Especially when the optical scanning means is both light in weight and compact, as in the case of a charge coupled device, for example, it is generally preferred to mount that scanning device on a platform that is movable in two dimensions over the stationary array of immunodiffusion or immunoelectrophoresis plates. That arrangement facilitates use of the scanning system as an aid to fill the antigen and antibody wells. For that purpose a digitally controlled dispenser is mounted in fixed but offset relation to the scanning device.Video signals from the scanner are supplied to the computer, which is suitably programmed to recognize the shapes of the individual wells, or to identify machine readable symbols associated with them. As the platform carrying scanner and dispenser is moved under computer control over the plates, it can then readily be stopped when the scanner axis is located accurately over a selected antigen well. The platform may then be moved a fixed increment to center the dispensing tip over the well so that the correct charge is dispensed accurately into the well.

When all wells are filled in this manner with the appropriate antigen solutions, with suitable washing of the dispenser between operations, electrophoresis is initiated. A similar loading operation is typically carried out to charge the antibody wells, either after completion of the electrophoresis, or immediately after filling the antigen wells if electrophoresis is omitted. After the desired time for diffusion the same scanning device is moved successively to all positions where zones of precipitation are to be scanned throughout the array of plates.

The scanning capability of the scanning device is preferably used also to scan the array of plates before the wells are loaded. The computer is programmed to compare the observed positions of the anitgen and antibody wells, recording any departure from uniformity of dimensions or relative positions among the plates o the array. Any plate or individual well found to depart too far from standard can then be automatically omitted during the loading operation; and smaller deviations can be compensated by automatically applying suitable corrections to the parameter values that are ultimately derived by the computer.

Initial zone scanning procedure typically comprises request by the computer for video information at successive exploration addresses with a selected x value and with y values shifting progressively by a specified interval over the region in which zone formation is anticipated. The video intensity values received from register 125 are stored with x, y and t data for each explored point. After each such y-scan, the x coordinate is shifted by a specified interval and a similar y-scan is performed and the results recorded, until the entire area of interest has been covered. Each received intensity value is typically compared with the previously received and stored values. If the observed intensity variations exceed a set threshold characteristic of a background area, indicating presence of a precipitation zone, each intensity value within the zone is so designated in memory.Also, for each y-scan, the intensity maximum as the zone is crossed is identified and recorded, establishing the zone axis in terms of a series of y values at uniformly spaced x values. The end points of each zone are identified typically as the terminal x values in each direction at which an intensity maximum was identified. If more precise location of the end points is desired, the computer is instructed to perform further scans in a defined region about each end point with reduced x and y intervals.

With the existing precipitation zones so located and recorded, straightforward arithmetic operations are performed on the stored data, yielding the described intensity sum parameter Ix for each cross-scan of a zone in the y direction; and addition of those values gives the total intensity parameter Iz for the zone. Substraction of the x values at the zone ends gives the parameter L. Additional parameters may be obtained by appropriate computation as desired.

The described measurements are preferably carried out as a unitary operation on a complete set of plates that includes one or more unknown samples and also a set of related standard solutions of the type discussed above and sufficiently complete to permit evaluation of the unknown samples. After each scan of such a plate set, the computer is preferably directed to derive concentration values for each of the proteins for which precipitation zones were found. If, as is preferred, multiple standard runs are included, the probable experimental error can be derived for each calculated concentration value.

Computer routines for such calculations are well known. Several independent determinations are preferably made of the concentration of each protein in the sample, typically by reference to different parameters; an average is then computed, with each value weighted in the usual way according to its probable error.

If the computed probable error for that average is within the specified range for the particular sample under study, no further measurements are needed. The computer then produces a conventional printout or other record of the final results, together with as much of the original data as may be requested by the program. For example, the physician for whom the analysis is being carried out may require a complete copy of all the data on magnetic tape or the like for possible future reference. Also, photographs of the immunodiffusion plate can be taken, either directly or using the monitor tube.

Ordinarily, one cycle of scanning does not permit an adequate quantitative determination of a protein unless it is present in high concentration. After a suitable length of supplementary incubation, which may vary from a few minutes, say, to an hour or more, the above described scanning process is repeated, typically for the entire set of test plates and the corresponding standards. The computer is typically instructed to explore plate areas where zones were expected but not found during the previous scan; and also areas corresponding to the previously found and recorded precipitation zones, preferably allowing for specified expansion or movement of those zones which are set into the program on the basis of previous experience. Thus, the comuter operations can maintain a continuity of treatment of the respective protein zones between one scan and the next.

Occasionally two precipitation zones due to distinct proteins are so close together that their normal growth ultimately produces overlap. For zones produced by immunoelectrophoresis, such overlap normally occurs at or near the zone ends, leading to crossing of the zones as shown typically for zones a and b of Figure 1D. When immunodiffusion is carred out without initial fractionation, the overlap tends to be confined to the zone center portions, as when zones d and d' of that figure grow together. The computer is typically programmed to anticipate such a zone overlap, to recognize it when it has occurred, and then to make suitable modifications in the procedure used for deriving the various parameters.

When the proteins under study are such that overlaps are expected, each plate is typically scanned at an early stage of zone development before zone overlaps have occurred (Figure 1B). The zone axis and end points can then be located without ambiguity. The computer is typically programmed as part of the regular processing of each scan to check for actual overlaps, for example by comparing the (x,y) coordinates for each measured point of a zone axis with those for the adjacent zones, coincidence within a specified threshold indicating an overlap. All scans preferably also include a check for potential overlaps. For example, the computer extrapolates each zone axis beyond the observed end points and compares the extrapolated axis points of adjacent zones.Such axis extrapolation may comprise a simple linear extension in the direction of the axis slope near each end. Alternatively, the observed zone axis may be fitted by a parabola or other curve, which can then be extrapolated accurately.

Zone overlap can also be detected by computing the rate of change of slope of the zone axis, or the radius of curvature of that axis, at a series of points along the zone. Any departure from the normally smooth variation of those functions indicates that the zone may be made up of two or more overlapping zones.

When an actual overlap is found on one half of a zone, sufficiently accurate compensation can often be made by simply replacing the observed intensity values in the area of overlap with the values observed at the corresponding points of the other half of the zone. Such correspondence between two points for zones due to immunoelectrophoresis is typically defined as equal x spacing on opposite sides of xg of the point of closest approach to the antibody well; and for zones due to direct immunodiffusion is defined as equal y spacing on opposite sides of the zone axis. The axis within the area of overlap can usually be located by extrapolation.

More elaborate and accurate compensation procedures are also available, if desired. For example, the computer is instructed to adjust the measurement at the selected symmetrically placed point of the unaffected half of the zone to take account of any actual lack of symmetry of the zone. That unsymmetry can be determined, for example, by comparing the values that were obtained for the two points during a previous scan prior to the overlap, and applying the ratio of those values as a correction factor.

A further illustrative compensation procedure takes account also of the intensity value actually measured at the region of overlap, which is due, of course, partly to one zone and partly to the other. The computer is instructed to divide that observed value at each point of overlap in an appropriate ratio. A suitable approximate ratio may be obtained by comparing measured values for the two zones at the respective symmetrically placed points already described, either with or without the described adjustment of those values.

Whenever possible, it is preferred to employ two or more distinct computation procedures for compensating for areas of overlap, averaging the results and determining the probable error. It is emphasized, however, that the region of overlap is ordinarily a small fraction of the entire area of a zone. Hence even quite appreciable errors in approximating the true value within the overlap still may not affect the final result significantly.

Also, for each of the decribed types of overlap, it is usually possible to make the determination of protein concentration in terms of parameters that are defined by portions of the zone not affected by the overlap. Thus, the measurements preferably employed are those near the zone center for processes similar to immunoelectrophoresis, and those near the zone ends for processes similar to direct immunodiffusion.

It will be recognized by those skilled in the art that many particulars of the described procedures can be replaced by their obvious equivalents without departing from the proper scope of the present inventon. For example, values for those parameters that do not require a series of measurements at successive incubation times can be obtained from measurements made after the incubation process has been allowed to reach substantial equilibrium.

Such measurements have the advantage that the zone configurations are virtually static, and the exact time of measurement is therefore not crucial. As a further example, at an desired stage of the incubation the zones of precipitation can be stained in known manner with an appropriate dye, and zone measurements can then be made using selectively transmitted light. That method of measurement is usually most useful after equilibrium has been reached, since it is then not necessary to assign a precise time of measurement.

In summary, methods have been described by which immunodiffusion, immunoelectrophoresis and analogous processes can be used to provide truly quantitative measurement of the concentrations of individual proteins in biological fluids. Actual use of the described methods has demonstrated that protein concentrations can be determined within five percent or better. The invention thus provides an effective and convenient method for making a wide variety of experimental and diagnostic determinations.

WHAT WE CLAIM IS: 1. The method of obtaining a quantitative measure of the concentration of a protein in an antigen sample which has been subjected to immunodiffusion with an antibody source contaning an antibody specific to said protein, the protein and antibody diffusing into reactive contact generally parallel to an axis and through an area of a supporting medium initially free of both reactants to produce at least one elongated precipitation zone of limited length; said method comprising scanning the zone by optical means to develop electrical signals responsive to precipitin concentration at a two-dimensional array of positions distributed partly within the zonre and partly outside the zone, deriving electronically from said electrical signals the corresponding value of a selected parameter of the zone which varies with said protein concentration, and comparing the resulting parameter value with a set of reference parameter values derived correspondingly from reference zones produced by equivalent immunodiffusion of a plurality of reference antigen solutions con taining respective known concentrations of said protein to provide a quantitative measure of the concentration of said protein in said sample.

2. The method according to claim 1 wherein, prior to said immunodiffusion, said antigen sample has been subjected to selective protein migration in a direction generally perpendicular to said axis.

3. Method according to claim 1 or 2 wherein said parameter comprises the distance between the ends of said precipitation zone.

4. Method according to claim 3 wherein said step of scanning the zone is carried out transversely of the zone along a plurality of mutually spaced parallel scan paths which include scan paths crossing the zone and scan paths beyond the respective ends of the zone, and said step of deriving a parameter value comprises comparing said electrical signals to select scan paths adjacent the respective zone ends and determining the mutual separation of the selected scan paths.

5. Method according to claim 1 or 2 wherein said parameter comprises the integrated zone intensity, said electrical signals represent light intensity at the respective positions of said array, and the step of deriving a parameter value comprises summation of signal values at positions within the zone corrected for values at generally adjacent positions outside the zone.

6. The method according to claim 1 or 2 wherein said parameter comprises the time of initial appearance of the zone.

7. The method according to claim 1 or 2 wherein said parameter comprises the time of initial appearance of the zone, said step of scanning the zone is repeated at a plurality of selected times as the zone develops, and said step of deriving said parameter value comprises deriving from the electrical signals for each said scan a value of a second parameter which varies in characteristic manner as a function of time and has a known value at the time of initial appearance of the zone, and extrapolating the resulting values of said second parameter back in time to said known value.

8. The method according to claim 7 wherein said second parameter comprises the distance between the ends of the precipitation zone.

9. The method according to claim 7 wherein said second parameter comprises the integrated zone intensity.

10. Method according to claim 1 or 2 wherein said step of scanning the zone to develop electrical signals comprises forming an optical image of the zone at the image surface of a video camera tube, representing the intensity of the image at each said position as a digital representation, and storing in a memory associated with a digital computer the digital representations and digital addresses identifying the respective image positions.

11. The method according to claim 1 or 2 including producing on the screen of a cathode ray tube an image responsive to said electrical signals and representing said zone.

12. The method according to clajim 1 or 2 wherein said step of scanning the zone is repeated at a plurality of times as the zone develops, and said electrical signals represent light intensity measurements at a plurality of positions, x and y coordinate values at the respective positions and time values corresponding to the respective measurements.

13. The method according to claim 1 or 2 wherein said step of deriving a parameter signal includes deriving a first signal which increases with increasing protein concentration, deriving a second signal which decreases with increasing protein concentration, and deriving as a combined parameter signal a third signal which represents a differential combination of said first and second signals.

14. The method according to claim 13 wherein said combined parameter signal represents the ratio of said first and second signals.

**WARNING** end of DESC field may overlap start of CLMS **.

Claims

**WARNING** start of CLMS field may overlap end of DESC **. making a wide variety of experimental and diagnostic determinations. WHAT WE CLAIM IS:

1. The method of obtaining a quantitative measure of the concentration of a protein in an antigen sample which has been subjected to immunodiffusion with an antibody source contaning an antibody specific to said protein, the protein and antibody diffusing into reactive contact generally parallel to an axis and through an area of a supporting medium initially free of both reactants to produce at least one elongated precipitation zone of limited length; said method comprising scanning the zone by optical means to develop electrical signals responsive to precipitin concentration at a two-dimensional array of positions distributed partly within the zonre and partly outside the zone, deriving electronically from said electrical signals the corresponding value of a selected parameter of the zone which varies with said protein concentration, and comparing the resulting parameter value with a set of reference parameter values derived correspondingly from reference zones produced by equivalent immunodiffusion of a plurality of reference antigen solutions con taining respective known concentrations of said protein to provide a quantitative measure of the concentration of said protein in said sample.

15. The method according to claim 1 or 2 wherein said array of positions include at least

one position within said zone on one side thereof and in a region of overlap with another zone, and said step of scanning the zone to develop electrical signals includes developing an electrical signal at a second position which is correspondingly placed on the other side of the first said zone and in a region free of overlap, and employing said second position signal for deriving said electrical signal for said one position.

16. The method according to claim 1 or 2 wherein respective portions of said antigen sample have been enriched by adding different known quantities of said protein prior to said immunodiffusion of each sample portion, said method including subjecting the resulting zone for each sample portion to said steps of scanning, deriving parameter values and comparing the resulting parameter values with said set of reference parameter values.

17. The method according to claim 1 or 2 wherein respective portions of said antigen sample have been enriched by adding different known quantities of said protein prior to said immunodiffusion of each sample portion, said method including subjecting the resulting zone for each sample portion to said steps of scanning and deriving parameter values, and effectively plotting said resulting parameter values for the respective sample portions against the protein concentrations added and extrapolating the resulting curve to obtain a value for said protein concentration in the original sample.

18. The method according to claim 1 or 2 including the step of staining the zone of precipitation with a dye prior to said step of scanning the zone.

19. Method according to claim 1 or 2 including the steps of deriving electronically from said electrical signals the corresponding value of a further parameter of the zone which differs characteristically from the first said parameter in presence of certain abnormalities of said protein, and comparing the values of the first said parameter and of said further parameter to detect presence of such abnormality.

20. Method according to claim 19 wherein one of said parameters comprises essentially the integrated zone intensity and the other said parameter comprises the time of initial appearance of the zone, and the parameter values are compared in terms of the respective protein concentrations derived from them.

21. Method according to claim 19 wherein said one parameter and said further parameter comprise intensity sums extending across the zone at respective spaced positions along the zone length.

22. Method according to claim 1 or 2 including the steps of deriving electronically from said electrical signals the corresponding values, at a plurality of positions along the zone, of a zone parameter which is responsive to irregular or unsymmetrical precipitate distribution along the zone, and comparing said parameter values as a function of the positions along the zone length to detect presence of protein abnormality.

23. Method according to claim 22 wherein said parameter values represent the zone curvatures at a plurality of zone segments.

24. Method according to claim 22 wherein each said parameter value represents a summation of signal values at a plurality of positions distributed transversely of the zone.

25. Apparatus for obtaining a quantitative measure of the concentration of antigen in an antigen containing sample which had been subjected to immuno-diffusion with an antibody specific to said antigen, the antigen and antibody diffusing through a supporting medium from respective sources mutually spaced along an axis of the supporting medium and reacting to form at least one elongated precipitation zone of limited length; said apparatus comprising means for automatically scanning the supporting medium to develop electrical signals representing light intensity at respective positions of a two-dimensional array distributed partly within the zone and partly outside the zone, each light representing signal being associated with electrical position representing signals identifying the corresponding position in said array, and means responsive jointly to position representing signals and the corresponding light representing signals for deriving the value of a selected parameter of the zone which varies with said antigen concentration.

26. Apparatus according to claim 25 including means for causing said scanning means to repeat said scanning action intermittently, and means for associating electrical time signals with said light representing signals.

27. Apparatus for obtaining a quantitative measure of the concentration of an antigen in an antigen containing sample by subjecting the same to immuno-diffusion with an antibody source containing an antibody specific to said antigen, said reactants diffusing in a supporting medium from respective antigen and antibody wells mutually spaced in a plane of the supporting medium to produce a precipitation zone; said apparatus comprising photoresponsive means for scanning the wells and the supporting medium to develop electrical signals representing light intensity at respective positions, each light representing signal being associated with electrical position representing signals identifying the corresponding position, dispensing means for dispensing antigen and antibody containing solutions selectively to the respective wells in response to electrical command signals, and control means responsive jointly to the light representing signals and the position representing signals and having one mode of operation for producing operative position of the dispensing means relative to each individual well and for supplying a command signal to the dispensing means to supply that well with a selected reagent, said control means having another mode of operation for deriving electronically from the light representing signals and the position representing signals the corresponding value of a selected zone parameter which varies with said antigen concentration.

28. Apparatus according to claim 27, said control means, when in said one mode, including means for deriving signals representing the position of each well, and being responsive to such well position signals for producing operative position of the dispensing means.

29. Apparatus according to claim 27 including means associated with said control means, when in said one mode, for comparing the relative positions of said wells with predetermined standard relative well positions to detect deviations therefrom, and means associated with said control means, when in said other mode, for compensating for such deviations during said parameter derivation.

30. A method of obtaining a quantitive measure of the concentration of a protein in an antigen sample, substantially as hereinbefore described with reference to the accompanying drawings.

31. An apparatus for obtaining a quantitive measure of the concentration of a protein in an antigen sample, substantially as hereinbefore described with reference to and as shown in the accompanying drawings.