Choline kinase and miR-32-5p: A crucial interaction promoting apoptosis and delaying wound repair in cervical cancer cells

Sweta Raikundalia; Ling Ling Few; Siti Asma' Hassan; Get Bee Yvonne-Τee; Wei Cun See Too; Sweta Raikundalia; Ling Ling Few; Siti Asma' Hassan; Get Bee Yvonne-Τee; Wei Cun See Too

doi:10.3934/biophy.2024016

AIMS Biophysics

2024, Volume 11, Issue 3: 281-295. doi: 10.3934/biophy.2024016

Previous Article Next Article

Research article

Choline kinase and miR-32-5p: A crucial interaction promoting apoptosis and delaying wound repair in cervical cancer cells

1.
School of Health Sciences, Health Campus, Universiti Sains Malaysia, 16150 Kubang Kerian, Kelantan, Malaysia
2.
School of Medical Sciences, Health Campus, Universiti Sains Malaysia, Kubang Kerian, Kelantan 16150, Malaysia

Received: 12 May 2024 Revised: 08 August 2024 Accepted: 19 August 2024 Published: 21 August 2024

Cervical cancer remains a significant global health challenge, especially in low- and middle-income countries. The elucidation of the molecular pathways underlying the carcinogenesis of cervical cancer is vital to develop targeted therapies. Our study investigates the regulatory effects of miR-32-5p on the choline kinase alpha (chka) gene in HeLa cells, a well-established cervical cancer cell line. Choline kinase alpha is recognized for its role in phosphatidylcholine biosynthesis, which is crucial for cell membrane formation, and is implicated in the oncogenic transformation of cells. Utilizing a combination of in silico prediction, luciferase assays, RT-qPCR, and Western blot analyses, we demonstrated that miR-32-5p directly targets the 3′ untranslated region (3′UTR) of chka mRNA, leading to a significant downregulation of chka expression. Our results demonstrate that miR-32-5p significantly downregulates chka at both the mRNA and protein levels, thus leading to decreased cellular proliferation and increased apoptosis. This was further confirmed by a cell cycle analysis, which showed a notable arrest in the G0/G1 phase. Additionally, scratch assays indicated a reduced migratory capacity in miR-32-5p-transfected cells, suggesting the potential anti-metastatic properties of miR-32-5p. These findings highlight the therapeutic potential of miR-32-5p as a biomarker and a target in cervical cancer treatment strategies. By modulating chka expression, miR-32-5p could serve as a novel approach to curb the progression and spread of cervical cancer, thus offering a promising avenue for future research and clinical applications. This study contributes to the growing understanding of miRNA-mediated gene regulation in cancer biology and underscores the importance of targeted genetic research in the development of cancer therapeutics.

Keywords:

Citation: Sweta Raikundalia, Ling Ling Few, Siti Asma' Hassan, Get Bee Yvonne-Τee, Wei Cun See Too. Choline kinase and miR-32-5p: A crucial interaction promoting apoptosis and delaying wound repair in cervical cancer cells[J]. AIMS Biophysics, 2024, 11(3): 281-295. doi: 10.3934/biophy.2024016

Related Papers:

[1]	Noman Zahid, Ali Hassan Sodhro, Usman Rauf Kamboh, Ahmed Alkhayyat, Lei Wang . AI-driven adaptive reliable and sustainable approach for internet of things enabled healthcare system. Mathematical Biosciences and Engineering, 2022, 19(4): 3953-3971. doi: 10.3934/mbe.2022182
[2]	Jianfei Cai, Guozheng Yang, Jingju Liu, Yi Xie . FastCAT: A framework for fast routing table calculation incorporating multiple protocols. Mathematical Biosciences and Engineering, 2023, 20(9): 16528-16550. doi: 10.3934/mbe.2023737
[3]	B. Kiruthika, Shyamala Bharathi P . Intelligent dynamic trust secure attacker detection routing for WSN-IoT networks. Mathematical Biosciences and Engineering, 2023, 20(2): 4243-4257. doi: 10.3934/mbe.2023198
[4]	Aida Ben Chehida Douss, Ryma Abassi, Damien Sauveron . State-of-the-art survey of in-vehicle protocols and automotive Ethernet security and vulnerabilities. Mathematical Biosciences and Engineering, 2023, 20(9): 17057-17095. doi: 10.3934/mbe.2023761
[5]	Xiaojun Liang . RETRACTED ARTICLE: A video images-aware knowledge extraction method for intelligent healthcare management of basketball players. Mathematical Biosciences and Engineering, 2023, 20(2): 1919-1937. doi: 10.3934/mbe.2023088
[6]	Suresh Y, Kalaivani T, Senthilkumar J, Mohanraj V . BF² VHDR based dynamic routing with hydrodynamics for QoS development in WSN. Mathematical Biosciences and Engineering, 2020, 17(1): 930-947. doi: 10.3934/mbe.2020050
[7]	Xiangyang Ren, Shuai Chen, Liyuan Ren . Optimization of regional emergency supplies distribution vehicle route with dynamic real-time demand. Mathematical Biosciences and Engineering, 2023, 20(4): 7487-7518. doi: 10.3934/mbe.2023324
[8]	Lingling Li, Congbo Li, Li Li, Ying Tang, Qingshan Yang . An integrated approach for remanufacturing job shop scheduling with routing alternatives. Mathematical Biosciences and Engineering, 2019, 16(4): 2063-2085. doi: 10.3934/mbe.2019101
[9]	Tianjun Lu, Xian Zhong, Luo Zhong, RuiqiLuo . A location-aware feature extraction algorithm for image recognition in mobile edge computing. Mathematical Biosciences and Engineering, 2019, 16(6): 6672-6682. doi: 10.3934/mbe.2019332
[10]	Pasqual Martí, Jaume Jordán, Pablo Chamoso, Vicente Julian . Taxi services and the carsharing alternative: a case study of valencia city. Mathematical Biosciences and Engineering, 2022, 19(7): 6680-6698. doi: 10.3934/mbe.2022314

