To assist in the efficient design of protein cavities, we have developed a minimization strategy that can predict with accuracy the fate of cavities created by mutation. We first modelled, under different conditions, the structures of six T4 lysozyme and cytochrome c peroxidase mutants of known crystal structure (where long, hydrophobic, buried side chains have been replaced by shorter ones) by minimizing the virtual structures derived from the corresponding wild-type co-ordinates. An unconstrained pathway together with an all-atom atom representation and a steepest descent minimization yielded modelled structures with lower root mean square deviations (r.m.s.d) from the crystal structures than other conditions. To test whether the method developed was generally applicable to other mutations of the kind, we have then modelled eighteen additional T4 lysozyme, barnase and cytochrome c peroxidase mutants of known crystal structure. The models of both cavity expanding and cavity collapsing mutants closely fit their crystal structures (average r.m.s.d. 0.33 ± 0.25 Å, with only one poorer prediction: L121A). The structure of protein cavities generated by mutation can thus be confidently simulated by energy minimization regardless of the tendency of the cavity to collapse or to expand. We think this is favoured by the fact that the typical response observed in these proteins to cavity-creating mutations is to experience only a limited rearrangement.