Development of an extended simulated annealing method: application to the modeling of complementary determining regions of immunoglobulins.

An extended simulated annealing process (ESAP) has been developed in order to obtain an ensemble of conformations of a peptide segment from a protein fluctuating at a given temperature. The annealing process was performed with a fast Monte Carlo method using the scaled collective variables developed by Noguti and Go. The system was divided into two parts: one consists of one or more peptide segments and is flexible around the main-chain and side-chain torsional angles; the other represents the rest of the molecule and was maintained fixed at the atomic positions determined by x-ray experiments. The target function included the nonbonding atomic interactions and a distance function to anchor the N and C terminal ends of each segment to the fixed part. Three systems of complementary determining regions (CDR) of antibodies were tested and compared to x-ray data: L2 loop (7 residues) of the light chain of lambda-type Bence-Jones protein, H1 and the H2 loops (14 residues) of McPC603, and H1 and H2 loops (12 residues) of HyHEL-5. Each state of CDR conformations was characterized at room temperature by the average of their coordinates (average conformation) and the internal energy. With a limited number of annealing processes (10), starting from the extended conformation, we have obtained states with conformations close to the observed x-ray structures, from 1.1 to 1.7 A root mean square deviation (rmsd) of main-chain atoms depending on the system. These states were identical or within 0.25 A rmsd of those of lowest internal energy. For unknown CDR structures the criteria of lowest internal energies from ESAP can be used to predict hypervariable loop structures in antibodies with an accuracy comparable to other methods.     

Structural studies of the streptavidin binding loop.

The streptavidin-biotin complex provides the basis for many important biotechnological applications and is an interesting model system for studying high-affinity protein-ligand interactions. We report here crystallographic studies elucidating the conformation of the flexible binding loop of streptavidin (residues 45 to 52) in the unbound and bound forms. The crystal structures of unbound streptavidin have been determined in two monoclinic crystal forms. The binding loop generally adopts an open conformation in the unbound species. In one subunit of one crystal form, the flexible loop adopts the closed conformation and an analysis of packing interactions suggests that protein-protein contacts stabilize the closed loop conformation. In the other crystal form all loops adopt an open conformation. Co-crystallization of streptavidin and biotin resulted in two additional, different crystal forms, with ligand bound in all four binding sites of the first crystal form and biotin bound in only two subunits in a second. The major change associated with binding of biotin is the closure of the surface loop incorporating residues 45 to 52. Residues 49 to 52 display a 3(10) helical conformation in unbound subunits of our structures as opposed to the disordered loops observed in other structure determinations of streptavidin. In addition, the open conformation is stabilized by a beta-sheet hydrogen bond between residues 45 and 52, which cannot occur in the closed conformation. The 3(10) helix is observed in nearly all unbound subunits of both the co-crystallized and ligand-free structures. An analysis of the temperature factors of the binding loop regions suggests that the mobility of the closed loops in the complexed structures is lower than in the open loops of the ligand-free structures. The two biotin bound subunits in the tetramer found in the MONO-b1 crystal form are those that contribute Trp 120 across their respective binding pockets, suggesting a structural link between these binding sites in the tetramer. However, there are no obvious signatures of binding site communication observed upon ligand binding, such as quaternary structure changes or shifts in the region of Trp 120. These studies demonstrate that while crystallographic packing interactions can stabilize both the open and closed forms of the flexible loop, in their absence the loop is open in the unbound state and closed in the presence of biotin. If present in solution, the helical structure in the open loop conformation could moderate the entropic penalty associated with biotin binding by contributing an order-to-disorder component to the loop closure.     

Alanine-scanning mutagenesis of the beta-sheet region of phage T4 lysozyme suggests that tertiary context has a dominant effect on beta-sheet formation

