ISPUB.com / IJGP/4/1/5941
  • Author/Editor Login
  • Registration
  • Facebook
  • Google Plus

ISPUB.com

Internet
Scientific
Publications

  • Home
  • Journals
  • Latest Articles
  • Disclaimers
  • Article Submissions
  • Contact
  • Help
  • The Internet Journal of Genomics and Proteomics
  • Volume 4
  • Number 1

Original Article

In Silico characterization of STK11- a protein involved in PJ Syndrome

M Reddy, B Bhanu, M Prasad, N Ramya, K Madhulika

Keywords

ampk, intestinal cancer, lkb1, peutz-jeghers syndrome, stk11, tumors

Citation

M Reddy, B Bhanu, M Prasad, N Ramya, K Madhulika. In Silico characterization of STK11- a protein involved in PJ Syndrome. The Internet Journal of Genomics and Proteomics. 2008 Volume 4 Number 1.

Abstract

Studies of hereditary cancer syndromes have contributed greatly to our understanding of molecular events involved in tumor genesis. One such syndrome of our interest is Peutz-Jeghers Syndrome, a hereditary disease, in which there is predisposition to benign and malignant tumors of many organ systems. Here, we investigate the genes responsible for the Peutz-Jeghers syndrome (PJS) and the pathway of disease causing mechanism. The STK11, Serine/Threonine Kinase 11, commonly known as LKB1 is the gene responsible for this syndrome. This gene is coding a protein of length 433 amino acids, 48.6 kDa. The protein sequence is retrieved from the protein database of NCBI. The analysis of this protein is performed using various in silico methods.

 

Introduction

Peutz-Jeghers syndrome is characterized by the development of growths called hamartomatous polyps in the gastrointestinal tract (particularly the stomach and intestines) [1]. Patients with this syndrome have a 15-fold increased risk of developing intestinal cancer compared to that of the general population [2,3]. Almost 50% of patients with Peutz-Jeghers syndrome develop and die from cancer around the age of 57 years. The cumulative risk for developing any cancers associated with Peutz-Jeghers syndrome in patients aged 15–64 years is 93% [4,5]. Peutz-Jeghers syndrome has been described in all races. The occurrence of cases in males and females is about equal. The average age at diagnosis is 23 years in men and 26 years in women. The gastrointestinal polyps found in this syndrome are typical hamartomas. Gastrointestinal polyps can result in chronic bleeding and anemia and cause recurrent obstruction and intussusceptions requiring repeated laparotomies and bowel resections. Mucocutaneous hyper pigmentation appears in childhood as dark blue to dark brown macules around the mouth, eyes, and nostrils, in the perianal area, and on the buccal mucosa [6,7]. The signaling pathway of the STK11 gene product currently is not identified, hence the mechanism of hamartomatous polyp formation and mucocutaneous pigmentation is not known. A locus for this condition (PJS) was assigned to chromosome 19p. The mutations of the STK11 gene mapped to chromosome 19p13.3 are responsible for PJS [8].

Chromosome: 19; Location: 19p13.3

Figure 1

The reference sequence has been curated by NCBI staff. The reference sequence was derived from AC011544.6 and AC004221.2. These sequences are submitted by DOE Joint Genome Institute, Stanford Human Genome Center & DOE Joint Genome Institute, Lawrence Livermore National Laboratory, Livermore, CA respectively.

Methodology

Identification of Mutations from the Uniprot database

The STK11 gene is a tumor suppressor gene, which means that it normally prevents cells from growing and dividing too rapidly or in an uncontrolled way. The mutations are found to occur in the domain region, disrupting its ability to restrain cell division. Thus, the mutations in the STK11 gene are responsible for the Peutz Jeghers Syndrome. These mutations are identified from the protein information given in Uniprot Database [9] listed in Table 1.

Sequence analysis of the STK11 protein

