Difference between revisions of "HMM PD00142"

From PhosphataseWiki
Jump to: navigation, search
m (Mark moved page HMM PD0142 to HMM PD00142)
 
(5 intermediate revisions by one user not shown)
Line 6: Line 6:
  
 
=== Why built the HMM? ===
 
=== Why built the HMM? ===
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, which lacks the C-terminal ~45 aa.
 
  
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>:
+
In brief, there is no available HMM profile for human MTMR14. Read [[HMM_PD0138|CC1_MTM_1_PD]] for more information.  
 
+
>Hsap_MTMR2:205-586 id=HsapP048_AA CC1:CC1:Myotubularin:MTMR1 [Homo sapiens]
+
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:
+
* MTM1, MTMR1, MTMR2: ~190-600, 100%
+
* 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)
+
* MTMR14: no hit
+
 
+
We just built the HMM followed the [[HMM#Guidance_for_HMM_building|general guidance]]. By using the protein phosphatases of the nine genomes, we found the HMM profile can detect the phosphatase domains of all subfamilies with good boundaries except MTMR14. The HMM has a length of 391 aa, all subfamilies except MTMR14 has the domains of > 380 aa long. Because the HMM does not capture MTMR14, we archive the HMM in a separate page [[HMM_PD0138|CC1_MTM_1_PD]].
+
 
+
The [[HMM_PD0138|CC1_MTM_1_PD]] HMM profile hit a region of 60 aa in human MTMR14. Because there is not structure of MTMR14, we looked into the conserved region of MTMR14 and investigated whether the region can align with [[HMM_PD0138|CC1_MTM_1_PD]].
+
  
 
=== How was the HMM built? ===
 
=== How was the HMM built? ===
==== Conserved region of [[PD0143|MTMR14]] ====
+
==== Conserved region of [[HMM_PD00143|MTMR14]] ====
 
We PSI-BLASTed human MTMR14 full sequence against NR database via NCBI BLAST server. Eyeballing the sequences, the query hit MTMR3 and other myotubularins. To build a HMM for full sequence MTMR14, we chose MTMR14 sequences from description. We removed redundant sequences above the threshold of 70% identity. We then aligned the sequences using Clustal Omega, followed by manual adjustment. We built the HMM (3.1b) using default parameter. By searching against the protein phosphatases of nine genomes, we found it matched to ~70 aa of other myotubularin subfamilies.
 
We PSI-BLASTed human MTMR14 full sequence against NR database via NCBI BLAST server. Eyeballing the sequences, the query hit MTMR3 and other myotubularins. To build a HMM for full sequence MTMR14, we chose MTMR14 sequences from description. We removed redundant sequences above the threshold of 70% identity. We then aligned the sequences using Clustal Omega, followed by manual adjustment. We built the HMM (3.1b) using default parameter. By searching against the protein phosphatases of nine genomes, we found it matched to ~70 aa of other myotubularin subfamilies.
  
 
==== Comparison with [[HMM_PD0138|CC1_MTM_1_PD]] ====
 
==== Comparison with [[HMM_PD0138|CC1_MTM_1_PD]] ====
 
We aligned the alignments of CC1_MTM_1_PD and the conserved region of MTMR14. We found the conserved region of MTMR14 is a subset of CC1_MTM_1_PD.
 
We aligned the alignments of CC1_MTM_1_PD and the conserved region of MTMR14. We found the conserved region of MTMR14 is a subset of CC1_MTM_1_PD.
 +
 +
==== Obtain phosphatase domain ====
 +
The conserved region contains N-termimal GRAM domain and C-terminal phosphatase domain. We define the phosphatase domain as the region adjacent to GRAM domain and retrieve the alignment block. We then built the HMM profile.
  
  

Latest revision as of 17:49, 28 October 2015

Back to List of HMMs

Symbol: CC1_MTM_MTMR14_PD

Name: MTMR14 phosphatase domain

Why built the HMM?

In brief, there is no available HMM profile for human MTMR14. Read CC1_MTM_1_PD for more information.

How was the HMM built?

Conserved region of MTMR14

We PSI-BLASTed human MTMR14 full sequence against NR database via NCBI BLAST server. Eyeballing the sequences, the query hit MTMR3 and other myotubularins. To build a HMM for full sequence MTMR14, we chose MTMR14 sequences from description. We removed redundant sequences above the threshold of 70% identity. We then aligned the sequences using Clustal Omega, followed by manual adjustment. We built the HMM (3.1b) using default parameter. By searching against the protein phosphatases of nine genomes, we found it matched to ~70 aa of other myotubularin subfamilies.

Comparison with CC1_MTM_1_PD

We aligned the alignments of CC1_MTM_1_PD and the conserved region of MTMR14. We found the conserved region of MTMR14 is a subset of CC1_MTM_1_PD.

Obtain phosphatase domain

The conserved region contains N-termimal GRAM domain and C-terminal phosphatase domain. We define the phosphatase domain as the region adjacent to GRAM domain and retrieve the alignment block. We then built the HMM profile.


Note: In addition to PSI-BLAST, we tried to obtain an alignment from HHPRED, but the quality is pretty bad.

References

  1. Begley MJ, Taylor GS, Kim SA, Veine DM, Dixon JE, and Stuckey JA. Crystal structure of a phosphoinositide phosphatase, MTMR2: insights into myotubular myopathy and Charcot-Marie-Tooth syndrome. Mol Cell. 2003 Dec;12(6):1391-402. DOI:10.1016/s1097-2765(03)00486-6 | PubMed ID:14690594 | HubMed [Begley03]