On the unfolding of α‐lytic protease and the role of the pro region

Molecular dynamics simulations of α‐lytic protease (αLP) alone and complexed with its pro region (PRO) are performed to understand the origin of its high unfolding (and folding) barrier when it is alone and how the pro region lowers this barrier. At room temperature, αLP exhibits lower dynamic fluctuations than α‐chymotrypsin. Simulation of PRO alone led to reorientation of its N terminal helix and collapse to a more compact state. A model for the uncleaved proenzyme was built and found to be stable in the time scale of the simulations. Energetic analysis suggests that the origin of strain in the uncleaved proenzyme compared with the cleaved complex is in the intramolecular backbone electrostatic interactions of the cleaved strand. In high temperature simulations, the interaction of the long beta hairpin of the enzyme with the C terminal beta sheet of PRO is among the most stable in the complex and a likely “nucleation site” for folding. In the course of unfolding, the C terminal tail of PRO is sometimes observed to intervene between the long hairpin and the aspartate loop of the enzyme, perhaps thereby lowering the energy barrier for separation of the two hairpins. Tighter interactions at the interface between the enzyme and its pro region are also occasionally observed, providing an additional mechanism for unfolding catalysis. Simulations of a mutant enzyme where the buried ion pair residues R102 and D142 were replaced by W and L, respectively, did not display any distinguishable behavior compared with the wild type. Proteins 2000;41:21–32. © 2000 Wiley‐Liss, Inc.