
Citation: Grujicic Mica, Ramaswami S., S. Snipes J., Yavari R.. Multi-Scale Computation-Based Design of Nano-Segregated Polyurea for Maximum Shockwave-Mitigation Performance[J]. AIMS Materials Science, 2014, 1(1): 15-27. doi: 10.3934/matersci.2014.1.15
[1] | Shubhra Singh, Douglas J. H. Shyu . Perspective on utilization of Bacillus species as plant probiotics for different crops in adverse conditions. AIMS Microbiology, 2024, 10(1): 220-238. doi: 10.3934/microbiol.2024011 |
[2] | Esther Menendez, Paula Garcia-Fraile . Plant probiotic bacteria: solutions to feed the world. AIMS Microbiology, 2017, 3(3): 502-524. doi: 10.3934/microbiol.2017.3.502 |
[3] | Esther Menendez, Paula Garcia-Fraile . Plant probiotic bacteria: solutions to feed the world. AIMS Microbiology, 2017, 3(4): 747-748. doi: 10.3934/microbiol.2017.4.747 |
[4] | Alexandra Díez-Méndez, Raul Rivas . Improvement of saffron production using Curtobacterium herbarum as a bioinoculant under greenhouse conditions. AIMS Microbiology, 2017, 3(3): 354-364. doi: 10.3934/microbiol.2017.3.354 |
[5] | Paola Ortiz-Ojeda, Katty Ogata-Gutiérrez, Doris Zúñiga-Dávila . Evaluation of plant growth promoting activity and heavy metal tolerance of psychrotrophic bacteria associated with maca (Lepidium meyenii Walp.) rhizosphere. AIMS Microbiology, 2017, 3(2): 279-292. doi: 10.3934/microbiol.2017.2.279 |
[6] | Luciana Porto de Souza Vandenberghe, Lina Marcela Blandon Garcia, Cristine Rodrigues, Marcela Cândido Camara, Gilberto Vinícius de Melo Pereira, Juliana de Oliveira, Carlos Ricardo Soccol . Potential applications of plant probiotic microorganisms in agriculture and forestry. AIMS Microbiology, 2017, 3(3): 629-648. doi: 10.3934/microbiol.2017.3.629 |
[7] | Alejandro Jiménez-Gómez, Lorena Celador-Lera, María Fradejas-Bayón, Raúl Rivas . Plant probiotic bacteria enhance the quality of fruit and horticultural crops. AIMS Microbiology, 2017, 3(3): 483-501. doi: 10.3934/microbiol.2017.3.483 |
[8] | Srirengaraj Vijayaram, Reshma Sinha, Caterina Faggio, Einar Ringø, Chi-Chung Chou . Biopolymer encapsulation for improved probiotic delivery: Advancements and challenges. AIMS Microbiology, 2024, 10(4): 986-1023. doi: 10.3934/microbiol.2024043 |
[9] | Ananya Mukherjee, Puja Bhattacharjee, Rituparna Das, Arundhati Pal, Amal K. Paul . Endophytic bacteria with plant growth promoting abilities from Ophioglossum reticulatum L.. AIMS Microbiology, 2017, 3(3): 596-612. doi: 10.3934/microbiol.2017.3.596 |
[10] | Sabrine Balti, Yassine Mabrouk, Mouna Souihi, Imen Hemissi, Ismail Amri, Ethan Humm, Noor Khan, Ann M. Hirsch . Combined inoculation of rhizobacteria with Mesorhizobium promotes growth, nutrient contents, and protects chickpea against Fusarium redolens. AIMS Microbiology, 2025, 11(2): 318-337. doi: 10.3934/microbiol.2025015 |
Mammalian cell division patterns are critical to understanding human tumor progression. It has drawn the attention of many researchers, and it has been the topic of intense research for a long time. Several research works utilize age-structured frameworks to investigate the cell cycle. Some examples of age-structured growth models include epidemic [1,2,3], microscopic virus [4,5] and cell population [6,7,8,9] models. The concealed molecular complexity of a tissue, on the other hand, demand a more thorough modeling framework that includes special cellular and molecular interactions.
Cells that divide in living tissues can be divided into two categories: proliferating and quiescent. The proliferating cells go through different phases in the cell cycle ($ G_1, S, G_2, M $) while dividing. Quiescent cells, however, do not grow or divide, but rather move to the $ G_0 $ phase, where they remain until they differentiate or undergo apoptosis. To preserve tissue homeostasis, cells must be capable to transition between proliferative and quiescent phases. The transition between proliferating and quiescent compartments, however, is dependent on signaling molecules called anti-growth or growth factors [10]. In a population of tumor cells, proliferating cells multiply till the tumor becomes active and aggressive. Additionally, to maintain homeostasis, the total number of cells in all sub-populations stays at equilibrium; hence, in a healthy and mutated cell populations, the proliferative compartment is constrained in size. Schematics in Figure 1 depicts multiscale-modeling framework used in this paper.
Mutations in genes that regulate cell growth and division can disturb the finely-tuned equilibrium that governs cell proliferation, which in turn can result in the emergence of cancer. Although other factors like environmental exposures and lifestyle choices may also play a role in cancer development, mutations in genes are a major contributor to this disease. Mutations can be of several types, those that are particularly significant for cancer involve increased potential for proliferation, decreased apoptosis, genetic instability, and reduced tumor suppression [11]. Furthermore, studies have shown that the transformation of a normal cell into a cancerous one usually requires the accumulation of one to ten mutations [11,12]. The mutated cell population also consists of its own progeny including different differentiation stages and quiescent cells.
This research focuses on the model development of the cell population in all healthy and mutated cell populations. We investigate the coupling dynamics of tissue cell density and cell cycle proteins. Previous studies have used age-structured models to investigate cell populations in quiescent phase [6] only, the proliferating phase [13,14] only, or both phases together [7,15,16,17,18,19]. Despite this, the impact of molecular interactions on the interplay between proliferative and quiescent phases at the subcellular level has yet to be examined. Thus, the primary aim of this paper is to develop a multiscale model utilizing mathematical techniques that can capture the intricacies of a complex system existing at the sub-cellular level.
The emphasis is laid upon the interaction between the macro- (population dynamics) and micro-scale (cell cycle dynamics) at two distinct levels. The concept of age is used to refer to the time since the last division, as described in previous studies [14,20]. In age-structured models, an additional variable, referred to as "age" ($ a $), is introduced, which has a physiological rather than physical interpretation. The idea of cell age reflects the variability (biological) within a population of dividing cells and is distinct from the physical time variable ($ t $). To model the behavior of cell populations in both the proliferative and quiescent stages at the macroscale, we use partial differential equations (PDEs). For predicting sub-cellular protein interactions related to cell cycle dynamics, we employ ordinary differential equations (ODEs). The two scales are connected via the feedback incorporated in both directions. Within a cycle of cell division ($ G_1 $, $ S $, $ G_2 $, $ M $), cells in the early proliferating phase ($ \mathrm{G}_1 $) can move to the quiescent-phase until the restriction point ($ \mathrm{R} $).
It is clear that the restriction point ($ \mathrm{R} $) cut the $ \mathrm{G}_1 $ phase into two portions, where before $ \mathrm{R} $, the cells become quiescent but once $ \mathrm{R} $ is passed, cells can no longer avoid division [21,22]. During the quiescent phase, cells do not grow or divide, but they still complete other cellular functions. The bidirectional transition of cells between quiescent and proliferation phases in both healthy and mutated cell populations plays a crucial role in maintaining tissue homeostasis, and it is controlled by the conditions of extracellular environmental [10]. In tumors, the transitional balance between two compartments is disturbed, and cells can grow uncontrollably [23]. Recent experimental results support that cyclin proteins are the most important to regulate any change in the cell-cycle phase [24]. Therefore, we utilize a vital aspect of cell cycle dynamics (i.e., the $ G_1-S $ phase-transition) to predict how the transitional balance between proliferating and quiescent subpopulations evolve, which is essential for maintaining homeostasis.
Several proteins are involved in regulating the transitions between different phases of the cell cycle at the microscale. The interactions between these proteins have been mathematically modeled and simulated using ODEs by various researchers, as documented in [25,26,27,28,29,30] and related references. However, to keep things simple, we focus only on four proteins: Cyclin $ \mathsf{D-CDK 4/6} $, $ \mathsf{p21} $, $ \mathsf{E2F} $, and $ \mathsf{Rb} $. These proteins are mainly involved in Cyclin $ \mathsf{D} $ activity and the transition of cells from $ G_1 $ phase to $ S $ phase. Experimental data supports this choice, as Cyclin $ \mathsf{D} $ has been shown to regulate the transition between the $ G_0 $ and $ G_1 $ phases [31,32,33]. In addition, over-expression of Cyclin $ \mathsf{D} $ leads to cell division commitment in the proliferative phase, while under-expression leads to the quiescent phase. It is essential to note that these molecular-interactions are considered to occur in a rapidly growing cell population rather than a single cell. Furthermore, we assume the average protein concentrations in both quiescent and proliferating cell subpopulations, thus omitting the variability between individual cells. The proteins involved in the cell cycle play a vital role in regulating its progression. Cyclin proteins, which are structural proteins, along with their inhibitors cyclin-dependent kinases ($ \mathsf{CDK} $), control the various phases of the cell cycle. Each phase of the cycle is regulated by a specific Cyclin/$ \mathsf{CDK} $ complex.
The initiation of a cell cycle is triggered when the microenvironment of a cell receives sufficient growth signals [34], which induce the activity of Cyclin-$ \mathsf{CDK} $ 4-6 complex. During the $ G_1 $ phase, Cyclin $ \mathsf{D} $ is activated only by growth signals, and in case of low or no growth factors, its concentration decreases, and the cell doesn't enter the cycle. Growth factors bind to specific receptors which are situated on the cell's external cytoplasmic membrane, activating intracellular signaling pathways (such as Ras/Map/Rap kinase), eventually which directing to the production of Cyclin $ \mathsf{D} $ (see [35,36,37], for more details). The active complex of Cyclin $ \mathsf{D} $ and $ \mathsf{CDK 4/6} $ is synthesized at maximum. Cyclin complex then triggers the transcription factor $ \mathsf{E2F} $'s activation by phosphorylating retinoblastoma protein $ \mathsf{Rb} $. Consequently, the transcription factor $ \mathsf{E2F} $ accumulates and triggers some other important cyclins. It is important to note that these interactions occur in a fast growing population of cells rather than in a single cell.
In summary, our contribution is to develop a multiscale model that describes the coupling between two predominant scales and focuses on the challenges related to the disruption of cell transitioning between proliferating and quiescent states, leading to uncontrolled tumor growth. Specifically, we investigate whether Cyclin $ \mathsf{D} $ complex is one of the important cause in creating a deregulation in cell transitioning between proliferating and quiescent cells.
The structure of rest of this paper is as follows. In Section 2, we delve into the details of multiscale mathematical modeling of proliferating and quiescent cells with regards to the dynamics of cell cycle. In Section 3, we investigate the existence and uniqueness of non-negative solutions utilizing spectral and semigroup theory. Moving on to Section 4, we start with deriving steady-state solutions and then establish spectral criteria for their local stability. Specifically, we show that if the linearised semigroup's growth bound is negative, the steady-state solution is locally asymptotically stable, while a positive growth bound implies instability of the steady-state solution. Lastly, Sections 5 and 6 offer results discussion and the concluding remarks of this paper, respectively.
The cell densities of healthy and mutated cells in the proliferative and quiescent compartments are described by nonlinear hyperbolic transport PDEs that relate the cell density distribution to both physiological age $ a $ and time $ t $. Specifically, the densities of healthy cells in the proliferative and quiescent phases are expressed as follows:
$ ∂∂tph(a,t)+∂∂a(gh(a)ph(a,t))=γh(N)qh(a,t)−(τh(a)+αh(a,x1)+μph(a))ph(a,t), $ | (2.1) |
$ ∂∂tqh(a,t)=αh(a,x1)ph(a,t)−(γh(N)+μqh(a))qh(a,t), $ | (2.2) |
where the rate evolution of a cell cycle is denoted by $ g_h(a) $ in the equation. The first term on the right side, $ \gamma_h(N)q_h(a, t) $, represents the transition to proliferating from quiescent cells, while the term $ \tau_h(a)p_h(a, t) $ represents the cell densities for completing cell division in some age of the proliferating phase. Cells which move to the quiescent phase without undergoing division are represented by $ \alpha_h(a, x_1)p_h(a, t) $. The loss in the proliferating cells due to apoptosis/necrosis is represented by the death rate $ \mu_{p_h}(a) $. The inflow from healthy proliferating cells, regulated by the microscale variable of Cyclin complex concentration $ x_1 $ for each age, is denoted by the first term in Eq (2.2): $ \alpha_h(a, x_1)p_h(a, t) $. The next term depicts a loss in the quiescent cells due to either by returning to proliferating phase at the rate $ \gamma_h(N) $ or by cell death due to apoptosis (or necrosis), as represented by the death rate $ \mu_{q_h}(a) $.
Next, the cell density of mutated cells in mutated proliferating $ (p_c) $ and quiescent $ (q_c) $ phases, respectively, is presented:
$ ∂∂tpc(a,t)+∂∂a(gc(a)pc(a,t))=γc(N)qc(a,t)−(τc(a)+αc(a,x1)+μpc(a))pc(a,t), $ | (2.3) |
$ ∂∂tqc(a,t)=αc(a,x1)pc(a,t)−(γc(N)+μqc(a))qc(a,t), $ | (2.4) |
where the terms used are similar to those in the case of healthier cells, as shown in Eqs (2.1) and (2.2). The total cell number in both healthy and mutated populations of cells in quiescent and proliferating phases is denoted by $ N(t) $ and is defined in Eq (2.5). In the case of quiescent cells, aging does not occur (i.e., the cells stop aging), so the convection term related to physiological age $ a $ is absent in Eqs (2.2) and (2.4). The total number of cells, represented by $ N(t) $, represents the sum of all cells in the proliferating and quiescent phases throughout all ages, and can be expressed as
$ N(t)=∫a⋆0(ph(a,t)+qh(a,t)+pc(a,t)+qc(a,t))da, $ | (2.5) |
where maximum age of the cells is given by $ a^\star $. The initial conditions are given below:
$ ph(a,0)=ph,0(a),qh(a,0)=qh,0(a),pc(a,0)=pc,0(a),qc(a,0)=qc,0(a),∀a≥0. $ | (2.6) |
The boundary conditions are given as follows:
$ gh(0)ph(0,t)=2(1−m)∫a⋆0τh(a)ph(a,t)da, $ | (2.7) |
$ gc(0)pc(0,t)=2∫a⋆0τc(a)pc(a,t)da+2m∫a⋆0τh(a)ph(a,t)da, $ | (2.8) |
for $ t > 0 $, where the number 2 shows the two newborn cells initializing in the proliferating phase, and the parameter $ m $ represents the mutation rate. Since healthy cell can acquire a mutation only during a division process, therefore, new born mutated cells start will start at age 0.
The function $ \tau_i(a) $ represents the cell number that finish dividing at a particular age in both healthy and mutated proliferating phases. Here, the index $ i $ indicates whether the compartment is healthy or cancerous, denoted by $ h $ and $ c $, respectively. The function $ \tau_i(a) $ is regulated by the age of the cell, denoted by $ a $, and is almost zero until a minimum cell age. Afterward, the function increases until it reaches the age of $ a^* $.
$ τi(a)=ρ1,iaγ1,iργ1,i2,i+aγ1,i, $ | (2.9) |
The maximum proliferation rate is represented by $ \rho_{1, i} $, while $ \rho_{2, i} $ is the age t achieve half-maximum. The Hill coefficient is represented by the exponent $ \gamma_{1, i} $.
Next, we establish the rate at which cells transition to the quiescent phase from the proliferating compartments, which depends on both the age of the cell ($ a $) and the quantity of the Cyclin $ \mathsf{D-CDK 4/6} $ complex ($ x_1 $).
$ αi(a,x1)=σ1,iσγ2,i2,i(σγ2,i2,i+xγ2,i1)σγ3,i3,i(σγ3,i3,i+aγ3,i). $ | (2.10) |
The function $ \alpha_i(a, x_1) $ depicts the number non-dividing cells due to anti-growth factors. The age-dependence of $ \alpha_i $ is motivated by the fact that cells transition to the quiescent phase from the proliferating phase only until they reach a specific age that marks a restriction point ($ \mathrm{R} $) in the cell cycle (which is also $ G_1-S $ phase transition). However, before the restriction point, the Cyclin complex's concentration $ x_1 $ must be below a certain value to enable cells to exit the proliferating phase. In Eq (2.10), the Hill coefficients are represented by $ \gamma_{2, i} $ and $ \gamma_{3, i} $, while $ \sigma_{2, i} $ and $ \sigma_{3, i} $ denote the concentration of the Cyclin $ \mathsf{D-CDK 4/6} $ complex $ x_1 $ and the age $ a $, respectively. After $ \gamma_{2, i} $ and $ \gamma_{3, i} $, the rate function $ \alpha $ decreases asymptotically to zero, preventing cells from transitioning to the quiescent phase. This implies that at age $ \sigma_{3, i} $, cells are inevitably committed to entering the proliferation phase. Finally, $ \sigma_{3, i} $ represents the threshold concentration of the Cyclin complex to determine the restriction point $ \mathrm{R} $.
The function $ \gamma_i(N) $, which determines the number of cells transitioning to the proliferating phase from the quiescent phase, is represented by a Hill function of $ N $ that decreases monotonically:
$ γi(N)=νiθκiiθκii+Nκi, $ | (2.11) |
where the Hill function are defined as follows: $ \nu_i $ specifies the maximum rate at which cells transition to proliferating from quiescent population, when there are no cells, i.e., $ N = 0 $; $ \kappa_i $ is the Hill coefficient, and $ \theta_i $ represents the proportion of the total cell population that reaches half the maximum value of $ \nu_i $. This implies that the number of quiescent cells transitioning to the proliferative compartment decreases to zero as the cell population increases, illustrating density inhibition.
Cell growth is controlled by proteins such as cytokines and other factors that regulate proliferation [38]. Cytokines bind to specific receptors, activating signaling pathways [39]. The cytokine signals that regulate cell numbers are reliant on the total population of cells, as demonstrated by various studies [40]. For a detailed explanation of cytokine signal dynamics, refer to [41,42]. Using the quasi-steady-state approximation, we can express the quantity of growth factors ($ g_f $) produced by the entire cell population ($ N $) as
$ gf=11+ktN. $ | (2.12) |
We restrict ourselves to only four microscale proteins (states of the model) in our cell cycle model, as mentioned earlier, that are sufficient to account for the reversible transitions between the quiescent and proliferating phases.We utilize Michaelis-Menten kinetics to depict the chemical reactions occurring between enzymes and substrates during the cell cycle. The basics of this approach are briefly outlined below. When adequate growth factors are present, Cyclin $ \mathsf{D} $ protein combines with its catalytic partner CDK4-6, resulting in the formation of Cyclin $ \mathsf{D-CDK 4/6} $ complex. As outlined in [31,43], this complex is responsible for phosphorylating various proteins involved in the cell cycle. These phosphorylated proteins are essential for advancing through the initial growth phase of the cell cycle and crossing the restriction point $ \mathrm{R} $. To provide further details, the Cyclin $ \mathsf{D-CDK 4/6} $ complex induces phosphorylation of the retinoblastoma protein $ \mathsf{Rb} $, resulting in its deactivation and subsequent release of the transcription factor $ \mathsf{E2F} $. This activates several signals that promote cell growth, facilitating the progression of the cell cycle. In contrast, $ \mathsf{p21} $ regulates the cell cycle by inhibiting the functions of various $ \mathsf{CDK} $ proteins. These proteins are described in Table 1.
Description | State |
Cyclin $ \mathsf{D-CDK 4/6} $ | $ x_1 $ |
$ \mathsf{E2F} $ | $ x_2 $ |
$ \mathsf{Rb} $ | $ x_3 $ |
$ \mathsf{p21} $ | $ x_4 $ |
In this study, we examine the behavior of a single cell to represent the dynamics of all cells within a population. Assuming that all cells behave similarly, we utilize a single ordinary differential equation (ODE) model with similar parameters for all cells in the population to describe the underlying cell cycle dynamics on a microscale. We also account for cells with shorter cycles on a macroscale using the function $ \tau_i(a) $. Our model considers a specific age $ a^\star $ at which the representative cell completes division. The cell cycle dynamics are described by the ODE system presented in [44].
$ dx1da=k1s(gfkgf+gf)−k14x4x1−k1d(x1k1+x1), $ | (2.13) |
$ dx2da=k21(x2t−x2k2+(x2t−x2))x1−k32x2x3−k2dx2, $ | (2.14) |
$ dx3da=k3s−k32x2x3−k31(x3k3+x3)x1−k3dx3, $ | (2.15) |
$ dx4da=k4s+k42(k34k34+x3)x2−k41(x4k4+x4)x1−k4dx4. $ | (2.16) |
A detailed description of all the terms and parameters involved in the model equations (2.13)–(2.16) can be found in [45]. Although we will not delve into the full derivation of these equations here, inquisitive readers can refer to [44] for a comprehensive explanation. To aid in understanding, we have included simulations of the four microscale states mentioned earlier in Figure 2.
This section presents the uniqueness of the solution to the initial-boundary value problem (2.1)–(2.7) and (2.13)–(2.16), which we will simplify by using the microscale model for the entire time $ t $, instead only until age $ a $. We introduce Banach spaces, $ X = L^1(0, a^\star)\times L^1(0, a^\star)\times L^1(0, a^\star)\times L^1(0, a^\star) $ and $ Y = L^1(0, a^\star)\times L^1(0, a^\star)\times L^1(0, a^\star)\times L^1(0, a^\star) $, with the norm $ |\phi| = \sum_{i = 1}^4|\phi_i|1 $ for $ \phi(a) = (\phi_1(a), \phi_2(a), \phi_3(a), \phi_4(a))^\mathrm{T}\in X $ and $ |\varphi| = \sum{i = 1}^4|\varphi_i|1 $ for $ \varphi(a) = (\varphi_1(a), \varphi_2(a), \varphi_3(a), \varphi_4(a))^\mathrm{T}\in Y $, where $ |\cdot|1 $ is the standard norm of $ L^1(0, a^\star) $. We first treat the initial-boundary value problem of system (2.1)–(2.7) as an abstract Cauchy problem on Banach space $ X $. We assume that $ g{h_a}, g{h_{aa}}, g_{c_a}, g_{c_{aa}}\in L^{\infty}((0, a^\star)\times\mathbb{R}^+) $, and non-negative death rates, that is, $ \mu_{p_h}(\cdot) = \mu_{q_h}(\cdot)\geq0 $ and $ \mu_{p_c}(\cdot) = \mu_{q_c}(\cdot)\geq0 $, and are locally integrable on $ [0, a^\star) $. The transition rate $ \alpha_i(a, x_1)\in L^{\infty}((0, a^\star)\times(0, a^\star)) $, and $ \tau_i(a)\in L^{1}(0, a^\star) $. We start by defining a linear operator $ \mathsf{A}_1 $ as follows:
$ (A1ϕ)(a)=(−∂(gh(a)ϕ1(a))∂a−(τh(a)+μph(a))ϕ1(a)−μqh(a)ϕ2(a)−∂(gc(a)ϕ3(a))∂a−(τc(a)+μpc(a))ϕ3(a)−μqc(a)ϕ4(a)),ϕ(a)=(ϕ1(a),ϕ2(a),ϕ3(a),ϕ4(a))T∈D(A1). $ |
The symbol $ \mathrm{T} $ denotes the transpose of the vector, and the domain $ \mathrm{D}(\mathsf{A}_1) $ is given by the following:
$ D(A1)={(ϕ1,ϕ2,ϕ3,ϕ4)|ϕi is absolute continuous on [0,a⋆),ϕ(0)=(2(1−m)∫a⋆0τh(a)ϕ1(a)da,0,2∫a⋆0τc(a)ϕ3(a)da+2m∫a⋆0τh(a)ϕ1(a)da,0)T}. $ |
The nonlinear operator $ \mathsf{F}_1:X \times Y\to X $ is given by
$ (F1(ϕ,φ))(a)=(νhθκhhϕ2(a)θκhh+(Nϕ)κh−αh(φ1,a)ϕ1(a)−νhθκhhϕ2(a)θκhh+(Nϕ)κh+αh(φ1,a)ϕ1(a)νcθκccϕ3(a)θκcc+(Nϕ)κc−αc(φ1,a)ϕ4(a)−νcθκccϕ3(a)θκcc+(Nϕ)κc+αc(φ1,a)ϕ4(a)),ϕ∈X,φ∈Y, $ |
where the linear operator $ N $ on $ L^1(0, a^\star)\times L^1(0, a^\star)\times L^1(0, a^\star)\times L^1(0, a^\star) $ is given by
$ Nϕ=∫a⋆0(ϕ1(a)+ϕ2(a)+ϕ3(a)+ϕ4(a))da. $ |
Consider $ \upsilon(t) = (p_h(\cdot, t), q_h(\cdot, t), p_c(\cdot, t), q_c(\cdot, t))^\mathrm{T}\in X $. We can define the initial-boundary value problem (2.1)–(2.7) as an abstract semilinear initial value problem (IVP) in $ X $, as shown below:
$ ddtυ(t)=A1υ(t)+F1(υ(t),v(t)),υ(0)=υ0∈X, $ | (3.1) |
where $ \upsilon_0(a) = (p_{h_0}(a), q_{h_0}(a), p_{c_0}(a), q_{c_0}(a)) $.
Next, we define IVP (2.13) as a Cauchy problem on the Banach space $ Y $. Suppose $ \mathsf{A}_2 $ is a linear operator which reads
$ (A2φ)(a)=(0−k2dφ2(a)k3s−k3dφ3(a)k4s−k4dφ4(a)),φ(a)=(φ1(a),φ2(a),φ3(a),φ4(a))T∈D(A2), $ |
where the domain $ \mathrm{D}(\mathsf{A}_2) $ is
$ \mathrm{D}(\mathsf{A}_2) = \{\varphi\in Y\, |\, \varphi_i\, \text{ is absolute continuous on }[0, a^\star), \, \, \varphi(0) = (0, 0, 0, 0)^\mathrm{T}\}. $ |
We define the nonlinear operator $ \mathsf{F}_2:X \times Y\to Y $ by
$ (F2(ϕ,φ))(a)=(k1s(gf(Nϕ)kgf+gf(Nϕ))−k14φ4(a)φ1(a)−k1d(φ1(a)k1+φ1(a)),k21(x2t−φ2(a)k2+(x2t−φ2(a)))φ1(a)−k32φ2(a)φ3(a)−k32φ2(a)φ3(a)−k31(φ3(a)k3+φ3(a))φ1(a)k42(k34k34+φ3(a))φ2(a)−k41(φ4(a)k4+φ4(a))φ1(a)), $ |
where $ \phi\in X, \, \varphi\in Y $. Take $ v(t) = (x_1(t), x_2(t), x_3(t), x_4(t))^\mathrm{T}\in Y $. Then (2.13)–(2.16) can be expressed as an abstract semilinear IVP in $ Y $:
$ ddtv(t)=A2v(t)+F2(υ(t),v(t)),v(0)=v0∈Y, $ | (3.2) |
where $ v_0(t) = (x_1^0, x_2^0, x_3^0, x_4^0) $, we can now establish a joint Cauchy problem for (3.1) and (3.2) as shown below:
$ ddt(υv)=(A100A2)(υv)+(F1(υ,v)F2(υ,v)),(υ(0)v(0))=(υ0v0)∈Z, $ |
$ ddtζ(t)=Aζ(t)+F(ζ(t)),ζ(0)=ζ0∈Z, $ | (3.3) |
where $ \zeta = (\upsilon, v) $, $ \zeta_0 = (\upsilon_0, v_0) $, $ \mathsf{A} = \left(A100A2\right) $, $ \mathsf{F} = \left(F1F2\right) $ and the Banach space is $ Z = \{X, Y\} $. Assuming that $ T(t) $ is a $ C_0 $-semigroup generated by $ \mathsf{A} $ for $ t\geq0 $, and the operator $ \mathsf{F} $ is continuously Fréchet differentiable on $ Z $ (specifically, both $ \mathsf{F}_1 $ and $ \mathsf{F}_2 $ are Fréchet differentiable on $ X $ and $ Y $), a continuous mild solution $ t\to \zeta(t, \zeta_0) $ exists and is unique for each $ \zeta_0\in Z $ on a maximal interval $ [0, t_1) $ in $ Z $.
$ ζ(t,ζ0)=T(t)ζ0+∫t0T(t−s)F(ζ(s,ζ0))ds,∀t∈[0,t1), $ | (3.4) |
where $ t_1 $ can be either $ +\infty $ or $ \lim_{t\to t_1^-}|\zeta(t, \zeta_0)| = \infty $. Moreover, if $ \zeta_0\in \mathrm{D}(\mathsf{A}) $, then $ \zeta(t, \zeta_0)\in \mathrm{D}(\mathsf{A}) $ for $ 0\leq t < t_1 $, and the function $ \zeta\to \zeta(t, \zeta_0) $ is continuously differentiable and satisfies (3.3) on $ [0, t_1) $. This has been established in Proposition 4.16 [46,47].
Remark 3.1. Let's take $ p_{h, \max} $, $ q_{h, \max} $, $ p_{c, \max} $, $ q_{c, \max} $, $ x_{1, \max} $, $ x_{2, \max} $, $ x_{3, \max} $ and $ x_{4, \max} $ to represent the maximum values of the solution variables. If we normalise the governing equations using $ N(a) = p_h(a, t)+q_h(a, t)+p_c(a, t)+q_c(a, t)+x_1(a)+x_2(a)+x_3(a)+x_4(a) $, then an a-priori estimate on these would lead to $ p_h(a, t)+q_h(a, t)+p_c(a, t)+q_c(a, t)+x_1(a)+x_2(a)+x_3(a)+x_4(a) = 1 $.
Lemma 3.1. Let $ \Omega = \{(p_h, q_h, p_c, q_c, x_1, x_2, x_3, x_4)\in Z|p_h\geq0, q_h\geq0, p_c\geq0, q_c\geq0, x_1\geq0, x_2\geq0, x_3\geq0, x_4\geq0\} $ and let $ \Omega_0 = \{(p_h, q_h, p_c, q_c, x_1, x_2, x_3, x_4)\in$ $ Z|0\leq p_h\leq p_{h, \max}, 0\leq q_h\leq q_{h, \max}, 0\leq p_c\leq p_{c, \max}, 0\leq q_c\leq q_{c, \max},$ $ 0\leq x_1\leq x_{1, \max}, 0\leq x_2\leq x_{2, \max}, 0\leq x_3\leq x_{3, \max}, 0\leq x_4\leq x_{4, \max}\} $. Then after a finite time, the mild solution $ \zeta(t, \zeta_0) $ of (3.3), where $ \zeta_0\in \Omega $, enters a positively invariant set $ \Omega_0 $.
Proof. To obtain the solution of Eq (2.1), we will begin by utilizing transformations $ \tilde{p}_h(a, t) = g_h(a)p_h(a, t) $ and $ \tilde{q}_h(a, t) = g_h(a)q_h(a, t) $ for $ t\in[0, t_1] $ and $ a\in[a_0, a^\star) $. Then, for $ t\in(0, t_1) $ and $ a\in(a_0, a^\star) $, we have from Eq (2.1)
$ ∂˜ph(a,t)∂t+gh(a)∂˜ph(a,t)∂a=γh(N(t))˜qh(a,t)−(τh(a)+αh(x1(a),a)+μph(a))˜ph(a,t), $ | (3.5) |
Next, we apply the parameter transform to remove the term $ g_h(a) $ and define a new age variable $ \eta $ for both $ p_h $ and $ q_h $, Lemma 3.1 [41]. This yields the expression:
$ \frac{\partial}{\partial \eta} \tilde{p}_h(a(\eta), t) = \frac{ \mathrm{d} a}{ \mathrm{d} \eta}\frac{\partial}{\partial a} \tilde{p}_h(a, t) = g_h(a)\frac{\partial}{\partial a} \tilde{p}_h(a, t), \:\: \text{where}\: \frac{ \mathrm{d} a}{ \mathrm{d} \eta} = g_h(a). $ |
Therefore, from Eq (3.5), it follows that
$ ∂˜ph(a(η),t)∂t+∂˜ph(a(η),t)∂η=γh(N(t))˜qh(a(η),t)−(τh(a(η))+αh(x1(a(η)),a(η))+μph(a(η)))˜ph(a(η),t). $ | (3.6) |
To obtain the explicit relation of $ \tilde{p}_h(a(\eta), t) $, we will utilize the method of characteristics (MOC). Specifically, we assume that $ \tilde{p}_h(a(\eta), t) $ is governed by an ODE along the curve $ (a(\psi_1(y)), \psi_2(y)) = \psi(y) $, and therefore, we have
$ ˙ψ1(y):=1⇒ψ1(y)=y+c1,˙ψ2(y):=1⇒ψ2(y)=y+c2,z(y):=˜ph(a(ψ1(y)),ψ2(y)), $ |
where $ c_1, c_2\in\mathbb{R} $ are constants. Then, it follows
$ dzdy=d˜ph(a(ψ1(y)),ψ2(y))dy=∂˜ph(a(ψ1(y)),ψ2(y))∂ada(ψ1(y))dψ1dψ1(y)dy+∂˜ph(a(ψ1(y)),ψ2(y))∂ψ2dψ2(y)dy=γh(N(ψ2(y)))˜qh(a(ψ1(y)),ψ2(y))−(τh(a(ψ1(y)))+αh(x1(a(ψ1(y))),a(ψ1(y)))+μph(a(ψ1(y))))˜ph(a(ψ1(y)),ψ2(y))=γh(N(ψ2(y)))˜qh(a(ψ1(y)),ψ2(y))−(τh(a(ψ1(y)))+αh(x1(a(ψ1(y))),a(ψ1(y)))+μph(a(ψ1(y))))z(y). $ | (3.7) |
We can now write $ \tilde{p}_h $ using an ODE (3.7) so that
$ ˜ph(a(y+c1),y+c2)=˜ph(a(ψ1(y)),ψ2(y))=z(y)=exp(−∫y0(τh(a(ψ1(ξ)))+αh(x1(a(ψ1(ξ))),a(ψ1(ξ)))+μph(a(ψ1(ξ))))dξ)[∫y0exp(∫ζ0(τh(a(ψ1(ξ)))+αh(x1(a(ψ1(ξ))),a(ψ1(ξ)))+μph(a(ψ1(ξ))))dξ)γh(N(ψ2(ζ)))˜qh(a(ψ1(ζ)),ψ2(ζ))dζ+˜ph(a(ψ1(0)),ψ2(0))]=exp(−∫y0(τh(a(ξ+c1))+αh(x1(a(ψ1(ξ+c1))),a(ξ+c1))+μph(a(ξ+c1)))dξ)[∫y0exp(∫ζ0(τh(a(ξ+c1))+αh(x1(a(ψ1(ξ+c1))),a(ξ+c1))+μph(a(ξ+c1)))dξ)γh(N(ζ+c2))˜qh(a(ζ+c1),ζ+c2)dζ+˜ph(a(c1),c2)]. $ |
Now, we define the boundary set $ \Gamma $ as $ {[a_0, a^\star)\times{0}}\cup{{0}\times[0, t_1]} $, which enables us to use the boundary condition to determine $ \tilde{p}_h(a(c_1), c_2) $ if a curve $ (a(\psi_1(y)), \psi_2(y)) $ begins in $ \Gamma $. In order for $ (a(y+c_1), y+c_2) $ to lie on $ \Gamma $, either $ c_1 = 0 $ or $ c_2 = 0 $. Therefore, we have the following two scenarios:
In the first scenario, we can randomly choose $ c_1 = 0 $ and $ c_2\in[0, t_1) $. In this case, we have
$ ˜ph(a(y),y+c2)=exp(−∫y0(τh(a(ξ))+αh(x1(a(ξ)),a(ξ))+μph(a(ξ)))dξ)[∫y0exp(∫ζ0(τh(a(ξ))+αh(x1(a(ξ)),a(ξ))+μph(a(ξ)))dξ)γh(N(ζ+c2))˜qh(a(ζ),ζ+c2)dζ+˜ph(a(0),c2)]. $ |
The solution in $ {(a(\eta), t)|t\in[0, t_1], \eta\in[0, \min(\eta^*, t))} $ can be obtained using the characteristic solution as follows:
$ \eta \stackrel{{{ ! }}}{ = }\psi_1(y) = y+c_1 = y\Rightarrow y = \eta\, \, \text{and}\, \, t \stackrel{{{ ! }}}{ = }\psi_2(y) = y+c_2\Rightarrow c_2 = t-y, $ |
which implies
$ ˜ph(a(η),t)=exp(−∫η0(τh(a(ξ))+αh(x1(a(ξ)),a(ξ))+μph(a(ξ)))dξ)[∫η0exp(∫ζ0(τh(a(ξ))+αh(x1(a(ξ)),a(ξ))+μph(a(ξ)))dξ)γh(N(ζ+t−η))˜qh(a(ζ),ζ+t−η)dζ+˜ph(a(0),t−η)]. $ |
Using the above equation, we can obtain the expression for $ g_h(a(\eta))p_h(a(\eta), t) $ when $ \eta < t $. By choosing an arbitrary $ c_1\in[0, \eta^*) $ and setting $ c_2 = 0 $, we obtain
$ ˜ph(a(y+c1),u)=exp(−∫y0(τh(a(ξ+c1))+αh(x1(a(ξ+c1)),a(ξ+c1))+μph(a(ξ+c1)))dξ)[∫y0exp(∫ζ0(τh(a(ξ+c1))+αh(x1(a(ξ+c1)),a(ξ+c1))μph(a(ξ+c1)))dξ)+γh(N(ζ))˜qh(a(ζ+c1),ζ)dζ+˜ph(a(c1),0)]. $ |
Using the characteristic solution, we can obtain a solution in the set $ {(a(\eta), t)|t\in[0, t_1], \eta\in[t, \eta^*)} $ as follows:
$ \eta \stackrel{{{ ! }}}{ = }\psi_1(y) = y+c_1\Rightarrow c_1 = \eta-y\, \, \text{and}\, \, t \stackrel{{{ ! }}}{ = }\psi_2(y) = y+c_2\Rightarrow y = t, $ |
which results into
$ ˜ph(a(η),t)=exp(−∫t0(τh(a(ξ+η−t))+αh(x1(a(ξ+η−t)),a(ξ+η−t))+μph(a(ξ+η−t)))dξ)[∫t0exp(∫ζ0(τh(a(ξ+η−t))+αh(x1(a(ξ+η−t)),a(ξ+η−t))+μph(a(ξ+η−t)))dξ)γh(N(ζ))˜qh(a(ζ+η−t),ζ)dζ+˜ph(a(η−t),0)]. $ |
Hence, the relation for $ g_h(a(\eta))p_h(a(\eta), t) $ is now established for $ \eta > t $. As a result, the ultimate solution for $ g_h(a(\eta))p_h(a(\eta), t) $ can be expressed as
$ ˜ph(a(η),t):={exp(−∫η0(τh(a(ξ))+αh(x1(a(ξ)),a(ξ))+μph(a(ξ)))dξ)[h(t−η)∫η0exp(∫ζ0(τh(a(ξ))+αh(x1(a(ξ)),a(ξ))+μph(a(ξ)))dξ)γh(N(ζ+t−η))˜qh(a(ζ),ζ+t−η)dζ],ˉa<t,exp(−∫t0(τh(a(ξ+η−t))+αh(x1(a(ξ+η−t)),a(ξ+η−t))+μph(a(ξ+η−t)))dξ)[p0(a(η−t))+∫t0exp(∫ζ0(τh(a(ξ+η−t))+αh(x1(a(ξ+η−t)),a(ξ+η−t))+μph(a(ξ+η−t)))dξ)γh(N(ζ))˜qh(a(ζ+η−t),ζ)dζ],ˉa≥t, $ |
where the boundary condition $ \tilde{p}_h(a(0), t-\eta) $ is denoted as $ h(t-\eta) $. Note that for positive initial data, the above expression is positive and for $ g_h(a)q_h(a, t)\geq0 $.
We then derive the solution expression from (2.2) as shown below:
$ qh(a,t):=exp(−∫t0μqh(a)+γh(N(t))dt){∫t0exp(−∫ξ0μqh(a)+γh(N(π))dπ)αh(x1(a),a)ph(a,ξ)dξ+qh,0(a)}. $ | (3.8) |
As a direct consequence, we observe that $ q_h(a, t) $ is non-negative for positive initial data and whenever $ g_h(a)q_h(a, t)\geq0 $. Similarly, we can obtain the solution expression for $ p_c $ and $ q_c $.
Next, to ensure the positivity of the coupled ODE model (2.13), we express the system of ODEs as follows:
$ {dx1da=F1(x1,x4),dx2da=F2(x1,x2,x3),dx3da=F3(x1,x2,x3),dx4da=F4(x1,x2,x3,x4), $ | (3.9) |
where $ \mathsf{F}_1 $, $ \mathsf{F}_2 $, $ \mathsf{F}_3 $ and $ \mathsf{F}_4 $ correspond to the vector fields of the microscale states $ x_1 $–$ x_4 $. It is worth noting that in (3.9), $ \mathsf{F}_1 $ does not depend on $ N $ (i.e., $ p_h $, $ q_h $, $ p_c $ and $ q_c $), as $ N $ changes with time and is a fixed constant at each time step, which determines the growth factors entire age range.
To ensure that the solutions of all ODEs is positive, it is essential to verify that the vector fields $ \mathsf{F}_1, \mathsf{F}_2, \mathsf{F}_3, \mathsf{F}_4 $ are smoothly differentiable and oriented in a direction that points away from the negative regions in the state space. Starting with the ODE for $ x_1 $ from (3.9), we set $ x_4 = 0 $ in $ \mathsf{F}_1(x_1, x_4) $ to obtain $ \dot{x}_1 = \mathsf{F}_1(x_1) $. We can observe that $ \mathsf{F}_1(x_1) = k_{1s}\bigg(\frac{g_f}{k_{gf}+g_f}\bigg)-k_{1d}\bigg(\frac{x_1}{k_{1}+x_1}\bigg) > 0 $ for all $ a > 0 $ when $ k_{1s}\bigg(\frac{g_f}{k_{gf}+g_f}\bigg) > k_{1d}\bigg(\frac{x_1}{k_{1}+x_1}\bigg) $. This implies that the concentration of $ x_1 $ consistently rises more than it falls over time, as the sole source of an increase in $ x_1 $ concentration is from growth factors. Consequently, during periods when growth factors are at their minimum, the concentration of $ x_1 $ is also at its minimum, meaning the amount of degradation or decrement cannot surpass the activation of the $ x_1 $ complex. Given that the solution to system (2.13)–(2.16) is unique for each initial condition, as can be observed from (3.3) and (3.4), we can infer that the solution remains in the first quadrant for any $ x_4 > 0 $. Therefore, the positivity of the solution for $ x_1 $ is guaranteed. To obtain an ODE for $ \dot{x}_2 $, we assume $ x_1 = 0 $ in $ \mathsf{F}_2(x_1, x_2, x_3) $. This results in $ \dot{x}_2 = \mathsf{F}_2(x_2, x_3) $, whose solution takes the form $ x_2(a) = x_2^0e^{-(k{32}x_3(a)-k{2d})a} $. Hence, for any positive initial data, $ x_2(a) $ remains positive for all ages and $ x_3(a) $ values. Similarly, we substitute $ x_3 = 0 $ in $ \mathsf{F}_2(x_1, x_2, x_3) $ to obtain a nonlinear ODE $ \dot{x}_2 = \mathsf{F}_2(x_1, x_2) $ for $ x_2 $. Although an explicit solution cannot be computed, the phase portrait of $ (x_1, x_2) $ reveals that the solution trajectories move away from the axis separating the positive and negative space for positive initial data. By following a similar procedure, we can establish sufficient conditions for the positivity of solutions for $ x_3(a) $ and $ x_4(a) $. Therefore, we conclude that if $ \zeta_0\in\Omega $, then $ \zeta(t, \zeta_0)\in\Omega $ for all $ t > 0 $.
The above analysis implies that the local solution $ \zeta(t, \zeta_0) $ of (3.3) with initial conditions $ \zeta_0\in \mathrm{D}(\mathsf{A})\cap \Omega $ has a well-defined and finite norm. Consequently, we obtain our final result.
Theorem 3.1. The abstract Cauchy problem (3.3) has a unique global classical solution on $ Z $ with respect to the initial data $ z_0\in\Omega\cap D(\mathsf{A}) $.
As a result of having positive initial data, the IVP (2.1)–(2.4) possesses a singular positive solution.
This section aims to determine the steady-state solution of the model and to present sufficient conditions for the existence of a positive steady-state. To this end, we specify some notation. Let $ X $ be a real or complex Banach space, and let $ X^\star $ denote its dual space. We denote the value of $ F\in X^\star $ at $ \psi\in X $ as $ \langle F, \psi\rangle $. Additionally, we define a cone $ X_+ $ as a non-zero set that satisfies $ X_+\cap(-X_+) = {0} $, $ \lambda X_+\subset X_+ $ for $ \lambda\geq0 $, and $ X_++X_+\subset X_+ $. Furthermore, we define the dual cone, denoted as $ X^\star_+ $, as the subset of the dual space.
The steady-states of the system (2.1)–(2.4) and (2.13)–(2.16) are denoted by $ \bar{p}_h(a) $, $ \bar{q}_h(a) $, $ \bar{p}_c(a) $, $ \bar{q}_c(a) $ and $ \bar{x}_1-\bar{x}_4 $. These steady-states must satisfy the following set of time-invariant ordinary differential equations:
$ ∂(gh(a)ˉph(a))∂a=ˉγhˉqh(a)−(τh(a)+αh(a,ˉx1)+μph(a))ˉph(a), $ | (4.1) |
$ 0=αh(a,ˉx1)ˉph(a)−(ˉγh+μqh(a))ˉqh(a), $ | (4.2) |
$ ∂(gc(a)ˉpc(a))∂a=ˉγcˉqc(a)−(τc(a)+αc(a,ˉx1)+μpc(a))ˉpc(a), $ | (4.3) |
$ 0=αc(a,ˉx1)ˉpc(a)−(ˉγc+μqc(a))ˉqc(a), $ | (4.4) |
$ ˉph(0)=2(1−m)∫a⋆0τh(a)ˉph(a)da, $ | (4.5) |
$ ˉpc(0)=2∫a⋆0τc(a)ˉpc(a)da+2∫a⋆0τh(a)ˉph(a)da, $ | (4.6) |
$ dˉx1da=k1s(ˉgfkgf+ˉgf)−k14ˉx4ˉx1−k1d(ˉx1k1+ˉx1), $ | (4.7) |
$ dˉx2da=k21(x2t−ˉx2k2+(x2t−ˉx2))ˉx1−k32ˉx2ˉx3−k2dˉx2, $ | (4.8) |
$ dˉx3da=k3s−k32ˉx2ˉx3−k31(ˉx3k3+ˉx3)ˉx2−k3dˉx3, $ | (4.9) |
$ dˉx4da=k4s+k42(k34k34+ˉx3)ˉx2−k41(ˉx4k4+ˉx4)ˉx1−k4dˉx4, $ | (4.10) |
where $ \bar{\gamma}_i = \gamma_i(\bar{N}) $, $ \bar{g}_f = g_f(\bar{N}) $ and $ \bar{N} = \int_0^{a^\star}(\bar{p}_h(a)+\bar{q}_h(a)+\bar{p}_c(a)+\bar{q}_c(a)) \mathrm{d} a $. Since the cell cycle model's ODEs are age-dependent and the system is in a steady-state due to the input of growth factors, all cell cycle states attain a steady-state. As a result, we can determine the steady-states of the quiescent and proliferating cell populations, represented by $ \bar{p}_h(a) $, $ \bar{q}_h(a) $, $ \bar{p}_c(a) $ and $ \bar{q}_c(a) $, without explicitly solving the equations of microscale model. Solving the system (4.10) for $ \bar{p}_h $, $ \bar{q}_h $, $ \bar{p}_c $ and $ \bar{q}_c $ allows us to obtain the values of $ \bar{q}_h $ and $ \bar{q}_c $:
$ ˉqh(a)=αh(a,ˉx1)ˉph(a)ˉγh+μqh(a),ˉqc(a)=αc(a,ˉx1)ˉpc(a)ˉγc+μqc(a), $ | (4.11) |
and substituting the aforementioned expressions for $ \bar{q}_h $ and $ \bar{q}_c $ into the equations for $ \bar{p}_h $ and $ \bar{p}_c $, respectively, results in the following expressions:
$ d(gh(a)ˉph(a))da+(αh(a,ˉx1)μqh(a)ˉγh+μqh(a)+τh(a)+μph(a))ˉph(a)=0, $ | (4.12) |
$ d(gc(a)ˉpc(a))da+(αc(a,ˉx1)μqc(a)ˉγc+μqc(a)+τc(a)+μpc(a))ˉpc(a)=0. $ | (4.13) |
Solving Eq (4.12) for $ \bar{p}_h(a) $ and $ \bar{p}_c(a) $, yields steady-state solutions for $ \bar{p}_h(a) $, $ \bar{q}_h(a) $, $ \bar{p}_c(a) $ and $ \bar{q}_c(a) $ as follows:
$ {ˉph(a)=ˉph(0)exp(−∫a01gh(a)(g′h(a)+αh(ˉx1,ξ)μqh(ξ)ˉγh+μqh(ξ)+τh(ξ)+μph(ξ))dξ),ˉqh(a)=αh(a,ˉx1)ˉph(0)ˉγh+μqh(a)exp(−∫a01gh(a)(g′h(a)+αh(ˉx1,ξ)μq,h(ξ)ˉγh+μqh(ξ)+τh(ξ)+μph(ξ))dξ),ˉpc(a)=ˉpc(0)exp(−∫a01gc(a)(g′c(a)+αc(ˉx1,ξ)μqc(ξ)ˉγc+μqc(ξ)+τc(ξ)+μpc(ξ))dξ),ˉqc(a)=αc(a,ˉx1)ˉpc(0)ˉγc+μqc(a)exp(−∫a01gc(a)(g′c(a)+αc(ˉx1,ξ)μq,c(ξ)ˉγc+μqc(ξ)+τc(ξ)+μpc(ξ))dξ). $ |
It is evident that the system described in Eqs (2.1)–(2.4), (2.13) always has a trivial steady-state.
Our next objective is to obtain the stability criteria for a positive steady-state solution. Suppose $ p_h(a, t) = \bar{p}_h $, $ q_h(a, t) = \bar{q}_h $, $ p_c(a, t) = \bar{p}_c $, $ q_c(a, t) = \bar{q}_c $, $ \forall t\geq0 $ represent equilibrium solutions to the PDE model (2.1)–(2.4) and $ p_h^*(a, t) $, $ q_h^*(a, t) $, $ p_c^*(a, t) $ and $ q_c^*(a, t) $ represent the corresponding perturbation terms:
$ p_h(a, t) = \bar{p}_h+\epsilon p_h^*(a, t), \: q_h(a, t) = \bar{q}_h+\epsilon q_h^*(a, t), \: p_c(a, t) = \bar{p}_c+\epsilon p_c^*(a, t), \: q_c(a, t) = \bar{q}_c+\epsilon q_c^*(a, t). $ |
After substituting the aforementioned expressions into the PDE model (2.1)–(2.4), we obtain
$ {ϵ∂∂tp∗h(a,t)+∂∂a(gh(a)(ˉph+ϵp∗h(a,t)))=(νhθκhhθκhh+(ˉN+ϵn(t))κh)(ˉph+ϵq∗h(a,t))−(τh(a)+αh(a,ˉx1)+μph(a))(ˉph+ϵp∗h(a,t)),ϵ∂∂tq∗h(a,t)=αh(a,ˉx1)(ˉph+ϵp∗h(a,t))−(νhθκhhθκhh+(ˉN+ϵn(t))κh+μqh(a))(ˉph+ϵq∗h(a,t)),ϵ∂∂tp∗c(a,t)+∂∂a(gc(a)(ˉpc+ϵp∗c(a,t)))=(νcθκccθκcc+(ˉN+ϵn(t))κc)(ˉpc+ϵq∗c(a,t))−(τc(a)+αc(a,ˉx1)+μpc(a))(ˉpc+ϵp∗c(a,t)),ϵ∂∂tq∗c(a,t)=αc(a,ˉx1)(ˉpc+ϵp∗c(a,t))−(νcθκhcθκcc+(ˉN+ϵn(t))κc+μqc(a))(ˉpc+ϵq∗c(a,t)),(ˉph(0)+ϵp∗h(0,t))=2(1−m)∫a⋆0τh(a)(ˉph+ϵp∗h(a,t))da,(ˉpc(0)+ϵp∗c(0,t))=2∫a⋆0τc(a)(ˉpc+ϵp∗c(a,t))da+2m∫a⋆0τh(a)(ˉph+ϵp∗h(a,t))da, $ |
where $ n(t): = \int_0^{a^\star}\left(p_h^*(a, t)+q_h^*(a, t)+p_c^*(a, t)+q_c^*(a, t)\right) \mathrm{d} a $. Then, take the derivative with respect to epsilon $ \epsilon $, leads to
$ {∂∂tp∗h(a,t)+∂∂a(gh(a)p∗h(a,t))=∂∂ϵ(νhθκhhϵθκhh+(ˉN+ϵn(t))κh)q∗h(a,t)−(τh(a)+αh(a,ˉx1)+μph(a))p∗h(a,t),∂∂tq∗h(a,t)=αh(a,ˉx1)p∗h(a,t)−(∂∂ϵ(νhθκhhϵθκhh+(ˉN+ϵn(t))κh)−μqh(a))q∗h(a,t),∂∂tp∗c(a,t)+∂∂a(gc(a)p∗c(a,t))=∂∂ϵ(νcθκccϵθκcc+(ˉN+ϵn(t))κc)q∗c(a,t)−(τc(a)+αc(a,ˉx1)+μpc(a))p∗c(a,t),∂∂tq∗c(a,t)=αc(a,ˉx1)p∗c(a,t)−(∂∂ϵ(νcθκccϵθκcc+(ˉN+ϵn(t))κc)−μqc(a))q∗c(a,t),p∗h(0,t)=2(1−m)∫a⋆0τh(a)p∗h(a,t)da,p∗c(0,t)=2∫a⋆0τc(a)p∗c(a,t)da+2m∫a⋆0τh(a)p∗h(a,t)da, $ |
which simplifies to
$ {∂∂tp∗h(a,t)+∂∂a(gh(a)p∗h(a,t))=νhθκhh(θκhh+(ˉN+ϵn(t))κh−κhϵn(t)(ˉN+ϵn(t))κh−1(θκhh+(ˉN+ϵn(t))κh)2)q∗h(a,t)−(αh(a,ˉx1)+τh(a)+μph(a))p∗h(a,t),∂∂tq∗h(a,t)=αh(a,ˉx1)p∗h(a,t)−(νhθκhh(θκhh+(ˉN+ϵn(t))κh−κhϵn(t)(ˉN+ϵn(t))κh−1(θκhh+(ˉN+ϵn(t))κh)2)−μqh(a))q∗h(a,t),∂∂tp∗c(a,t)+∂∂a(gc(a)p∗c(a,t))=νcθκcc(θκcc+(ˉN+ϵn(t))κc−κcϵn(t)(ˉN+ϵn(t))κc−1(θκcc+(ˉN+ϵn(t))κc)2)q∗c(a,t)−(αc(a,ˉx1)+τc(a)+μpc(a))p∗c(a,t),∂∂tq∗c(a,t)=αc(a,ˉx1)p∗c(a,t)−(νcθκcc(θκcc+(ˉN+ϵn(t))κc−κcϵn(t)(ˉN+ϵn(t))κc−1(θκcc+(ˉN+ϵn(t))κc)2)−μqc(a))q∗c(a,t),p∗h(0,t)=2(1−m)∫a⋆0τh(a)p∗h(a,t)da,p∗c(0,t)=2∫a⋆0τc(a)p∗c(a,t)da+2m∫a⋆0τh(a)p∗h(a,t)da. $ |
In the limit as $ \epsilon $ approaches zero, we arrive at a linear system of partial differential equations:
$ {∂∂tp∗h(a,t)+∂∂a(gh(a)p∗h(a,t))=γh(ˉN)q∗h(a,t))−(αh(a,ˉx1)+τh(a)+μph(a))p∗h(a,t),∂∂tq∗h(a,t)=αh(a,ˉx1)p∗h(a,t)−(μqh(a)+γh(ˉN))q∗h(a,t),∂∂tp∗c(a,t)+∂∂a(gc(a)p∗c(a,t))=γc(ˉN)q∗c(a,t))−(αc(a,ˉx1)+τc(a)+μpc(a))p∗c(a,t),∂∂tq∗c(a,t)=αc(a,ˉx1)p∗c(a,t)−(μqc(a)+γc(ˉN))q∗c(a,t),p∗h(0,t)=2(1−m)∫a⋆0τh(a)p∗h(a,t)da,p∗c(0,t)=2∫a⋆0τc(a)p∗c(a,t)da+2m∫a⋆0τh(a)p∗h(a,t)da, $ | (4.14) |
where $ \gamma_i(\bar{N}) = \nu_i\theta_i^{\kappa_i}/ (\theta_i^{\kappa_i}+\bar{N}^{\kappa_i}) $, where $ i = \{h, c\} $. Next, we formulate (4.14) as semilinear problem:
$ ddtω(t)=Cω(t),ω(0)=ω0∈X, $ | (4.15) |
where the generator $ \mathrm{C} $ is defined on the Banach space $ X $ as follows:
$ (Cϕ)(a)=(−(∂∂a+1gh(a)(τh(a)+αh(a,ˉx1)+μph(a)))gh(a)ϕ1(a)+γh(ˉN)ϕ2(a)αh(a,ˉx1)ϕ1(a)−(γh(ˉN)+μqh(a))ϕ2(a)−(∂∂a+1gc(a)(τc(a)+αc(a,ˉx1)+μpc(a)))gc(a)ϕ1(a)+γc(ˉN)ϕ2(a)αc(a,ˉx1)ϕ1(a)−(γc(ˉN)+μqc(a))ϕ2(a)), $ |
where
$ \phi(a) = (\phi_1(a), \phi_2(a), \phi_3(a), \phi_4(a))^\mathrm{T}\in \mathrm{D}(\mathrm{C}), $ |
where $ \mathrm{D}(\mathrm{C}) $ is defined below:
$ D(C)={(ϕ1,ϕ2)|ϕi is absolute continuous on [0,a⋆),ϕ(0)=(2(1−m)∫a⋆0τh(a)ϕ1(a)da,0,2∫a⋆0τc(a)ϕ2(a)da+2m∫a⋆0τh(a)ϕ1(a)da,0)T}. $ |
Next, we take the resolvent equation for the operator $ \mathrm{C} $:
$ (λI−C)ϕ=ψ,ϕ∈D(C),ψ∈X,λ∈C, $ | (4.16) |
which leads to
$ −γh(ˉN)ϕ2(a)+∂∂a(gh(a)ϕ1(a))+(λ+τh(a)+αh(a,ˉx1)+μph(a))ϕ1(a)=ψ1(a), $ | (4.17) |
$ (λ+γh(ˉN)+μqh(a))ϕ2(a)−αh(a,ˉx1)ϕ1(a)=ψ2(a), $ | (4.18) |
$ −γc(ˉN)ϕ4(a)+∂∂a(gc(a)ϕ3(a))+(λ+τc(a)+αc(a,ˉx1)+μpc(a))ϕ3(a)=ψ3(a), $ | (4.19) |
$ (λ+γc(ˉN)+μqc(a))ϕ4(a)−αc(a,ˉx1)ϕ3(a)=ψ4(a), $ | (4.20) |
and
$ \phi_1(0) = 2(1-m)\int_0^{a^\star}\tau_h(a)\phi_1(a) \mathrm{d} a, \:\:\:\phi_3(0) = 2\int_0^{a^\star}\tau_c(a)\phi_3(a) \mathrm{d} a+2m\int_0^{a^\star}\tau_h(a)\phi_1(a) \mathrm{d} a. $ |
By solving (4.18) and (4.20), we get
$ ϕ2(a)=ψ2(a)+αh(a,ˉx1)ϕ1(a)λ+γh(ˉN)+μqh(a),ϕ4(a)=ψ4(a)+αc(a,ˉx1)ϕ3(a)λ+γc(ˉN)+μqc(a), $ | (4.21) |
which after substituting in Eqs. (4.17) and (4.19) gives
$ ϕ1(a)=exp(−∫a0τh(ξ)+αh(ˉx1,ξ)+λ+μph(ξ)−γh(ˉN)αh(ˉx1,ξ)gh(ξ)(λ+γh(ˉN)+μqh(ξ))dξ)[∫a0exp(∫ζ0τh(ξ)+αh(ˉx1,ξ)+λ+μph(ξ)−γh(ˉN)αh(ˉx1,ξ)gh(ξ)(λ+γh(ˉN)+μqh(ξ))dξ)1gh(ζ){ψ1(ζ)+γh(ˉN)ψ2(ζ)λ+γh(ˉN)+μqh(ζ)}dζ+ϕ1(0)],ϕ3(a)=exp(−∫a0τc(ξ)+αc(ˉx1,ξ)+λ+μpc(ξ)−γc(ˉN)αc(ˉx1,ξ)gc(ξ)(λ+γc(ˉN)+μqc(ξ))dξ)[∫a0exp(∫ζ0τc(ξ)+αc(ˉx1,ξ)+λ+μpc(ξ)−γc(ˉN)αc(ˉx1,ξ)gc(ξ)(λ+γc(ˉN)+μqc(ξ))dξ)1gc(ζ){ψ3(ζ)+γc(ˉN)ψ4(ζ)λ+γc(ˉN)+μqc(ζ)}dζ+ϕ3(0)]. $ |
Substituting $ \phi_1(a) $ and $ \phi_3(a) $ back in Eq (4.21) yields
$ ϕ2(a)=1λ+γh(ˉN)+μqh(a)[exp(−∫a0τh(ξ)+αh(ˉx1,ξ)+λ+μph(ξ)−γh(ˉN)αh(ˉx1,ξ)gh(ξ)(λ+γh(ˉN)+μqh(ξ))dξ){∫a0exp(∫ζ0τh(ξ)+αh(ˉx1,ξ)+λ+μph(ξ)−γh(ˉN)αh(ˉx1,ξ)gh(ξ)(λ+γh(ˉN)+μqh(ξ))dξ)1gh(ζ){ψ1(ζ)+γh(ˉN)ψ2(ζ)λ+γh(ˉN)+μqh(ζ)}dζ+ϕ1(0)}αh(a,ˉx1)+ψ2(a)],ϕ4(a)=1λ+γc(ˉN)+μqc(a)[exp(−∫a0τc(ξ)+αc(ˉx1,ξ)+λ+μpc(ξ)−γc(ˉN)αc(ˉx1,ξ)gc(ξ)(λ+γc(ˉN)+μqc(ξ))dξ){∫a0exp(∫ζ0τc(ξ)+αc(ˉx1,ξ)+λ+μpc(ξ)−γc(ˉN)αc(ˉx1,ξ)gc(ξ)(λ+γc(ˉN)+μqc(ξ))dξ)1gc(ζ){ψ3(ζ)+γc(ˉN)ψ4(ζ)λ+γc(ˉN)+μqc(ζ)}dζ+ϕ3(0)}αc(a,ˉx1)+ψ4(a)]. $ |
Lemma 4.1. The resolvent of operator $ \mathrm{C} $ is compact and its spectrum, denoted by $ \sigma(\mathrm{C}) $, satisfies the condition:
$ σ(C)=σP(C)={λ∈C|1∈σp(Uλ)}. $ | (4.22) |
Here, $ \sigma_P(\mathrm{C}) $ refers to the point spectrum of $ \mathrm{C} $, and $ U_\lambda $ is an operator dependent on $ \lambda $.
Proof. The expression of $ \phi_1(a) $ and $ \phi_3(a) $ can be re-written as
$ ϕ1(a)=1αh(a,ˉx1){(λ+γh(ˉN)+μqh)(Uh,λψ2)(a)+γh(ˉN)(Uh,λψ1)(a)},ϕ3(a)=1αc(a,ˉx1){(λ+γc(ˉN)+μqc)(Uc,λψ4)(a)+γc(ˉN)(Uc,λψ3)(a)}, $ |
where the linear operator on Banach space, $ U_{i, \lambda} $ is given as
$ (Ui,λψ)(a)=∫a⋆0Hi,λ(ζ,a)ψ(ζ)dζ,i={h,c}, $ | (4.23) |
where
$ Hλ(ζ,a)=αi(a,ˉx1)gi(ζ)(λ+γi(ˉN)+μqi(a))exp(−∫a0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)g(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)exp(∫ζ0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)gi(ξ)(λ+γi(ˉN)+μqi(ξ))dξ). $ | (4.24) |
Similarly, we rewrite $ \phi_2(a) $ and $ \phi_4(a) $ as
$ ϕ2(a)=(Ui,λψ2)(a)+(Vi,λψ1)(a),ϕ4(a)=(Ui,λψ4)(a)+(Vi,λψ3)(a), $ |
where the linear operator on Banach space, $ V_{i, \lambda} $ is given as
$ (Vi,λψ)(a)=∫a⋆0Gi,λ(ζ,a)ψ(ζ)dζ,Gi,λ(ζ,a)=1gi(ζ)(λ+γi(ˉN)+μqi(ξ))(γi(ˉN)Hi,λ(ζ,a)+gi(ζ)a⋆). $ |
Let $ \Lambda = {\lambda\in\mathbb{C}, |1\in \sigma(U_\lambda)} $. For $ \lambda\in\mathbb{C}\backslash \Lambda $, the operators $ U_{i, \lambda} $ and $ V_{i, \lambda} $ are compact operators from $ X $ to $ L^1(0, a^\star) $, implying that $ \phi_1(a) $ and $ \phi_3(a) $ are represented by compact operators, and similarly, $ \phi_2(a) $ and $ \phi_4(a) $ are also represented by compact operators. As a result, the operator $ \mathrm{C} $ has a compact resolvent, which confirms that its spectrum $ \sigma(\mathrm{C}) $ constitutes only isolated eigenvalues, i.e., $ \sigma(\mathrm{C}) = \sigma_P(\mathrm{C}) $ (see Theorem 6.29 on page 187 in [48]). Hence, $ \mathbb{C}\backslash\Lambda\subset\rho(\mathrm{C}) $, where $ \rho(\mathrm{C}) $ is the resolvent of operator $ \mathrm{C} $. Therefore, $ \sigma_P(\mathrm{C}) = \sigma(\mathrm{C})\subset\Lambda $. Since $ U_\lambda $ is a compact operator, we have $ \sigma(U_\lambda)\backslash{0} = \sigma_P({U_\lambda)\backslash{0}} $. If $ \lambda\in\Lambda $, there exists an eigenfunction $ \psi_\lambda $ such that $ U_\lambda\psi_\lambda = \psi_\lambda $. It is easy to see that $ (\phi_1(a), \phi_2(a), \phi_3(a), \phi_4(a))^\mathrm{T} $ provides an eigenvector of $ \mathrm{C} $ for the eigenvalue $ \lambda $. Thus, we have $ \Lambda\subset\sigma_P(\mathrm{C}) $, and we can conclude that (4.22) is satisfied.
Lemma 4.2. Consider the operator $ {\mathrm{C}} $ which generates $ C_0 $-semigroup for $ t\geq0 $. Then, $ {\bf T}(t) $ is eventually norm continuous (ENC), and we have
$ ω0(C)=s(C)=supReλ,|λ∈σ(C), $ | (4.25) |
where $ s(\mathrm{C}) $ represents the spectral bound of the operator $ \mathrm{C} $, and $ \omega_0(\mathrm{C}) $ denotes the growth bound of the semigroup $ {\bf T}(t) $.
Proof. To begin, we express the bounded operator $ \mathrm{C} $ as
$ Cϕ=(kh(a)γh(ˉN)00αh(a,ˉx1)−γh(ˉN)−μqh(a)0000kc(a)γc(ˉN)00αc(a,ˉx1)−γc(ˉN)−μqc(a))(ϕ1(a)ϕ2(a)ϕ3(a)ϕ4(a)), $ |
for $ \phi\in X $, $ k_h(a) = -\frac{\partial g_h(a)}{\partial a}-\tau_h(a)-\alpha_h(a, \bar{x}_1)-\mu_{p_h}(a) $ and $ k_c(a) = -\frac{\partial g_c(a)}{\partial a}-\tau_c(a)-\alpha_c(a, \bar{x}_1)-\mu_{p_c}(a) $. To establish the compactness of $ \mathrm{C} $, our strategy is to demonstrate that for any bounded sequence $ (\phi^n){n\in \mathbb{N}} $ in $ X $, the sequence $ (\mathrm{C}\phi^n){n\in \mathbb{N}} $ contains a subsequence that converges uniformly. To accomplish this, we invoke the Arzelà-Ascoli Theorem, which requires us to verify that $ (\mathrm{C}\phi^n){n\in \mathbb{N}} $ is uniformly bounded and uniformly equicontinuous. To prove boundedness, since we assume $ (\phi^n){n\in \mathbb{N}} $ to be bounded, we get
$ \|\mathrm{C}\phi^n\|_1\leq\|\mathrm{C}\|\, \|\phi^n\|_1\leq \|\mathrm{C}\|\, \sup\limits_{n\in \mathbb{N}}\|\phi^n\|_1, $ |
which determines that $ (\mathrm{C}\phi^n){n\in \mathbb{N}} $ is also bounded. For uniform equicontinuity, we consider
$ ∫R|(Cϕ)(a+h)−(Cϕ)(a)|da=∫R|C(a+h)−C(a)||ϕ(a)|da≤∫R|(kh(a+h)γh(ˉN)00αh(a+h,ˉx1)−γh(ˉN)−μqh(a+h)0000kc(a+h)γc(ˉN)00αc(a+h,ˉx1)−γc(ˉN)−μqc(a+h))−(kh(a)γh(ˉN)00αh(a,ˉx1)−γh(ˉN)−μqh(a)0000kc(a)γc(ˉN)00αc(a,ˉx1)−γc(ˉN)−μqc(a))||ϕ1(a)ϕ2(a)ϕ3(a)ϕ4(a)|da=∫R|kh(a+h)−kh(a)000αh(a+h,ˉx1)−αh(a,ˉx1)−μqh(ah)−μqh(a)0000kc(a+h)−kc(a)000αc(a+h,ˉx1)−αc(a,ˉx1)−μqc(a+h)−μqc(a)||ϕ1(a)ϕ2(a)ϕ3(a)ϕ4(a)|da≤‖ϕ‖∫R|kh(a+h)−kh(a)000αh(a+h,ˉx1)−αh(a,ˉx1)−μqh(ah)−μqh(a)0000kc(a+h)−kc(a)000αc(a+h,ˉx1)−αc(a,ˉx1)−μqc(a+h)−μqc(a)|da. $ |
Hence, we have shown that $ (\mathrm{C}\phi^n){n\in \mathbb{N}} $ is equicontinuous, and by the Arzelà-Ascoli Theorem, we can conclude that there exists a uniformly convergent subsequence of $ (\mathrm{C}\phi^n){n\in \mathbb{N}} $. Hence, $ \mathrm{C} $ is compact, which in turn implies that $ {\bf T} $ is an ENC semigroup. Since the spectral mapping theorem can be applied to ENC semigroups, we have the spectral determined growth condition given by $ \omega_0(\mathrm{C}) = s(\mathrm{C}) $. Thus, we obtain (4.25).
The local exponential asymptotic stability of the steady-state solution $ \omega = 0 $ of (4.15) is established when $ \omega_0(\mathrm{C}) < 0 $. Specifically, there exist constants $ \epsilon > 0 $, $ M\geq1 $, and $ \gamma < 0 $ such that if $ x\in X $ and $ |x|\leq\epsilon $, then the solution $ \omega(t, x) $ of (4.15) exists globally and satisfies $ |\omega(t, x)|\leq M \exp(\gamma t)|x| $ for all $ t > 0 $. In order to examine the stability of steady states, it is necessary to identify the dominant singular point within the set $ \Lambda $, which corresponds to the element with the highest real value. By utilizing (4.22) and (4.25), we can then determine the growth bound of the semigroup $ {\bf T} $.
Lemma 4.3. For any $ \lambda\in \mathbb{R} $, the operator $ U_{i, \lambda} $ is nonsupporting with respect to $ X_+ $ and
$ limλ→+∞r(Ui,λ)=0 $ | (4.26) |
holds.
Proof. By Eqs. (4.23) and (4.24), we can conclude that the operator $ U_{i, \lambda}, \lambda\in \mathbb{R} $ is strictly positive. To prove that $ U_\lambda, \lambda\in \mathbb{R} $ is non-supporting, we can easily demonstrate the inequality
$ Ui,λψ≥⟨fi,λ,ψ⟩c,c=1∈X+,ψ∈X+, $ | (4.27) |
where the linear function $ f_{i, \lambda} $, is
$ ⟨fi,λ,ψ⟩=∫a⋆0[si(ζ)gi(ζ)(λ+γi(ˉN)+μqi(a))exp(−∫a0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)g(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)exp(∫ζ0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)gi(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)]ψ(ζ)dζ. $ | (4.28) |
Thereby, it leads us to $ U_{i, \lambda}^{n+1}\psi\geq\langle f_{i, \lambda}, \psi\rangle\langle f_{i, \lambda}, c\rangle^n c, \, \forall n. $ which holds for all $ \psi\in X_+ $, where $ f_{i, \lambda} $ is strictly positive and the constant function $ c = 1 $ is a quasi-interior point of $ L^1(0, a^\star) $. This implies that $ \langle F, U_{i, \lambda}^n\rangle > 0 $ for every pair $ \psi\in X_+\backslash{0} $, $ F\in X_+^{*}\backslash{0} $, and therefore $ U_{i, \lambda}, , \lambda\in \mathbb{R} $ is non-supporting. We then use inequality (4.27) and take the duality pairing with the eigenfunctional $ F_{i, \lambda} $ of $ U_{i, \lambda} $ corresponding to $ r(U_{i, \lambda}) $, yielding
$ r(U_{i, \lambda})\langle F_\lambda, \psi\rangle\geq\langle F_\lambda, e\rangle\langle f_{i, \lambda}, \psi\rangle. $ |
Assuming $ \psi = c $, we obtain the inequality:
$ r(U_{i, \lambda})\geq\langle f_{i, \lambda}, c\rangle, $ |
where
$ ⟨fi,λ,c⟩=∫a⋆0si(ζ)gi(ζ)(λ+γi(ˉN)+μqi(ζ))exp(−∫a0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)g(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)exp(∫ζ0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)gi(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)dζ. $ | (4.29) |
It follows that
$ ⟨fλ,c⟩≥ϵ∫a⋆01gi(ζ)(λ+γi(ˉN)+μqi(ζ))exp(−∫a0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)g(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)exp(∫ζ0τi(ξ)+αi(ˉx1,ξ)+λ+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)gi(ξ)(λ+γi(ˉN)+μqi(ξ))dξ)dζ. $ | (4.30) |
By using the positivity of $ \gamma_i(\bar{N}), \, \mu_{p_i}, \, \mu_{q_i}, \, \alpha_i $ and $ \tau_i $, we conclude the following:
$ \lim\limits_{\lambda\to+\infty}r(U_{i, \lambda}) = 0. $ |
Hence proved. The previous lemma implies that the function $ \lambda\to r(U_{i, \lambda}) $ is decreasing for all $ \lambda\in \mathbb{R} $. Moreover, if there exists a $ \lambda\in \mathbb{R} $ such that $ r(U_{i, \lambda}) = 1 $, then it follows that $ \lambda\in \Lambda $, as $ r(U_{i, \lambda})\in\sigma_P(U_{i, \lambda}) $. Combining this with the monotonicity property of $ r(U_{i, \lambda}) $ and inequality (4.26), we obtain the following result.
Lemma 4.4. There exists a unique $ \lambda_0\in\mathbb{R}\cap\Lambda $ such that $ r(U_{i, \lambda_{0}}) = 1 $, and $ \lambda_0 > 0 $ if $ r(U_{0}) > 1 $; $ \lambda_0 = 0 $ if $ r(U_{0}) = 1 $; $ \lambda_0 < 0 $ if $ r(U_{0}) < 1 $.
We will demonstrate that $ \lambda_0 $ is a dominant singular point, utilizing Theorem 6.13 in [49].
Lemma 4.5. If there exists a $ \lambda\in\Lambda, \lambda\neq\lambda_0 $, then $ \mathit{\text{Re}}\lambda < \lambda_0 $.
Proof. Let $ \lambda\in \Lambda $ and $ U_{i, \lambda}\psi = \psi $. Then $ |\psi|(a) = |\psi(a)| $, and we have $ |U_{i, \lambda}\psi| = |\psi| $. Therefore, we obtain $ U_{i, \text{Re}\lambda}\psi\geq\psi $. By taking the duality pairing with $ F_{\text{Re}\lambda}\in X^\star_+ $, we get $ r(U_{i, \text{Re}\lambda})\langle F_{\text{Re}\lambda}, |\psi|\rangle\geq\langle F_{\text{Re}\lambda}, |\psi|\rangle $. We have $ r(U_{i, \text{Re}\lambda})\geq 1 $, as $ F_{\text{Re}\lambda} $ is strictly positive. Since $ r(U_{i, \lambda}), \lambda\in\mathbb{R} $ is a declining function, we conclude that $ \text{Re}\lambda\leq\lambda_0 $. Suppose $ \text{Re}\lambda = \lambda_0 $. Then $ U_{i, \lambda_0}|\psi| = |\psi| $. If we assume $ U_{i, \lambda_0}|\psi| > |\psi| $, then taking the duality pairing with the eigenfunctional $ F_0 $ corresponding to $ r(U_{i, \lambda_0}) = 1 $ results in $ \langle F_{0}, |\psi|\rangle > \langle F_{0}, |\psi|\rangle $, which is a contradiction. Therefore, we have $ U_{i, \lambda_0}|\psi| = |\psi| $, and we can deduce that $ |\psi| = c\psi_{0} $, where constant $ c $ is assumed to be $ 1 $, and $ \psi_0 $ is the eigenfunction relating to $ r(U_{i, \lambda_0}) = 1 $. Hence, we have $ \psi(a) = \psi_{0}(a)\exp(iv(a)) $ for a real-valued function $ v(a) $. Substituting this into $ U_{i, \lambda_0}\psi_0 = |U_{i, \lambda}\psi| $ yields
$ αi(a,ˉx1)gi(a)(λ0+γi(ˉN)+μqi(a))∫a⋆0exp(∫ζaτi(ξ)+αi(ˉx1,ξ)+λ0+μpi(ξ)−γi(ˉN)αi(ˉx1,ξ)λ0+γi(ˉN)+μqi(ξ)dξ)ψ0(ζ)dζ=|αi(a,ˉx1)gi(a)(λ0+iImλ+γi(ˉN)+μqi(a))∫a⋆0exp(∫ζaτi(ξ)+αi(ˉx1,ξ)+λ0+iImλ+μpi(ξ)−γ(iˉN)αi(ˉx1,ξ)λ0+iImλ+γi(ˉN)+μqi(ξ)dξ)exp(iv(ζ))ψ0(ζ)dζ|. $ |
Lemma 6.12 [49] implies that $ \text{Im}\lambda+v(\zeta) $ equals a constant $ \Theta $. Using the fact that $ U_{i, \lambda}\psi = \psi $, we obtain the equation $ \exp(i\Theta)U_{i, \lambda_0}\psi_{\lambda_0} = \psi_{\lambda_0}\exp(iv(\zeta)) $. This equation shows that if $ \Theta = v(\zeta) $, then $ \text{Im}\lambda = 0 $. Therefore, the proof is complete.
Theorem 4.1. The equilibrium state $ (\bar{p}_h(a), \bar{q}_h(a), \bar{p}_c(a), \bar{q}_c(a))^{\mathrm{T}}, $ for (2.1)–(2.4), is locally asymptotically stable if $ r(U_{0}) < 1 $ and locally unstable if $ r(U_{0}) > 1 $.
Proof. Lemmas 4.4 and 4.5 suggests that $ \sup{\text{Re}\lambda:1\in\sigma_P(U_{i, \lambda})} = \lambda_0 $. This implies that if $ r(U_{0}) < 1 $, then $ s(\mathrm{C}) = \sup{\text{Re}\lambda:1\in\sigma_P(U_{i, \lambda})} < 0 $. Conversely, if $ r(U_{0}) > 1 $, then $ s(\mathrm{C}) = \sup{\text{Re}\lambda:1\in\sigma_P(U_{i, \lambda})} > 0 $. Therefore, the proof is complete.
In this section, we show simulations of the model to investigate the evolution of healthy and mutated sub-populations of quiescent and proliferative cells in relation to the cell cycle dynamics. These simulations are performed in the Matlab and we have used finite volume method with discretization central upwind scheme. Table 2 displays the model parameters that were used in the simulations. A maximum cell age of $ 50 $ is assumed, and the time step and spatial step size are set to $ \Delta t = 0.02 $ and $ \Delta a = 0.5 $, respectively. Additionally, we use unit speed, such that $ g_h(a) = g_c(a) = 1 $.
Para. | Description | Healthy | Mutated | Unit |
$ m $ | Mutation rate | 0.2 | - | day$ ^{-1} $ |
$ \nu_i $ | Maximum transition rate from quiescent to proliferation phase | 0.6 [50] | 0.6 | day$ ^{-1} $ |
$ \theta_i $ | Total cell population beyond which $ \Gamma $ is zero | $ 0.095\times{10}^{6} $ [50] | $ 0.095\times{10}^{6} $ | - |
$ \kappa_i $ | Hill coefficient | 1 [50] | 1 | - |
$ \rho_{1, i} $ | Maximal effect of Cyclin $ \mathsf{D-CDK 4/6} $ complex on the division of cell | 0.7 | 0.7 | - |
$ \rho_{2, i} $ | Value of Cyclin $ \mathsf{D-CDK 4/6} $ complex to achieve half maximum effect | 0.35 | 0.35 | - |
$ \gamma_{1, i} $ | Hill coefficient | 8 | 8 | - |
$ \sigma_{1, i} $ | Maximum rate of switching cells from proliferating to quiescent phase | 0.01 | 0.01 | - |
$ \sigma_{2, i} $ | Switching Cyclin $ \mathsf{D-CDK 4/6} $ complex value, after that $ \alpha $ is close to zero | 0.5[45] | 0.45 | |
$ \sigma_{3, i} $ | Switching age value beyond which $ \alpha $ is close to zero | 14 | 15 | h |
$ \gamma_{2, i} $ | Hill coefficient | 7 | 7 | - |
$ \gamma_{3, i} $ | Hill coefficient | 7 | 7 | - |
$ k_t $ | Rate constant which measures the effect of total population on growth factors | $ 1.80\times{10}^{-9} $ [51] | $ 1.80\times{10}^{-9} $ | - |
Steady-state dynamics of healthy and mutated cell populations:
This case study aims to study the steady-state dynamics of both healthy and mutated cell populations, where the death rates are set to $ \mu_{p_h} = \mu_{q_h} = 0.0014 $ and $ \mu_{p_c} = \mu_{q_c} = 0.0014 $. Additionally, $ \gamma(N) = 6.8964\times10^{-6} $ and $ \rho_1 = 1.0 $, while the mutation rate is set to $ m = 0.1 $. The initial conditions for the healthy cell populations are defined as normal distributions $ p_h(a, 0) = q_h(a, 0) = \frac{k_0}{\sqrt{2\pi\sigma^2}}\exp\big(-\frac{(a-\mu)^2}{2\sigma^2}\big) $, where $ k_0 = 10^6 $, $ \mu = 2 $, and $ \sigma^2 = 200 $. Similar initial distributions are used for the mutated cell populations, but with $ k_0 = 10^2 $. A normal distribution is used because it provides a good approximation of the cell distribution within a population by incorporating heterogeneity with respect to cell age in a population.
Figure 3 illustrates the number density distribution of different cell populations, namely healthy (a) proliferating, (b) quiescent, and mutated (c) proliferating, (d) quiescent, as they evolve and eventually reach a steady-state. Here, the cell age is measured in time and cell density is measured as cells per cubic millimeter. The mutation rate is set to $ m = 0.1 $. Over time, the population of mutated proliferating and quiescent cells gradually declines as they grow faster and occupy the tissues completely. However, the total population of cells (consisting of both healthy and mutated proliferating and quiescent cells) denoted by N(t) increases rapidly, as shown in Figure 4(a), before eventually reaching a steady-state. Furthermore, Figure 4(b) demonstrates that the growth factors, which are influenced by the cell population N(t), initially increase due to the low cell count and subsequently decrease until reaching an equilibrium. Finally, Figure 4(c) portrays the transition rate of cells, denoted by $ \gamma(N) $, from quiescent phase to proliferating phase.
Exponential growth of mutated cell populations:
For this case study, we selected a mutation rate of $ m = 0.2 $ and used death rates of $ \mu_{p_h} = \mu_{q_h} = 0.0014 $ and $ \mu_{p_c} = \mu_{q_c} = 0.0010 $. Additionally, we altered several other parameters, including $ \nu_{1, c} = 0.045 $, $ \rho_{1, c} = 1.0 $, $ \rho_{2, c} = 30 $, $ \sigma_{1, c} = 0.040 $ and $ \sigma_{2, c} = 0.45 $. Figure 5 shows the cell density distribution of healthy and mutated proliferating and quiescent cells, respectively. Both healthy subpopulations exhibit the trends of reaching a trivial steady-state with time. However, the mutated cell populations $ p_c(a, t) $ and $ q_c(a, t) $ increase exponentially, mimicking cancerous behavior. Figure 6 illustrate the total number of cells, growth factors, and the switching function from quiescent to proliferating phase, $ \gamma_i(N) $, respectively. The total cell population, which comprises healthy and mutated proliferating and quiescent cell populations, exhibits an exponential increase in cell number over time. The growth factors are initially maximum due to the low cell count and gradually decline until reaching extremely low levels. Finally, the transition functions $ \gamma_h $ and $ \gamma_c $ also decline as the cell population grows.
The proposed model has some limitations that should be considered. Firstly, exclusion of cell heterogeneity, that is a essential aspect to account for cellular noise. Additionally, the feedback model that is comprised of growth factors is fairly simple straight forward, and to characterize the activation of the Cyclin $ \mathsf{D} $ complex, all signaling pathways should be considered. Moreover, at the microscale, it would have been necessary to model cell cycle dynamics separately for healthy and mutated cell populations to investigate their respective compartments' more heterogeneous behavior. Moreover, while the Cyclin-CDK$ \mathsf{CDK 4/6} $ are crucial for the $ G_1 $ to $ S $ phase transition, the model overlooks the other restriction point that detects DNA damage in the $ S $ phase.
This reserch presents non-linear, multiscale modeling of physiologically-structured healthy and mutated quiescent and proliferating cells in relation to cell cycle dynamics. We modeled reversible transitioning from quiescent to proliferating cells and vice versa. We checked the wellposedness of the model, derive non-trivial equilibrium solutions and find spectral criteria for local stability. We also performed numerical simulations to study the impact of Cyclin $ \mathsf{D-CDK 4/6} $ complex on the transition between two sub-populations. Furthermore, we predict that the Cyclin' complex plays an important role in the reversible transition, and any abnormality in this transition can result in cancer.
All authors declare no conflicts of interest in this paper.
[1] | Castagna AM, Pangon A, Choi T, et al. (2012) The role of soft segment molecular weight on microphase separation and dynamics of bulk polymerized polyureas.Macromol 45: 8438-8444. |
[2] | Grujicic M, Yavari R, Snipes JS, et al. (2013) Discrete Element Modeling and Analysis of Structural Collapse/Survivability of a Building Subjected to Improvised Explosive Device (IED) Attack.Adv Mater Sci Appl 2: 9-24. |
[3] | Grujicic M, Pandurangan B, He T, et al. (2010) Computational Investigation of Impact Energy Absorption Capability of Polyurea Coatings via Deformation-Induced Glass Transition.Mater Sci Eng A, 527: 7741-7751. |
[4] | Grujicic M, Arakere A, Pandurangan B, et al. (2012) Computational Investigation of Shock-Mitigation Efficacy of Polyurea when used in a Combat Helmet: A Core Sample Analysis.Multidisc Model Mater Struc 8: 297-331. |
[5] | Grujicic M, Snipes JS, Ramaswami S, et al.Meso-Scale Computational Investigation of Shock-Wave Attenuation by Trailing Release-Wave in Different Grades of Polyurea.J Mater Eng Perform, in press . |
[6] | Settles G. (2013) Pennsylvania State University.work in progress . |
[7] | Grujicic M, Yavari R, Snipes JS, et al. (2013) Molecular-Level Computational Investigation of Shock-Wave Mitigation Capability of Polyurea.J Mater Sci 47: 8197-8215. |
[8] | Grujicic M, Snipes JS, Ramaswami S, et al. (2013) Coarse-Grained Molecular-Level Analysis of Polyurea Properties and Shock-Mitigation Potential.J Mater Eng Perform 22: 1964-1981. |
[9] | Sun H, Ren P, Fried JR (1998) The Compass Force-field: Parameterization and Validation for Phosphazenes.Comput Theor Polym Sci 8: 229-246. |
[10] | Grujicic M, Bell WC, Pandurangan B, et al. (2010) Blast-wave Impact-Mitigation Capability of Polyurea When Used as Helmet Suspension Pad Material.Mater Des 31: 4050-4065. |
[11] | Grujicic A, LaBerge M, Grujicic M, et al. (2012) Potential Improvements in Shock-Mitigation Efficacy of A Polyurea-Augmented Advanced Combat Helmet: A Computational Investigation.J Mater Eng Perform 21: 1562-1579. |
[12] | Grujicic M, Olson GB, Owen WS (1985) Mobility of the β1-γ1' Martensitic Interface in Cu-Al-Ni. I. Experimental Measurements.Metall Trans 16A: 1723-1734. |
[13] | Grujicic M, Bell WC, Pandurangan B, et al. (2012) Inclusion of Material Nonlinearity and Inelasticity into a Continuum-Level Material Model for Soda-lime Glass.J Mater Des 35: 144-155. |
[14] | Grujicic M, Snipes JS, Ramaswami S, et al. (2013) Molecular- and Domain-Level Microstructure Dependent Material Model for Nano-Segregated Polyurea.Multidisc Model Mater Struc 9: 548-578. |
[15] | Grujicic M, He T, Pandurangan B, et al. (2011) Development and Parameterization of a Time-Invariant (Equilibrium) Material Model for Segmented Elastomeric Polyureas.J Mater: Des Appl 225: 182-194. |
[16] | Grujicic M, d'Entremont BP, Pandurangan B, et al. (2012) A Study of the Blast-induced Brain White-Matter Damage and the Associated Diffuse Axonal Injury.Multidisc Model Mater Struc 8: 213-245. |
[17] | Grujicic M, d'Entremont BP, Pandurangan B, et al. (2012) Concept-Level Analysis and Design of Polyurea for Enhanced Blast-Mitigation Performance.J Mater Eng Perform 21: 2024-2037. |
[18] | Grujicic M, Pandurangan B (2012) Meso-Scale Analysis of Segmental Dynamics in Micro-phase Segregated Polyurea.J Mater Sci 47: 3876-3889. |
1. | W. K. Sexton, M. Fidero, J. C. Spain, L. Jiang, K. Bucalo, J. M. Cruse-Sanders, G. S. Pullman, Characterization of endophytic bacterial communities within greenhouse and field-grown rhizomes of three rare pitcher plant species (Sarracenia oreophila, S. leucophylla, and S. purpurea spp. venosa) with an emphasis on nitrogen-fixing bacteria, 2020, 447, 0032-079X, 257, 10.1007/s11104-019-04372-8 | |
2. | Esther Menendez, Lorena Carro, 2019, Chapter 15, 978-3-030-18932-7, 333, 10.1007/978-3-030-18933-4_15 | |
3. | Jocelyn T. Zarate, Genevieve Mae B. Aquino, Joana Marie C. Cruz, Neilyn O. Villa, Albert Remus R. Rosana, J. Cameron Thrash, Draft Genome Sequence of Enterobacter sp. Strain AD2-3, Isolated from a Postmining Site in Benguet, Philippines, 2019, 8, 2576-098X, 10.1128/MRA.00563-19 | |
4. | Lorena Carro, Imen Nouioui, Vartul Sangal, Jan P. Meier-Kolthoff, Martha E. Trujillo, Maria del Carmen Montero-Calasanz, Nevzat Sahin, Darren Lee Smith, Kristi E. Kim, Paul Peluso, Shweta Deshpande, Tanja Woyke, Nicole Shapiro, Nikos C. Kyrpides, Hans-Peter Klenk, Markus Göker, Michael Goodfellow, Genome-based classification of micromonosporae with a focus on their biotechnological and ecological potential, 2018, 8, 2045-2322, 10.1038/s41598-017-17392-0 | |
5. | Lorena Carro, Esther Menéndez, 2020, 9780128184691, 169, 10.1016/B978-0-12-818469-1.00014-6 | |
6. | Imen Nouioui, Carlos Cortés-albayay, Lorena Carro, Jean Franco Castro, Maher Gtari, Faten Ghodhbane-Gtari, Hans-Peter Klenk, Louis S. Tisa, Vartul Sangal, Michael Goodfellow, Genomic Insights Into Plant-Growth-Promoting Potentialities of the Genus Frankia, 2019, 10, 1664-302X, 10.3389/fmicb.2019.01457 | |
7. | Rachael E. Antwis, Ellen L. Fry, Chloë E. James, Natalie Ferry, 2020, Chapter Nine, 9781108654418, 182, 10.1017/9781108654418.009 | |
8. | Natalya Harutyunyan, Almagul Kushugulova, Narine Hovhannisyan, Astghik Pepoyan, One Health Probiotics as Biocontrol Agents: One Health Tomato Probiotics, 2022, 11, 2223-7747, 1334, 10.3390/plants11101334 | |
9. | Imen Nouioui, Sung-min Ha, Inwoo Baek, Jongsik Chun, Michael Goodfellow, Genome insights into the pharmaceutical and plant growth promoting features of the novel species Nocardia alni sp. nov, 2022, 23, 1471-2164, 10.1186/s12864-021-08257-y | |
10. | Federica Massa, Roberto Defez, Carmen Bianco, Exploitation of Plant Growth Promoting Bacteria for Sustainable Agriculture: Hierarchical Approach to Link Laboratory and Field Experiments, 2022, 10, 2076-2607, 865, 10.3390/microorganisms10050865 | |
11. | O.V. Ryabova, A.A. Gagarina, Actinomycetes as the basis of Probiotics for Plants, 2021, 37, 2500-2341, 3, 10.21519/0234-2758-2021-37-2-3-19 | |
12. | Celia Alvarez-Sastre, Lorena Carro, 2022, Desert Actinobacterial Strains Increase Salt Stress Resilience in Crops, 17, 10.3390/environsciproc2022016017 | |
13. | Noureddine Bouras, Amine Yekkour, Slim Smaoui, Lotfi Mellouli, Mustapha Barakate, 2021, Chapter 16, 978-981-16-5616-3, 337, 10.1007/978-981-16-5617-0_16 | |
14. | Maria Tartaglia, Rosaria Sciarrillo, Daniela Zuzolo, Alessia Postiglione, Antonello Prigioniero, Pierpaolo Scarano, Valentino Ruggieri, Carmine Guarino, Exploring an enhanced rhizospheric phenomenon for pluricontaminated soil remediation: Insights from tripartite metatranscriptome analyses, 2022, 428, 03043894, 128246, 10.1016/j.jhazmat.2022.128246 | |
15. | Lorena Carro, Álvaro Peix, Encarna Velázquez, 2021, Chapter 12, 978-3-030-77594-0, 289, 10.1007/978-3-030-77595-7_12 | |
16. | Wei Zhao, Hong Jiang, Xiao-Wan Liu, Jian Zhou, Bin Wu, Polyene Macrolactams from Marine and Terrestrial Sources: Structure, Production Strategies, Biosynthesis and Bioactivities, 2022, 20, 1660-3397, 360, 10.3390/md20060360 | |
17. | Sigita Jurkonienė, Rima Mockevičiu̅tė, Jurga Jankauskienė, Elžbieta Jankovska-Bortkevič, Gabija Armalytė, Virgilija Gavelienė, Application of Commercial Plant Probiotics Improves the Berry Yield and Quality of Field-Grown Blackcurrant, 2021, 1, 2692-1952, 615, 10.1021/acsagscitech.1c00115 | |
18. | Saurabh Gangola, Samiksha Joshi, Divya Joshi, Jyoti Rajwar, Shilippreet Kour, Jyoti Singh, Saurabh Kumar, 2022, Chapter 18, 978-3-030-86168-1, 385, 10.1007/978-3-030-86169-8_18 | |
19. | Martino Schillaci, Aida Raio, Fabiano Sillo, Elisa Zampieri, Shahid Mahmood, Muzammil Anjum, Azeem Khalid, Mauro Centritto, Pseudomonas and Curtobacterium Strains from Olive Rhizosphere Characterized and Evaluated for Plant Growth Promoting Traits, 2022, 11, 2223-7747, 2245, 10.3390/plants11172245 | |
20. | Juan Sanjuán, Maria Caridad Nápoles, Daniel Pérez-Mendoza, María J. Lorite, Dulce N. Rodríguez-Navarro, Microbials for Agriculture: Why Do They Call Them Biostimulants When They Mean Probiotics?, 2023, 11, 2076-2607, 153, 10.3390/microorganisms11010153 | |
21. | O. V. Ryabova, A. A. Gagarina, Actinomycetes as the Basis of Probiotics for Plants, 2022, 58, 0003-6838, 827, 10.1134/S0003683822070055 | |
22. | Tunjung Mahatmanto, Nabilah Rizka Estiningtyas, Syafira Kintan Maharani, Ahmad Sultan Shalahuddin, Muhammad Yuda Dwitama Agustian, Untung Murdiyatmo, Feather Waste Biorefinery using Chryseobacterium sp. A9.9 Adapted to Feathers as its Sole Carbon and Nitrogen Source, 2022, 13, 1877-2641, 4137, 10.1007/s12649-022-01766-1 | |
23. | Nattakorn Kuncharoen, Masahiro Yuki, Takuji Kudo, Moriya Okuma, Auttaporn Booncharoen, Wuttichai Mhuantong, Somboon Tanasupawat, Comparative genomics and proposal of Streptomyces radicis sp. nov., an endophytic actinomycete from roots of plants in Thailand, 2022, 254, 09445013, 126889, 10.1016/j.micres.2021.126889 | |
24. | Pooja Bhadrecha, Shilpy Singh, Vagish Dwibedi, ‘A plant’s major strength in rhizosphere’: the plant growth promoting rhizobacteria, 2023, 205, 0302-8933, 10.1007/s00203-023-03502-2 | |
25. | Viabhav Kumar Upadhayay, Manoj Kumar Chitara, Dhruv Mishra, Manindra Nath Jha, Aman Jaiswal, Geeta Kumari, Saipayan Ghosh, Vivek Kumar Patel, Mayur G. Naitam, Ashish Kumar Singh, Navneet Pareek, Gohar Taj, Damini Maithani, Ankit Kumar, Hemant Dasila, Adita Sharma, Synergistic impact of nanomaterials and plant probiotics in agriculture: A tale of two-way strategy for long-term sustainability, 2023, 14, 1664-302X, 10.3389/fmicb.2023.1133968 | |
26. | Rhydum Sharma, Richa Salwan, Vivek Sharma, 2023, 9780323917346, 271, 10.1016/B978-0-323-91734-6.00016-8 | |
27. | Phumzile Mkhize, Josephine Malatji, Phetole Mangena, 2023, 2, 9789815165319, 37, 10.2174/9789815165319123020007 | |
28. | Valeria Razmilic, Imen Nouioui, Andrey Karlyshev, Rana Jawad, Martha E. Trujillo, Jose M. Igual, Barbara A. Andrews, Juan A. Asenjo, Lorena Carro, Michael Goodfellow, Micromonospora parastrephiae sp. nov. and Micromonospora tarensis sp. nov., isolated from the rhizosphere of a Parastrephia quadrangularis plant growing in the Salar de Tara region of the Central Andes in Chile, 2023, 73, 1466-5026, 10.1099/ijsem.0.006189 | |
29. | Shubhra Singh, Douglas J. H. Shyu, 2024, 978-1-83767-271-4, 481, 10.1039/BK9781837673131-00481 | |
30. | Kavitha Komire, Mamta Tiwari, Prakash Chandra Gupta, Theeshan Bahorun, Nisha Sharma, 2024, Chapter 9, 978-3-031-68137-0, 197, 10.1007/978-3-031-68138-7_9 | |
31. | Wuyu Liu, Guoqing Wang, Shiming Wen, Yiwen Zhao, Yuxin Ding, Baihui Yao, Zhelin Wang, Duntao Shu, Gehong Wei, Juan Chen, Zhouping Shangguan, Microbiome-Mediated Mechanisms Regulating Adaptability to Iron Deficiency in the Intercropping System of Soybean and Maize, 2025, 15, 2073-4395, 286, 10.3390/agronomy15020286 | |
32. | Yifei Huang, Yuankun Zhou, Hengyi Xu, Unlocking Wellness: Probiotics as Key Drivers in Functional Food Innovation and Health Promotion, 2025, 15, 2076-3417, 6498, 10.3390/app15126498 |
Description | State |
Cyclin $ \mathsf{D-CDK 4/6} $ | $ x_1 $ |
$ \mathsf{E2F} $ | $ x_2 $ |
$ \mathsf{Rb} $ | $ x_3 $ |
$ \mathsf{p21} $ | $ x_4 $ |
Para. | Description | Healthy | Mutated | Unit |
$ m $ | Mutation rate | 0.2 | - | day$ ^{-1} $ |
$ \nu_i $ | Maximum transition rate from quiescent to proliferation phase | 0.6 [50] | 0.6 | day$ ^{-1} $ |
$ \theta_i $ | Total cell population beyond which $ \Gamma $ is zero | $ 0.095\times{10}^{6} $ [50] | $ 0.095\times{10}^{6} $ | - |
$ \kappa_i $ | Hill coefficient | 1 [50] | 1 | - |
$ \rho_{1, i} $ | Maximal effect of Cyclin $ \mathsf{D-CDK 4/6} $ complex on the division of cell | 0.7 | 0.7 | - |
$ \rho_{2, i} $ | Value of Cyclin $ \mathsf{D-CDK 4/6} $ complex to achieve half maximum effect | 0.35 | 0.35 | - |
$ \gamma_{1, i} $ | Hill coefficient | 8 | 8 | - |
$ \sigma_{1, i} $ | Maximum rate of switching cells from proliferating to quiescent phase | 0.01 | 0.01 | - |
$ \sigma_{2, i} $ | Switching Cyclin $ \mathsf{D-CDK 4/6} $ complex value, after that $ \alpha $ is close to zero | 0.5[45] | 0.45 | |
$ \sigma_{3, i} $ | Switching age value beyond which $ \alpha $ is close to zero | 14 | 15 | h |
$ \gamma_{2, i} $ | Hill coefficient | 7 | 7 | - |
$ \gamma_{3, i} $ | Hill coefficient | 7 | 7 | - |
$ k_t $ | Rate constant which measures the effect of total population on growth factors | $ 1.80\times{10}^{-9} $ [51] | $ 1.80\times{10}^{-9} $ | - |
Description | State |
Cyclin $ \mathsf{D-CDK 4/6} $ | $ x_1 $ |
$ \mathsf{E2F} $ | $ x_2 $ |
$ \mathsf{Rb} $ | $ x_3 $ |
$ \mathsf{p21} $ | $ x_4 $ |
Para. | Description | Healthy | Mutated | Unit |
$ m $ | Mutation rate | 0.2 | - | day$ ^{-1} $ |
$ \nu_i $ | Maximum transition rate from quiescent to proliferation phase | 0.6 [50] | 0.6 | day$ ^{-1} $ |
$ \theta_i $ | Total cell population beyond which $ \Gamma $ is zero | $ 0.095\times{10}^{6} $ [50] | $ 0.095\times{10}^{6} $ | - |
$ \kappa_i $ | Hill coefficient | 1 [50] | 1 | - |
$ \rho_{1, i} $ | Maximal effect of Cyclin $ \mathsf{D-CDK 4/6} $ complex on the division of cell | 0.7 | 0.7 | - |
$ \rho_{2, i} $ | Value of Cyclin $ \mathsf{D-CDK 4/6} $ complex to achieve half maximum effect | 0.35 | 0.35 | - |
$ \gamma_{1, i} $ | Hill coefficient | 8 | 8 | - |
$ \sigma_{1, i} $ | Maximum rate of switching cells from proliferating to quiescent phase | 0.01 | 0.01 | - |
$ \sigma_{2, i} $ | Switching Cyclin $ \mathsf{D-CDK 4/6} $ complex value, after that $ \alpha $ is close to zero | 0.5[45] | 0.45 | |
$ \sigma_{3, i} $ | Switching age value beyond which $ \alpha $ is close to zero | 14 | 15 | h |
$ \gamma_{2, i} $ | Hill coefficient | 7 | 7 | - |
$ \gamma_{3, i} $ | Hill coefficient | 7 | 7 | - |
$ k_t $ | Rate constant which measures the effect of total population on growth factors | $ 1.80\times{10}^{-9} $ [51] | $ 1.80\times{10}^{-9} $ | - |