Unfolding and refolding of the native structure of bovine pancreatic trypsin inhibitor studied by computer simulations.

A new procedure for studying the folding and unfolding of proteins, with an application to bovine pancreatic trypsin inhibitor (BPTI), is reported. The unfolding and refolding of the native structure of the protein are characterized by the dimensions of the protein, expressed in terms of the three principal radii of the structure considered as an ellipsoid. A dynamic equation, describing the variations of the principal radii on the unfolding path, and a numerical procedure to solve this equation are proposed. Expanded and distorted conformations are refolded to the native structure by a dimensional-constraint energy minimization procedure. A unique and reproducible unfolding pathway for an intermediate of BPTI lacking the [30,51] disulfide bond is obtained. The resulting unfolded conformations are extended; they contain near-native local structure, but their longest principal radii are more than 2.5 times greater than that of the native structure. The most interesting finding is that the majority of expanded conformations, generated under various conditions, can be refolded closely to the native structure, as measured by the correct overall chain fold, by the rms deviations from the native structure of only 1.9-3.1 A, and by the energy differences of about 10 kcal/mol from the native structure. Introduction of the [30,51] disulfide bond at this stage, followed by minimization, improves the closeness of the refolded structures to the native structure, reducing the rms deviations to 0.9-2.0 A. The unique refolding of these expanded structures over such a large conformational space implies that the folding is strongly dictated by the interactions in the amino acid sequence of BPTI. The simulations indicate that, under conditions that favor a compact structure as mimicked by the volume constraints in our algorithm, the expanded conformations have a strong tendency to move toward the native structure; therefore, they probably would be favorable folding intermediates. The results presented here support a general model for protein folding, i.e., progressive formation of partially folded structural units, followed by collapse to the compact native structure. The general applicability of the procedure is also discussed.

[1]  H A Scheraga,et al.  Assessment of some problems associated with prediction of the three-dimensional structure of a protein from its amino-acid sequence. , 1975, Proceedings of the National Academy of Sciences of the United States of America.

[2]  T. Creighton,et al.  Kinetic role of a meta-stable native-like two-disulphide species in the folding transition of bovine pancreatic trypsin inhibitor. , 1984, Journal of molecular biology.

[3]  J. Skolnick,et al.  Discretized model of proteins. I. Monte Carlo study of cooperativity in homopolypeptides , 1992 .

[4]  M. Levitt.  Protein folding by restrained energy minimization and molecular dynamics. , 1983, Journal of molecular biology.

[5]  H. Scheraga,et al.  Energy parameters in polypeptides. VII. Geometric parameters, partial atomic charges, nonbonded interactions, hydrogen bond interactions, and intrinsic torsional potentials for the naturally occurring amino acids , 1975 .

[6]  Y Okamoto,et al.  Beta-sheet folding of fragment (16-36) of bovine pancreatic trypsin inhibitor as predicted by Monte Carlo simulated annealing. , 1992, Protein engineering.

[7]  K. Dill,et al.  Intrachain loops in polymers: Effects of excluded volume , 1989 .

[8]  Folding of model heteropolymers by configurational-bias Monte Carlo , 1992 .

[9]  H. Scheraga,et al.  Statistical mechanics of noncovalent bonds in polyamino acids. VIII. Covalent loops in proteins , 1965 .

[10]  M Levitt,et al.  A model of the molten globule state from molecular dynamics simulations. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[11]  N. Go.  Thermodynamics of Small-Amplitude Conformational Fluctuations in Native Globular Proteins , 1980 .

[12]  A Kolinski,et al.  Dynamic Monte Carlo simulations of a new lattice model of globular protein folding, structure and dynamics. , 1991, Journal of molecular biology.

[13]  T. Creighton,et al.  Conformational restrictions on the pathway of folding and unfolding of the pancreatic trypsin inhibitor. , 1977, Journal of molecular biology.

[14]  W. Kauzmann.  Some factors in the interpretation of protein denaturation. , 1959, Advances in protein chemistry.

[15]  H. Scheraga,et al.  Introduction of short-range restrictions in a protein-folding algorithm involving a long-range geometrical restriction and short-, medium-, and long-range interactions. , 1981, Proceedings of the National Academy of Sciences of the United States of America.

