Global computational mutagenesis provides a critical stability framework in protein structures

A protein’s amino acid sequence dictates the folds and final structure the macromolecule will form. We propose that by identifying critical residues in a protein’s atomic structure, we can select a critical stability framework within the protein structure essential to proper protein folding. We use global computational mutagenesis based on the unfolding mutation screen to test the effect of every possible missense mutation on the protein structure to identify the residues that cannot tolerate a substitution without causing protein misfolding. This method was tested using molecular dynamics to simulate the stability effects of mutating critical residues in proteins involved in inherited disease, such as myoglobin, p53, and the 15(th) sushi domain of complement factor H. In addition we prove that when the critical residues are in place, other residues may be changed within the structure without a stability loss. We validate that critical residues are conserved using myoglobin to show that critical residues are the same for crystal structures of 6 different species and comparing conservation indices to critical residues in 9 eye disease-related proteins. Our studies demonstrate that by using a selection of critical elements in a protein structure we can identify a critical protein stability framework. The frame of critical residues can be used in genetic engineering to improve small molecule binding for drug studies, identify loss-of-function disease-causing missense mutations in genetics studies, and aide in identifying templates for homology modeling.