Loading [MathJax]/jax/output/SVG/jax.js
Research article

Mathematical analysis of nonlinear combination drug delivery

  • Received: 19 January 2025 Revised: 13 March 2025 Accepted: 20 March 2025 Published: 28 March 2025
  • This study examined a single-compartment pharmacokinetic model with intravenous and oral administration. It investigated the trends in steady-state drug exposure and average steady-state plasma drug concentration under multiple dosing regimens. First, the combined drug administration model was a non-autonomous system, and we approximated the solution of the model by estimating the exact upper and lower bounds in conjunction with the comparison theorem for differential equations. We also demonstrated the existence, uniqueness, and stability of the solution. Second, we derived the steady-state drug exposure and compared it with the solution of the intravenous drug delivery model alone. The results indicate that the combined drug delivery scheme offers superior performance. Finally, we theoretically proved the change rule of average steady-state blood drug concentration under different dosing regimens, and verified its feasibility and rationality by combining the numerical simulation results of Phenytoin sodium.

    Citation: Jiali Shi, Zhongyi Xiang. Mathematical analysis of nonlinear combination drug delivery[J]. Electronic Research Archive, 2025, 33(3): 1812-1835. doi: 10.3934/era.2025082

    Related Papers:

    [1] Martin Burger, Marco Di Francesco . Large time behavior of nonlocal aggregation models with nonlinear diffusion. Networks and Heterogeneous Media, 2008, 3(4): 749-785. doi: 10.3934/nhm.2008.3.749
    [2] Iryna Pankratova, Andrey Piatnitski . Homogenization of convection-diffusion equation in infinite cylinder. Networks and Heterogeneous Media, 2011, 6(1): 111-126. doi: 10.3934/nhm.2011.6.111
    [3] Alexander Mielke, Sina Reichelt, Marita Thomas . Two-scale homogenization of nonlinear reaction-diffusion systems with slow diffusion. Networks and Heterogeneous Media, 2014, 9(2): 353-382. doi: 10.3934/nhm.2014.9.353
    [4] Patrick Henning . Convergence of MsFEM approximations for elliptic, non-periodic homogenization problems. Networks and Heterogeneous Media, 2012, 7(3): 503-524. doi: 10.3934/nhm.2012.7.503
    [5] Feiyang Peng, Yanbin Tang . Inverse problem of determining diffusion matrix between different structures for time fractional diffusion equation. Networks and Heterogeneous Media, 2024, 19(1): 291-304. doi: 10.3934/nhm.2024013
    [6] Elisabeth Logak, Isabelle Passat . An epidemic model with nonlocal diffusion on networks. Networks and Heterogeneous Media, 2016, 11(4): 693-719. doi: 10.3934/nhm.2016014
    [7] Patrick Henning, Mario Ohlberger . The heterogeneous multiscale finite element method for advection-diffusion problems with rapidly oscillating coefficients and large expected drift. Networks and Heterogeneous Media, 2010, 5(4): 711-744. doi: 10.3934/nhm.2010.5.711
    [8] José Antonio Carrillo, Yingping Peng, Aneta Wróblewska-Kamińska . Relative entropy method for the relaxation limit of hydrodynamic models. Networks and Heterogeneous Media, 2020, 15(3): 369-387. doi: 10.3934/nhm.2020023
    [9] Martin Heida, Benedikt Jahnel, Anh Duc Vu . Regularized homogenization on irregularly perforated domains. Networks and Heterogeneous Media, 2025, 20(1): 165-212. doi: 10.3934/nhm.2025010
    [10] Mogtaba Mohammed, Mamadou Sango . Homogenization of nonlinear hyperbolic stochastic partial differential equations with nonlinear damping and forcing. Networks and Heterogeneous Media, 2019, 14(2): 341-369. doi: 10.3934/nhm.2019014
  • This study examined a single-compartment pharmacokinetic model with intravenous and oral administration. It investigated the trends in steady-state drug exposure and average steady-state plasma drug concentration under multiple dosing regimens. First, the combined drug administration model was a non-autonomous system, and we approximated the solution of the model by estimating the exact upper and lower bounds in conjunction with the comparison theorem for differential equations. We also demonstrated the existence, uniqueness, and stability of the solution. Second, we derived the steady-state drug exposure and compared it with the solution of the intravenous drug delivery model alone. The results indicate that the combined drug delivery scheme offers superior performance. Finally, we theoretically proved the change rule of average steady-state blood drug concentration under different dosing regimens, and verified its feasibility and rationality by combining the numerical simulation results of Phenytoin sodium.



    NICM: nonischaemic cardiomyopathy; PBMCs: peripheral blood mononuclear cells; DEGs: differentially expressed genes; GO: gene ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; PPI: protein-protein interaction; HF: heart failure; IL-6:interleukin-6; IFN-γ: interferon-gamma; GEO: Gene Expression Omnibus; BP: biological process; MF: molecular function; CC: cellular component; NCBI: National Center for Biotechnology Information; CT: cycle threshold; IBD: inflammatory bowel disease

    Heart failure (HF), a chronic cardiac syndrome, could give rise to a mean survival of 5 years once it is diagnosed, placing millions of individuals globally at risk of death at present, despite the vital progress that has been made in its diagnosis and treatment [1,2]. Nonischaemic cardiomyopathy (NICM) is one of the most common causes of HF [3]. At present, patients diagnosed with HF due to NICM, especially if conventional therapy shows no effect for them, require further diagnostic and therapeutic procedures involving endomyocardial biopsy. However, endomyocardial biopsy is invasive and thus may easily cause adverse complications [4]. Hence, it is of great necessity to establish a series of non-invasive approaches that could be operated by routine assays in the clinic.

    Increasing evidence has demonstrated that dysregulated levels of circulating pro-inflammatory cytokines in peripheral blood mononuclear cells (PBMCs) from patients with NICM exhibit a strong relationship with the disease severity and/or prognosis [5]. Pressure and volume overload trigger profound alterations in cardiac extracellular matrix components, which eventually induce cardiac remodelling by affecting the removal of pro-inflammatory cytokines [6]. Pro-inflammatory monocytes and T-cell infiltration are also observed when myocardial dysfunction and cardiac remodelling occur [7]. These immune cells could mediate hypertension, cardiac hypertrophy and fibrosis by releasing certain cytokines, such as interleukin-6 (IL-6) and interferon-gamma (IFN-γ) [8,9]. PBMCs are a group of immune cells, mainly including T cells, B cells and monocytes/macrophages, the profiling of which has drawn particular attention in a number of diseases due to the convenience and noninvasiveness of PBMCs assessments in recent years [10,11].

    Bioinformatic analysis and gene expression profiling analysis could help us screen molecular markers between patients and healthy individuals, thereby providing novel insights into diseases at multiple levels ranging from the alterations of copy number at the genome level to gene expression at the transcriptome level and even epigenetic alterations [12]. In this study, we sought to determine the differentially expressed genes (DEGs) in PBMCs between patients with NICM and healthy controls based on the gene expression profiles from GSE9128, which was downloaded from the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/). Meanwhile, we analysed gene ontology (GO) involving biological process (BP), molecular function (MF), cellular component (CC), and the KEGG pathways of these DEGs. Furthermore, we constructed a protein-protein interaction (PPI) network of these DEGs and identified the top 15 hub genes with a high degree of connectivity. Finally, we detected the mRNA levels of the top 5 upregulated and downregulated DEGs in PBMCs from a mouse model induced by isoproterenol.

    Gene expression profiles from the GSE9128 dataset were downloaded from the National Center for Biotechnology Information (NCBI) GEO database, which is free and publicly available [13]. The dataset, based on the GPL96 platform [HG-U133A] Affymetrix Human Genome U133A Array, contains samples of PBMCs from 4 NICM patients and 3 normal individuals. Meanwhile, the Series Matrix File of GSE9128 from the GEO database in PubMed was also downloaded. The peripheral blood samples in this dataset were from Caucasian patients with NICM and healthy controls. Chronic HF was classified according to New York Heart Association (NYHA) functional class Ⅲ/Ⅳ as lasting for more than 6 months; the left ventricular ejection fraction (LVEF) of these patients was less than 36%. All of these NICM patients had no history of angiographically abnormal coronary arteries and myocardial infarction. Gene expression profiling was performed by Affymetrix microarrays.

    We screened the DEGs between NICM and healthy samples using GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r), an interactive analysis tool for the GEO database based on the R language. We modified the DEGs as differentially expressed with a logFC < 1 (upregulated genes) or logFC > 1 (downregulated genes) according to the criteria described in the references ADDIN EN.CITE [1,2]. Meanwhile, a P value < 0.05 was regarded as significantly different, aiming to reduce the false positive rate. Finally, after downloading the relatively raw TXT data, we employed visual hierarchical cluster analysis to display the volcano plot and heat map of these DEGs in PBMCs between the NICM and healthy samples using ImageGP (http://www.ehbio.com/ImageGP/index.php/Home/Index/index.html).

    Gene ontology (GO) analysis is a common framework that can annotate genes and gene products involving functions of cellular components (CC), biological pathways (BP) and molecular function (MF) [14]. The Kyoto Encyclopedia of Genes and Genomes (KEGG) is a collection of databases that help to handle genomes, biological pathways, diseases, chemical substances and drugs [15]. In our study, GO and KEGG analysis was performed by the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.ncifcrf.gov) (version 6.7), which is an online biological function database integrating considerable biological data and analysis tools [16]. P < 0.05 was regarded as the cut-off criterion for significant difference. After the key BP, MF, CC and KEGG pathways of these DEGs were determined using DAVID, a scatterplot was generated by ImageGP according to the results of the GO and KEGG pathway analyses.

    The Search Tool for the Retrieval of Interacting Genes (STRING) is an online tool to assess and integrate protein-protein interaction (PPI) information, containing physical and functional associations [17]. The PPI network can then identify the core hub genes and key gene modules between patients and healthy individuals. To find the potential correlations of the DEGs, we used STRING to map the PPI network of these DEGs. Subsequently, we used the Cytoscape software platform on the basis of the PPI associations to select 15 hub genes according to their ranked order of connectivity degree. Moreover, the Molecular Complex Detection (MCODE) app was also utilized to screen the typical modules of the PPI network in Cytoscape in line with degree cutoff = 2, k-core = 2, node score cutoff = 0.2 and max depth = 100.

    All animal experimental procedures in this study were approved by the Animal Care and Use Committee of Renmin Hospital of Wuhan University and were performed in accordance with the Care and Use of Laboratory Animals published by the US National Institute of Health (Revised 2011) (approval number: WDRX-2017K012). C57/B6 male mice (8-10 weeks old) weighing 25.1 ± 1.6 g were obtained from the Institute of Laboratory Animal Science of the Chinese Academy of Medical Sciences (Beijing, China). The animals were randomly divided into two groups (20 per group) for treatment, the control group and the NICM group. For the mouse model of NICM, isoprenaline (dissolved in sterile saline) was injected (60 mg/kg/day) intraperitoneally once daily for 2 consecutive weeks [18]. Success for HF model creation was assessed by serial echocardiography demonstrating dilation, hypertrophy, and ventricular dysfunction [19].

    Blood samples were collected by intracardiac puncture (3–5 mL per animal). As described previously [20], the PBMCs of the mice were isolated from whole blood by density gradient centrifugation with Ficoll-Paque Plus reagent and suspended in QIAzol Lysis Reagent (Qiagen, Dusseldorf, Germany). Subsequently, the layer with the PBMCs was removed and washed twice with saline (0.9%), and the recovered cells were used for total RNA extraction.

    Total RNA was isolated using a TRIzol (Invitrogen, Carlsbad, CA, USA) assay, the concentrations and purities of which were quantified using an ultraviolet spectrophotometer. The RNA was then reverse transcribed according to the previous description [21]. The expression levels of the top 5 upregulated genes and the top 5 downregulated genes in PBMCs were normalized to GAPDH. Relative mRNA expression levels were analysed by the 2-∆∆ cycle threshold (CT) method. Apart from detecting the mRNA levels of the top 5 upregulated genes and the top 5 downregulated genes in the control group and NICM group, we also measured the mRNA levels of these genes in the isolated different cell types, the T cells (CD4+ and CD8+), B cells (CD20+) and monocytes (CD14+), from the obtained PBMCs of the mice by positive selection (Anti-PE MicroBeads UltraPure, Miltenyi Biotec, Teterow, Germany).

    The obtained data are presented as the mean ± SD (standard deviation) and were assessed by a two-tailed Student's t-test. P < 0 05 was considered statistically significant.

    In the present study, a comparison of 3 healthy samples with 4 NICM samples was executed by employing the GEO2R online analysis tool based on a P value < 0.05 and logFC ≤ −1 or logFC ≥ 1 criteria. After analysing the GSE9128 microarray, a total of 371 DEGs were selected, consisting of 288 upregulated DEGs and 83 downregulated DEGs. The expression proportion of these DEGs is displayed in the volcano plot (Figure 1A). The expression levels of the top 10 upregulated and top 10 downregulated DEGs are displayed in a heat map to visualize the results (Figure 1B). Among these 371 DEGs, the top 5 upregulated DEGs were EGR3, MBNL1, PTGS2, EGR2 and TXNIP, while the top 5 downregulated DEGs were GABARAPL3, NPIPB3, RGS1, CREM and JUN. The biological functions and gene titles of these 10 DEGs are presented in Table 1.

    Figure 1.  Volcano plot and heat map of the differentially expressed genes (DEGs) in peripheral blood mononuclear cells (PBMCs) between normal samples and patients with nonischaemic cardiomyopathy (NICM). (A) Volcano plot of genes detected in PBMCs from patients with NICM. Red indicates upregulated DEGs; blue indicates downregulated DEGs; green indicates no difference. (B) Heat map of the top 10 upregulated DEGs and top 10 downregulated DEGs in PBMCs.
    Table 1.  The top 5 up-regulated and down-regulated differentially expressed genes in PBMCs from patients with NICM.
    DEGs Gene title Gene symbol LogFC Biological function
    Up-regulated early growth response 3 EGR3 3.25296217 Regulating systemic autoimmunity
    muscleblind like splicing regulator 1 MBNL1 3.22538031 Regulating pre-mRNA alternative splicing
    prostaglandin-endoperoxide synthase 2 PTGS2 3.21406038 Catalyzings the conversion of arachidonic acid to prostaglandins; Immune response suppression
    early growth response 2 EGR2 2.96759548 Intrinsic negative regulator of DC immunogenicity; Mmonocyte/macrophage cell fate determination
    thioredoxin interacting protein TXNIP 2.76192581 Linking oxidative stress to inflammasome activation
    Down-regulated GABA type A receptor associated protein like 3 pseudogene GABARAPL3 −1.78817386 Inhibiting inflammasome
    nuclear pore complex interacting protein family member B3 NPIPB3 −1.82835573 Modulating Type I interferon
    regulator of G-protein signaling 1 RGS1 −1.87910047 Regulating chemokine signalling
    cAMP responsive element modulator CREM −1.93814192 Suppresseing IL-2 and TCRζ chain gene transcription
    Jun proto-oncogene, AP-1 transcription factor subunit JUN −2.12183042 Modulating inflammation

     | Show Table
    DownLoad: CSV

    To gain a more extensive and in-depth knowledge of the selected DEGs, we used DAVID to analyse significantly enriched GO functions. By analysing GO enrichment of all upregulated and downregulated DEGs via DAVID, we discovered that the upregulated DEGs in BP were principally enriched in inflammatory response, immune response, response to lipopolysaccharide, cellular response to mechanical stimulus and negative regulation of translation, while the downregulated DEGs were enriched in circadian regulation of gene expression, negative regulation of sequence-specific DNA binding transcription factor activity, innate immune response, cellular response to hydrogen peroxide and negative regulation of I-kappaB kinase/NF-kappaB signaling. For CC, the upregulated DEGs were mainly enriched in extracellular exosome, nucleus, membrane and nucleoplasm and cytosol, while the downregulated DEGs were enriched in nucleus, nucleoplasm, centrosome and transcription factor complex. MF analysis revealed that the upregulated DEGs were mainly responsible for ATP binding, poly(A) RNA binding, chromatin binding, nucleotide binding and GTP binding. By contrast, the downregulated DEGs may function in protein binding, mitochondrial uncoupling and cytokine activity. The details of the GO analysis are presented in Figure 2 AC and Table 2.

    Figure 2.  Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in PBMCs of NICM. (A) The enriched GO terms in the biological process (BP); (B) the enriched GO terms in the cellular component (CC); (C) the enriched GO terms in the molecular function (MF); (D) the enriched KEGG pathways in PBMCs of NICM.
    Table 2.  Gene ontology analysis of differentially expressed genes in PBMCs from patients with NICM.
    Expression Category Term Count PValue
    Up-regulated GOTERM_BP_DIRECT GO:0006954~inflammatory response 9 0.008
    GOTERM_BP_DIRECT GO:0006955~immune response 8 0.012
    GOTERM_BP_DIRECT GO:0032496~response to lipopolysaccharide 7 0.001
    GOTERM_BP_DIRECT GO:0071260~cellular response to mechanical stimulus 6 2.66E-04
    GOTERM_BP_DIRECT GO:0017148~negative regulation of translation 5 0.001
    GOTERM_CC_DIRECT GO:0070062~extracellular exosome 51 1.41E-04
    GOTERM_CC_DIRECT GO:0005634~nucleus 49 0.011
    GOTERM_CC_DIRECT GO:0016020~membrane 27 5.87E-05
    GOTERM_CC_DIRECT GO:0005654~nucleoplasm 27 0.021
    GOTERM_CC_DIRECT GO:0005829~cytosol 19 0.049
    GOTERM_MF_DIRECT GO:0005524~ATP binding 25 0.039
    GOTERM_MF_DIRECT GO:0044822~poly(A) RNA binding 23 0.006
    GOTERM_MF_DIRECT GO:0003682~chromatin binding 11 0.008
    GOTERM_MF_DIRECT GO:0000166~nucleotide binding 10 0.016
    GOTERM_MF_DIRECT GO:0005525~GTP binding 10 0.044
    Down-regulated GOTERM_BP_DIRECT GO:0032922~circadian regulation of gene expression 3 0.008
    GOTERM_BP_DIRECT GO:0043433~negative regulation of sequence-specific DNA binding transcription factor activity 3 0.011
    GOTERM_BP_DIRECT GO:0045087~innate immune response 3 0.001
    GOTERM_BP_DIRECT GO:0070301~cellular response to hydrogen peroxide 2 0.003
    GOTERM_BP_DIRECT GO:0043124~negative regulation of I-kappaB kinase/NF-kappaB signaling 2 ≤0.001
    GOTERM_CC_DIRECT GO:0005634~nucleus 14 0.032
    GOTERM_CC_DIRECT GO:0005654~nucleoplasm 13 7.15E-04
    GOTERM_CC_DIRECT GO:0005813~centrosome 5 0.019
    GOTERM_CC_DIRECT GO:0005667~transcription factor complex 3 0.007
    GOTERM_MF_DIRECT GO:0005515~protein binding 58 0.003
    GOTERM_MF_DIRECT GO:0008083~mitochondrial uncoupling 5 0.031
    GOTERM_MF_DIRECT GO:0005125~cytokine activity 5 0.032
    GO: Gene Ontology.

     | Show Table
    DownLoad: CSV

    To obtain more integrated information regarding the crucial pathways of these selected DEGs, KEGG pathways were also analysed via DAVID, the result of which is presented in Table 3 and Figure 2D. The downregulated DEGs were enriched in protein processing in the endoplasmic reticulum, legionellosis, NF-kappa B signaling pathway, NOD-like receptor signaling pathway, estrogen signaling pathway while the upregulated DEGs were enriched in inflammatory bowel disease (IBD), TGF-beta signaling pathway, neurotrophin signaling pathway and Epstein-Barr virus infection.

    Table 3.  KEGG pathway analysis of differentially expressed PBMCs from patients with NICM.
    Category Term Count % P-Value Genes
    Down-regulated ptr04141:Protein processing in endoplasmic reticulum 11 0.030 5.49E-04 SEC23A, HSP90AB1, HSP90B1, HSP90AA1, UBE2J1, MAN1A1, EDEM3, UBE2D1, SAR1B, TRAM1, SSR1
    ptr05134:Legionellosis 7 0.020 5.25E-04 CXCL2, CXCL8, IL1B, TLR4, HSPD1, CASP1, SAR1B
    ptr04064:NF-kappa B signaling pathway 7 0.020 0.004304358 PTGS2, BCL2A1, CXCL8, IL1B, TLR4, TAB2, ATM
    ptr04621:NOD-like receptor signaling pathway 6 0.017 0.002195688 HSP90AB1, HSP90AA1, CXCL8, IL1B, CASP1, TAB2
    ptr04915:Estrogen signaling pathway 6 0.017 0.023154474 HSP90AB1, HSP90B1, HSP90AA1, FKBP5, CREB1, PRKACB
    Up-regulated cjc05321:Inflammatory bowel disease (IBD) 3 0.030 0.000636646 JUN, SMAD3, RORA
    cjc04350:TGF-beta signaling pathway 3 0.030 0.03910693 SMAD7, ID1, SMAD3
    cjc04722:Neurotrophin signaling pathway 3 0.030 0.000546883 NFKBIE, JUN, MAPKAPK2
    cjc05169:Epstein-Barr virus infection 3 0.030 0.000983148 NFKBIE, JUN, TNFAIP3
    KEGG: Kyoto Encyclopedia of Genes and Genomes; FDR: False Discovery Rate.

     | Show Table
    DownLoad: CSV

    Furthermore, 15 hub genes identified from the Cytoscape software, which occupy the core position in the PPI network, were found in accordance with a high degree of connectivity (Table 4 & Figure 3A). We built the PPI network from the top 15 hub genes via the information from the STRING protein query of public databases. The top 15 hub genes, possessing a high degree of connectivity in NICM, are as follows: JUN, MYC, HSP90AA1, PCNA, CREB1, IL1B, IL8, SMARCA2, TLR4, RB1, RANBP2, EGR1, PTGS2, ENO1 and XPO1. Among the top 15 hub genes, only one gene (JUN) was significantly downregulated in the PBMCs from the patients with NICM, while the others were significantly increased in the PBMCs. Intriguingly, JUN is not only the DEG with the highest connectivity but is also the downregulated DEG with the highest fold change. Additionally, the top 15 hub genes could interact with 149 DEGs directly. These hub genes could also interact with each other and exhibited a strong correlation. For example, JUN could directly interact with IL8, RB1, EGR1, PTGS2, IL1B, SMARCA2, CREB1, HSP90AA1, TLR4 and MYC, while CREB1 interacts with IL1B, PTGS2, JUN, SMARCA2, RB1, MYC, XPO1, TLR4 and HSP90AA1. Moreover, we used the MCODE plug-in to detect the highest modules in the PPI network. We identified the top 2 modules of the network; Module 1 (Figure 3B) was associated with immune response and protein processing in the endoplasmic reticulum, while Module 2 (Figure 3C) was associated with the signaling pathway and cellular response to inflammation. Collectively, the PPI network suggested that the top 15 hub genes, especially JUN and CREB, may play essential roles in the immune response of NICM.

    Table 4.  Top 15 hub genes with higher degree of connectivity.
    Gene Degree of connectivity P value Log FC
    JUN 47 ≤0.001 -2.121
    MYC 44 ≤0.001 1.932
    HSP90AA1 40 0.040 1.155
    PCNA 35 0.004 1.649
    CREB1 33 0.006 1.617
    IL1B 33 0.009 1.351
    IL8 33 ≤0.001 1.241
    SMARCA2 30 0.002 1.131
    TLR4 30 0.008 1.821
    RB1 28 0.003 1.176
    RANBP2 28 0.002 1.855
    EGR1 24 0.003 1.721
    PTGS2 24 0.001 3.214
    ENO1 23 0.013 1.953
    XPO1 22 0.015 1.067
    Top 15 hub genes with higher degree of connectivity

     | Show Table
    DownLoad: CSV
    Figure 3.  Protein-protein interaction (PPI) network of the top 15 hub genes with high connectivity degrees and module analysis. (A) PPI network of the top 15 hub genes. (B) Module 1. (B) Module 2.

    To ensure the credibility of the GSE9128 microarray and obtain further credible analysis, we re-identified the top 5 upregulated genes and the top 5 downregulated genes via qPCR in a mouse model induced by isoprenaline. The results from PCR (Figure 4) demonstrated that the mRNA expression levels of EGR3, MBNL1, PTGS2, EGR2 and TXNIP in PBMCs from mice with NICM were significantly higher than those in control mice (P < 0.05), while the mRNA expression levels of GABARAPL3, NPIPB3, RGS1, CREM and JUN in the NICM group were statistically lower than those in the control group (P < 0.05).

    Figure 4.  Validation of the top 5 upregulated and top 5 downregulated DEGs in a mouse model induced by isoprenaline.*P < 0.05, **P < 0.001, ***P < 0.0001 compared with the control group.

    Furthermore, we detected the mRNA levels of the top 5 upregulated and top 5 downregulated genes in the subtypes of the PBMCs. To our knowledge, the main subtypes of PBMCs consists of T cells, B cells and monocytes. Therefore, we principally explored the mRNA levels of these 10 genes in CD4+ T cells, CD8+ T cells, CD14+ monocytes and CD20+ B cells (Figure 5). EGR3 was mainly expressed in T cells in mice with NICM. MBNL, PTGS2 and EGR2 were consistently higher in monocytes; however, their expression was significantly higher in mice with NICM. TXNIP was mainly expressed in CD4+ T cells, CD8+ T cells and CD20+ B cells, the expression of which was significantly increased in the presence of NICM. Regarding the downregulated genes, we found that JUN was highly expressed in T cells in normal subjects, while its level was reduced in T cells but increased in B cells in mice with NICM. CREM was principally expressed in T cells, while its level was decreased in T cells in mice with NICM. RGS1 was primarily expressed in B cells in normal mice. The level of NPIPB3 was higher in T cells and monocytes of normal mice compared with the mice with NICM. The mRNA levels of GABARAPL3 in the four types of immune cells were consistently higher in normal mice than in mice with NICM.

    Figure 5.  Expression levels of the top 5 upregulated and top 5 downregulated DEGs in the different isolated cell types (CD4+ T cells, CD8+ T cells, CD14+ monocytes and CD20+ B cells) of PBMCs from control mice and NICM mice.

    Patients with left ventricular dysfunction resulting from NICM possess an annual rate of death from any cause of approximately 7 percent when they received standard medical therapy for HF [22]. The aetiology of NICM is complicated and multifactorial, and at present, it is recognized that NICM is usually an integrated pathologic change implicated with genetic factors, mechanical stress, metabolic disorders, and infectious (mainly viral) causes, as well as toxicity-related causes [23]. The potential pathophysiological mechanisms underlying these alterations in NICM remain incompletely explored, and it is believed that the interplay among altered sarcomeric and cytoskeletal proteins, autoimmune mechanisms, postinfection immune mechanisms, direct pathogen damage, and free oxygen radical species are associated with the development and progression of NICM [24]. In particular, immune mechanisms in NICM have received increasing attention and have been unveiled in recent years. Some therapeutic strategies have also been proposed to indirectly modulate the inflammatory and immune status of patients with NICM. However, to date, our knowledge of immune mediators and mechanisms in NICM is still too superficial. Hence, some key immune-related diagnostic biomarkers and therapeutic targets in PBMCs should be verified as early as possible.

    In this study, we first performed a comprehensive investigation on the expression profile of PBMCs obtained from patients with NICM and healthy individuals. Our study included 3 normal individuals and 4 patients with NICM from the GSE9128 GEO database. In our analysis, a total of 371 DEGs (accounting for 2.1 % of all detected genes) were found: 288 upregulated genes and 83 downregulated genes. Among the 371 DEGs, we identified the top 5 upregulated (EGR3, MBNL1, PTGS2, EGR2 and TXNIP) and top 5 downregulated (GABARAPL3, NPIPB3, RGS1, CREM and JUN) genes between healthy individuals and patients with NICM. Furthermore, we detected their mRNA levels in PBMCs in a mouse model induced by isoprenaline, and the results based on our experiments validated the bioinformatics analysis. Thus, we hypothesize that these 10 DEGs are expected to be promising biomarkers for the diagnosis of NICM. Meanwhile, based on a PPI network, we selected 15 hub genes with high connectivity, which all occupied the core node of the network, indicating that these hub genes may exert pivotal influence in the immune response of NICM. Furthermore, by annotating and analysing all DEGs in the PBMCs, we found that immune response, inflammatory response, extracellular exosome and cell differentiation were enriched in the progression of NICM. In fact, most of the DEGs, especially the top 5 upregulated/downregulated genes and the top 15 hub genes, perform vital immunoregulatory functions.

    EGR3 is a transcription factor of the Cys2His2-type zinc finger transcription factor family that can inhibit IL-2 and interferon (IFN)-γ secretion in T cells and upregulate the level of the E3 ligase Cbl-b, which is important for the modulation of T cell energy and tolerance [25]. In a NICM model induced by transverse aortic constriction (TAC), IL-2 was significantly downregulated, and pretreatment with IL-2 was effective in the selective induction of regulatory T cells (Tregs) and in alleviating the development of cardiac remodelling through decreasing cardiac systemic inflammation [26]. On the other hand, the frequency of T helper (Th) 17 cells, which are derived from naive Th cells, was significantly elevated in patients with HF [27]. Meanwhile, EGR3 could trigger a Th17 response by enhancing the development of γδ T cells, suggesting that EGR3 may mediate the activation of Th17 cells and induce HF [28]. In our study, EGR3 was significantly increased in T cells of patients and mice with NICM. Thus, we speculate that EGR3 may suppress the secretion of IL2 by inhibiting Tregs and promoting Th17 cells, thereby deteriorating cardiac injury and dysfunction.

    Excessive inflammation and monocyte/macrophage activation are featured in cardiac remodelling associated with NICM [29]. However, the roles of monocytes/macrophages in cardiac remodelling and HF remain controversial. For example, inhibition of macrophages in hypertensive rats resulted in an acceleration of cardiomyopathy, which was evidenced by left ventricle systolic dysfunction [30]. Other studies have provided primary evidence that the localization of macrophages may contribute to inflammation, fibrosis and cardiac injury induced by TAC [31,32]. RGS1 was reported to inactivate G-protein signalling and suppress the response to continuous chemokine stimulation. A recent study reported that in non-activated monocytes, RGS1 was low, the expression of which was immediately upregulated upon the activation of monocytes by pro-inflammatory signals during the recruitment phase. Subsequently, monocytes began to differentiate into macrophages, resulting in the increased expression of RGS1. The increased RGS1 further terminated chemokine signalling, eventually reducing the capacity for inflammatory cell migration and increasing the accumulation of macrophages [33]. In our study, we found that RGS1 was highly expressed in monocytes in the control group and that the mRNA level of RGS1 was significantly decreased in patients and mice with NICM, which hinted to us that RGS1 may exert a protective role in NICM by blocking the differentiation of monocytes/macrophages. In fact, another top 5 upregulated gene in our study, EGR2 was also uncovered to regulate the differentiation of monocytes into macrophages by modulating the expression of colony-stimulating factor-1 (CSF-1). In the inflammatory response, the hub gene (TLR4) has been reported to be significantly activated in monocytes and macrophages in atherosclerosis [34]. Collectively, our study suggested that the DEGs we screened participated in the differentiation of monocytes into macrophages in the development of NICM from the perspective of bioinformatics analysis.

    Another finding in our study was that the inflammatory response played an essential role in the development of NICM [35]. Although a vast number of studies have proven the biological status of inflammation in NICM, few studies have provided evidence using bioinformatics analysis. The GO and KEGG pathway analysis in our study revealed that many inflammatory pathways in patients with NICM were activated, including the NF-kappa B signaling pathway, mitochondrial uncoupling, Toll-like receptor signaling, NOD-like receptor signaling pathway and TGF-beta signaling pathway. Additionally, many of the screened DEGs and hub genes were also associated with inflammation, including TXNIP, GABARAPL3, CREM, JUN, PCNA, TLR4, IL8, IL1B and HSP90AA1. To our knowledge, activation of the immune system is not independent of inflammation in the progression of HF. In chronic HF, activating the immune system can usually contribute to the activation of the complement system, secretion of pro-inflammatory cytokines, and production and release of autoantibodies, indicating that immune-inflammatory activation in patients with NICM may occur synchronously but not in isolation [36].

    Taken together, using a series of bioinformatics analyses, we performed an integrated and novel illustration of mRNA expression profiles to identify DEGs expressed in PBMCs, which may play critical roles in the occurrence and development in patients with NICM. Genes and pathways implicated in inflammation and the immune response were significantly altered in the PBMCs of NICM. These findings will greatly contribute to unveiling the molecular mechanisms of NICM and promoting a non-invasive diagnosis. To allow these biomarkers and targets to be used more routinely in the clinic, further investigations into these genes in PBMCs should be performed.

    All procedures performed in study involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. All applicable international, national, and/or institutional guidelines for the care and use of animals were followed.

    This work was supported by grants from the National Natural Science Foundation of China (No. 81470516, 81530012, 81300070, and 81770399), the Key Project of the National Natural Science Foundation (No.81530012) and Fundamental Research Funds for the Central Universities of China (2042018kf0121).

    The data used to support the findings of this study are available from the corresponding author upon request.

    ZJ conceived the experiments. LN conducted the experiments. LN analyzed the data. ZJ wrote the paper.

    The authors declare that they have no competing interests.



    [1] R. Panchagnula, N. S. Thomas, Biopharmaceutics and pharmacokinetics in drug research, Int. J. Pharm., 201 (2000), 131–150. https://doi.org/10.1016/s0378-5173(00)00344-6 doi: 10.1016/s0378-5173(00)00344-6
    [2] M. A. Hedaya, Basic Pharmacokinetics, Routledge, 2023. https://doi.org/10.1016/b978-0-12-815499-1.00050-8
    [3] X. H. Huang, Q. S. Zheng, Pharmacokinetic and pharmacodynamic data analysis: concepts and applications, Am. J. Pharm. Educ., 74 (2010). https://doi.org/10.5688/aj740353b
    [4] C. Csajka, D. Verotta, Pharmacokinetic–pharmacodynamic modelling: history and perspectives, J. Pharmacokinet. Pharmacodyn., 33 (2006), 227–279. https://doi.org/10.1007/s10928-005-9002-0 doi: 10.1007/s10928-005-9002-0
    [5] W. Wang, E. Q. Wang, J. P. Balthasar, Monoclonal antibody pharmacokinetics and pharmacodynamics, Clin. Pharmacol. Ther., 84 (2008), 548–558. https://doi.org/10.1038/clpt.2008.170 doi: 10.1038/clpt.2008.170
    [6] Y. M. C. Wang, B. Sloey, T. Wong, P. Khandelwal, R. Melara, Y. N. Sun, Investigation of the pharmacokinetics of romiplostim in rodents with a focus on the clearance mechanism, Pharm. Res., 28 (2011), 1931–1938. https://doi.org/10.1007/s11095-011-0420-y doi: 10.1007/s11095-011-0420-y
    [7] S. Kozawa, N. Yukawa, J. Liu, A. Shimamoto, E. Kakizaki, T. Fujimiya, Effect of chronic ethanol administration on disposition of ethanol and its metabolites in rat, Alcohol, 41 (2007), 87–93. https://doi.org/10.1016/j.alcohol.2007.03.002 doi: 10.1016/j.alcohol.2007.03.002
    [8] P. H. van der Graaf, N. Benson, L. A. Peletier, Topics in mathematical pharmacology, J. Dyn. Differ. Equations, 28 (2016), 1337–1356. https://doi.org/10.1007/s10884-015-9468-4 doi: 10.1007/s10884-015-9468-4
    [9] H. Moore, R. Allen, What can mathematics do for drug development?, Bull. Math. Biol., 81 (2019), 3421–3424. https://doi.org/10.1007/s11538-019-00632-x doi: 10.1007/s11538-019-00632-x
    [10] S. Tang, Y. Xiao, One-compartment model with Michaelis-Menten elimination kinetics and therapeutic window: an analytical approach, J. Pharmacokinet. Pharmacodyn., 34 (2007), 807–827. https://doi.org/10.1007/s10928-007-9070-4 doi: 10.1007/s10928-007-9070-4
    [11] A. Dokoumetzidis, R. Magin, P. Macheras, Fractional kinetics in multi-compartmental systems, J. Pharmacokinet. Pharmacodyn., 37, (2010), 507–524. https://doi.org/10.1007/s10928-010-9170-4 doi: 10.1007/s10928-010-9170-4
    [12] X. Wu, J. Li, F. Nekka, Closed form solutions and dominant elimination pathways of simultaneous first-order and Michaelis–Menten kinetics, J. Pharmacokinet. Pharmacodyn., 42 (2015), 151–161. https://doi.org/10.1007/s10928-015-9407-3 doi: 10.1007/s10928-015-9407-3
    [13] X. Wu, F. Nekka, J. Li, Analytical solution and exposure analysis of a pharmacokinetic model with simultaneous elimination pathways and endogenous production: The case of multiple dosing administration, Bull. Math. Biol., 81 (2019), 3436–3459. https://doi.org/10.1007/s11538-019-00651-8 doi: 10.1007/s11538-019-00651-8
    [14] N. A. Daryakenari, M. De Florio, K. Shukla, G. E. Karniadakis, Al-Aristotle: A physics-informed framework for systems biologygray-box identification, PLoS Comput. Biol., 20 (2024), e1011916. https://doi.org/10.1371/journal.pcbi.1011916 doi: 10.1371/journal.pcbi.1011916
    [15] N. A. Daryakenari, S. Wang, G. Karniadakis, CMINNs: Compartment model informed neural networks-Unlocking drug dynamics, Comput. Biol. Med., 184, (2025), 109392. https://doi.org/10.1016/j.compbiomed.2024.109392 doi: 10.1016/j.compbiomed.2024.109392
    [16] R. Panaccione, S. Ghosh, S. Middleton, J. R. Márquez, B. B. Scott, L. Flint, et al., Combination therapy with infliximab and azathioprine is superior to monotherapy with either agent in ulcerative colitis, Gastroenterology, 146 (2014), 392–400. https://doi.org/10.1053/j.gastro.2013.10.052 doi: 10.1053/j.gastro.2013.10.052
    [17] N. M. Maruthur, E. Tseng, S. Hutfless, L. M. Wilson, C. Suarez-Cuervo, Z. Berger, et al., Diabetes medications as monotherapy or metformin-based combination therapy for type 2 diabetes: a systematic review and meta-analysis, Ann. Intern. Med., 164 (2016), 740–751. https://doi.org/10.7326/M15-2650 doi: 10.7326/M15-2650
    [18] Y. Ma, Q. Chen, Y. Zhang, J. Xue, Q. Liu, Y. Zhao, et al., Pharmacokinetics, safety, tolerability, and feasibility of apatinib in combination with gefitinib in stage IIIB-IV EGFR-mutated non-squamous NSCLC: a drug-drug interaction study, Cancer Chemother. Pharmacol., 92 (2023), 411–418. https://doi.org/10.1007/s00280-023-04563-2 doi: 10.1007/s00280-023-04563-2
    [19] L. Hanum, Q. Ertiningsih, N. Susyanto, Sensitivity analysis unveils the interplay of drug-sensitive and drug-resistant Glioma cells: Implications of chemotherapy and anti-angiogenic therapy, Electron. Res. Arch., 32 (2024), 72–89. https://doi.org/10.3934/era.2024004 doi: 10.3934/era.2024004
    [20] X. Tan, S. Fan, K. Duan, M. Xu, J. Zhang, P. Sun, et al., A novel drug-drug interactions prediction method based on a graph attention network, Electron. Res. Arch., 31 (2023), 5632–5648. https://doi.org/10.3934/era.2023286 doi: 10.3934/era.2023286
    [21] N. Erawaty, N. Aris, A mathematical study of effects of delays arising from the interaction of anti-drug antibody and therapeutic protein in the immune response system, AIMS Math., 5 (2020), 7176–7199. https://doi.org/10.3934/math.2020460 doi: 10.3934/math.2020460
    [22] T. Gould, R. J. Roberts, Therapeutic problems arising from the use of the intravenous route for drug administration, J. Pediatr., 95 (1979), 465–471. https://doi.org/10.1016/s0022-3476(79)80538-7 doi: 10.1016/s0022-3476(79)80538-7
    [23] P. Fasinu, V. Pillay, V. M. K. Ndesendo, L. C. du Toit, Y. E. Choonara, Diverse approaches for the enhancement of oral drug bioavailability, Biopharm. Drug Dispos., 32 (2011), 185–209. https://doi.org/10.1002/bdd.750 doi: 10.1002/bdd.750
    [24] R. Selimov, E. Goncharova, P. Koriakovtsev, D. Gabidullina, J. Karsakova, S. Kozlov, et al., Comparative pharmacokinetics of florfenicol in heifers after intramuscular and subcutaneous administration, J. Vet. Pharmacol. Ther., 46 (2023), 177–184. https://doi.org/10.1111/jvp.13110 doi: 10.1111/jvp.13110
    [25] N. Z. Kerin, R. D. Blevins, H. Frumin, K. Faitel, M. Rubenfire, Intravenous and oral loading versus oral loading alone with amiodarone for chronic refractory ventricular arrhythmias, Am. J. Cardiol., 55 (1985), 89–91. https://doi.org/10.1016/0002-9149(85)90305-4 doi: 10.1016/0002-9149(85)90305-4
    [26] H. U. Xiao-hu, T. San-yi, Approximate Solutions to the Nonlinear Compartmental Model for Extravascular Administration, Appl. Math. Mech., 35 (2014), 1033–1045. https://doi.org/10.3879/j.issn.1000-0887.2014.09.009 doi: 10.3879/j.issn.1000-0887.2014.09.009
    [27] S. H. Jang, P. M. Colangelo, J. V. S. Gobburu, Exposure-response of posaconazole used for prophylaxis against invasive fungal infections: evaluating the need to adjust doses based on drug concentrations in plasma, Clin. Pharmacol. Ther., 88 (2010), 115–119. https://doi.org/10.1038/clpt.2010.64 doi: 10.1038/clpt.2010.64
    [28] M. Prasad, P. R. Krishnan, R. Sequeira, Khaldoon. Al-Roomi, Anticonvulsant therapy for status epilepticus, Cochrane Database Syst. Rev., 2014. https://doi.org/10.1002/14651858.CD003723.pub3
    [29] E. H. Reynolds, D. Chadwick, A. W. Galbraith, One drug(phenytoin) in the treatment of epilepsy, The lancet, 307 (1976), 923–926. https://doi.org/10.1016/s0140-6736(76)92881-6 doi: 10.1016/s0140-6736(76)92881-6
    [30] J. Patocka, Q. Wu, E. Nepovimova, K. Kuca, Phenytoin–An anti-seizure drug: Overview of its chemistry, pharmacology and toxicology, Food Chem. Toxicol., 142 (2020), 111393. https://doi.org/10.1016/j.fct.2020.111393 doi: 10.1016/j.fct.2020.111393
    [31] W. J. Jusko, J. R. Koup, G. Alván, Nonlinear assessment of phenytoin bioavailability, J. Pharmacokinet. Biopharm., 4 (1976), 327–336. https://doi.org/10.1007/bf01063122 doi: 10.1007/bf01063122
    [32] R. Gugler, C. V. Manion, D. L. Azarnoff, Phenytoin: pharmacokinetics and bioavailability, Clin. Pharmacol. Ther., 19 (1976), 135–142. https://doi.org/10.1002/cpt1976192135 doi: 10.1002/cpt1976192135
    [33] L. Osterberg, T. Blaschke, Adherence to medication, N. Engl. J. Med., 353 (2005), 487–497. https://doi.org/10.7748/ns.25.2.59.s53 doi: 10.7748/ns.25.2.59.s53
    [34] B. Vrijens, S. De Geest, D. A. Hughes, K. Przemyslaw, J. Demonceau, T. Ruppar, et al., A new taxonomy for describing and defining adherence to medications, Br. J. Clin. Pharmacol., 73 (2012), 691–705. https://doi.org/10.1111/j.1365-2125.2012.04167.x doi: 10.1111/j.1365-2125.2012.04167.x
    [35] M. E. Csete, J. C. Doyle, Reverse engineering of biological complexity, Science, 295 (2002), 1664–1669. https://doi.org/10.1126/science.1069981 doi: 10.1126/science.1069981
    [36] B. Munos, Lessons from 60 years of pharmaceutical innovation, Nat. Rev. Drug Discov., 8 (2009), 959–968. https://doi.org/10.1038/nrd2961 doi: 10.1038/nrd2961
  • This article has been cited by:

    1. Hongru Ma, Yanbin Tang, Homogenization of a semilinear elliptic problem in a thin composite domain with an imperfect interface, 2023, 46, 0170-4214, 19329, 10.1002/mma.9628
    2. Yongqiang Zhao, Yanbin Tang, Approximation of solutions to integro-differential time fractional order parabolic equations in Lp-spaces, 2023, 2023, 1029-242X, 10.1186/s13660-023-03057-2
    3. Caihong Gu, Yanbin Tang, Well-posedness of Cauchy problem of fractional drift diffusion system in non-critical spaces with power-law nonlinearity, 2024, 13, 2191-950X, 10.1515/anona-2024-0023
    4. Junlong Chen, Yanbin Tang, Homogenization of non-local nonlinear p-Laplacian equation with variable index and periodic structure, 2023, 64, 0022-2488, 10.1063/5.0091156
    5. Feiyang Peng, Yanbin Tang, Inverse problem of determining diffusion matrix between different structures for time fractional diffusion equation, 2024, 19, 1556-1801, 291, 10.3934/nhm.2024013
    6. Hongru Ma, Yanbin Tang, Homogenization of a nonlinear reaction‐diffusion problem in a perforated domain with different size obstacles, 2024, 104, 0044-2267, 10.1002/zamm.202300333
    7. Tomasz Klimsiak, Tomasz Komorowski, Lorenzo Marino, Homogenization of stable-like operators with random, ergodic coefficients, 2025, 00220396, 10.1016/j.jde.2025.02.054
  • Reader Comments
  • © 2025 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(328) PDF downloads(37) Cited by(0)

Figures and Tables

Figures(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog