Journal:Proteins:2

From Proteopedia

Jump to: navigation, search

Protein Stability and in Vivo Concentration of Missense Mutations in Phenylalanine Hydroxylase

Zhen Shi, Jenn Sellers, and John Moult [1]


Molecular Tour
A previous computational analysis of missense mutations linked to monogenic disease found a high proportion of missense mutations affect protein stability, rather than other aspects of protein structure and function. The purpose of the present study is to relate the presence of such stability damaging missense mutations to the levels of a particular protein present under in vivo like conditions, and to test the reliability of the computational methods. Experimental data on a set of missense mutations of the enzyme phenylalanine hydroxylase (PAH) associated with the monogenic disease phenylketonuria (PKU) have been compared with the expected in vivo impact on protein function, obtained using SNPs3D, an in silico analysis package. A high proportion of the PAH mutations are predicted to be destabilizing. The overall agreement between predicted stability impact and experimental evidence for lower protein levels is in accordance with the estimated error rates of the methods. For these mutations, destabilization of protein three dimensional structure is the major molecular mechanism leading to PKU, and results in a substantial reduction of in vivo PAH protein concentration. Although of limited scale, the results support the view that destabilization is the most common mechanism by which missense mutations cause monogenic disease. In turn, this conclusion suggests the general therapeutic strategy of developing drugs targeted at restoring wild type stability.

Under physiological conditions, human PAH is a homo-tetramer, with each subunit composed of three domains. From N terminal to C terminal these are the regulatory, catalytic and tetramerization domains. To date, no experimentally determined structure of the complete human molecule is available. Three PDB structures were selected to model specific mutations in different domains based on crystal structure resolution, structure quality, and coverage: 1j8u, human PAH structure containing mainly the catalytic domain (monomeric form); 2pah, human PAH structure covering the catalytic and tetramerization domains (tetrameric complex); and 1phz, rat PAH structure covering the regulatory and catalytic domains (dimeric complex). The high resolution human 1j8u structure was used to model catalytic domain mutations. Regulatory domain mutations were modeled using a homology model of the human domain, based on the rat 1phz structure, as were three catalytic domain mutations, R261Q, R413P, and Y414C, that are in contact with the regulatory domain across a subunit interface. Rat PAH protein has 93% sequence identity with human PAH. There are no insertions or deletions in sequence between the two proteins. Main chain coordinates were taken directly from the rat structure. Side chains conformations were optimized using SCRWL. Catalytic domain mutations R408W and R408Q are in contact with the tetramerization domain of another subunit and were modeled using 2pah.

Composite model is build from PDB structures 1j8u, 2pah, and 1phz. Domains are: regulatory (yellow); catalytic (green); tetramerization (blue). The substrate L-Phe (substrate analog beta(2-thienyl) alanine (TIH) is shown) and cofactor tetrahydrobiopterin (BH4 or H4B) both have binding sites in the catalytic domain. The TIH, Fe (++) ion and cofactor Tetrahydrobiopterin (BH4) are shown space filled.

Category 1: 28 missense mutations are expected to affect stability only

28 of the 35 mutations with destabilization assignments are remote from any known ligand binding or the catalytic site, and so are expected to have a low experimental protein level, and wild type specific activity.

16 of the 28 (F39L, G46S, L48S, I65T, A104D, P122Q, R157N, F161S, R243Q, R252G, R252Q, R252W, A259T, A259V, L311P, R408W, in orangered) have protein levels less than 50% wild type, as expected. Of these, all but two have wild type specific activity. The two exceptions, F39L and L48S, have approximately three fold higher specific activities than the wild type. These mutations lie in the regulatory domain, suggesting a possible explanation for the high activity level. The 16 mutants are classified into clinical categories of mild PKU (A104D), moderate PKU (F39L, L48S, I65T), and classic PKU (G46S, R243Q, R252G/Q/W, A259V, L311P, R408W).

Selected mutations from this group:

  • A104D, this mutation probably decreased hydrophobic interaction and formed new hydrogen bonds. and the . Click here to see animation of this scene.
  • Mutation R157N caused saltbridge lost and hydrogen bonds lost. and the . Click here to see animation of this scene. Water molecules shown as red spheres.
  • Mutation R243Q caused saltbridge lost. and the . Click here to see animation of this scene.
  • Mutation A259V caused overpacking. and the . Click here to see animation of this scene.
  • Mutation A259T caused overpacking and buried polar. and the . Click here to see animation of this scene.
  • Mutation R408W caused hydrogen bonds lost. and the . Click here to see animation of this scene.

Mutations R252G/Q/W caused saltbridge lost and hydrogen bonds lost:

  • Mutation R252G. and the . Click here to see animation of this scene.
  • Mutation R252Q. and the . Click here to see animation of this scene.
  • Mutation R252W. and the . Click here to see animation of this scene.

Nine of remaining mutations expected to affect stability only (L41F, R68G, R68S, E76G, G218V, P244L, A309V, A403V, R408Q, in blueviolet) have reported experimental protein levels greater than 50% of wild type (all 100%, except one of the R408Q experiments with 70%), inconsistent with the computational assignment.

