Bayesian and deep-learning models applied to the early detection of ovarian cancer using multiple longitudinal biomarkers.

BACKGROUND: Ovarian cancer is the most lethal of all gynecological cancers. Cancer Antigen 125 (CA125) is the best-performing ovarian cancer biomarker which however is still not effective as a screening test in the general population. Recent literature reports additional biomarkers with the potential to improve on CA125 for early detection when using longitudinal multimarker models. METHODS: Our data comprised 180 controls and 44 cases with serum samples sourced from the multimodal arm of UK Collaborative Trial of Ovarian Cancer Screening (UKCTOCS). Our models were based on Bayesian change-point detection and recurrent neural networks. RESULTS: We obtained a significantly higher performance for CA125-HE4 model using both methodologies (AUC 0.971, sensitivity 96.7% and AUC 0.987, sensitivity 96.7%) with respect to CA125 (AUC 0.949, sensitivity 90.8% and AUC 0.953, sensitivity 92.1%) for Bayesian change-point model (BCP) and recurrent neural networks (RNN) approaches, respectively. One year before diagnosis, the CA125-HE4 model also ranked as the best, whereas at 2 years before diagnosis no multimarker model outperformed CA125. CONCLUSIONS: Our study identified and tested different combination of biomarkers using longitudinal multivariable models that outperformed CA125 alone. We showed the potential of multivariable models and candidate biomarkers to increase the detection rate of ovarian cancer.

Technology of the photobiostimulation of the brain's drainage system during sleep for improvement of learning and memory in male mice.

In this study on healthy male mice using confocal imaging of dye spreading in the brain and its further accumulation in the peripheral lymphatics, we demonstrate stronger effects of photobiomodulation (PBM) on the brain's drainage system in sleeping vs. awake animals. Using the Pavlovian instrumental transfer probe and the 2-objects-location test, we found that the 10-day course of PBM during sleep vs. wakefulness promotes improved learning and spatial memory in mice. For the first time, we present the technology for PBM under electroencephalographic (EEG) control that incorporates modern state of the art facilities of optoelectronics and biopotential detection and that can be built of relatively cheap and commercially available components. These findings open a new niche in the development of smart technologies for phototherapy of brain diseases during sleep.

Prediction of spontaneous preterm birth in women with previous full dilatation cesarean delivery.

BACKGROUND: A previous term (≥37 weeks' gestation), full-dilatation cesarean delivery is associated with an increased risk for a subsequent spontaneous preterm birth. The mechanism is unknown. We hypothesized that the cesarean delivery scar characteristics and scar position relative to the internal cervical os may compromise cervical function, thereby leading to shortening of the cervical length and spontaneous preterm birth. OBJECTIVE: This study aimed to determine the relationship of cesarean delivery scar characteristics and position, assessed by transvaginal ultrasound, in pregnant women with previous full-dilatation cesarean delivery with the risk of shortening cervical length and spontaneous preterm birth. STUDY DESIGN: This was a single-center, prospective cohort study of singleton pregnant women (14 to 24 weeks' gestation) with a previous term full-dilatation cesarean delivery who attended a high-risk preterm birth surveillance clinic (2017-2021). Women underwent transvaginal ultrasound assessment of cervical length, cesarean delivery scar distance relative to the internal cervical os, and scar niche parameters using a reproducible transvaginal ultrasound technique. Spontaneous preterm birth prophylactic interventions (vaginal cervical cerclage or vaginal progesterone) were offered for short cervical length (≤25 mm) and to women with a history of spontaneous preterm birth or late miscarriage after full-dilatation cesarean delivery. The primary outcome was spontaneous preterm birth; secondary outcomes included short cervical length and a need for prophylactic interventions. A multivariable logistic regression analysis was used to develop multiparameter models that combined cesarean delivery scar parameters, cervical length, history of full-dilatation cesarean delivery, and maternal characteristics. The predictive performance of models was examined using the area under the receiver operating characteristics curve and the detection rate at various fixed false positive rates. The optimal cutoff for cesarean delivery scar distance to best predict a short cervical length and spontaneous preterm birth was analyzed. RESULTS: Cesarean delivery scars were visualized in 90.5% (220/243) of the included women. The spontaneous preterm birth rate was 4.1% (10/243), and 12.8% (31/243) of women developed a short cervical length. A history- (n=4) or ultrasound-indicated (n=19) cervical cerclage was performed in 23 of 243 (9.5%) women; among those, 2 (8.7%) spontaneously delivered prematurely. A multiparameter model based on absolute scar distance from the internal os best predicted spontaneous preterm birth (area under the receiver operating characteristics curve, 0.73; 95% confidence interval, 0.57-0.89; detection rate of 60% for a fixed 25% false positive rate). Models based on the relative anatomic position of the cesarean delivery scar to the internal os and the cesarean delivery scar position with niche parameters (length, depth, and width) best predicted the development of a short cervical length (area under the receiver operating characteristics curve, 0.79 [95% confidence interval, 0.71-0.87]; and 0.81 [95% confidence interval, 0.73-0.89], respectively; detection rate of 73% at a fixed 25% false positive rate). Spontaneous preterm birth was significantly more likely when the cesarean delivery scar was <5.0 mm above or below the internal os (adjusted odds ratio, 6.87; 95% confidence interval, 1.34-58; P =.035). CONCLUSION: In pregnancies following a full-dilatation cesarean delivery, cesarean delivery scar characteristics and distance from the internal os identified women who were at risk for spontaneous preterm birth and developing short cervical length. Overall, the spontaneous preterm birth rate was low, but it was significantly increased among women with a scar located <5.0 mm above or below the internal cervical os. Shortening of cervical length was strongly associated with a low scar position. Our novel findings indicate that a low cesarean delivery scar can compromise the functional integrity of the internal cervical os, leading to cervical shortening and/or spontaneous preterm birth. Assessment of the cesarean delivery scar characteristics and position seem to have use in preterm birth clinical surveillance among women with a previous, full-dilatation cesarean delivery and could better identify women who would benefit from prophylactic interventions.

