On the multiple-minima problem in the conformational analysis of polypeptides. V. Application of the self-consistent electrostatic field and the electrostatically driven Monte Carlo methods to bovine pancreatic trypsin inhibitor.
In connection with the accompanying paper to test various models for the hydration of polypeptides, we have explored a limited portion of the conformational energy hyperspace of the small protein bovine pancreatic trypsin inhibitor (BPTI) with the aid of two search methods developed in this laboratory. A series of low-energy conformations was obtained as a result of this study. These conformations constitute a set of local minima in the conformational energy space of the molecule as described by the ECEPP/2 (Empirical Conformational Energy Program for Peptides) potential energy function, without the inclusion of hydration. Five different initial conformations were used in this exploration: the first corresponds to an energy-refined structure based on the crystallographic coordinates (4PTI) provided by Deisenhofer and Steigemann and reported previously by Meirovitch and Scheraga. The remaining four initial conformations were obtained by using a Variable-Target-Function procedure, applied to the experimental Cartesian coordinates (5PTI) reported by Wlodawer et al. The self-consistent electrostatic field (SCEF) and the electrostatically driven Monte Carlo (EDMC) methods were used to search the conformational space. The SCEF and EDMC methodologies assume that a polypeptide or protein molecule is driven toward the native structure mainly by the action of the electrostatic interactions. Application of these methodologies led to a set of conformations (up to 50 kcal/mol lower than the starting ones) with ECEPP/2 energies lower than any of those that we had previously found. Application of both methods to the initial conformation generated from 4PTI led to a series of low-energy conformations exhibiting similar rms deviations with respect to the experimental data (4PTI) as did the starting conformation. However, statistical analysis of the runs that had started from the conformations generated by using the variable-target-function procedure (and applying the EDMC method) indicated that the rms deviations of the atomic positions of the new low-energy conformations tended to increase as the energy improved, when compared with the X-ray data from which the starting conformations had been generated. The structures with the lowest energies also had radii of gyration smaller than the experimentally observed one. These results indicated a need to include hydration in the potential function, and provided the conformations used in the accompanying paper to test various hydration models.
Ripoll
D R
DR
Baker Laboratory of Chemistry, Cornell University, Ithaca, New York 14853-1301.
Piela
L
L
Vásquez
M
M
Scheraga
H A
HA
D000818:Animals; D007611:Aprotinin / Q000737:chemistry*; D002417:Cattle; D008432:Mathematical Computing; D008956:Models, Chemical; D009010:Monte Carlo Method; D011487:Protein Conformation; D011489:Protein Denaturation; D013499:Surface Properties; D014961:X-Ray Diffraction
D007611:Aprotinin