In this study we evaluate, at full atomic detail, the folding processes of two small helical proteins, the B domain of protein A and the Villin headpiece. Folding kinetics are studied by performing a large number of ab initio Monte Carlo folding simulations using a single transferable all-atom potential. Using these trajectories, we examine the relaxation behavior, secondary structure formation, and transition-state ensembles (TSEs) of the two proteins and compare our results with experimental data and previous computational studies. To obtain a detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Moreover, rigorous pfold analysis is used to obtain representative samples of the TSEs and a good quantitative agreement between experimental and simulated ⌽ values is obtained for protein A. ⌽ values for Villin also are obtained and left as predictions to be tested by future experiments. Our analysis shows that the two-helix hairpin is a common partially stable structural motif that gets formed before entering the TSE in the studied proteins. These results together with our earlier study of Engrailed Homeodomain and recent experimental studies provide a comprehensive, atomic-level picture of folding mechanics of three-helix bundle proteins.transition state ensemble ͉ Villin ͉ protein A A n eventual solution to the protein-folding problem will involve a close calibration of theoretical methods to experimental data (1-4). In the endeavor of obtaining a quantitative agreement between theory and experiments, two small ␣-helical proteins have played a central role, namely the B domain of protein A from Staphylococcus aurues and the Villin headpiece subdomain from chicken. Although these proteins belong to different SCOP fold classes (5), both have simple three-helix bundle native topologies and fold autonomously on the microsecond time scale (6, 7), which makes them ideal test cases for protein simulations and numerous simulation studies, ranging from simple C ␣ Go-type to all-atom models with explicit water, have been undertaken for both protein A (8-21) and Villin (16,17,[22][23][24][25][26][27][28][29].Important advances have been made toward agreements with experiments for both proteins, but several key issues remain unresolved (6, 30, 31). The need for additional studies also is emphasized by recent experiments. Fersht et al. (31,32) performed a comprehensive mutational analysis on protein A by obtaining ⌽ values at Ͻ30 aa positions, providing an important benchmark for simulation studies. The obtained ⌽ values suggest that the transition-state ensemble (TSE) is characterized mainly by a well formed H2 (we denote the three individual helices from N-to C-terminal by H1, H2, and H3, following previous convention) stabilized by hydrophobic interactions with H1. Recent experimental studies of Villin (6, 33) have focused mainly on achieving fast-folding mutants, although new biophysical characterization of wild-type Villin a...