Situation-Based Neuromorphic Memory in Spiking Neuron-Astrocyte Network.

Mammalian brains operate in very special surroundings: to survive they have to react quickly and effectively to the pool of stimuli patterns previously recognized as danger. Many learning tasks often encountered by living organisms involve a specific set-up centered around a relatively small set of patterns presented in a particular environment. For example, at a party, people recognize friends immediately, without deep analysis, just by seeing a fragment of their clothes. This set-up with reduced "ontology" is referred to as a "situation." Situations are usually local in space and time. In this work, we propose that neuron-astrocyte networks provide a network topology that is effectively adapted to accommodate situation-based memory. In order to illustrate this, we numerically simulate and analyze a well-established model of a neuron-astrocyte network, which is subjected to stimuli conforming to the situation-driven environment. Three pools of stimuli patterns are considered: external patterns, patterns from the situation associative pool regularly presented to the network and learned by the network, and patterns already learned and remembered by astrocytes. Patterns from the external world are added to and removed from the associative pool. Then, we show that astrocytes are structurally necessary for an effective function in such a learning and testing set-up. To demonstrate this we present a novel neuromorphic computational model for short-term memory implemented by a two-net spiking neural-astrocytic network. Our results show that such a system tested on synthesized data with selective astrocyte-induced modulation of neuronal activity provides an enhancement of retrieval quality in comparison to standard spiking neural networks trained via Hebbian plasticity only. We argue that the proposed set-up may offer a new way to analyze, model, and understand neuromorphic artificial intelligence systems.

An efficient reduced-lattice model of IP3R for probing Ca2+ dynamics.

Numerous cellular processes are regulated by Ca2+ signals, and the endoplasmic reticulum (ER) membrane's inositol triphosphate receptor (IP3R) is critical for modulating intracellular Ca2+ dynamics. The IP3Rs are seen to be clustered in a variety of cell types. The combination of IP3Rs clustering and IP3Rs-mediated Ca2+-induced Ca2+ release results in the hierarchical organization of the Ca2+ signals, which challenges the numerical simulation given the multiple spatial and temporal scales that must be covered. The previous methods rather ignore the spatial feature of IP3Rs or fail to coordinate the conflicts between the real biological relevance and the computational cost. In this work, a general and efficient reduced-lattice model is presented for the simulation of IP3Rs-mediated multiscale Ca2+ dynamics. The model highlights biological details that make up the majority of the calcium events, including IP3Rs clustering and calcium domains, and it reduces the complexity by approximating the minor details. The model's extensibility provides fresh insights into the function of IP3Rs in producing global Ca2+ events and supports the research under more physiological circumstances. Our work contributes to a novel toolkit for modeling multiscale Ca2+ dynamics and advances knowledge of Ca2+ signals.

