Loading [MathJax]/jax/element/mml/optable/MathOperators.js
Research article Special Issues

Dynamics of T cell receptor distributions following acute thymic atrophy and resumption

  • Naive human T cells are produced and developed in the thymus, which atrophies abruptly and severely in response to physical or psychological stress. To understand how an instance of stress affects the size and "diversity" of the peripheral naive T cell pool, we derive a mean-field autonomous ODE model of T cell replenishment that allows us to track the clone abundance distribution (the mean number of different TCRs each represented by a specific number of cells). We identify equilibrium solutions that arise at different rates of T cell production, and derive analytic approximations to the dominant eigenvalues and eigenvectors of the mathematical model linearized about these equilibria. From the forms of the eigenvalues and eigenvectors, we estimate rates at which counts of clones of different sizes converge to and depart from equilibrium values-that is, how the number of clones of different sizes "adjusts" to the changing rate of T cell production. Under most physiological realizations of our model, the dominant eigenvalue (representing the slowest dynamics of the clone abundance distribution) scales as a power law in the thymic output for low output levels, but saturates at higher T cell production rates. Our analysis provides a framework for quantitatively understanding how the clone abundance distribution evolves under small changes in the overall T cell production rate.

    Citation: Stephanie M. Lewkiewicz, Yao-Li Chuang, Tom Chou. Dynamics of T cell receptor distributions following acute thymic atrophy and resumption[J]. Mathematical Biosciences and Engineering, 2020, 17(1): 28-55. doi: 10.3934/mbe.2020002

    Related Papers:

    [1] Hengki Tasman, Edy Soewono, Kuntjoro Adji Sidarto, Din Syafruddin, William Oscar Rogers . A model for transmission of partial resistance to anti-malarial drugs. Mathematical Biosciences and Engineering, 2009, 6(3): 649-661. doi: 10.3934/mbe.2009.6.649
    [2] Alberto Gambaruto, João Janela, Alexandra Moura, Adélia Sequeira . Shear-thinning effects of hemodynamics in patient-specific cerebral aneurysms. Mathematical Biosciences and Engineering, 2013, 10(3): 649-665. doi: 10.3934/mbe.2013.10.649
    [3] Yuyang Xiao, Juan Shen, Xiufen Zou . Mathematical modeling and dynamical analysis of anti-tumor drug dose-response. Mathematical Biosciences and Engineering, 2022, 19(4): 4120-4144. doi: 10.3934/mbe.2022190
    [4] Mario Grassi, Giuseppe Pontrelli, Luciano Teresi, Gabriele Grassi, Lorenzo Comel, Alessio Ferluga, Luigi Galasso . Novel design of drug delivery in stented arteries: A numerical comparative study. Mathematical Biosciences and Engineering, 2009, 6(3): 493-508. doi: 10.3934/mbe.2009.6.493
    [5] Alberto M. Gambaruto, João Janela, Alexandra Moura, Adélia Sequeira . Sensitivity of hemodynamics in a patient specific cerebral aneurysm to vascular geometry and blood rheology. Mathematical Biosciences and Engineering, 2011, 8(2): 409-423. doi: 10.3934/mbe.2011.8.409
    [6] Weirui Lei, Jiwen Hu, Yatao Liu, Wenyi Liu, Xuekun Chen . Numerical evaluation of high-intensity focused ultrasound- induced thermal lesions in atherosclerotic plaques. Mathematical Biosciences and Engineering, 2021, 18(2): 1154-1168. doi: 10.3934/mbe.2021062
    [7] Dehua Feng, Xi Chen, Xiaoyu Wang, Xuanqin Mou, Ling Bai, Shu Zhang, Zhiguo Zhou . Predicting effectiveness of anti-VEGF injection through self-supervised learning in OCT images. Mathematical Biosciences and Engineering, 2023, 20(2): 2439-2458. doi: 10.3934/mbe.2023114
    [8] Gesham Magombedze, Winston Garira, Eddie Mwenje . Modelling the immunopathogenesis of HIV-1 infection and the effect of multidrug therapy: The role of fusion inhibitors in HAART. Mathematical Biosciences and Engineering, 2008, 5(3): 485-504. doi: 10.3934/mbe.2008.5.485
    [9] Avner Friedman, Kang-Ling Liao . The role of the cytokines IL-27 and IL-35 in cancer. Mathematical Biosciences and Engineering, 2015, 12(6): 1203-1217. doi: 10.3934/mbe.2015.12.1203
    [10] Xiaoying Wang, Xingfu Zou . Pattern formation of a predator-prey model with the cost of anti-predator behaviors. Mathematical Biosciences and Engineering, 2018, 15(3): 775-805. doi: 10.3934/mbe.2018035
  • Naive human T cells are produced and developed in the thymus, which atrophies abruptly and severely in response to physical or psychological stress. To understand how an instance of stress affects the size and "diversity" of the peripheral naive T cell pool, we derive a mean-field autonomous ODE model of T cell replenishment that allows us to track the clone abundance distribution (the mean number of different TCRs each represented by a specific number of cells). We identify equilibrium solutions that arise at different rates of T cell production, and derive analytic approximations to the dominant eigenvalues and eigenvectors of the mathematical model linearized about these equilibria. From the forms of the eigenvalues and eigenvectors, we estimate rates at which counts of clones of different sizes converge to and depart from equilibrium values-that is, how the number of clones of different sizes "adjusts" to the changing rate of T cell production. Under most physiological realizations of our model, the dominant eigenvalue (representing the slowest dynamics of the clone abundance distribution) scales as a power law in the thymic output for low output levels, but saturates at higher T cell production rates. Our analysis provides a framework for quantitatively understanding how the clone abundance distribution evolves under small changes in the overall T cell production rate.



    Abdominal aortic aneurysm (AAA) is a leading cause of death globally and involves complex and multi-factorial pathogenesis including factors such as hypertension, dyslipidemia, and tobacco use. Half of patients with untreated aneurysms larger than 5.5 cm will die of rupture within five years, and there is currently no recognized treatment for these AAAs. Patients with AAAs frequently have atherosclerosis, and numerous studies show the association of coronary heart disease and peripheral atherosclerosis with AAA [1]. Whether this association between AAA and atherosclerosis is causal or simply due to common risk factors is unknown. However, the term "atherosclerotic aneurysms" is commonly used and prevalent today. As one of most common causes of AAA, atherosclerosis originates with an accumulation of blood-borne lipids in the inner-most arterial wall (intima), followed by a detailed immune response and the progressive formation of a lesion in the arterial wall [2]. Mature lesions then form atherosclerotic plaques, then arterial walls lose the elasticity via arterial remodeling and develop the AAA. Thus, it is vital to understand the pathophysiology of atherosclerotic aneurysms whose pathology is largely unknown and likely multifactorial [3,4]. Though the pathophysiology is hard to study experimentally due to its network complexity, a computational model could integrate this complex network and include all the possible factors in the pathogenesis. Moreover, computational modeling techniques, such as fluid-structure interactions or computational fluid dynamics, has been successfully implemented in the study of CVD and provides many insightful suggestions for clinical practice [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21].

    In this paper, we develop a mathematical model that accounts not just for the biological networks associated with the pathophysiology of atherosclerotic aneurysms, but also for the multi-layered structure of the arterial wall. Modeling one or two layers would not be sufficient for capturing the progression of the atherosclerotic aneurysm. Therefore our model accounts for the multi-layered structure of the arterial walls, with attention paid to the distinct physical properties of each layer. Biologically speaking, atherosclerosis silently and slowly hardens, narrows, and blocks the arteries, putting blood flow at risk [22]. Then atherosclerotic plaques restrict blood vessels (a process called stenosis) and lead to the destruction of smooth muscle cells (SMCs), depletion of elastin, and disruption of the extracellular matrix (ECM) in the media/adventitia which causes the AAA. Our mathematical model, by including both the mechanism of the arterial wall and the biology of the atherosclerotic aneurysm, can help us to construct the cardiovascular environment virtually and to mimic plaque/aneurysm growth. This is an innovative tool for potential biomarkers' detection, drug efficacy evaluation, and reducing the risk of misdiagnosis due to patient variations. More specifically, we use our model to assess the risk of atherosclerotic aneurysms on a 2D biomarker, the cholesterol ratio and the level of the receptor tyrosine kinase, discoidin domain receptor 1 (DDR1). We further explore the efficacy of three treatments, anti-cholesterol drugs, the anti-DDR1 injection, and the combined therapy. Our simulation results suggest that the dosage of anti-cholesterol drugs is critical for treating atherosclerotic aneurysm patients who can also benefit from additional anti-DDR1 injections.

    In order to model the atherosclerotic aneurysm, we present a hybrid computational model which consists of a mechanical model, including the psychical properties and multi-layer structure of arterial wall, and a biological model, including the molecular mechanism of etiology and pathophysiology of atherosclerotic aneurysm.

    The arterial wall has a multi-layer structure (see Figure 1), comprising intima, media, and adventitia. [23]. This multi-layered structure is central to the progression atherosclerotic aneurysm which involves changes in all three arterial layers, including changes to their geometry and composition. Modeling the multi-layered structure makes our model possible to render an accurate virtual cardiovascular environment, even recreating and observing the pathophysiological networks of both atherosclerosis (affecting the intima) and aneurysm formation (affecting the media/adventitia).

    Figure 1.  Right: Detail of the arterial wall's three concentric layers: adventitia, media, and intima. Left: Schematic representation of a 2D section of the computational domain: ΩI represents the intima and ΩA is the media/adventitia. There are three free boundaries: The inner and outer surfaces of the arterial wall (ΓI and ΓA), and the surface between the intima and media/adventitia (ΓM).

    Porous medium for intima: In light of the intima's high permeability to white cells and platelets [24,25,26], we propose to treat the intima as a porous medium and all the cells in the intima as a low speed flow [27,28] moving with a common velocity, v, which is the result of the movement of macrophages, T cells and SMCs in the intima. The pressure, p, results from the cells' movement and proliferation. This leads to the following Darcy's law and mass conservation equation [27,28]:

    {div(v)=f(C,S,M)andv+Kμp=0inΩI(t)p=pB+γκonΓI(t),pn=σn|ΩA+γκnonΓM(t), (2.1)

    where

    f(C,S,M)=λMM+λMCC+λSMMKM+M+λCMMdMMdSSdCCKC+KS+KM,

    which is derived by the biological model for which variables M, S, and C refer to macrophages, SMCs, and collagen, respectively. See next section for more details. Here K is the Darcy permeability of the intima layer, μ is the viscosity of the cell fluid in the intima, pB is the blood pressure (we treat as a constant since we do not include the blood fluid dynamics), γ is a dimensionless surface tension coefficient, κ is the mean curvature of the corresponding free boundaries and σ is the stress tensor determined by the media/adventitia in ΩA(t).

    Elasticity for media/adventitia: From the continuum mechanics point of view, mechanical response of the media/adventitia layers of an arterial wall, is a elastic material [29,30]. Specifically, the Cauchy stress tensor is given by

    σ=λtr(ε)I+2μsε, (2.2)

    where μs and λ are material constants, called Lame or elastic moduli coefficients, and ε is the so-called strain stressor defined as

    ε=sym(u)=12(u+uT), (2.3)

    where u is the displacement. Therefore,

    σ=λuI+μs(u+uT). (2.4)

    The coefficients are directly related to Young's modulus, E, and Poisson's ratio, ν, namely, λ=Eν(12ν)(1+ν) and μs=E1+ν. In our model, E=E0β(SS01) where S is the SMCs density and produces the elastic modulus. Then the momentum equation is given by

    {dDttuσ=0 in ΩA(t)v=tuonΓM(t)σn=0onΓA(t), (2.5)

    where d is the tissue density. Here we impose the stress-free boundary condition on ΓA and the velocity continuous condition on ΓM(t).

    Pathophysiology of atherosclerosis (intima): The atherosclerosis involves a complex series of biological processes with the formation of plaque inside the intima layer. This process begins with a lesion, allowing low-density lipoprotein (LDL), to move from the blood into the intima, where LDL then becomes oxidized (ox-LDL) by free radicals [31]. Because ox-LDL is toxic, macrophages are recruited and differentiated in the intima [32,28]. The ingestion of large amounts of ox-LDL transforms the fatty macrophages into foam cells, which contribute to plaque formation [33,34,35,36]. Moreover, The presence of collagen in the intima also causes activation of macrophages via the DDR1, which is a member of the receptor tyrosine kinase family and binds to several collagen(s) as its ligand, including collagen I, III and IV present in the vasculature [37]. With regard to vascular diseases, it is important to note that DDR1 is expressed in monocytes and macrophages, endothelial cells, SMCs as well as in adventitial fibroblasts. It can thus have a multifaceted role in regulation of vascular diseases such as atherosclerosis and aneurysms. Studies using vessel wall or bone-marrow specific deletion of DDR1 in mice models have elucidated on how DDR1 expression on macrophages mediates their infiltration into the developing plaque, whereas DDR1 expression in resident vessel wall SMCs limits SMC infiltration and lesion fibrosis [38]. The biological network is shown in Figure 2.

    Figure 2.  The schematic network of atherosclerotic aneurysm. In atherosclerosis, foam cells and ECM in the intima contribute to the buildup of plaque. Meanwhile, SMCs and ECM in the media/adventitia play major roles in aneurysm formation.

    Pathophysiology of aneurysm growth (media/adventitia): Aneurysm is due to the destruction of smooth muscle cells (SMCs), depletion of elastin in the media, and disruption of the extracellular matrix (ECM) in the adventitia. Foam cells in the plaque [39,40], degrade elastin and ECM to weaken the pressure of the adventitia layer to withstand blood stress [41]. Foam cells are also known to cause apoptosis in SMCs, which leads to a reduction in elastin and subsequent overall decline of the media's strength [39,40,42].

    Equation for Macrophages (M): Activated ECs secrete cytokines like MCP-1 to recruit monocytes to the fatty streak area [43]. The collagens in the intima also causes activation of macrophages. Then the equation of M is written as

    {Mt+(vM)DMΔM=λMM+λMCCdMM in ΩI(t)Mnα(C)M0=0onΓI(t)Mn=0onΓM(t), (2.6)

    where λMM is the macrophages proliferation due to the scavenger receptor class A (SR-A) expressed by tissue macrophages [44,45], λMCC accounts for the macrophages activation by the DDR1 [46], dMC represents the natural death of macrophages, α(C)=α0+αCKC+C for which α0 indicates the MCP-1 induced by the oxidized LDL [47] and αCKC+C describes the production of MCP-1 activated by collagens via the DDR1 [48], and M0 is the monocytes density in the blood.

    Equation for SMCs (S): Foam cells secrete large amounts of cytokines and growth factors which stimulate the migration of SMCs from the medial layer of the artery into the intima [49]. Then the equation of S is written as

    {St+(tuS)DSΔS=λSχΩA+λSMMKM+MχΩIdSS in Ω(t)Sn=0onΓIΓASn+α(M)S=0onΓM(t), (2.7)

    where χΩA=1 in ΩA, χΩA=0 in ΩI. S=S0=λSdS represents the SMCs density in health, λSMMKM+M describes the SMCs proliferation due to the stimulation from platelet-derived growth factor (PDGF) produced by activated macrophages [50], and α(M)=αMKM+M accounts for the migration of SMCs due to the foam cells. Since no migration of SMCs occurs on the boundaries ΓI and ΓA, nonflux boundary conditions are specified on them.

    Equation for collagen (C): Foam cells form extracellular deposits of lipid material, known as the lipid or necrotic core in the intima while SMCs produce elastin. Thus we have

    {Ct+(tuC)DCΔC=λCSSχΩA+λCMMχΩIdCC in Ω(t)[Cn]=α(S0S)onΓM(t)Cn=0onΓAΓI(t), (2.8)

    where on ΓM(t), since the 'synthetic' SMCs secrete a matrix rich elastin in the initima [51], we use S0S account for the 'synthetic' SMCs and model the production through a flux boundary condition.

    We cover the healthy case in our model by setting α=0 and α0=0 in Eq (2.6) which means that there is no inflammation in the arterial wall. Moreover, we assume α is proportional to the ratio of the cholesterol levels, LDL/HDL, and refer α=0 in our model to the cholesterol ratio less than 5 [52,53]. Thus we have the density of macrophages in Eq (2.6) be zero. Therefore, there is no migration of SMCs on the interface between the intima and the media in Eq (2.7) which implies that the density of SMCs achieves the steady-state in health, namely, S=S0. Similarly, the density of collagen in Eq (2.8) has steady states C=KC in ΩA and C=0 in ΩI. Thus we have v=0 in Eq (2.1) and u=0 in Eq (2.5). Therefore, there is no plaque formation and no deformation of the arterial wall for the healthy case.

    We simulate the hybrid model in the 2D radially symmetric case with finite difference method and moving mesh method [54]. More specifically, for given a function f, we have the following formulas in the polar coordinate with radial symmetry:

    f=1r(rfr)r and f=frr,

    where fr is the projection of the function f on er. By denoting grid points as r0<r1<rN and the numerical approximation of f(ri) as fi, we employ the following numerical schemes to approximate the first and second derivatives with respect to r:

    frfi(h2i+1h2i1)+fi+1h2i1fi1h2i+1hi+1(h2i1+hi+1hi1),2fr22fi(hi+1hi1)+2fi+1hi12fi1hi+1hi+1(h2i1+hi+1hi1),

    where hi+1=ri+1ri and hi1=ri1ri. If the grid point is uniformly distributed, these numerical schemes recover the central difference schemes with the second order accuracy by setting hi+1=hi1. Due to the free boundaries, we use the moving mesh method to update the grid points [54], namely, ri(t+Δt)=ri(t)+vrΔt. We outline the algorithm of solving the hybrid system with free boundaries in Algorithm 1.

    Algorithm 1: The algorithm to solve the hybrid coupled system (Eqs 2.1-2.8)

     | Show Table
    DownLoad: CSV

    In order to ensure the numerical stability, we need to satisfy the CFL condition [55], namely, max|v|ΔtΔhC, where C is the so-called Courant number between 0 and 1. For simplicity, we used a small fixed Δt in the numerical computation to ensure this stability condition. More specifically, we choose the time step Δt=0.01 day and the initial grid size Δh=0.01 cm, for which the choice has been tested in the simulations to satisfy the CFL condition. The initial conditions have also been specified as:

    M=0g/cm3,S=S0, and C=KC.

    Since the thickness of the abdominal aorta is measured 2.77 cm ± 0.20 cm [56], we set the diameter of lumen, the thickness of intima layer, and the thickness of media layer to be 2 cm, 0.1 cm, and 0.2 cm, respectively, and simulate our model by setting α=1 and α0=0.1 in Eq (2.6) for 2 years. Figure 3 shows the density distributions of macrophages in the intima and SMCs in the media every three months. The macrophage density starts increasing from the boundary ΓI and then spreads to the boundary ΓM in the first six months. Then SMCs start decreasing from month 9-15 and triggers the expansion of the aneurysm after 1 year while the plaque growth starts after 6 months. This is consistent with the pathology of atherosclerotic aneurysm: the inflammation associated with atherosclerosis first leads to destruction and plaque formation then weakening of vascular wall, specifically the tunica media, which leads to aneurysmal dilation of the vessel [57]. In Figure 4 Left, the diameter of the aneurysm grows from 2.60 to 3.48 cm, thus the average aneurysm growth rate is 4.4 mm/year which is in the 95% range of the aneurysm diameter growth rate, 1.0 to 6.1 mm/year [58]. The plaque area is computed by π(R2MR2I), where RI and RM represents the radius of intima and media, respectively, and grows from 2.72 to 20.37 cm2, thus the growth rate from Y1 to Y2 is around 100%.

    Figure 3.  The density distributions of macrophage/SMCs for every three months over 2 years. The color is associated with the density in terms of g/cm3. We take α=1 and α0=0.1 in Eq (2.6) and other parameters in Table 1. The red and blue circles represent the free boundaries ΓA(t) and ΓI(t) respectively.
    Figure 4.  Left: The change of the aneurysm diameter (red) and the total plaque area (blue) v.s. years; Right: The effect of blood pressure on the aneurysm diameter and the plaque area. We take α=1 and α0=0.1 in Eq (2.6) and blood pressure pB=120 mmHg in the left figure and vary pB from 90 mmHg to 200 mmHg in the right figure.
    Table 1.  Parameters' description and value: we take λMC=102 which is twice the production rate in [28] and take λM=0.1λMC. We also assume that the production rate of SMCs by PDGF is five times the production rate of macrophages, namely, λSM=5λM and that KS=13S0. All the other parameters are taken from related references.
    ParameterDescriptionValue
    DMdiffusion coefficient of cells 8.64×104 cm2 day1
    [27,62,28,63]
    DCdiffusion coefficient of T-cell 8.64×104 cm2 day1
    [27,62,28,63]
    DSdiffusion coefficient of SMCs 8.64×104 cm2 day1
    [27,28]
    λMproduction rate of macrophages by SR-A 1×103 day 1
    λMCactivation rate of macrophages by DDR1 1×102 day1
    [27,28]
    λSproduction rate of SMCs 5.16×101 g/cm3 day1
    [54]
    λSMactivation rate of SMCs by PDGF 5×103 g/cm3 day1
    λCMactivation rate of ECM due to macrophages 9×102 day1
    [28]
    λCSactivation rate of ECM due to SMCs 6.17×102 day1
    [54]
    dM death rate of macrophages 1.5×102 day1
    [27,62,28]
    dC degradation rate of ECM 3.7×101 day1
    [27,62,28]
    dS death rate of SMCs 8.6×101 day1
    [27,62,28]
    M0source/influx of macrophages from blood 5×103 g/cm3
    [28]
    S0initial steady state of SMCs 6×101 g/cm3
    [28]
    KC ECM saturation 1×102 g/cm3
    [28]
    KS SMCs saturation 2×101 g/cm3
    KM macrophage saturation 5×101 g/cm3
    [28]
    E0Young's modulus in health 5×102 kpa
    [65,64]
    β Coefficient of aneurysm tissue material properties 8×103 kpa
    [54,66]
    νPoisson's ratio 4.9×101
    [65]

     | Show Table
    DownLoad: CSV

    Hypertension is a risk factor for the development of atherosclerosis [59] and also has long been suspected to increase the growth rate of abdominal aortic aneurysms [60], although the mechanisms have not been well elucidated. In our model, we studied the effect of blood pressure on the growth of the atherosclerotic aneurysm. The systolic blood pressure corresponds to pB in our model via Eq (2.1) which is taken to be 120 mmHg as the normal blood pressure in Figure 3. We simulate our model for different blood pressure with pB ranging from 90 mmHg (hypotension) to 200 mmHg (hypertensive crisis) [61] and show the results in Figure 4 Right. The effect of blood pressure on the aneurysm diameter varies from 3.46 to 3.51 cm (1.45%), while the effect on the plaque area varies from 20.30 to 20.53 cm2 (1.13%). The effect of blood pressure in our model is not significant due to the early stage of the atherosclerotic aneurysm and also gets larger as the time evolves shown in Figure 4. Moreover, hypertension also induces oxidative stress on the arterial wall [59] and the inflammatory cell infiltration [60]. However, we do not include these mechanisms in the current model so we might underestimate the effect of blood pressure on atherosclerotic aneurysm.

    The LDL-cholesterol to HDL-cholesterol (LDL/HDL-cholesterol) ratio is recognized as a good biomarker in atherosclerosis [67]. Although there are no effective biomarkers for aneurysm [68], the level of DDR1 is highly associated with aneurysm growth in mice experiments [69,70]. We will explore the possibility of using the cholesterol ratio and DDR1 as a 2D biomarkers for atherosclerotic aneurysm via our model. As in section 3.1, we refer the cholesterol ratio as the flux rate, α, in Eq (2.6).

    The DDR1 is deleted exclusively in vessel wall cells, macrophages can colonize and initiate an inflammatory response in the arterial intima of both Ddr1+/+ and Ddr1-/- hosts. However, vessel-wall specific deletion of DDR1 results in (over 3-fold) increase in the size of aortic sinus plaques due to enhanced collagen accumulation. The cellular composition of the lesions in DDR1-/- hosts shifted toward increased numbers of vessel wall-derived SMCs compared to bone marrow-derived macrophages. In vitro studies confirmed that Ddr1(-/-) SMCs expressed more matrix, proliferated more, and migrated farther than Ddr1(+/+) SMCs. [38]. One plausible cause is that in the absence of DDR1, N-cadherin mediated cell-cell contacts in vSMCs are mislocalised from cell-cell junctions [71], thus facilitating cell invasion. As Ddr1-/- SMCs invade the plaque, their enhanced expression of interstitial matrix genes results in extracellular matrix synthesis and assembly which contributes directly to plaque expansion.

    Moreover, DDR1 is also reported to be a negative regulator of collagen(s) in SMCs. The mRNA for collagen type I (and III) was significantly increased in DDR1-/- SMC as compared to DDR1+/+ cells [38,70]. Consistent with these results, overexpression by transient transcription of DDR1 in human SMC was found to induce a significant decrease in collagen type I mRNA, which was accompanied by a similar decrease at the protein level [72]. Besides decreasing the mRNA levels, expression of DDR1 on the cell surface is also reported to inhibit collagen. These multiple effects of DDR1 in suppressing collagen deposition account for increased collagen deposition observed in the adventitial layer of DDR1 KO mice [69]. Therefore we model the level of DDR1 as follows:

    ● Since DDR1 limits SMC migration [38] and regulates the inflammation [70], α(M) in Eq (2.7) becomes α(M)=αDDR1(MKM+MdDDR1), where αDDR1 is the SMC migration and dDDR1 accounts for the inflammation regulation;

    ● Since DDR1 increases SMC-mediated collagen and elastin degradation [73] and is also a central mediator of SMC migration [74], α in Eq (2.8) becomes α=αDDR1.

    Thus we define the risk of atherosclerotic aneurysm associated with these two biomarkers as

    R(L/H,DDR1)=The outer diameter of the arteryThe inner diameter of the artery,

    where the LDL/HDL ratio, L/H, varies from 5 to 15 [67], the level of DDR1 varies from 5 to 15 μm/ml [70], and the outer/inner diameters of the artery is measured at year 2. The risk, R, combines both the aneurysm growth (the outer diameter) and the plaque size (the inner diameter) and is a good indicator of the progression of the combined disease. The contour plot of the risk map is shown in Figure 5 with two contour curves 1.5 and 2. We have a high-risk region (above 2), a medium-risk region (between 1.5 and 2), and a low-risk region (below 1.5). The DDR1 level contributes significantly to the risk map: if the DDR1 level is low, no matter how large the cholesterol ratio is, the risk remains low; if the DDR1 level is high, even if the cholesterol ratio is low, the risk could be very high. Considering that the DDR1 level remains the same for one patient in the early stage of the aneurysm, thus controlling the cholesterol ratio can still minimize the risk of atherosclerotic aneurysm.

    Figure 5.  Risk map for atherosclerotic aneurysm over 2 years associated with the LDL/HDL-cholesterol ratio and the level of DDR1. The risk map is a contour plot of R(L/H,DDR1) with three regions divided by two contour curves 2 and 1.5: high-risk, medium-risk, and low-risk.

    We next test some potential intervention strategies for atherosclerotic aneurysm via our model, e.g. anti-cholesterol drugs, anti-DDR1 injection, and the combined therapy. Statins, only cholesterol-lowering drug class, are recommended for most patients with atherosclerosis to reduce the risk of heart attack or stroke [75]. In our study, we tested the efficacy of Statins on the atherosclerotic aneurysm. Although there are some medications for aneurysm patients whose size is less than 5 cm such as beta blockers and calcium channel blockers [76], these drugs are to lower your blood pressure and relax your blood vessels not related directly to the pathology of aneurysm [77]. Anti-DDR1 has been used in the mice experiments to slow the progression of the aneurysm [78,69]. Therefore, we will test the feasibility of the anti-DDR1 as a treatment for the atherosclerotic aneurysm. More specifically, we will test the three potential intervention strategies as follows:

    ● Anti-cholesterol drugs, statins, help lower LDL/HDL-cholesterol ratio in our model. We take the flux rate, α, in Eq (2.6) as a control function

    α(t)=106.5min(t/60,1)10,

    where the initial LDL/HDL-cholesterol ratio is assumed to be 10. One patient starts taking the medications to achieve the ideal cholesterol ratio, 3.5, by 2 months (60 days);

    ● The anti-DDR1 injection affects the flux rates, α, in Eqs (2.7) and (2.8) and is modeled as follows:

    α(t)=10(1et/100),

    where the DDR1 concentration is assumed to be 10 μm/ml initially and returns back to the initial value after one year since the injection;

    ● The combined therapy of the first two treatments is modeled by changing the flux rates, α, in Eqs (2.6)-(2.8) as the control functions mentioned above.

    Figure 6 shows the control functions for both the anti-cholesterol drugs and the anti-DDR1 injection in our study: anti-cholesterol drugs reduces the flux rate consistently after 2 months while the anti-DDR1 injection reduces the flux rate immediately and returns the original flux rate after around 1 year. To explore the three potential treatment strategies, we first let the atherosclerotic aneurysm develop to a certain size, e.g., running our model for 2 years, and then start three treatments for three years. The anit-DDR1 injection is given annually. Figure 7 shows the effects of three different treatments on both the aneurysm diameter and the plaque area. The anti-cholesterol drugs reduce the growth rate of plaque area from 7.19 cm2/year to 2.33 cm2/year and the growth rate of aneurysm diameter from 3.2 mm/year to 1.04 mm/year while the anti-DDR1 injection reduces the growth rates of plaque area and aneurysm diameter to 5.8 cm2/year and 2.56 mm/year, respectively. The combined therapy has the best efficacy by reducing the growth rates of plaque area and aneurysm diameter to 1.42 cm2/year and 0.64 mm/year, respectively. To further explore the efficacy of the combined therapy, we model the dosage of two drugs as follows

    α(t)=10γAmin(t/60,1)10 and α(t)=10(1γBet/100),
    Figure 6.  Flux rates α(t) are control functions for anti-cholesterol drugs (left) and the anti-DDR1 injection (right).
    Figure 7.  Three potential intervention strategies are given for three years after atherosclerotic aneurysm developed. Anti-cholesterol drugs are consistently given from year 2 to year 5 while the anti-DDR1 injection is given annually. Solid curves are for the control case (without treatments) while dash curves are for three different treatments.

    where γA and γB is the dosage of anti-cholesterol drugs and the anti-DDR1 injection, respectively. Then we define the efficacy effectiveness as

    E1(γA,γB)=d(γA,γB)d(0,0) and E2(γA,γB)=Area(γA,γB)Area(0,0),

    where d(γA,γB) and Area(γA,γB) are the diameter of aneurysm and the plaque area at Year 5 for given dosage γA and γB, respectively. Thus d(0,0) and Area(0,0) are the control values of no drug treatment. By varying γA from 0.3 to 0.9 with a step size 0.05 and γB from 0.1 to 1 with a stepsize 0.05, we have the efficacy maps for both E1 and E2 plotted in Figure 8 and show that the dosage of anti-cholesterol drugs affects significantly both plaque area and aneurysm diameter. Moreover, for any given dosage of anti-cholesterol drugs, the anti-DDR1 injection can also improve the efficacy of the combined therapy.

    Figure 8.  Efficacy maps of aneurysm diameter (left) and plaque area (right) for the combined therapy.

    More than 80% of aortic aneurysms are associated with atherosclerosis which develops firstly when cholesterol and fat build up inside the arteries. Atherosclerosis accelerates the breakdown of collagen and SMCs, structure and elasticity to the arterial wall. Over time, this causes the arterial walls to weaken and become damaged. As the diameter of the vessel grows, the wall tension increases, leading to even more dilation. The end result is an aneurysm. In this paper, we account for both cardiovascular diseases in the model and develop a comprehensive mathematical model to mimic the progression of the atherosclerotic aneurysm. Our model includes both the multi-layer structure of the arterial wall and the pathophysiology of the atherosclerotic aneurysm. More specifically, we model the plaque formed in the intima layer as a porous medium and use the linear elastic model to describe the media layer; we also include the role of macrophage, SMCs migration, and collagen proliferation into our hybrid model. Then our model is used to predict the progression of atherosclerotic aneurysm, the growth rate of the plaque and the aneurysm diameter. A 2D biomarker, the cholesterol ratio and the level of DDR1, is revealed to evaluate the risk of atherosclerotic aneurysm via our model shown in Figure 5. Then the efficacy of three clinical treatments is also illustrated via our model: the anti-cholesterol treatment is better than the anti-DDR1 injection but the combined therapy is optimal. We accordingly explore the efficacy of the combined therapy with different dosages of the anti-cholesterol treatment and the anti-DDR1 injection. Based on our results shown in Figure 8, we suggest that the anti-cholesterol treatment significantly decreases the risk of the atherosclerotic aneurysm while patients can also benefit from the additional anti-DDR1 injection.

    In the current study, we use a simplified circular geometry of the arterial wall with the radially symmetric assumption and do not consider the effect of complex arterial shapes. In the future study, we intend to incorporate the CT scan data of atherosclerotic aneurysm patients into our model and give a personalized prediction for each patient. This would also help us to evaluate the personalized efficacy of different treatments. Another limitation of the current work is the effect of blood pressure via the force balance condition in our model. However, we do not take the biological effects of the blood pressure into consideration, e.g. the inflammatory cell infiltration, therefore, might underestimate the effect of hypertension on the atherosclerotic aneurysm. In the future, we will include a more complex biological network to evaluate the effect of blood pressure. Moreover, the parameters in the current model will be calibrated when clinical data become available. Another future direction is to incorporate a better clinical serum biomarker than DDR1 which is not easily obtainable from a biomarker standpoint. One possible serum biomarker is platelet counts or activated platelet counts since platelets play a significant role in atherosclerosis and have been proved the cross-talk with other inflammatory cells. Furthermore, we plan to validate our suggested treatments to aneurysm patients and quantify the progression of atherosclerosis by using clinical data.

    GK and WH carried out all simulations and analyzed the data. GK, CH, GA, KO, MG and WH discussed the project and wrote the article.

    Data sharing is not applicable to this article as no new data were created or analyzed in this study.

    This work is supported by National Institutes of Health (Grant R01HL124155, CH) and by the American Heart Association (Grant 17SDG33660722, WH).

    The authors declare there is no conflict of interest.



    [1] K. Murphy, Immunobiology, Garland Science, 2012.
    [2] J. Gameiro, P. Nagib and L. Verinaud, The thymus microenvironment in regulating thymocyte differentiation, Cell Adhes. Migr., 4 (2010), 382–390.
    [3] B. Alberts, A. Johnson, J. Lewis, et al., Molecular Biology of the Cell, Garland Science, 2002.
    [4] A. Corthday, How do regulatory T cells work?, Scand. J. Immunol., 70 (2009), 326–336.
    [5] D. L. Farber, N. A. Yudanin and N. P. Restifo, Human memory T cells: Generation, compartmentalization and homeostasis, Nat. Rev. Immunol., 14 (2014), 24–35.
    [6] Y. Takahama, Journey through the thymus: Stromal guides for T-cell development and selection, Nat. Rev. Immunol., 6 (2006), 127–135.
    [7] C. H. Bassing, W. Swat and F. W. Alt, The mechanism and regulation of chromosomal V(D)J recombination, Cell, 109 (2002), S45–S55.
    [8] D. Mason, A very high level of crossreactivity is an essential feature of the T-cell receptor, Trends Immunol., 19 (1998), 395–404.
    [9] J. Nicolić-Žugić, M. K. Slifka and I. Messaoudi, The many important facets of T-cell repertoire diversity, Nat. Rev. Immunol., 4 (2004), 123–132.
    [10] D. J. Laydon, C. R. M. Bangham and B. Asquith, Estimating T-cell repertoire diversity: Limitations of classical estimators and a new approach, Philos. Trans. R. Soc. B, 370 (2015), 1–11.
    [11] M. S. Chaudhry, E. Velardi, J. A. Dudakov, et al., Thymus: The next (re)generation, Immunol. Rev., 271 (2016), 56–71.
    [12] G. G. Steinmann, B. Klaus and H. K. Müller-Hermelink, The involution of the ageing human thymic epithelium is independent of puberty, Scand. J. Immunol., 22 (1985), 563–575.
    [13] A. Globerson and R. B. Effros, Aging of lymphocytes and lymphocytes in the aged, Immunol. Today, 21 (2000), 515–521.
    [14] A. L. Gruver, L. L. Hudson and J. D. Sempowski, Immunosenescence of aging, J. Pathol., 211 (2007), 144–156.
    [15] D. P. Shanley, D. Aw, N. R. Manley, et al., An evolutionary perspective on the mechanisms of immunosenescence, Trends Immunol., 30 (2009), 374–381.
    [16] A. L. Gruver and G. D. Sempowski, Cytokines, leptin, and stress-induced thymic atrophy, J. Leukocyte Biol., 84 (2008), 915–923.
    [17] J. Dooley and A. Liston. Molecular control over thymic involution: From cytokines and microRNA to aging and adipose tissue, Eur. J. Immunol., 42 (2012), 1073–1079.
    [18] H. Selye, Thymus and the adrenals in the response of the organ to injuries and intoxications, Br. J. Exper. Pathol., 17 (1936), 234–248.
    [19] W. Savino, The thymus is a common target organ in infectious diseases, PLoS Pathog., 2 (2006), e62.
    [20] S. D. Wang, K. J. Huang, Y. S. Lin, et al., Sepsis-induced apoptosis of the thymocytes in mice, J. Immunol., 152 (1994), 5014–5021.
    [21] W. W Grody, S. Fliegiel and F. Naeim, Thymus involution in the acquired immunodeficiency syndrome, Am. J. Clin. Pathol., 84 (1985), 85–95.
    [22] W. Savino, M. Dardenne, L. A. Velloso, et al., The thymus is a common target in malnutrition and infection, Br. J. Nutr., 98 (2007), S11–S16.
    [23] C. L. Mackall, T. A. Fleischer, M. R. Brown, et al., Age, thymopoiesis, and CD4+ T-lymphocyte regeneration after intensive chemotherapy, New Engl. J. Med., 332 (1995), 143–149.
    [24] J. Storek, R. P. Witherspoon and R. Storb, T cell reconstitution after bone marrow transplantation into adult patients does not resemble T cell development in early life, Bone Marrow Transplant., 16 (1995), 413–425.
    [25] A. L. Zoller, F. J. Schnell and G. J. Kersh, Murine pregnancy leads to reduced proliferation of maternal thymocytes and decreased thymic emigration, Immunology, 121 (2007), 207–215.
    [26] T. A. Tibbetts, F. DeMayo, S. Rich, et al., Progesterone receptors in the thymus are required for thymic involution during pregnancy and for normal fertility, PNAS, 96 (1999), 12021–12026.
    [27] A. G. Rijhsinghani, K. Thompson and S. K. Bhatia, Estrogen blocks early T cell development in the thymus, Am. J. Reprod. Immunol., 36 (1996), 269–277.
    [28] J. D. Ashwell, F. W. M. Lu and M. S. Vacchio, Glucocorticoids in T cell development and function, Annu. Rev. Immunol., 18 (2000), 309–345.
    [29] D. A. da Cruz, J. S. Silva, V. C. de Almeida, et al., Altered thymocyte migration during experimental acute trypanosoma cruzi infection: Combined role of fibronectin and the chemokines CXCL12 and CCL4, Eur. J. Immunol., 36 (2006), 1486–1493.
    [30] S. K. Stanley, J. M. McCune, H. Kaneshima, et al., Human immunodeficiency virus infection of the human thymus and disruption of the thymic microenvironment in the SCID-hu mouse, J. Exper. Med., 178 (1993), 1151–1163.
    [31] M. G. Durdov, O. Springer, V. Ćapkun, et al., The grade of acute thymus involution in neonates correlates with the duration of acute illness and with the percentage of lymphocytes in peripheral blood smear, Biol. Neonate, 83 (2003), 229–234.
    [32] J. van Baarlen, H. J. Schuurman, R. Reitsma, et al., Acute thymuc involution during infancy and childhood: Immunohistology of the thymus and peripheral lymphoid tissues after acute illness, Pediat. Pathol., 19 (1989), 261–275.
    [33] E. Juretić, A. Juretić, B. Užarević, et al., Alterations in lymphocyte phenotype of preterm infected newborns, Biol. Neonate, 80 (2001), 223–227.
    [34] S. M. Falkenberg, C. Johnson, F.V. Bauermann, et al., Changes observed in the thymus and lymph nodes 14 days after exposure to BVDV field strains of enhanced or typical virulence in neonatal calves, Vet. Immunol. Immunopathol., 160 (2014), 70–80.
    [35] J. P. L. Rangel, C. S. Galan-Enriquez, M. López-Medina, et al., Bacterial clearance reverses a skewed T-cell repertoire induced by Salmonella infection, Immun., Inflammation Dis., 3 (2015), 209–223.
    [36] S. Yovino, L. Kleinberg, S. A. Grossman, et al., The etiology of treatment-related lymphopenia in patients with malignant gliomas: Modeling radiation dose to circulating lymphocytes explains clinical observations and suggests methods of modifying the impact of radiation on immune cells, Cancer Invest., 31 (2013), 140–144.
    [37] J. S. Mendez, A. Govindan, J. Leong, et al., Association between treatment-related lymphopenia and overall survival in elderly patients with newly diagnosed glioblastoma, J. Neuro-Oncol., 127 (2016), 329–335.
    [38] J. L. Campian, X. Ye, M.Brock, et al., Treatment-related lymphopenia in patients with stage Ⅲ non-small-cell lung cancer, Cancer Invest., 31 (2013), 183–188.
    [39] S. S. Long, Laboratory manifestations of infectious disease, in Principles and Practice of Pediatric Infectious Diseases, 4th edition, Elsevier, (2012), 1400–1412.
    [40] S. M. Ciupe, B. H. Devlin, M. L. Markert, et al., The dynamics of T-cell receptor repertoire diversity following thymus transplantation for digeorge anomaly, PLoS Comput. Biol., 5 (2009), e1000396.
    [41] J. F. Purton, J. A. Monk, D. R. Liddicoat, et al., Expression of the glucocorticoid receptor from the 1A promotor correlates with T lymphocyte sensitivity to glucocorticoid–induced cell death, J. Immunol., 173 (2004), 3816–3824.
    [42] F. K. Kong, C. L. H. Chen and M. D. Cooper, Reversible disruption of thymic function by steroid treatment, J. Immunol., 168 (2002), 6500–6505.
    [43] J. A. Guevara Patiño, M. W. Marino, V. N. Ivanov, et al., Sex steroids induce apoptosis of CD8+ CD4+ double positive thymocytes via TNF-α, Eur. J. Immunol., 30 (2000), 2586–2592.
    [44] P. L. Choyke, R. K. Zemon, J. E. Gootenberg, et al., Thymic atrophy and regrowth in response to chemotherapy: CT evaluation, Am. J. Roentgenol., 149 (1987), 269–272.
    [45] M. Cohen, C. A. Hill, A. Cangir, et al., Thymic rebound after treatment of childhood tumors, Am. J. Roentgenol., 135 (1980), 151–156.
    [46] D. E. DeFriend, J. M. Coote, M. P. Williams, et al., Thymic enlargement in an adult following a severe infection, Clin. Radiol., 56 (2001), 331–333.
    [47] D. W. Gelfand, A. S. Goldman, E. J. Law, et al., Thymic hyperplasia in children recovering from thermal burns, J. Trauma, 12 (1972), 813–817.
    [48] L. M. Bradley, L. Haynes and S. L. Swain, IL-7: Maintaining T-cell memory and achieving homeostasis, Trends Tmmunol., 26 (2005), 172–176.
    [49] J. T. Tan, E. Dudl, E. LeRoy, et al., IL-7 is critical for homeostatic proliferation and survival of naive T cells, Proc. Natl. Acad. Sci., 98 (2001), 8732–8737.
    [50] L. Vivien, C. Benoist and D. Mathis, T lymphocytes need IL-7 but not IL-4 or IL-6 to survive in vivo, Int. Immunol., 13 (2001), 763–768.
    [51] T. J. Fry and C. L. Mackall, The many faces of IL-7: From lymphopoesis to peripheral T cell maintenance, J. Immunol., 174 (2005), 6571–6576.
    [52] S. Xu and T. Chou, Immigration-induced phase transition in a regulated multispecies birth-death process, J. Phys. A: Math. Theor., 51 (2018), 425602.
    [53] J. Hataye, J. J. Moon, A. Khoruts, et al., Naïve and memory CD4+ T cell survival controlled by clonal abundance, Science, 312 (2006), 114–116.
    [54] R. Dessalles, M. D'Orsogna and T. Chou, Exact steady-state distributions of multispecies birth–death–immigration processes: Effects of mutations and carrying capacity on diversity, J. Stat. Phys., 173 (2018), 182–221.
    [55] L. Westera, V. van Hoeven, J. Drylewicz, et al., Lymphocyte maintenance during healthy aging requires no substantial alterations in cellular turnover, Aging Cell, 14 (2015), 219–227.
    [56] H. S. Robins, P. V. Campregher, S. K. Srivastava, et al., Comprehensive assessment of T-cell receptor β-chain diversity in αβ T cells, Blood, 114 (2009), 4099–4107.
    [57] G. Lythe, R. E. Callard, R. L. Hoare, et al., How many TCR clonotypes does a body maintain?, J. Theor. Biol., 389 (2016), 214–224.
    [58] S. Lewkiewicz, Y. L. Chuang and T. Chou, A mathematical model of the effects of aging on naive T cell populations and diversity, Bull. Math. Biol., 81 (2019), 2783–2817.
    [59] N. Vrisekoop, I. den Braber, A. B. de Boer, et al., Sparse production but preferential incorporation of recently produced naive T-cells in the human peripheral pool, Proc. Natl. Acad. Sci., 105 (2008), 6115–6120.
    [60] R. J. de Boer and A. S. Perelson, Quantifying T lymphocyte turnover, J. Theor. Biol., 327 (2013), 45–87.
  • This article has been cited by:

    1. Shankar Narayan S, Sunanda Saha, Time-dependent study of blood flow in an aneurysmic stenosed coronary artery with inelastic walls, 2021, 47, 22147853, 4718, 10.1016/j.matpr.2021.05.608
    2. Yatao Liu, Mathematical modeling and simulation of atherosclerotic formation and progress at arterial bifurcations, 2023, 80, 17468094, 104266, 10.1016/j.bspc.2022.104266
    3. Ali A. Rostam-Alilou, Hamid R. Jarrah, Ali Zolfagharian, Mahdi Bodaghi, Fluid–structure interaction (FSI) simulation for studying the impact of atherosclerosis on hemodynamics, arterial tissue remodeling, and initiation risk of intracranial aneurysms, 2022, 21, 1617-7959, 1393, 10.1007/s10237-022-01597-y
    4. Shankar Narayan S., Anuradha Bhattacharjee, Sunanda Saha, CFD analysis of the hyper-viscous effects on blood flow across abdominal aortic aneurysm in COVID patients: multiphysics approach, 2023, 1025-5842, 1, 10.1080/10255842.2023.2194474
    5. Egle Kavaliunaite, Thomas Emil Andersen, Jes S. Lindholt, Jane Stubbe, Daily fenugreek intake does not attenuate abdominal aortic aneurysm growth in rats, 2025, 0301-1526, 10.1024/0301-1526/a001185
  • Reader Comments
  • © 2020 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(4398) PDF downloads(480) Cited by(2)

Figures and Tables

Figures(7)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog