Difference between revisions of "HMM PD00008"

From PhosphataseWiki
Jump to: navigation, search
(Created page with "Back to '''List of HMMs''' '''Symbol''': CC1_Sac_PD '''Name''': Sac phosphatase domain === Description === We built the profile because the boundary of phosphatase ...")
 
(Description)
 
(8 intermediate revisions by one user not shown)
Line 6: Line 6:
  
 
=== Description ===
 
=== Description ===
We built the profile because the boundary of phosphatase domain defined by Pfam profile is not the same as the one reported in crystal structure study <cite>Begley03</cite> and SMART does not have the profile for the myotubularin phosphatase domain. The crystal structure shows the phosphatase domain ranges from 205-586 in human MTMR2 (Figure 6A in <cite>Begley03</cite>). The Pfam profile [http://pfam.xfam.org/family/PF06602.10 Myotub-related (PF06602)] predicted the range from 192-529.
+
Based upon the crystal structure of the yeast member, SAC1 has two structural domains: SacN and the catalytic phosphatase domain <cite>Manford10</cite> (see the sequences in supplementary materials). The SacN ranges from 1-181; the catalytic phosphatase domain ranges from 182-504. The SacN domain mediates the interaction with VPS74, which is proposed to mediate packaging of medial Golgi glycosyltransferases into coatomer (COP1)-coated vesicles, thereby maintaining Golgi residence <cite>Cai14</cite>.  
  
=== How the HMM was built ===
+
When using the HMM/PSSM profiles of public databases to determine the region of phosphatase domain, they gave different regions from the region determined from the first crystal structure:
We PSI-BLASTed the region from 190 to 600 of human MTMR2. The region is the phosphatase domain determined in crystal structure study with flanking regions <cite>Begley03</cite>:
+
* Pfam database: Syja_N (PF02383), 56-343
 +
* CDD database: 1) Cdd:pfam02383, 55-348, 2) Cdd:COG5329, 3-578
 +
* SMART: no hit
  
>Hsap_MTMR2:205-586 id=HsapP048_AA CC1:CC1:Myotubularin:MTMR1 [Homo sapiens]
+
Thus, Pfam profile contains part of SacN domain and part of the phosphatase domain, and it lacks the region of 451-511 is required for the PtdIns4P phosphatase activity of yeast SAC1 <cite>Cai14</cite>. COG profile contains both SacN and the catalytic phosphatase domain, but also the trans-membrane–spanning domains (523-578) (figure 4A in <cite>Cai14</cite>). We therefore built our own profiles for SacN and catalytic phosphatase domain.
GWKLYDPLLEYRRQGIPNESWRITKINERYELCDTYPALLVVPANIPDEELKRVASFRSRGRIPVLSWIHPESQATITRCSQPMVGVSGKRSKEDEKYLQAIMDSNAQSHKIFIFDARPSVNAVANKAKGGGYESEDAYQNAELVFLDIHNIHVMRESLRKLKEIVYPNIEETHWLSNLESTHWLEHIKLILAGALRIADKVESGKTSVVVHSSDGWDRTAQLTSLAMLMLDGYYRTIRGFEVLVEKEWLSFGHRFQLRVGHGDKNHADADRSPVFLQFIDCVWQMTRQFPTAFEFNEYFLITILDHLYSCLFGTFLCNSEQQRGKENLPKRTVSLWSYINSQLEDFTNPLYGSYSNHVLYPVASMRHLELWVGYYIRWNPR
+
  
We first PSI-BLASTed the region against RefSeq database via NCBI BLAST server. We found 5000+ hits in 1st round, which is difficult to look into. We then  PSI-BLASTed the region against RefSeq database via NCBI BLAST server. The coverage of human myotubularins:
+
=== How the HMM was built ===
* MTM1, MTMR1, MTMR2: ~190-600, 100%
+
We PSI-BLASTed the full sequence of yeast SAC1 (see supplementary materials) against SwissProt via NCBI BLAST server. The search converged right after 1st round. The query yeast SAC1 hit the SACs from animals and plants. We downloaded the aligned sequences, aligned them with Clustal Omega, manually adjusted the alignment, built the HMM and validated the HMM using the protein phosphatases of the nine genomes.
* MTMR3: ~214-595 (human, rat, mouse); MTMR4: ~225-588 (human, mouse), ~213-598 (xenopus)
+
* MTMR5: ~317-535 (human), ~317-579 (mouse); MTMR13 (not in swissprot?)
+
* MTMR6: ~195-595 (human, mouse); MTMR7: 192-592 (human), 203-593 (mouse); 192-593 (human, chicken)
+
* MTMR9: 192-588 (human, mouse), 204-592 (cow)
+
* MTMR10 (MTMRA): 214-532 (human, mouse, xenopus, zebra fish); MTMR11 (MTMRB): 214-561 (human, mouse); MTMR12 (MTMRC): 202-556 (human), 225-556 (mouse, rat, zebra fish),
+
* MTMR13 (MTMRD): 223-581 (human), 345-581 (xenopus)
+
  
The result shows we can build four individual HMMs for the subfamilies:
+
We found the alignment contains both SacN and the catalytic phosphatase domain. We split them based upon the ranges of the two domains of yeast SAC1 derived from the crystal structure study mentioned above.
* MTMR1, MTMR6, MTMR9 (perhaps MTMR3)
+
* MTMR10
+
* MTMR13
+
* MTMR14
+
 
+
Note: In addition to PSI-BLAST, we tried to obtain an alignment from HHPRED, but the quality is pretty bad.
+
  
 
=== References ===
 
=== References ===
 
<biblio>
 
<biblio>
 +
#Cai14 pmid=25113029
 +
#Manford10 pmid=20389282
 
</biblio>
 
</biblio>
 +
 +
=== Supplementary materials ===
 +
The full sequence of yeast SAC1:
 +
 +
>yeast_SAC1
 +
MTGPIVYVQNADGIFFKLAEGKGTNDAVIHLANQDQGVRVLGAEEFPVQGEVVKIASLMGFIKLKLNRYAIIANTVEETGRFNGHVFYRVLQHSIVSTKFNSRIDSEEAEYIKLLELHLKNSTFYFSYTYDLTNSLQRNEKVGPAASWKTADERFFWNHYLTEDLRNFAHQDPRIDSFIQPVIYGYAKTVDAVLNATPIVLGLITRRSIFRAGTRYFRRGVDKDGNVGNFNETEQILLAENPESEKIHVFSFLQTRGSVPIYWAEINNLKYKPNLVLGENSLDATKKHFDQQKELYGDNYLVNLVNQKGHELPVKEGYESVVHALNDPKIHYVYFDFHHECRKMQWHRVKLLIDHLEKLGLSNEDFFHKVIDSNGNTVEIVNEQHSVVRTNCMDCLDRTNVVQSVLAQWVLQKEFESADVVATGSTWEDNAPLLTSYQNLWADNADAVSVAYSGTGALKTDFTRTGKRTRLGAFNDFLNSASRYYQNNWTDGPRQDSYDLFLGGFRPHTASIKSPFPDRRPVYIQLIPMIICAALTVLGATIFFPKDRFTSSKNLLYFAGASIVLALSTKFMFKNGIQFVNWPKLVDVGFLVVHQTHDKEQQFKGLKYAQSPKFSKPDPLKRD
 +
 +
The sequence of SacN domain of SAC1:
 +
>yeast_SAC1_SacN
 +
MTGPIVYVQNADGIFFKLAEGKGTNDAVIHLANQDQGVRVLGAEEFPVQGEVVKIASLMGFIKLKLNRYAIIANTVEETGRFNGHVFYRVLQHSIVSTKFNSRIDSEEAEYIKLLELHLKNSTFYFSYTYDLTNSLQRNEKVGPAASWKTADERFFWNHYLTEDLRNFAHQDPRIDSFIQP
 +
 +
The sequence of catalytic domain of yeast SAC1:
 +
>yeast_SAC1_Catalytic_Domain
 +
VIYGYAKTVDAVLNATPIVLGLITRRSIFRAGTRYFRRGVDKDGNVGNFNETEQILLAENPESEKIHVFSFLQTRGSVPIYWAEINNLKYKPNLVLGENSLDATKKHFDQQKELYGDNYLVNLVNQKGHELPVKEGYESVVHALNDPKIHYVYFDFHHECRKMQWHRVKLLIDHLEKLGLSNEDFFHKVIDSNGNTVEIVNEQHSVVRTNCMDCLDRTNVVQSVLAQWVLQKEFESADVVATGSTWEDNAPLLTSYQNLWADNADAVSVAYSGTGALKTDFTRTGKRTRLGAFNDFLNSASRYYQNNWTDGPRQDSYDLFLGG

Latest revision as of 18:22, 26 March 2017

Back to List of HMMs

Symbol: CC1_Sac_PD

Name: Sac phosphatase domain

Description

Based upon the crystal structure of the yeast member, SAC1 has two structural domains: SacN and the catalytic phosphatase domain [1] (see the sequences in supplementary materials). The SacN ranges from 1-181; the catalytic phosphatase domain ranges from 182-504. The SacN domain mediates the interaction with VPS74, which is proposed to mediate packaging of medial Golgi glycosyltransferases into coatomer (COP1)-coated vesicles, thereby maintaining Golgi residence [2].

When using the HMM/PSSM profiles of public databases to determine the region of phosphatase domain, they gave different regions from the region determined from the first crystal structure:

  • Pfam database: Syja_N (PF02383), 56-343
  • CDD database: 1) Cdd:pfam02383, 55-348, 2) Cdd:COG5329, 3-578
  • SMART: no hit

