Packing as a structural basis of protein stability: understanding mutant properties from wildtype structure

Modeling of protein core mutations using sidechain packing can forecast their effects on stability. We have assessed the structural basis of this approach, by evaluating the accuracy of our 1991 model of a three-site mutant of lambda repressor (V36L/M40L/V47I), against the recently reported crystal structure. The three mutated residues matched the crystal structure to within 0.89A (1.11A for sidechain atoms), giving fairly accurate sidechain placement and packing (81-99th percentile rank in coordinate accuracy). However, the model used different sidechain torsional angles than seen in the crystal structure at residues 36 and 40, apparently to compensate for the backbone shifts present in the actual mutant structure, but not accounted for in our modeling method. To understand the structural basis of stability across a set of lambda repressor core mutants, we have analyzed the mutant models, revealing several simple packing effects: V36I, predicted to be stabilized by filling a hydrophobic cavity; M40V, destabilized by a steric clash with the unusual structural demands of a helix-turn transition. These effects illustrate how mutant stability can often be understood directly from scrutiny of wildtype structure. Simply adding the calculated energies of neighboring point mutations predicts the stability effect of the combined mutant relatively well, with little apparent cooperativity, yielding simple rules for each site's amino acid preferences. Our treatment of core packing indicates that it can permit a large fraction of sequences to fit the native fold, as observed experimentally, far more than indicated by rotamer hard-sphere models.