The sequence analysis is performed using blastp [10] from NCBI. Sequences with high similarity are structurally conserved and are functionally related. Hence, the similar protein is identified to analyze if it has any role in the syndrome. Similarity search is performed against the non redundant database taking the expect threshold as 1 and blosum matrix 80. The query protein is compared with various model organisms. The protein having good score, high identities & similarities and fewer gaps is selected from the results of the similarity search in all the model organisms. This resulted in the identification of STK11 proteins in all the organisms selected (Table 2) and also the AMPK (Activated AMP Protein Kinase) as the closest protein to STK11 (Table 3). The AMPK has high sequence similarity and a common domain (kinase) with STK11.

Phylogenetic Studies using Biology Workbench

Evolutionary studies are performed on the sequences taken from the similarity search results. A multiple sequence alignment is built and the dendrogram is constructed (Figure 1) by Clustal W [11] to find the evolutionary relationship between STK11 and AMPK. The conservations in these proteins are analyzed from the Multiple Sequence Alignment built using both the local (Figures 2&3) and the global alignments (Figure 4) by Texshade [12] and Boxshade respectively. The distance matrix (Figure 5) is constructed using ClustalDist [13] to confirm the results of the dendrogram.

Structural Comparison of STK11 & AMPK

The structures of STK11 and the AMPK are analyzed to find the structural similarity. The secondary structure elements, alpha helices, beta sheets and coils are predicted using the tool Hierarchical Neural Network (HNN) [14] from the Expasy proteomics server. The composition of the secondary structure elements are listed in table 5. The 3D structures are retrieved using comparative modeling (CPH Models) [15] and threading methods (3DPSSM [16] server & HHPred [17]). The same structural entry, 2H6D is obtained from both the principles for the two proteins. The cartoon view of the structure is given in Figure 6. Also, the other properties of the protein viz, the molecular weight, atomic composition, iso electric point, half-life, extinction co efficient, Gravy, amino acid composition, instability index and aliphatic index are retrieved using the Protparam [18] tool, the results of which are listed in the table 4.