Thus, Pfam profile contains part of SacN domain and part of the phosphatase domain, and it lacks the region of 451-511 is required for the PtdIns4P phosphatase activity of yeast SAC1 [2]. COG profile contains both SacN and the catalytic phosphatase domain, but also the trans-membrane–spanning domains (523-578) (figure 4A in [2]). We therefore built our own profiles for SacN and catalytic phosphatase domain.

How the HMM was built

We PSI-BLASTed the full sequence of yeast SAC1 (see supplementary materials) against SwissProt via NCBI BLAST server. The search converged right after 1st round. The query yeast SAC1 hit the SACs from animals and plants. We downloaded the aligned sequences, aligned them with Clustal Omega, manually adjusted the alignment, built the HMM and validated the HMM using the protein phosphatases of the nine genomes.

We found the alignment contains both SacN and the catalytic phosphatase domain. We split them based upon the ranges of the two domains of yeast SAC1 derived from the crystal structure study mentioned above.

References

  1. Manford A, Xia T, Saxena AK, Stefan C, Hu F, Emr SD, and Mao Y. Crystal structure of the yeast Sac1: implications for its phosphoinositide phosphatase function. EMBO J. 2010 May 5;29(9):1489-98. DOI:10.1038/emboj.2010.57 | PubMed ID:20389282 | HubMed [Manford10]
  2. Cai Y, Deng Y, Horenkamp F, Reinisch KM, and Burd CG. Sac1-Vps74 structure reveals a mechanism to terminate phosphoinositide signaling in the Golgi apparatus. J Cell Biol. 2014 Aug 18;206(4):485-91. DOI:10.1083/jcb.201404041 | PubMed ID:25113029 | HubMed [Cai14]
All Medline abstracts: PubMed | HubMed

Supplementary materials

The full sequence of yeast SAC1:

>yeast_SAC1 MTGPIVYVQNADGIFFKLAEGKGTNDAVIHLANQDQGVRVLGAEEFPVQGEVVKIASLMGFIKLKLNRYAIIANTVEETGRFNGHVFYRVLQHSIVSTKFNSRIDSEEAEYIKLLELHLKNSTFYFSYTYDLTNSLQRNEKVGPAASWKTADERFFWNHYLTEDLRNFAHQDPRIDSFIQPVIYGYAKTVDAVLNATPIVLGLITRRSIFRAGTRYFRRGVDKDGNVGNFNETEQILLAENPESEKIHVFSFLQTRGSVPIYWAEINNLKYKPNLVLGENSLDATKKHFDQQKELYGDNYLVNLVNQKGHELPVKEGYESVVHALNDPKIHYVYFDFHHECRKMQWHRVKLLIDHLEKLGLSNEDFFHKVIDSNGNTVEIVNEQHSVVRTNCMDCLDRTNVVQSVLAQWVLQKEFESADVVATGSTWEDNAPLLTSYQNLWADNADAVSVAYSGTGALKTDFTRTGKRTRLGAFNDFLNSASRYYQNNWTDGPRQDSYDLFLGGFRPHTASIKSPFPDRRPVYIQLIPMIICAALTVLGATIFFPKDRFTSSKNLLYFAGASIVLALSTKFMFKNGIQFVNWPKLVDVGFLVVHQTHDKEQQFKGLKYAQSPKFSKPDPLKRD

The sequence of SacN domain of SAC1: >yeast_SAC1_SacN MTGPIVYVQNADGIFFKLAEGKGTNDAVIHLANQDQGVRVLGAEEFPVQGEVVKIASLMGFIKLKLNRYAIIANTVEETGRFNGHVFYRVLQHSIVSTKFNSRIDSEEAEYIKLLELHLKNSTFYFSYTYDLTNSLQRNEKVGPAASWKTADERFFWNHYLTEDLRNFAHQDPRIDSFIQP

The sequence of catalytic domain of yeast SAC1: >yeast_SAC1_Catalytic_Domain VIYGYAKTVDAVLNATPIVLGLITRRSIFRAGTRYFRRGVDKDGNVGNFNETEQILLAENPESEKIHVFSFLQTRGSVPIYWAEINNLKYKPNLVLGENSLDATKKHFDQQKELYGDNYLVNLVNQKGHELPVKEGYESVVHALNDPKIHYVYFDFHHECRKMQWHRVKLLIDHLEKLGLSNEDFFHKVIDSNGNTVEIVNEQHSVVRTNCMDCLDRTNVVQSVLAQWVLQKEFESADVVATGSTWEDNAPLLTSYQNLWADNADAVSVAYSGTGALKTDFTRTGKRTRLGAFNDFLNSASRYYQNNWTDGPRQDSYDLFLGG