In silico mutational study reveal improved interaction between Beta-Hexosaminidase A and GM2 activator essential for the breakdown of GM2 and GA2 Gangliosides on Tay-Sachs disease
Kritika Kapruwan, Rahul Parida, Rajiniraja Muniyan*
School of Bio-Sciences and Technology, VIT University, Vellore-632014, Tamil Nadu, India
*Corresponding Author E-mail: rajiniraja.m@vit.ac.in
ABSTRACT:
Gene Hex-A on chromosome 15 in human provides instructions for making subunit of an enzyme called Beta-Hexosaminidase A (BHA).Within lysosomes BHA forms part of a complex that breaks down a fatty substance called GM2 gangliosides. The mutation in BHA of Hex-A gene is unable to degrade GM2 gangliosides, causes Tay-Sach’s disease.GM2 Activator (GM2A) a large binding pocket binds gangliosides and stimulates the breakdown of GM2 and GA2 gangliosides by BHA. The GSEP283 alpha loop of BHAis involved in GM2A binding. The aim of the present study is to represent and compare the interaction between BHA and GM2A proteins before and after mutation. The mutation Ser279Pro was performed using tool PyMOL and optimized using Arguslab with UFF force field. We used HADDOCK 2.2 server for protein-protein docking, SPDBV and PyMOL viewer software were used for understanding and analyzing the interaction. There were slight changes in the amino acid interactions at GSEP283 alpha loop region after mutation with decreasing binding energy when compared with unmutated protein. The docked results provided an evidence that the GSEP283 loop of BHA was binding to GM2A and we have represented for the first time a model on how the mutation on Hex-A affects GSEP283 interaction with GM2A activator. Our study also highlights that the TYR137 of GM2A is present in the surface that play role in lipid transfer. This in-silico study is a contribution for the research related to the interaction between BHA and GM2A proteins.
KEYWORDS: A mutation, GM2 Activator, Haddock, protein-protein interaction, GSEP283 alpha loop, Beta-Hexosaminidase A.
INTRODUCTION:
Lysosomes organelles which have acidified cytoplasm consist of degradative or hydrolytic enzymes specific for the degradation of cellular components1. The group of enzymes which are involved in the degradation of terminal, non-reducing sugar units from oligosaccharides molecules are exo-glycosidases. Gene Hex-A on chromosome 15 in human provides instructions for making subunit of an enzyme called Beta-Hexosaminidase A (BHA).
The BHA are lysosomal isozymes that separate the glycosidic linkages of the non-reducing, terminal L-DN-acetylglucosamine (GlcNAc) or L-D-N-acetyl galactosamine (GalNAc) residues on glycolipids, glycoproteins, also glycosaminoglycans. The Activator is a substrate specific (GM2 ganglioside) co-factor for one of the hexosaminidase isozymes, hexosaminidase - A. The presence of an N-acetyl-L-D-glucosaminidase (NAG) was first reported in 19362.
Specifically the protein produced from the Hex-A gene forms the alpha subunit of the enzymes BHA3. One alpha subunit interacts with one beta subunit produced from Hex-B gene to form a functioning enzyme. This complex formed interacts with GM2 Activator to prevent the accumulation of GM2 gangliosides1.Tay-Sachs is autosomal recessive genetic disorder which is caused by mutation in a vital enzyme BHA. Classical Tay-Sachs patients usually do not survive beyond 4 years of age4. There were many mutations listed out of which we chose the Ser279 to Pro2795. Tay-sach’s is neurodegenerative lysosomal storage diseases associated with excessive accumulation of GM26. The history of TaySach’s diseases was reported are listed(Table 1). Due to Hex-A genedysfunctioning there is accumulation of a fatty substance called GM2 gangliosides especially in the nerve cells of the brain. This progressive accumulation of GM2 gangliosides leads to damage in the nerve cells1.GM2A is an activator that initiates degradation of GM2 gangliosides3. It acts as substrate specific co-factor for BHA as it is a large binding pocket that can hold single chain phospholipids and fatty acids. It has TYR137 on the surface that plays a major role in lipid transfer from the membrane to the BHA7. It works by binding gangliosides and stimulates GM2 degradation by calcium independent phospholipase activity4. It only stimulates breakdown of GM2 and GA2 by the beta Hex-A. It extracts single GM2 units from membrane which is presented in soluble form to BHA for cleavage of N-acteyl D-galactosamine and conversion to GM3.Since the major interaction is between GM2A and BHA, the GSEP283 loop(alpha loop) in Hex-A protein is reported as to be responsible in binding to GM2A3.
Table 1.TaySach’s history
|
Year |
Event |
Ref. |
|
1881 |
Warren Tay: initial description of clinical characteristics, cherry red spot |
20 |
|
1896 |
Bernard Sachs: `Amaurotic family idiocy' prevalent in the Jews |
21 |
|
1962 |
GM2 ganglioside identi¢ed as the compound stored in TSD |
22 |
|
1968 |
Demonstration of two major Hex isozymes in normal human tissue (Hex A and Hex B) |
10 |
|
1969 |
Hex A, but not Hex B is de¢cient in Jewish TSD |
12 |
|
1969 |
Isozyme patterns in non-Jewish TSD patients: (a) deficiency of both Hex A and B (SD); (b) no Hex deficiency AB-variant and B1-variant |
13 |
|
1973 |
Immunological basis for the subunit compositions of Hex A (KL) and B (LL): deficiency of the K results in TSD, deficiency of L results in SD |
23 |
|
1979 |
Purification and characterization of the GM2 Activator protein: its deficiency results in the AB-variant form |
24 |
|
1984-1985 |
Isolation of cDNA clones encoding the K- and L-subunits |
15,25 |
|
1986 |
Identification of partial HEXA and HEXB gene deletions resulting in infantile TSD and SD |
14,17 |
|
1988 |
Identification of two major HEXA mutations responsible for TSD in Ashkenazi Jews |
26-28 |
|
1989-1992 |
Isolation of a cDNA clone encoding the GM2 activator |
16 |
|
1991-1992 |
Association of a GM2A mutation with the AB-variant form |
18,19 |
GM2 is an intermediate in both the amalgamation (ER) also, degradation (lysosome) of the higher brain gangliosides, e.g. GM1. So as to hydrolyze GM2, Hex A needs the little, monomeric, heat stable GM2 activator protein. The Activator connects with both the carbohydrate and lipid part of the ganglioside and it takes out steric hindrance from the membrane by solubilizing a ganglioside. Also, Activator binding may interfere with a specific hydrogen bond between the acetamido-NH of the GalNAc residue and the carboxylic group of the NeuAc residue of GM2, making the GalNAc accessible for hydrolysis by Hex A8. The relative significance of these two factions is right now a matter of level headed discussion in the writing9.
Our main focus is to provide evidence that GSEP283 alpha loop of BHA interacts with GM2 Activator for the degradation of GM2 gangliosides, the amino acids involved in the interaction and comparative study of binding energy before and after mutation. We also provide our own hypothesis related to functioning of GM2 Activator with BHA.
MATERIALS AND METHODS:
The software and server used are Arguslab 4.0.1, PyMOL, and SPDBV.HADDOCK2.2: High Ambiguity Driven protein-protein DOCKing- is information operated flexible docking method for the modelling of biomolecular complexes. This docking online application differs itself from ab-initio docking methods encodes information from identified or predicted protein interfaces in ambiguous interaction restraints to drive the docking process. HADDOCK can deal with many classes of modelling problems such as protein-protein, protein-nucleic acids and protein-ligand complexes11.
The protein3Dstructures were obtained from Protein Databank (PDB). BHA protein (PDB ID:2GK1_A) and GM2A (PDB ID:2AF9) were taken as receptor and ligand respectively. We preceded our studies in two ways, firstly the unmutated A chain of 2GK1 was docked with 2AF9 by using HADDOCK2.2 server. Secondly, the 2GK1_A was mutated (S279P) using SPDBV and PyMol. This mutated Hex-A protein (2GK1_A) was then docked with GM2A (2AF9) and the amino acid interactions were studied. The docked clusters were then viewed in PyMol and SPDBV.
RESULTS AND DISCUSSION:
The HADDOCK results provided the best structures or solution. The Z score indicates how many deviations from the average a particular cluster is located in terms of score (the more the negative the better). The lower the energy is better.
Table 2.Energy profile of before and after mutation (Ser279Pro) of BHA.
|
Energies/parameters |
Binding energy before mutation (kcal/Mol) |
Binding energy after mutation (kcal/Mol) |
|
Z score |
-1.2 |
-1.5 |
|
RMSD from the overall lowest energy structure |
7.3 +/-0.1 |
1.2 +/-0.8 |
|
Van der Waals energy |
-29.4 +/-4.4 |
-34.8 +/-4.9 |
|
Electrostatic energy |
-127.4 +/-26.5 |
-184.1 +/-6.3 |
|
Restraints violation energy |
58.4 +/- 3.00 |
39.5 +/-10.53 |
|
Haddock score |
-66.6 +/-2.7 |
-89.9 +/- 2.4 |
Figure 1.The graphical representation of the energy before and after mutation
Unmutated A chain of BHA (2GK1_A) docking with GM2A (2AF9):
HADDOCK clustered 397 structures in3 cluster(s) out of which top 10 were shown, which represents 99.25% of water-refined models generated. Figure 1 shows the cluster 2 was the best solution with Z score of -1.2 out of 3 clusters. The electrostatic and Vander Waals energy were recorded as -127.4 and -29.4 respectively and the lowest RMSD value 7.3 (Table 2).
The unmutated A chain of Hex-A (2GK1_A) interaction with the GM2A (2AF9) shows that, the GSEP283 alpha loop of 2GK1_A is interacting with 2AF9. We also found the amino acids Met0, Glu58, Asn136 and Ala60 of 2AF9 were interacting with the 2GK1_A GSEP283 loop (Figure 2). The TYR137 of GM2A (2AF9) plays a major role in lipid transfer from the membrane to the Hex-A (2GK1_A), confirmed from the result obtained and also highlighted that its present on the surface.
Figure2. Partial GSEP283 loop(right) interaction with 2AF9 (left)
Mutated A chain of BHA (2GK1_A) docking with GM2A (2AF9)
HADDOCK clustered 379 structures in 4 cluster(s) out of which top 10 were shown, which represents 94.75% of water-refined models generated. Figure 1 shows that the cluster 1 was the best solution with Z score of -1.5 out of 4 clusters. The electrostatic and Vander Waals energies were-184.1 and -34.8 respectively (Table 2).
Figure 3.Partial GSEP283 loop (right) interaction with 2AF9 (left)
The mutated 2GK1_A interaction with the 2AF9 when viewed in SPDBV and PyMol we observed that the GSEP283 loop is partially interacting with the 2AF9 (the GM2A) as compare to the unmutated 2GK1_A docking results. We also analysed that this docking result are better than the unmutated 2GK1_A docking with 2AF9 as the Z score is -1.5. We found that the Gly280 and Ser281 of GSEP283 loop are the only one interacting with 2AF9 amino acids Asn136 and Glu58 and Tyr137 respectively (Figure 3). The Tyr137 interaction with the Ser281 is a weak hydrogen bond.
CONCLUSION:
This study is the evidence to research done earlier reporting that GSEP283 alpha loop of Hex-A protein binds to the GM2Activator. From the above results we can hypothesise that Glu282 of GSEP283 alpha loop plays a crucial role in the functioning of Hex-A and GM2Actiavtor in degradation of gangliosides and that the Tyr137 is restricted after mutation to play its role in lipid transfer since it is interacting with Ser281(BHA) and results in TaySach’s disease. This in silico model also demonstrates for the first time that the GSEP283 loop of Hex-A protein interact with Met0, Asn136, Glu58, Ala60 of GM2Activator. Our study opens the gate for further research regarding mutation, interactions and functioning involved in Tay-sach’s disease. It can also help in the study of drug related to TaySach’s disease
REFERENCES:
1. Mahuran DJ. Biochemical consequences of mutations causing the GM2 gangliosidoses. Biochimicaet Biophysica Acta (BBA) - Molecular Basis of Disease. 1455(2-3); 1999: 105-138.
2. Watanabe K. Biochemical Studies On Carbohydrates xxii. On Animal β-N-Monoacetylglucosaminidase. Preliminary Report. J. Biochem. 24; 1936: 297.
3. Zarghoonia M, Bukovaca S, Tropaka M, Callahana J, and Mahuran DJ. An α-subunit loop structure is required for GM2 activator protein binding by β-hexosaminidase A ,Research Institute, The Hospital for Sick Children, Toronto, Ont., Canada. 324(3); 2004:1048-1052
4. GravelRA, Kaback MM, Proia RL et al. The Metabolic and Molecular Bases of Inherited Disease, McGraw-Hill, New York 2(7); 1995: 2839-2879.
5. Drucker L, Hemli JA, Navon R. Two mutated HEXA alleles in a Druze patient with late-infantile Tay-Sachs disease. Hum Mutat. 10; 1997: 451-457.
6. KitakazeK, Mizutani Y, Sugiyama E et al. Protease-resistant modified human β-hexosaminidase B amelioratessymptoms in GM2 gangliosidosis model. Journal of Clinical Investigation.126(5);2016:1691–1703
7. Wright CS, Mi LZ, Lee S, et al. Crystal Structure Analysis of Phosphatidylcholine-GM2-Activator Product Complexes: Evidence for Hydrolase Activity. Biochemistry. 44(41); 2005:13510–13521
8. Wu YY, Lockyer JM, Sugiyama E, et al. Expression and specificity of human GM2 activator protein. Biol. Chem. 269;1994:16276-16283.
9. Mahuran DJ. The GM2 activator protein, its roles as a co-factor in GM2 hydrolysis and as a general glycolipid transport protein, Biochimica et BiophysicaActa.1393(1); 1998: 1-18.
10. Robinson D and Stirling JL. N-Acetyl-β-glucosaminidases in human spleen. Biochem. J. 107;1968: 321-327.
11. Dominguez C, Boelens R, Alexandre MJJB,HADDOCK: A Protein−Protein Docking Approach Based on Biochemical or Biophysical Information. Contribution from the Department of NMR Spectroscopy, BijvoetCenter for Biomolecular Research, Utrecht University. 125 (7); 2003:1731–1737
12. Okada S and O'Brien JS.Tay-Sachs Disease: Generalized Absence of a Beta-D-N-Acetylhexosaminidase Component Science 165;1969: 698-700.
13. Sandhoff K, Variation of β-N-acetylhexosaminidase-pattern in Tay-Sachs disease. FEBS Lett. 4;1969: 351-354.
14. O'Dowd BF, Klavins MH, Willard HF,et al. Molecular heterogeneity in the infantile and juvenile forms of Sandhoff disease (O-variant GM2 gangliosidosis). J. Biol. Chem. 261;1986: 12680-12685.
15. O'Dowd BF, Quan F, Willard HF et al. Isolation of cDNA clones coding for the beta subunit of human beta-hexosaminidase. Proc. Natl. Acad. Sci. USA 82;1985: 1184-1188.
16. Nagarajan S, Chen HC, Li SC et al. Evidence for two cDNA clones encoding human GM2-activator protein. Biochem. J. 282;1992: 807-813.
17. Myerowitz R andHogikyan ND. Different mutations in Ashkenazi Jewish and non-Jewish French Canadians with Tay-Sachs disease. Science 232;1986: 1646-1648.
18. Schröder M, Schnabel D, Suzuki K et al. A mutation in the gene of a glycolipid-binding protein (GM2 activator) that causes GM2-gangliosidosis variant AB. FEBS Lett. 290;1991: 1-3.
19. Xie B, Wang W, Mahuran DJ. A Cys138-to-Arg substitution in the GM2 activator protein is associated with the AB variant form of GM2 gangliosidosis.Am J Hum Genet. 50; 1992: 1046-1052.
20. Tay W. Symmetrical changes in the region of the yellow spot in each eye op an infant Trans. Ophthalmol. Soc. UK 1;1881:155-157.
21. Sachs B.A Family Form of Idiocy: Generally Fatal, and Associated with Early Blindness (amaurotic Family Idiocy) Dis. J. Nerv. Ment.21; 1896: 475-479.
22. Svennerholm L. The chemical structure of normal human brain and Tay-Sachs gangliosides. Biochem. Biophys. Res. Commun. 9;1962:436.
23. Srivastava SK and Beutler E. Hexosaminidase-A and Hexosamini-dase-B: Studies in Tay-Sachs' and Sandhoff's Disease. Nature 241;1973 :463-463.
24. Conzelmann E and S and hoff K. Purification and characterization of an activator protein for the degradation of glycolipids GM2 and GA2 by hexosaminidase A. Physiol. Chem.360;1979:1837-1849.
25. Myerowitz R and Proia RL. cDNA clone for the alpha-chain of human beta-hexosaminidase: deficiency of alpha-chain mRNA in Ashkenazi Tay-Sachs fibroblasts. Proc. Natl. Acad. Sci. USA 81;1984: 5394-5398.
26. Myerowitz R and Costigan FC. The major defect in Ashkenazi Jews with Tay-Sachs disease is an insertion in the gene for the alpha-chain of beta-hexosaminidase. J. Biol. Chem. 263;1988:18587-18589.
27. Myerowitz R. Splice junction mutation in some Ashkenazi Jews with Tay-Sachs disease: evidence against a single defect within this ethnic group. Proc. Natl. Acad. Sci. USA 85;1988:3955-3959.
28. Arpaia E, Dumbrille-Ross A, Maler T et al. Identification of an altered splice site in Ashkenazi Tay-Sachs disease. Nature 333;1988:85-86.
Received on 11.06.2017 Modified on 07.07.2017
Accepted on 20.09.2017 © RJPT All right reserved
Research J. Pharm. and Tech 2017; 10(11): 3899-3902.
DOI: 10.5958/0974-360X.2017.00708.9