Inhibition of Neuronal Necroptosis Mediated by RIPK1 Provides Neuroprotective Effects on Hypoxia and Ischemia In Vitro and In Vivo.

Ischemic brain injury is a widespread pathological condition, the main components of which are a deficiency of oxygen and energy substrates. In recent years, a number of new forms of cell death, including necroptosis, have been described. In necroptosis, a cascade of interactions between the kinases RIPK1 and RIPK3 and the MLKL protein leads to the formation of a specialized death complex called the necrosome, which triggers MLKL-mediated destruction of the cell membrane and necroptotic cell death. Necroptosis probably plays an important role in the development of ischemia/reperfusion injury and can be considered as a potential target for finding methods to correct the disruption of neural networks in ischemic damage. In the present study, we demonstrated that blockade of RIPK1 kinase by Necrostatin-1 preserved the viability of cells in primary hippocampal cultures in an in vitro model of glucose deprivation. The effect of RIPK1 blockade on the bioelectrical and metabolic calcium activity of neuron-glial networks in vitro using calcium imaging and multi-electrode arrays was assessed for the first time. RIPK1 blockade was shown to partially preserve both calcium and bioelectric activity of neuron-glial networks under ischemic factors. However, it should be noted that RIPK1 blockade does not preserve the network parameters of the collective calcium dynamics of neuron-glial networks, despite the maintenance of network bioelectrical activity (the number of bursts and the number of spikes in the bursts). To confirm the data obtained in vitro, we studied the effect of RIPK1 blockade on the resistance of small laboratory animals to in vivo modeling of hypoxia and cerebral ischemia. The use of Necrostatin-1 increases the survival rate of C57BL mice in modeling both acute hypobaric hypoxia and ischemic brain damage.

Novel Algorithm of Network Calcium Dynamics Analysis for Studying the Role of Astrocytes in Neuronal Activity in Alzheimer's Disease Models.

Accumulated experimental data strongly suggest that astrocytes play an important role in the pathogenesis of neurodegeneration, including Alzheimer's disease (AD). The effect of astrocytes on the calcium activity of neuron-astroglia networks in AD modelling was the object of the present study. We have expanded and improved our approach's capabilities to analyze calcium activity. We have developed a novel algorithm to construct dynamic directed graphs of both astrocytic and neuronal networks. The proposed algorithm allows us not only to identify functional relationships between cells and determine the presence of network activity, but also to characterize the spread of the calcium signal from cell to cell. Our study showed that Alzheimer's astrocytes can change the functional pattern of the calcium activity of healthy nerve cells. When healthy nerve cells were cocultivated with astrocytes treated with Aß42, activation of calcium signaling was found. When healthy nerve cells were cocultivated with 5xFAD astrocytes, inhibition of calcium signaling was observed. In this regard, it seems relevant to further study astrocytic-neuronal interactions as an important factor in the regulation of the functional activity of brain cells during neurodegenerative processes. The approach to the analysis of streaming imaging data developed by the authors is a promising tool for studying the collective calcium dynamics of nerve cells.

Synchronization in multiplex models of neuron-glial systems: Small-world topology and inhibitory coupling.

In this work, we investigate the impact of mixed coupling on synchronization in a multiplex oscillatory network. The network mimics the neural-glial systems by incorporating interacting slow ("glial") and fast ("neural") oscillatory layers. Connections between the "glial" elements form a regular periodic structure, in which each element is connected to the eight other neighbor elements, whereas connections among "neural" elements are represented by Watts-Strogatz networks (from regular and small-world to random Erdös-Rényi graph) with a matching mean node degree. We find that the random rewiring toward small-world topology readily yields the dynamics close to that exhibited for a completely random graph, in particular, leading to coarse-graining of dynamics, suppressing multi-stability of synchronization regimes, and the onset of Kuramoto-type synchrony in both layers. The duration of transient dynamics in the system measured by relaxation times is minimized with the increase of random connections in the neural layer, remaining substantial only close to synchronization-desynchronization transitions. "Inhibitory" interactions in the "neural" subnetwork layer undermine synchronization; however, the strong coupling with the "glial" layer overcomes this effect.

Double-Edged Sword of Vitamin D3 Effects on Primary Neuronal Cultures in Hypoxic States.

