Analysis of the factors and implications of an empirical method for estimating the stability of mutant human haemoglobins

Thermitase, a thermostable subtilisin: comparison of predicted and experimental structures and the molecular cause of thermostability

The subtilisin family of proteases has four members of known sequence and structure: subtilisin Carlsberg, subtilisin novo, proteinase K, and thermitase. Using thermitase as a test case, we ask two questions. How good are methods for model building a three-dimensional structure of a protein based on sequence homology to a known structure? And what are the molecular causes of thermostability? First, we compare predicted models of thermitase, refined by energy minimization and varied by molecular dynamics, with the preliminary crystal structure. The predictions work best in the conserved structural core and less well in seven loop regions involving insertions and deletions relative to subtilisin. Here, variation of loop regions by molecular dynamics simulation in vacuo followed by energy minimization does not improve the prediction since we find no correlation between in vacuo energy and correctness of structure when comparing local energy minima. Second, in order to identify the molecular cause of thermostability we confront hypotheses derived by calculation of the details of interatomic interactions and estimates of hydrophobic interactions with inactivation experiments. As a result, we can exclude salt bridges and hydrophobic interactions as main causes of thermostability. Based on a combination of theoretical and experimental evidence, the unusually tight binding of calcium by thermitase emerges as the most likely single influence responsible for its increased thermostability.

Single amino acid substitutions producing instability of globular proteins. Calculation of their frequencies in the entire mutational spectra of the alpha- and beta-subunits of human hemoglobin

The frequencies of substitutions resulting in protein instability were calculated by a method estimating changes in stability produced by amino acid substitutions. The method takes into account the accessibility of an amino acid position to a solvent and changes in the specificity of amino acid interactions. When tested on human mutant hemoglobins, the method yielded predictions with a preciseness of 80%. The consideration of the evolutionary homologous proteins in the analysis allowed us to estimate the evolutionary constraints imposed on stability of their spatial structure. With these limitations, approximately 50% of amino acid substitutions in the entire mutational spectra of the alpha- and beta-subunits of human hemoglobin were found to damage the spatial structure of the globular proteins.