In general, alpha-helical conformations in proteins depend in large part on the amino acid residues within the helix and their proximal interactions. For example, an alanine residue has a high propensity to adopt an alpha-helical conformation, whereas that of a glycine residue is low. The sequence preferences for beta-sheet formation are less obvious. To identify the factors that influence beta-sheet conformation, a series of scanning polyalanine mutations were made within the strands and associated turns of the beta-sheet region in T4 lysozyme. For each construct the stability of the folded protein was reduced substantially, consistent with removal of native packing interactions. However, the crystal structures showed that each of the mutants retained the beta-sheet conformation. These results suggest that the structure of the beta-sheet region of T4 lysozyme is maintained to a substantial extent by tertiary interactions with the surrounding parts of the protein. Such tertiary interactions may be important in determining the structures of beta-sheets in general.     

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Determination of structurally conservative amino acids of the HIV-1 protein gp120 V3 loop as promising targets for drug design by protein engineering approaches

Based on the published NMR spectroscopy data, three-dimensional structures of the HIV-1 gp120 protein V3 loop were obtained by computer modeling in the viral strains HIV-Haiti and HIV-MN. In both cases, the secondary structure elements and conformations of irregular stretches were determined for the fragment representing the principal antigenic determinant of the virus, as well as determinants of the cellular tropism and syncytium formation. Notwithstanding the high variability of the amino acid sequence of gp120 protein, more than 50% of the V3 loop residues retained their conformations in the different HIV-1 virions. The combined analysis of the findings and the literature data on the biological activity of the individual residues of the HIV-1 V3 loop resulted in identification of its structurally conservative amino acids, which seem to be promising targets for antiviral drug design by protein engineering approaches.     

Three-dimensional structure of apotransketolase. Flexible loops at the active site enable cofactor binding.

The structure determination of apotransketolase and the comparison of its three-dimensional structure with that of the holoenzyme has revealed that no large conformational changes are associated with cofactor binding. Two loops at the active site are flexible in the apoenzyme which enables ThDP to reach its binding site. Binding of the cofactor induces defined conformations for these two loops at the active site. One of these loops is directly involving in binding of the cofactors, Ca2+ and ThDP. This loop acts like a flap which closes off the diphosphate binding site. After binding of the cofactor, residues of this loop form interactions to residues of loop 383-398 from the second subunit. These interactions stabilize the conformation of the two loops from a flexible to a 'closed' conformation.     

Protein loops on structurally similar scaffolds: database and conformational analysis

A general problem in comparative modeling and protein design is the conformational evaluation of loops with a certain sequence in specific environmental protein frameworks. Loops of different sequences and structures on similar scaffolds are common in the Protein Data Bank (PDB). In order to explore both structural and sequential diversity of them, a data base of loops connecting similar secondary structure fragments is constructed by searching the data base of families of structurally similar proteins and PDB. A total of 84 loop families having 2-13 residues are found among the well-determined structures of resolution better than 2.5 A. Eight alpha-alpha, 20 alpha-beta, 19 beta-alpha, and 37 beta-beta families are identified. Every family contains more than 5 loop motifs. In each family, no loops share same sequence and all the frameworks are well superimposed. Forty-three new loop classes are distinguished in the data base. The structural variability of loops in homologous proteins are examined and shown in 44 families. Motif families are characterized with geometric parameters and sequence patterns. The conformations of loops in each family are clustered into subfamilies using average linkage cluster analysis method. Information such as geometric properties, sequence profile, sequential and structural variability in loop, structural alignment parameters, sequence similarities, and clustering results are provided. Correlations between the conformation of loops and loop sequence, motif sequence, and global sequence of PDB chain are examined in order to find how loop structures depend on their sequences and how they are affected by the local and global environment. Strong correlations (R > 0.75) are only found in 24 families. The best R value is 0.98. The data base is available through the Internet.     

Predicting antibody hypervariable loop conformations. II: Minimization and molecular dynamics studies of MCPC603 from many randomly generated loop conformations