The use of vitamin D3 along with traditional therapy opens up new prospects for increasing the adaptive capacity of nerve cells to the effects of a wide range of stress factors, including hypoxia-ischemic processes. However, questions about prophylactic and therapeutic doses of vitamin D3 remain controversial. The purpose of our study was to analyze the effects of vitamin D3 at different concentrations on morpho-functional characteristics of neuron-glial networks in hypoxia modeling in vitro. We showed that a single administration of vitamin D3 at a high concentration (1 µM) in a normal state has no significant effect on the cell viability of primary neuronal cultures; however, it has a pronounced modulatory effect on the functional calcium activity of neuron-glial networks and causes destruction of the network response. Under hypoxia, the use of vitamin D3 (1 µM) leads to total cell death of primary neuronal cultures and complete negation of functional neural network activity. In contrast, application of lower concentrations of vitamin D3 (0.01 µM and 0.1 µM) caused a pronounced dose-dependent neuroprotective effect during the studied post-hypoxic period. While the use of vitamin D3 at a concentration of 0.1 µM maintained cell viability, preventive administration of 0.01 µM not only partially preserved the morphological integrity of primary neuronal cells but also maintained the functional structure and activity of neuron-glial networks in cultures. Possible molecular mechanisms of neuroprotective action of vitamin D3 can be associated with the increased expression level of transcription factor HIF-1α and maintaining the relationship between the levels of BDNF and TrkB expression in cells of primary neuronal cultures.

Neuroprotective Effect of Kinase Inhibition in Ischemic Factor Modeling In Vitro.

The contribution of many neuronal kinases to the adaptation of nerve cells to ischemic damage and their effect on functional neural network activity has not yet been studied. The aim of this work is to study the role of the four kinases belonging to different metabolic cascades (SRC, Ikkb, eEF2K, and FLT4) in the adaptive potential of the neuron-glial network for modeling the key factors of ischemic damage. We carried out a comprehensive study on the effects of kinases blockade on the viability and network functional calcium activity of nerve cells under ischemic factor modeling in vitro. Ischemic factor modelling was performed on day 14 of culturing primary hippocampal cells obtained from mouse embryos (E18). The most significant neuroprotective effect was shown in the blockade of FLT4 kinase in the simulation of hypoxia. The studies performed revealed the role of FLT4 in the development of functional dysfunction in cerebrovascular accidents and created new opportunities for the study of this enzyme and its blockers in the formation of new therapeutic strategies.

Quasi-stationary states of game-driven systems: A dynamical approach.

Evolutionary game theory is a framework to formalize the evolution of collectives ("populations") of competing agents that are playing a game and, after every round, update their strategies to maximize individual payoffs. There are two complementary approaches to modeling evolution of player populations. The first addresses essentially finite populations by implementing the apparatus of Markov chains. The second assumes that the populations are infinite and operates with a system of mean-field deterministic differential equations. By using a model of two antagonistic populations, which are playing a game with stationary or periodically varying payoffs, we demonstrate that it exhibits metastable dynamics that is reducible neither to an immediate transition to a fixation (extinction of all but one strategy in a finite-size population) nor to the mean-field picture. In the case of stationary payoffs, this dynamics can be captured with a system of stochastic differential equations and interpreted as a stochastic Hopf bifurcation. In the case of varying payoffs, the metastable dynamics is much more complex than the dynamics of the means.

Signatures of the Consolidated Response of Astrocytes to Ischemic Factors In Vitro.

Whether and under what conditions astrocytes can mount a collective network response has recently become one of the central questions in neurobiology. Here, we address this problem, investigating astrocytic reactions to different biochemical stimuli and ischemic-like conditions in vitro. Identifying an emergent astrocytic network is based on a novel mathematical approach that extracts calcium activity from time-lapse fluorescence imaging and estimates the connectivity of astrocytes. The developed algorithm represents the astrocytic network as an oriented graph in which the nodes correspond to separate astrocytes, and the edges indicate high dynamical correlations between astrocytic events. We demonstrate that ischemic-like conditions decrease network connectivity in primary cultures in vitro, although calcium events persist. Importantly, we found that stimulation under normal conditions with 10 µM ATP increases the number of long-range connections and the degree of corresponding correlations in calcium activity, apart from the frequency of calcium events. This result indicates that astrocytes can form a large functional network in response to certain stimuli. In the post-ischemic interval, the response to ATP stimulation is not manifested, which suggests a deep lesion in functional astrocytic networks. The blockade of Connexin 43 during ischemic modeling preserves the connectivity of astrocytes in the post-hypoxic period.

Transforming Lindblad Equations into Systems of Real-Valued Linear Equations: Performance Optimization and Parallelization of an Algorithm.

