Calculation of protein conformation as an assembly of stable overlapping segments: application to bovine pancreatic trypsin inhibitor.

Abstract

Conformations of bovine pancreatic trypsin inhibitor were calculated by assuming that the final structure as well as properly chosen overlapping segments thereof are simultaneously in low-energy (not necessarily the lowest-energy) conformational states. Therefore, the whole chain can be built up from building blocks whose conformations are determined primarily by short-range interactions. Our earlier buildup procedure was modified by taking account of a statistical analysis of known amino acid sequences that indicates that there is nonrandom pairing of amino acid residues in short segments along the chain, and by carrying out energy minimization on only these segments and on the whole chain [without minimizing the energies of intermediate-size segments (20-30 residues long)]. Results of this statistical analysis were used to determine the variable sizes of the overlapping oligopeptide building blocks used in the calculations; these varied from tripeptides to octapeptides, depending on the amino acid sequence. Successive stages of approximations were used to combine the low-energy conformations of these building blocks in order to keep the number of variables in the computations to a manageable size. The calculations led to a limited number of conformations of the protein (only two different groups, with very similar structure within each group), most residues of which were in the same conformational state as in the native structure.