Using a Recurrent Neural Network To Inform the Use of Prostate-specific Antigen (PSA) and PSA Density for Dynamic Monitoring of the Risk of Prostate Cancer Progression on Active Surveillance.

The global uptake of prostate cancer (PCa) active surveillance (AS) is steadily increasing. While prostate-specific antigen density (PSAD) is an important baseline predictor of PCa progression on AS, there is a scarcity of recommendations on its use in follow-up. In particular, the best way of measuring PSAD is unclear. One approach would be to use the baseline gland volume (BGV) as a denominator in all calculations throughout AS (nonadaptive PSAD, PSADNA), while another would be to remeasure gland volume at each new magnetic resonance imaging scan (adaptive PSAD, PSADA). In addition, little is known about the predictive value of serial PSAD in comparison to PSA. We applied a long short-term memory recurrent neural network to an AS cohort of 332 patients and found that serial PSADNA significantly outperformed both PSADA and PSA for follow-up prediction of PCa progression because of its high sensitivity. Importantly, while PSADNA was superior in patients with smaller glands (BGV ≤55 ml), serial PSA was better in men with larger prostates of >55 ml. Patient summary: Repeat measurements of prostate-specific antigen (PSA) and PSA density (PSAD) are the mainstay of active surveillance in prostate cancer. Our study suggests that in patients with a prostate gland of 55 ml or smaller, PSAD measurements are a better predictor of tumour progression, whereas men with a larger gland may benefit more from PSA monitoring.

Integrated Placental Modelling of Histology with Gene Expression to Identify Functional Impact on Fetal Growth.

Fetal growth restriction (FGR) is a leading cause of perinatal morbidity and mortality. Altered placental formation and functional capacity are major contributors to FGR pathogenesis. Relating placental structure to function across the placenta in healthy and FGR pregnancies remains largely unexplored but could improve understanding of placental diseases. We investigated integration of these parameters spatially in the term human placenta using predictive modelling. Systematic sampling was able to overcome heterogeneity in placental morphological and molecular features. Defects in villous development, elevated fibrosis, and reduced expression of growth and functional marker genes (IGF2, VEGA, SLC38A1, and SLC2A3) were seen in age-matched term FGR versus healthy control placentas. Characteristic histopathological changes with specific accompanying molecular signatures could be integrated through computational modelling to predict if the placenta came from a healthy or FGR pregnancy. Our findings yield new insights into the spatial relationship between placental structure and function and the etiology of FGR.

Time series radiomics for the prediction of prostate cancer progression in patients on active surveillance.

Serial MRI is an essential assessment tool in prostate cancer (PCa) patients enrolled on active surveillance (AS). However, it has only moderate sensitivity for predicting histopathological tumour progression at follow-up, which is in part due to the subjective nature of its clinical reporting and variation among centres and readers. In this study, we used a long short-term memory (LSTM) recurrent neural network (RNN) to develop a time series radiomics (TSR) predictive model that analysed longitudinal changes in tumour-derived radiomic features across 297 scans from 76 AS patients, 28 with histopathological PCa progression and 48 with stable disease. Using leave-one-out cross-validation (LOOCV), we found that an LSTM-based model combining TSR and serial PSA density (AUC 0.86 [95% CI: 0.78-0.94]) significantly outperformed a model combining conventional delta-radiomics and delta-PSA density (0.75 [0.64-0.87]; p = 0.048) and achieved comparable performance to expert-performed serial MRI analysis using the Prostate Cancer Radiologic Estimation of Change in Sequential Evaluation (PRECISE) scoring system (0.84 [0.76-0.93]; p = 0.710). The proposed TSR framework, therefore, offers a feasible quantitative tool for standardising serial MRI assessment in PCa AS. It also presents a novel methodological approach to serial image analysis that can be used to support clinical decision-making in multiple scenarios, from continuous disease monitoring to treatment response evaluation. KEY POINTS: •LSTM RNN can be used to predict the outcome of PCa AS using time series changes in tumour-derived radiomic features and PSA density. •Using all available TSR features and serial PSA density yields a significantly better predictive performance compared to using just two time points within the delta-radiomics framework. •The concept of TSR can be applied to other clinical scenarios involving serial imaging, setting out a new field in AI-driven radiology research.

Serum biomarker-based early detection of pancreatic ductal adenocarcinomas with ensemble learning.