With their constantly increasing peak performance and memory capacity, modern supercomputers offer new perspectives on numerical studies of open many-body quantum systems. These systems are often modeled by using Markovian quantum master equations describing the evolution of the system density operators. In this paper, we address master equations of the Lindblad form, which are a popular theoretical tools in quantum optics, cavity quantum electrodynamics, and optomechanics. By using the generalized Gell-Mann matrices as a basis, any Lindblad equation can be transformed into a system of ordinary differential equations with real coefficients. Recently, we presented an implementation of the transformation with the computational complexity, scaling as O(N5logN) for dense Lindbaldians and O(N3logN) for sparse ones. However, infeasible memory costs remains a serious obstacle on the way to large models. Here, we present a parallel cluster-based implementation of the algorithm and demonstrate that it allows us to integrate a sparse Lindbladian model of the dimension N=2000 and a dense random Lindbladian model of the dimension N=200 by using 25 nodes with 64 GB RAM per node.

Spatiotemporal dynamics of distributed synthetic genetic circuits.

We propose and study models of two distributed synthetic gene circuits, toggle-switch and oscillator, each split between two cell strains and coupled via quorum-sensing signals. The distributed toggle switch relies on mutual repression of the two strains, and oscillator is comprised of two strains, one of which acts as an activator for another that in turn acts as a repressor. Distributed toggle switch can exhibit mobile fronts, switching the system from the weaker to the stronger spatially homogeneous state. The circuit can also act as a biosensor, with the switching front dynamics determined by the properties of an external signal. Distributed oscillator system displays another biosensor functionality: oscillations emerge once a small amount of one cell strain appears amid the other, present in abundance. Distribution of synthetic gene circuits among multiple strains allows one to reduce crosstalk among different parts of the overall system and also decrease the energetic burden of the synthetic circuit per cell, which may allow for enhanced functionality and viability of engineered cells.

eXplainable Artificial Intelligence (XAI) in aging clock models.

XAI is a rapidly progressing field of machine learning, aiming to unravel the predictions of complex models. XAI is especially required in sensitive applications, e.g. in health care, when diagnosis, recommendations and treatment choices might rely on the decisions made by artificial intelligence systems. AI approaches have become widely used in aging research as well, in particular, in developing biological clock models and identifying biomarkers of aging and age-related diseases. However, the potential of XAI here awaits to be fully appreciated. We discuss the application of XAI for developing the "aging clocks" and present a comprehensive analysis of the literature categorized by the focus on particular physiological systems.

Map of epigenetic age acceleration: A worldwide analysis.

We present a systematic analysis of epigenetic age acceleration based on by far the largest collection of publicly available DNA methylation data for healthy samples (93 datasets, 23 K samples), focusing on the geographic (25 countries) and ethnic (31 ethnicities) aspects around the world. We employed the most popular epigenetic tools for assessing age acceleration and examined their quality metrics and ability to extrapolate to epigenetic data from different tissue types and age ranges different from the training data of these models. In most cases, the models proved to be inconsistent with each other and showed different signs of age acceleration, with the PhenoAge model tending to systematically underestimate and different versions of the GrimAge model tending to systematically overestimate the age prediction of healthy subjects. Referring to data availability and consistency, most countries and populations are still not represented in GEO, moreover, different datasets use different criteria for determining healthy controls. Because of this, it is difficult to fully isolate the contribution of "geography/environment", "ethnicity" and "healthiness" to epigenetic age acceleration. Among the explored metrics, only the DunedinPACE, which measures aging rate, appears to adequately reflect the standard of living and socioeconomic indicators in countries, although it has a limited application to blood methylation data only. Invariably, by epigenetic age acceleration, males age faster than females in most of the studied countries and populations.

Organ-specific biological clocks: Ageotyping for personalized anti-aging medicine.

Aging is a complex multidimensional, progressive remodeling process affecting multiple organ systems. While many studies have focused on studying aging across multiple organs, assessment of the contribution of individual organs to overall aging processes is a cutting-edge issue. An organ's biological age might influence the aging of other organs, revealing a multiorgan aging network. Recent data demonstrated a similar yet asynchronous inter-organs and inter-individuals progression of aging, thereby providing a foundation to track sources of declining health in old age. The integration of multiple omics with common clinical parameters through artificial intelligence has allowed the building of organ-specific aging clocks, which can predict the development of specific age-related diseases at high resolution. The peculiar individual aging-trajectory, referred to as ageotype, might provide a novel tool for a personalized anti-aging, preventive medicine. Here, we review data relative to biological aging clocks and omics-based data, suggesting different organ-specific aging rates. Additional research on longitudinal data, including young subjects and analyzing sex-related differences, should be encouraged to apply ageotyping analysis for preventive purposes in clinical practice.

