It is generally accepted that a protein's primary sequence determines its three-dmenslonal structure. It has proved difficult, however, to obtain detailed structural information about the actual protein folding process and intermediate states. We present the results of molular dynamics simulations of the unfoldin of reduced bovine pancreatic trypsin inhibitor. The resulting partially "denatured" state was compact but expanded relative to the native state (11-25%); the expansion was not caused by an influx of water molecules. The structures were mobile, with overal secondary structure contents comparable to those of the native protein.The protein experienced relatively local unfolding, with the largest changes in the structure occurring in the loop regions.A hydrophobic core was maintained althongh ping of the side chains was compromised. The properties displayed in the simulation are consistent with unfolding to a molten globule state. Our smulations provide an in-depth view of this state and details of water-protein interactions that cannot yet be obtained experimentally.A protein must assume a stable and precisely ordered conformation to perform its biological function properly. Although much is known about these conformations and their synthesis, little is known about the detailed structures and dynamic transitions that occur during protein folding. Determination of the structural characteristics of partially folded intermediate states are crucial for understanding the mechanism of folding. Unfortunately, the cooperative nature of folding results in only minute amounts of these intermediates at equilibrium. So, instead, kinetic means have been used to characterize transiently populated intermediates (1-3). However, it has been shown that a few proteins adopt stable, partially folded equilibrium intermediates under unfolding conditions (4). There is some evidence that these intermediates, termed "molten globules," also occur along the kinetic pathway (ref. 4 and references therein). However, the molecular details of even the equilibrium folding intermediates remain difficult to characterize experimentally.Here we present the results of molecular dynamics (MD) simulations of protein unfolding and the resulting partially unfolded structures. The approach of studying unfolding instead of folding has computational advantages, because simulations proceed from a well-defined starting structure. In principle, MD is capable of providing detailed information about unfolding and the resulting unfolded state that cannot be obtained experimentally. This method has yielded realistic representations ofa native protein in water (5, 6) and has been useful in understanding the nature of less-structured peptides (7-9) and the unfolding-folding transitions that they undergo (7,10,11 (14). All atoms were explicitly present during the simulations, and the potential energy function and associated parameters are based on those described earlier (15, 16), except for the parameters for chloride (R* = 2.35 A; e* = 0.10 kcal/mol...