- Research note
- Open Access
A mathematical model relates intracellular TLR4 oscillations to sepsis progression
BMC Research Notes volume 11, Article number: 462 (2018)
Oscillations of physiological parameters describe many biological processes and their modulation is determinant for various pathologies. In sepsis, toll-like receptor 4 (TLR4) is a key sensor for signaling the presence of Gram-negative bacteria. Its intracellular trafficking rates shift the equilibrium between the pro- and anti-inflammatory downstream signaling cascades, leading to either the physiological resolution of the bacterial stimulation or to sepsis. This study aimed to evaluate the effects of TLR4 increased expression and intracellular trafficking on the course and outcome of sepsis.
Using a set of three differential equations, we defined the TLR4 fluxes between relevant cell organelles. We obtained three different regions in the phase space: (1) a limit-cycle describing unstimulated physiological oscillations, (2) a fixed-point attractor resulting from moderate LPS stimulation that is resolved and (3) a double-attractor resulting from sustained LPS stimulation that leads to sepsis. We used this model to describe available hospital data of sepsis patients and we correctly characterize the clinical outcome of these patients.
The immune system is replete with oscillations of various parameters needed for mounting an appropriate response upon stimulation. TLR4 is an important bacterial recognition that serves as a link between innate and adaptive immunity . In particular, TLR4 experiences a significant upregulation in mRNA production and presentation to the cell surface at the initial stages of sepsis in both humans and experimental models . An emerging theme in TLR4 signaling posits that its cellular localization is determinant for its functions . Throughout the continuum of sepsis, complete TLR4 signaling includes not only the initial surface-bound pro-inflammatory signaling, but also its subsequent endocytosis and intracellular trafficking. This results in competing endosomal anti-inflammatory cytokine production and further into either receptor recycling to cell membrane or signal termination within endolysosomes. Initial responsiveness to LPS is therefore regulated by the concentration of cell surface TLR4 that depends in turn on TLR4 expression, on TLR4 trafficking from the Golgi apparatus to the plasma membrane and on the amount of TLR4 already internalized in endosomes .
Mechanistic model of TLR4 trafficking
An overview of the known TLR4 intracellular trafficking routes that influence its signaling is presented in Fig. 1.
Upon endotoxin stimulation, initial TLR4 immobilization (step 1) may lead to monomeric LPS being internalized and trafficked to the Golgi apparatus within seconds of stimulation, without activating TLR4 (step 2). This is followed by TLR4 clustering (step 3) with monomeric LPS . Internalization by either clathrin- and dynamin-mediated processes  results in a switch in TLR4 signaling pathways by means of different adaptors (step 4). Provided that TLR4 endocytosis has occurred, signaling continues induction of type-I interferons  from early endosomes (step 5). For the signal to be terminated (step 6), the TLR4 complex is ubiquitinated and marked for lysosomal degradation. Alternatively, TLR4 can be recycled for new signaling cycles back to the cell surface via the endosomal recycling compartment (step 7). Under limit-cycle unstimulated physiological oscillations, cell surface TLR4 are present in low concentrations in macrophages or are undetectable in dendritic cells, with most resident TLR4 being distributed in the Golgi apparatus . Rapid TLR4 mobilization to cell membrane follows LPS activation , canceling the downregulation present in physiological conditions that serves to desensitize cells to low endotoxin levels. While the overall sequence of TLR4 activation has been elucidated, the rates of TLR4 trafficking are not quantified, nor are available absolute numbers for TLR4 expression on cell surfaces.
Simulation of TLR4 trafficking routes
In order to simulate in silico the initial TLR4 trafficking events between the endosomal recycling compartment (ERC) and the trans-Golgi network (TGN) to and from cell surface and within the early endosomes–endolysosome (EE) system, we have constructed a dynamic model (Additional file 1), based on the three ordinary differential equations presented below:
where: x = concentration of TLR4 in TGN and ERC, y = concentration of TLR4 in endosomes/endolysosomes (EE), z = concentration of TLR4 on cell surface, ϕ = rate of TLR4 mRNA production, β = rate of TLR4 trafficked to lysosomes from endosomes, α = rate of TLR4 retroactively trafficked to ERC from endosomes, γ = rate of TLR4 on cell surface trafficked to TGN, σ = rate of TLR4 on cell surface trafficked to endosomal system.
The TLR4 flux in the system as indicated by Eq. (1) is influenced by the TRAM distribution within ERC that shifts onto the enlarged CD14/LPS-positive endosomes upon TLR4 activation . The adaptor TRAM is also constitutively present at the plasma membrane anchored at a N-terminal myristoylation site and traffics concomitantly the TLR4 signaling complex unidirectionaly to the endosomal system . This synergy allows for the anti-inflammatory signaling phase to take preponderance, possibly due to unique TLR4 conformation brought on by the endosomal acidic environment, as previously proposed . These events are dominant after about 30 min upon LPS stimulation , allowing for TLR4 to traffic, in a first stage mostly bidirectionally from the ERC to EE (Eq. 2). The small GTPase Rab7b is upregulated upon LPS exposure in the early endosomes and is a key regulator of intracellular trafficking of the TLR4 signaling complex to either late endosomes/lysosomes for signal termination (parameter β), or to ERC (parameter α) . In Rab7b-silenced macrophages, after LPS stimulation, continued TLR4 presence only in the EE system has adverse effects as to its prolonged anti-inflammatory signaling . Equation (3) describes the TLR4 cell surface concentration changes as the difference between the pool of available TLR4 in TGN + ERC and in EE, and the TLR4 that is actively being prevented from clustering on the cell surface (so as to increase downstream signaling), be it directly from the surface towards TGN (parameter γ) or towards EE (parameter σ).
Results and discussion
In the absence of available data from the literature, the parameter values for ϕ, β, γ, σ and ϕ were varied until a stable limit cycle was obtained, a characteristic of protein expression, as previously shown . The first four parameters were kept constant to reflect the steady-state, non-stimulated oscillations in the TLR4 intracellular trafficking, while the ϕ-parameter that was determined experimentally elsewhere  was allowed to vary. Using Eqs. 1–3, we sought to model the cellular regimes that are impacted by the overall TLR4 sensitivity to LPS, as reflected by the initial rate of tlr4 mRNA synthesis, upon sepsis diagnosis and prior to clinical intervention. The ϕ-parameter augments markedly in experimental models of sepsis and directly correlates with mortality, with peak increases between 1 and 3 h post sepsis induction . As such, its oscillations will determine TLR4 changes within various relevant cell compartments and dictate the timing and preponderance of the pro- and anti-inflammatory responses. We defined three regions in the phase space for the plasma membrane and intracellular TLR4 distribution, based on the variations of ϕ-parameter (Additional file 2): (i) a steady-state with TLR4 expression and concentration oscillating within a narrow margin throughout the relevant cell compartments, (ii) a low to medium tlr4 mRNA production following LPS stimulation that results in an initial increase of TLR4 concentration on the cell surface and subsequently in the endosomal system, followed by a regulated decrease, (iii) a third, high tlr4 mRNA output matching increasing LPS stimulation where TLR4 concentrations oscillate stably and irreversibly on the cell surface and within the EE. The variations in tlr4 mRNA measured in the patients served as the first parameter (ϕ) to be changed, responsible for initial TLR4 distribution within the relevant cell compartments. TLR4 is unique among other pathogen-recognition receptors in that its intracellular trafficking is determinant for the inflammatory signaling it initiates . Depending on initial conditions and rate changes, the ensuing orbits either approach stable fixed points or undergo variations, each having a different physiological interpretation, as presented in Fig. 2.
a. Physiological variations in TLR4 concentrations
In all simulations, we assumed that initial expression levels on the cell surface, TGN/ERC and EE are low. For steady-state conditions, we proposed that TLR4 concentration oscillations are of low amplitude, reflecting the experimental data on tlr4 mRNA in human monocytes in vitro . A stable limit cycle is achieved with ϕ = 1.2, β = 3.6, α = 1.2, γ = 2.4, σ = 1.3 (Fig. 2a).
b. Sepsis progression and resolution
We propose that following a moderate LPS stimulation, TLR4 levels initially increase in order to proportionally signal the Gram-negative bacterial presence . A fixed-point attractor is obtained with ϕ < 1.2, β = 3.6, α = 1.2, γ = 2.4, σ = 1.3 (Fig. 2b).
c. Sepsis progression and mortality
Upon increasing LPS stimulation, we assumed that tlr4 mRNA rates are amplified proportionally, the result of which leads to the system moving to a double-attractor. TLR4 concentrations oscillate with highest amplitude and indefinitely between cell surface and EE compartments, with no signal resolution, using ϕ > 1.2, β = 3.6, α = 1.2, γ = 2.4, σ = 1.3 (Fig. 2c).
We note that this model, when using the initial, pre-treatment rates of tlr4 mRNA from the patient cohort yielded appropriate descriptions of both the clinical outcome in 8 out of 10 patients , and the category of attractor each patient belongs to (Additional file 2). Patients whose TLR4 concentrations changes evolved towards one attractor survived sepsis (patients #1, 4, 5, and 8). In contrast, those patients that presented a double-attractor state for TLR4 died within 3 days after ICU admission (patients #3, 6, 7, and 10). As a test to the sensitivity and specificity of our model, patient #2 died 9 days after ICU admittance due to Candida albicans infection, a pathogen known to stimulate both TLR2 and TLR4 and is commonly associated with severe immunosuppression . Patient #9 survived with negative microbiological cultures, a result of false-negative cultures or sepsis without infection, a situation not accounted for in our model. We have used initial tlr4 mRNA expression levels from sepsis patients in a dynamic model in order to describe the distribution of TLR4 within the cell surface compartment (pro-inflammatory role), or intracellularly (anti-inflammatory and signal termination functions). We discriminated Gram-negative infections from the overall cohort and correctly described the clinical outcome of 8/10 patients. Confirming this model with in vivo measurements of TLR4 intracellular trafficking rates would provide further insight into their contribution to sepsis onset and progression.
The study was unable to account for phagosome signaling of whole Gram negative bacteria, nor for the additional TLR4 subpopulation that trafficks from ERC to phagosome after LPS stimulation.
endosomal recycling compartment
Medzhitov R, Preston-Hurlburt P, Janeway CA. A human homologue of the Drosophila Toll protein signals activation of adaptive immunity. Nature. 1997;388:394–7. https://doi.org/10.1038/41131.
Williams DL, Ha T, Li C, Kalbfleisch JH, Schweitzer J, Vogt W, Browder IW. Modulation of tissue Toll-like receptor 2 and 4 during the early phases of polymicrobial sepsis correlates with mortality. Crit Care Med. 2003;31:1808–18. https://doi.org/10.1097/01.CCM.0000069343.27691.F3.
Chaturvedi A, Pierce SK. How location governs toll-like receptor signaling. Traffic. 2009;10:621–8. https://doi.org/10.1111/j.1600-0854.2009.00899.x.
Tsujimoto H, Ono S, Efron PA, Scumpia PO, Moldawer LL, Mochizuki H. Role of Toll-like receptors in the development of sepsis. Shock. 2008;29:315–21. https://doi.org/10.1097/SHK.0b013e318157ee55.
Klein DC, Skjesol A, Kers-Rebel ED, Sherstova T, Sporsheim B, Egeberg KW, Stokke BT, Espevik T, Husebye H. CD14, TLR4 and TRAM show different trafficking dynamics during LPS stimulation. Traffic. 2015;16:677–90. https://doi.org/10.1111/tra.12274.
Husebye H, Halaas Ø, Stenmark H, Tunheim G, Sandanger Ø, Bogen B, Brech A, Latz E, Espevik T. Endocytic pathways regulate Toll-like receptor 4 signaling and link innate and adaptive immunity. EMBO J. 2006;25:683–92. https://doi.org/10.1038/sj.emboj.7600991.
Kagan JC, Su T, Horng T, Chow A, Akira S, Medzhitov R. TRAM couples endocytosis of Toll-like receptor 4 to the induction of interferon-beta. Nat Immunol. 2009;9:361–8. https://doi.org/10.1038/ni1569.
Uronen-Hansson H, Allen J, Osman M, Squires G, Klein N, Callard RE. Toll-like receptor 2 (TLR2) and TLR4 are present inside human dendritic cells, associated with microtubules and the Golgi apparatus but are not detectable on the cell surface: integrity of microtubules is required for interleukin-12 production in response to internalized bacteria. Immunology. 2004;111:173–8. https://doi.org/10.1111/j.0019-2805.2003.01803.x.
Rocuts F, Ma Y, Zhang X, Gao W, Yue Y, Vartanian T, Wang H. Carbon monoxide suppresses membrane expression of TLR4 via myeloid differentiation factor-2 in betaTC3 cells. J Immunol. 2010;185:2134–9. https://doi.org/10.3389/fimmu.2014.00473.
Brubaker SW, Bonham KS, Zanoni I, Kagan JC. Innate immune pattern recognition: a cell biological perspective. Annu Rev Immunol. 2015;33:257–90. https://doi.org/10.1146/annurev-immunol-032414-112240.
Verstak B, Stack J, Ve T, Mangan M, Hjerrild K, Jeon J, Stahl R, Latz E, Gay N, Kobe B, Bowie AG, Mansell A. The TLR signaling adaptor TRAM interacts with TRAF6 to mediate activation of the inflammatory response by TLR4. J Leukoc Biol. 2014;96:427–36. https://doi.org/10.1189/jlb.2A0913-487R.
Gangloff M. Different dimerisation mode for TLR4 upon endosomal acidification? Trends Biochem Sci. 2012;37:92–8. https://doi.org/10.1016/j.tibs.2011.11.003.
Distefano MB, Kjos I, Bakke O, Progida C. Rab7b at the intersection of intracellular trafficking and cell migration. Commun Integr Biol. 2015;8(6):e1023492. https://doi.org/10.1080/19420889.2015.1023492.
Wang Y, Chen T, Han C, He D, Liu H, An H, Cai Z, Cao X. Lysosome-associated small Rab GTPase Rab7b negatively regulates TLR4 signaling in macrophages by promoting lysosomal degradation of TLR4. Blood. 2007;110:962–71. https://doi.org/10.1182/blood-2007-01-066027.
Swain PS, Elowitz MB, Siggia ED. Intrinsic and extrinsic contributions to stochasticity in gene expression. Proc Natl Acad Sci. 2002;99:12795–800. https://doi.org/10.1073/pnas.162041399.
Stan RS, Bonin C, Porto R, Soriano F, Camargo MM. Negative correlation between tlr4 and grp78 expression is characteristic of sepsis onset and progression. bioRxiv. 2017. https://doi.org/10.1101/133264.
Togbe D, Schnyder-Candrian S, Schnyder B, Couillin I, Maillet I, Bihl F, Malo D, Ryffel B, Quesniaux VF. TLR4 gene dosage contributes to endotoxin-induced acute respiratory inflammation. J Leukoc Biol. 2006;80:451–7. https://doi.org/10.1189/jlb.0206099.
Kollef M, Micek S, Hampton N, Doherty JA, Kumar A. Septic shock attributed to Candida infection: importance of empiric therapy and source control. Clin Infect Dis. 2012;54:739–1746. https://doi.org/10.1093/cid/cis305.
MMDC and FGS performed the original experimental work whose data this study is partly built on. MMDC and RCS obtained funding for the work. RCS and MMDC participated in the study design and drafted the manuscript. RCS prepared the Mathematica files. All authors read and approved the final manuscript.
We thank Dr. Grégoire Altan-Bonnet for critical reading of the manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Consent for publication
Ethics approval and consent to participate
This study was supported by CNPq (400662/2014-0 for RCS, 309041/2012-0 for MMDC), FAPESP (11/51778-6 for MMDC).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.