Small immunological clocks identified by deep learning and gradient boosting.

Background: The aging process affects all systems of the human body, and the observed increase in inflammatory components affecting the immune system in old age can lead to the development of age-associated diseases and systemic inflammation. Results: We propose a small clock model SImAge based on a limited number of immunological biomarkers. To regress the chronological age from cytokine data, we first use a baseline Elastic Net model, gradient-boosted decision trees models, and several deep neural network architectures. For the full dataset of 46 immunological parameters, DANet, SAINT, FT-Transformer and TabNet models showed the best results for the test dataset. Dimensionality reduction of these models with SHAP values revealed the 10 most age-associated immunological parameters, taken to construct the SImAge small immunological clock. The best result of the SImAge model shown by the FT-Transformer deep neural network model has mean absolute error of 6.94 years and Pearson ρ = 0.939 on the independent test dataset. Explainable artificial intelligence methods allow for explaining the model solution for each individual participant. Conclusions: We developed an approach to construct a model of immunological age based on just 10 immunological parameters, coined SImAge, for which the FT-Transformer deep neural network model had proved to be the best choice. The model shows competitive results compared to the published studies on immunological profiles, and takes a smaller number of features as an input. Neural network architectures outperformed gradient-boosted decision trees, and can be recommended in the further analysis of immunological profiles.

Differential protease content of mast cells and the processing of IL-33 in Alternaria alternata induced allergic airway inflammation in mice.

Background: Recent in vitro studies strongly implicated mast cell-derived proteases as regulators of IL-33 activity by enzymatic cleavage in its central domain. A better understanding of the role of mast cell proteases on IL-33 activity in vivo is needed. We aimed to compare the expression of mast cell proteases in C57BL/6 and BALB/c mice, their role in the cleavage of IL-33 cytokine, and their contribution to allergic airway inflammation. Results: In vitro, full-length IL-33 protein was efficiently degraded by mast cell supernatants of BALB/c mice in contrast to the mast cell supernatants from C57BL/6 mice. RNAseq analysis indicated major differences in the gene expression profiles of bone marrow-derived mast cells from C57BL/6 and BALB/c mice. In Alternaria alternata (Alt) - treated C57BL/6 mice the full-length form of IL-33 was mainly present, while in BALB/c mice, the processed shorter form of IL-33 was more prominent. The observed cleavage pattern of IL-33 was associated with a nearly complete lack of mast cells and their proteases in the lungs of C57BL/6 mice. While most inflammatory cells were similarly increased in Alt-treated C57BL/6 and BALB/c mice, C57BL/6 mice had significantly more eosinophils in the bronchoalveolar lavage fluid and IL-5 protein levels in their lungs than BALB/c mice. Conclusion: Our study demonstrates that lung mast cells differ in number and protease content between the two tested mouse strains and could affect the processing of IL-33 and inflammatory outcome of Alt -induced airway inflammation. We suggest that mast cells and their proteases play a regulatory role in IL-33-induced lung inflammation by limiting its proinflammatory effect via the IL-33/ST2 signaling pathway.

Epigenetics of the far northern Yakutian population.

BACKGROUND: Yakuts are one of the indigenous populations of the subarctic and arctic territories of Siberia characterized by a continental subarctic climate with severe winters, with the regular January average temperature in the regional capital city of Yakutsk dipping below - 40 °C. The epigenetic mechanisms of adaptation to such ecologies and environments and, in particular, epigenetic age acceleration in the local population have not been studied before. RESULTS: This work reports the first epigenetic study of the Yakutian population using whole-blood DNA methylation data, supplemented with the comparison to the residents of Central Russia. Gene set enrichment analysis revealed, among others, geographic region-specific differentially methylated regions associated with adaptation to climatic conditions (water consumption, digestive system regulation), aging processes (actin filament activity, cell fate), and both of them (channel activity, regulation of steroid and corticosteroid hormone secretion). Further, it is demonstrated that the epigenetic age acceleration of the Yakutian representatives is significantly higher than that of Central Russia counterparts. For both geographic regions, we showed that epigenetically males age faster than females, whereas no significant sex differences were found between the regions. CONCLUSIONS: We performed the first study of the epigenetic data of the Yakutia cohort, paying special attention to region-specific features, aging processes, age acceleration, and sex specificity.

Controlling synchronization of gamma oscillations by astrocytic modulation in a model hippocampal neural network.