We describe a method for predicting the conformations of loops in proteins and its application to four of the complementarity determining regions [CDRs] in the crystallographically determined structure of MCPC603. The method is based on the generation of a large number of randomly generated conformations for the backbone of the loop being studied, followed by either minimization or molecular dynamics followed by minimization starting from these random structures. The details of the algorithm for the generation of the loops are presented in the first paper in this series (Shenkin et al. [submitted]). The results of minimization and molecular dynamics applied to these loops is presented here. For the two shortest CDRs studied (H1 and L2, which are five and seven amino acids long), minimizations and dynamics simulations which ignore interactions of the loop amino acids beyond the carbon beta replicate the conformation of the crystal structure closely. This suggests that these loops fold independently of sequence variation. For the third CDR (L3, which is nine amino acids), those portions of the CDR near its base which are hydrogen bonded to framework are well replicated by our procedures, but the top of the loop shows significant conformational variability. This variability persists when side chain interactions for the MCPC603 sequence are included. For a fourth CDR (H3, which is 11 amino acids long), new low-energy backbone conformations are found; however, only those which are close to the crystal are compatible with the sequence when side chain interactions are taken into account. Results from minimization and dynamics on single CDRs with all other CDRs removed are presented. These allow us to explore the extent to which individual CDR conformations are determined by interactions with framework only.     

Local propensities and statistical potentials of backbone dihedral angles in proteins

The following three issues concerning the backbone dihedral angles of protein structures are presented. (1) How do the dihedral angles of the 20 amino acids depend on the identity and conformation of their nearest residues? (2) To what extent are the native dihedral angles determined by local (dihedral) potentials? (3) How to build a knowledge-based potential for a residue's dihedral angles, considering the identity and conformation of its nearest residues? We find that the dihedral angle distribution for a residue can significantly depend on the identity and conformation of its adjacent residues. These correlations are in sharp contrast to the Flory isolated-pair hypothesis. Statistical potentials are built for all combinations of residue triplets and depend on the dihedral angles between consecutive residues. First, a low-resolution potential is obtained, which only differentiates between the main populated basins in the dihedral angle density plots. Minimization of the dihedral potential for 125 test proteins reveals that most native alpha-helical residues (89%) and a large fraction of native beta-sheet residues (47%) adopt conformations close to their native one. For native loop residues, the percentage is 48%. It is also found that this fraction is higher for residues away from the ends of alpha or beta secondary structure elements. In addition, a higher resolution potential is built as a function of dihedral angles by a smoothing procedure and continuous functions interpolations. Monte Carlo energy minimization with this potential results in a lower fraction for native beta-sheet residues. Nevertheless, because of the higher flexibility and entropy of beta structures, they could be preferred under the influence of non-local interactions. In general, most alpha-helices and many beta-sheets are strongly determined by the local potential, while the conformations in loops and near the end of beta-sheets are more influenced by non-local interactions.     

Browsing the SLoop database of structurally classified loops connecting elements of protein secondary structure

We describe a web server, which provides easy access to the SLoop database of loop conformations connecting elements of protein secondary structure. The loops are classified according to their length, the type of bounding secondary structures and the conformation of the mainchain. The current release of the database consists of over 8000 loops of up to 20 residues in length. A loop prediction method, which selects conformers on the basis of the sequence and the positions of the elements of secondary structure, is also implemented. These web pages are freely accessible over the internet at http://www-cryst.bioc.cam.ac.uk/ approximately sloop.     

Comparison of Molecular Dynamics Averaged Structures for Complexes of Normal and Oncogenic ras-p21 with SOS Nucleotide Exchange Protein, Containing Computed Conformations for Three Crystallographically Undefined Domains, Suggests a Potential Role of These Domains in ras Signaling

ras-p21 protein binds to the son-of-sevenless (SOS) guanine nucleotide-exchange promoter that allows it to exchange GDP for GTP. Previously, we performed molecular dynamics calculations on oncogenic (Val 12-) and wild-type ras-p21 bound to SOS. By superimposing the average structures of these two complexes, we identified four domains (residues 631-641, 676-691, 718-729, and 994-1004) in SOS that change conformation and were candidates for being effector domains. These calculations were performed in the absence of three crystallographically undefined loops (i.e., residues 591-596, 654-675, and 742-751). We have now modeled these loops into the SOS structure and have re-performed the dynamics calculations. We find that all three loop domains undergo large changes in conformation that involve mostly changes in their positioning and not their individual conformations. We have also identified another potential effector domain (i.e., residues 980-989). Overall, our current results suggest that SOS interactions with oncogenic ras-p21 may enhance ras-p21 mitogenic signaling through prolonging its activation by maintaining its binding to GTP and by allowing its effector domains to interact with intracellular targets.     

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.     

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.     

