To be of any value, a predicted model of an antibody combining site should have an accuracy approaching that of antibody structures determined by x-ray crystallography (1.6-2.7 A). A number of modeling protocols have been proposed, which fall into two main categories-those that adopt a knowledge-based approach and those that attempt to construct the hypervariable loop regions of the antibody ab initio. Here we present a combined algorithm requiring no arbitrary decisions on the part of the user, which has been successfully applied to the modeling of the individual loops in two systems: the anti-lysozyme antibody HyHel-5, the crystal structure of which is as a complex with lysozyme [Sheriff, S., Silverton, E. W., Padlan, E. A., Cohen, G. H., Smith-GM, S. The wide range of specificities exhibited by antibodies is a function of the sequence and length variability of six hypervariable loops or complementarity-determining regions (CDRs) (1), which form the antigen combining site. These six CDRs supported on a highly conserved framework region constitute the variable region of the antigen binding fragment (Fab). A knowledge of antibody structure is essential for intelligent design ofantibody enzymes (2), tailoring of affinity (3), and CDR replacement strategies (4). However, sequence information vastly exceeds structural information from x-ray crystallography and, until crystallographic structure determination becomes no less routine than sequencing, modeling of structures is necessary. Since the framework region is conserved, it has proved relatively easy to model, whereas the CDRs, by their very nature (5), present a more challenging problem since accuracy in their modeling is of paramount importance.The approaches taken to modeling the antibody combining site, so far, fall into two groups: knowledge based and ab initio. Knowledge-based approaches have been used to model a number of antibodies, including J539 (6), GLOOP1-5 (5), HyHel-10 (7), and D1.3 (8). Although the methods differ in their detail, the common feature of all the approaches has been to examine only the known antibody crystal structures and select CDRs from these on the basis of length and/or sequence. Although most methods use simple sequence homology to select model loop conformations, Chothia and Lesk (9) have obtained better results by selecting conformations on the basis of the conservation of "key" residues that affect loop packing or conformation. However, the chief problem with any such method is the limited size of the knowledge base: while this has been extended to include all protein loops in the broader protein modeling field (10), a general data base has not previously been used in antibody modeling.The second approach has been to use ab initio conformational search algorithms to saturate the conformational space available to a loop and select an appropriate structure on the basis of its energy, calculated by using an empirical energy function (11). Whereas this overcomes the limited size of the knowledge base, it fails to make us...