Abstract

1. Introduction

Biological systems are permeated with deterministic laws and randomness. Noise are inevitable at every level of biology, from the most essential molecular, sub-cellular processes to the kinetics of tissues, organs, organisms and populations ^[1]. According to different origination of noise, noise can be divided into two terms: intrinsic and extrinsic noise ^[2,3,4,5]. Intrinsic noise originates from the stochasticity of biochemical events, while extrinsic noise is caused by status differences in individual cells, e.g., ribosomes and metabolites, the concentration of RNA polymerases or local environmental conditions ^[6]. Some experimental and theoretical studies have shown that noise may lead to genetic diseases, gene damage, etc. ^[7,8]. However, over the past few decades, it has been revealed that noise may also play a positive role in various dynamic systems and dynamic processes ^[9,10], for example, stochastic transitions ^[11] and switch processes ^[3,12].

Experimental studies have exposed that fluctuations occur in various stages of gene expression, including transcription, degradation, translation, binding, and so on ^[6,13,14] and noise affects gene expression ^[15]. Smolen et al. ^[16] proposed a simple model of gene transcription system dynamics, which is a bistable model with a positive feedback loop. They pointed out that the simple model manifested multiple stable state and brief perturbations could switch the model between those states. Based on this model, Liu et al. showed that the switch process could be successfully induced by the fluctuation of the rate of degradation and the synthesis of transcription factors ^[17]. The state transition of a rapid response to environmental changes is a common and important phenomenon in biological systems. Noise can also be used to regulate the dynamics of bistable genetic regulatory systems ^[18]. Wang et al. ^[19] further considered the effect of a time delay on the switch process. The related studies have focused on white Gaussian noise and colored noise ^{[3,17,19,20,21,22,23,24]}. However, the range of the Gaussian noise is unbounded and large values can be obtained, but the nature of real physical quantities is always bounded. The use of Gaussian noise may lead to unreasonable results in some cases ^[25,26]. Therefore, researchers have explored the modelling of bounded noise in dynamic systems ^[27]. Among them, sine-Wiener noise is the most widely used form of bounded noise ^[28] that is generated by a sinusoidal function with a constant amplitude, constant frequency and random change of phase with the Wiener process. Sinusoidal functions are bounded, resulting in their noise values not exceeding a fixed amplitude. In recent years, the influence of cross-correlated sine-Wiener (CCSW) noises on nonlinear dynamic systems has attracted wide attention ^[29,30]. Furthermore, our previous studies ^{[31,32,33,34]} revealed temporal coherent resonance phenomenon induced by CCSW noises in regular and small-world neuronal networks and indicated the occurrence of CCSW noise-induced coherent resonance. The increasing degree of correlation of CCSW noises may enhance or impair temporal regularity, which is dependent on the cross-correlation time ^[34]. CCSW noises in the gene transcriptional regulatory system has not been reported to date. Therefore, this study focuses on the influence of CCSW noises on the switch process in the gene transcriptional regulatory system.

2. Model and algorithm

Smolen et al. ^[16] have developed a simple kinetic model to examine the complex dynamic activity of genetic regulatory systems. The dynamics of the system are governed by a single ordinary differential equation to determine the concentration of the transcription factor activator (TF-A):

$\frac{{dx}}{{dt}} = \frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - {k_d}x + {R_{bas}},$

(1)

where x(t) denotes the concentration of TF-A. k_d and R_bas represent the decomposition rate and basic synthesis rate of TF-A, respectively. K_d describes the dissociation concentration of the TF-A dimer from the responsive elements, TF-Res. Under the following parameters:

${\left[ { - {{\left( {\frac{{{k_f} + {R_{bas}}}}{{3{k_d}}}} \right)}^3} + \frac{{{K_d}\left( {{k_f} + {R_{bas}}} \right)}}{{6{k_d}}} - \frac{{{K_d}{R_{bas}}}}{{2{k_d}}}} \right]^2} + {\left[ {\frac{{{K_d}}}{3} - {{\left( {\frac{{{k_f} + {R_{bas}}}}{{3{k_d}}}} \right)}^2}} \right]^3} \lt 0,$

(2)

the potential

${U_0}{\rm{(}}x{\rm{)}} = {k_f}\sqrt {{K_d}} {\rm{arctan}}\frac{x}{{\sqrt {{K_d}} }} + \frac{{{k_d}}}{2}{x^2} - {\rm{(}}{R_{bas}}{\rm{ + }}{k_f}{\rm{)}}x.$

(3)