Modeling the anti-CEA antibody combining site by homology and conformational search.

A model for an antibody specific for the carcinoembryonic antigen (CEA) has been constructed using a method which combines the concept of canonical structures with conformational search. A conformational search technique is introduced which couples random generation of backbone loop conformations to a simulated annealing method for assigning side chain conformations. This technique was used both to verify conformations selected from the set of known canonical structures and to explore conformations available to the H3 loop in CEA ab initio. Canonical structures are not available for H3 due to its variability in length, sequence, and observed conformation in known antibody structures. Analysis of the results of conformational search resulted in three equally probable conformations for H3 loop in CEA. Force field energies, solvation free energies, exposure of charged residues and burial of hydrophobic residues, and packing of hydrophobic residues at the base of the loop were used as selection criteria. The existence of three equally plausible structures may reflect the high degree of flexibility expected for an exposed loop of this length. The nature of the combining site and features which could be important to interaction with antigen are discussed.     

Unbiased,scalable sampling of protein loop conformations from probabilistic priors

Background

Protein loops are flexible structures that are intimately tied to function, but understanding loop motion and generating loop conformation ensembles remain significant computational challenges. Discrete search techniques scale poorly to large loops, optimization and molecular dynamics techniques are prone to local minima, and inverse kinematics techniques can only incorporate structural preferences in adhoc fashion. This paper presents Sub-Loop Inverse Kinematics Monte Carlo (SLIKMC), a new Markov chain Monte Carlo algorithm for generating conformations of closed loops according to experimentally available, heterogeneous structural preferences.

Results

Our simulation experiments demonstrate that the method computes high-scoring conformations of large loops (> 10 residues) orders of magnitude faster than standard Monte Carlo and discrete search techniques. Two new developments contribute to the scalability of the new method. First, structural preferences are specified via a probabilistic graphical model (PGM) that links conformation variables, spatial variables (e.g., atom positions), constraints and prior information in a unified framework. The method uses a sparse PGM that exploits locality of interactions between atoms and residues. Second, a novel method for sampling sub-loops is developed to generate statistically unbiased samples of probability densities restricted by loop-closure constraints.

Conclusion

Numerical experiments confirm that SLIKMC generates conformation ensembles that are statistically consistent with specified structural preferences. Protein conformations with 100+ residues are sampled on standard PC hardware in seconds. Application to proteins involved in ion-binding demonstrate its potential as a tool for loop ensemble generation and missing structure completion.

     

Multiple conformations of amino acid residues in ribonuclease A

The highly refined 1.26 A structure (R = 0.15) of phosphate-free bovine pancreatic ribonuclease A was modeled with 13 residues having discrete multiple conformations of side chains. These residues are widely distributed over the protein surface, but only one of them, Lys 61, is involved in crystal packing interactions. The discrete conformers have no unusual torsion angles, and their interactions with the solvent and with other atoms of the protein are similar to those residues modeled with a single conformation. For three of the residues--Val 43, Asp 83, and Arg 85--two correlated conformations are found. The observed multiple conformations on the protein surfaces will be of significance in analyzing structure-function relationships and in performing protein engineering.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Common features of the conformations of antigen-binding loops in immunoglobulins and application to modeling loop conformations.

Using database screening techniques we have examined the relationship between antigen-binding loops in immunoglobulins, and regions of similar conformation in other protein families. The conformations of most antigen-binding loops are not unique to immunoglobulins. But in many cases, the geometrical relationship between the loop and the peptides flanking it differs between the immunoglobulins and other structures with the same loop. We assess model building by data base screening, compared with that based on canonical structures.