BACKGROUND: Earlier detection of pancreatic ductal adenocarcinoma (PDAC) is key to improving patient outcomes, as it is mostly detected at advanced stages which are associated with poor survival. Developing non-invasive blood tests for early detection would be an important breakthrough. METHODS: The primary objective of the work presented here is to use a dataset that is prospectively collected, to quantify a set of cancer-associated proteins and construct multi-marker models with the capacity to predict PDAC years before diagnosis. The data used is part of a nested case-control study within the UK Collaborative Trial of Ovarian Cancer Screening and is comprised of 218 samples, collected from a total of 143 post-menopausal women who were diagnosed with pancreatic cancer within 70 months after sample collection, and 249 matched non-cancer controls. We develop a stacked ensemble modelling technique to achieve robustness in predictions and, therefore, improve performance in newly collected datasets. RESULTS: Here we show that with ensemble learning we can predict PDAC status with an AUC of 0.91 (95% CI 0.75-1.0), sensitivity of 92% (95% CI 0.54-1.0) at 90% specificity, up to 1 year prior to diagnosis, and at an AUC of 0.85 (95% CI 0.74-0.93) up to 2 years prior to diagnosis (sensitivity of 61%, 95% CI 0.17-0.83, at 90% specificity). CONCLUSIONS: The ensemble modelling strategy explored here outperforms considerably biomarker combinations cited in the literature. Further developments in the selection of classifiers balancing performance and heterogeneity should further enhance the predictive capacity of the method.

Pancreatic cancers are most frequently detected at an advanced stage. This limits treatment options and contributes to the dismal survival rates currently recorded. The development of new tests that could improve detection of early-stage disease is fundamental to improve outcomes. Here, we use advanced data analysis techniques to devise an early detection test for pancreatic cancer. We use data on markers in the blood from people enrolled on a screening trial. Our test correctly identifies as positive for pancreatic cancer 91% of the time up to 1 year prior to diagnosis, and 78% of the time up to 2 years prior to diagnosis. These results surpass previously reported tests and should encourage further evaluation of the test in different populations, to see whether it should be adopted in the clinic.

An Account of Models of Molecular Circuits for Associative Learning with Reinforcement Effect and Forced Dissociation.

The development of synthetic biology has enabled massive progress in biotechnology and in approaching research questions from a brand-new perspective. In particular, the design and study of gene regulatory networks in vitro, in vivo, and in silico have played an increasingly indispensable role in understanding and controlling biological phenomena. Among them, it is of great interest to understand how associative learning is formed at the molecular circuit level. Mathematical models are increasingly used to predict the behaviours of molecular circuits. Fernando's model, which is one of the first works in this line of research using the Hill equation, attempted to design a synthetic circuit that mimics Hebbian learning in a neural network architecture. In this article, we carry out indepth computational analysis of the model and demonstrate that the reinforcement effect can be achieved by choosing the proper parameter values. We also construct a novel circuit that can demonstrate forced dissociation, which was not observed in Fernando's model. Our work can be readily used as reference for synthetic biologists who consider implementing circuits of this kind in biological systems.

A mathematical model of in vitro hepatocellular cholesterol and lipoprotein metabolism for hyperlipidemia therapy.

Cardiovascular diseases associated with high cholesterol (hypercholesterolemia) and low-density lipoproteins (LDL) levels are significant contributors to total mortality in developing and developed countries. Mathematical modeling of LDL metabolism is an important step in the development of drugs for hypercholesterolemia. The aim of this work was to develop and to analyze an integrated mathematical model of cholesterol metabolism in liver cells and its interaction with two types of drugs, statins and PCSK9 inhibitors. The model consisted of 21 ordinary differential equations (ODE) describing cholesterol biosynthesis and lipoprotein endocytosis in liver cells in vitro. The model was tested for its ability to mimic known biochemical effects of familial hypercholesterolemia, statin therapy, and PCSK9 inhibitors. The model qualitatively reproduced the well-known biology of cholesterol regulation, which confirms its potential for minimizing cellular research in initial testing of new drugs for cardiology.

Unraveling Aß-Mediated Multi-Pathway Calcium Dynamics in Astrocytes: Implications for Alzheimer's Disease Treatment From Simulations.

