3.1.1 Figure 1 and Table 1 shows the results of the Proteomic characteristics of TP53 and functional partners genes associated with Hepato carcinogenesis and liver cancer patients. The circos visualization of the results revealed a great diversity in proteomic properties among TP53 gene and the other functional partner genes. The dimensions of each strand linkage is proportional to the weight and measure of the features under consideration. The molecular weight (MW) for some of the functional partners genes in green colour were very high and consequently represented and linked with very large strands to ATM, CREBBP and EP300. When compared to other genes with small strands, it implies that the molecular weights for these genes were high and least in the gens with small strands. The extinction coefficient of the gene when exposed to UVR property and the hydropathicity attributes of the functional partner genes presented the least measurements and consequently represented with the smallest black-coloured strands. Hence the size of the strand per linkage is proportional to size and impact on the subject.
Table 1 presents the results of proteomic properties and figure 2 presents the results of interactive intensity of TP53 gene and associated functional partner genes in Hepato carcinoma and liver cancer patients. The number of amino acids varied from 3065 aa in ATM to 248 aa in SFN. The molecular weight of the genes was proportional to the number of amino acids present in the gene. ATM showed the highest molecular weight of 350,687 mg/mol while SFN had the least molecular weight of 27774 mg/mol. The isoelectric and electrode potentials of the genes varied from moderately acidic (4.64) to low alkalinity (8.83). The total number of negatively charged amino acid residues (Asp + Glu) ranges from 48 in SFN gene to 389 in ATM functional gene partner. The total number of positively charged amino acids residues (Arg + Lys) varied from 30 to 366 in SFN. The number of atoms in the gene molecules range from only 3864 in SFN to 49436 in ATM. Extinction coefficients among the genes also varied from 0.315 g/L in DAXX to 1.02 g/L in ATM. The estimated half-life of the TP53 and associated genes were 30 hrs when estimated using the human recticulyte cell. All TP53 and associated genes showed genetic and thermal instability except for RPA1 gene which was thermally and genetically stable with instability index of 35.83%. All the evaluated genes showed high aliphatic indices and ranged from 59.74% in HSP90AA1 to 94.61% in SFN gene. The solvent affinity index (hydropathicity) results revealed that TP53 and all the associated genes were hydrophilic with positive values of GRAVY ranging from 0.362 in RPA1 gene to 0.786 in TP53 gene. As revealed in the interactive heatmap, the intensity of the proteomic properties among the genes was more in ATM (deep brown colour) gene and least intense in DAXX gene (Deep blue colour) as shown in Figure 2.
Table 1: Proteomic characteristics of TP53 gene and functional partner genes associated with hepatocellular carcinoma and liver cancer patients.
S/N
|
Proteomic Characteristics
|
DAXX
|
PRA1
|
SFN
|
ATM
|
HS90AA1
|
EP300
|
CREBBP
|
SIRT1
|
MDM2
|
TP53
|
1
|
No. of Amino acid
|
740
|
616
|
248
|
3065
|
865
|
2414
|
2442
|
747
|
479
|
1134
|
2
|
Molecular weight
|
81372
|
68138
|
27774
|
350687
|
102343
|
201454
|
232655
|
82197
|
56328
|
132654
|
3
|
Theoretical pl.
|
4.79
|
6.92
|
4.68
|
6.39
|
5.07
|
8.81
|
8.83
|
4.55
|
4.64
|
5.88
|
4
|
Negative charged amino acid residue
|
135
|
71
|
48
|
389
|
160
|
178
|
183
|
134
|
52
|
146
|
5
|
Positively charged amino acid residue
|
87
|
71
|
30
|
366
|
122
|
208
|
214
|
75
|
54
|
132
|
6
|
Total no. of atoms
|
11264
|
9536
|
3864
|
49436
|
13769
|
36568
|
36792
|
11314
|
7682
|
17652
|
7
|
Extinction coefficient
|
0.315
|
0.931
|
0.937
|
1.02
|
0.824
|
0.561
|
0.543
|
0.575
|
0.765
|
0.655
|
8
|
Estimated half life
|
30
|
30
|
30
|
30
|
30
|
30
|
30
|
30
|
30
|
30
|
9
|
Instability index
|
73.37
|
35.83
|
48.59
|
45.21
|
47.71
|
67.04
|
65.91
|
53.88
|
58.71
|
58.22
|
10
|
Aliphatic index
|
69.14
|
78.96
|
94.61
|
78.09
|
59.74
|
60.95
|
78.41
|
65.68
|
70.75
|
78.17
|
11
|
Hydropathicity
|
0.848
|
0.362
|
0.599
|
0.247
|
0.699
|
0.725
|
0.69
|
0.551
|
0.786
|
0.76
|
Biological functions of Network functional partner genes of TP53 Associated With Hepatocellular Carcinoma and Liver Cancer
TP53
TP53 (Fig. 3 and Fig. 3b) is cellular tumor antigen p53; the gene acts as a tumor suppressor in many tumor types including hepato cellular carcinoma and liver cirrhosis. TP53 gene induces growth arrest or apoptosis depending on the physiological circumstances and cell type. The gene is also involved in cell cycle regulation as a trans-activator that acts to negatively regulate cell division by controlling a set of functional partner genes required for this process. One of the activated genes is an inhibitor of cyclin-dependent kinases. Apoptosis induction seems to be mediated either by stimulation of BAX and FAS antigen expression, or by repression of Bcl-2 expression. It has 393 amino acid residues.
SIRT1
SIRT 1 (Fig.3) is NAD-dependent protein deacetylase sirtuin-1 gene. It is a NAD-dependent protein deacetylase that links transcriptional regulation directly to intracellular energetics and participates in the coordination of several separated cellular functions such as cell cycle, response to DNA damage, metabolism, apoptosis and autophagy. It modulates chromatin function through deacetylation of histones and promote alterations in the methylation of histones and DNA, leading to transcriptional repression. This gene deacetylates a broad range of transcription factors and coregulators, thereby regulating target gene expression in cancer cells.
RPA1
RPA 1 (Fig.3) is replication protein A with 70 kDa DNA-binding subunit and an N-terminally processed amino acid. As part of the heterotrimeric replication protein A complex (RPA/RP-A), it binds and stabilizes single-stranded DNA intermediates, that form during DNA replication or upon DNA stress. It prevents their reannealing and in parallel, recruits and activates different proteins and complexes involved in DNA metabolism. Thereby plays an essential role both in DNA replication and the cellular response to DNA damage. In the cellular response to DNA damage, the RPA complex controls DNA repair and DNA damage checkpoint activation.
MDM2
MDM2 is (Fig.3) E3 ubiquitin-protein ligase Mdm2; E3 ubiquitin-protein ligase that mediates ubiquitination of p53/TP53, leading to its degradation by the proteasome. It inhibits p53/TP53- and p73/TP73-mediated cell cycle arrest and apoptosis by binding its transcriptional activation domain. MDM2 also acts as ubiquitin ligase E3 toward itself and ARRB1. It further permits the nuclear export of p53/TP53. The gene promotes proteasome-dependent ubiquitin-independent degradation of retinoblastoma RB1 protein and also inhibits DAXX-mediated apoptosis by inducing its ubiquitination and degradation.
CREBBP
CREBBP is (Fig.3) CREB-binding protein. This gene acetylates histones by giving a specific tag for transcriptional activation. It also acetylates non- histone proteins, like DDX21, FBL, IRF2, MAFG, NCOA3, POLR1E/PAF53 and FOXO1. It binds specifically to phosphorylated CREB and enhances its transcriptional activity toward cAMP-responsive genes. The gene acts as a coactivator of ALX1 and as a circadian transcriptional coactivator which enhances the activity of the circadian transcriptional activators: NPAS2-ARNTL/BMAL1 and CLOCK-ARNTL/BMAL1 heterodimers.
EP300
EP300 (Fig.3) is Histone acetyltransferase p300 gene. It functions as histone acetyltransferase and regulates transcription via chromatin remodeling. It also acetylates all four core histones in nucleosomes. Histone acetylation gives an epigenetic tag for transcriptional activation. EP300 mediates cAMP-gene regulation by binding specifically to phosphorylated CREB protein. It further mediates acetylation of histone H3 at 'Lys-122' (H3K122ac), a modification that localizes at the surface of the histone octamer and stimulates transcription, possibly by promoting nucleosome instability.
DAXX
DAXX (Fig. 3) is Death domain-associated protein 6 gene. The gene is a transcription corepressor known to repress transcriptional potential of several sumoylated transcription factors. The gene also down-regulates basal and activated transcription. Its transcription repressor activity is modulated by recruiting it to subnuclear compartments like the nucleolus or PML/POD/ND10 nuclear bodies through interactions with MCSR1 and PML, respectively. DAXX helps to regulate transcription in PML/POD/ND10 nuclear bodies together with PML and may influence TNFRSF6-dependent apoptosis and thereby inhibits transcriptional activation of PAX3 and ETS1.
ATM
ATM (Fig.3) is Serine-protein kinase. The Serine/threonine protein kinase activates checkpoint signaling upon double strand breaks (DSBs), apoptosis and genotoxic stresses such as ionizing ultraviolet A light (UVA), thereby acting as a DNA damage sensor. ATM gene recognizes the substrate consensus sequence [ST]-Q .It phosphorylates 'Ser-139' of histone variant H2AX at double strand breaks (DSBs), thereby regulating DNA damage response mechanism. It also plays a role in pre-B cell allelic exclusion, a process leading to expression of a single immunoglobulin heavy chain allele to enforce clonality and monospecific reactions.
HSP90AA1
HSP90AA1 (Fig. 3) is Heat shock protein HSP 90-alpha gene. It is a molecular chaperone that promotes the maturation, structural maintenance and proper regulation of specific target proteins involved for instance in cell cycle control and signal transduction. The gene undergoes a functional cycle that is linked to its ATPase activity which is essential for its chaperone activity. This cycle probably induces conformational changes in the client proteins, thereby causing their activation. The gene interacts dynamically with various co-chaperones that modulate its substrate recognition, ATPase cycle and chaperone function.
SFN
SFN (Fig.3) is 14-3-3 protein sigma gene. The adapter protein is implicated in the regulation of a large spectrum of both general and specialized signaling pathways. It binds to a large number of functional partners, usually by recognition of a phosphoserine or phosphothreonine motif. It’s binding generally results in the modulation of the activity of the binding partner. When bound to SIRT17, SFN regulates protein synthesis and epithelial cell growth by stimulating Akt/mTOR pathway. SFN also regulate MDM2 autoubiquitination and degradation and thereby activating p53/TP53 tumour protein.
TP53BP2
TP53BP2 (Fig.3) is Apoptosis-stimulating of p53 protein 2 gene. It is a regulator that plays a central role in regulation of apoptosis and cell growth via its interactions with proteins such as TP53. It regulates TP53 by enhancing the DNA binding and transactivation function of TP53 on the promoters of proapoptotic genes in vivo. It also inhibits the ability of APPBP1 to conjugate NEDD8 to CUL1, and thereby decreases APPBP1 ability to induce apoptosis. The gene impedes cell cycle progression at G2/M. Its apoptosis-stimulating activity is inhibited by its interaction with DDX42.
Genetic correlation analysis of TP53 and functional partner genes associated with hepatocellular carcinoma and liver cancer in normal and disease conditions.
Results of analysis of genetic correlation for TP53 and SIRT1 genes (Fig. 4a and Table 2) using the Genetic Correlation R Analyzer programme showed positive slope values (rate of mutation in genes), positive intercept values, positive determination coefficients R and weak, non-significant strength of associations (<0.5) in both cancerous and normal liver conditions. However, the degree of association was higher (37%) in liver cancer than in normal liver conditions (25%).
Results of analysis of genetic correlation for TP53 and RPA1 genes (figure 4b) using the Genetic Correlation R Analyzer programme showed very small negative slope values (rate of mutation in genes), positive intercept values, positive determination coefficients R and very strong, significant strength of associations (>0.5) in both cancerous and normal liver conditions. The degree of association was however higher (93%) in liver cancer than in normal liver conditions (87%).
Genetic correlation analysis for TP53 and MDM2 genes (Fig. 4c and Table 2) using the Genetic Correlation R Analyzer programme showed very small positive slope values (rate of mutation in genes), positive intercept values, positive determination coefficients R and very weak, non-significant strength of associations (>0.5) in both cancerous and normal liver conditions. The degree of association was very low (4%) in liver cancer and a bit higher in normal liver conditions (16%).
Genetic correlation analysis for TP53 and CREBBP genes (Fig. 4d and Table 2) using the Genetic Correlation R Analyzer programme showed very small positive slope values (rate of mutation of genes), positive intercept values, positive determination coefficients R and very weak, non-significant strength of association (>0.5) in cancerous condition but strong, significant strength of association (>0.5) in normal liver conditions. The degree of association was very low (5.9%) in liver cancer but very high in normal liver conditions (82%).
Results of analysis of genetic correlation for TP53 and EP300 genes (Fig. 4e and Table 2) using the Genetic Correlation R Analyzer programme showed very small positive slope values (rate of mutation of genes), positive intercept values, positive determination coefficients R and a weak, non-significant strength of association (>0.5) in cancerous liver and a strong, significant strength of association (>0.5) in normal liver conditions. The degree of association was however low (24%) in liver cancer but high in normal liver conditions (75%).
Results of analysis of genetic correlation for TP53 and DAXX genes (Fig. 4f and Table 2) using the Genetic Correlation R Analyzer programme showed very small negative slope values (rate of mutation of genes), positive intercept values, positive determination coefficients R and strong, significant strength of associations (>0.5) in both cancerous and normal liver conditions. The degree of association was however higher (89%) in liver cancer than in normal liver conditions (88%).
Results of analysis of genetic correlation for TP53 and ATM genes (Fig. 4g and Table 2) using the Genetic Correlation R Analyzer programme showed very small positive slope values (rate of mutation of genes), positive intercept values, positive determination coefficients R and very weak, non-significant strength of associations (>0.5) in both cancerous and normal liver conditions. The degree of association was (10%) in liver cancer and 12% in normal liver conditions.
Results of analysis of genetic correlation for TP53 and HSP90AA1 genes (Fig. 4h and Table 2) using the Genetic Correlation R Analyzer programme showed very small positive slope value (rate of mutation in genes) in liver cancer and a very small negative slope value (rate of mutation of genes) in normal liver conditions, positive intercept values, positive determination coefficients R and strong, significant strength of associations (>0.5) in both cancerous and normal liver conditions. The degree of association was however higher (78%) in liver cancer than in normal liver conditions (50%).
Results of analysis of genetic correlation for TP53 and SFN genes (Fig. 4i and Table 2) using the Genetic Correlation R Analyzer programme showed very small negative slope value (rate of mutation in genes) in liver cancer and very small positive slope value (rate of mutation of genes) in normal liver condition, positive intercept values, positive determination coefficients R and a strong, significant strength of associations (>0.5) in both cancerous and normal liver cases. The degree of association was however higher (77%) in normal liver than in liver cancer conditions (73%).
Results of analysis of genetic correlation for TP53 and TP53BP2 genes (Fig. 4j and Table) using the Genetic Correlation R Analyzer programme showed very small positive slope values (rate of mutation of genes), positive intercept values, positive determination coefficients R and very weak, non-significant strength of association (<0.5) in liver cancer and a strong, significant strength of association (>0.5) in normal liver condition. The degree of association was higher (83%) in normal liver than in liver cancer conditions (5.2%).
Table 2: Genetic correlation and P- values of TP53 and related functional partner genes associated with Hepatic carcinoma and Liver cancer among patients under consideration.
S/N
|
Gene name
|
Description
|
TP53 Correlation in Liver Cancer
|
TP53 Correlation in Normal Liver
|
Rate of change (mutation of genes)
|
P-Value
|
1
|
TP53
|
Tumour protein 53
|
|
|
|
|
2
|
SIRT 1
|
Sirtuin 1
|
Y=0.028 +0.34x; R= 0.37
|
Y=0.018 +0.14x; R= 0.25
|
Small and positive for both conditions
|
0.00124
|
3
|
ATM
|
ATM serine/Threonine kinase
|
Y=0.043 +0.13x; R= 0.1
|
Y=0.083 +0.11x; R= 0.12
|
Small and positive for both conditions
|
0.001687
|
4
|
TP53BP2
|
Tumour protein P53 binding protein 2
|
Y=0.045 +0.047x; R= 0.052
|
Y=0.016 + 0.66x; R= 0.83
|
Small and positive for both conditions
|
0.001792
|
5
|
MDM2
|
MDM2 proto-oncogene
|
Y=0.032 +0.05x; R= 0.44
|
Y=0.04 +0.11x; R= 0.16
|
Small and positive for both conditions
|
0.002179
|
6
|
RPA 1
|
Replication protein A1
|
Y=-0.0059 +1.3x; R= 0.93
|
Y=-0.022 +0.95x; R= 0.87
|
Small and negative for both conditions
|
0.000001465
|
7
|
DAXX
|
Death domain associated protein
|
Y=-0.0097 +1.2x; R= 0.89
|
Y=-0.0042 +0.91x; R= 0.88
|
Small and negative for both conditions
|
0.000002221
|
8
|
HSP90AA1
|
Heat shock protein 90 alpha family class A member 1
|
Y=0.014 +0.77x; R= 0.78
|
Y=-0.08 +0.52x; R= 0.5
|
Small and positive in cancer but negative in normal condition
|
0.00001064
|
9
|
SFN
|
Stratifin
|
Y=-0.018 +0.82x; R= 0.72
|
Y=0.0066 +0.63x; R= 0.77
|
Small and negative for cancer but positive for normal condition
|
0.0000265
|
10
|
CREBBP
|
CREB binding protein
|
Y=0.041 +0.049x; R= 0.059
|
Y=0.051 +0.069x; R= 0.82
|
Small and positive for both conditions
|
0.00002827
|
11
|
EP300
|
E1A binding protein P300
|
Y=0.033 +0.22x; R= 0.24
|
Y=0.061 +0.68x; R= 0.75
|
Small and positive for both conditions
|
0.00005194
|
The results of determination coefficients R for TP53 and functional partner genes is presented in figure 5. The coefficient depicts the strength of association between TP53 and each partner gene in normal liver and liver cancer conditions. The degree or strength of association was higher in normal liver for 6 of the 10 evaluations while 4 showed high strength of associations in liver cancer conditions.
The result of functional partner gene reaction hotspots of TP53 and other genes are presented in Table 3. The prosites in silico tool of expasy.org suite was used to obtain the sequence for each functional partner gene from the protein – protein interaction (Figure 2) network and used for the scanning, analysis, and identification of the reaction prosites of these genes in liver cancer cells. The results revealed that the protein kinase II phosphorylation site II was the most functional protein reactional site for the TP53 and associated genes in liver cancer initiation and progression. Other functional protein domain hotspot sites for TP53 and associated genes included the N-glycosylation protein site, N -Myristylation protein site, N-Amidation protein sites, Tyrosine kinase 1 and II phosphorylation sites, cAMP and cGMP phosphorylation sites as well as the casein kinase II phosphorylation sites (Table 3).
Table 3: Functional Domain Reaction Hotspots for TP53 and functional gene partners in hepatocellular carcinoma and Liver Cancer conditions
Gene
Name
|
Description
|
Molecular protein Family signature
|
PKC
|
cAMP/cGMP site
|
N-Glycosylation site
|
N-Myristylation site
|
Tyrosine Phosphorylation site 1
|
Tyrosine phosphorylation site 2
|
CK2
|
N-Amidation site
|
Rich region
|
TP53
|
Tumour protein 53
|
P53 family signature (237 – 249 bp)
|
99 – 101 bp
|
-
|
239 – 242 bp
|
117 – 122 bp
|
-
|
320 – 327 bp
|
18 -21 bp
|
-
|
Proline
(58 -92 bp
|
SIRT 1
|
Sirtuin 1
|
Sirtuin Catalytic domain signature 236 – 496 bp
|
177- 179
|
-
|
342 – 345 bp
|
12 – 17 bp
|
635 – 642 bp
|
516 – 523 bp
|
47 – 50 bp
|
-
|
Alanine (54 -–98 bp)
|
ATM
|
ATM serine/Threonine kinase
|
FAT domain family signature (1940 – 2566 KDa)
|
21 – 23 TeR
|
1920 - 1923 KRpS
|
81 - 84 NVSA
|
134 - 139 GAIYGA
|
116 – 123 bp
|
2160 – 2167 KRSLFSVY
|
127 - 130 TVKD
|
|
Leucine (1217 – 1238 KDa)
|
TP53BP2
|
Tumour protein P53 binding protein 2
|
P53 protein family signature
|
7 -9 TaR
|
9-12 RKsT
116 – 119 KRvT
|
|
|
|
|
|
|
Proline
(76 – 79)
VaTR
|
MDM2
|
MDM2 proto-oncogene
|
Zinc finger domain family signature 309 – 328 KDa
|
55 – 57 TmK
|
-
|
11-14 NMSV
|
18 – 23 GAVtTS
|
|
|
28 – 31 SeqE
|
|
-
|
RPA 1
|
Replication protein A1
|
Replication factor protein family signature
|
43 -45 SrK
|
AsTR
129 -132
|
|
|
|
|
|
|
-
|
DAXX
|
Death domain associated protein
|
Protein signature
|
138 – 140
SaK
|
590 -593
RKqS
|
364 – 367 NRSL
|
8-13 GAiaAI
|
595 -602 KPvDyreY
|
327 – 334 RrhlDliY
|
40 -43 SssE
|
54 – 57 gGKK
|
Glutamic acid 434 – 572 KDa
|
HSP90AA1
|
Heat shock protein 90 alpha family class A member 1
|
Heat shock protein family signature 160 – 169 KDa
|
16 -18 SIR
|
330 – 333 KHhS
|
173 -176 NSSD
|
6 -11 GGdqsT
|
275 – 282 KhnDdegT
|
635 – 642 KhglEviY
|
5 – 8 SqqD
|
|
Glutamic acid 345 – 390 KDa
|
SFN
|
Stratifin
|
Protein signature 1 41 – 51 KDa
|
146 – 148 SaR
|
|
175 -178 NFSV
|
53 – 58 GGqrAA
|
77 – 84 kqPEVreY
|
|
37 – 40 SceE
|
|
Proline 213 – 232 KDa
|
CREBBP
|
CREB binding protein
|
ABC Transporter family signature 123 – 137 bp
|
387 – 389 bp
|
668 – 671 bp
|
25 – 28 bp
|
48 -53 bp
|
634 – 641 bp
|
|
23 – 26 bp
|
1200 – 1203 bp
|
Glutamine (683 -960 bp)
|
EP300
|
E1A binding protein P300
|
Bromodomain family signature 1067 – 1139 KDa
|
12 – 14 SaK
|
647 – 650 KRrT
|
49 – 52 NSTE
|
30 -35 GTdfGS
|
613 – 620 RKVEgdmY
|
|
26 – 29 SasD
|
1164 – 1167 cGRK
|
Proline 666 – 943 KDa
|
Single Linkage genetic cluster and principal component analysis of TP53 and functional partner genes in hepatocellular carcinoma and liver cancer patient results are presented in Table 4. Two major principal components were used to explain the dimensionality of variations in gene action in liver cancer conditions. The Serine/Threonine kinase ATM gene showed the highest impact (0.46) on the first principal component while the heat shock protein 90 alpha family class (HSP90AA1) gene impacted the least (-0.0822) effect on PC1. The contributions of the different genes to PC2 showed all negative dimensionality with SFN gene contributing the highest loading value of -0.0958 to the dimensionality observed in PC2 while the least loading value of -0.434 to the overall variation was contributed by RPA1 gene. The TP53 and associated gens were delineated into two major cluster groups based on their dimensionality. Cluster one had TP53, RPA1, DAXX and SFN genes. Cluster 2 has ATM, MDM2, EP300, SIRT 1, CREBBP and HSP90AA1 genes as sub cluster members (Table 4). The red-coloured site (Fig. 6) revealed axis with high reaction intensity among the TP53 and associated genes. The green-coloured axis portray interactions that are very weak and less intense while the black-coloured axis reveals areas of moderate to average interactions among TP53 and associated functional partner genes (Fig. 6)
Table 4: Single linkage genetic cluster and principal component analysis of TP53 and other functional partner genes in hepatocellular carcinoma and liver cancer patients