Advances in simulation techniques and computing hardware have created a substantial overlap between the timescales accessible to atomic-level simulations and those on which the fastest-folding proteins fold. Here we demonstrate, using simulations of four variants of the human villin headpiece, how simulations of spontaneous folding and unfolding can provide direct access to thermodynamic and kinetic quantities such as folding rates, free energies, folding enthalpies, heat capacities, Φ-values, and temperaturejump relaxation profiles. The quantitative comparison of simulation results with various forms of experimental data probing different aspects of the folding process can facilitate robust assessment of the accuracy of the calculations while providing a detailed structural interpretation for the experimental observations. In the example studied here, the analysis of folding rates, Φ-values, and folding pathways provides support for the notion that a norleucine double mutant of villin folds five times faster than the wild-type sequence, but following a slightly different pathway. This work showcases how computer simulation has now developed into a mature tool for the quantitative computational study of protein folding and dynamics that can provide a valuable complement to experimental techniques.Amber ff99SB*-ILDN | enthalpy | heat capacity | pre-exponential factor | transition path time P roteins are synthesized in the cell or in vitro as unstructured polypeptide chains that, in most cases, self-assemble into their functionally active three-dimensional shapes. This process, called protein folding, occurs on a broad range of timescales ranging from microseconds to seconds and higher. From a purely physical-chemical perspective, it should be possible in principle to characterize the folding mechanism of a given protein at atomistic resolution and to reconstruct its free-energy landscape, given only its primary sequence, through molecular dynamics (MD) simulations based on elementary physical principles. This direct approach has been rarely pursued because even the simplest systems representing a protein immersed in water consist of several thousand atoms, and simulating their behavior on the timescales typical of protein folding is computationally extremely demanding. The discovery and design of fast-folding proteins (1) significantly narrowed the timescale gap between simulations and experiments, making such simulations feasible, at least for the fastest-folding proteins.The C-terminal fragment of the villin headpiece [referred to in the remainder of this paper simply as "villin" (2)], one of the fastest-folding protein domains known (3), has proven to be an excellent target for folding simulations with physics-based force fields and an atomistically detailed representation of both the solute and the surrounding solvent (4-7). Until recently, the length of such simulations was limited to a few microseconds-a timescale sufficient to capture, at best, a single folding event (8,9). With this limitation, it has bee...