Next Article in Journal
Experimental Morogoro Virus Infection in Its Natural Host, Mastomys natalensis
Next Article in Special Issue
Will SARS-CoV-2 Become Just Another Seasonal Coronavirus?
Previous Article in Journal
Convalescent Plasma Transfusion for the Treatment of COVID-19 in Adults: A Global Perspective
Previous Article in Special Issue
Endogenously Produced SARS-CoV-2 Specific IgG Antibodies May Have a Limited Impact on Clearing Nasal Shedding of Virus during Primary Infection in Humans
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling the Effect of MUC1 on Influenza Virus Infection Kinetics and Macrophage Dynamics

1
School of Mathematics and Statistics, The University of Melbourne, Parkville, VIC 3010, Australia
2
Peter Doherty Institute for Infection and Immunity, The Royal Melbourne Hospital and The University of Melbourne, Parkville, VIC 3010, Australia
3
Melbourne School of Population and Global Health, The University of Melbourne, Parkville, VIC 3010, Australia
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(5), 850; https://doi.org/10.3390/v13050850
Submission received: 25 March 2021 / Revised: 20 April 2021 / Accepted: 4 May 2021 / Published: 7 May 2021
(This article belongs to the Special Issue Mathematical Modeling of Viral Infection)

Abstract

:
MUC1 belongs to the family of cell surface (cs-) mucins. Experimental evidence indicates that its presence reduces in vivo influenza viral infection severity. However, the mechanisms by which MUC1 influences viral dynamics and the host immune response are not yet well understood, limiting our ability to predict the efficacy of potential treatments that target MUC1. To address this limitation, we use available in vivo kinetic data for both virus and macrophage populations in wildtype and MUC1 knockout mice. We apply two mathematical models of within-host influenza dynamics to this data. The models differ in how they categorise the mechanisms of viral control. Both models provide evidence that MUC1 reduces the susceptibility of epithelial cells to influenza virus and regulates macrophage recruitment. Furthermore, we predict and compare some key infection-related quantities between the two mice groups. We find that MUC1 significantly reduces the basic reproduction number of viral replication as well as the number of cumulative macrophages but has little impact on the cumulative viral load. Our analyses suggest that the viral replication rate in the early stages of infection influences the kinetics of the host immune response, with consequences for infection outcomes, such as severity. We also show that MUC1 plays a strong anti-inflammatory role in the regulation of the host immune response. This study improves our understanding of the dynamic role of MUC1 against influenza infection and may support the development of novel antiviral treatments and immunomodulators that target MUC1.

1. Introduction

Influenza is a contagious respiratory disease. It remains a major public health burden that affects and threatens millions of people each year [1]. Influenza virus (IV) primarily attacks the epithelial cells that line the upper respiratory tract (URT) of the host, causing an acute infection [2]. The host immune response has been shown to play an important role against influenza infection [3,4]. As part of the innate immune response, macrophages that reside in airways limit viral dissemination through phagocytosis of viral particles and prevent the virus from spreading to the lungs [5,6]. Activated macrophages produce inflammatory molecules, such as TNF- α , which stimulates recruitment of additional immune cells, such as monocyte-derived macrophages (MDMs) to the site of infection. These molecules also facilitate the activation of adaptive immune responses, such as maturation of B cells and effector CD8 + T cells [7]. Thus, macrophages play a critical role against influenza viral infection [8,9,10].
However, recruited macrophages also amplify inflammation. Overstimulation of the host immune response can lead to pathology, indicating that there is a subtle balance between a protective and a destructive response [1,11]. A dysregulated immune response, often marked by an excessive recruitment of macrophages to the site of infection and a high level of cytokine production, can lead to lung pathology, causing serious and sometimes fatal infection outcomes [12,13,14,15].
MUC1 belongs to the family of cell surface (cs-) mucin and is constitutively expressed at the surface of respiratory epithelial cells and macrophages, as reviewed in [16,17,18]. It was shown to be capable of modulating cytokine production in vitro viral infection [19,20,21] and in vivo bacterial infection [22,23]. More recently, McAuley and colleagues investigated the in vivo effects of MUC1 on influenza viral infection [24]. They first intranasally infected wildtype (WT) and MUC1-knockout (KO) mice with influenza A virus, then measured and compared the kinetic time-series data of viral load as well as different immune cells between the two groups. They found that the virus grows more quickly and reaches a peak earlier in MUC1-KO mice. Mice displayed a more enhanced inflammatory response, dominated by a higher number of macrophages and a high level of cytokine production. Based on these observations, they hypothesized that MUC1 acts as physical barrier to prevent virus from infecting epithelial cells and contribute to regulation of the host immune response. However, the potential effects of MUC1 in vivo are poorly quantified, limiting our ability to predict the efficacy of potential treatments that target MUC1. To address this limitation, we incorporated the hypothesized effects of MUC1 into mathematical models of influenza viral dynamics and applied Bayesian inference to estimate key parameter values and provided new quantitative insight into the role of MUC1 in shaping influenza virus infection and the host immune response.
Influenza viral dynamics models have been used to study many aspects of influenza infection and the host immune response, as reviewed in [25]. Studies focusing on the immune system have used viral dynamics models to study various types of immunological data, sharpening to our understanding of the contribution of different immunological components to influenza viral infection [26,27,28].
In this work, we use available in vivo kinetic data for both virus and macrophage populations in wildtype and MUC1 knockout mice. We analyze the data with two mathematical models of influenza viral dynamics under a Bayesian framework, quantifying the potential effects of cs-mucin MUC1 in influenza infection. The two models differ in how they categorise mechanisms of viral control. We also use the data-calibrated models to evaluate and analyze the dependence of various infection-related quantities on MUC1 expression. Finally, we discuss the biological implications of our results.

2. Results

2.1. Model Fitting

In vivo viral load and macrophage data in WT and MUC1-KO mice were used in model fitting. We fitted a Target cell-Infected-cell-Virus (TIV) model (Equations (1)–(4) in Materials and Methods) and an Immune Response (IR) model (Equations (5)–(18)) to the data, respectively. MUC1 has been suggested to prevent virus from infecting epithelial cells. It also has been implicated in the regulation of the host innate immune response, associated with macrophage recruitment [24]. As detailed in Materials and Methods, both models capture these effects. The reduction in susceptibility of target cells is captured by a parameter ε 1 , modulating viral infectivity to the target cells in d T / d t = ( 1 ε 1 ) β T V (Equation (1)). The effect of the limitation of macrophage recruitment is captured by a parameter ε 2 and is modelled in d M / d t = s + ( 1 ε 2 ) ϕ I δ M M (Equation (4)). In the absence of MUC1 expression, e.g., in MUC1-KO mice, we set ε 1 = ε 2 = 0 to represent a complete knockout effect.
The fitting results are shown in Figure 1. The median of the posterior prediction (solid line) and a 95% predict interval (PI, shaded area) were computed from 4000 model simulation based on 4000 samples from the posterior distribution of model parameters (provided in Supplementary Figures). The trend for both the viral kinetics (Figure 1A,B) and macrophages dynamics (Figure 1C,D) is well captured by the median prediction in both models, suggesting that both models are able to explain the data. Moreover, the narrow 95% PI indicates a relatively high certainty level for model predictions.

