Mathematical model for glutathione dynamics in the retina.

The retina is highly susceptible to the generation of toxic reactive oxygen species (ROS) that disrupt the normal operations of retinal cells. The glutathione (GSH) antioxidant system plays an important role in mitigating ROS. To perform its protective functions, GSH depends on nicotinamide adenine dinucleotide phosphate (NADPH) produced through the pentose phosphate pathway. This work develops the first mathematical model for the GSH antioxidant system in the outer retina, capturing the most essential components for formation of ROS, GSH production, its oxidation in detoxifying ROS, and subsequent reduction by NADPH. We calibrate and validate the model using experimental measurements, at different postnatal days up to PN28, from control mice and from the rd1 mouse model for the disease retinitis pigmentosa (RP). Global sensitivity analysis is then applied to examine the model behavior and identify the pathways with the greatest impact in control compared to RP conditions. The findings underscore the importance of GSH and NADPH production in dealing with oxidative stress during retinal development, especially after peak rod degeneration occurs in RP, leading to increased oxygen tension. This suggests that stimulation of GSH and NADPH synthesis could be a potential intervention strategy in degenerative mouse retinas with RP.

Insights into pathological mechanisms and interventions revealed by analyzing a mathematical model for cone metabolism.

This work analyzes a mathematical model for the metabolic dynamics of a cone photoreceptor, which is the first model to account for energy generation from fatty acids oxidation of shed photoreceptor outer segments (POS). Multiple parameter bifurcation analysis shows that joint variations in external glucose, the efficiency of glucose transporter 1 (GLUT1), lipid utilization for POS renewal, and oxidation of fatty acids affect the cone's metabolic vitality and its capability to adapt under glucose-deficient conditions. The analysis further reveals that when glucose is scarce, cone viability cannot be sustained by only fueling energy production in the mitochondria, but it also requires supporting anabolic processes to create lipids necessary for cell maintenance and repair. In silico experiments are used to investigate how the duration of glucose deprivation impacts the cell without and with a potential GLUT1 or oxidation of fatty acids intervention as well as a dual intervention. The results show that for prolonged duration of glucose deprivation, the cone metabolic system does not recover with higher oxidation of fatty acids and requires greater effectiveness of GLUT1 to recover. Finally, time-varying global sensitivity analysis (GSA) is applied to assess the sensitivity of the model outputs of interest to changes and uncertainty in the parameters at specific times. The results reveal a critical temporal window where there would be more flexibility for interventions to rescue a cone cell from the detrimental consequences of glucose shortage.

Data assimilation of synthetic data as a novel strategy for predicting disease progression in alopecia areata.

The goal of patient-specific treatment of diseases requires a connection between clinical observations with models that are able to accurately predict the disease progression. Even when realistic models are available, it is very difficult to parameterize them and often parameter estimates that are made using early time course data prove to be highly inaccurate. Inaccuracies can cause different predictions, especially when the progression depends sensitively on the parameters. In this study, we apply a Bayesian data assimilation method, where the data are incorporated sequentially, to a model of the autoimmune disease alopecia areata that is characterized by distinct spatial patterns of hair loss. Using synthetic data as simulated clinical observations, we show that our method is relatively robust with respect to variations in parameter estimates. Moreover, we compare convergence rates for parameters with different sensitivities, varying observational times and varying levels of noise. We find that this method works better for sparse observations, sensitive parameters and noisy observations. Taken together, we find that our data assimilation, in conjunction with our biologically inspired model, provides directions for individualized diagnosis and treatments.

A physiological model of the inflammatory-thermal-pain-cardiovascular interactions during an endotoxin challenge.

KEY POINTS: Inflammation in response to bacterial endotoxin challenge impacts physiological functions, including cardiovascular, thermal and pain dynamics, although the mechanisms are poorly understood. We develop an innovative mathematical model incorporating interaction pathways between inflammation and physiological processes observed in response to an endotoxin challenge. We calibrate the model to individual data from 20 subjects in an experimental study of the human inflammatory and physiological responses to endotoxin, and we validate the model against human data from an independent study. Using the model to simulate patient responses to different treatment modalities reveals that a multimodal treatment combining several therapeutic strategies gives the best recovery outcome. ABSTRACT: Uncontrolled, excessive production of pro-inflammatory mediators from immune cells and traumatized tissues can cause systemic inflammatory conditions such as sepsis, one of the ten leading causes of death in the USA, and one of the three leading causes of death in the intensive care unit. Understanding how inflammation affects physiological processes, including cardiovascular, thermal and pain dynamics, can improve a patient's chance of recovery after an inflammatory event caused by surgery or a severe infection. Although the effects of the autonomic response on the inflammatory system are well-known, knowledge about the reverse interaction is lacking. The present study develops a mathematical model analyzing the inflammatory system's interactions with thermal, pain and cardiovascular dynamics in response to a bacterial endotoxin challenge. We calibrate the model with individual data from an experimental study of the inflammatory and physiological responses to a one-time administration of endotoxin in 20 healthy young men and validate it against data from an independent endotoxin study. We use simulation to explore how various treatments help patients exposed to a sustained pathological input. The treatments explored include bacterial endotoxin adsorption, antipyretics and vasopressors, as well as combinations of these. Our findings suggest that the most favourable recovery outcome is achieved by a multimodal strategy, combining all three interventions to simultaneously remove endotoxin from the body and alleviate symptoms caused by the immune system as it fights the infection.

Toward Predicting the Spatio-Temporal Dynamics of Alopecia Areata Lesions Using Partial Differential Equation Analysis.

Hair loss in the autoimmune disease, alopecia areata (AA), is characterized by the appearance of circularly spreading alopecic lesions in seemingly healthy skin. The distinct spatial patterns of AA lesions form because the immune system attacks hair follicle cells that are in the process of producing hair shaft, catapults the mini-organs that produce hair from a state of growth (anagen) into an apoptosis-driven regression state (catagen), and causes major hair follicle dystrophy along with rapid hair shaft shedding. In this paper, we develop a model of partial differential equations (PDEs) to describe the spatio-temporal dynamics of immune system components that clinical and experimental studies show are primarily involved in the disease development. Global linear stability analysis reveals there is a most unstable mode giving rise to a pattern. The most unstable mode indicates a spatial scale consistent with results of the humanized AA mouse model of Gilhar et al. (Autoimmun Rev 15(7):726-735, 2016) for experimentally induced AA lesions. Numerical simulations of the PDE system confirm our analytic findings and illustrate the formation of a pattern that is characteristic of the spatio-temporal AA dynamics. We apply marginal linear stability analysis to examine and predict the pattern propagation.

Analysing the dynamics of a model for alopecia areata as an autoimmune disorder of hair follicle cycling.

Alopecia areata (AA) is a CD8$^{+}$ T cell-dependent autoimmune disease that disrupts the constantly repeating cyclic transformations of hair follicles (HFs). Among the three main HF cycle stages-growth (anagen), regression (catagen) and relative quiescence (telogen)-only anagen HFs are attacked and thereby forced to prematurely enter into catagen, thus shortening active hair growth substantially. After having previously modelled the dynamics of immune system components critically involved in the disease development (Dobreva et al., 2015), we here present a mathematical model for AA which incorporates HF cycling and illustrates the anagen phase interruption in AA resulting from an inflammatory autoimmune response against HFs. The model couples a system describing the dynamics of autoreactive immune cells with equations modelling the hair cycle. We illustrate states of health, disease and treatment as well as transitions between them. In addition, we perform parameter sensitivity analysis to assess how different processes, such as proliferation, apoptosis and input from stem cells, impact anagen duration in healthy versus AA-affected HFs. The proposed model may help in evaluating the effectiveness of existing treatments and identifying new potential therapeutic targets.

Mathematical model for alopecia areata.

Alopecia areata (AA) is an autoimmune disease, and its clinical phenotype is characterized by the formation of distinct hairless patterns on the scalp or other parts of the body. In most cases hair falls out in round patches. A well-established hypothesis for the pathogenesis of AA states that collapse of hair follicle immune privilege is one of the essential elements in disease development. To investigate the dynamics of alopecia areata, we develop a mathematical model that incorporates immune system components and hair follicle immune privilege agents whose involvement in AA has been confirmed in clinical studies and experimentally. We perform parameter sensitivity analysis in order to determine which inputs have the greatest effect on outcome variables. Our findings suggest that, among all processes reflected in the model, immune privilege guardians and the pro-inflammatory cytokine interferon-γ govern disease dynamics. These results agree with the immune privilege collapse hypothesis for the development of AA.