SURROGATES FOR TARGETS AND IMPROVED REFERENCE PANELS

Abstract: A method to determine reactivity of a candidate compound with a target which method does not require the physical presence of the target is disclosed. By providing a formula for treating data obtained from a reference set of target substitutes which formula is predictive of reactivity with the target, the compound to be tested can be physically assessed with respect to the reference panels, the formula applied, and reactivity with the actual target may be predicted. Panels which consist of individual members, said members comprising proteins, wherein at least one of the members of the panel is a protein other than an immunoglobulin (Ig) or fragment thereof and wherein the presence of said non-Ig protein enriches the panel are described. These panels can be tested for reactivity with an analyte to create a profile. Such profiles can be used in pattern matching, analysis of samples and other analyses.

Description: SURROGATES FOR TARGETS AND IMPROVED REFERENCE PANELS Technical Field The invention relates to identification of compounds that are useful in analysis, therapy and other applications where it is desirable to provide a substance which binds specifically to a target molecule, i.e., a specific pattern- matching technique which permits candidate binding substances to be screened in the absence of the target molecule. The invention also relates to an improvement in the construction of reference panels for use in profiling and pattern matching. Specifically, the invention concerns reference panels for the production of cross-reaction fingerprints which comprise enzymes and/or other nonimmunoglobulin proteins as affinity targets. Background Art There are numerous instances in which it is desirable to find a ligand that specifically binds a receptor or other target. To cite the most obvious examples, if a receptor is responsible for activation of a particular type of cell, ligands which bind the receptor may find therapeutic use in either activating or preventing the activation of the receptor, with a corresponding physiological effect on the cell. If the cell is contained in an animal or a plant, the effect may be felt by the entire organism. Thus, a very popular approach to designing new drugs rests on finding appropriate binding agents for these receptors. Ligands that bind specific targets can also find applications in analytical contexts. For example, antibodies are useful components in immunoassay procedures. All of these procedures rely on the specific interaction between an antigen and an antibody; either partner may be the analyte. In addition, separation procedures and other processes with industrial application may take advantage of specific binding. To take a very straightforward illustration, an impurity may effectively be removed from a composition by treating the composition with a solid support to which is bound a "receptor" capable of binding the impurity to the relative exclusion of the other components of the composition, provided the affinity of the receptor for the impurity is sufficiently greater than for the desired components. In all of the above cases, the amount of affinity that characterizes the specific binding and the degree of specificity required depends on the circumstances. Some applications are benefited by a relatively weak interaction, whereas others require a high affinity. Some applications are more demanding of specificity than others. The obvious brute force method to find a ligand that will bind a target of interest is physically to test the capability of a large number of compounds which are potential ligands with respect to their ability to bind the target itself. This method would no doubt eventually lead to finding a successful ligand in virtually every case but is clearly more time consuming and labor-intensive than would be desirable for practical utility. First, the target, e.g, . a receptor must be produced in some physical form that can be tested and sufficient quantities must be provided to test the range of compounds that are candidates. Second, if compounds are tested in just random order, a large quantity of target will be needed. This, especially in the case of cellular receptors, may be prohibitively expensive. Several approaches have been suggested to minimize these difficulties. First, rather than testing compounds at random, a systematically varied panel of compounds could be used. Such systematically varied panels can conveniently be constructed by forming polymers from monomer units of predetermined characteristics. The most convenient such polymers are peptides, but polysaccharides, polynucleotides and the like could also be used. The parameters that are important and the manner of constructing such panels are described in U.S. patents 4,963,263, 5,133,866 and 5,340,474 the contents of which are incorporated herein by reference. In addition to, or instead of, using systematically varied panels of compounds as candidates, the screening itself can be conducted in such a way as to minimize the number of physical measurements that are required. For example, as set forth in U.S. patent 5,217,869, which is incorporated herein by reference, a reactivity profile for a ligand known to react with a target can be established by providing a standard panel of binding agents. The profile obtained characterizes this particular ligand known to bind the receptor. The candidate compounds can then be tested against the same panel to obtain their corresponding profiles. When a corresponding profile matches that of a ligand known to be a successful binder to the target, the compound which generated the matching profile will have a high probability of binding the target. In an alternative, inverse image panels are prepared with varying characteristics, and profiles obtained for the receptor and ligand against opposite panels are matched. Various other technologies are directed to methods to improve the ease with which the physical binding of receptor to candidate ligand can be measured, such as the use of robotics, fluorescence detection of reactivity, physical arrangements of the panels, and so forth. Other methods which seek to find specific binding pair members include computer based methods such as three dimensional database searching, x-ray crystallography, molecular modeling, and the like. Other methods employ antibodies as surrogate targets or simply rely on the behavior of the compound with respect to related target receptors. For example, the behavior of a compound as an inhibitor of a particular serine protease, or of a number of serine proteases, might lead one to assume that it will be a useful inhibitor of an additional serine protease for which its inhibition activity has not yet been determined. The validity of this last mentioned method relies on the similarity of the serine proteases that are the "reference receptors" for which the binding characteristics of the test compound are known to the target receptor (serine protease) for which the binding characteristics are not known. The present invention provides another method to match ligand with a target. It is especially helpful when limited supplies of the target are available. The invention method is especially useful in drug design projects where the target has never been fully purified, is unstable or otherwise not available in adequate quantities for large-scale screening, or when the assay procedure for the target is complex and costly. Further, the method minimizes consumption of receptor in a program of screening against many potential ligands. U.S. Patent No. 5,300,425, incorporated herein by reference, describes methods of preparing characteristic profiles of a particular analyte, matching similar profiles to correlate binding properties among various analytes, and the use of inverse image panels to create profiles for this purpose. In the methods described in the '425 patent, immunoglobulins or their immunologically reactive fragments were used as members of panels of binding ligands to obtain the characteristic profiles used in characterization and correlation. A modification of this technology, described in U.S. Patent No. 5,340,474, mentioned above substitutes panels of diverse paralogs for the antibodies and fragments used in the profiling panels. Paralogs are defined as polymeric moieties preferably of MW less than 7.5 kD composed of monomers with characteristics such that maximal diversity could be obtained across the panel members with a minimum number of paralogs. By maximizing diversity, the range of space/charge contours that characterize "chemical space" can then be achieved with relatively small numbers of compounds. As described in the above-referenced patents, such reference panels are useful in a number of contexts. The panel can be used to obtain a "fingerprint" that characterizes a particular analyte. The fingerprint can be used as an analytical tool to identify a particular substance much in the same way that an IR spectrum or NMR spectrum could be used. In addition, it was recognized that analytes that have similar fingerprints or similar features contained in their fingerprints have similar binding or reactivity properties in general or with respect to the property associated with the similar feature. Therefore, if, for example, a receptor of interest has a known ligand, other compounds that will bind to the receptor can be found by matching their fingerprints against the reference panel with the fingerprint obtained from the known ligand. Similar matching of complementary members of a binding pair can be obtained using inverse image sets wherein a fingerprint for a ligand against a reference panel will match the fingerprint of the receptor against a set of compounds which is an inverse image of that reference panel. Still another application for which panels of reagents are useful is in determining analyte composition of a sample. This application is described in U.S. 5,338,659 incorporated herein by reference. The fingerprint obtained for an unknown sample is matched with predetermined fingerprints or profiles determined on standard known compositions. Certain computational techniques can be employed to facilitate this comparison as described in this patent. In this case, however, it is not generally thought that a wide range of binding capabilities will be required since the application is focused on compositions which contain analytes, generally with related structures, and means for correlating the fingerprints with the other inherent properties of the analytes themselves are not needed. Thus, in this case, it might be considered logical to use panel members which are not necessarily antibodies nor maximally diverse paralogs. In the present application, an additional method of identifying binding partners is described using a computational combination of results against a reference panel as a surrogate for a desired target. The reference panel illustrated is comprised of enzymes. It is thus found, surprisingly, that sufficient diversity of reactivity can be obtained to achieve meaningful results, even though the enzymes used in the illustrative reference panel were not designed by nature to have a vast multiplicity of binding activities (as are antibodies) . Neither were the enzymes expected to have the maximal diversity ascribable to a small number of panel members that was achieved through the design of paralogs. Nonetheless, by utilizing enzymes, even isoenzymes with similar activities, as members of the reference panels, a satisfactory surrogate can be achieved to predict binding of candidate ligands to targets, including targets entirely unrelated by any similarity of amino acid sequence to the enzymes that are panel members. It has thus been found that such enzymes or reference molecules in general should also be useful in the profiling and pattern-matching can thus be searched computationally with results comparable to a direct physical screen of the ligands. Thus, for each compound subsequently tested, reactivity against each member of the reference panel is obtained and the formula derived from the training set is applied to obtain a predicted value with respect to the target. Rather than directly testing the reactivity of a candidate compound with a target, it is possible instead to test its reactivity with respect to a panel of readily available reference receptors, apply the formula to the results, and predict what would have happened had the target itself been used. The larger the library of stored ligand profiles against a reference set, the larger the increase in efficiency for screening by surrogate. It has also now been found that nonimmunoglobulin proteins (some of which are naturally occurring, but regardless of how they are actually produced) can be used successfully to constitute a reference panel for use in profiling analytes, predicting binding capabilities of candidate compounds with respect to targets, as well as for the analytical purposes described in U.S. 5,338,659. Thus, panels useful in the methods of the present invention can be comprised entirely of such proteins as enzymes, T cell receptors, olfactory receptors, lectins, and artificially modified proteins containing arbitrary binding sites. The panels may also include antibodies or fragments thereof, or paralogs as members; however, in the panels useful in the methods of the present invention the non-Ig/nonparalog members must "enrich" the panel beyond the contribution of any immunoglobulin proteins and paralogs also contained in the panels, as described hereinbelow. In one aspect, the invention is directed to a method to determine the ability of a candidate compound to react with a target which method comprises providing a surrogate for the target. The surrogate is that formula representing a computational combination, preferably a linear combination, of at least 2 reference reactivity profiles, which best agrees with the empirical binding data of the target against the training set of compounds. The reference reactivity profiles represent the reaction of each member of a panel of reference receptors with respect to a set of compounds, which set of compounds can be designated a "training set". The formula is then applied to the reactivities with respect to each of the members of the panel of reference receptors that is obtained for each candidate compound. The outcome of applying this formula mimics what would be found had the compound been tested directly with the target receptor. This aspect of the invention thus relates to a method to identify a candidate reactive with a target, which method comprises: (a) providing a formula that represents a combination of the reactivity profiles of at least two members of a reference panel with respect to a first set of compounds, which formula calculates a predicted profile that best matches the reactivity profile of the target itself with respect to said first set of compounds; (b) testing the reactivity of said at least two proteins of the reference panel with respect to a candidate; and (c) calculating a predicted reactivity with respect to the target for said candidate by applying said formula to the reactivities determined in step (b) to estimate the reactivity of the candidate with respect to the target. A successful candidate is then identified and synthesized from the appropriate starting materials. Another aspect of the invention is a particularly preferred combination of a training set and panel. In this preferred matrix, each member of the reference panel has effectively an inverse image member in the training set of compounds. In this way, the number of reference panel members and training compounds is minimized by removing redundant overlaps. The invention is also directed to a database of fingerprints obtained with respect to a reference panel. The database can be used for a variety of purposes as described below. In still another aspect, the invention is directed to methods to construct the reference panels of the invention. In another aspect, the invention is directed to a method to characterize a single analyte, which method comprises contacting said analyte with each member of a panel enriched by or constituted by the above-described non Ig proteins which react in a multiplicity of differing degrees with said single analyte; detecting the degree of reactivity of said analyte to each of said members; recording said degree of reactivity of said analyte to each of said panel members; and arranging said recorded degrees of reactivity so as to provide a characteristic profile of said analyte. In another aspect, the invention is directed to panels containing non Ig proteins useful in the various pattern- matching methods and to physical embodiments of the fingerprints obtained by the four methods. In another aspect, the invention is directed to a method to identify a candidate, which candidate will be effective in reacting with a target, wherein said target has a known ligand with which it reacts, which method comprises: contacting said candidate with each member of the panel enriched by or constituted by the above-described non Ig proteins which react in a multiplicity of differing degrees with said candidate; detecting the degree of reactivity of said candidate to each of said panel members; recording each said degree of reactivity of said candidate to each of said panel members; arranging said recorded degrees of reactivity so as to provide a characteristic profile of said candidate; comparing said profile to a profile analogously obtained of said ligand with respect to said multiplicity of panel members; wherein similarity of the profile of said candidate to the profile of said ligand indicates the ability of the candidate to react with said target. A substance identified as a successful candidate is then identified and synthesized from the appropriate starting materials. In a third aspect, the invention is directed to a method to select a candidate, from a multiplicity of candidates, that reacts specifically with a known target, which method comprises: providing a profile of reactivity of said target against a maximally diverse set of compounds; providing a panel including non Ig proteins as described above which is an inverse image of said maximally diverse set; preparing a profile of the reactivity of the candidate to the inverse image panel; comparing the maximally diverse set profile of the target with the inverse image panel profile of the candidate; and wherein similarity of the inverse image panel profile with the diverse set profile indicates the probability that the candidate will bind to the target. A successful candidate is then identified and synthesized from the appropriate starting materials. This method can be "reversed" in that the choice of which substance is considered a "candidate" and which a "target" is arbitrary -- i.e., the target can be profiled vs. the inverse image panel and the candidate vs. the maximally diverse panel. In addition, the invention is directed to a method to determine the ability of a candidate to react with a target which method comprises providing a surrogate for the target, as described above, and including non Ig proteins in the reference panel. Brief Description of the Drawings Figure 1 is a flow diagram of the method to calculate the probability of a candidate binding to target using a surrogate. Figure 2 shows a preferred embodiment of the training set/reference matrix. Figure 3 is a flow diagram of the method to determine the profile of an analyte. Figures 4a and 4b show typical embodiments of fingerprints obtained by various invention methods. Figure 5 is a flow diagram of the method for comparing profiles of a ligand with a candidate compound. Figure 6 is a flow diagram of the method to compare inverse image profiles. Figures 7a, 7b and 7c represent distance distributions for profiles of 800 compounds determined with respect to reference panels of 5, 7 and 10 reference proteins, respectively. Figures 8a-8c are distance distributions for points in 10-dimensional space representing profiles of 50, 100 and 1000 compounds, respectively, with respect to a panel of 10 reference proteins. Figures 9a-9c show distributions for the profiles with respect to 10 reference proteins of various collections of compounds. Figure 9a is the same as Figure 8a which shows the distance distribution representing profiles of 50 random compounds. Figure 9b shows the distance distribution for profiles of 50 known pharmaceutically active compounds. Figure 9c represents a similar distribution for 50 peptides of varying biological activity. Figure 10 shows the results obtained when a training set of compounds is tested with respect to a panel of reference GST isozymes to generate a surrogate for a target receptor. The results of testing a multiplicity of additional compounds against the panel of reference enzymes and applying the formula defining the surrogate is compared to testing the additional compounds directly against the target receptor. Gray scale indicates IC 50 values. Figure 11a shows the predictions and actual empirical data from Figure 10 as a scatter plot indicating high degree of correlation. Figure lib shows the residuals from Figure 11a. Figure 12 shows a list of 122 compounds and their symbols used as the compound library in the results obtained against an enriched reference panel of eight selected enzymes. Figure 13 shows the experimental and predicted ability of the compounds of Figure 12 to bind GRd and AdDH, as well as the characteristic profiles of these compounds against a reference panel where the first 12 compounds listed are the initial training set. An additional set of 10 training compounds used in the second iteration predictions are denoted by adjacent black bars, with a different set of 10 for each target. Figures 14a and 14b show correlation plots of predicted and experimental values according to the results shown in Figure 13. Figure 15 shows the correlation between fitted and experimental binding of a multiplicity of compounds against nine different targets. Modes of Carrying Out the Invention The method of the invention permits a large number of candidate compounds to be tested for their ability to react with, and in particular to bind to, a target without necessity for large amounts of the target per se. The target itself is required only in sufficient quantity and purity to generate the formula which creates the surrogate. As used herein, the term "target" includes, for example, molecules that reside on the surface of cells and mediate activation of the cells by activating ligands, but also is used generically to mean any molecule that binds specifically to a counterpart. One member of a specific binding pair could arbitrarily be called a "receptor" or "target" and the other a "ligand". No particular physiological function need be associated with this specific binding. Thus, for example, a "target" might include antibodies, immunologically reactive portions of antibodies, molecules that are designed to complement other molecules, and so forth. Indeed, in the context of the present invention, the distinction between "target" and "ligand" is entirely irrelevant; the invention concerns pairs of molecules which specifically bind each other noncovalently with greater affinity than either binds other molecules. However, for ease of explanation, the invention methods will often be discussed in terms of target, such as an enzyme (again, simply a molecule for which a counterpart is sought that will react or bind with it) and "ligand" simply represents that counterpart (such as a low molecular weight inhibitor) . The Use of Surrogates In order to practice surrogate method of the present invention, the following elements are needed: First, a reference set of model targets against which measurable activity can be assessed. Various techniques for determining reactivity of compounds with this set of reference targets are possible, and within the skill of the art as described above. It is important to emphasize that it is unnecessary that the reference panel members be in any way related by primary amino acid sequence or empirical chemical structure or by known biological function to the target for which they provide a model. For example, in the illustration below, various enzymes, including glutathione S-transferase (GST) , are used as the reference receptors while the actual target is glutathione reductase (GRd) , aldehyde dehydrogenase, or a variety of other proteins. There is no previously discernible similarity between the enzymes of the panel and any of the targets at the levels of primary structure or of known enzymatic function. One of the advantages of the present invention is that the reference proteins can be quite different in known reactivity and in primary structure from the target, because the predictive information is present in their relative correlations with the target, not their homology. The reference panel may contain as few as 1, but preferably 2-50 and more preferably 8-25 non-Ig proteins; the total number of panel members can also be similarly described. Second, a training set of ligands representative of the compounds desired to be further tested with respect to their reactivities with the reference panel is required. If there is a library of compounds to be further tested, a multivariate clustering method can be used to determine representative compounds from the library, or similar to those in the library, for use in the training set. Similarly, compounds with maximally systematically varying properties can also be used. In general, this training set of compounds should include at least as many compounds as the number of reference proteins and preferably about 3 times that number. Third, there must be enough target available to test the training set empirically, although the target need not necessarily be pure. The target must be free of undesired interfering impurities, however. With these compounds and reference panels in hand, the profiles of each reference panel member with respect to the training set and the profile of the target with respect to the training set can be obtained by physical measurement. A fourth requirement then is a fitting procedure to match the target's profile with a combination of the reference panel member profiles. In addition to techniques for linear regression, nonlinear regression methods can also be used for this purpose, including partly linear models as well as rule-based methods such as clustering by recursive partitioning. Indeed, any algorithms used in hemometric analysis or pattern recognition generally can be combined with the physical assay data, represented by fingerprints prepared as taught here, in order to classify compounds. Such mathematical techniques are well understood in the art, and result in the formula which serves as a surrogate for testing of further compounds. Application of the formula to the profile obtained for a newly tested compound with respect to the reference panel results in an estimate of the ability of the newly tested compounds to bind target. Of course, this represents a probability and not an absolute. The predicted result amounts to a screening procedure to identify compounds with a high probability of binding the target (or not binding the target) . While one compound at a time can be tested with respect to the reference panel and the formula applied to estimate a target reactivity value, the most useful application of the method of the invention pertains to screening libraries of candidate compounds. Thus, quite frequently, a large number of candidate compounds is available and the method of the invention can be used to select those which do and those which do not bind the target. When the method is thus applied to libraries, the results from the newly screened candidates can be added, if desired, to the training set and the process repeated in an iterative loop. Thus, the original training set could be supplemented with selected compounds which are estimated to bind the target strongly and selected compounds which are estimated to bind the target only weakly or undetectably and these compounds used in addition to, or instead of, certain members of the training set to obtain the profiles with respect to reference panel members and actual targets. The formula can then be recalculated taking account of these additional members. Further, not all profiles of the reference panel proteins with regard to the training set need be included, in the end, in the formula. That is, some of the coefficients for model receptor profiles in the linear combination may be zero or negative. The general approach to the use of surrogates is outlined in Figure 1. In Figure 1, a fingerprint database is first assembled according to the procedure shown in Figure 3 described below for a multiplicity of compounds against a representative reference panel. The reference panel itself will have been selected using preliminary data to include members that have the ability to, collectively, react with a wide range of compounds but wherein each panel member reacts with different sets of such compounds. When a suitable panel has been chosen, a training set is also selected from among the profiles for testing against the target. Each of the members of the training set is thus tested and the resultant with respect to target is obtained for each member of the training set. This amounts to a profile of the target using the training set as panel members. The fingerprints of the training sets can then be inverted conceptually, since the same data points are involved, to provide a profile of each member of the panel with respect to the compounds of the training set. These conceptually inverted profiles can be analyzed mathematically, for example, using linear regression analysis, to obtain a mathematical surrogate as shown in Figure 1. The profile of any candidate, including candidates for which profiles are already available in the database, can be mathematically treated according to the surrogate to predict the reactivity with the target. Successful candidates can be identified using the surrogate-generated predictions and the successful candidates synthesized using the relevant starting materials. There is also a feedback loop which permits such predictions to be tested and revisions to the training set made on the basis of these predictions leading to modifications of the surrogate. Illustrative Method The method of the invention can be further illustrated using a simplified hypothetical matrix, and a linear regression method of combination. The matrix set forth below represents a hypothetical matrix used to illustrate the generation of the relevant formula as surrogate. Across the top labeled MR1-MR5 are five panel members which represent panel members, such as enzymes used as reference model targets for the actual target receptor TR. Along the side, labeled TC1-TC5 are five training compounds which bind or otherwise react in varying degrees with each of the reference panel members The degree of reactivity is arbitrarily assigned a value on a scale of 1-10 where 10 indicates high reactivity and 1 indicates low reactivity. Generally, a logarithmic scale of measured values is used. Sample Matrix MR1 MR2 MR3 MR4 MR5 TCI 6 1 1 7 2 TC2 2 4 2 6 2 TC3 1 3 8 1 5 TC4 5 9 10 10 1 TC5 9 1 10 5 9 In these hypothetical results, profiles for each of the set of training compounds with respect to the reference panel are shown in the horizontal rows and profiles for each reference enzyme with respect to the training set of compounds are shown in the vertical columns. Thus, for example, for MR1, there is a moderately high level of reactivity with TCI, low reactivity with TC2, very low with TC3, moderate reactivity with TC4 and very high reactivity with TC5. Thus, each of MR1-MR5 has a particular profile of reactivity with regard to the training set. On the right, marked TR, the target receptor shows a profile against the training set with monotonically increasing reactivities over the TC1-TC5 range, a pattern grossly different from any of the reference profiles. A formula is then generated by assigning weights to each of the elements of the five MR1-MR5 profiles to obtain a predicted target receptor profile that matches that actually obtained for the target. The weighing values will need to be the same for each element of the profiles. Thus, the weights applied to the TCI element with respect to how the values from MR1-MR5 are counted have to be the same as those applied to TC2. Ultimately the algorithm will be of the form A(MR1) + B(MR2) + C(MR3) + D(MR4) + E(MR5) = the value assigned to the predicted value according to the surrogate, shown in the table as PR. Each of the coefficients A-E will have a numerical value; some of the coefficients may be zero. This same equation, with the same values of A-E will be used to calculate the predicted reactivity with the target receptor for any individual candidate compound. In the above example, A=+2; B=+3; C=-l; D=-2; E= + l. Here the coefficients allow a perfect match between the Predicted Receptor (PR) profile and the target receptor (TR) profile with respect to the training set. In general, and if more compounds are included in the training set a perfect match may not be possible; but the closest approximation obtainable is useful to the same end. Thus, for any new compound, a prediction for reactivity with target is obtained as follows: A profile that provides reactivity values for MR1-MR5 is obtained. The values obtained are then substituted into the formula set forth above, with the predetermined values of A-E. A predicted value is calculated. Thus, a new candidate compound, which gives a profile with values of MR1=8, MR2=9, MR3=4, MR4=7 and MR5=5, will be evaluated according to the formula: (+2) (8) + (+3) (9) + (-1) (4) + (-2) (7) + (+1) (5) = PR to provide a predicted reactivity value of 30. This demonstrates that the method can predict higher reactivity than available in the training set. Confirmed high reactivity compounds can be added to the training set to refine the formula. Examples 3 and 4 set forth below indicate that this general approach is successful in predicting the reactivity of any candidate compound with a target; accordingly, no further supplies of target receptor are required in order to test an arbitrary number of compounds. In a preferred embodiment of the original matrix, both the reference panel and the training set are maximally diverse and represent inverse images. This is illustrated in Figure 2 which shows a hypothetical matrix of reference panel members and reference binding agents. As illustrated in the figure, reference panel member 1 and set member 1' interact strongly; reference panel member 2 and set member 2' do so; reference panel member 3 and set member 3', etc. There is relatively weak interaction between, say, set member 3' and reference panel member 2 or reference panel member 1. In effect, the reference panel and the training set represent inverse images. Kits can be prepared which include, in separate containers, each of the members of the training set, each of the members of the reference panel, and the target, along with reagents for testing their reactivity. Inclusion of Non Ig Proteins Performance of the above surrogate method resulted in the surprising discovery that fingerprint matching to identify compounds with desirable properties, such as the ability to bind to a desired target, the ability behave as an enzyme inhibitor, a specific pharmacological activity, and so forth, can be based on panels which are substantially enriched by proteins that are neither immunoglobulins and their fragments nor specifically designed maximally diverse paralogs. Surprisingly, a range of complementarily or other interactive ability sufficient to cover substantially all of "chemical space" can be achieved by employing naturally occurring proteins such as enzymes, lectins, T cell receptors, olfactory receptors and the like, or by employing proteins which are modified forms of naturally occurring proteins. By choosing a suitable set of these proteins, a sufficient range of reactivity can be obtained to provide enhanced fingerprints in these contexts. Thus, panels enriched by or constituted by nonimmunoglobulin proteins serve to provide suitable reference sets of data points for obtaining a characteristic profile of an individual substance. The profiles can be manipulated in a number of ways as further described below. It would be possible and is within the scope of the invention to construct panels which contain as members not only these proteins but also antibodies and/or paralogs or other arbitrarily chosen quantitative reactivity events. When the word "reactivity" is used in the present application, it refers to noncovalent interaction between the stated participants. In a sense, then, "reactivity" is substantially similar to noncovalent binding. Such binding may or may not be coupled with catalytic or allosteric responses. However, the panel must at least be enriched by the alternative proteins. A protein "enriches" the panel if its membership in the panel does any of the following or some combination: (a) expands the coverage of the panel over chemical space (see below) ; (b) increases the average distance between fingerprints of different compounds in the library (see below) ; (c) decreases the number of reference panel members required to obtain a given number of principal components (see below) . (a) It is, of course, desired to cover all of chemical space. However, 90%, but preferably 95%, coverage is generally satisfactory. "Covering" chemical space means that all compounds tested against the panel show at least some reactivity with at least one panel member, and preferably 3-5. (b) The "distance" between fingerprints or profiles can be best understood by the device of assigning each profile The proteins in the panel can preferably be chosen as follows: An iterative process is used to select the members of any panel for use in fingerprinting. A few candidate panel members including non-Ig proteins are arbitrarily chosen and fingerprints for any arbitrary set of compounds are obtained. Comparisons are made between the fingerprints. Any method of comparison could be used, but some particularly effective methods are described hereinbelow. Whatever the method of comparison, compounds which have very similar fingerprints are clearly redundant members of the library of compounds for this purpose and only one of the compounds in such a group should be retained in the selection set. The remaining fingerprints are then again compared for similarity, only this time an inverse profile is obtained for each of the reference panel members with respect to the remaining compounds in the selection set . Now it becomes possible to discard panel members which provide similar inverse profiles with respect to the compound library. Thus, if three candidate members in the panel seem to provide similar reaction patterns across the compound library tested, only one of the members is retained in the panel. If the panel, including non-Ig proteins, which has been thus reduced for redundancy continues to generate viable fingerprints for all new compounds and if the new compounds do not reveal any further redundancy in the panel, then the panel is satisfactory. However, if the panel fails to provide a meaningful fingerprint for any new compound, additional members need to be added to the panel, although it becomes harder and harder to find a new member which provides distinct patterns as compared to those already present. The screening for new members in the panels is preferably conducted on compounds that were not detected with the members already present. The new member candidates are then evaluated on a maximally diverse set of the compounds already tested. The ideal panel provides high coverage with high independence and a small number of members, preferably under 100, more preferably under 25 and most preferably under 15. It has been found that among 100 enzymes of widely varying function, 12 of them provide 95% coverage against 1000 compounds from a wide variety of chemical classes. The 12 enzymes are independent since about 9 statistically meaningful principal components are needed to describe the 12; if they were totally independent, 12 would be needed. Arrangement of the Panel The members, including non-Ig proteins, that comprise the panel must be physically embodied in such a way that an individual result for each member can be retrieved and recorded so as to construct the profile. Of course, it is possible simply to react each member independently in an individual reaction container with the relevant analyte; to record the results of each container individually; and to manually construct the profile that results. More convenient alternative approaches involve displaying the panel members in an orderly fashion on some type of solid support, such as a microtiter plate or other support with multiple test regions and to scan the regions for the individual results. The scan can assess the results in each region sequentially or simultaneously using known technology. In general, the reactivity of the analyte with each test region or container is assessed in terms of the binding affinity of the analyte to the panel member contained therein. The art is replete with the methodologies for detecting the degree of binding of one substance to another. In a prototypic approach, one partner, in this case the panel member, is bound to solid support and the other partner, in this case the analyte, is labeled using radioisotopes, fluorescence, enzymes and the like, and after contact of the analyte with the supported panel member, the support is washed free, if necessary, of unbound analyte and the amount of label measured. In the alternative, the binding affinity can be measured by competition between the analyte and a labeled competitor. One method of such competitive binding described in the above- referenced patents involves competition between the analyte and a diverse mixture of labeled compounds which mixture is sufficiently diverse that the mixture binds uniformly to every member of the test panel so that the diminution in label directly gives a measure of the degree of binding for the competitor analyte. Methods are also available to detect the degree of binding between two substances in homogeneous media as in, for example, the EMIT technology. In all of these methods, any conventional method of labeling may be used. Among preferred methods is the use of fluorescent label competition, for example using fluorescence polarization. The invention does not concern specific methods of detecting the degree of binding, and any conventionally used procedure for measuring the binding affinity between analyte and the member of the panel can be used. It is preferable to use assay methods with a wide dynamic range. Quantitation of affinity by IC 50 for inhibition of substrate turnover, or other competitive binding events, can often be measured over more than five log units of potency, for example. Profile Determination Determination of a characteristic profile provides the basic tool for the matching techniques of the invention. Each profile or fingerprint is determined by measuring the individual reactivities, such as binding affinities of the analyte for each member of the panel. The reactivities are then recorded in an orderly arrangement so as to provide this characteristic profile. Figure 3 is a flow chart showing the steps in obtaining the characteristic profile for an analyte. First, the analyte is contacted with each panel member (panel member i) in a panel of n members. For each of these contacts, the reaction of the analyte with the panel member is detected and measured. Then, the extent of reaction is recorded to obtain a data point for reactivity associated with each of the n members of the panel . Then the recorded data points are arranged in an orderly manner to obtain the profile. One convenient way of arranging these data points is to plot each reactivity in one of the dimensions of n-dimensional space. 00 © o m -5 ∞ m o unnecessary to know what the complementary motifs are in terms of their molecular structure. U.S. 5,338,659, referenced above and incorporated herein, discloses a particularly efficient approach to making comparisons between profiles. This approach is to plot the obtained profiles or fingerprints in n-dimensional space, wherein n is the number of members of the relevant panel and the location of the point in each dimension is a function of its reactivity with each panel member. The proximity of the points representing the unknown and any of the predetermined profiles in n-dimensional space represents the similarity of their compositions. Multiparametric statistical techniques can also be employed to define which of the n dimensions have the greatest information content relative to the assay so as to permit a selection of the minimum number of characteristics or dimensions to be measured. In order to use the profiles as tools in predicting properties of test substances or in other pattern-matching applications, regardless of the specific pattern-matching techniques used, the reference panel should be capable of covering at least 90% of the chemical space, and should provide an average distance between fingerprints of all pairs of at least about three times the noise level generated by replicate determination of profiles for a single compound. In addition, the fingerprint provided by the reference panel should provide at least five principal components with respect to the range of sms.ll organic compounds that are available commercially. For example, this range is typified by any set of approximately 1,000 compounds among those available from the Aldrich Catalog of Fine Chemicals. Use of Surrogates The panels are also used to create surrogates for a desired target in order to evaluate binding of candidate compounds as described above. Applications Applications of these pattern-matching processes with respect to the profiles or fingerprints are manifold. For example, it is possible, because of their ease of synthesis or because of their native occurrence, to obtain peptides or proteins that behave in biologically important ways. However, peptides and proteins are not attractive as drugs as they cannot easily be orally administered and metabolized and present problems in manufacturing and storage; small molecules are preferred. By matching profiles, either by direct pattern- matching, inverse panels, or surrogates, suitable small molecule substitutes can be found. Another important application is the prediction of toxicity in candidate drugs. Comparing the appropriate aspects of the fingerprint of the candidate drug with features of fingerprints of known toxins permits such prediction. Likewise, construction of surrogates for proteins similar in sequence or function to target allows side effects due to cross-reactions to be estimated in advance of animal testing using only trace amounts of the related protein. Still another application of the profiles of the invention and their correlation relates to providing parameters for improving the three-dimensional models of spatial arrangement of pharmacophores obtained by conventional computer modeling. Comparison of the fingerprint for a particular candidate compound, whose three-dimensional structure is to be compared with an idealized description of an appropriate ligand (the pharmacophore) , to fingerprints of compounds having related activities provides substantial additional empirical information which can permit construction of more accurate three-dimensional representations of peptides or other macromolecules subject to conformational variation. The techniques of the invention also permit the reduction of large libraries of compounds to smaller sets that will, nevertheless, contain the compounds most likely to have a desired biological activity. The reduced size of the library permits more sophisticated tools to be applied to prediction of the affinity of the compounds in the reduced library for a target. Thus, because the size of a library is reduced, extensive conformational analysis of the ligand in the active site as well as of conformational changes of the active site in the presence of the ligand can be studied for the library members. This also permits a more accurate analysis of the electrostatic interactions between the ligand and the binding site, including solvation effects which are related to desolvation of the binding cavity and the ligand when these interact. The reduced library enabled by the invention is considerably smaller than that generally used in three- dimensional databases, allowing proportionately more computational effort to be expended on each compound. The most general application is simply to provide the maximum functional diversity for a given size of chemical library; this chemical library provides a core set for screening, a core set for computer screening, training sets, generally, and chromatographic ligands. This application is especially useful applied to a combinatorial library, in which large numbers of quite similar compounds are typically found. The utility of the pattern-comparison approach has been successfully shown to identify nonsteroidal antiinflammatory drugs (NSAIDs) . Many NSAIDs have been selected based on their ability to inhibit cyclooxygenase (COX, also known as prostaglandin synthase) which catalyzes the first step in the synthesis of prostaglandins, as well as on their activity in animal models. A second cyclooxygenase, COX-II has recently been discovered which is an isoenzyme of the originally known COX-I. COX-II is largely restricted to cells of the immune system and is believed more important than COX-I in inflammation. As set forth in more detail in Example 2, fingerprints were obtained for several hundred compounds, including two NSAIDs using the protein panels of the invention (containing 8-10 proteins) . Examination of the fingerprints of these two compounds showed a common feature which proved to be shared with several additional known NSAIDs for which profiles or fingerprints were subsequently obtained. The fingerprints for the several hundred compounds already tested were then searched for the presence or absence of this feature. Twelve compounds were found and these were tested for their ability to inhibit COX-I. Two compounds showed moderate and one measurable but low ability in this respect, although no NSAID activity had previously been reported for these compounds. The panel of proteins was then optimized as generally described above and used to evaluate a group of structurally diverse compounds containing seven known COX inhibitors and six inhibitors of other targets. The fingerprints obtained permitted completely accurate prediction of whether or not the compound was a COX inhibitor, although the proteins in the panel did not represent any proteins which were related to COX either by homology or by enzyme activity. Fingerprint Databases The reference panels of the invention, and reference panels generally, can be used to generate a fingerprint database which contains the fingerprints of a library of compounds in physically stored form to permit their retrieval. This form may be either a "paper" database, or is preferably in computer readable form. The database will contain the fingerprints of generally over 1,000 compounds with respect to a panel of proteins or other members wherein the number of panel members is less than three times the number of principal components represented in the panel. The compounds will represent a range of binding affinities for the panel members which is greater than three logs represented as IC 50 s. In the selected database, more than 95% of the compounds will provide fingerprints which are visible — i.e., are greater than the noise distance from the origin and will have mean separation from their nearest neighbors of more than three times the noise distance. These databases are useful in a variety of contexts. By applying multivariate statistical methods, equally diverse subsets can be obtained so that it can be verified that a subset selected from the database is of equal interest to the diversity represented by an alternative subset obtained by another method. Multivariate statistics can also be used to select a subset of maximal diversity for a defined size of the library,* for example, if the defined size is five times the number of members of the reference panel, it can be used as a training set, as described above. The database can also be used as a source for a diverse set of chromatography ligands. The following examples are intended to illustrate but not to limit the invention. Example 1 Factors Determining Minimum Requirements for Surrogate Construction In order to construct a surrogate, both the reference panel and the training set must be adequate. In order to obtain successful candidate compounds for a desired property, the library must be adequate as well. Confirmation that the reference panel contains an adequate number of properly chosen proteins can be accomplished by obtaining an X-Y plot of the distance between points in n-dimensional space (X axis) versus the frequency of this observed distance (Y axis) (distance distribution) . It will be recalled that each point in n-dimensional space represents the profile obtained for a single compound from the compound library with respect to a reference panel of n members. The height, shape, and maximum span of this distance distribution provides information as to the adequacy of the panel and the library. Ideally, a Poisson distribution should be obtained where the maximum of the distribution is at a high value of the distance between pairs. Figures 7a, 7b and 7c represent the distance distributions for the same set of compounds with respect to reference panels containing 5, 7 and 10 proteins, respectively. It is seen that when only five proteins are used in the panel, the shape of the distribution is somewhat irregular and the most frequent distance between points is relatively low. However, when the number of proteins in the panel is increased, a more regularly shaped Poisson distribution emerges with a larger distance between points at its maximum. The number of members in the panel is adequate when further addition of members fails to improve the position and shape of this distribution. Conversely, Figures 8a-8c reflect progress toward achieving an ideal distribution by a simple increase in the number of randomly chosen compounds in the compound library. The plot of pair-wise distances among compounds in a chemical library should provide a random distribution of distances if the collection of compounds is complete. If there are discontinuities, the collection is incomplete. In addition, large values of the maximum distance between members of a pair indicate more diversity in a set of compounds. This is illustrated in Figures 8a-8c. Figure 8a shows the frequency vs. distance plot for points representing fingerprints determined against a set of ten reference proteins for 50 compounds selected at random. The data do not result in a Poisson distribution and the maximum span of the distance is slightly over eight units. Figure 8b shows similar results when the fingerprints of 100 compounds are included; the distribution has become more regular and the maximum span has increased to approximately 12 units. When fingerprints for 1000 compounds were obtained and compared, the maximum separation between the points in n-dimensional space reaches 15 units and the distribution assumes the typical Poisson shape (Figure 8c) . Similar comparisons can be used to evaluate the adequacy of smaller numbers of putatively more representative compounds, for example, to evaluate the adequacy of combinatorial libraries comprised entirely of peptides. Figure 9b shows the distance distribution for 50 commercially available drugs. Comparing this distribution with that shown in Figure 9a (same as Figure 8a) for 50 structurally diverse random compounds reveals that the distributions are quite similar. However, when these distributions are compared to that obtained for a library of peptides ranging from dipeptides to 32-mers, as shown in Figure 9c, the portion of the space spanned is more than a unit smaller. This leads to the conclusion that peptide libraries per se may be inadequate to represent all of chemical space. The character of the distance distribution can also be used as a measure of the diversity of a particular set of candidate compounds, for example substances available as chromatographic ligands. Using the distance distribution as a criterion, a minimum number of ligands may be supplied to offer the widest possible spectrum of separation efficiency. In other words, such distance distributions can be used to verify the maximal diversity of panels of chromatographic ligands constructed as described in U.S. Patent 4,963,263, incorporated herein by reference, or to select nonpolymeric compounds to serve as diverse chromatographic ligands. Example 2 Discovery of Additional NSAIDs A data base of fingerprinted compounds which included fenoprofen, flufenamic acid, ibuprofen, endoprofen, ketoprofen, mefenamic acid, naproxen, piroxicam, and sulindac was prepared. A panel of proteins was prepared which were commercially obtained or expressed recombinantly in E. coli and purified. All of the proteins were enzymes in the initial panel and an IC 50 was determined in an enzymatic assay. A revised panel included other proteins and binding could be determined by fluorescence polarization. None of these proteins had any homology to the target of the NSAIDS, cyclooxygenase. Verification of predicted COX inhibitors was done by assessing COX activity in the presence and absence of the fingerprinted compound. Both COX-I from ram and COX-II from sheep were tested. The assays were conducted by incubating the enzyme at 37°C in 0.1 mM arachidonic acid contained in 0.1 M TRIS, pH 8.0 with 20 mM phenol. The reaction was stirred vigorously to maintain a significant concentration of dissolved oxygen for 3 minutes. The reaction was then stopped by addition of 5 mM citric acid and samples were diluted. The PGE 2 concentration was measured by EIA using a standard kit from Caymen Chemicals. The first several hundred compounds tested against the initial panel of 10 enzymes contained ibuprofen and indomethacin, two known NSAIDs. The 10 enzymes in this panel are those shown herein in Figure 10 and listed in Example 3 below. Both shared common features in their fingerprints that were tentatively considered to be diagnostic of an NSAID. In evaluating the remaining compound fingerprints, 12 additional compounds were selected which shared this feature. These were evaluated for their ability to inhibit COX-I and COX-II. Two COX-I inhibitors were found with moderate affinity and one COX-I inhibitor of low affinity. Thus, nine compounds which had been tentatively identified did not inhibit these enzymes, but two new leads were found without screening the whole library against COX. These novel leads are significantly different in structure from the known NSAIDS: ibuprofen and indomethacin. The reference panel was then revised to include enzymes which enriched the panel by expanding the range of chemicals that can be fingerprinted, and by increasing the average and maximal distances between fingerprints. The proteins in the panel are those listed on Figure 13 and are listed in Example 4 below. The compounds were refingerprinted. Of the nine compounds originally selected that then did not inhibit COX, only two remained putatively similar to the NSAID profile against this revised panel. The revised panel was used to fingerprint a group of 13 unidentified compounds and the fingerprints were compared to the NSAID consensus fingerprint obtained from ibuprofen and indomethacin. When compared, seven of the compounds exhibited features which predicted they would inhibit COX and six led to the prediction that they would not. The fingerprints accurately identified flosulide, phenylbutazone, pirprofen, prinomid, oxindanac, oxindanac analog, and diclofenac as inhibitors of COX and the compounds chlordiazepoxide, maprotiline, imipramine, metoprolol, and pentopril as noninhibitors. Among the predicted noninhibitors was also included a diclofenac prodrug which itself does not inhibit this enzyme. Example 3 Construction of a Surrogate In this example, the reference panel members whose profiles will be obtained with respect to a training set of compounds were isoenzymes of glutathione-S-transferase (GST) . The reference panel containing ten such isoenzymes is shown at the top of Figure 10. The target in this example was glutathione reductase (GRd) shown at the right. The first 20 80 test compounds of Figure 10 not used in the fitting procedure and Figure lib shows the residuals (experimental-predicted) from Figure 11a. The mathematical form for the linear regression is: n log ( IC 50 ) , tT = ∑ C Rj log ( IC 50 ) Rj . ( 1 ) j =l As shown in this formula, the IC 50 of compound i is measured against target T or reference protein R- weighted by a fitted coefficient C Rj . The successful correlation obtained above is surprising since the GRd derived from yeast is a NADPH dependent protein which has a different enzymatic function from GST. These enzymes share no sequence homology, and comparison of the crystal structures of GST and GRd reveals no tertiary structural similarities. The common use of glutathione does not appear to contribute to the correlation since the six peptide variants of glutathione which do bind various GSTs do not bind particularly well to GRd. Example 4 Improved Reference Panel/Compound Library Combinations The general procedure set forth in Example 3 was followed but using a different reference panel and expanded compound library. An initial set of eight proteins was chosen by preliminary screening of about 100 proteins generally expected to display a broad cross-reactivity against small organic molecules. The eight panel members were chosen based on enriching the panel of GSTs used in Example 3, as described above. Four of the final panel members were glutathione-S- transferase (GST) isoenzymes: human Al, rat R8, housefly HF2 and schistosome SI. The remaining panel members were D-amino acid oxidase (DAO) from porcine kidney (EC1.4.3.3) ; butyryl cholinesterase (BCh) from horse serum (EC3.1.1.8); papain (Pap) (EC3.4.22.2) and snake venom phosphodiesterase I (PDE) from Crotalus adamantaeus (EC3.1.4.1) . Cross-reactivity profiles were obtained with respect to this panel of eight proteins for a representative sample of 122 diverse compounds listed, along with their identification codes in Figure 12. For convenience, in determining the fingerprints, the binding of each compound to each protein was quantified as the concentration needed to inhibit 50% of the protein's activity (IC 50 ) . The IC 50 values ranged over more than four log units from 1 mM to less than 0.05 μM. A subset of 12 of the 122 compounds initially tested was chosen based on a high selectivity of these compounds for one or another of the proteins in the reference panel. This initial training set of 12 compounds was assayed for inhibitory activity with respect to two target enzymes: glutathione reductase (GRd) and aldehyde dehydrogenase (AdDH) . These two proteins are not related to each other and are not related by amino acid homology or activity to any of the reference proteins in the panel. The 12 compounds selected for the training set are the first 12 compounds for which results are shown in Figure 13. A surrogate was obtained based on this training set by applying a linear regression to the data to obtain the coefficients in Equation (1) above. This resulted in the following regression equations for this iteration: For glutathione reductase: 0.11 BCh + 0.19 HF2 + 1.79; For aldehyde dehydrogenase: 0.55 PDE + 1.35. The resulting surrogate was used to choose (for each target) a second set of 10 compounds (from the remaining 110) that were expected to be more representative of the range of potencies for the targets. These 10 compounds (marked by vertical bars in Figure 13, and different in most instances for each target) were then tested directly against the target compounds and the data obtained from these tests were used to supplement the results from the first 12 compounds, providing a total training set of 22 compounds for each target. Linear regression applied to this newly defined training set yielded the following forms of Equation (1) for the two targets: For glutathione reductase: 0.21 BCh + 0.72 HF2 + 0.24S1 - 0.05; © © in H υ A. oo ro 00 O fashion that Example 4 used a panel enriched with respect to that of Example 3. Surrogates were constructed against the additional targets: estrogen receptor, glycerol kinase, schistosome GST, nucleoside 5' -diphosphate kinase, human Factor Xa, trypsin and glyoxalase I. In each instance, a diverse set of 15-50 compounds, drawn from a database catalog of over 1,000 compounds, was used for the fitting. For each determination, the panel included at least the following enzymes: GST Al-l; acid 0.-1 glycoprotein; GST Pl-1; human serum albumin; papain; GST Rat 12:12(0); GST Housefly 3; butyryl cholinesterase; GST Rat 8:8; trypsin; and alcohol dehydrogenase. Of course, trypsin was not included in the panel for which trypsin was the counterpart target. In some instances, plasmin was substituted for GST Rat 8:8 and/or antitrypsin was substituted for alcohol dehydrogenase. These surrogates were correlated with experimentally determined binding as shown in Figure 15. Correlations generally showed a good match between the surrogate and the actual target. In each case, a different linear combination of the reference proteins provided the best fit. In all cases, there is no sequence homology between target and fitting proteins.