Selected mutations from this group:

  • Mutation L41F caused overpacking; gain of hydrophobic interaction. and the . Click here to see animation of this scene.
  • Mutations R68G/S caused saltbridge lost and hydrogen bond lost; hydrophobic interaction decreased (Mutation R68S is shown). and the . Click here to see animation of this scene.
  • Mutation E76G caused hydrogen bonds lost. and the . Click here to see animation of this scene.
  • Mutation P244L caused overpacking; gain of hydrophobic interaction. and the . Click here to see animation of this scene.
  • Mutation R408Q caused hydrogen bonds lost. and the . Click here to see animation of this scene.

Category 2: Seven missense mutations are expected to affect both stability and molecular function

There are seven mutations (G247V, L255S, R270S, E280K, S349L, S349P and Y277D, in magenta) with atomic contacts of 6.5 Å or less to the phenylalanine substrate (substrate analog beta(2-thienyl) alanine (TIH) is shown), the BH4 cofactor or the Fe++ ion, and that are assigned as destabilizing by the structure SVM. Residues interacting with TIH, BH4 and the Fe++ ion are in green. Water molecules shown as red spheres. These mutant proteins are therefore expected to exhibit a combination of lower specific activity and a lower total protein level. Six of the seven (G247V, L255S, R270S, E280K, S349L, S349P) have protein levels less than half or in one case close to half (G247V, 56%) that of wild type, and very low protein activity, consistent with expectations. Clinical categories are available for E280K, S349L, and S349P, and are all “classic PKU”, consistent with the results and with experiment. The remaining mutant in this category, Y277D, has an experimental activity of zero, and is classified as mild or classic PKU, consistent with the profile SVM assignments. But the measured protein level is reported as 99% of wild type, inconsistent with a modest confidence stability assignment. This may be a computational false positive with respect to stability.

  • Mutation G247V caused overpacking. and the . Click here to see animation of this scene.
  • Mutation L255S caused decrease of hydrophobic interaction.
  • Mutation R270S caused saltbridge lost and hydrogen bond lost; hydrophobic interaction decreased. and the . Click here to see animation of this scene.
  • Mutation E280K caused saltbridge lost. and the . Click here to see animation of this scene.
  • Mutation S349L caused hydrogen bond lost. and the . Click here to see animation of this scene.
  • Mutation S349P caused hydrogen bond lost. and the . Click here to see animation of this scene
  • Mutation Y277D caused decrease of hydrophobic interaction and probably formed new hydrogen bond. and the . Click here to see animation of this scene.

Category 3: Nine mutations are expected to impact molecular function only

A total of nine mutations are classified as high impact by the sequence conservation method, classified as not destabilizing by the stability method, and so are expected to impact molecular function but not stability, implying wild type protein levels and lower activity. Four of these, L255V, P281L, A322G, and L348V (in cyan) have atomic contacts of 6.5 Å or less to a ligand. Experimental data for two, A322G and L348V, are consistent with expectations, with low activity and normal protein levels:

  • Mutation A322G. and the . Click here to see animation of this scene.
  • Mutation L348V caused cavity creation; hydrophobic interaction decreased. and the . Click here to see animation of this scene.

The remaining two, L255V and P281L, have low activity, but also low protein level. Both are in direct contact with the BH4 cofactor, and would disrupt binding substantially:

  • Mutation L255V. and the . Click here to see animation of this scene.
  • Mutation P281L, this mutation probably caused clashes with ligands BH4 and TIH. and the . Click here to see animation of this scene.

The other five mutations in this category, K42I, D59Y, D143G, V388M and R413P (in deepskyblue), are not near to any known ligand binding or catalytic site.

Mutations from this group:

  • Mutation K42I caused saltbridge lost; is located on surface of regulatory domain. and the . Click here to see animation of this scene.
  • Mutation D59Y is located on the surface of the regulatory domain. and the . Click here to see animation of this scene.
  • Mutation D143G caused hydrogen bond lost and disruption of hydrogen bonds network. and the . Click here to see animation of this scene.
  • Mutation V388M probably cause overpacking. and the . Click here to see animation of this scene.
  • Mutation R413P caused saltbridge lost. and the . Click here to see animation of this scene.

Category 4: Two mutations are assigned low impact by both the sequence conservation and stability methods

Two mutations, T92I, and P211T (in darkmagenta), are assigned low impact by both computational methods. Both sets of experimental results show close to normal activity and protein levels, consistent with the analysis results. Also reasonably consistent, T92I is assigned to the mild MHP category of disease, suggesting a subtle effect on protein function. Inconsistent with both experiment and computational analysis, P211T is assigned to the “classic PKU” category, based on a single functionally hemizygous patient genotype.

Evolutionary conservation

Evolutionary conservation

Image:Consurf_key_small.gif


PDB ID 2pah

Drag the structure with the mouse to rotate
  1. Shi Z, Sellers J, Moult J. Protein stability and in vivo concentration of missense mutations in phenylalanine hydroxylase. Proteins. 2012 Jan;80(1):61-70. doi: 10.1002/prot.23159. Epub 2011 Sep 21. PMID:21953985 doi:http://dx.doi.org/10.1002/prot.23159

Proteopedia Page Contributors and Editors (what is this?)

Alexander Berchansky, Jaime Prilusky, Joel L. Sussman

This page complements a publication in scientific journals and is one of the Proteopedia's Interactive 3D Complement pages. For aditional details please see I3DC.
Personal tools