A new approach in implementing classical molecular dynamics simulation for parallel computers has enabled a simulation to be carried out on a protein with explicit representation of water an order of magnitude longer than previously reported and will soon enable such simulations to be carried into the microsecond time range. We have used this approach to study the folding of the villin headpiece subdomain, a 36-residue small protein consisting of three helices, from an unfolded structure to a molten globule state, which has a number of features of the native structure. The time development of the solvation free energy, the radius of gyration, and the mainchain rms difference from the native NMR structure showed that the process can be seen as a 60-nsec ''burst'' phase followed by a slow ''conformational readjustment'' phase. We found that the burial of the hydrophobic surface dominated the early phase of the folding process and appeared to be the primary driving force of the reduction in the radius of gyration in that phase.Understanding the mechanism of protein folding has been a grand challenge in protein chemistry for a few decades. Important advances in this area can have a profound impact in many biologically relevant fields. A direct benefit from the understanding of the mechanism should be the ability to predict protein structures more accurately, which in turn should enable the pharmaceutical industry to improve drug discovery. The recent hypothesis of folding-related diseases is another example of the significance of folding (1, 2). Despite great progress made by a variety of experimental and theoretical approaches, it has been difficult to establish detailed descriptions of the folding process and mechanism at an atomic level.Computer simulation has been a powerful tool that can provide rich information at various levels of resolution. Simulation has been instrumental in understanding protein folding mechanisms. One such approach is lattice model simulations (3-5). These types of models use reduced representations by treating the residues as (one or two) linked beads. A more detailed approach is the atomic level model, either unitedatom or all-atom, which represents all or most of the atoms of the protein explicitly, with an implicit representation of the solvent (6, 7). Molecular dynamics (MD) simulations with full atomic representation of both protein and solvent possess a unique advantage to study protein folding because of their atomic level resolution and accuracy. This method with associated simulation parameters (i.e., the force field) derived from experiments and from gas-phase quantum mechanical calculations has been tested by using smaller molecular systems in many comparisons with experimental results.Because of the complexity of this approach, the large number of atoms, and the need to take time steps of 1 to 2 fsec, such simulations have to date been limited to a few nanoseconds; the longest single MD trajectories of proteins with explicit water have been Ͻ10 nsec (8). This has pre...