Homology modeling of divergent proteins

A method is presented for homology modeling of proteins bearing weak sequence identity to proteins of known tertiary structure. To accommodate non-identical amino acids in the core region, the backbone of the structurally conserved core of the model protein is allowed to deviate from that of the template protein. We have expanded FOLDER, a distance geometry-based homology modeling method, to allow for such displacements in the structurally conserved core. Models are built by rigidly constraining the interatomic distances within a structurally conserved segment and by allowing the interatomic distances between these segments to vary by a "divergence factor". We test this method by simulating models of the beta-barrel domain D1 of CD4 and a four-helix bundle protein cytochrome b562 using the crystal structures of Bence-Jones protein and cytochrome c' as templates, respectively. In both cases, previously published structure-based sequence alignments were used for simulating models. The root-mean-square (r.m.s.) deviation of the backbone atoms in the common core between the templates and models was found to be a function of the imposed divergence factor. Our results demonstrate that this r.m.s. deviation results from the relative displacements of structurally conserved segments to accommodate the amino acid replacements in the core of the model protein. To test the integrity of the simulated structures we compared them with their respective crystal structures. The r.m.s. deviation of the backbone atoms in the core regions of the simulated models and their respective crystal structures is approximately 1.4 A. The r.m.s. deviation for all the backbone atoms in the models, including those in the structurally variable regions, which are modeled de novo, is 2.4 A for CD4 and 3.2 A for cytochrome b562 when compared with their respective X-ray structures.