On the computation of the tertiary structure of globular proteins--IV. Use of secondary structure information.

Abstract

The distance geometry approach for computing the tertiary structure of globular proteins emphasized in this series of papers (Goelet al., J. theor. Biol. 99, 705–757, 1982) is developed further. This development includes incorporation of some secondary structure information—the location of alpha helices in the primary sequence—in the algorithm to compute the tertiary structure of alpha helical globular proteins. An algorithm is developed which estimates the interresidue distances between chain-proximate helices. These distances, in conjunction with the global statistical average distances obtainable from a database of real proteins and determined by the primary sequence of the protein under study, are used to determine the tertiary structure. Five proteins, parvalbumin, hemerythrin, human hemoglobin, lamprey hemoglobin, and sperm whale myoglobin, are investigated. The root mean square (RMS) errors between the calculated structures and those determined by X-ray diffraction range from 4.78 to 7.56 A. These RMSs are 0.21–2.76 A lower than those estimated without the secondary structure information. Contact maps and three-dimensional backbone representations also show considerable improvements with the introduction of secondary structure information.