As the studies reveal good similarity of AMPK with STK11, the KEGG [19] pathway (Figure 7) is obtained to gain more insights into the disease causing mechanism. It is observed that AMPK is involved in the mTOR signaling pathway and that it is activated by its phosphorylation through STK11. Therefore, the phosporylation sites of AMPK are predicted using the tool NetPhos [20]. Also, the Prosite [21] pattern (Table 6), domain (ProDom [22] (Table 7) and family (Pfam [23]) of the proteins (Tables 8 & 9) is retrieved. The protein localization and the solubility are identified from Psort [24] (Table 10) & Sosui [25] (Table 11) respectively.

Results

Figure 2
Table 1: List of Mutations from the Uniprot Database

The sequence of STK11 is retrieved from the protein database of NCBI. The accession number of the protein is NP_000446.1. The protein is 433aa long. The FASTA format of the protein is given below and the similarity search results are listed in Table 2 & 3.

FASTA FORMAT of the query sequence (STK11):

>gi|4507271|ref|NP_000446.1| serine/threonine protein kinase 11 [Homo sapiens]

MEVVDPQQLGMFTEGELMSVGMDTFIHRIDSTEVIYQPRRKRAKLIGKYLMGDLLGEGSYGKVKEVLDSETLCRRAVKILKKKKLRRIPNGEANVKKEIQLLRRLRHKNVIQLVDVLYNEEKQKMYMVMEYCVCGMQEMLDSVPEKRFPVCQAHGYFCQLIDGLEYLHSQGIVHKDIKPGNLLLTTGGTLKISDLGVAEALHPFAADDTCRTSQGSPAFQPPEIANGLDTFSGFKVDIWSAGVTLYNITTGLYPFEGDNIYKLFENIGKGSYAIPGDCGPPLSDLLKGMLEYEPAKRFSIRQIRQHSWFRKKHPPAEAPVPIPPSPDTKDRWRSMTVVPYLEDLHGADEDEDLFDIEDDIIYTQDFTVPGQVPEEEASHNGQRRGLPKAVCMNGTEAAQLSTKSRAEGRAPNPARKACSASSKIRRLSACKQQ

Figure 3
Table 2: STK11 Sequences – Results of Similarity Searches

Figure 4
Table 3: AMPK Sequences – Results of Similarity Searches

The similarity search result has given significant similarity between the query and the AMPK in all the model organisms listed above in the tables. The studies of the proteins in the results indicate that the STK11 protein and the AMPK are having a kinase domain in common.

Figure 5
Figure 1: Dendrogram (Rooted Tree) constructed using Clustal W

Figure 6
Figure 2: Multiple Sequence Alignment built on local alignment to identify the conserved regions (Pg 1).

Figure 7
Figure 3: Multiple Sequence Alignment built on local alignment to identify the conserved regions (Pg2).

Figure 8
Figure 4: Multiple Sequence Alignment built on global alignment to identify the conserved regions

Violet and green colors are observed dominantly in the results of TEXSHADE and BOXSHADE respectively. These colors indicate the residues conserved in all the species taken.

Figure 9
Figure 5: Clustal Distance Matrix

Figure 10
Table 4: Physico-chemical properties of STK11 & AMPK

# Extinction coefficients are in units of M-1 cm-1, at 280 nm measured in water.

* Abs 0.1% (=1 g/l) 0.812, assuming ALL Cys residues appear as half cystines

** Abs 0.1% (=1 g/l) 0.799, assuming NO Cys residues appear as half cystines

~ The N-terminal of the sequence considered is M (Met). The estimated half-life is: 30 hours (mammalian reticulocytes, in vitro), >20 hours (yeast, in vivo), >10 hours (Escherichia coli, in vivo).

` The instability index provides an estimate of the stability of a protein in a test tube. Statistical analysis of 12 unstable and 32 stable proteins has revealed that there are certain dipeptides, the occurrence of which is significantly different in the unstable proteins compared with those in the stable ones. The authors of this method have assigned a weight value of instability to each of the 400 different dipeptides (DIWV). Using these weight values it is possible to compute an instability index (II) which is defined as:

Figure 11

Where: L is the length of sequence

DIWV(xx [i+1]) is the instability weight value for the dipeptide starting in position i.

A protein whose instability index is smaller than 40 is predicted as stable, a value above 40 predicts that the protein may be unstable.

Both these proteins exhibit similar physico-chemical properties, the properties that are similar are highlighted in the table above.

Figure 12

Figure 13
Table 5: %Composition of the Secondary Structure Elements in STK11 & AMPK

Figure 14
Figure 6: The 3d structure of the PDB entry 2H6D – Rasmol View

The cartoon view of the secondary structural elements, alpha helices, beta sheets & coils are colored red, green & yellow respectively. The amino acid Leucine contributing the highest percent in composition is highlighted using the ball& stick rasmol view (orange color).

Figure 15
Table 6: Prosite Pattern of STK11 & AMPK

The PROSITE pattern observed in the results of these proteins is the signature pattern of the Serine – Threonine protein kinase family. The kinase domain is shared in common for both the proteins. This confirms the hypothesis made in the similarity search.

Figure 16
Table 7: Domain identification using ProDom

Figure 17
Table 8: Pfam Results of AMPK

Figure 18
Table 9: Pfam Result for STK11

Both these proteins belong to a family of PKinase. However, the protein AMPK is a member of many other families as seen in the result. This is because this protein is involved in many metabolic pathways.

Figure 19
Table 10: Identification of Protein Localization using Psort

These results confirm that these are mostly cytoplasmic proteins.

Figure 20
Table 11: Confirmation of the protein's solubility from Sosui

Figure 21
Figure 7: KEGG Metabolic Pathway

Figure 22

Discussion

The Insilico analysis of the STK11 reveals the role of AMPK in the Peutz-Jeghers syndrome. AMP-activated protein kinase (AMPK) is a primary regulator of the cellular response to lowered ATP levels in eukaryotic cells [26,27]. AMPK is a serine/threonine protein kinase [26].

The pathway also suggests involvement of both these proteins in the same pathway of cell growth. However, AMPK is involved in many other metabolic processes. AMPK acts as a metabolic master switch regulating several intracellular systems including the cellular uptake of glucose, the β-oxidation of fatty acids and the biogenesis of glucose transporter 4 (GLUT4). STK11 has been shown to phosphorylate AMPK on Thr-172 and serves as an essential component of physiological AMPK activation. The identified human mutations of STK11 occur in the catalytic domain and cause a loss of its kinase activity. This failure of STK11 phosphorylation impairs its downstream signaling of AMPK [28], thus disturbing the normal cycle of cell growth to abnormal growth leading to the formation of hamartomatous, cancerous polyps, i.e. the Peutz-Jeghers Syndrome. The role of STK11/AMPK in the survival of a subset of genetically defined tumor cells may provide opportunities for cancer therapeutics.

Acknowledgements

We would like to thank Mr. G. Ashok Kumar, Director, Ventura Institute of Biosciences, Hyderabad, for providing necessary facilities during the period of the project work.

References

1. Peter Beighton & Gretha Beighton. The Man behind the Syndrome. Springer-Verlag, Berlin-Heidelberg, 1986.
2. Jeghers, H, Mckusick, VA, KATZ, KH.Generalized intestinal polyposis and melanin spots of the oral mucosa, lips and digits; a syndrome of diagnostic significance. N Engl J Med 1949; 241:993.
3. Peutz, JL. Over een zeer merkwaardige, gecombineerde familiaire pollyposis van de sligmliezen van den tractus intestinalis met die van de neuskeelholte en gepaard met eigenaardige pigmentaties van huid-en slijmvliezen. Ned Maandschr v Gen 1921; 10:134.
4. Giardiello FM, Brensinger JD, Tersmette AC, Goodman SN, Petersen GM, Booker SV, Cruz-Correa M, Offerhaus JA (2000) Very high risk of cancer in familial Peutz-Jeghers syndrome. Gastroenterology Med 119:1447-53 [Medline]
5. Giardiello FM, Welsh SB, Hamilton SR, Offerhaus GJ, Gittelsohn AM, Booker SV, Krush AJ, Yardley JH, Luk GD (1987) Increased risk of cancer in the Peutz-Jeghers syndrome. N Engl J Med 316:1511-4 [Medline]
6. Bruwer A, Bargen Ja, Kierland Rr (1954) Surface pigmentation and generalized intestinal polyposis; (Peutz-Jeghers syndrome). Mayo Clin Proc 29:168-71 [Medline]
7. Burke CA, Santisi J, Church J, Levinthal G (2005) the utility of capsule endoscopy small bowel surveillance in patients with polyposis. Am J Gastroenterol 100:1498-502 [Medline]
8. Jenne DE, Reimann H, Nezu J, Friedel W, Loff S, Jeschke R, Müller O, Back W, Zimmer M. Peutz-Jeghers syndrome is caused by mutations in a novel serine threonine kinase. Nature Genetics. 1998 Jan; 18(1):38-43.
9. Bairoch A., Apweiler R., Wu C.H., Barker W.C., Boeckmann B., Ferro S.,Gasteiger E., Huang H., Lopez R., Magrane M., Martin M.J.,Natale D.A., O'Donovan C., Redaschi N., YehL.S. The Universal Protein Resource (UniProt), Nucleic Acids Res. 33:D154-D159 (2005).
10. Altschul, S.F., Gish, W., Miller, W., Myers, E.W. & Lipman, D.J. (1990) "Basic local alignment search tool." J. Mol. Biol. 215:403-410. PubMed
11. Thompson J.D., Higgins D.G., Gibson T.J. "CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice." Nucleic Acids Res. 22:4673-4680(1994).
12. Beitz, E. (2000), TeXshade: shading and labeling of multiple sequence alignments using LaTeX2e. Bioinformatics 16: 135-139.
13. Thompson J.D., Higgins D.G., Gibson T.J. "CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice." Nucleic Acids Res. 22:4673-4680(1994).
14. Combinaison de classifieurs statistiques, Application a la prediction de structure secondaire des proteins, PhD Thesis , Guermeur, Y
15. O. Lund, M. Nielsen, C. Lundegaard, P. Worning, CPHmodels 2.0: X3M a Computer Program to extract 3D Models. Abstract at the CASP5 conferenceA102, 2002.
16. http://www.sbg.bio.ic.ac.uk/~3dpssm/index2.html
17. Söding J, Biegert A, and Lupas AN. (2005) The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Research 33, W244--W248 (Web Server issue).doi:10.1093/nar/gki40.
18. Gasteiger E., Hoogland C., Gattiker A., Duvaud S., Wilkins M.R., Appel R.D., Bairoch A.; Protein Identification and Analysis Tools on the ExPASy Server;(In) John M. Walker (ed): The Proteomics Protocols Handbook, Humana Press (2005) pp. 571-607
19. http://www.genome.jp/kegg/pathway.html
20. Blom, N., Gammeltoft, S., and Brunak, S. Sequence- and structure-based prediction of eukaryotic protein phosphorylation sites. Journal of Molecular Biology: 294(5): 1351-1362, 1999.
21. De Castro E., Sigrist C.J.A., Gattiker A., Bulliard V., Langendijk-Genevaux P.S., Gasteiger E., Bairoch A., Hulo N. ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006 Jul 1; 34(Web Server issue):W362-5.
22. Catherine Bru, Emmanuel Courcelle, Sébastien Carrère, Yoann Beausse, Sandrine Dalmar, and Daniel Kahn (2005) The ProDom database of protein domain families: more emphasis on 3D. Nucleic Acids Res. 33: D212-D215
23. Robert D. Finn, Jaina Mistry, Benjamin Schuster-Böckler, Sam Griffiths-Jones, Volker Hollich, Timo Lassmann, Simon Moxon, Mhairi Marshall, Ajay Khanna, Richard Durbin, Sean R. Eddy, Erik L. L. Sonnhammer and Alex Bateman, Pfam: clans, web tools and services. Nucleic Acids Research (2006) Database Issue 34:D247-D251
24. Paul Horton, Keun-Joon Park, Takeshi Obayashi, Naoya Fujita, Hajime Harada, C.J. Adams-Collier, & Kenta Nakai, "WoLF PSORT: Protein Localization Predictor", Nucleic Acids Research, doi:10.1093/nar/gkm259, 2007.
25. Hirokawa T., Boon-Chieng S., and Mitaku S., SOSUI: classification and secondary structure prediction system for membrane proteins. Bioinformatics, 14 378-9 (1998)
26. Hardie DG, Scott JW, Pan DA, Hudson ER: Management of cellular energy by the AMP-activated protein kinase system. FEBS Lett 2003, 546:113-20.
27. Kemp BE, Stapleton D, Campbell DJ, Chen ZP, Murthy S, Walter M, Gupta A, Adams JJ, Katsis F, van Denderen B, Jennings IG, Iseli T, Michell BJ, Witters LA: AMP-activated protein kinase, super metabolic regulator. Biochem Soc Trans 2003, 31:162-8.
28. Ojuka EO. Role of calcium and AMP kinase in the regulation of mitochondrial biogenesis and GLUT4 levels in muscle. Proc Nutr Soc. 63: 275-278, 2004.

Author Information

M. Vishwanath Reddy, Ph.D.
Department of Bioinformatics, Ventura Institute of Biosciences

B. Divya Bhanu, M.Sc
Department of Bioinformatics, Ventura Institute of Biosciences

MPSNR Prasad, M.tech
Department of Bioinformatics, Ventura Institute of Biosciences

N. Ramya, B.Tech
Department of Bioinformatics, Ventura Institute of Biosciences

K. Madhulika, B.Tech
Department of Bioinformatics, Ventura Institute of Biosciences

Download PDF

Your free access to ISPUB is funded by the following advertisements:

 

 

BACK TO TOP
  • Facebook
  • Google Plus

© 2013 Internet Scientific Publications, LLC. All rights reserved.    UBM Medica Network Privacy Policy

Close

Enter the site

Login

Password

Remember me

Forgot password?

Login

SIGN IN AS A USER

Use your account on the social network Facebook, to create a profile on BusinessPress