2.2. Estimates of MUC1 Parameters

The marginal posterior densities for ε 1 and ε 2 provide insight into the role of MUC1. The median parameter estimates and their associated 95% credible intervals (CIs) are given in Table 1. The median estimate of MUC1 on reduction of viral infectivity ( ε 1 ) is 0.44 (95% CI: 0.23–0.71) in the TIV model and 0.42 (95% CI: 0.22–0.58) in the IR model. Furthermore, the estimated median values of MUC1 on regulation of macrophage recruitment ( ε 2 ) are 0.45 (95% CI: 0.18–0.64) and 0.38 (95% CI: 0.06–0.63) in the TIV and IR models, respectively. Biologically, the median estimate for ε 1 indicates that the presence of MUC1 reduces the rate of virus infection to epithelial cells by 44% (for the TIV model) or 42% (for the IR model). In addition, the median estimate for ε 2 indicates that the presence of MUC1 reduces the recruitment rate of macrophages induced by infected cells by 45% ( for the TIV model) or 38% (for the IR model).
The posterior-median estimates are qualitatively consistent between the two models. The median estimates of ε 1 and ε 2 both exclude 0 (in the 95% credible interval) for the WT group (for both the TIV and IR models), indicating a reduced viral infectivity (i.e., ( 1 ε 1 ) β ) and a reduced rate of macrophage recruitment induced by infected cells (i.e., ( 1 ε 2 ) ϕ ) in the presence of MUC1. The results support the experimental hypothesis [24] and provide quantitative evidence that the presence of MUC1 reduces viral infectivity to epithelial cells. They also provide evidence that MUC1 reduces macrophage recruitment and thus regulates the host innate immune response.
Detailed posteriors of model parameters are provided in Supplementary Figures S1–S10, and correlation maps of the estimated parameters for the TIV and IR models are given in Figures S11 and S12. There is a low correlation coefficient between ε 1 and ε 2 for the TIV ( R = 0.08 ) and IR ( R = 0.22 ) models, suggesting the two parameters have a weak relationship. In particular, we found that the posterior-median estimate of the phagocytosis rate of virus by macrophages ( κ M ) is approximately 10 8 for the TIV model, and the estimate is in agreement with the estimate for the IR model (Figure S9). We used the median estimates of model parameters to compute the ratio of macrophage-mediated viral decay ( κ M M ( t ) ) to overall viral decay rate in the TIV ( κ M M ( t ) + δ V ) and IR ( κ M M ( t ) + δ V + κ A S A S ( t ) + κ A L A L ( t ) ) models as a time-series, respectively. We found that κ M M ( t ) only has a minor contribution to viral clearance (Figure S14). The result suggests that macrophages, although important to maintain gas exchange in lungs and reduce infection severity, are not directly involved in limiting viral replication, as shown in [29,30].

2.3. Prediction of Infection-Related Quantities

Influenza pathogenesis is often associated with a high viral load and an overstimulated immune response [15]. In the absence of MUC1, mice showed a significantly high mortality rate [24]. Here, we use the 4000 joint posterior distributions to predict the impact of MUC1 on some key infection-related quantities that likely influence infection severity. We then compare these quantities between the two models.
The basic reproduction number of viral replication ( R 0 ) is defined as the average number of secondary infected cells that are produced by an initially infected cell when the target cell population is not depleted and is fully susceptible [31]. An infection may be established only if R 0 > 1 . It is a critical indicator that quantifies how fast an infection is established and evolved.
Figure 2A,B show the R 0 between WT and MUC1-KO groups in the TIV and IR models, respectively. Both models predict a significantly higher median value of R 0 (dashed line) in the MUC1-KO group (20 in MUC1-KO group versus 11.1 in WT group for the TIV model, and 45.6 versus 26.4 for the IR model). The estimates of R 0 are comparable to previous estimates from fitting viral dynamics models to viral kinetic data in humans [32] and mice [33].
To assess the impact of MUC1 on viral dynamics, we compute the area under the viral load (without log-transformation) curve, which is often used as a marker for infectiousness (shown in Equation (21) in Materials and Methods). Both the TIV (Figure 2C) and IR (Figure 2D) models predict very similar log 10 (AUC V ) in WT and MUC1-KO mice. This implies that a paucity of MUC1 expression has little, if any, effect on the cumulative viral load. This observation is supported by data in [24] in which they found that MUC1-KO mice were still capable of clearing virus after day 7 post infection.
An excessive accumulation of macrophages is considered to be a hallmark for severe infection, often observed in highly pathogenic influenza viral infection [14]. We use the area under the macrophage time-series curve (without log-transformation; Equation (22) in Materials and Methods) as a surrogate for the strength of immune response stimulation and explore the dependence of the AUC M on MUC1. As shown in Figure 2E,F, both models predict a higher median value of log 10 (AUC M ) in MUC1-KO mice compared to WT mice. This suggests that MUC1 reduces the accumulation of macrophages and thus contributes to the regulation of the host immune response.
We also assessed the influence of MUC1 on peak viral load (Figure S13A,B in Supplementary Figures) and peak viral load time (Figure S13C,D) for the two models. Both models predict that the presence of MUC1 delays the time at which viral load peaks but only has a subtle influence on the magnitude of peak viral load, as shown in [24].
In summary, both models predict a higher value of R 0 (Figure 2A,B) and increased macrophage accumulation (Figure 2E,F) in the absence of MUC1 expression. The results emphasise the dual roles for MUC1 in reducing viral infectivity and limiting macrophage recruitment. Furthermore, they suggest that the absence of MUC1, while not driving a significant increase in cumulative viral load, facilitates viral replication and dissemination within the host in the early stages of infection. More epithelial cells are infected in a short time interval, heightening macrophage recruitment, likely contributing to lung pathology and providing an explanation for the heightened mortality rate in MUC1 KO mice.

2.4. Delineation the Effects of MUC1 on Macrophage Recruitment

We showed that the presence of MUC1 reduces AUC M (Figure 2E,F), which may alleviate infection severity. The accumulation of macrophages is not only directly impacted by the regulatory effect of MUC1, (i.e., ε 2 ), but is also indirectly affected by antigen levels, which are influenced by ε 1 through modulating dynamics for infected cells (I). Here, we analyze the relative contribution of the two effects of MUC1 on the AUC M . We use the macrophage reduction efficiency, defined as the decrease in the area under the macrophage curve in wild type mice (AUC M , W T ) relative to the AUC of the macrophage curve in MUC1 knockout mice (AUC M , K O ):
Macrophage Reduction Efficiency = 1 A U C M , W T A U C M , K O
Figure 3 shows the estimated marginal posterior density of ε 1 and ε 2 for the TIV model (top panel) and a heatmap of the dependence of macrophage reduction efficiency on ε 1 and ε 2 (bottom panel). The heatmap predicts the dependence of the macrophage reduction efficiency for various values of ε 1 and ε 2 within the 95% CI. We observe that a higher ε 1 or ε 2 leads to a higher macrophage reduction level, suggesting that both effects contribute to reduce the accumulation of macrophages. However, the macrophage reduction efficiency is notably more sensitive to changes in ε 2 . In particular, taking the median parameter values as a reasonable prediction point (black circle), the rate of change in the macrophage reduction efficiency is strongly dependent on ε 2 and only weakly dependent on ε 1 . (indicated by the arrow line). The result suggests that the reduction in macrophage recruitment (i.e., via an increase in ε 2 ) is more influential than the reduction of infected cells (i.e., via a reduction in ε 1 ) in determining AUC M . We also assess the macrophage reduction efficiency as a function of ε 1 and ε 2 for the IR model. As shown in Figure 4, the results are qualitatively consistent—the macrophage reduction efficiency is strongly influenced by ε 2 .
Both models predict a strong effect for ε 2 and relatively small effect for ε 1 on the AUC M . This is understood by recalling that the presence of MUC1 does not significantly influence the cumulative viral load, as shown in Figure 2C,D. Thus, a change in the reduction of viral infectivity to target cells ( ε 1 ) has only a minor effect on the AUC M . The results emphasise a strong regulatory effect of MUC1 on macrophage accumulation.

3. Discussion

In this work, we studied the in vivo immunological effects of cs-mucin MUC1 in influenza viral infection. To the best of our knowledge, this is the first study to incorporate the dynamical roles of MUC1 into models of influenza virus dynamics. In our models, MUC1 regulates both viral replication and macrophage recruitment. We incorporated the experimentally hypothesized roles of MUC1 into two mathematical models and fitted kinetic data of both virus and macrophage populations to the models in a Bayesian framework. Our estimation results (Table 1) provide evidence that MUC1 reduces the susceptibility of epithelial cells to viral infection. They also provide evidence that MUC1 limits the recruitment of macrophages and thus regulates the host immune response. Both models predict the influence of MUC1 on various infection-related quantities (Figure 2). While the expression of MUC1 has little impact on the cumulative viral load (AUC V ), it delays viral infection by reducing the basic reproduction number of viral replication (Figure 2A,B) and delaying viral load peak time (Figure S11C,D). More importantly, we found that the presence of MUC1 significantly reduces the accumulation of macrophages (Figure 2E,F). The decreased level of macrophages is primarily driven by the direct regulatory effect ( ε 2 ) of MUC1 on macrophage recruitment (Figure 3 and Figure 4).
Our model-based analyses provide new insight into the mechanisms by which MUC1 influences viral dynamics and the host immune response. This is also the first study that we are aware of that provides quantitative estimates of the in vivo effects of cs-mucin MUC1 on influenza infection. Our analyses enhance our ability to predict the efficacy of potential treatments that target MUC1. Influenza pathogenesis is often marked by a high viral load, and infection of epithelial cells is a key determinant of the level of viral load [11,15,34]. MUC1 is rapidly stimulated at the surface of epithelial cells and macrophages on infection, and is thought to act as a “releasable decoy”, preventing virus from attaching and infecting the cells, thereby reducing viral infectivity [17]. Regardless of the specific mechanism, our model predictions suggest that MUC1 only effectively reduce R 0 (Figure 2B) but not the AUC V (Figure 2C,D). The biological implications of this are two-fold. Firstly, MUC1, as part of the innate immune response, has been shown to be rapidly up-regulated within a few hours post in vitro infection [21]. The decreased R 0 suggests that MUC1 expression contributes to limit and delay viral infection, and more importantly, to prevent viral dissemination within the host. This provides strong protection to the host and reduce infection severity. Viral spread to the lower respiratory tract (LRT) is known to cause complications, leading to more severe infection outcomes [34]. Secondly, the comparable AUC V between WT and MUC1-KO group implies that a lack of cs-mucin MUC1 protection have a subtle influence on other immunological components that are responsible for viral clearance, such as the host adaptive immune response. This may be partially supported in [24], where MUC1-KO mice were shown to clear virus from the lungs at day 7 post infection. A more comprehensive dataset that captures the dynamics of antibodies or effector CD8 + T cells would greatly improve our understanding of the impact on MUC1 to the adaptive immune response.
Beyond these virological indicators, viral pathogenesis is also associated with the strength of the host immune response induced by influenza infection. An excessive recruitment of macrophages to the site of infection is a hallmark of overstimulated immune responses [14,34]. The anti-inflammatory role of MUC1 has been shown to inhibit activation of Toll-like receptors (TLRs) in macrophages and infected cells [24]. In MUC1-KO mice, our models predicted a significantly enhanced level of AUC M (Figure 2E,F), which may reflect the high mortality rate in the group. This finding emphasises the importance of quantities related to the immune response, which can be critical indicators for predicting the severity of infection and facilitating the assessment of antiviral therapies, as suggested in [34,35]. Furthermore, we showed that the decreased AUC M is primarily due to the direct regulatory effect of MUC1 on macrophages (i.e., ε 2 ), which highlights a strong anti-inflammatory effect for MUC1. This may support the development of novel immunomodulators that target cs-mucin MUC1.
In conducting this study, we applied two mathematical models to the kinetic data. The models models differ in how they model adaptive immunity. We compared the key estimation results of MUC1’s effects and model predictions of infection-related quantities between the two models. We found that both models fit the in vivo viral load and macrophage data well (Figure 1), giving comparable parameter estimates and consistent biological insights.
One of the most important applications of viral dynamic models is to estimate key kinetic parameters, as reviewed in [25]. Model selection for data fitting is an important but unresolved challenge in influenza dynamics modelling due to limited time-series data on numerous quantities of interested. Parameter estimates vary substantially between different studies, and the predictive power of any given model is influenced by the selection of model components, as showed in previous work by our group [27,35] and others [36]. In our study, there are advantages and disadvantages in applying the TIV and IR models. Due to its simple model structure, the TIV model is more computationally efficient. However, its lack of a detailed characterization of adaptive immunity makes the model difficult to use to explore potential interactions between different immunological components, e.g., interactions between macrophages and CD8 + T cells. The IR model, on the other hand, is more computationally intensive and has far more parameters to either estimate or determine from the literature. However, it is more suitable for explaining in vivo kinetic viral load data to which adaptive immunity has been shown to have an influence. It also provides a platform to study more complicated virus-immunity dynamics and interaction between different components of immune responses.
Neither the TIV nor IR models consider the full spectrum of host immune response which are known to contribute to viral control and that have been included in other modelling works, e.g., interferon dynamics [27,28]. Regardless, we argue our two models are sufficient for this study in which we focus on the influence of MUC1 on viral dynamics and macrophage kinetics, which are both explicitly considered in the models. Furthermore, there is no evidence to suggest that MUC1 has an impact on the adaptive immune response. Combined with the observation that MUC1-KO mice clear virus after day 7 post infection [24], the effects of MUC1 may be minimally influenced by the detailed dynamics of adaptive immunity.
Our study has some limitations. We only incorporated the two hypothesized effects of cs-mucin MUC1 on influenza viral infection into our mathematical models, but did not consider the detailed dynamics of MUC1 itself due to a lack of MUC1 kinetic data. As a result, the critical timing at which MUC1 starts to take effect has not been estimated. This could be an important factor that influences disease severity [17]. In future work, explicitly modelling the time dependent MUC1 effects would be of interest given availability of time-series data of MUC1 expression. Another limitation is that we assumed a fixed adaptive immune response, such that the adaptive immune responses dominate viral clearance at day 5 post infection regardless of MUC1 expression [27,37]. Though there is no evidence so far that MUC1 would affect the magnitude and/or timing of the adaptive immune response, extension of the IR model to allow for such an effect may be of interest.

