Immunotherapy is a targeted therapy that can be applied to cervical cancer patients to prevent DNA damage caused by human papillomavirus (HPV). The HPV infects normal cervical cells withing a specific cell age interval, i.e., between the G1 to S phase of the cell cycle. In this study, we developed a new mathematical model of age-dependent immunotherapy for cervical cancer. The model is a four-dimensional first-order partial differential equation with time- and age-independent variables. The cell population is divided into four sub-populations, i.e., susceptible cells, cells infected by HPV, precancerous cells, and cancer cells. The immunotherapy term has been added to precancerous cells since these cells can experience regression if appointed by proper treatments. The immunotherapy process is closely related to the rate of T-cell division. The treatment works in the same cell cycle that stimulates and inhibits the immune system. In our model, immunotherapy is represented as a periodic function with a small amplitude. It is based on the fluctuating interaction between T-cells and precancerous cells. We have found that there are two types of steady-state conditions, i.e., infection-free and endemic. The local and global stability of an infection-free steady-state has been analyzed based on basic reproduction numbers. We have solved the Riccati differential equation to show the existence of an endemic steady-state. The stability analysis of the endemic steady-state has been determined by using the perturbation approach and solving integral equations. Some numerical simulations are also presented in this paper to illustrate the behavior of the solutions.
Citation: Eminugroho Ratna Sari, Lina Aryati, Fajar Adi-Kusumo. An age-structured SIPC model of cervical cancer with immunotherapy[J]. AIMS Mathematics, 2024, 9(6): 14075-14105. doi: 10.3934/math.2024685
[1] | Qiuyan Wang, Weixin Liu, Jianming Wang, Yang Yan . A class of nearly optimal codebooks and their applications in strongly regular Cayley graphs. AIMS Mathematics, 2024, 9(7): 18236-18246. doi: 10.3934/math.2024890 |
[2] | Nuttawoot Nupo, Sayan Panma . Certain structural properties for Cayley regularity graphs of semigroups and their theoretical applications. AIMS Mathematics, 2023, 8(7): 16228-16239. doi: 10.3934/math.2023830 |
[3] | Yang Yan, Xingguo Zhang, Rize Jin, Limin Zhou . A construction of strongly regular Cayley graphs and their applications to codebooks. AIMS Mathematics, 2024, 9(2): 2672-2683. doi: 10.3934/math.2024132 |
[4] | Jiangmin Pan, Yingnan Zhang . On cubic semisymmetric bi-Cayley graphs on nonabelian simple groups. AIMS Mathematics, 2022, 7(7): 12689-12701. doi: 10.3934/math.2022702 |
[5] | Tariq A. Alraqad, Hicham Saber . On the structure of finite groups associated to regular non-centralizer graphs. AIMS Mathematics, 2023, 8(12): 30981-30991. doi: 10.3934/math.20231585 |
[6] | He Song, Longmin Wang, Kainan Xiang, Qingpei Zang . The continuity of biased random walk's spectral radius on free product graphs. AIMS Mathematics, 2024, 9(7): 19529-19545. doi: 10.3934/math.2024952 |
[7] | Yan Wang, Kai Yuan, Ying Zhao . Perfect directed codes in Cayley digraphs. AIMS Mathematics, 2024, 9(9): 23878-23889. doi: 10.3934/math.20241160 |
[8] | Igal Sason . Observations on graph invariants with the Lovász ϑ-function. AIMS Mathematics, 2024, 9(6): 15385-15468. doi: 10.3934/math.2024747 |
[9] | Piyapat Dangpat, Teerapong Suksumran . Regularity of extended conjugate graphs of finite groups. AIMS Mathematics, 2022, 7(4): 5480-5498. doi: 10.3934/math.2022304 |
[10] | Bilal A. Rather, Fawad Ali, Nasim Ullah, Al-Sharef Mohammad, Anwarud Din, Sehra . Aα matrix of commuting graphs of non-abelian groups. AIMS Mathematics, 2022, 7(8): 15436-15452. doi: 10.3934/math.2022845 |
Immunotherapy is a targeted therapy that can be applied to cervical cancer patients to prevent DNA damage caused by human papillomavirus (HPV). The HPV infects normal cervical cells withing a specific cell age interval, i.e., between the G1 to S phase of the cell cycle. In this study, we developed a new mathematical model of age-dependent immunotherapy for cervical cancer. The model is a four-dimensional first-order partial differential equation with time- and age-independent variables. The cell population is divided into four sub-populations, i.e., susceptible cells, cells infected by HPV, precancerous cells, and cancer cells. The immunotherapy term has been added to precancerous cells since these cells can experience regression if appointed by proper treatments. The immunotherapy process is closely related to the rate of T-cell division. The treatment works in the same cell cycle that stimulates and inhibits the immune system. In our model, immunotherapy is represented as a periodic function with a small amplitude. It is based on the fluctuating interaction between T-cells and precancerous cells. We have found that there are two types of steady-state conditions, i.e., infection-free and endemic. The local and global stability of an infection-free steady-state has been analyzed based on basic reproduction numbers. We have solved the Riccati differential equation to show the existence of an endemic steady-state. The stability analysis of the endemic steady-state has been determined by using the perturbation approach and solving integral equations. Some numerical simulations are also presented in this paper to illustrate the behavior of the solutions.
The directed strongly regular graph (or DSRG for short) is a potential generalization of the well-established strongly regular graphs. It was introduced by Duval in [4] in 1988. Although it has received less attention compared to its undirected counterparts, DSRGs have gained popularity in recent years.
A DSRG can be interpreted in the framework of adjacency matrices. The adjacency matrix of a directed graph X of order n is denoted by A=A(X)=(aij)n×n. We employ the notation I=In and J=Jn to represent the identity matrix and all-one matrix of size n, respectively. A directed graph X is a DSRG with parameters (n,k,μ,λ,t) if and only if it satisfies:
(ⅰ) JA=AJ=kJ,
(ⅱ) A2=tI+λA+μ(J−I−A).
When t=k, X becomes undirected strongly regular graph. Duval [4] demonstrated that DSRGs with t=0 correspond to the doubly regular tournaments. Consequently, it is customary to assume that 0<t<k.
From history, many infinite families of DSRGs were constructed in light of several parameters of DSRGs as well as some sporadic examples. Despite the extensive literature on the existence, structure, and construction of DSRGs for various parameter values [8,10], a significant number of DSRGs remain shrouded in mystery, with their existence yet to be determined. In fact, it is a challenging problem for the complete characterization of DSRGs. By using representation theory as a powerful tool, He and Zhang [6] obtained a large family of DSRGs on dihedral groups, which extended certain findings presented in [10]. S. Hayat, J. H. Koolen and M. Riaz gave a similar conclusion for undirected strongly regular graphs in [5]. For more results, one may refer to [1,9].
Our objective is to derive novel infinite families of DSRGs. Inspired by the methodology outlined in [6], we explore directed strongly regular Cayley graphs derived from Dic4n, where n is odd. As previously noted in[3], the dicyclic group exhibits a distinct nature compared to the dihedral group. Consequently, it would be intriguing to investigate the potential applications of this group.
This paper is structured as follows. Initially, we introduce several classes of directed strongly regular Cayley graphs (or DSRCGs for short) over dicyclic groups. Subsequently, we provide a criterion for the Cayley graph C(Dic4n,M∪Mb∪Mb2∪Mb3) with M∩M(−1)=∅ to be directed strongly regular.
For comprehensive insights into representation theory and associated concepts, we follow [7].
Let G be a finite group. We denote by IRR(G) (resp. Irr(G)) the set of all non-equivalent irreducible representations (resp. irreducible characters) of G. Our subsequent discussion needs the characters associated with cyclic groups.
Lemma 2.1 ([7]). For a cyclic group Cn=⟨ν⟩ of order n, IRR(Cn)={ℓs|0≤s≤n−1}, with ℓs(νk)=ωskn (0≤s,k≤n−1), where ωn=e2πin represents the n-th primitive root of unity.
We denote the group algebra of a group G over the complex field C as CG. It contains all formal sums of the form ∑g∈Gmgg, where mg∈C. The multiplication is defined as:
(∑g∈Gmgg)(∑h∈Gnhh)=∑g∈G∑h∈Gmgnhgh. |
Lemma 2.2. [7] For an element A=∑g∈Gmgg∈CG with G being abelian group, we have
mg=1|G|∑χ∈Irr(G)χ(A)¯χ(g),∀g∈G. |
We now introduce some notations regarding multisets. Let A represent a multiset characterized by a multiplicity function, denoted as δA:S→N. Here, δA(x) denotes the frequency of occurrence of x in A. We define x is an element of A (i.e., x∈A) if and only if the multiplicity of x, as determined by δA(x), is greater than zero. For two multisets A and B, with multiplicity functions δA and δB, respectively, their union, represented as A⊎B, is determined by the function δA⊎B=δA+δB. For instance, consider A={1,1,2,3}, and B={1,2,3,3,4,}, then the union A⊎B={1,1,1,2,2,3,3,3,4}.
Let ¯X be an element in CG corresponding to any multisubset X of the group G. ¯X may be expressed as
¯X=∑x∈XδX(x)x. |
Let S be a non-empty subset of a finite group G. We denote the set {s−1∣s∈S} as S(−1). If S∩S(−1)=∅, then it is called antisymmetric. Let us assume that the identity element e∉S. In this case, the graph Γ=C(G,S) is a directed Cayley graph, its vertex set is V(Γ)=G, and there is an arc from x to y, represented by x→y, if yx−1∈S.
Utilizing the group algebra, the lemma presented below establishes a criterion for a Cayley graph to be DSRG.
Lemma 2.3. [4] C(G,S) is a DSRG with parameters (n,k,μ,λ,t) if and only if |G|=n, |S|=k, and
¯S2=te+λ¯S+μ(¯G−e−¯S). |
Our primary focus is on a non-abelian group–the dicyclic group. It is denoted as Dic4n, and is typically presented with the following group presentation:
Dic4n=⟨x,y|x2n=1,xn=y2,y−1xy=x−1⟩. |
For n odd, let x2=α and y=β. Then, we have
Dic4n=⟨α,β|αn=β4=1,β−1αβ=α−1⟩={αk,αkβ,αkβ2,αkβ3|0≤k≤n−1}. | (2.1) |
The following relationships will be commonly referenced within the context of our discussion.
Lemma 2.4. Given that n is odd, for the dicyclic group Dic4n, we can observe the following properties:
(ⅰ) αkβ=βα−k; αkβ2=β2αk; αkβ3=β3α−k;
(ⅱ) (αkβ)−1=αkβ3; (αkβ2)−1=α−kβ2.
Proof: Using (2.1), the conclusions are immediate.
We aim to present various constructions of DSRCGs derived over Dic4n, with n being an odd integer.
Any subset S of Dic4n can be expressed as S=S0∪S1b∪S2b2∪S3b3 with S0,S1,S2,S3⊆Cn. Our first main result is:
Theorem 3.1. Γ=C(Dic4n,S0∪S1b∪S2b2∪S3b3) is a DSRG with parameters (4n,|S0|+|S1|+|S2|+|S3|),μ,λ,t) if and only if the following four statements hold:
(ⅰ) ¯S02+¯S1¯S(−1)3+¯S22+¯S3¯S(−1)1=(t−μ)e+(λ−μ)¯S0+μ¯Cn;
(ⅱ) ¯S0¯S1+¯S1¯S(−1)0+¯S2¯S3+¯S3¯S(−1)2=(λ−μ)¯S1+μ¯Cn;
(ⅲ) 2¯S0¯S2+¯S1¯S(−1)1+¯S3¯S(−1)3=(λ−μ)¯S2+μ¯Cn;
(ⅳ) ¯S0¯S3+¯S1¯S(−1)2+¯S2¯S1+¯S3¯S(−1)0=(λ−μ)¯S3+μ¯Cn.
Proof: According to Lemma 2.4, the following holds:
(¯S0∪¯S1b∪¯S2b2∪¯S3b3)2=¯S02+¯S0¯S1b+¯S0¯S2b2+¯S0¯S3b3+¯S1¯S(−1)0b+¯S1¯S(−1)1b2+¯S1¯S(−1)2b3+¯S1¯S(−1)3+¯S2¯S0b2+¯S2¯S1b3+¯S22+¯S2¯S3b+¯S3¯S(−1)0b3+¯S3¯S(−1)1+¯S3¯S(−1)2b+¯S3¯S(−1)3b2=¯S02+¯S1¯S(−1)3+¯S22+¯S3¯S(−1)1+(¯S0¯S1+¯S1¯S(−1)0+¯S2¯S3+¯S3¯S(−1)2)b+(2¯S0¯S2+¯S1¯S(−1)1+¯S3¯S(−1)3)b2+(¯S0¯S3+¯S1¯S(−1)2+¯S2¯S1+¯S3¯S(−1)0)b3. |
From Lemma 2.3, the Cayley graph Γ=C(Dic4n,S0∪S1b∪S2b2∪S3b3) is recognized as a DSRG with parameters (4n,|S0|+|S1|+|S2|+|S3|),μ,λ,t) if and only if
(¯S0∪¯S1b∪¯S2b2∪¯S3b3)2=te+λ(¯S0∪¯S1b∪¯S2b2∪¯S3b3)+μ(¯Cn+¯Cnb+¯Cnb2+¯Cnb3)−μe−μ(¯S0∪¯S1b∪¯S2b2∪¯S3b3)=((t−μ)e+(λ−μ)¯S0+μ¯Cn)+((λ−μ)¯S1+μ¯Cn)b+((λ−μ)¯S2+μ¯Cn)b2+((λ−μ)¯S3+μ¯Cn)b3. |
From the above two equations, we complete the proof.
Let S0=S2=M and S1=S3=N in Theorem 3.1, then we have:
Theorem 3.2. Γ=C(Dic4n,M∪Nb∪Mb2∪Nb3) is a DSRG with parameters (4n,2(|M|+|N|),μ,λ,t) if and only if the following two conditions hold for t=μ and M, N:
(ⅰ) 2(¯M2+¯N¯N(−1))=(λ−μ)¯M+μ¯Cn;
(ⅱ) 2¯N(¯M+¯M(−1))=(λ−μ)¯N+μ¯Cn.
Proof: As S0=S2=M and S1=S3=N, from (i), (iii) of Theorem 3.1, we derive
2(¯M2+¯N¯N(−1))=(t−μ)e+(λ−μ)¯M+μ¯Cn, |
and
2(¯M2+¯N¯N(−1))=(λ−μ)¯M+μ¯Cn. |
Therefore, comparing the above two equations, we have t=μ and 2(¯M2+¯N¯N(−1))=(λ−μ)¯M+μ¯Cn.
Similarly, as S0=S2=M and S1=S3=N, from (ⅱ), (ⅳ) of Theorem 3.1, we have 2¯N(¯M+¯M(−1))=(λ−μ)¯N+μ¯Cn.
This completes the proof.
Setting N=M in Theorem 3.2, we have:
Theorem 3.3. C(Dic4n,M∪Mb∪Mb2∪Mb3) is a DSRG with parameters (4n,4|M|,μ,λ,t) if and only if the following two conditions hold for t,μ and M:
(ⅰ) t=μ;
(ⅱ) 2¯M(¯M+¯M(−1))=(λ−μ)¯M+μ¯Cn.
Remark 3.1. Theorem 3.3 (ii) also implied that
2¯M(−1)(¯M+¯M(−1))=(λ−μ)¯M(−1)+μ¯Cn. |
Thus, by Theorem 3.3 (ⅱ) and Remark 3.1, we derive
2(¯M+¯M(−1))2=(λ−μ)(¯M+¯M(−1))+2μ¯Cn. |
Next, we can get several classes of DSRGs from the above results.
Corollary 3.1. For an odd number n, suppose that the two conditions hold for M,N⊆Cn:
(ⅰ) ¯M+¯M(−1)=¯Cn−e;
(ⅱ) ¯N¯N(−1)−¯M¯M(−1)=ε¯Cn,ε=0 or 1.
Then, Γ=C(Dic4n,M∪Nb∪Mb2∪Nb3) is a DSRG with parameters (4n,2n−2+2ε,n−1+2ε,n−3+2ε,n−1+2ε).
Furthermore, if M satisfies (i), and N=Mh or N=M(−1)h for h∈Cn, then Γ is a DSRG with parameters (4n,2n−2,n−1,n−3,n−1).
Proof: By (i), we have |M|=n−12. By (ii), we have |N|2=|M|2+εn=(|M|+ε)2, because ε=0 or 1, and |M|=n−12. Thus, we obtain |N|=|M|+ε=n−12+ε. By (ⅰ) and (ⅱ), we obtain:
2(¯M2+¯N¯N(−1))=2(¯M2+¯M¯M(−1)+ε¯Cn)=2¯M(¯M+¯M(−1))+2ε¯Cn=2¯M(¯Cn−e)+2ε¯Cn=−2¯M+2(|M|+ε)¯Cn=−2¯M+(n−1+2ε)¯Cn, |
and
2¯N(¯M+¯M(−1))=2¯N(¯Cn−e)=2|N|Cn−2¯N=−2¯N+(n−1+2ε)¯Cn. |
Therefore, the desired result is obtained through Theorem 3.2.
In the following, we denote l=nv where v|n. Given that ⟨av⟩ is a subgroup of Cn, a corresponding transversal within Cn consists of the elements {e,a,a2,⋯,av−1}. Let T⊆{e,a,a2,⋯,av−1} as a subset. We define:
T⟨av⟩=⋃at∈Tat⟨av⟩, |
where at⟨av⟩ are coset of ⟨a⟩ in Cn, for at∈T. Then we have:
Corollary 3.2. Let T⊆{e,a,a2,⋯,av−1} as a subset, with v∣n, and the following two conditions hold for M,N⊆Cn:
(ⅰ) N=T⟨av⟩=M∪⟨av⟩;
(ⅱ) N⨄N(−1)=Cn⨄⟨av⟩.
Then, Γ=C(Dic4n,M∪Nb∪Mb2∪Nb3) is a DSRG with parameters (4n,2n,n+l,n−l,n+l).
Proof: According to (i) and (ii), we derive |N|=|M|+l=n+l2, and ¯M+¯M(−1)=¯Cn−¯⟨av⟩. Therefore,
2(¯M2+¯N¯N(−1))=2(¯M2+(¯M+¯⟨av⟩)(¯M(−1)+¯⟨av⟩))=2(¯M2+¯M¯M(−1)+(¯M+¯M(−1))¯⟨av⟩+l¯⟨av⟩)=2(¯M+¯⟨av⟩)(¯M+¯M(−1))+2l¯⟨av⟩=2(¯M+¯⟨av⟩)(¯Cn−¯⟨av⟩)+2l¯⟨av⟩=(n+l)¯Cn−2¯M¯⟨av⟩−2¯⟨av⟩¯⟨av⟩+2l¯⟨av⟩=(n+l)¯Cn−2l¯M−2l¯⟨av⟩+2l¯⟨av⟩=−2l¯M+(n+l)¯Cn, |
and
2¯N(¯M+¯M(−1))=2¯N(¯Cn−¯⟨av⟩)=2|N|Cn−2¯N¯⟨av⟩=−2l¯N+(n+l)¯Cn. |
Therefore, the result is obtained through Theorem 3.2.
Corollary 3.3. Let T⊆{e,a,a2,⋯,av−1} be a subset, with v∣n. Assume that the following conditions hold for M⊆Cn:
(ⅰ) M=T⟨av⟩;
(ⅱ) M∪M(−1)=Cn∖⟨av⟩.
Then, Γ=C(Dic4n,M∪Mb∪Mb2∪Mb3) is a DSRG with parameters (4n,4|M|,n−l,n−3l,n−l).
Proof: According to (i) and (ii), we derive
2¯M(¯M+¯M(−1))=2¯M(¯Cn−¯⟨Mv⟩)=2|M|¯Cn−2l¯M=−2l¯M+(n−l)¯Cn. |
By Theorem 3.3, we obtain the result.
Remark 3.2. If Γ=C(Dic4n,M∪Mb∪Mb2∪Mb3) is a DSRG in Corollary 3.3, then we derive M∩M(−1)=∅.
Now, we give an example to illustrate our results, whose parameters are listed in [2].
Example 3.1. Let Dic36=⟨a,b|a9=b4=1,b−1ab=a−1⟩ be the dicyclic group of order 36 and υ=3, then l=93=3 and ⟨a3⟩={e,a3,a6}. Let T={a2}⊆{e,a,a2}. Then, we have M={a2}⟨a3⟩={a,a4,a7} and M(−1)={a2,a5,a8}. Thus, M∪M(−1)={a,a2,a4,a5,a7,a8}=C9∖⟨a3⟩ satisfies the condition of Corollary 3.3. By direct computation, we have
¯S2=6e+0¯S+6(¯Dic36−e−¯S), |
where ¯S=¯M∪¯Mb∪¯Mb2∪¯Mb3. Therefore, by Lemma 2.3, we have Γ=C(Dic36,M∪Mb∪Mb2∪Mb3) is a DSRG with parameters (36,12,6,0,6).
Suppose that M∩M(−1)=∅, i.e., M is an antisymmetric subset of Cn. We end this paper with a criterion for certain Cayley graph to be a DSRG.
Theorem 3.4. Γ=C(Dic4n,M∪Mb∪Mb2∪Mb3) with M∩M(−1)=∅ is a DSRG with parameters (4n,4|M|,μ,λ,t) if and only if the following conditions hold for a subset T of {a,a2,⋯,av−1}:
(ⅰ) M=T⟨av⟩;
(ⅱ) M∪M(−1)=Cn∖⟨av⟩, where v=2nμ−λ is an odd positive integer.
Proof: Let W=M∪M(−1)⊆Cn. Then, δW(h)=0 or 1 for any h∈W. Hence, ¯W=¯M+¯M(−1). By Corollary 3.3, if Γ=C(Dic4n,M∪Mb∪Mb2∪Mb3) satisfying (i) and (ii) is a DSRG, then M∩M(−1)=∅.
Now, we consider the converse part. Suppose that Γ=C(Dic4n,M∪Mb∪Mb2∪Mb3) is a DSRG with parameters (4n,4|M|,μ,λ,t), where n is odd. By Theorem 3.3 and Remark 3.1, we have
2¯W2=(λ−μ)¯W+2μ¯(Cn). | (3.1) |
Therefore, we have χ(W)∈{0,λ−μ2} for any non-principal characters χ(W) of Cn. Now, we define
W={j|j=1,2,⋯,n−1,χj(¯W)=λ−μ2}. |
By Lemma 2.2, we have that
δW(h)=1n∑χ∈Irr(Cn)χ(¯W)¯χ(h)=λ−μ2n∑j∈W¯χj(h)+2|M|n. | (3.2) |
As e∉W, hence, we obtain δW(e)=0, and therefore,
δW(e)=λ−μ2n|W|+2|M|n=0. |
Then, we have 4|M|=(μ−λ)|W|. Thus, the expression (3.2) becomes
δW(h)=μ−λ2n(|W|−∑j∈W¯χj(h)). | (3.3) |
By Eq (3.3), we have |W|−∑j∈W¯χj(h)∈Q. Note that |W|−∑j∈W¯χj(h)∈Z[ωn]. Therefore,
2nμ−λ∈Z. |
By Eq (3.3), we also have
δW(h)=0⇔|W|−∑j∈W¯χj(h)=0⇔χj(h)=1⇔g∈⋂j∈WKχj, |
where j∈W). Let Rdef=⋂j∈WKχj is some subgroup of Cn. Thus, |R|||Cn|. Since |Cn|=n is odd, |R| is odd too. Thus, we have ¯W=¯Cn−¯R. Moreover,
2¯W2=2(¯Cn−¯R)2=2(n−2|R|)¯Cn+2|R|¯R=−2|R|¯W+2(n−|R|)¯Cn, |
then, we have |R|=μ−λ2, n−|R|=μ, and |M|=n−|R|2=μ2. Therfore, R=⟨a2nμ−λ⟩=⟨av⟩. Since |Cn|=n and |R|=μ−λ2 are all odds, we have μ=2nμ−λ is odd too. Thus, we proved (ii). In this case, by Theorem 3.3, we have
(λ−μ)¯M+μ¯Cn=2¯M¯W=2¯M(¯Cn−¯⟨av⟩)=μ¯Cn−2¯M¯⟨av⟩, |
i.e., μ−λ2¯M=¯M¯⟨av⟩. Thus, we have M=T⟨av⟩ for some subset T of {a,a2,⋯,av−1}; therefore, we proved (ⅰ).
Tao Cheng: conceptualization, writing review and editing, data curation, writing original draft preparation, supervision; Junchao Mao: conceptualization, writing review and editing, investigation, project administration.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors would like to express their gratitude to the referee for his or her careful reading and valuable suggestions. This work was supported by National Natural Science Foundation of China (Grant No. 12271311).
The authors declare no conflict of interest.
[1] |
P. Zibako, N. Tsikai, S. Manyame, T. G. Ginindza, Cervical cancer management in Zimbabwe (2019–2020), Plos one, 17 (2022), e0274884. https://doi.org/10.1371/journal.pone.0274884 doi: 10.1371/journal.pone.0274884
![]() |
[2] |
D. S. Chen, I. Mellman, Oncology meets immunology: the cancer-immunity cycle, Immunity, 39 (2013), 1–10. https://doi.org/10.1016/j.immuni.2013.07.012 doi: 10.1016/j.immuni.2013.07.012
![]() |
[3] |
L. Ferrall, K. Y. Lin, R. B. S Roden, C. F. Hung, T. C. Wu, Cervical cancer immunotherapy: facts and hopes, Clin. Cancer Res., 27 (2021), 4953–4973. https://doi.org/10.1158/1078-0432.CCR-20-2833 doi: 10.1158/1078-0432.CCR-20-2833
![]() |
[4] |
T. S. N. Asih, S. Lenhart, S. Wise, L. Aryati, F. Adi-Kusumo, M. S. Hardianti, et al., The dynamics of HPV infection and cervical cancer cells, Bull. Math. Biol., 78 (2016), 4–20. https://doi.org/10.1007/s11538-015-0124-2 doi: 10.1007/s11538-015-0124-2
![]() |
[5] |
L. Aryati, T. S. Noor-Asih, F. Adi-Kusumo, M. S. Hardianti, Global stability of the disease free equilibrium in a cervical cancer model: a chance to recover, FJMS, 103 (2018), 1535–1546. https://doi.org/10.17654/MS103101535 doi: 10.17654/MS103101535
![]() |
[6] |
K. Allali, Stability analysis and optimal control of HPV infection model with early-stage cervical cancer, Biosystems, 199 (2021), 104321. https://doi.org/10.1016/j.biosystems.2020.104321 doi: 10.1016/j.biosystems.2020.104321
![]() |
[7] |
V. V. Akimenko, F. Adi-Kusumo, Stability analysis of an age-structured model of cervical cancer cells and HPV dynamics, Math. Biosci. Eng., 18 (2021), 6155–6177. https://doi.org/10.3934/mbe.2021308 doi: 10.3934/mbe.2021308
![]() |
[8] |
V. V. Akimenko, F. Adi-Kusumo, Age-structured delayed SIPCV epidemic model of HPV and cervical cancer cells dynamics I. Numerical method, Biomath, 10 (2021), 1–23. https://doi.org/10.11145/j.biomath.2021.10.027 doi: 10.11145/j.biomath.2021.10.027
![]() |
[9] |
E. R. Sari, F. Adi-Kusumo, L. Aryati, Mathematical analysis of a SIPC age-structured model of cervical cancer, Math. Biosci. Eng., 19 (2022), 6013–6039. https://doi.org/10.3934/mbe.2022281 doi: 10.3934/mbe.2022281
![]() |
[10] |
V. V. Akimenko, F. Adi-Kusumo, Age-structured delayed SIPCV epidemic model of HPV and cervical cancer cells dynamics II. Convergence of numerical solution, Biomath, 11 (2022), 1–20. https://doi.org/10.55630/j.biomath.2022.03.278 doi: 10.55630/j.biomath.2022.03.278
![]() |
[11] |
N. Cifuentes-Muñoz, R. E. Dutch, R. Cattaneo, Direct cell-to-cell transmission of respiratory viruses: The fast lanes, PLoS pathog., 14 (2018), e1007015. https://doi.org/10.1371/journal.ppat.1007015 doi: 10.1371/journal.ppat.1007015
![]() |
[12] |
D. Kirschner, J. C. Panetta, Modeling immunotherapy of the tumor-immune interaction, J. Math Biol., 37 (1998), 235–252. https://doi.org/10.1007/s002850050127 doi: 10.1007/s002850050127
![]() |
[13] |
E. Nikolopoulou, L. R. Johnson, D. Harris, J. D. Nagy, E. C. Stites, Y. Kuang, Tumour-immune dynamics with an immune checkpoint inhibitor, Lett. Biomath., 5 (2018), 137–159. https://doi.org/10.30707/LiB5.2Nikolopoulou doi: 10.30707/LiB5.2Nikolopoulou
![]() |
[14] |
F. Adi-Kusumo, R. S. Winanda, Bifurcation analysis of the cervical cancer cells, effector cells, and IL-2 compounds interaction model with immunotherapy, FJMS, 99 (2016), 869–883. http://dx.doi.org/10.17654/MS099060869 doi: 10.17654/MS099060869
![]() |
[15] |
A. Pavan, I. Attili, G. Pasello, V. Guarneri, P. F. Conte, L. Bonanno, Immunotherapy in small-cell lung cancer: from molecular promises to clinical challenges, J. Immunother. Cancer, 7 (2019), 205. https://doi.org/10.1186/s40425-019-0690-1 doi: 10.1186/s40425-019-0690-1
![]() |
[16] |
H. Tabbal, A. Septier, M. Mathieu, C. Drelon, S. Rodriguez, C. Djari, et al., EZH2 cooperates with E2F1 to stimulate expression of genes involved in adrenocortical carcinoma aggressiveness, Br. J. Cancer, 121 (2019), 384–394. https://doi.org/10.1038/s41416-019-0538-y doi: 10.1038/s41416-019-0538-y
![]() |
[17] |
M. Barberis, T. Helikar, P. Verbruggen, Simulation of stimulation: Cytokine dosage and cell cycle crosstalk driving timing-dependent T cell differentiation, Front. Physiol., 9 (2018), 879. https://doi.org/10.3389/fphys.2018.00879 doi: 10.3389/fphys.2018.00879
![]() |
[18] |
Y. Pala, M. O. Ertas, An analytical method for solving general Riccati equation, Int. J. Math. Comput. Sci., 11 (2017), 125–130. https://doi.org/10.5281/zenodo.1340124 doi: 10.5281/zenodo.1340124
![]() |
[19] |
Z. Ahmed, M. Kalim, A new transformation technique to find the analytical solution of general second order linear ordinary differential equation, Int. J. Adv. Appl. Sci., 5 (2018), 109–114. https://doi.org/10.21833/ijaas.2018.04.014 doi: 10.21833/ijaas.2018.04.014
![]() |
[20] |
D. P. Singh, A. Ujlayan, An alternative approach to write the general solution of a class of second-order linear differential equations, Resonance, 26 (2021), 705–714. https://doi.org/10.1007/s12045-021-1170-8 doi: 10.1007/s12045-021-1170-8
![]() |
[21] |
A. Wilmer Ⅲ, G. B. Costa, Solving second-order differential equations with variable coefficients, Int. J. Math. Educ. Sci., 39 (2008), 238–243. https://doi.org/10.1080/00207390701464709 doi: 10.1080/00207390701464709
![]() |
[22] |
A. Khan, G. Zaman, Global analysis of an age-structured SEIR endemic model, Chaos Soliton. Fract., 108 (2018), 154–165. https://doi.org/10.1016/j.chaos.2018.01.037 doi: 10.1016/j.chaos.2018.01.037
![]() |
[23] |
H. Inaba, Mathematical analysis of an age-structured SIR epidemic model with vertical transmission, Discrete Cont. Dyn. B, 6 (2006), 69–96. https://doi.org/10.3934/dcdsb.2006.6.69 doi: 10.3934/dcdsb.2006.6.69
![]() |
[24] | A. H. Siddiqi, Functional analysis and applications, Singapore: Springer, 2018. https://doi.org/10.1007/978-981-10-3725-2 |
[25] | A. Pazy, Semigroups of linear operators and applications to partial differential equations, New York: Springer, 1983. https://doi.org/10.1007/978-1-4612-5561-1 |
[26] | G. F. Webb, Theory of nonlinear age-dependent population dynamics, Marcel Dekker, 1985. |
[27] | L. C. Evans, Partial differential equations, American Mathematical Society, 2010. |
[28] |
O. Diekmann, J. A. P. Heesterbeek, J. A. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365–382. https://doi.org/10.1007/BF00178324 doi: 10.1007/BF00178324
![]() |
[29] |
C. Castillo-Chavez, H. W. Hethcote, V. Andreasen, S. A. Levin, W. M. Liu, Epidemiological models with age structure, proportionate mixing, and cross-immunity, J. Math. Biology, 27 (1989), 233–258. https://doi.org/10.1007/BF00275810 doi: 10.1007/BF00275810
![]() |
[30] |
A. K. Miller, K. Munger, F. R. Adler, A mathematical model of cell cycle dysregulation due to Human Papillomavirus infection, Bull. Math. Biol., 79 (2017), 1564–1585. https://doi.org/10.1007/s11538-017-0299-9 doi: 10.1007/s11538-017-0299-9
![]() |
[31] |
H. Cho, Z. Wang, D. Levy, Study of dose-dependent combination immunotherapy using engineered T cells and IL-2 in cervical cancer, J. Theor. Biol., 505 (2020), 110403. https://doi.org/10.1016/j.jtbi.2020.110403 doi: 10.1016/j.jtbi.2020.110403
![]() |
[32] |
S. Patil, R. S. Rao, N. Amrutha, D. S. Sanketh, Analysis of human papillomavirus in oral squamous cell carcinoma using p16: An immunohistochemical study, J. Int. Soc. Prev. Commu., 4 (2014), 61–66. https://doi.org/10.4103/2231-0762.131269 doi: 10.4103/2231-0762.131269
![]() |
[33] | X. Z. Li, J. Yang, M. Martcheva, Age structured epidemic modeling, Switzerland: Springer Nature, 2020. https://doi.org/10.1007/978-3-030-42496-1 |
[34] |
K. Hattaf, Y. Yang, Global dynamics of an age-structured viral infection model with general incidence function and absorption, Int. J. Biomath., 11 (2018), 1850065. https://doi.org/10.1142/S1793524518500651 doi: 10.1142/S1793524518500651
![]() |
[35] |
K. Hattaf, A new class of generalized fractal and fractal-fractional derivatives with non-singular kernels, Fractal Fract., 7 (2023), 395. https://doi.org/10.3390/fractalfract7050395 doi: 10.3390/fractalfract7050395
![]() |
[36] |
K. Hattaf, A new mixed fractional derivative with applications in computational biology, Computation, 12 (2024), 7. https://doi.org/10.3390/computation12010007 doi: 10.3390/computation12010007
![]() |