Development of simulation models for protein folding in a thermal annealing process -I: a simulation of BPTI folding by the pearl necklace model

A model system is proposed to simulate the folding processesof proteins during thermal annealing. This system consists offour subsystems: (i) the pearl necklace model with isotropicinter-residue interactions; (ii) the extended pearl necklacemodel with anisotropic interaction potentials; (iii) moltenglobule phase dynamics; and (iv) final generation of the three-dimensionalstructure of a given protein. In this paper results obtainedwith the pearl necklace model are reported. This model consistsof spherical elements and virtual bonds of 3.8 Å in lengthand is intended to sinudate dynamical processes at relativelyhigh temperature where entropic terms play a dominant role.Inter-residue interactions are composed of spherical soft repulsivepotentials and hydrophobic interactions inherent to respectiveresidues. A simulation of folding processes of BPTI startingfrom the fully extended conformation indicated that intermediates,even at early stages of folding, are not randomly coiled butassume organized structures that resemble, to some extent, thenative conformation.