4. Materials and Methods

4.1. Mathematical Models

In this study, we considered two mathematical models that are often used to study within-host influenza dynamics, but which differ in how they categorise the mechanisms of viral control.

4.1.1. The TIV Model

The Target cell-Infected cell-Virus (TIV) model depicts a simple but fundamental interaction between target cells and influenza virus, as originally presented in [32]. To estimate the in vivo impacts of MUC1, we incorporate the two hypothesized effects of MUC1 on viral infectivity and innate immune responses into the TIV model. We also consider a component of macrophage dynamics and critical interactions between macrophages and virus. The model is described by a set of ordinary differential equations (ODEs):
d T d t = g T 1 T + I T m a x ( 1 ε 1 ) β T V ,
d I d t = ( 1 ε 1 ) β T V δ I I ,
d V d t = p I δ V V κ M M V ,
d M d t = s + ( 1 ε 2 ) ϕ I δ M M .
Equations (1)–(3) describe the interaction between virus and epithelial cells. In detail, epithelial cells (T), the target cells for influenza virus, are infected with virus (V) and become infected cells (I) at an infectivity rate β V per day. Target cells are replenished at a rate g T ( 1 ( T + I ) / T m a x ) , where T m a x is the maximal number of epithelial cells that line the upper respiratory tract (URT). The infectivity rate is modified by MUC1, parameterised by ε 1 . Infected cells produce free virus at a rate p per day. Apoptosis occurs at a rate δ I per day. We do not explicitly model the role of macrophages in removing apoptotic infected cells [3,6,38]. While an established role of macrophages, it is not required as we have no data on dead cell dynamics and so our model does not include the dynamics of the dead cell population. The decrease of free virus is either due to natural decay at a constant rate δ V per day, or internalization by macrophages (M) at a rate κ M M .
Equation (4) models the dynamics of macrophages. We assume a constant supplementary rate and a decay rate of macrophages at s and δ M per day, respectively. Upon infection, monocytes are recruited from peripheral blood to the site of infection and become monocyte-derived macrophages (MDMs) in the presence of cytokines. We assume the recruitment rate is proportional to the level of infected cells, ϕ I , as infected cells contribute to cytokines production. The cs-mucin MUC1 regulates the recruitment rate of macrophages, parameterised by ε 2 .

4.1.2. The IR Model

The immune response (IR) model is based on the TIV model and includes a detailed adaptive immune response, which contributes to viral clearance over a distinct timescale [28]. The model is formulated by a system of ODEs:
d T d t = g T 1 T + I T m a x ( 1 ε 1 ) β T V ,
d I d t = ( 1 ε 1 ) β T V δ I I κ E E I ,
d V d t = p I δ V V κ M M V κ A S A S V κ A L A L V ,
d M d t = s + ( 1 ε 2 ) ϕ I δ M M ,
d E 0 d t = γ E V V + E 50 E 0 ,
d E 1 d t = γ E V V + E 50 E 0 n E τ E E 1 ,
d E i d t = n E τ E ( E i 1 E i ) , i = 2 , , n E
d E d t = ϕ E n E τ E E n E δ E E ,
d B 0 d t = γ B V V + B 50 B 0 ,
d B 1 d t = γ B V V + B 50 B 0 n B τ B B 1 ,
d B i d t = n B τ B ( B i 1 B i ) , i = 2 , , n B
d P d t = ϕ p n B τ B B n B δ p P ,
d A S d t = μ S P δ A S A S ,
d A L d t = μ L P δ A L A L .
Equations (5)–(8) retain the skeleton of the TIV model, describing the essential target cell-virus dynamics, except for additional components in d I / d t and d V / d t related to adaptive immune responses. κ E E in Equation (6) represents the rate of infected cells lysis by effector CD8 + T cells. The extra terms κ A S A S and κ L S A L in Equation (7) represent virus clearance mediated by a short-lived ( A S , e.g., IgM) and a long-lasting antibody ( A L , e.g., IgG), respectively.
Equations (9)–(12) describe a major component of the cellular adaptive immune response mediated by CD8 + T cells. Naïve CD8 + T cells ( E 0 ) initiate proliferation and differentiate into effector cells E 1 on stimulation via antigen-presentation at a rate γ E V / ( V + E 50 ) , where γ E is the maximal stimulation rate, and E 50 is a half saturation level at which half of the stimulation rate is obtained (as shown in Equation (9)). Effector cells E 1 perform programmed proliferation to E i where i denotes proliferation stages (Equations (10) and (11)) for τ E days, experience through n E stages [39], finally becoming mature effector cytotoxic T lymphocytes (E) at a rate ϕ E at the final stage. The decay rate of E is δ E , as shown in Equation (12).
Similarly, the dynamics of the humoral adaptive immune response are described by Equations (13)–(16). Naïve B cells ( B 0 ) start to proliferate and differentiate into plasma cells ( B 1 ) once stimulated by virus at a rate γ B V / ( V + B 50 ) , where γ B is the maximal stimulation rate and B 50 is a half-saturation level, as shown in Equation (13). Equations (14) and (15) capture how plasma cells ( B 1 ) undergo programmed proliferation through n B stages into B i , where i denotes proliferation stages, for τ B days [39]. Finally, mature plasma cells P (Equation (16)) are produced at a rate ϕ B and decay at a rate δ p .
Equations (17) and (18) describe the dynamics of a short-lived antibody ( A S ) and a long-lived antibody ( A L ). A S and A L are produced by plasma cells (P) at rates μ S and μ L and decay at rates δ A S and δ A L , respectively.