Claims: Claims 1. A method to identify a candidate reactive with a target, which method comprises: (a) providing a formula that represents a combination of the reactivity profiles with respect to a first set of candidates of at least two members of a reference panel, which formula calculates a predicted profile that best matches the reactivity profile of the target with respect to said first set of candidates; (b) testing the reactivity of said at least two members of the panel with respect to a candidate; and (c) calculating a predicted reactivity with respect to the target for said candidate by applying said formula to the reactivities determined in step (b) to estimate the reactivity of the candidate with respect to the target . 2. The method of claim 1 wherein the reactivities in step (b) are tested with respect to each member of a library of candidates, and/or wherein the combination in (a) is a linear combination. 3. The method of claim 1 wherein application of step (c) results in compounds which are estimated to react well and compounds that are estimated to react poorly with a target, and at least some of the compounds which are estimated to react well and at least some of the compounds which are estimated to react poorly with the target are added to the first set of candidates, to generate a second set of candidates and step (a) is repeated with said second set of candidates to obtain an improved formula, and which further includes (d) testing the reactivity of said at least two members of said reference panel with respect to a candidate; and (e) calculating a predicted reactivity with respect to the target for said candidate by applying said formula to the reactivities determined in step (d) to estimate the reactivity of the candidate with respect to the target. 4. A method to prepare a substance that reacts with a desired target which method comprises: (a) providing a formula that represents a combination of the reactivity profiles with respect to a first set of candidates of at least two members of a reference panel which formula calculates a predicted profile that best matches the reactivity profile of the target with respect to said first set of candidates; (b) testing the reactivity of said at least two members of the panel with respect to a candidate; (c) calculating a predicted reactivity with respect to the target for said candidate by applying said formula to the reactivities determined in step (b) to estimate the reactivity of the candidate with respect to the target; (d) identifying a substance as being a candidate predicted to react with the target, and (e) assembling the identified substance from starting materials appropriate to said substance. 5. A method to construct a reference panel for predicting reactivity of a candidate for a target which method comprises: arbitrarily identifying an initial set of panel members; obtaining profiles of reactivity for an initial set of arbitrarily chosen compounds with respect to said initial set of panel members; comparing the profiles obtained; discarding compounds and panel members which result in redundant profiles; substituting additional provisional panel members and compounds for the panel members and compounds discarded to obtain a second set of panel members and a second set of compounds; obtaining profiles for the second set of compounds with respect to said second set of panel members; again comparing the profiles obtained and discarding compounds and panel members that result in redundant profiles; and repeating the foregoing steps until a panel which covers at least 90% of chemical space is obtained; and/or which provides at least 5 principal components with respect to the range of compounds marketed as small organic molecules; and/or wherein for said panel, the average of the differences between a profile for any first compound from that of any second compound is at least three times the differences observed for repeated determinations of the profile of said first compound. 6. A reference panel of members comprising proteins, wherein at least one of the members of the panel is a protein other than an immunoglobulin (Ig) or fragment thereof; and wherein said panel provides members binding in a multiplicity of differing degrees with respect to a population of compounds; and wherein the presence of said non-Ig protein enriches the panel. 7. The panel of claim 6 wherein said panel covers 90% of chemical space; and/or wherein said panel provides at least 5 principal components with respect to the range of compounds marketed as small organic molecules; and/or wherein for said panel, the average of the differences between a profile for any first compound from that of any second compound is at least three times the differences observed for repeated determinations of the profile of said first compound. 8. The panel of claim 6 wherein said panel comprises at least 2 enzymes or at least 2 lectins or at least 2 T cell receptors or at least 2 olfactory receptors, and/or wherein said panel comprises at least 10% non-Ig proteins. 9. A method to characterize a single analyte, which method comprises: contacting said analyte with each member of a panel of claim 8, wherein said compound is the analyte; detecting the degree of reactivity of said analyte to each of said members; recording said degree of reactivity of said analyte to each of said members; and arranging said recorded degrees of reactivity so as to provide a characteristic profile of said analyte. 10. A method to identify a candidate which candidate will be effective in reacting with a target, wherein said target has a known ligand with which it reacts, which method comprises: contacting said candidate with each member of the panel of claim 8 wherein said compound is the candidate; detecting the degree of reactivity of said candidate to each of said members; recording each said degree of reactivity of said candidate to each of said members; arranging said recorded degrees of reactivity so as to provide a characteristic profile of said candidate; comparing said profile to a profile analogously obtained of said ligand with respect to said multiplicity of members; wherein similarity of the profile of said candidate to the profile of said ligand indicates the probability that the candidate will react with said target. 11. The method of claim 10 wherein said comparing includes the steps of: determining a point obtained by plotting, in n-dimensional space, the profile of reactivity of the candidate for each member of the panel, wherein each n dimension represents a different member of the panel and the reactivity of the candidate with said each member is plotted in said each n dimension; and comparing the position of said point to the point in said n-dimensional space determined for the profile representing the reactivity of the known ligand for each member of the panel wherein proximity of the points indicates the degree of binding of the candidate to the target. 12. A method to prepare a substance that reacts with a desired target wherein the target has a known ligand with which it reacts which method comprises: contacting a candidate with each member of the panel of claim 8 wherein said compound is said candidate; detecting the degree of reactivity of said candidate to each of said members; recording each said degree of reactivity of said candidate to each of said members; arranging said recorded degrees of reactivity so as to provide a characteristic profile of said candidate; comparing said profile to a profile analogously obtained of said ligand with respect to said multiplicity of members; wherein similarity of the profile of said candidate to the profile of said ligand indicates the probability that the candidate will react with said target; identifying a substance as a candidate having a high probability of reacting with target that has a profile similar to that of the ligand, and assembling the identified substance from the starting materials appropriate to said substance. 13. A method to select from a multiplicity of candidates a candidate that reacts specifically with a known target, which method comprises: providing a profile of reactivity of said target against a maximally diverse set; preparing a profile of the reactivity of the candidate with respect to the panel of claim 8 which panel is the inverse image of the maximally diverse set; comparing the maximally diverse set profile of the target with the inverse image panel profile of the candidate; and wherein similarity of the inverse image panel profile with maximally diverse set profile indicates the probability that the candidate will bind to the target . 14. The method of claim 13 wherein said comparing includes the steps of : determining a point obtained by plotting, in n-dimensional space, the profile of reactivity of the candidate for each member of the inverse image panel, wherein each n dimension represents a different member of the inverse image panel and the reactivity of the candidate with said each member is plotted in said each n dimension; and comparing the position of said point to the point in said n-dimensional space determined for the profile representing the reactivity of the target for each member of the maximally diverse set wherein proximity of the points indicates the degree of binding of the candidate to the target . 15. A method to prepare a substance that reacts with a desired target which method comprises: providing a profile of reactivity of said target against a maximally diverse set; preparing a profile of the reactivity of the candidate tc the panel of claim 8 which is the inverse image of the maximally diverse set; comparing the maximally diverse set profile of the target with the inverse image panel profile of the candidate; and wherein similarity of the inverse image panel profile with the maximally diverse set profile indicates the probability that the candidate will bind to the target; identifying a substance as being a candidate having a high probability of reacting with target that has a profile similar to that of the receptor, and assembling the identified substance from the starting materials appropriate to said substance. 16. A method to select from a multiplicity of candidates a candidate that reacts specifically with a known target, which method comprises: providing a profile of reactivity of said candidate against a maximally diverse set; preparing a profile of the reactivity of the target with respect to the panel of claim 8, which panel is the inverse image of said maximally diverse set; comparing the maximally diverse panel profile of the candidate with the inverse image panel profile of the target; and wherein similarity of the inverse image panel profile with maximally diverse set profile indicates the probability that the candidate will bind to the target. 17. The method of claim 16 wherein said comparing includes the steps of: determining a point obtained by plotting, in n-dimensional space, the profile of reactivity of the candidate for each member of the maximally diverse set, wherein each n dimension represents a different member of the panel and the reactivity of the candidate with said each member is plotted in said each n dimension; and comparing the position of said point to the point in said n-dimensional space determined for the profile representing the reactivity of the target for each member of the inverse image panel wherein proximity of the points indicates the degree of binding of the candidate to the target. 18. A method to prepare a substance that reacts with a desired target which method comprises: providing a profile of reactivity of a candidate against a maximally diverse set; preparing a profile of the reactivity of the target to panel of claim 8 which is the inverse image of the maximally diverse set; comparing the maximally diverse set profile of the candidate with the inverse image panel profile of the target; and wherein similarity of the inverse image panel profile with maximally diverse set profile indicates the probability that the candidate will bind to the target; identifying a substance as being a candidate having a high probability of reacting with target that has a profile similar to that of the target, and assembling the identified substance from the starting materials appropriate to said substance. 19. A fingerprint for a compound which is a characteristic profile obtained by the method of claim 9. 20. A database comprising multiple fingerprints of claim 19. 21. A subset of the database of claim 20 selected to provide maximal diversity. 22. The subset of claim 21 for use as a core screening set, or for use as a source of chromatographic ligands. 23. A database which comprises the physical embodiment of fingerprints for more than 1000 compounds against a reference panel, said compounds representing a range of binding with the panel members of at least three log units, and wherein at least 95% of the compounds generate visible fingerprints. 24. A subset of the database of claim 23 selected to provide maximal diversity. 25. The subset of claim 24 for use as a core screening set, or for use as a source of chromatographic ligands.

Inventors: KAUVAR LAWRENCE M, VILLAR HUGO O, KAUVAR, LAWRENCE, M., VILLAR, HUGO, O.
Applicants: TERRAPIN TECH INC, TERRAPIN TECHNOLOGIES, INC.
Application year: 1995
Main classification: G01N33
IPC codes: G01N33