The accumulation of amyloid ß peptide (Aß) in the brain is hypothesized to be the major factor driving Alzheimer's disease (AD) pathogenesis. Mounting evidence suggests that astrocytes are the primary target of Aß neurotoxicity. Aß is known to interfere with multiple calcium fluxes, thus disrupting the calcium homeostasis regulation of astrocytes, which are likely to produce calcium oscillations. Ca2+ dyshomeostasis has been observed to precede the appearance of clinical symptoms of AD; however, it is experimentally very difficult to investigate the interactions of many mechanisms. Given that Ca2+ disruption is ubiquitously involved in AD progression, it is likely that focusing on Ca2+ dysregulation may serve as a potential therapeutic approach to preventing or treating AD, while current hypotheses concerning AD have so far failed to yield curable therapies. For this purpose, we derive and investigate a concise mathematical model for Aß-mediated multi-pathway astrocytic intracellular Ca2+ dynamics. This model accounts for how Aß affects various fluxes contributions through voltage-gated calcium channels, Aß-formed channels and ryanodine receptors. Bifurcation analysis of Aß level, which reflected the corresponding progression of the disease, revealed that Aß significantly induced the increasing [Ca2+] i and frequency of calcium oscillations. The influence of inositol 1,4,5-trisphosphate production (IP3) is also investigated in the presence of Aß as well as the impact of changes in resting membrane potential. In turn, the Ca2+ flux can be considerably changed by exerting specific interventions, such as ion channel blockers or receptor antagonists. By doing so, a "combination therapy" targeting multiple pathways simultaneously has finally been demonstrated to be more effective. This study helps to better understand the effect of Aß, and our findings provide new insight into the treatment of AD.

Parenclitic and Synolytic Networks Revisited.

Parenclitic networks provide a powerful and relatively new way to coerce multidimensional data into a graph form, enabling the application of graph theory to evaluate features. Different algorithms have been published for constructing parenclitic networks, leading to the question-which algorithm should be chosen? Initially, it was suggested to calculate the weight of an edge between two nodes of the network as a deviation from a linear regression, calculated for a dependence of one of these features on the other. This method works well, but not when features do not have a linear relationship. To overcome this, it was suggested to calculate edge weights as the distance from the area of most probable values by using a kernel density estimation. In these two approaches only one class (typically controls or healthy population) is used to construct a model. To take account of a second class, we have introduced synolytic networks, using a boundary between two classes on the feature-feature plane to estimate the weight of the edge between these features. Common to all these approaches is that topological indices can be used to evaluate the structure represented by the graphs. To compare these network approaches alongside more traditional machine-learning algorithms, we performed a substantial analysis using both synthetic data with a priori known structure and publicly available datasets used for the benchmarking of ML-algorithms. Such a comparison has shown that the main advantage of parenclitic and synolytic networks is their resistance to over-fitting (occurring when the number of features is greater than the number of subjects) compared to other ML approaches. Secondly, the capability to visualise data in a structured form, even when this structure is not a priori available allows for visual inspection and the application of well-established graph theory to their interpretation/application, eliminating the "black-box" nature of other ML approaches.

Modeling Working Memory in a Spiking Neuron Network Accompanied by Astrocytes.

We propose a novel biologically plausible computational model of working memory (WM) implemented by a spiking neuron network (SNN) interacting with a network of astrocytes. The SNN is modeled by synaptically coupled Izhikevich neurons with a non-specific architecture connection topology. Astrocytes generating calcium signals are connected by local gap junction diffusive couplings and interact with neurons via chemicals diffused in the extracellular space. Calcium elevations occur in response to the increased concentration of the neurotransmitter released by spiking neurons when a group of them fire coherently. In turn, gliotransmitters are released by activated astrocytes modulating the strength of the synaptic connections in the corresponding neuronal group. Input information is encoded as two-dimensional patterns of short applied current pulses stimulating neurons. The output is taken from frequencies of transient discharges of corresponding neurons. We show how a set of information patterns with quite significant overlapping areas can be uploaded into the neuron-astrocyte network and stored for several seconds. Information retrieval is organized by the application of a cue pattern representing one from the memory set distorted by noise. We found that successful retrieval with the level of the correlation between the recalled pattern and ideal pattern exceeding 90% is possible for the multi-item WM task. Having analyzed the dynamical mechanism of WM formation, we discovered that astrocytes operating at a time scale of a dozen of seconds can successfully store traces of neuronal activations corresponding to information patterns. In the retrieval stage, the astrocytic network selectively modulates synaptic connections in the SNN leading to successful recall. Information and dynamical characteristics of the proposed WM model agrees with classical concepts and other WM models.

Estimating integrated information in bidirectional neuron-astrocyte communication.