4.2. Statistical Inference

We extracted the kinetic data of both virus and macrophage population in wild type (WT) and MUC1 knockout mice using WebPlotDigitizer (version 4.4) from [24]. In the study, groups of wild type and MUC1-KO mice were intranasally infected with influenza A virus (PR8). There were five mice in each group. We assumed the variability of virus and macrophage data between different mice within the same group was due to measurement error, so that the data from different mice were pooled together for analysis.
We took a Bayesian inference approach to fit the TIV and IR model (detailed in Model) to the log-transformed kinetic data. In detail, our model has 10 parameters to estimate, and the parameter space is denoted as Φ = ( ε 1 , β , δ I , p , δ V , s , δ M , ε 2 , κ M , ϕ ) . Upon calibrating the IR model, we fixed all parameters of the adaptive immune responses (e.g., all parameters in Equations (9)–(18)) to previous estimated values in the literature [27,39]. We fixed the parameters because estimating the immunological effects of adaptive immunity is not a focus of this study, and [24] does not provide sufficient data for estimation of these parameters. We chose the number of effector T cell and B cell division cycle (i.e., n E and n B ) and the total proliferation time of the cells (i.e., τ E and τ B ) from [39], such that the adaptive immune responses only become activated five days post infection. The fixed parameter values are given in Table S1 in Supplementary Materials.
Furthermore, we assumed WT and MUC1-KO mice only differ in ε 1 and ε 2 , a reasonable assumption given inbred mice and use of the same virus for all experiment. We fitted log-transformed WT and MUC1-KO data simultaneously to the models with the same parameter vector set, only differing except for ε 1 and ε 2 , which were set to ε 1 = ε 2 = 0 for MUC1-KO mice. The prior distribution for model parameters ( Φ ) is given in Table S2 in Supplementary Materials. The distribution of the observed log-transformed viral load and macrophage measurement is assumed to be a normal distribution with a mean value given by the model simulation results and standard deviation (SD) parameter with prior distribution of a normal distribution with a mean of 0 and a SD of 1.
Model fitting was performed in R (version 4.0.2) and Stan (Rstan 2.21.0). Hamiltonian Monte Carlo (HMC) optimized by the No-U-Turn Sampler (NUTS) [40] was implemented to draw samples from the joint posterior distribution of the model parameters. A detailed theoretical foundation of HMC can be found in [41]. In particular, we used four chains with different starting points and ran 2000 iterations for each chain, discarding the first 1000 iterations as burn-in. We retained 4000 samples in total from 4 chains (1000 for each) after the burn-in iterations. The marginal posterior and prior density for all parameters are shown in Supplementary Materials. We calculated the median and quantiles of 2.5% and 97.5% of the 4000 model outputs at each time for posterior prediction and a 95% prediction interval (PI), respectively (e.g., Figure 2).

4.3. Infection-Related Quantities

The basic reproduction number of viral replication ( R 0 ) is given by
R 0 = ( 1 ε 1 ) β T 0 V ( δ I + κ E E ( 0 ) ) ( δ V + κ M M 0 + κ A S A S ( 0 ) + κ A L A L ( 0 ) ) ,
where T 0 is the initial number of epithelial cells, and M 0 is the number of macrophages in a disease-free equilibrium, given by s / δ M . E ( 0 ) , A S ( 0 ) and A L ( 0 ) are initial values of effector CD8 + T cells, a short-term antibody and a long-term antibody, respectively, which are set to zero at the beginning of infection as mice are naive (i.e., not previously exposed to influenza). Therefore, R 0 simplifies to
R 0 = ( 1 ε 1 ) β T 0 V δ I ( δ V + κ M M 0 ) .
Please note that ε 1 = 0 in MUC1-KO group. The area under the viral load time-series curve (AUC V ) and under the macrophage time-series curve (AUC M ) are given by
A U C V = 0 τ V ( t ) d t ,
A U C M = 0 τ M ( t ) d t ,
where τ is a cut-off day for calculation. We set τ = 14 , which covers the duration of viral infection, macrophage dynamics and clinical dynamics in [24]. V ( t ) and M ( t ) are simulated time series of viral load and macrophages, respectively.
The estimates of the infection-related quantities were computed using the 4000 posterior samples by solving the ode solver ode15s in MATLAB R2019b with a relative tolerance of 1 × 10 5 and an absolute tolerance of 1 × 10 10 . The initial values for different model components in the TIV model is ( T , I , V , M ) = ( 1 × 10 7 , 0 , 30 , s / δ M ) , where s and δ M are estimated from fitting the macrophage data to the model. For the IR model, the initial values were ( T , I , V , M , E 0 , E 1 E , B 0 , B 1 P , A S , A L ) = ( 1 × 10 7 , 0 , 30 , s / δ M , 100 , 0 , 0 , 100 , 0 , 0 , 0 , 0 ) . The values of fixed parameters are given in Supplementary Materials (Table S2). All visualization was performed in R (version 4.0.2). Computer codes to produce all the figures in this study can be found at https://github.com/keli5734/MUC1 (accessed on 5 May 2021).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13050850/s1, Figures S1–S10: Posterior distribution for estimated model parameters Φ . Figure S11–S12: Correlation map of estimated parameters in the TIV and IR models. Figure S13: Comparison of model predictions for peak viral load, viral load peak time and initial viral growth rate in the TIV and IR models. Figure S14: Relative contribution of macrophage-mediated viral clearance in the TIV and IR models. Table S1: Parameter values for fixed parameters. Table S2: Priors for estimated model parameters.

Author Contributions

Conceptualization, K.L., P.C. and J.M.M.; Formal analysis, K.L., P.C. and J.M.M.; Methodology, K.L., P.C. and J.M.M.; Supervision, P.C. and J.M.M.; Validation, K.L., P.C. and J.M.M.; Visualization, K.L., P.C. and J.M.M.; Writing—original draft, K.L.; Writing—review & editing, K.L., P.C. and J.M.M. All authors have read and agreed to the published version of the manuscript.

Funding

Ke Li is supported by a Melbourne Research Scholarship. This work was supported by an Australian Research Council (ARC) Discovery Project (DP170103076) and a National Health and Medical Research Council (NHMRC) funded Centre for Research Excellence in Infectious Diseases Modelling to Inform Public Health Policy (1078068).

Acknowledgments

We would like to thank Julie McAuley and Lorena Brown for helpful conversations. This research was supported by use of the Nectar Research Cloud, a collaborative Australian research platform supported by the National Collaborative Research Infrastructure Strategy (NCRIS).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Taubenberger, J.K.; Morens, D.M. The pathology of influenza virus infections. Annu. Rev. Pathol. Mech. Dis. 2008, 3, 499–522. [Google Scholar] [CrossRef]
  2. Rosen, D.G.; Lopez, A.E.; Anzalone, M.L.; Wolf, D.A.; Derrick, S.M.; Florez, L.F.; Gonsoulin, M.L.; Hines, M.O.; Mitchell, R.A.; Phatak, D.R.; et al. Postmortem findings in eight cases of influenza A/H1N1. Mod. Pathol. 2010, 23, 1449–1457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Iwasaki, A.; Pillai, P.S. Innate immunity to influenza virus infection. Nat. Rev. Immunol. 2014, 14, 315–328. [Google Scholar] [CrossRef] [PubMed]
  4. Koutsakos, M.; Kedzierska, K.; Subbarao, K. Immune responses to avian influenza viruses. J. Immunol. 2019, 202, 382–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Fujisawa, H.; Tsuru, S.; Taniguchi, M.; Zinnaka, Y.; Nomoto, K. Protective mechanisms against pulmonary infection with influenza virus. I. Relative contribution of polymorphonuclear leukocytes and of alveolar macrophages to protection during the early phase of intranasal infection. J. Gen. Virol. 1987, 68, 425–432. [Google Scholar] [CrossRef]
  6. Hashimoto, Y.; Moki, T.; Takizawa, T.; Shiratsuchi, A.; Nakanishi, Y. Evidence for phagocytosis of influenza virus-infected, apoptotic cells by neutrophils and macrophages in mice. J. Immunol. 2007, 178, 2448–2457. [Google Scholar] [CrossRef] [Green Version]
  7. Swain, S.L.; McKinstry, K.K.; Strutt, T.M. Expanding roles for CD4+ T cells in immunity to viruses. Nat. Rev. Immunol. 2012, 12, 136–148. [Google Scholar] [CrossRef] [PubMed]
  8. Davies, L.C.; Jenkins, S.J.; Allen, J.E.; Taylor, P.R. Tissue-resident macrophages. Nat. Immunol. 2013, 14, 986. [Google Scholar] [CrossRef]
  9. Hussell, T.; Bell, T.J. Alveolar macrophages: Plasticity in a tissue-specific context. Nat. Rev. Immunol. 2014, 14, 81–93. [Google Scholar] [CrossRef]
  10. Kim, H.M.; Lee, Y.W.; Lee, K.J.; Kim, H.S.; Cho, S.W.; Van Rooijen, N.; Guan, Y.; Seo, S.H. Alveolar macrophages are indispensable for controlling influenza viruses in lungs of pigs. J. Virol. 2008, 82, 4265–4274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. La Gruta, N.L.; Kedzierska, K.; Stambas, J.; Doherty, P.C. A question of self-preservation: Immunopathology in influenza virus infection. Immunol. Cell Biol. 2007, 85, 85–92. [Google Scholar] [CrossRef] [PubMed]
  12. Wendy, C.; Chan, R.W.; Wang, J.; Travanty, E.A.; Nicholls, J.M.; Peiris, J.M.; Mason, R.J.; Chan, M.C. Viral replication and innate host responses in primary human alveolar epithelial cells and alveolar macrophages infected with influenza H5N1 and H1N1 viruses. J. Virol. 2011, 85, 6844–6855. [Google Scholar]
  13. Cheung, C.; Poon, L.; Lau, A.; Luk, W.; Lau, Y.; Shortridge, K.; Gordon, S.; Guan, Y.; Peiris, J. Induction of proinflammatory cytokines in human macrophages by influenza A (H5N1) viruses: A mechanism for the unusual severity of human disease? Lancet 2002, 360, 1831–1837. [Google Scholar] [CrossRef]
  14. Perrone, L.A.; Plowden, J.K.; García-Sastre, A.; Katz, J.M.; Tumpey, T.M. H5N1 and 1918 pandemic influenza virus infection results in early and excessive infiltration of macrophages and neutrophils in the lungs of mice. PLoS Pathog. 2008, 4, e1000115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. De Jong, M.D.; Simmons, C.P.; Thanh, T.T.; Hien, V.M.; Smith, G.J.; Chau, T.N.B.; Hoang, D.M.; Chau, N.V.V.; Khanh, T.H.; Dong, V.C.; et al. Fatal outcome of human influenza A (H5N1) is associated with high viral load and hypercytokinemia. Nat. Med. 2006, 12, 1203–1207. [Google Scholar] [CrossRef] [PubMed]
  16. Apostolopoulos, V.; Stojanovska, L.; Gargosky, S.E. MUC1 (CD227): A multi-tasked molecule. Cell. Mol. Life Sci. 2015, 72, 4475–4500. [Google Scholar] [CrossRef] [PubMed]
  17. Dhar, P.; McAuley, J. The role of the cell surface mucin MUC1 as a barrier to infection and regulator of inflammation. Front. Cell. Infect. Microbiol. 2019, 9, 117. [Google Scholar] [CrossRef] [PubMed]
  18. Kim, K.C.; Lillehoj, E.P. MUC1 mucin: A peacemaker in the lung. Am. J. Respir. Cell Mol. Biol. 2008, 39, 644–647. [Google Scholar] [CrossRef] [PubMed]
  19. Baños-Lara, M.; Rocío, D.; Piao, B.; Guerrero-Plata, A. Differential mucin expression by respiratory syncytial virus and human metapneumovirus infection in human epithelial cells. Mediat. Inflamm. 2015, 2015. [Google Scholar] [CrossRef] [Green Version]
  20. Kato, K.; Lillehoj, E.P.; Kim, K.C. MUC1 regulates epithelial inflammation and apoptosis by PolyI: C through Inhibition of Toll/IL-1 receptor-domain–containing adapter-inducing IFN-β (TRIF) recruitment to toll-like receptor 3. Am. J. Respir. Cell Mol. Biol. 2014, 51, 446–454. [Google Scholar] [CrossRef] [Green Version]
  21. Li, Y.; Dinwiddie, D.L.; Harrod, K.S.; Jiang, Y.; Kim, K.C. Anti-inflammatory effect of MUC1 during respiratory syncytial virus infection of lung epithelial cells in vitro. Am. J. Physiol.-Lung Cell. Mol. Physiol. 2010, 298, L558–L563. [Google Scholar] [CrossRef]
  22. Dhar, P.; Ng, G.Z.; Dunne, E.M.; Sutton, P. Mucin 1 protects against severe Streptococcus pneumoniae infection. Virulence 2017, 8, 1631–1642. [Google Scholar] [CrossRef] [Green Version]
  23. Lu, W.; Hisatsune, A.; Koga, T.; Kato, K.; Kuwahara, I.; Lillehoj, E.P.; Chen, W.; Cross, A.S.; Gendler, S.J.; Gewirtz, A.T.; et al. Cutting edge: Enhanced pulmonary clearance of Pseudomonas aeruginosa by Muc1 knockout mice. J. Immunol. 2006, 176, 3890–3894. [Google Scholar] [CrossRef] [Green Version]
  24. McAuley, J.; Corcilius, L.; Tan, H.; Payne, R.; McGuckin, M.; Brown, L. The cell surface mucin MUC1 limits the severity of influenza A virus infection. Mucosal Immunol. 2017, 10, 1581–1593. [Google Scholar] [CrossRef] [Green Version]
  25. Handel, A.; Liao, L.E.; Beauchemin, C.A. Progress and trends in mathematical modelling of influenza A virus infections. Curr. Opin. Syst. Biol. 2018, 12, 30–36. [Google Scholar] [CrossRef]
  26. Miao, H.; Hollenbaugh, J.A.; Zand, M.S.; Holden-Wiltse, J.; Mosmann, T.R.; Perelson, A.S.; Wu, H.; Topham, D.J. Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza A virus. J. Virol. 2010, 84, 6687–6698. [Google Scholar] [CrossRef] [Green Version]
  27. Cao, P.; Yan, A.W.; Heffernan, J.M.; Petrie, S.; Moss, R.G.; Carolan, L.A.; Guarnaccia, T.A.; Kelso, A.; Barr, I.G.; McVernon, J.; et al. Innate immunity and the inter-exposure interval determine the dynamics of secondary influenza virus infection and explain observed viral hierarchies. PLoS Comput. Biol. 2015, 11, e1004334. [Google Scholar] [CrossRef]
  28. Cao, P.; Wang, Z.; Yan, A.W.; McVernon, J.; Xu, J.; Heffernan, J.M.; Kedzierska, K.; McCaw, J.M. On the role of CD8+ T cells in determining recovery time from influenza virus infection. Front. Immunol. 2016, 7, 611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kim, H.M.; Kang, Y.M.; Ku, K.B.; Park, E.H.; Yum, J.; Kim, J.C.; Jin, S.Y.; Lee, J.S.; Kim, H.S.; Seo, S.H. The severe pathogenicity of alveolar macrophage-depleted ferrets infected with 2009 pandemic H1N1 influenza virus. Virology 2013, 444, 394–403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Schneider, C.; Nobs, S.P.; Heer, A.K.; Kurrer, M.; Klinke, G.; Van Rooijen, N.; Vogel, J.; Kopf, M. Alveolar macrophages are essential for protection from respiratory failure and associated morbidity following influenza virus infection. PLoS Pathog. 2014, 10, e1004053. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A. On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Baccam, P.; Beauchemin, C.; Macken, C.A.; Hayden, F.G.; Perelson, A.S. Kinetics of influenza A virus infection in humans. J. Virol. 2006, 80, 7590–7599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Smith, A.M.; Adler, F.R.; McAuley, J.L.; Gutenkunst, R.N.; Ribeiro, R.M.; McCullers, J.A.; Perelson, A.S. Effect of 1918 PB1-F2 expression on influenza A virus infection kinetics. PLoS Comput. Biol. 2011, 7, e1001081. [Google Scholar] [CrossRef] [Green Version]
  34. Hayden, F.G.; Fritz, R.; Lobo, M.C.; Alvord, W.; Strober, W.; Straus, S.E. Local and systemic cytokine responses during experimental human influenza A virus infection. Relation to symptom formation and host defense. J. Clin. Investig. 1998, 101, 643–649. [Google Scholar] [CrossRef] [PubMed]
  35. Cao, P.; McCaw, J.M. The mechanisms for within-host influenza virus control affect model-based assessment and prediction of antiviral treatment. Viruses 2017, 9, 197. [Google Scholar] [CrossRef]
  36. Pawelek, K.A.; Huynh, G.T.; Quinlivan, M.; Cullinane, A.; Rong, L.; Perelson, A.S. Modeling within-host dynamics of influenza virus infection including immune responses. PLoS Comput. Biol. 2012, 8, e1002588. [Google Scholar] [CrossRef] [Green Version]
  37. Yan, A.W.; Zaloumis, S.G.; Simpson, J.A.; McCaw, J.M. Sequential infection experiments for quantifying innate and adaptive immunity during influenza infection. PLoS Comput. Biol. 2019, 15, e1006568. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, J.; Nikrad, M.P.; Travanty, E.A.; Zhou, B.; Phang, T.; Gao, B.; Alford, T.; Ito, Y.; Nahreini, P.; Hartshorn, K.; et al. Innate immune response of human alveolar macrophages during influenza A infection. PLoS ONE 2012, 7, e29879. [Google Scholar] [CrossRef] [Green Version]
  39. Yan, A.W.; Cao, P.; Heffernan, J.M.; McVernon, J.; Quinn, K.M.; La Gruta, N.L.; Laurie, K.L.; McCaw, J.M. Modelling cross-reactivity and memory in the cellular adaptive immune response to influenza infection in the host. J. Theor. Biol. 2017, 413, 34–49. [Google Scholar] [CrossRef] [Green Version]
  40. Hoffman, M.D.; Gelman, A. The No-U-Turn sampler: Adaptively setting path lengths in Hamiltonian Monte Carlo. J. Mach. Learn. Res. 2014, 15, 1593–1623. [Google Scholar]
  41. Chatzilena, A.; van Leeuwen, E.; Ratmann, O.; Baguelin, M.; Demiris, N. Contemporary statistical inference for infectious disease models using Stan. Epidemics 2019, 29, 100367. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Results of model fitting for WT and MUC1-KO mice. Data are presented by solid circles. Panels (A,B) show the median of posterior predictions (solid line) and a 95% prediction interval (shaded area) of viral load data for both TIV (red) and IR (yellow) models for WT and MUC1-KO mice, respectively. Panels (C,D) show the model predictions of macrophage data in the two models for WT and MUC1-KO mice, respectively. The priors of model parameters are given in Supplementary Materials. The posteriors of estimated model parameters are given in Supplementary Figures.
Figure 1. Results of model fitting for WT and MUC1-KO mice. Data are presented by solid circles. Panels (A,B) show the median of posterior predictions (solid line) and a 95% prediction interval (shaded area) of viral load data for both TIV (red) and IR (yellow) models for WT and MUC1-KO mice, respectively. Panels (C,D) show the model predictions of macrophage data in the two models for WT and MUC1-KO mice, respectively. The priors of model parameters are given in Supplementary Materials. The posteriors of estimated model parameters are given in Supplementary Figures.
Viruses 13 00850 g001
Figure 2. Comparison of model predictions for selected key biological quantities. Distributions are calculated using the 4000 joint posterior distributions. Panels (A,B) show the distribution of the basic reproduction number of viral replication in wildtype (purple) and MUC1-knockout (green) group in TIV (left panel) and IR models (right panel), respectively. Panels (C,D) show the distribution of the cumulative viral load in different mice groups in the two models. Panels (E,F) show the accumulative macrophages in WT and MUC1-KO mice group in the two models.
Figure 2. Comparison of model predictions for selected key biological quantities. Distributions are calculated using the 4000 joint posterior distributions. Panels (A,B) show the distribution of the basic reproduction number of viral replication in wildtype (purple) and MUC1-knockout (green) group in TIV (left panel) and IR models (right panel), respectively. Panels (C,D) show the distribution of the cumulative viral load in different mice groups in the two models. Panels (E,F) show the accumulative macrophages in WT and MUC1-KO mice group in the two models.
Viruses 13 00850 g002
Figure 3. The dependence of the AUC M on the effects of MUC1 for the TIV model. The upper panel shows the marginal posterior distribution of ε 1 (left) and ε 2 (right). The 95% credible interval (CI) for the parameters is indicated between the two red-dashed lines, and the red-solid line indicates parameters’ median value. The heatmap shows dependence of macrophage reduction efficiency on ε 1 and ε 2 . The black circle indicates the pair of median values of ε 1 and ε 2 , and the arrow indicates the direction of the rate of change in macrophage reduction efficiency at that point.
Figure 3. The dependence of the AUC M on the effects of MUC1 for the TIV model. The upper panel shows the marginal posterior distribution of ε 1 (left) and ε 2 (right). The 95% credible interval (CI) for the parameters is indicated between the two red-dashed lines, and the red-solid line indicates parameters’ median value. The heatmap shows dependence of macrophage reduction efficiency on ε 1 and ε 2 . The black circle indicates the pair of median values of ε 1 and ε 2 , and the arrow indicates the direction of the rate of change in macrophage reduction efficiency at that point.
Viruses 13 00850 g003
Figure 4. The dependence of the AUC M on the effects of MUC1 for the IR model. The upper panel shows the marginal posterior distribution of ε 1 (left) and ε 2 (right). The 95% credible interval (CI) for the parameters is indicated between the two red-dashed lines, and the red-solid line indicates parameters’ median value. The heatmap shows dependence of macrophage reduction efficiency on ε 1 and ε 2 . The black circle indicates the pair of median values of ε 1 and ε 2 , and the arrow indicates the direction of the rate of change in macrophage reduction efficiency at that point.
Figure 4. The dependence of the AUC M on the effects of MUC1 for the IR model. The upper panel shows the marginal posterior distribution of ε 1 (left) and ε 2 (right). The 95% credible interval (CI) for the parameters is indicated between the two red-dashed lines, and the red-solid line indicates parameters’ median value. The heatmap shows dependence of macrophage reduction efficiency on ε 1 and ε 2 . The black circle indicates the pair of median values of ε 1 and ε 2 , and the arrow indicates the direction of the rate of change in macrophage reduction efficiency at that point.
Viruses 13 00850 g004
Table 1. Estimates of MUC1 parameters and comparison between models. The estimates of MUC1 on reduction of target cell susceptibility to influenza virus ( ε 1 ) and on reduction of macrophage recruitment rate induced by infected cells ( ε 2 ). The lower and upper boundary of the 95% credible interval (CI) of the parameter is given by calculating the 2.5% and 97.5% quantile of the marginal posterior parameter distribution.
Table 1. Estimates of MUC1 parameters and comparison between models. The estimates of MUC1 on reduction of target cell susceptibility to influenza virus ( ε 1 ) and on reduction of macrophage recruitment rate induced by infected cells ( ε 2 ). The lower and upper boundary of the 95% credible interval (CI) of the parameter is given by calculating the 2.5% and 97.5% quantile of the marginal posterior parameter distribution.
ParameterDescriptionMedian (95% CI)
TIVIR
ε 1 The reduction in target cell susceptibility to infection due to MUC10.44 (0.23, 0.71)0.42 (0.22, 0.58)
ε 2 The reduction in recruitment rate of macrophages due to MUC10.45 (0.18, 0.64)0.38 (0.06, 0.63)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, K.; Cao, P.; McCaw, J.M. Modelling the Effect of MUC1 on Influenza Virus Infection Kinetics and Macrophage Dynamics. Viruses 2021, 13, 850. https://doi.org/10.3390/v13050850

AMA Style

Li K, Cao P, McCaw JM. Modelling the Effect of MUC1 on Influenza Virus Infection Kinetics and Macrophage Dynamics. Viruses. 2021; 13(5):850. https://doi.org/10.3390/v13050850

Chicago/Turabian Style

Li, Ke, Pengxing Cao, and James M. McCaw. 2021. "Modelling the Effect of MUC1 on Influenza Virus Infection Kinetics and Macrophage Dynamics" Viruses 13, no. 5: 850. https://doi.org/10.3390/v13050850

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop