Ab initio simulations of the folding pathways are currently limited to very small proteins. For larger proteins, some approximations or simplifications in protein models need to be introduced. Protein folding and unfolding are among the basic processes in the cell and are very difficult to characterize in detail by experiment or simulation. Chymotrypsin inhibitor 2 (CI2) and barnase are probably the best characterized experimentally in this respect. For these model systems, initial folding stages were simulated by using CA-CB-side chain (CABS), a reduced-space protein-modeling tool. CABS employs knowledge-based potentials that proved to be very successful in protein structure prediction. With the use of isothermal Monte Carlo (MC) dynamics, initiation sites with a residual structure and weak tertiary interactions were identified. Such structures are essential for the initiation of the folding process through a sequential reduction of the protein conformational space, overcoming the Levinthal paradox in this manner. Furthermore, nucleation sites that initiate a tertiary interactions network were located. The MC simulations correspond perfectly to the results of experimental and theoretical research and bring insights into CI2 folding mechanism: unambiguous sequence of folding events was reported as well as cooperative substructures compatible with those obtained in recent molecular dynamics unfolding studies. The correspondence between the simulation and experiment shows that knowledge-based potentials are not only useful in protein structure predictions but are also capable of reproducing the folding pathways. Thus, the results of this work significantly extend the applicability range of reduced models in the theoretical study of proteins.protein structure prediction ͉ Monte Carlo simulations ͉ protein denatured state ͉ folding nucleus ͉ residual structure A large number of folded protein structures was determined by x-ray crystallography or NMR. For a few proteins, the folding intermediates were characterized by using protein engineering (1) and NMR techniques (2, 3). The major transition states (TS) of chymotrypsin inhibitor 2 (CI2) and barnase were mapped at the level of individual residues by protein engineering (4, 5). Much less is known, however, about early folding events. It is very important to understand how protein folding is initiated and how the native structure emerges subsequently. The denatured state, an ensemble of partially folded, highly mobile conformations, is very difficult to study, although there have been recent reports of NMR studies of residual structure in denatured proteins. Such structures, along with hydrophobic clusters, were discovered even under highly denaturing conditions (6, 7). Moreover, denatured proteins can exhibit a longrange ordering of native-like topology (8). Therefore, the folding process can be directed from the very beginning when starting from a specific structure (9, 10). It becomes evident that the denatured state plays a crucial role in all aspects of protein ...