An effort to combine theoretical analyses and protein engineering methods has been made to probe the folding mechanism of SH3 by using Energy Landscape Theory and a -value analysis. Particular emphasis was given to core residues and the effect of desolvation during the folding event by replacing the core valines with isosteric threonines. These mutations have the advantage of keeping the core structurally invariant while affecting core stability relative to the unfolded state. Although the valines that form the core appear spatially invariant, the folding kinetics of their threonine mutants varies, indicating their different extent of solvation in the transition-state ensemble. Theoretical studies predicted the distribution of folding kinetics of threonine mutants without previous knowledge of the measured rates. This initial success encourages further investigations of the molecular details behind these macroscopic phenomena and of the role of solvation in the folding mechanism.A large body of recent data suggests that proteins, especially small fast-folding (submillisecond) proteins (1-3), have sequences with a level of energetic frustration sufficiently reduced that their overall energy landscape (4-8) resembles a moderately rough funnel (9). The landscape roughness corresponds to local free energy minima arising from geometric (topological) and energetic traps (10). Geometric or topological traps are associated with chain connectivity and the shape of the native fold and occur when correct contacts form prematurely (10). Because the energetic roughness is minimal, the structural heterogeneity observed in the transitionstate ensemble (TSE) is strongly influenced by the topological effects, which, to a large extent, may be inferred from the native structure (11-16). This explains why simple energetically unfrustrated (Go -like) models (17) reproduce nearly all experimental results for the global geometrical features of the TSE and͞or intermediates of a large number of real proteins, which are two-or three-state folders (13). These models have been recently generalized to include the ability to capture effects such as the desolvation of the hydrophobic core (14).Using this theoretical framework, the funnel landscape coupled with microscopic desolvation can be pictured by using two ordering parameters: Q, a fraction of the total native contact formation, and pseudo Q, a fraction of the total single-water separated native contact formation (14). The latter is a direct consequence of using a pairwise potential with a characteristic desolvation barrier, which accounts for entropy costs to expel a single water molecule between two interacting apolar groups before a full contact can be made (14,(18)(19)(20). Although no explicit water molecule has been included in the simulations, the emergence of this desolvation barrier emulates the granular properties of solvents. By using this model, we have recently predicted that the folding of SH3 is controlled by a structural-search collapse followed by desolvation of the hyd...