Theoretical determination of sterically allowed conformations of a polypeptide chain by a computer method

Abstract

Abstract A mathematical method has been developed, together with a digital computer program for its evaluation, for the determination of possible conformations of a polypeptide chain for any given amino acid sequence. Allowed conformations were required to satisfy predetermined criteria regarding the locations of potential energy minima for rotation around single bonds, steric restrictions arising from van der Waals contacts between atoms, and the formation of a closed polypeptide loop by means of a peptide linkage or a disulfide bridge. Requirements for specific noncovalent interactions can also be considered. Various assumptions regarding the selection of the minima for rotational potentials and of the van der Waals contact distances are discussed. According to the assumptions adopted in the paper, 21 conformations are allowed for a glycine residue which is part of a dipeptide. These conformations are assumed to lie at minima of the rotational potential function, and satisfy simultaneously the steric requirements. Dipeptides containing alanine or other residues are restricted to 7 such conformations. As an example of the applicability of the method, the computation was performed for the smallest loop containing a disulfide bridge in the ribonuclease molecule, i.e., the cyclic octapeptide formed by residues 65 to 72. Only 15 allowed conformations were found for the loop as a whole. Parameters characterizing these allowed conformations are tabulated. Approximate calculations indicate the importance of steric restrictions and of the presence of side chains in reducing the conformational freedom of this peptide. Further potentialities of the method are discussed.