An analysis of incorrectly folded protein models. Implications for structure predictions.

Proteins with homologous amino acid sequences have similar folds and it has been assumed that an unknown three-dimensional structure can be obtained from a known homologous structure by substituting new side-chains into the polypeptide chain backbone, followed by relatively small adjustment of the model. To examine this approach of structure prediction and, more generally, to isolate the characteristics of native proteins, we constructed two incorrectly folded protein models. Sea-worm hemerythrin and the variable domain of mouse immunoglobulin K-chain, two proteins with no sequence homology, were chosen for study; the former is composed of a bundle of four alpha-helices and the latter consists of two 4-stranded beta-sheets. Using an automatic computer procedure, hemerythrin side-chains were substituted into the immunoglobulin domain and vice versa. The structures were energy-minimized with the program CHARMM and the resulting structures compared with the correctly folded forms. It was found that the incorrect side-chains can be incorporated readily into both types of structures (alpha-helices, beta-sheets) with only small structural adjustments. After constrained energy-minimization, which led to an average atomic co-ordinate shift of no more than 0.7 to 0.9 A, the incorrectly folded models arrived at potential energy values comparable to those of the correct structures. Detailed analysis of the energy results shows that the incorrect structures have less stabilizing electrostatic, van der Waals' and hydrogen-bonding interactions. The difference is particularly pronounced when the electrostatic and van der Waals' energy terms are calculated by modified equations that include an approximate representation of solvent effects. The incorrectly folded structures also have a significantly larger solvent-accessible surface and a greater fraction of non-polar side-chain atoms exposed to solvent. Examination of their interior shows that the packing of side-chains at the secondary structure interfaces, although corresponding to sterically allowed conformations, deviates from the characteristics found in normal proteins. The analysis of incorrectly folded structures has made it clear that the absence of bad non-bonded contacts, though necessary, is not sufficient to demonstrate the validity of model-built structures and that modeling of homologous structures has to be accompanied by a thorough quantitative evaluation of the results. Further, certain features that characterize native proteins are made evident by their absence in misfolded models.