Machine Learning Discoveries of HVCN1-X Synergy in ETC-1922159 Treated Colorectal Cancer Cells

Machine learning discoveries of HVCN1-X synergy in ETC-1922159 treated colorectal cancer cells shriprakash sinha Independent Researcher; Orcid ID : orcid.org/0000-0001-7027-5788 104-Madhurisha Heights Phase 1, Risali, Bhilai-490006, India Abstract HVCN1 is the gene that encodes for voltage gated proton channel. It works along with other genes/proteins, however, many of the synergistic combinations are not known. In colorectal cancer (CRC) cells treated with ETC-1922159, HVCN1 was found to be dow regulated along with other genes. A recently developed search engine ranked combinations of hydrogen voltage gated channel HVCN1-X (X, a particular gene/protein) at 2nd order level after drug administration. These rankings reveal which HVCN1-X combinations might be working synergistically in CRC. In this research work, I cover combinations of HVCN1 with NADPH oxidase 1 (NOX1), ecto-NOX disulfide-thiol exchanger 2 (ENOX2), mitogen-activated protein kinase (MAPK), matrix metallopeptidase (MMP), tumor protein p53 (TP53), marker of proliferation Ki-67 (MKI67), glutathione s-transferases (GST) and phospholipase C (PLC) family. Keywords: HVCN1, Porcupine inhibitor ETC-1922159, Sensitivity analysis, Colorectal cancer. 1. Introduction 1.1. Ion channels Ion channels, as the name suggests, as we deconstruct the terminology, in a lay man’s term means, a channel or a passage that helps pass the ions. When we talk of passage, then the flow of the ions has to happen through something permeable, a wall, a barrier,5 an obstruction, etc. Also, the channel should open and close with some purpose and there has to be some definite aim, for which the channel behind which the passage will work. For other conditions, the channel has to remain closed. Investigating these conditions, will help in understanding the ambinence of the channel itself. Another aspect of particular importance is the fact regarding the kind of ions that will pass through10 IML dicoveries of HVCN1-X synergy in ETC-1922159 treated CRC cells Email address: [email protected] (shriprakash sinha) 1Aspects of unpublished work were presented in a poster session at (1) the recently concluded first ever Wnt Gordon Conference, from 6-11 August 2017, held in Stowe, VT 05672, USA. Preprint submitted to Preprint December 23, 2024 the passage and why such ions need to pass through such a channel. Along with the study of these channels, comes the study of the varied structure of these channels which help the passage of these ions. All these issues gives rise to the field of simulation and modeling of these channels also. Broadly speaking, ion channel is one of the two classes of an ionophoric pro-15 teins, that will let the ions pass downhill along with the electrochemical gradient in the channel without the input of any metabolic activity; the other being ion transporters. Ionophoric proteins are structures that reversibly binds with ions. Ion transporters are proteins that will carry the ions along with them through a permeable wall or barrier or an obstruction against a concentration gradient. This passage of the ions requires20 openning and closing mechanism in the channel and in electrophysiology, gating is the term used for such a mechanism. It is the process were the channel transforms from a conducting mode where ions flow to a non-conducting mode were ions do not pass through. The rate of closing and openning is referred to as kinetics of gating. There are two major kinds of ion channels namely, (a) Voltage gated ion channel25 and (b) Ligand gated ion channel. The former operates via the difference in barrier or membrane voltage potential while the latter works via the binding of ligand to the channel. Also, there are different kinds of ions that go through the membrane. These include chloride, potassium, sodium, calcium and proton, to name a few. Depending on the kind of ions, certain functionalities are characteristics of a respective ion channel30 under study. But before I delve deeper into the specific area of interest, it is important to ask, what is the purpose of such a channel and why this gating process in helpful? The answer to this question will also address the ambience in which the channel is working or not working and thus build the context for the above posed question. 1.2. Voltage gated proton ion channel - HVCN135 The voltage gated proton ion channel is a unique channel unlike other ion channels, as it has some special properties. The gated channel is not a pump as it works without the help of metabolic activity of ATP. It is a higly selective channel and has a perfect selectivity in comparison to the other ion channels which are sometime imperfect in nature. With strong dependence on the temperature and the lack of an aqueous pore that40 probably makes the selectivity extremely perfect, the channel distinguishes in making the flow of the hydrogen proton outward. Also, since the flow of protons can happen in various ways, the route of which the other ions cannot take, voltage gated ion channels form a separate class in themselves. That, HVCN1 is the gene that encodes for voltage gated proton channel was dis-45 covered recently in 2006 by two different group simulataneously (i.e Ramsey et al. [1] and Sasaki et al. [2]). The channels resembles the voltage sensing domain VSD of ordinary voltage gated cation channels. Structurally, the proton channel constitutes four membrane-spanning region S1-S4 and lacks the separate pore domain. The human Hv1 forms a dimer in the membrane. Lee et al. [3] shows this dimeric structure in figure 150 as a probable model. Advanced review on the topic can be found in Rebolledo et al. [4]. Earlier we talked about the regulation of the pH. This regulation happens by balancing the acidity, i.e the concentration of the hydrogen protons in the cells. What we 2 Figure 1: Dimeric structures at extracellular level mediated by S1 and a colied coil at the intracellular domain. Reprinted with permission from Lee et al. [3]. Copyright 2008 National Academy of Sciences, USA. find is that, if there is excess of concentration of the hydrogen ions in the cell, then55 the proton channel works by gating the protons in an outward direction to stabilize the concentration to a physiological optimum for proper functioning of the cell. However, when the protons are being directed outside, the balance of the concentration in the thriving environment of the cell is disturbed. To compensate for this, simultaneously, a flux of the electrons is made possible via NADPH oxidase. Thus the balance60 in concentration is maintained. As a consequence of the above phenomena, one of the jobs of proton channels is to regulate the acidity by extrusion of the acidic material at intracellular level and thus increase the pHi. Simultaneously, from the extracellular perspective the acid secretion leads to the regulation of acidity outside the cell and hence the decrease of pHo. Finally, the extrusion of protons has a controlled effect on65 the osmotic swelling in comparison to other protons. 1.3. HVCN1 in colorectal cancer Wang et al. [5] have recently shown that HVCN1 was found to be highly expressed in colorectal cancer cell lines. Suppression of Hv1 in highly metastatic colorectal cancer restrained the extrusion of protons and thus leading to low pHi. Since the expression70 of Hv1 was found to be very high in colorectal cancer cells, the extrusion of protons is expected to be high. What this means is that there is an increase in the pHiand cancer cells are able to thrive in hypoxic environments where the acidity of the envirnoment is beyond the normal level. A more recent role of voltage-gated proton channel (Hv1) in cancer biology has been studied by Alvear-Arias et al. [6]. El Chemaly et al. [7]75 in their study, show that the genetic deletion of Hv1 delays tumor development in a mouse model of granulocytic sarcoma and report the discovery and characterization of two novel bioavailable inhibitors of Hv1 channels that they validate by orthogonal assays and electrophysiological recordings. HVCN1 was found to be highly suppressed in colorectal cancer cell lines after the treatment of ETC-1922159 drug as observed80 in Madan et al. [8]. Most studies till now, have dealt with how the HVCN1 works and there is very less informnation reagarding which gene/protein combinations might 3 be working synergistically along with it. It would be nice to observe if there is any connection between the independently observed factors in the form of unknown biological hypotheses. To solve the issue, the next section addresses poses a solution to85 the problem. 1.4. Combinatorial search problem and a possible solution In a recently published work Sinha [9], a frame work of a search engine was developed which can rank combinations of factors (genes/proteins) in a signaling pathway. Readers are requested to go through the adaptation of the above mentioned work for90 gaining deeper insight into the working of the pipeline and its use of published data set generated after administration of ETC-1922159, Sinha [10]. The work uses SVM package by Joachims [11] in https://www.cs.cornell.edu/people/tj/svm_light/ svm_rank.html. I use the adaptation to rank 2nd order gene combinations. 2. Results & Discussion95 2.1. HVCN1 related synergies 2.1.1. HVCN1-NOX1/ENOX2 The reactive oxygen species (ROS)-producing enzyme NADPH oxidase (NOX) was first identified in the membrane of phagocytic cells. Its discovery and molecular mechanisms has been described in detail in Vermot et al. [12]. Hv1 channels sustain high-100 level production of ROS by NOX2 in phagocytes (Ramsey et al. [13], El Chemaly et al. [14] and El Chemaly et al. [7]). In colorectal cancer cells treated with ETC-1922159, NOX1 and HVCN1, were found to be down regulated and recorded independently. Ecto-NOX disulfide-thiol exchanger 2 (ENOX2), is a member of the NOX family of NADPH oxidases. ENOX2 proteins are restricted to cancer (Hostetler et al. [15],105 Chueh et al. [16]). In colorectal cancer cells treated with ETC-1922159, ENOX2 and HVCN1, were found to be down regulated and recorded independently. I was able to rank 2nd order combination of NOX1/ENOX2 and HVCN1, that were down regulated. Table 1 shows rankings of these combinations. Followed by this is the unexplored combinatorial hypotheses in table 2 generated from analysis of the ranks in table 1.110 The table 1 shows rankings of NOX family w.r.t HVCN1. NOX1 - HVCN1 shows low ranking of 349 (laplace) and 855 (linear). ENOX2 - HVCN1 shows low ranking of 1404 (laplace), 601 (linear) and 1503 (rbf). These rankings point to the synergy existing between the two components, which have been down regulated after the drug treatment.115 One can also interpret the results of the table 1 graphically, with the following influences - •NOX family w.r.t HVCN1 with HVCN1 −>NOX1; and HVCN1 −> ENOX2. 4 RANKING NOX FAMILY VS HVCN1 RANKING OF NOX FAMILY W.R.THVCN1 laplace linear rbf NOX1 - HVCN1 349 855 2050 ENOX2 - HVCN1 1404 601 1503 Table 1: 2nd order interaction ranking between HVCN1 VS NOX family members. UNEXPLORED COMBINATORIAL HYPOTHESES NOX family w.r.t HVCN1 NOX1 HVCN1 ENOX2 HVCN1 Table 2: 2nd order combinatorial hypotheses between HVCN1 and NOX family members. 5 2.1.2. HVCN1-MAPK/MMP Banskota et al. [17] show that NOX2-derived ROS regulates NOX1, NOX2, and matrix120 metalloproteinase-7 (MMP7) expression through the mitogen-activated protein kinase (MAPK) pathway, in colon cancer. Based on this experimentally tested phenomena, there might be possible connection between HVCN1 and MMP/MAPK, though the current understanding of how the synergy works might be lacking. In ETC-1922159 treated colorectal cancer cells, MMP11 and MAPK15 were both down regulated and125 their regulation was recorded individually. I was able to rank 2nd order combination of MMP11/MAPK15 and HVCN1, that were down regulated. Table 3 shows rankings of these combinations. Followed by this is the unexplored combinatorial hypotheses in table 4 generated from analysis of the ranks in table 3. The table 3 shows rankings of MMP11/MAPK15 w.r.t HVCN1. MAPK15 - HVCN1130 shows low ranking of 1612 (linear) and 60 (rbf). The ranking for MAPK15 points to the possible synergy existing between the two components, which have been down regulated after the drug treatment. However, the ranking of MMP11 does not point to possible synergy with HVCN1 (majority getting high numerical valued rank for down regulation, see table).135 RANKING MMP11/MAPK15 VS HVCN1 RANKING OF MMP11/MAPK15 W.R.THVCN1 laplace linear rbf MMP11 - HVCN1 2297 618 2123 MAPK15 - HVCN1 2314 1612 60 Table 3: 2nd order interaction ranking between HVCN1 VS MMP/MAPK members. One can also interpret the results of the table 3 graphically, with the following influences - •MAPK15 w.r.t HVCN1 with HVCN1 −>MAPK15. UNEXPLORED COMBINATORIAL HYPOTHESES MMP/MAPK w.r.t HVCN1 MAPK15 HVCN1 Table 4: 2nd order combinatorial hypotheses between HVCN1 and MMP/MAPK members. 6 2.1.3. HVCN1clinicopathological parameters TP53/MKI67/GST In colorectal cancer, while studying Hv1 as a new potential biomarker for diagnosis and prognosis, Wang et al. [5] reported clinicopathological parameters like TP53, MKI67140 and GST. Correlation between Hv1 expression levels in colorectal cancer and clinicopathological parameters, were recorded. In ETC-1922159 treated colorectal cancer cells, TP53, MKI67 and members of GST family were found down regulated and their regulation was recorded individually. I was able to rank 2nd order combination of TP53/MKI67/GST and HVCN1, that were down regulated.145 Table 5 shows rankings of these combinations. Followed by this is the unexplored combinatorial hypotheses in table 6 generated from analysis of the ranks in table 5. The table 5 shows rankings of TP53/MKI67/GST w.r.t HVCN1. TP53 - HVCN1 shows low ranking of 1532 (laplace) and 1640 (rbf). MKI67 - HVCN1 shows low ranking of 184 (linear) and 572 (rbf). MKI67 - HVCN1 shows low ranking of 184 (linear) and 572150 (rbf). GSTO2 - HVCN1 shows low ranking of 689 (laplace), 1341 (linear) and 1591 (rbf). GSTM4 - HVCN1 shows low ranking of 741 (laplace) and 593 (rbf). RANKING CLINICOPATHOLOGICAL PARAMETERS TP53/MKI67/GST VS HVCN1 RANKING OF TP53/MKI67/GST FAMILY W.R.THVCN1 laplace linear rbf TP53 - HVCN1 1532 2541 1640 TP53BP1 - HVCN1 1259 2552 1972 MKI67 - HVCN1 2713 184 572 GSTM3 - HVCN1 23 2431 2531 GSTA4 - HVCN1 201 1767 2047 MGST2 - HVCN1 527 2163 2409 MGST1 - HVCN1 627 1980 2361 GSTO2 - HVCN1 689 1341 1591 GSTM4 - HVCN1 741 1716 593 Table 5: 2nd order interaction ranking between HVCN1 VS TP53/MKI67/GST members. One can also interpret the results of the table 5 graphically, with the following influences - •TP53 w.r.t HVCN1 with HVCN1 −>TP53; •MKI67 w.r.t HVCN1 with HVCN1 −>MKI67; and •GST w.r.t HVCN1 with HVCN1 −>GSTO2; HVCN1155 −>GSTM4. 2.1.4. HVCN1-PLC Alvear-Arias et al. [6], state that the effect called ”enhanced gating mode” is a phenomenon that occurs after stimulation with phorbol 12-meristate-13-acetate (PMA). PMA addition to cells, promotes a shift of activation curves of Hv1 channel currents160 to negative voltages and faster opening kinetics (B´ anfi et al. [18], DeCoursey et al. [19]). PMA operates by activating the phospholipase C (PLC) that leads to the activation of protein kinase C (PKC). PKC phosphorylates Threonine 29 in the N-terminal of hHv1, leading to the enhancement of its gating Musset et al. [20]. These experi7 UNEXPLORED COMBINATORIAL HYPOTHESES Clinicopathological parameters w.r.t HVCN1 TP53 HVCN1 MKI67 HVCN1 GST-O2/M4 HVCN1 Table 6: 2nd order combinatorial hypotheses between HVCN1 and Clinicopathological parameters. mental findings point to the synergy that might be existing between PLC and HVCN1.165 In ETC-1922159 treated colorectal cancer cells, PLC family and HVCN1 were found down regulated and their regulation was recorded individually. I was able to rank 2nd order combination of PLC family and HVCN1, that were down regulated. Table 7 shows rankings of these combinations. Followed by this is the unexplored combinatorial hypotheses in table 8 generated from analysis of the ranks in table 7.170 The table 7 shows rankings of PLC family w.r.t HVCN1. PLCH1 - HVCN1 shows low ranking of 494 (laplace), 929 (linear) and 922 (rbf). PLCB4 - HVCN1 shows low ranking of 1279 (laplace), 1248 (linear) and 677 (rbf). PLCB2 - HVCN1 shows low ranking of 988 (linear) and 1501 (rbf). RANKING PLC FAMILY VS HVCN1 RANKING OF PLC FAMILY W.R.THVCN1 laplace linear rbf PLCH1 - HVCN1 494 929 922 PLCG1 - HVCN1 1275 2682 1958 PLCB4 - HVCN1 1279 1248 677 PLCB2 - HVCN1 2324 988 1501 Table 7: 2nd order interaction ranking between HVCN1 VS PLC members. One can also interpret the results of the table 7 graphically, with the following175 influences - •PLC family w.r.t HVCN1 with HVCN1 −>PLCH1; HVCN1 −> PLCB4; and HVCN1 −>PLCB2. . 8 UNEXPLORED COMBINATORIAL HYPOTHESES PLC family w.r.t HVCN1 PLC-H1/B4/B2 HVCN1 Table 8: 2nd order combinatorial hypotheses between HVCN1 and PLC family. 3. Conclusion Presented here are a range of multiple synergistic HVCN1 2nd order combinations that180 were ranked via a machine learning based search engine. Via majority voting across the ranking methods, it was possible to find plausible unexplored synergistic combinations of HVCN1-X that might be prevalent in CRC cells after treatment with ETC-1922159 drug. Conflict of interest185 There are no conflicts to declare. Author’s contributions Concept, design, in silico implementation - SS. Analysis and interpretation of results - SS. Manuscript writing - SS. Manuscript revision - SS. Approval of manuscript - SS Acknowledgements190 Special thanks to Mrs. Rita Sinha and Mr. Prabhat Sinha for supporting the author financially, without which this work could not have been made possible. Source of Data Data used in this research work was released in a publication in Madan et al. [8]. 4. References195 References [1] I. S. Ramsey, M. M. Moran, J. A. Chong, D. E. Clapham, A voltage-gated proton-selective channel lacking the pore domain, Nature 440 (2006) 1213–1216. 9