Toward mechanistic models for genotype-phenotype correlations in phenylketonuria using protein stability calculations

Phenylketonuria (PKU) is a genetic disorder caused by variants in the gene encoding phenylalanine hydroxylase (PAH), resulting in accumulation of phenylalanine to neurotoxic levels. Here, we analyzed the cellular stability, localization, and interaction with wild-type PAH of 20 selected PKU-linked PAH protein missense variants. Several were present at reduced levels in human cells, and the levels increased in the presence of a proteasome inhibitor, indicating that proteins are proteasome targets. We found that all the tested PAH variants retained their ability to associate with wild-type PAH, and none formed aggregates, suggesting that they are only mildly destabilized in structure. In all cases, PAH variants were stabilized by the cofactor tetrahydrobiopterin (BH₄ ), a molecule known to alleviate symptoms in certain PKU patients. Biophysical calculations on all possible single-site missense variants using the full-length structure of PAH revealed a strong correlation between the predicted protein stability and the observed stability in cells. This observation rationalizes previously observed correlations between predicted loss of protein destabilization and disease severity, a correlation that we also observed using new calculations. We thus propose that many disease-linked PAH variants are structurally destabilized, which in turn leads to proteasomal degradation and insufficient amounts of cellular PAH protein.