[16]  D Thirumalai,et al.  The nature of folded states of globular proteins , 1992, Biopolymers.

[17]  H. Scheraga,et al.  Formation of Local Structures in Protein Folding , 1989 .

[18]  P. Flory.  Principles of polymer chemistry , 1953 .

[19]  M. Levitt,et al.  Computer simulation of protein folding , 1975, Nature.

[20]  H. Scheraga,et al.  A Method for Predicting Nucleation Sites for Protein Folding Based on Hydrophobic Contacts , 1978 .

[21]  M Karplus,et al.  Conformations of intermediates in the folding of the pancreatic trypsin inhibitor. , 1987, Journal of molecular biology.

[22]  Robert Huber,et al.  Structure of bovine pancreatic trypsin inhibitor , 1984 .

[23]  H. Scheraga,et al.  Local interactions favor the native 8-residue disulfide loop in the oxidation of a fragment corresponding to the sequence Ser-50-Met-79 derived from bovine pancreatic ribonuclease A , 1988, Journal of protein chemistry.

[24]  R. L. Baldwin.  Seeding protein folding , 1986 .

[25]  Harold A. Scheraga,et al.  Consideration of the Hydrodynamic Properties of Proteins1,2 , 1953 .

[26]  C. Tanford,et al.  Interfacial free energy and the hydrophobic effect. , 1979, Proceedings of the National Academy of Sciences of the United States of America.

[27]  S. Hawley,et al.  An electrophoretic study of reversible protein denaturation: chymotrypsinogen at high pressures. , 1975, Biochemistry.

[28]  Homer Jacobson,et al.  Intramolecular Reaction in Polycondensations. I. The Theory of Linear Systems , 1950 .

[29]  E. Freire,et al.  Structural energetics of the molten globule state , 1993, Proteins.

[30]  P. Flory,et al.  Macrocyclization Equilibrium Constants and the Statistical Configuration of Poly(dimethylsiloxane) Chains , 1966 .

[31]  I. Kuntz,et al.  Disulfide bond-coupled folding of bovine pancreatic trypsin inhibitor derivatives missing one or two disulfide bonds. , 1992, Biochemistry.

[32]  S W Englander,et al.  Stable submolecular folding units in a non-compact form of cytochrome c. , 1991, Journal of molecular biology.

[33]  H. Scheraga,et al.  Energy parameters in polypeptides. 9. Updating of geometrical parameters, nonbonded interactions, and hydrogen bond interactions for the naturally occurring amino acids , 1983 .

[34]  H A Scheraga,et al.  Variable-target-function and build-up procedures for the calculation of protein conformation. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data. , 1988, Journal of biomolecular structure & dynamics.

[35]  T Schlick,et al.  Supercoiled DNA energetics and dynamics by computer simulation. , 1992, Journal of molecular biology.

[36]  Terrence G. Oas,et al.  A peptide model of a protein folding intermediate , 1988, Nature.

[37]  M. Hao,et al.  Effects of compact volume and chain stiffness on the conformations of native proteins. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[38]  T. Schlick,et al.  Molecular dynamics by the backward-Euler method , 1989 .

[39]  H A Scheraga,et al.  Folding of the twisted beta-sheet in bovine pancreatic trypsin inhibitor. , 1985, Biochemistry.

[40]  D. Covell,et al.  Conformations of folded proteins in restricted spaces. , 1990, Biochemistry.

[41]  F E Cohen,et al.  Protein folding. Effect of packing density on chain conformation. , 1991, Journal of molecular biology.

[42]  H. Scheraga,et al.  Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions for use in the ECEPP/3 algorithm, with application to proline-containing peptides , 1994 .

[43]  P. S. Kim,et al.  Reexamination of the folding of BPTI: predominance of native intermediates , 1991, Science.

[44]  A. Kossiakoff,et al.  Structural effects induced by removal of a disulfide-bridge: the X-ray structure of the C30A/C51A mutant of basic pancreatic trypsin inhibitor at 1.6 A. , 1990, Protein engineering.

[45]  H. Scheraga,et al.  Calculation of protein conformation by the build-up procedure. Application to bovine pancreatic trypsin inhibitor using limited simulated nuclear magnetic resonance data. , 1988, Journal of biomolecular structure & dynamics.

[46]  K. Dill,et al.  Protein core assembly processes , 1993 .