There is growing evidence that suggests the importance of astrocytes as elements for neural information processing through the modulation of synaptic transmission. A key aspect of this problem is understanding the impact of astrocytes in the information carried by compound events in neurons across time. In this paper, we investigate how the astrocytes participate in the information integrated by individual neurons in an ensemble through the measurement of "integrated information." We propose a computational model that considers bidirectional communication between astrocytes and neurons through glutamate-induced calcium signaling. Our model highlights the role of astrocytes in information processing through dynamical coordination. Our findings suggest that the astrocytic feedback promotes synergetic influences in the neural communication, which is maximized when there is a balance between excess correlation and spontaneous spiking activity. The results were further linked with additional measures such as net synergy and mutual information. This result reinforces the idea that astrocytes have integrative properties in communication among neurons.

Impact of modular mitochondrial epistatic interactions on the evolution of human subpopulations.

Investigation of human mitochondrial (mt) genome variation has been shown to provide insights to the human history and natural selection. By analyzing 24,167 human mt-genome samples, collected for five continents, we have developed a co-mutation network model to investigate characteristic human evolutionary patterns. The analysis highlighted richer co-mutating regions of the mt-genome, suggesting the presence of epistasis. Specifically, a large portion of COX genes was found to co-mutate in Asian and American populations, whereas, in African, European, and Oceanic populations, there was greater co-mutation bias in hypervariable regions. Interestingly, this study demonstrated hierarchical modularity as a crucial agent for these co-mutation networks. More profoundly, our ancestry-based co-mutation module analyses showed that mutations cluster preferentially in known mitochondrial haplogroups. Contemporary human mt-genome nucleotides most closely resembled the ancestral state, and very few of them were found to be ancestral-variants. Overall, these results demonstrated that subpopulation-based biases may favor mitochondrial gene specific epistasis.

Integrated Information in the Spiking-Bursting Stochastic Model.

Integrated information has been recently suggested as a possible measure to identify a necessary condition for a system to display conscious features. Recently, we have shown that astrocytes contribute to the generation of integrated information through the complex behavior of neuron-astrocyte networks. Still, it remained unclear which underlying mechanisms governing the complex behavior of a neuron-astrocyte network are essential to generating positive integrated information. This study presents an analytic consideration of this question based on exact and asymptotic expressions for integrated information in terms of exactly known probability distributions for a reduced mathematical model (discrete-time, discrete-state stochastic model) reflecting the main features of the "spiking-bursting" dynamics of a neuron-astrocyte network. The analysis was performed in terms of the empirical "whole minus sum" version of integrated information in comparison to the "decoder based" version. The "whole minus sum" information may change sign, and an interpretation of this transition in terms of "net synergy" is available in the literature. This motivated our particular interest in the sign of the "whole minus sum" information in our analytical considerations. The behaviors of the "whole minus sum" and "decoder based" information measures are found to bear a lot of similarity-they have mutual asymptotic convergence as time-uncorrelated activity increases, and the sign transition of the "whole minus sum" information is associated with a rapid growth in the "decoder based" information. The study aims at creating a theoretical framework for using the spiking-bursting model as an analytically tractable reference point for applying integrated information concepts to systems exhibiting similar bursting behavior. The model can also be of interest as a new discrete-state test bench for different formulations of integrated information.

Brain aging and garbage cleaning : Modelling the role of sleep, glymphatic system, and microglia senescence in the propagation of inflammaging.

Brain aging is a complex process involving many functions of our body and described by the interplay of a sleep pattern and changes in the metabolic waste concentration regulated by the microglial function and the glymphatic system. We review the existing modelling approaches to this topic and derive a novel mathematical model to describe the crosstalk between these components within the conceptual framework of inflammaging. Analysis of the model gives insight into the dynamics of garbage concentration and linked microglial senescence process resulting from a normal or disrupted sleep pattern, hence, explaining an underlying mechanism behind healthy or unhealthy brain aging. The model incorporates accumulation and elimination of garbage, induction of glial activation by garbage, and glial senescence by over-activation, as well as the production of pro-inflammatory molecules by their senescence-associated secretory phenotype (SASP). Assuming that insufficient sleep leads to the increase of garbage concentration and promotes senescence, the model predicts that if the accumulation of senescent glia overcomes an inflammaging threshold, further progression of senescence becomes unstoppable even if a normal sleep pattern is restored. Inverting this process by "rejuvenating the brain" is only possible via a reset of concentration of senescent glia below this threshold. Our model approach enables analysis of space-time dynamics of senescence, and in this way, we show that heterogeneous patterns of inflammation will accelerate the propagation of senescence profile through a network, confirming a negative effect of heterogeneity.