Recent in vitro and in vivo experiments demonstrate that astrocytes participate in the maintenance of cortical gamma oscillations and recognition memory. However, the mathematical understanding of the underlying dynamical mechanisms remains largely incomplete. Here we investigate how the interplay of slow modulatory astrocytic signaling with fast synaptic transmission controls coherent oscillations in the network of hippocampal interneurons that receive inputs from pyramidal cells. We show that the astrocytic regulation of signal transmission between neurons improves the firing synchrony and extends the region of coherent oscillations in the biologically relevant values of synaptic conductance. Astrocyte-mediated potentiation of inhibitory synaptic transmission markedly enhances the coherence of network oscillations over a broad range of model parameters. Astrocytic regulation of excitatory synaptic input improves the robustness of interneuron network gamma oscillations induced by physiologically relevant excitatory model drive. These findings suggest a mechanism, by which the astrocytes become involved in cognitive function and information processing through modulating fast neural network dynamics.

Role of mitochondrial genetic interactions in determining adaptation to high altitude human population.

Physiological and haplogroup studies performed to understand high-altitude adaptation in humans are limited to individual genes and polymorphic sites. Due to stochastic evolutionary forces, the frequency of a polymorphism is affected by changes in the frequency of a near-by polymorphism on the same DNA sample making them connected in terms of evolution. Here, first, we provide a method to model these mitochondrial polymorphisms as "co-mutation networks" for three high-altitude populations, Tibetan, Ethiopian and Andean. Then, by transforming these co-mutation networks into weighted and undirected gene-gene interaction (GGI) networks, we were able to identify functionally enriched genetic interactions of CYB and CO3 genes in Tibetan and Andean populations, while NADH dehydrogenase genes in the Ethiopian population playing a significant role in high altitude adaptation. These co-mutation based genetic networks provide insights into the role of different set of genes in high-altitude adaptation in human sub-populations.

Disease classification for whole-blood DNA methylation: Meta-analysis, missing values imputation, and XAI.

BACKGROUND: DNA methylation has a significant effect on gene expression and can be associated with various diseases. Meta-analysis of available DNA methylation datasets requires development of a specific workflow for joint data processing. RESULTS: We propose a comprehensive approach of combined DNA methylation datasets to classify controls and patients. The solution includes data harmonization, construction of machine learning classification models, dimensionality reduction of models, imputation of missing values, and explanation of model predictions by explainable artificial intelligence (XAI) algorithms. We show that harmonization can improve classification accuracy by up to 20% when preprocessing methods of the training and test datasets are different. The best accuracy results were obtained with tree ensembles, reaching above 95% for Parkinson's disease. Dimensionality reduction can substantially decrease the number of features, without detriment to the classification accuracy. The best imputation methods achieve almost the same classification accuracy for data with missing values as for the original data. XAI approaches have allowed us to explain model predictions from both populational and individual perspectives. CONCLUSIONS: We propose a methodologically valid and comprehensive approach to the classification of healthy individuals and patients with various diseases based on whole-blood DNA methylation data using Parkinson's disease and schizophrenia as examples. The proposed algorithm works better for the former pathology, characterized by a complex set of symptoms. It allows to solve data harmonization problems for meta-analysis of many different datasets, impute missing values, and build classification models of small dimensionality.

Accelerated epigenetic aging and inflammatory/immunological profile (ipAGE) in patients with chronic kidney disease.

Chronic kidney disease (CKD) is defined by a reduced estimated glomerular filtration rate (eGFR). This failure can be related to a phenotype of accelerated aging. In this work, we considered 76 patients with end-stage renal disease (ESRD) and 83 healthy controls. We concomitantly evaluated for the first time two measures that can be informative of the rate of aging, i.e., whole blood DNA methylation using the Illumina Infinium EPIC array and plasma levels of a selection of inflammatory/immunological proteins using multiplex immunoassays. First of all, we demonstrated accelerated aging in terms of the most common epigenetic age estimators in CKD patients. Moreover, we developed a new clock/predictor of age based on the inflammatory/immunological profile (ipAGE) and identified the inflammatory/immunological biomarkers differentially expressed between cases and controls. IpAGE appeared to be more sensitive than epigenetic clocks in quantifying the accelerated aging phenotype of ESRD patients. Interestingly, we did not find any correlation between the age acceleration evaluated according to the epigenetic clocks and ipAGE in either the ESRD group or the control group. On the whole, our data show a consistent accelerated aging phenotype in ESRD patients, which is better appreciated by quantifying the underlying inflammatory processes (inflammaging) by ipAGE than by using epigenetic clocks.