Figure 1(a) shows that when K_d = 10, k_d = 1, R_bas = 0.4, and k_f is in the interval [5.45, 6.68], each value of k_f corresponds to three TF-A, representing two steady states and one unstable state. The corresponding bistable potential U₀ is shown in Figure 1b. When k_f = 6, the stable steady states are x_- ≈ 0.62685 and x₊ ≈ 4.28343 and the unstable steady state is x_u ≈ 1.48971, as shown in Figure 1. To simulate the stochastic effects of the biochemical reaction rates R_bas and k_d, it is assumed that the stochasticity is added to the reaction rates as k_d→k_d + η₁(t) and R_bas→R_bas +η₂(t). Here, multiplicative and additive cross-correlated sine-Wiener noises η₁(t) and η₂(t) are introduced into the original deterministic kinetics. Thus,

$\frac{{dx}}{{dt}} = \frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - \left[ {{k_d} + {\eta _1}\left( t \right)} \right]x + {R_{bas}} + {\eta _2}\left( t \right).$

(4)

Figure 1. (a) Bifurcation plot for the steady state of TF-A on the control parameter of the transcription rate k_f. (b) The bistable potential of Eq (3). k_f = 6, K_d = 10, k_d = 1, and R_bas = 0.4.

DownLoad: Full-Size Img PowerPoint

The explicit definitions of CCSW noises are represented as follows:

${\eta _1}\left( t \right) = A\sin \left( {\sqrt {\frac{2}{{{\tau _1}}}} {\omega _1}\left( t \right)}\right) ,$

(5)

${\eta _2}\left( t \right) = B\sin \left( {\sqrt {\frac{2}{{{\tau _2}}}} {\omega _2}\left( t \right)} \right).$

(6)

The noise strength and self-correlation times of η₁(t) and η₂(t) are indicated as A, B, τ₁ and τ₂. ω₁ and ω₂ are two cross-correlated standard Wiener processes, i.e., $\left\langle {{\omega _1}\left(t \right){\omega _2}\left({t'} \right)} \right\rangle = \left\langle {{\omega _1}\left({t'} \right){\omega _2}\left(t \right)} \right\rangle = \lambda \cdot \min \left({t, t'} \right)$ . $\left\langle { \cdot \cdot \cdot } \right\rangle$ denotes an ensemble average and $\min \left({t, t'} \right)$ denotes taking the smaller value between t and $t'$ . As previously mentioned, the statistical properties of ${\eta _1}\left(t \right)$ and ${\eta _2}\left(t \right)$ are given by:

$\left\langle {{\eta _1}\left( t \right)} \right\rangle = \left\langle {{\eta _2}\left( t \right)} \right\rangle = 0,$

(7)

$\left\langle {{\eta _1}\left( t \right){\eta _1}\left( {t'} \right)} \right\rangle = \frac{{A_{}^2}}{2}\exp \left( { - \frac{{t - t'}}{{{\tau _1}}}} \right)\left[ {1 - \exp \left( { - \frac{{4t'}}{{{\tau _1}}}} \right)} \right],$

(8)

$\left\langle {{\eta _2}\left( t \right){\eta _2}\left( {t'} \right)} \right\rangle = \frac{{B_{}^2}}{2}\exp \left( { - \frac{{t - t'}}{{{\tau _2}}}} \right)\left[ {1 - \exp \left( { - \frac{{4t'}}{{{\tau _2}}}} \right)} \right],$

(9)

$\left\langle {{\eta _1}\left( t \right){\eta _2}\left( {t'} \right)} \right\rangle = \left\langle {{\eta _1}\left( {t'} \right){\eta _2}\left( t \right)} \right\rangle = \lambda \frac{{AB}}{2}\exp \left( { - \frac{{t - t'}}{\tau }} \right)\left[ {1 - \exp \left( { - \frac{{4t'}}{\tau }} \right)} \right], \;t \gt t', {\tau _1} = {\tau _2} = \tau.$

(10)

Here, $\lambda$ and τ stand for the cross-correlation intensity (0 ≤ $\lambda$ ≤ 1) and cross-correlation time between η₁(t) and η₂(t), respectively. The coupling between CCSW noises results in a difficulty of directly integrating the dynamic equation (4). However, CCSW noises can be decoupled and transformed as below ^[32,35,36]:

${\xi _1}\left( t \right) = A\sin \left( {\sqrt {\frac{2}{\tau }} {\omega _\alpha }\left( t \right)} \right),$

(11)

${\xi _2}\left( t \right) = B\lambda \sin \left( {\sqrt {\frac{2}{\tau }} {\omega _\alpha }\left( t \right)} \right) + B\sqrt {1 - {\lambda ^2}} \sin \left( {\sqrt {\frac{2}{\tau }} {\omega _\beta }\left( t \right)} \right).$

(12)

Here, ω_α and ω_β denote two independent standard Wiener processes. Note that the transformations do not change the statistical properties of Eqs 7–10. After the above-mentioned decoupling operations, the dynamic equation (4) is rewritten as follows:

$\frac{{dx}}{{dt}} = \frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - \left[ {{k_d} + {\xi _1}\left( t \right)} \right]x + {R_{bas}} + {\xi _2}\left( t \right).$

(13)

According to the Euler forward procedure, Eq (13) is discretized in our numerical simulations as follows ^[37]:

$\begin{gathered} x\left( {t + \Delta t} \right) = x\left( t \right) + \left[ {\frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - {k_d}x + {R_{bas}}} \right]\;\Delta t \\ - Ax\left( t \right){X_1}\left( t \right)\Delta t + B\left[ {\lambda {X_1}\left( t \right) + \sqrt {1 - {\lambda ^2}} {X_2}\left( t \right)} \right]\;\Delta t \\ - \frac{1}{2}A{X_1}\left( t \right)\;\left\{ {\; - Ax\left( t \right){X_1}\left( t \right) + B\left[ {\lambda {X_1}\left( t \right) + \sqrt {1 - {\lambda ^2}} {X_2}\left( t \right)} \right]\;} \right\}\;\Delta {t^2}, \\ \end{gathered}$

(14)

with

${X_1}\left( t \right) = \sin \left( {\sqrt {{2 \mathord{\left/ {\vphantom {2 \tau }} \right. } \tau }} {\omega _\alpha }\left( t \right)} \right),$

(15)

${X_2}\left( t \right) = \sin \left( {\sqrt {{2 \mathord{\left/ {\vphantom {2 \tau }} \right. } \tau }} {\omega _\beta }\left( t \right)} \right),$

(16)

${\omega _\alpha }\left( t \right) = {\omega _\alpha }\left( {t - \Delta t} \right) + \sqrt { - 2\Delta t\ln {\chi _1}} \cos \left( {2{\rm{ \mathsf{ π} }}{\chi _2}} \right),$

(17)

${\omega _\beta }\left( t \right) = {\omega _\beta }\left( {t - \Delta t} \right) + \sqrt { - 2\Delta t\ln {\chi _3}} \cos \left( {2{\rm{ \mathsf{ π} }}{\chi _4}} \right),$

(18)

Here, ${\chi _1}$ , ${\chi _2}$ , ${\chi _3}$ and ${\chi _4}$ are four independent random numbers uniformly distributed on the unit interval. Δt = 0.001 is the time step.

3. Results and discussion

3.1. Steady state characteristics

We first investigated the effect of the cross-correlation time τ on the steady-state probability distribution (SPD) and switching process. For a given time series of a concentration of TF-A, the total number of time points is N. The observed number n_i of x within the distance bin from x_i to x_i +^△x is counted for a long x (t) sequence that discards transient processing, 0 < x < 10, ^△x = 0.1. The probability n_i/N of each distance is calculated, and thus, the steady-state probability distribution (SPD) is obtained. As shown in Figure 2, the novel phenomenon is that the high concentration state of SPD undergoes a transition from unimodal to bimodal, and therefore, the overall SPD has a three-peak structure with a larger noise cross-correlation time τ. With the increase of the noise cross-correlation time τ, the peak of the probability distribution in the low steady state increases, which corresponds to a decrease of the TF-A monomer concentration and the occurrence of a switch process (i.e., "on"→"off"). That is, the cross-correlation time τ not only causes the on-off switching process but also leads to the formation of double peaks of a high steady state and a triple peaks structure in the overall SPD. Tristability itself may be valuable; for example, a genetic tristate can be found or used in various biochemical networks. This triple peaks structure has only ever been observed in an experiment in synthetic biology by Rui Ma et al. ^[38]. They noted that the triple steady state phenomenon resulted from the discrete and fluctuating nature of a small system. Stem cells and other kinds of multipotent cells might differentiate into two states, and then, one of them redifferentiates into three states according to Waddington's epigenetic landscape ^[39]. Our work also provides a possible mechanism for the formation of the triple peaks structure.

Figure 2. (a) Probability distribution of x(t) for different cross-correlation time τ with A = 0.1, B = 0.1 and λ = 0.5. (b). Sample paths of x(t) with τ = 0.5 and τ = 20. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

Based on the different effects of the noise cross-correlation time τ on the system, the influence of the other parameters of noise is considered when different values of the cross-correlation time τ are selected. First, we consider the effect of multiplicative noise, which is shown in Figure 3. It can be seen that the increase of the multiplicative noise intensity A leads the system to be in a low steady state, corresponding to the decrease of the TF-A monomer concentration and the switch to an off position when τ is smaller and larger. However, multiplicative noise has a strong influence on the bimodal structure of a high steady state. With the increase of the multiplicative noise intensity A, the distance between the two peaks increases. Finally, the left peak of the two peaks weakens and disappears, and the whole peak structure returns to the double peak structure. The extreme value of the peak of the high steady state deviates from the steady-state point to the right.

Figure 3. Probability distribution and sample paths of x(t) for different multiplicative noise intensities A. The additive noise intensity B = 0.2, the cross-correlation intensity λ = 0.3, the cross-correlation intensity τ = 0.1 and τ = 20. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

Then, we investigated the effect of additive noise, which is shown in Figure 4. When τ is small, the increase of the additive noise intensity B leads to a decrease of the low concentration and an increase of the high concentration (Figure 4a). Correspondingly, the time series diagrams of the TF-A monomer concentrations for additive noise intensities B = 0.1 and B = 0.6 are given in Figure 4b. It can be seen that the increase of the additive noise intensity B drives the system to the high steady state, which corresponds to an increase of the TF-A monomer concentration, and the system switches to the "on" position. However, when τ is large, the distance between the two peaks decreases, and they become a single peak with the increase of additive noise B, as shown in Figure 4c-d. At the same time, switching of the process from "on" to "off" is induced.

Figure 4. Probability distribution and sample paths of x(t) for different additive noise intensities B. (a-b) The multiplicative noise intensity A = 0.6, the cross-correlation intensity λ = 0.1, the cross-correlation intensity τ = 0.1. (c-d) the additive noise intensity A = 0.1, the cross-correlation intensity λ = 0.3, the cross-correlation intensity τ = 20. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

Figure 5a shows the effect of the cross-correlation intensity on SPD when τ is small. It can be found that an increase of the noise cross-correlation intensity λ causes the SPD to shift from a high steady state to a low concentration state. However, it is interesting to note that the shift between the different concentration states does not mean the process switches from "on" → "off". As can be seen from Figure 5b, the mean first passage time (MFPT), which is the mean time spent in one state before switching to another, monotonously increases with the increase of the cross-correlation intensity λ. The mean first passage time reflects the switching time between on and off states. That is, the smaller MFPT is, the more frequent the jump between the two steady states. On the contrary, the longer MFPT is, the longer the system maintains the current state. The corresponding time series of the TF-A monomer concentration, as shown in Figure 5c, indicates that a large cross-correlation intensity λ can maintain the TF-A monomer in the current concentration state, that is, keep in the "on" or "off" state. When τ is large, the increase of the cross-correlation intensity λ makes the high-concentration binary states more distinct (Figure 5d-e) and also causes the process to switch from "off" to "on". Contrary to the case of a smaller τ, for a large cross-correlation intensity λ, the system cannot maintain the current state when τ is larger.

Figure 5. (a, d) Probability distribution, (b) MFPT and (c, e) sample paths of x(t) for different cross-correlation intensities λ. (a) A = 0.6, B = 0.2, τ = 0.1. (d) A = 0.1, B = 0.1 and τ = 10. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

In addition, there is feedback delay in the real system, so equation (1) becomes:

$\frac{{dx\left( t \right)}}{{dt}} = \frac{{{k_f}x{{\left( {t - {t_d}} \right)}^2}}}{{x{{\left( {t - {t_d}} \right)}^2} + {K_d}}} - \left[ {{k_d} + {\eta _1}\left( t \right)} \right]x\left( t \right) + {R_{bas}} + {\eta _2}\left( t \right),$

(19)

${t_d}$ is the delay time, and the other parameters are the same as in Eq (1).

Figure 6 shows that the feedback time delay plays an important role in random systems. With the increase of the delay time, the peak of the probability distribution at a low steady state decreases, while those of the high steady state increase, corresponding to an increase of the TF-A monomer concentration, and the system switches from "off" to "on".

Figure 6. Probability distribution and sample paths of x(t) for different delay times t_d. A = 0.8, B = 0.1, λ = 0.5 and τ = 0.1. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

3.2. Mean first passage time (MFPT)

When the system is stochastic bistable, the quantity of interest is the time from one state to another. This time is a random variable, often referred to as the first passage time. Here, we consider the mean first passage time (MFPT) ^[40]. We record the time of the process x(t) to reach the low concentration state x_- from the initial condition x (t = 0) = x₊ (the high concentration state) for the first time and repeat the process 1000 times to obtain the mean first pass time (MFPT). The MFPT can quantify the effect of noise on the transition between two steady states. The main results are shown in Figures 7-10.

Figure 7. The MFPT as a function of the cross-correlation τ (a) for different cross-correlation intensities, λ = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6, B = 0.2; (b) for different multiplicative noise intensities, A = 0.6, 0.7, 0.8, 0.9 with B = 0.2, λ = 0.5; (c) for different additive noise intensities, B = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6, λ = 0.5.

DownLoad: Full-Size Img PowerPoint

Figure 8. The MFPT as a function of the multiplicative noise intensity A. (a) for different cross-correlation intensities, λ = 0.1, 0.3, 0.5, 0.7, 0.9 with B = 0.2 and τ = 0.1; (b) for different additive noise intensities, B = 0.1, 0.3, 0.5, 0.7, 0.9 with λ = 0.5 and τ = 0.1.

DownLoad: Full-Size Img PowerPoint

Figure 9. The MFPT as a function of the additive noise intensity B; (a) for different cross-correlation intensities, λ = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6 and τ = 0.1; (b) for different multiplicative noise intensities, A = 0.6, 0.7, 0.8, 0.9 with λ = 0.5 and τ = 0.1.

DownLoad: Full-Size Img PowerPoint

Figure 10. The MFPT as a function of the cross-correlation intensity λ; (a) for different additive noise intensities, B = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6 and τ = 0.1; (b) for different multiplicative noise intensities, A = 0.6, 0.7, 0.8, 0.9 with λ = 0.5 and τ = 0.1.

DownLoad: Full-Size Img PowerPoint

The MFPT as a function of the cross-correlation time τ for different cross-correlation intensities λ, multiplicative noise intensities A and additive noise intensities B is shown in Figure 7. It can be seen that the MFPT always decreases monotonously with the increase of the cross-correlation time τ. These results indicate that the switch of the process from "on" to "off" becomes easier when the cross-correlation time τ is larger. In addition, the MFPT increases with the increase of the cross-correlation intensity λ, as illustrated in Figures 7a and 10, while the MFPT decreases with the increase of the multiplicative noise intensity A, as shown in Figures 7b and 8. Then, the MFPT varies nonmonotonically with the increase of the additive noise intensity B when λ = 0.5, as shown in Figures 7c, 9 and 10a.

Figure 8 shows the MFPT as a function of multiplicative noise A for different cross-correlation intensities λ and additive noise intensities B. The MFPT always decreases monotonously with the increase of the multiplicative noise intensity A, which is similar to the effect of the cross-correlation time τ. In addition, there are both monotonic and non-monotonic changes with the increase of λ or B when A is fixed in the MFPT as shown in Figure 8a, b. Further results are presented in Figures 9 and 10.

Then, the MFPT as a function of the additive noise intensity B for different cross-correlation intensities and multiplicative noise intensities is shown in Figure 9. The MFPT slightly decreases when the cross-correlation intensity λ is small (λ = 0.1). Moreover, the MFPT first increases and then decreases when λ reaches a moderate value (λ = 0.5). Furthermore, the MFPT increases monotonously and rapidly with a large λ (λ = 0.9). The MFPT showed different trends with the enhancement of B for different values of λ (monotonic decrease, non-monotonic, monotonic increase), which is also shown in Figure 10a. The larger λ is, the more TF-A is able to maintain its current state and the larger MFPT is (Figure 5b). Therefore, it can be inferred that the monotonous decrease of the black dotted line in Figure 9a is due to the role of B when λ is small. This result suggests that the increasing additive noise intensity B also makes the switching process easier. When the value of λ is large, the rapid increase of the MFPT is caused by the enhancement of the ability of λ to maintain the current state. Therefore, when the value of λ is set moderately, the enhancement of B reduces the MFPT, while a moderate value of λ allows the MFPT to be increased, which results in an increase and then a decrease of the MFPT. Moreover, Figure 9b shows that the above nonmonotonicity is destroyed with the increase of A when the value of λ is moderate (λ = 0.5), and the MFPT then undergoes a monotonous increase.

Figure 10 shows the MFPT as a function of the cross-correlation intensity λ for different additive noise intensities B and multiplicative noise intensities A. The MFPT always increases monotonously when an increase of the cross-correlation intensity λ occurs. The results imply that the switch of the process from "on" to "off" becomes difficult. The MFPT decreases monotonically with the increase of B in the region of 0 < λ < 0.1. The MFPT increases monotonically with the increase of B in the region of 0.6 < λ < 1 and presents non-monotonicity in the region of 0.1 < λ < 0.6 (Figure 10a). In Figure 10b, when A increases, the trend of the MFPT increases slowly with the increase of λ. That is, the role of A is opposite to that of λ.

4. Conclusion

In this study, we investigated switch processes in a transcriptional regulatory system driven by CCSW noises. It was found that the increase of the multiplicative noise intensity A and cross-correlation time τ could reduce the concentration of the TF-A monomer and switch to an "off" state. In addition, the effect of the cross-correlation time on SPD presents a novel phenomenon. When the cross-correlation time τ is small, the enhancement of the cross-correlation intensity λ presents a seemingly contradictory result, that is, an increase of the MFPT and an increase of the low steady state, which indicates that the role of the cross-correlation intensity λ is to maintain the current state. The combined effects of the additive noise intensity B and cross-correlation intensity λ show that when λ is small, the effect of maintaining the current state is weak. The enhancement of B leads the SPD to shift from a high steady state to a low concentration one, a switch to the "on" state, and a monotonous decrease of MFPT. Next, when λ is large, the ability to maintain the current state is enhanced and the MFPT increases rapidly. Third, when λ is moderate, the enhancement of B decreases the MFPT and the increase of cross-correlation intensity λ increases the MFPT, which results in the non-monotonicity of the MFPT (first increases and then decreases).

When the cross-correlation time τ is large, the single peak of the high steady state splits into two peaks and the SPD has a three-peak structure. The effect of the cross-correlation intensity λ is different from that of the small cross-correlation time τ. The enhancement of the cross-correlation intensity λ leads to a decrease of the low steady state and an increase of the high steady state, and the bimodal structure of the high steady state is sharper. The SPD is shifted from a high steady state to a low one when the additive noise intensity is decreased. However, the effect of multiplicative noise and additive noise is opposite to that of the double peaks in the high concentration state. The increase of the multiplicative noise intensity leads to an increase of the distance between two peaks in the high concentration state, while the increase of the intensity of additive noise presents an inverse process, that is, the distance between the two peaks decreases until they transition into a single peak. The above results show different characteristics from white Gaussian noise, which lays a foundation for further research on the influence of CCSW noises on nonlinear systems.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No.11675060, 91730301, 11804111), the Huazhong Agricultural University Scientific and the Fundamental Research Funds for the Central University (Grant No. 26622018JC017).

Conflict of Interest

The authors declare that they have no competing interests.

Conflict of interest

The authors declare no conflict of interest.

Author contributions

Wei Cun See Too, Ling Ling Few, Siti Asma' Hassan and Get Bee Yvonne-Τee contributed to the study conception and design. Material preparation, data collection and analysis were performed by Sweta Raikundalia. The first draft of the manuscript was written by Ling Ling Few and Sweta Raikundalia. The funding was procured by Wei Cun See Too, Siti Asma' Hassan, Get Bee Yvonne-Τee and Ling Ling Few. All authors read and approved the final manuscript.

References

[1]	World Health Organization Global Health Estimates 2020: Deaths by Cause, Age, sex, by Country and by Region, 2000−2019, 2023. Available from: https://www.who.int/data/gho/data/themes/mortality-and-global-health-estimates/ghe-leading-causes-of-death
[2]	Sung H, Ferlay J, Siegel RL, et al. (2021) Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 71: 209-249. https://doi.org/10.3322/caac.21660
[3]	Ramírez de Molina A, Gallego-Ortega D, Sarmentero J, et al. (2005) Choline kinase is a novel oncogene that potentiates RhoA-induced carcinogenesis. Cancer Res 65: 5647-5653. https://doi.org/10.1158/0008-5472.CAN-04-4416
[4]	Lacal Sanjuán JC, Zimmerman T, Campos Rosa JM (2021) Choline kjinase: An unexpected journey for a precision medicine strategy in human diseases. Pharmaceutics 13: 788. https://doi.org/10.3390/pharmaceutics13060788
[5]	Rizzo A, Satta A, Garrone G, et al. (2021) Choline kinase alpha impairment overcomes TRAIL resistance in ovarian cancer cells. J Exp Clin Cancer Res 40: 5. https://doi.org/10.1186/s13046-020-01794-6
[6]	Falcon SC, Hudson CS, Huang Y, et al. (2013) A non-catalytic role of choline kinase alpha is important in promoting cancer cell survival. Oncogenesis 2: e38. https://doi.org/10.1038/oncsis.2013.2
[7]	Trousil S, Kaliszczak M, Schug Z, et al. (2016) The novel choline kinase inhibitor ICL-CCIC-0019 reprograms cellular metabolism and inhibits cancer cell growth. Oncotarget 7: 37103-37120. https://doi.org/10.18632/oncotarget.9466
[8]	Sola-Leyva A, López-Cara LC, Ríos-Marco P, et al. (2019) Choline kinase inhibitors EB-3D and EB-3P interferes with lipid homeostasis in HepG2 cells. Sci Rep 9: 5109. https://doi.org/10.1038/s41598-019-40885-z
[9]	Shah V, Shah J (2020) Recent trends in targeting miRNAs for cancer therapy. J Pharm Pharmacol 72: 1732-1749. https://doi.org/10.1111/jphp.13351
[10]	Lan H, Lu H, Wang X, et al. (2015) MicroRNAs as potential biomarkers in cancer: opportunities and challenges. Biomed Res Int 2015: 125094. https://doi.org/10.1155/2015/125094
[11]	Arlauckas SP, Popov AV, Delikatny EJ (2016) Choline kinase alpha—Putting the ChoK-hold on tumor metabolism. Prog Lipid Res 63: 28-40. https://doi.org/10.1016/j.plipres.2016.03.005
[12]	Gruber J, See Too WC, Wong MT, et al. (2012) Balance of human choline kinase isoforms is critical for cell cycle regulation: implications for the development of choline kinase-targeted cancer therapy. FEBS J 279: 1915-1928. https://doi.org/10.1111/j.1742-4658.2012.08573.x
[13]	Ramírez de Molina A, Gutiérrez R, Ramos MA, et al. (2002) Increased choline kinase activity in human breast carcinomas: clinical evidence for a potential novel antitumor strategy. Oncogene 21: 4317-4322. https://doi.org/10.1038/sj.onc.1205556
[14]	Clem BF, Clem AL, Yalcin A, et al. (2011) A novel small molecule antagonist of choline kinase-α that simultaneously suppresses MAPK and PI3K/AKT signaling. Oncogene 30: 3370-3380. https://doi.org/10.1038/onc.2011.51
[15]	Gokhale S, Xie P (2021) ChoK-full of potential: Choline kinase in B cell and T cell malignancies. Pharmaceutics 13: 911. https://doi.org/10.3390/pharmaceutics13060911
[16]	Crilly KS, Tomono M, Kiss Z (1998) The choline kinase inhibitor hemicholinium-3 can inhibit mitogen-induced DNA synthesis independent of its effect on phosphocholine formation. Arch Biochem Biophys 352: 137-143. https://doi.org/10.1006/abbi.1998.0601
[17]	Schiaffino-Ortega S, Baglioni E, Mariotto E, et al. (2016) Design, synthesis, crystallization and biological evaluation of new symmetrical biscationic compounds as selective inhibitors of human Choline Kinase α1 (ChoKα1). Sci Rep 6: 23793. https://doi.org/10.1038/srep23793
[18]	Lin XM, Hu L, Gu J, et al. (2017) Choline kinase α mediates interactions between the epidermal growth factor receptor and mechanistic target of rapamycin complex 2 in hepatocellular carcinoma cells to promote drug resistance and xenograft tumor progression. Gastroenterology 152: 1187-1202. https://doi.org/10.1053/j.gastro.2016.12.033
[19]	Barthe A, Bodaar K, Issaian A, et al. (2022) Choline kinase a is required for alkylating agent resistance in T-cell acute lymphoblastic leukemia. Blood 140: 2258. https://doi.org/10.1182/blood-2022-162508
[20]	Qin SY, Li B, Chen M, et al. (2022) MiR-32-5p promoted epithelial-to-mesenchymal transition of oral squamous cell carcinoma cells via regulating the KLF2/CXCR4 pathway. Kaohsiung J Med Sci 38: 120-128. https://doi.org/10.1002/kjm2.12450
[21]	Zeng S, Liu S, Feng J, et al. (2020) MicroRNA-32 promotes ovarian cancer cell proliferation and motility by targeting SMG1. Oncol Lett 20: 733-741. https://doi.org/10.3892/ol.2020.11624
[22]	Zhang JX, Yang W, Wu JZ, et al. (2021) MicroRNA-32-5p inhibits epithelial-mesenchymal transition and metastasis in lung adenocarcinoma by targeting SMAD family 3. J Cancer 12: 2258-2267. https://doi.org/10.7150/jca.48387
[23]	Liu YJ, Zhou HG, Chen LH, et al. (2019) MiR-32-5p regulates the proliferation and metastasis of cervical cancer cells by targeting HOXB8. Eur Rev Med Pharmacol Sci 23: 87-95. https://doi.org/10.26355/eurrev_201901_16752

This article has been cited by:

1.	P. Dilli Prasad, R. Balamanigandan, 2024, 3193, 0094-243X, 020209, 10.1063/5.0233224
2.	Chanchal Maurya, Vijay Kumar Chaurasiya, Efficient Anonymous Batch Authentication Scheme With Conditional Privacy in the Internet of Vehicles (IoV) Applications, 2023, 24, 1524-9050, 9670, 10.1109/TITS.2023.3271355
3.	Jili Tao, Zejiang Xu, Longhua Ma, Guanzhong Tian, Chengyu Wu, 2024, Neural Network and Reinforcement Learning Based Energy Management Strategy for Battery/Supercapacitor HEV, 979-8-3503-6086-8, 1, 10.1109/ICIEA61579.2024.10665301
4.	Yingjie Xia, Xuejiao Liu, Jing Ou, Oubo Ma, RLID-V: Reinforcement Learning-Based Information Dissemination Policy Generation in VANETs, 2023, 24, 1524-9050, 14151, 10.1109/TITS.2023.3300948
5.	Faten S. Alamri, Khalid Haseeb, Tanzila Saba, Jaime Lloret, Jose M. Jimenez, Multimedia IoT-surveillance optimization model using mobile-edge authentic computing, 2023, 20, 1551-0018, 19174, 10.3934/mbe.2023847
6.	Behnam Ojaghi, Ferran Adelantado, Christos Verikoukis, Enhancing V2X QoS: A Dual Approach of Packet Duplication and Dynamic RAN Slicing in B5G, 2024, 25, 1524-9050, 7848, 10.1109/TITS.2024.3350236
7.	Muhammad Irfan Younas, Muhammad Jawed Iqbal, Abdul Aziz, Ali Hassan Sodhro, Toward QoS Monitoring in IoT Edge Devices Driven Healthcare—A Systematic Literature Review, 2023, 23, 1424-8220, 8885, 10.3390/s23218885
8.	Marcel Vološin, Eugen Šlapak, Zdenek Becvar, Taras Maksymyuk, Adam Petík, Madhusanka Liyanage, Juraj Gazda, Blockchain-Based Route Selection With Allocation of Radio and Computing Resources for Connected Autonomous Vehicles, 2023, 24, 1524-9050, 7230, 10.1109/TITS.2023.3255301
9.	Zhi Cai, Qian Dong, Meihui Shi, Xing Su, Limin Guo, Zhiming Ding, VOLTCom: A Novel Online Trajectory Compression Method Based on Vector Processing, 2023, 24, 1524-9050, 14982, 10.1109/TITS.2023.3294941
10.	Jingyi Zhang, Jili Tao, Yuanmin Hu, Longhua Ma, An Energy Management Strategy Based on DDPG With Improved Exploration for Battery/Supercapacitor Hybrid Electric Vehicle, 2024, 25, 1524-9050, 3999, 10.1109/TITS.2023.3327105
11.	Ruiqi Liu, Meng Hua, Ke Guan, Xiping Wang, Leyi Zhang, Tianqi Mao, Di Zhang, Qingqing Wu, Abbas Jamalipour, 6G Enabled Advanced Transportation Systems, 2024, 25, 1524-9050, 10564, 10.1109/TITS.2024.3362515
12.	Yazhou Zhang, Prayag Tiwari, Qian Zheng, Abdulmotaleb El Saddik, M. Shamim Hossain, A Multimodal Coupled Graph Attention Network for Joint Traffic Event Detection and Sentiment Classification, 2023, 24, 1524-9050, 8542, 10.1109/TITS.2022.3205477
13.	Muhammad Ayzed Mirza, Yu Junsheng, Salman Raza, Manzoor Ahmed, Muhammad Asif, Azeem Irshad, Neeraj Kumar, MCLA Task Offloading Framework for 5G-NR-V2X-Based Heterogeneous VECNs, 2023, 24, 1524-9050, 14329, 10.1109/TITS.2023.3292140
14.	Hana Elhachi, Farouk Boumehrez, Mohamed Aymen Labiod, Salah Redadaa, Abdelhamid Mellouk, Smart cross‐layer approach to multi‐access terrestrial and non‐terrestrial networks (NTNs): Real‐time mobile‐health use case, 2024, 37, 1074-5351, 10.1002/dac.5941
15.	Liangliang Lou, Mingxin Song, Xinquan Chen, Xiaoming Zhao, Shiqing Zhang, Optimized Wireless Sensing and Deep Learning for Enhanced Human-Vehicle Recognition, 2024, 25, 1524-9050, 7508, 10.1109/TITS.2024.3352820
16.	Tao Wang, Jiaping Liao, Jing Jin, Robust Decoupled Cruise Control Design of the High-Speed Train as Uncertain Large-Scale System, 2024, 25, 1524-9050, 16096, 10.1109/TITS.2024.3443287
17.	Debadreeta Das, Manoj Kumar Rana, Bhaskar Sardar, Tommaso Pecorella, Debashis Saha, A Comparative Analysis of Distributed Mobility Management Schemes for 5G-Based Intelligent Transportation Systems, 2025, 26, 1524-9050, 2434, 10.1109/TITS.2024.3506662
18.	Sonika Bhardwaj, Ramesh Saha, 2024, A Comprehensive Review of Routing Protocols in Vehicular Ad Hoc Networks (VANETs): Challenges, Solutions and Future Directions, 979-8-3503-7871-9, 1, 10.1109/UPCON62832.2024.10982762

Reader Comments

Your name:*

Email:*
© 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)