We study the folding mechanism of a three-helix bundle protein at atomic resolution, including effects of explicit water. Using replica exchange molecular dynamics we perform enough sampling over a wide range of temperatures to obtain the free energy, entropy, and enthalpy surfaces as a function of structural reaction coordinates. Simulations were started from different configurations covering the folded and unfolded states. Because many transitions between all minima at the free energy surface are observed, a quantitative determination of the free energy barriers and the ensemble of configurations associated with them is now possible. The kinetic bottlenecks for folding can be determined from the thermal ensembles of structures on the free energy barriers, provided the kinetically determined transition-state ensembles are similar to those determined from free energy barriers. A mechanism incorporating the interplay among backbone ordering, sidechain packing, and desolvation arises from these calculations. Large ⌽ values arise not only from native contacts, which mostly form at the transition state, but also from contacts already present in the unfolded state that are partially destroyed at the transition. T he energy landscape theory and the funnel concept, along with a new generation of experiments, have created a ''new view'' of protein folding (1-4). Folding is best described as an ensemble of converging pathways toward the native structure. For good folding sequences, these paths are all energetically very similar with barriers between them on the k B T energy scale. Minor fluctuations due to environmental or mutational changes may vary the path probabilities, but the overall global landscape flow will be only weakly altered. This leads to a key result of the energy landscape theory: protein folding dynamics can be properly understood as the diffusion of an ensemble of protein configurations over a low-dimensional free energy surface, which may be constructed by using many different order parameters.In this article, we describe the free energy landscape for protein folding simulated with full atomic details of both protein and solvent, over a broad range of temperatures. For our studies, a natural choice is the 10-55 helical fragment B of protein A from Staphylococcus aureus, using the replica exchange molecular dynamics (REMD) algorithm described by Sugita and Okamoto (5). Protein A folds into a simple three-helix bundle (6) whose folding has been widely studied by using minimalist and all-atom simulations (7-17) as well as experiments (18,19). All of these different studies provided some information about the folding mechanism and the nature of the transition-state ensemble (TSE). The study presented here aims to establish a quantitative picture by integrating the earlier qualitative findings with detailed simulations. We simulated 82 replicas of the water-protein system, with temperature, T, ranging from 277 to 548 K, and each replica was started from a different configuration spanning the space of folde...