Search

The field of precision radiation therapy has seen remarkable advancements in both experimental and computational methods. Recent literature has introduced various approaches such as Spatially Fractionated Radiation Therapy (SFRT). This unconventional treatment, demanding high-precision radiotherapy, has shown promising clinical outcomes. A comprehensive computational scheme for SFRT, extrapolated from a case report, is proposed. This framework exhibits exceptional flexibility, accommodating diverse initial conditions (shape, inhomogeneity, etc.) and enabling specific choices for sub-volume selection with administrated higher radiation doses. The approach integrates the standard linear quadratic model and, significantly, considers the activation of the immune system due to radiotherapy. This activation enhances the immune response in comparison to the untreated case. We delve into the distinct roles of the native immune system, immune activation by radiation, and post-radiotherapy immunotherapy, discussing their implications for either complete recovery or disease regrowth. [ABSTRACT FROM AUTHOR]

Background: The aim of the present paper is to construct an emulator of a complex biological system simulator using a machine learning approach. More specifically, the simulator is a patient-specific model that integrates metabolic, nutritional, and lifestyle data to predict the metabolic and inflammatory processes underlying the development of type-2 diabetes in absence of familiarity. Given the very high incidence of type-2 diabetes, the implementation of this predictive model on mobile devices could provide a useful instrument to assess the risk of the disease for aware individuals. The high computational cost of the developed model, being a mixture of agent-based and ordinary differential equations and providing a dynamic multivariate output, makes the simulator executable only on powerful workstations but not on mobile devices. Hence the need to implement an emulator with a reduced computational cost that can be executed on mobile devices to provide real-time self-monitoring. Results: Similarly to our previous work, we propose an emulator based on a machine learning algorithm but here we consider a different approach which turn out to have better performances, indeed in terms of root mean square error we have an improvement of two order magnitude. We tested the proposed emulator on samples containing different number of simulated trajectories, and it turned out that the fitted trajectories are able to predict with high accuracy the entire dynamics of the simulator output variables. We apply the emulator to control the level of inflammation while leveraging on the nutritional input. Conclusion: The proposed emulator can be implemented and executed on mobile health devices to perform quick-and-easy self-monitoring assessments. [ABSTRACT FROM AUTHOR]

Background: Macrophages cover a major role in the immune system, being the most plastic cell yielding several key immune functions. Methods: Here we derived a minimalistic gene regulatory network model for the differentiation of macrophages into the two phenotypes M1 (pro-) and M2 (anti-inflammatory). Results: To test the model, we simulated a large number of such networks as in a statistical ensemble. In other words, to enable the inter-cellular crosstalk required to obtain an immune activation in which the macrophage plays its role, the simulated networks are not taken in isolation but combined with other cellular agents, thus setting up a discrete minimalistic model of the immune system at the microscopic/intracellular (i.e., genetic regulation) and mesoscopic/intercellular scale. Conclusions: We show that within the mesoscopic level description of cellular interaction and cooperation, the gene regulatory logic is coherent and contributes to the overall dynamics of the ensembles that shows, statistically, the expected behaviour. [ABSTRACT FROM AUTHOR]

Duchenne muscular dystrophy (DMD) is a genetic disease that results in the death of affected boys by early adulthood. The genetic defect responsible for DMD has been known for over 25 years, yet at present there is neither cure nor effective treatment for DMD. During early disease onset, the mdx mouse has been validated as an animal model for DMD and use of this model has led to valuable but incomplete insights into the disease process. For example, immune cells are thought to be responsible for a significant portion of muscle cell death in the mdx mouse; however, the role and time course of the immune response in the dystrophic process have not been well described. In this paper we constructed a simple mathematical model to investigate the role of the immune response in muscle degeneration and subsequent regeneration in the mdx mouse model of Duchenne muscular dystrophy. Our model suggests that the immune response contributes substantially to the muscle degeneration and regeneration processes. Furthermore, the analysis of the model predicts that the immune system response oscillates throughout the life of the mice, and the damaged fibers are never completely cleared. [ABSTRACT FROM AUTHOR]

In this article we present evidence for a relationship between chromosome gene loci and the topological properties of the protein-protein interaction network corresponding to the set of genes under consideration. Specifically, for each chromosome of the Saccharomyces cerevisiae genome, the distribution of the intra-chromosome inter-gene distances was analyzed and a positive correlation with the distance among the corresponding proteins of the protein-protein interaction network was found. In order to study this relationship we used concepts based on non-parametric statistics and information theory. We provide statistical evidence that if two genes are closely located, then it is likely that their protein products are closely located in the protein-protein interaction network, or in other words, that they are involved in the same biological process. [ABSTRACT FROM AUTHOR]

Continuous antiretroviral therapy is currently the most effective way to treat HIV infection. Unstructured interruptions are quite common due to side effects and toxicity, among others, and cannot be prevented. Several attempts to structure these interruptions failed due to an increased morbidity compared to continuous treatment. The cause of this failure is poorly understood and often attributed to drug resistance. Here we show that structured treatment interruptions would fail regardless of the emergence of drug resistance. Our computational model of the HIV infection dynamics in lymphoid tissue inside lymph nodes, demonstrates that HIV reservoirs and evasion from immune surveillance themselves are sufficient to cause the failure of structured interruptions. We validate our model with data from a clinical trial and show that it is possible to optimize the schedule of interruptions to perform as well as the continuous treatment in the absence of drug resistance. Our methodology enables studying the problem of treatment optimization without having impact on human beings. We anticipate that it is feasible to steer new clinical trials using computational models. [ABSTRACT FROM AUTHOR]

The time of initiation of antiretroviral therapy in HIV-1 infected patients has a determinant effect on the viral dynamics. The question is, how far can the therapy be delayed? Is sooner always better? We resort to clinical data and to microsimulations to forecast the dynamics of the viral load at therapy interruption after prolonged antiretroviral treatment. A computational model previously evaluated, produces results that are statistically adherent to clinical data. In addition, it allows a finer grain analysis of the impact of the therapy initiation point to the disease course. We find a swift increase of the viral density as a function of the time of initiation of the therapy measured when the therapy is stopped. In particular there is a critical time delay with respect to the infection instant beyond which the therapy does not affect the viral rebound. Initiation of the treatment is beneficial because it can down-regulate the immune activation, hence limiting viral replication and spread. [ABSTRACT FROM AUTHOR]

Abstract: We present a mathematical model to investigate the role of the immune system in the Duchenne muscular dystrophy disease. It is based on the assumption that the immune system contributes to the tissue damage and indeed its interaction with the muscle tissue after an initial endogenous damage can be described as a predator–prey system showing typical oscillations. In this article we investigate the dynamical properties of the system. We find that, for a biologically relevant parameters range, it shows two phase transitions between qualitative different behaviors corresponding to complete recover or to a state where muscle regeneration and degeneration coexist. [Copyright &y& Elsevier]

We have developed a rat brain organotypic culture model, in which tissue slices contain cortex-subventricular zone-striatum regions, to model neuroblast activity in response to in vitro ischemia. Neuroblast activation has been described in terms of two main parameters, proliferation and migration from the subventricular zone into the injured cortex. We observed distinct phases of neuroblast activation as is known to occur after in vivo ischemia. Thus, immediately after oxygen/glucose deprivation (6-24 hours), neuroblasts reduce their proliferative and migratory activity, whereas, at longer time points after the insult (2 to 5 days), they start to proliferate and migrate into the damaged cortex. Antagonism of ionotropic receptors for extracellular ATP during and after the insult unmasks an early activation of neuroblasts in the subventricular zone, which responded with a rapid and intense migration of neuroblasts into the damaged cortex (within 24 hours). The process is further enhanced by elevating the production of the chemoattractant SDf-1α and may also be boosted by blocking the activation of microglia. This organotypic model which we have developed is an excellent in vitro system to study neurogenesis after ischemia and other neurodegenerative diseases. Its application has revealed a SOS response to oxygen/ glucose deprivation, which is inhibited by unfavorable conditions due to the ischemic environment. Finally, experimental quantifications have allowed us to elaborate a mathematical model to describe neuroblast activation and to develop a computer simulation which should have promising applications for the screening of drug candidates for novel therapies of ischemia-related pathologies. [ABSTRACT FROM AUTHOR]

Background: Protein interactions support cell organization and mediate its response to any specific stimulus. Recent technological advances have produced large data-sets that aim at describing the cell interactome. These data are usually presented as graphs where proteins (nodes) are linked by edges to their experimentally determined partners. This representation reveals that protein-protein interaction (PPI) networks, like other kinds of complex networks, are not randomly organized and display properties that are typical of "hierarchical" networks, combining modularity and local clustering to scale free topology. However informative, this representation is static and provides no clue about the dynamic nature of protein interactions inside the cell. Results: To fill this methodological gap, we designed and implemented a computer model that captures the discrete and stochastic nature of protein interactions. In ProtNet, our simplified model, the intracellular space is mapped onto either a two-dimensional or a three-dimensional lattice with each lattice site having a linear size (5 nm) comparable to the diameter of an average globular protein. The protein filled lattice has an occupancy (e.g. 20%) compatible with the estimated crowding of proteins in the cell cytoplasm. Proteins or protein complexes are free to translate and rotate on the lattice that represents a sort of naïve unstructured cell (devoid of compartments). At each time step, molecular entities (proteins or complexes) that happen to be in neighboring cells may interact and form larger complexes or dissociate depending on the interaction rules defined in an experimental protein interaction network. This whole procedure can be seen as a sort of "discrete molecular dynamics" applied to interacting proteins in a cell. We have tested our model by performing different simulations using as interaction rules those derived from an experimental interactome of Saccharomyces cerevisiae (1378 nodes, 2491 edges) and we have compared the dynamics of complex formation in a two and a three dimensional lattice model. Conclusion: ProtNet is a cellular automaton model, where each protein molecule or complex is explicitly represented and where simple interaction rules are applied to populations of discrete particles. This tool can be used to simulate the dynamics of protein interactions in the cell. [ABSTRACT FROM AUTHOR]

We construct a population dynamics model of the competition among immune system cells and generic tumor cells. Then, we apply the theory of optimal control to find the optimal schedule of injection of autologous dendritic cells used as immunotherapeutic agent. The optimization method works for a general ODE system and can be applied to find the optimal schedule in a variety of medical treatments that have been described by a mathematical model. [ABSTRACT FROM AUTHOR]

Abstract: Vaccination protocols designed to elicit anti-cancer immune responses have, many times, failed in producing tumor eradication and in prolonging patient survival. Usually in cancer vaccination, epitopes from one organism are included in the genome or linked with some protein of another in the hope that the immunogenic properties of the latter will boost an immune response to the former. However, recent results have demonstrated that injections of two different vectors encoding the same recombinant antigen generate high levels of specific immunity. Systematic comparison of the efficacy of different vaccination protocols has been hampered by technical limitations, and clear evidence that the use of multiple vectors has advantages over single carrier injections is lacking. We used a computational model to investigate the dynamics of the immune response to different anti-cancer vaccines based on randomly generated antigen/carrier compounds. The computer model was adapted for simulations to this new area in immunology research and carefully validated to the purpose. As a matter of fact, it reproduces a relevant number of experimental observations. The model shows that when priming and boosting with the same construct, competition rather than cooperation develops amongst T cell clones of different specificities. Moreover, from the simulations, it appears that the sequential use of multiple carriers may generate more robust anti-tumor immune responses and may lead to effective tumor eradication in a higher percentage of cases. Our results provide a rational background for the design of novel strategies for the achievement of immune control of cancer. [Copyright &y& Elsevier]

Abstract: Agent-based models, an emerging paradigm of simulation of complex systems, appear very suitable to parallel processing. However, during the parallelization of a simulator of financial markets, we found that some features of these codes highlight non-trivial issues of the present hardware/software platforms for parallel processing. Here we present the results of a series of tests, on different platforms, of simplified codes that reproduce such problems and can be used as a starting point in the search of a possible solution. [Copyright &y& Elsevier]

Mathematical and computational models are designed to improve our understanding of biological phenomena, to confirm/reject hypotheses, and to find points of intervention by altering the behavior of the studied systems. Here we describe the role of mathematical/computational models of the immune system. In particular, we analyze some examples of how mathematical modeling can contribute to finding optimal vaccination strategies. Indeed, computational modeling offers an intriguing opportunity from the theoretical point of view, and it will be of interest for clinically oriented investigators who wish to find optimal therapeutic strategies and for pharmaceutical industries that want to produce effective and successful drugs. [ABSTRACT FROM AUTHOR]

Financial forecasting is a difficult task due to the intrinsic complexity of the financial system. A simplified approach in forecasting is given by “black box” methods like neural networks that assume little about the structure of the economy. In the present paper we relate our experience using neural nets as financial time series forecast method. In particular we show that a neural net able to forecast the sign of the price increments with a success rate slightly above 50% can be found. Target series are the daily closing price of different assets and indexes during the period from about January 1990 to February 2000. [ABSTRACT FROM AUTHOR]

We present a minimal regulatory model for the dynamics of the humoral immune response of two lymphocytes populations (B and T helper) interacting with a specific antigen pool (bacterium). Stability analysis reveals that the system accounts for the occurrence of multiple steady states in the absence as well as in the presence of the antigen population. The model exhibits (i) a state of immune memory, (ii) one state with high antigen and low helper concentration (disease), and (iii) one state with low antigen and high helper concentration (tolerance). The latter state allows oscillatory behavior. Injection of high antigen doses as well as minimal changes of structural parameters provoke the system to jump from the state of disease to the state of tolerance. This is reminiscent of therapies where the patient is treated with allergen, immuno-suppressants or drugs. [ABSTRACT FROM AUTHOR]

Focuses on a study which proposed a semi-automatic procedure to compute the keratoconus index (KI) from ultrasound biomicroscope images. Cubic spline interpolation of corneal contours; Details on compute thickness and KI; Results of the study.

We describe a new model to simulate the dynamic interactions between market price and the decisions of two different kind of traders. They possess spatial mobility allowing to group together to form coalitions. Each coalition follows a strategy chosen from a proportional voting "dominated" by a leader's decision. The interplay of both kind of agents gives rise to complex price dynamics that is consistent with the main stylized facts of financial time series. The present model incorporates many features of other known models and is meant to be the first step toward the construction of an agent-based model that uses more realistic markets rules, strategies, and information structures. [ABSTRACT FROM AUTHOR]

Medical digital twins are computational models of human biology relevant to a given medical condition, which are tailored to an individual patient, thereby predicting the course of disease and individualized treatments, an important goal of personalized medicine. The immune system, which has a central role in many diseases, is highly heterogeneous between individuals, and thus poses a major challenge for this technology. In February 2023, an international group of experts convened for two days to discuss these challenges related to immune digital twins. The group consisted of clinicians, immunologists, biologists, and mathematical modelers, representative of the interdisciplinary nature of medical digital twin development. A video recording of the entire event is available. This paper presents a synopsis of the discussions, brief descriptions of ongoing digital twin projects at different stages of progress. It also proposes a 5-year action plan for further developing this technology. The main recommendations are to identify and pursue a small number of promising use cases, to develop stimulation-specific assays of immune function in a clinical setting, and to develop a database of existing computational immune models, as well as advanced modeling technology and infrastructure. [ABSTRACT FROM AUTHOR]

Background: Neutrophils are one of the key players in the human innate immune system (HIIS). In the event of an insult where the body is exposed to inflammation triggering moieties (ITMs), neutrophils are mobilized towards the site of insult and antagonize the inflammation. If the inflammation is cleared, neutrophils go into a programmed death called apoptosis. However, if the insult is intense or persistent, neutrophils take on a violent death pathway called necrosis, which involves the rupture of their cytoplasmic content into the surrounding tissue that causes local tissue damage, thus further aggravating inflammation. This seemingly paradoxical phenomenon fuels the inflammatory process by triggering the recruitment of additional neutrophils to the site of inflammation, aimed to contribute to the complete neutralization of severe inflammation. This delicate balance between the cost and benefit of the neutrophils’ choice of death pathway has been optimized during the evolution of the innate immune system. The goal of our work is to understand how the tradeoff between the cost and benefit of the different death pathways of neutrophils, in response to various levels of insults, has been optimized over evolutionary time by using the concepts of evolutionary game theory. Results: We show that by using evolutionary game theory, we are able to formulate a game that predicts the percentage of necrosis and apoptosis when exposed to various levels of insults. Conclusion: By adopting an evolutionary perspective, we identify the driving mechanisms leading to the delicate balance between apoptosis and necrosis in neutrophils’ cell death in response to different insults. Using our simple model, we verify that indeed, the global cost of remaining ITMs is the driving mechanism that reproduces the percentage of necrosis and apoptosis observed in data and neutrophils need sufficient information of the overall inflammation to be able to pick a death pathway that presumably increases the survival of the organism. [ABSTRACT FROM AUTHOR]

It is coming nowadays more clear that in order to obtain a unified description of the different mechanisms governing the behavior and causality relations among the various parts of a living system, the development of comprehensive computational and mathematical models at different space and time scales is required. This is one of the most formidable challenges of modern biology characterized by the availability of huge amount of high throughput measurements. In this paper we draw attention to the importance of multiscale modeling in the framework of studies of biological systems in general and of the immune system in particular. [ABSTRACT FROM AUTHOR]

One of the greatest challenges in biomedicine is to get a unified view of observations made from the molecular up to the organism scale. Towards this goal, multiscale models have been highly instrumental in contexts such as the cardiovascular field, angiogenesis, neurosciences and tumour biology. More recently, such models are becoming an increasingly important resource to address immunological questions as well. Systematic mining of the literature in multiscale modelling led us to identify three main fields of immunological applications: host-virus interactions, inflammatory diseases and their treatment and development of multiscale simulation platforms for immunological research and for educational purposes. Here, we review the current developments in these directions, which illustrate that multiscale models can consistently integrate immunological data generated at several scales, and can be used to describe and optimize therapeutic treatments of complex immune diseases. [ABSTRACT FROM AUTHOR]

Background: Cell organization is governed and maintained via specific interactions among its constituent macromolecules. Comparison of the experimentally determined protein interaction networks in different model organisms has revealed little conservation of the specific edges linking ortholog proteins. Nevertheless, some topological characteristics of the graphs representing the networks - namely non-random degree distribution and high clustering coefficient - are shared by networks of distantly related organisms. Here we investigate the role of the topological features of the protein interaction network in promoting cell organization. Methods: We have used a stochastic model, dubbed ProtNet representing a computer stylized cell to answer questions about the dynamic consequences of the topological properties of the static graphs representing protein interaction networks. Results: By using a novel metrics of cell organization, we show that natural networks, differently from random networks, can promote cell self-organization. Furthermore the ensemble of protein complexes that forms in pseudocells, which self-organize according to the interaction rules of natural networks, are more robust to perturbations. Conclusions: The analysis of the dynamic properties of networks with a variety of topological characteristics lead us to conclude that self organization is a consequence of the high clustering coefficient, whereas the scale free degree distribution has little influence on this property. [ABSTRACT FROM AUTHOR]

Background: Despite progress in conventional cancer therapies, cancer is still one of the leading causes of death in industrial nations. Therefore, an urgent need of progress in fighting cancer remains. A promising alternative to conventional methods is immune therapy. This relies on the fact that low-immunogenic tumours can be eradicated if an immune response against them is induced. Peptide vaccination is carried out by injecting tumour peptides into a patient to trigger a specific immune response against the tumour in its entirety. However, peptide vaccination is a highly complicated treatment and currently many factors like the optimal number of epitopes are not known precisely. Therefore, it is necessary to evaluate how certain parameters influence the therapy. Results: We present the VaccImm Server that allows users to simulate peptide vaccination in cancer therapy. It uses an agent-based model that simulates peptide vaccination by explicitly modelling the involved cells (immune system and cancer) as well as molecules (antibodies, antigens and semiochemicals). As a new feature, our model uses real amino acid sequences to represent molecular binding sites of relevant immune cells. The model is used to generate detailed statistics of the population sizes and states of the single cell types over time. This makes the VaccImm web server well suited to examine the parameter space of peptide vaccination in silico. VaccImm is publicly available without registration on the web at http://bioinformatics.charite.de/vaccimm; all major browsers are supported. Conclusions: The VaccImm Server provides a convenient way to analyze properties of peptide vaccination in cancer therapy. Using the server, we could gain interesting insights into peptide vaccination that reveal the complex and patient-specific nature of peptide vaccination. [ABSTRACT FROM AUTHOR]

Background: The immune contribution to cancer progression is complex and difficult to characterize. For example in tumors, immune gene expression is detected from the combination of normal, tumor and immune cells in the tumor microenvironment. Profiling the immune component of tumors may facilitate the characterization of the poorly understood roles immunity plays in cancer progression. However, the current approaches to analyze the immune component of a tumor rely on incomplete identification of immune factors. Methods: To facilitate a more comprehensive approach, we created a ranked immunological relevance score for all human genes, developed using a novel strategy that combines text mining and information theory. We used this score to assign an immunological grade to gene expression profiles, and thereby quantify the immunological component of tumors. This immunological relevance score was benchmarked against existing manually curated immune resources as well as high-throughput studies. To further characterize immunological relevance for genes, the relevance score was charted against both the human interactome and cancer information, forming an expanded interactome landscape of tumor immunity. We applied this approach to expression profiles in melanomas, thus identifying and grading their immunological components, followed by identification of their associated protein interactions. Results: The power of this strategy was demonstrated by the observation of early activation of the adaptive immune response and the diversity of the immune component during melanoma progression. Furthermore, the genome-wide immunological relevance score classified melanoma patient groups, whose immunological grade correlated with clinical features, such as immune phenotypes and survival. Conclusions: The assignment of a ranked immunological relevance score to all human genes extends the content of existing immune gene resources and enriches our understanding of immune involvement in complex biological networks. The application of this approach to tumor immunity represents an automated systems strategy that quantifies the immunological component in complex disease. In so doing, it stratifies patients according to their immune profiles, which may lead to effective computational prognostic and clinical guides. [ABSTRACT FROM AUTHOR]

Abstract: Vaccines represent a special class of drugs, capable of stimulating immune system responses against pathogens and tumors. Vaccine development is a lengthy process that includes expensive laboratory experiments in order to assess safety and effectiveness. As the efficacy of a vaccine was demonstrated by biological/chemical investigations and pre-clinical studies, then a major problem is represented by the search for an optimal vaccination dosage. Optimality here assumes the meaning of assuring a high degree of efficacy and safety (lack of toxic or side effects). In lack of quantitative methods, this is usually achieved by a consensus technique, a public statement on a particular aspect of medical knowledge available at the time it was written, and that is generally agreed upon as the evidence-based, state-of-the-art (or state-of-science) knowledge by a representative group of experts in that area. In this article, we focus on the difficult problem of the search for an optimal vaccination dosage in the field of tumor immunology, that is a major issue in biomedical research. This, indeed, represents a first step toward a personalized medicine approach. [Copyright &y& Elsevier]

Clostridium difficile is a spore-forming gram-positive bacterium, recognized as the primary cause of antibiotic-associated nosocomial diarrhoea. Clostridium difficile infection (CDI) has emerged as a major health-associated infection with increased incidence and hospitalization over the years with high mortality rates. Contamination and infection occur after ingestion of vegetative spores, which germinate in the gastro-intestinal tract. The surface layer protein and flagellar proteins are responsible for the bacterial colonization while the spore coat protein, is associated with spore colonization. Both these factors are the main concern of the recurrence of CDI in hospitalized patients. In this study, the CotE, SlpA and FliC proteins are chosen to form a multivalent, multi-epitopic, chimeric vaccine candidate using the immunoinformatics approach. The overall reliability of the candidate vaccine was validated in silico and the molecular dynamics simulation verified the stability of the vaccine designed. Docking studies showed stable vaccine interactions with Toll‐Like Receptors of innate immune cells and MHC receptors. In silico codon optimization of the vaccine and its insertion in the cloning vector indicates a competent expression of the modelled vaccine in E. coli expression system. An in silico immune simulation system evaluated the effectiveness of the candidate vaccine to trigger a protective immune response. [ABSTRACT FROM AUTHOR]

Systems medicine (SM) has emerged as a powerful tool for studying the human body at the systems level with the aim of improving our understanding, prevention and treatment of complex diseases. Being able to automatically extract relevant features needed for a given task from high-dimensional, heterogeneous data, deep learning (DL) holds great promise in this endeavour. This review paper addresses the main developments of DL algorithms and a set of general topics where DL is decisive, namely, within the SM landscape. It discusses how DL can be applied to SM with an emphasis on the applications to predictive, preventive and precision medicine. Several key challenges have been highlighted including delivering clinical impact and improving interpretability. We used some prototypical examples to highlight the relevance and significance of the adoption of DL in SM, one of them is involving the creation of a model for personalized Parkinson's disease. The review offers valuable insights and informs the research in DL and SM. [ABSTRACT FROM AUTHOR]

Background: The aim of a recent research project was the investigation of the mechanisms involved in the onset of type 2 diabetes in the absence of familiarity. This has led to the development of a computational model that recapitulates the aetiology of the disease and simulates the immunological and metabolic alterations linked to type-2 diabetes subjected to clinical, physiological, and behavioural features of prototypical human individuals. Results: We analysed the time course of 46,170 virtual subjects, experiencing different lifestyle conditions. We then set up a statistical model able to recapitulate the simulated outcomes. Conclusions: The resulting machine learning model adequately predicts the synthetic dataset and can, therefore, be used as a computationally-cheaper version of the detailed mathematical model, ready to be implemented on mobile devices to allow self-assessment by informed and aware individuals. The computational model used to generate the dataset of this work is available as a web-service at the following address: http://kraken.iac.rm.cnr.it/T2DM. [ABSTRACT FROM AUTHOR]

In the past two decades, 7 coronaviruses have infected the human population, with two major outbreaks caused by SARS-CoV and MERS-CoV in the year 2002 and 2012, respectively. Currently, the entire world is facing a pandemic of another coronavirus, SARS-CoV-2, with a high fatality rate. The spike glycoprotein of SARS-CoV-2 mediates entry of virus into the host cell and is one of the most important antigenic determinants, making it a potential candidate for a vaccine. In this study, we have computationally designed a multi-epitope vaccine using spike glycoprotein of SARS-CoV-2. The overall quality of the candidate vaccine was validated in silico and Molecular Dynamics Simulation confirmed the stability of the designed vaccine. Docking studies revealed stable interactions of the vaccine with Toll-Like Receptors and MHC Receptors. The in silico cloning and codon optimization supported the proficient expression of the designed vaccine in E. coli expression system. The efficiency of the candidate vaccine to trigger an effective immune response was assessed by an in silico immune simulation. The computational analyses suggest that the designed multi-epitope vaccine is structurally stable which can induce specific immune responses and thus, can be a potential vaccine candidate against SARS-CoV-2. [ABSTRACT FROM AUTHOR]

In the past two decades, 7 coronaviruses have infected the human population, with two major outbreaks caused by SARS-CoV and MERS-CoV in the year 2002 and 2012, respectively. Currently, the entire world is facing a pandemic of another coronavirus, SARS-CoV-2, with a high fatality rate. The spike glycoprotein of SARS-CoV-2 mediates entry of virus into the host cell and is one of the most important antigenic determinants, making it a potential candidate for a vaccine. In this study, we have computationally designed a multi-epitope vaccine using spike glycoprotein of SARS-CoV-2. The overall quality of the candidate vaccine was validated in silico and Molecular Dynamics Simulation confirmed the stability of the designed vaccine. Docking studies revealed stable interactions of the vaccine with Toll-Like Receptors and MHC Receptors. The in silico cloning and codon optimization supported the proficient expression of the designed vaccine in E. coli expression system. The efficiency of the candidate vaccine to trigger an effective immune response was assessed by an in silico immune simulation. The computational analyses suggest that the designed multi-epitope vaccine is structurally stable which can induce specific immune responses and thus, can be a potential vaccine candidate against SARS-CoV-2. [ABSTRACT FROM AUTHOR]

Systems medicine holds many promises, but has so far provided only a limited number of proofs of principle. To address this road block, possible barriers and challenges of translating systems medicine into clinical practice need to be identified and addressed. The members of the European Cooperation in Science and Technology (COST) Action CA15120 Open Multiscale Systems Medicine (OpenMultiMed) wish to engage the scientific community of systems medicine and multiscale modelling, data science and computing, to provide their feedback in a structured manner. This will result in follow-up white papers and open access resources to accelerate the clinical translation of systems medicine. [ABSTRACT FROM AUTHOR]

The beneficial effects of physical activity for the prevention and management of several chronic diseases are widely recognized. Mathematical modeling of the effects of physical exercise in body metabolism and in particular its influence on the control of glucose homeostasis is of primary importance in the development of eHealth monitoring devices for a personalized medicine. Nonetheless, to date only a few mathematical models have been aiming at this specific purpose. We have developed a whole-body computational model of the effects on metabolic homeostasis of a bout of physical exercise. Built upon an existing model, it allows to detail better both subjects’ characteristics and physical exercise, thus determining to a greater extent the dynamics of the hormones and the metabolites considered. [ABSTRACT FROM AUTHOR]

Motivation: A computational model equipped with the main immunological features of the sea bass (Dicentrarchus labrax L.) immune system was used to predict more effective vaccination in fish. The performance of the model was evaluated by using the results of two in vivo vaccinations trials against L. anguillarum and P. damselae. Results: Tests were performed to select the appropriate doses of vaccine and infectious bacteria to set up the model. Simulation outputs were compared with the specific antibody production and the expression of BcR and TcR gene transcripts in spleen. The model has shown a good ability to be used in sea bass and could be implemented for different routes of vaccine administration even with more than two pathogens. The model confirms the suitability of in silico methods to optimize vaccine doses and the immune response to them. This model could be applied to other species to optimize the design of new vaccination treatments of fish in aquaculture. [ABSTRACT FROM AUTHOR]

Interleukin-6 (IL-6) has been recently shown to play a central role in glucose homeostasis, since it stimulates the production and secretion of Glucagon-like Peptide-1 (GLP-1) from intestinal L-cells and pancreas, leading to an enhanced insulin response. In resting conditions, IL-6 is mainly produced by the adipose tissue whereas, during exercise, skeletal muscle contractions stimulate a marked IL-6 secretion as well. Available mathematical models describing the effects of exercise on glucose homeostasis, however, do not account for this IL-6 contribution. This study aimed at developing and validating a system model of exercise’s effects on plasma IL-6 dynamics in healthy humans, combining the contributions of both adipose tissue and skeletal muscle. A two-compartment description was adopted to model plasma IL-6 changes in response to oxygen uptake’s variation during an exercise bout. The free parameters of the model were estimated by means of a cross-validation procedure performed on four different datasets. A low coefficient of variation (<10%) was found for each parameter and the physiologically meaningful parameters were all consistent with literature data. Moreover, plasma IL-6 dynamics during exercise and post-exercise were consistent with literature data from exercise protocols differing in intensity, duration and modality. The model successfully emulated the physiological effects of exercise on plasma IL-6 levels and provided a reliable description of the role of skeletal muscle and adipose tissue on the dynamics of plasma IL-6. The system model here proposed is suitable to simulate IL-6 response to different exercise modalities. Its future integration with existing models of GLP-1-induced insulin secretion might provide a more reliable description of exercise’s effects on glucose homeostasis and hence support the definition of more tailored interventions for the treatment of type 2 diabetes. [ABSTRACT FROM AUTHOR]

Immunosenescence is thought to result from cellular aging and to reflect exposure to environmental stressors and antigens, including cytomegalovirus (CMV). However, not all of the features of immunosenescence are consistent with this view, and this has led to the emergence of the sister theory of “inflammaging”. The recently discovered diffuse tissue distribution of resident memory T cells (T RM ) which don't recirculate, calls these theories into question. These cells account for most T cells residing in barrier epithelia which sit in and travel through the extracellular matrix (ECM). With almost all studies to date carried out on peripheral blood, the age-related changes of the ECM and their consequences for T cell mobility, which is crucial for the function of these cells, have been largely ignored. We propose an update of the theoretical framework of immunosenescence, based on a novel hypothesis: the increasing stiffness and cross-linking of the senescent ECM lead to a progressive immunodeficiency due to an age-related decrease in T cell mobility and eventually the death of these cells. A key element of this mechanism is the mechanical stress to which the cell cytoplasm and nucleus are subjected during passage through the ECM. This hypothesis is based on an “evo-devo” perspective bringing together some major characteristics of aging, to create a single interpretive framework for immunosenescence. [ABSTRACT FROM AUTHOR]

Objectives: Identifying the predictive factors of Sustained Virological Response (SVR) represents an important challenge in new interferon-based DAA therapies. Here, we analyzed the kinetics of antiviral response associated with a triple drug regimen, and the association between negative residual viral load at different time points during treatment. Methods: Twenty-three HCV genotype 1 (GT 1a n = 11; GT1b n = 12) infected patients were included in the study. Linear Discriminant Analysis (LDA) was used to establish possible association between HCV RNA values at days 1 and 4 from start of therapy and SVR. Principal component analysis (PCA) was applied to analyze the correlation between HCV RNA slope and SVR. A ultrasensitive (US) method was established to measure the residual HCV viral load in those samples which resulted “detected <12IU/ml” or undetectable with ABBOTT standard assay, and was retrospectively used on samples collected at different time points to establish its predictive power for SVR. Results: According to LDA, there was no association between SVR and viral kinetics neither at time points earlier than 1 week (days 1 and 4) after therapy initiation nor later. The slopes were not relevant for classifying patients as SVR or no-SVR. No significant differences were observed in the median HCV RNA values at T0 among SVR and no-SVR patients. HCV RNA values with US protocol (LOD 1.2 IU/ml) after 1 month of therapy were considered; the area under the ROC curve was 0.70. Overall, PPV and NPV of undetectable HCV RNA with the US method for SVR was 100% and 46.7%, respectively; sensitivity and specificity were 38.4% and 100% respectively. Conclusion: HCV RNA “not detected” by the US method after 1 month of treatment is predictive of SVR in first generation Protease inhibitor (PI)-based triple therapy. The US method could have clinical utility for advanced monitoring of virological response in new interferon based DAA combination regimens. [ABSTRACT FROM AUTHOR]

Motivation: The recognition of antigenic peptides is a major event of an immune response. In current mesoscopic-scale simulators of the immune system, this crucial step has been modeled in a very approximated way.Results: We have equipped an agent-based model of the immune system with immuno-informatics methods to allow the simulation of the cardinal events of the antigenic recognition, going from single peptides to whole proteomes.The recognition process accounts for B cell-epitopes prediction through Parker-scale affinity estimation, class I and II HLA peptide prediction and binding through position-specific scoring matrices based on information from known HLA epitopes prediction tools, and TCR binding to HLA–peptide complex calculated as the averaged sum of a residue–residue contact potential.These steps are executed for all lymphocytes agents encountering the antigen in a wide-reaching Monte Carlo simulation.Availability: http://www.cbs.dtu.dk/services/C-ImmSim-10.1/Contact: f.castiglione@iac.cnr.it [ABSTRACT FROM AUTHOR]

Background Direct-acting antiviral drugs (DAA) regimen improve the SVR rate. However, adverse effects often lead to therapy interruption. This underlines the importance to find some predictive parameters of response in order to consider the possibility of a shorter time of antiviral treatment in the appearance of adverse effects without affecting the success of the therapy. Objectives We aimed to examine the HCVAg kinetics in the early phase of treatment and its predictive value of SVR in patients undergoing TPV/Peg-IFN/RBV treatment. Study design Twenty-three patients infected by HCV genotype 1 (1a n = 11; 1b n = 12) were included in this prospective study. Results At baseline the median Log of HCVAg concentration in RVR and EVR patients were 3.15 fmol/L and 3.45 fmol/L, respectively with no significant differences. The baseline median HCV-RNA to HCVAg ratio was 233.77, this ratio was significantly lower when measured on day 1 (27.52) and on day 6 (24.84) ( p < 0.001). The two-tailed Fisher’s exact test indicated that the SVR response is statistically significantly different in patients with detected HCVAg at week1 compared to patients with no detectable HCVAg ( p = 0.05). The sensitivity, specificity, and negative and positive predictive values (NPV, PPV) were 53.8, 87.5, 53.8 and 87.5%, respectively. The area under the ROC curve was 0.71 at day T6, the best cut-off of 3 fmol/L when evaluated with the HCVAg plasma concentration at dayT6. Conclusion Undetectable HCVAg in the early phase of TPV/Peg-IFN/RBV treatment could represent an important parameter for predicting SVR. [ABSTRACT FROM AUTHOR]

We developed an approach that integrates different network-based methods to analyze the correlation network arising from large-scale gene expression data. By studying grapevine (Vitis vinifera) and tomato (Solanum lycopersicum) gene expression atlases and a grapevine berry transcriptomic data set during the transition from immature to mature growth, we identified a category named "fight-club hubs" characterized by a marked negative correlation with the expression profiles of neighboring genes in the network. A special subset named "switch genes" was identified, with the additional property of many significant negative correlations outside their own group in the network. Switch genes are involved in multiple processes and include transcription factors that may be considered master regulators of the previously reported transcriptome remodeling that marks the developmental shift from immature to mature growth. All switch genes, expressed at low levels in vegetative/green tissues, showed a significant increase in mature/woody organs, suggesting a potential regulatory role during the developmental transition. Finally, our analysis of tomato gene expression data sets showed that wild-type switch genes are downregulated in ripening-deficient mutants. The identification of known master regulators of tomato fruit maturation suggests our method is suitable for the detection of key regulators of organ development in different fleshy fruit crops. [ABSTRACT FROM AUTHOR]

Peptide vaccination in cancer therapy is a promising alternative to conventional methods. However, the parameters for this personalized treatment are difficult to access experimentally. In this respect, in silico models can help to narrow down the parameter space or to explain certain phenomena at a systems level. Herein, we develop two empirical interaction potentials specific to B-cell and T-cell receptor complexes and validate their applicability in comparison to a more general potential. The interaction potentials are applied to the model VaccImm which simulates the immune response against solid tumors under peptide vaccination therapy. This multi-agent system is derived from another immune system simulator (CImmSim) and now includes a module that enables the amino acid sequence of immune receptors and their ligands to be taken into account. The multi-agent approach is combined with approved methods for prediction of major histocompatibility complex (MHC)-binding peptides and the newly developed interaction potentials. In the analysis, we critically assess the impact of the different modules on the simulation with VaccImm and how they influence each other. In addition, we explore the reasons for failures in inducing an immune response by examining the activation states of the immune cell populations in detail. In summary, the present work introduces immune-specific interaction potentials and their application to the agent-based model VaccImm which simulates peptide vaccination in cancer therapy. [ABSTRACT FROM AUTHOR]

Background: When to initiate antiretroviral therapy in HIV infected patients is a diffcult clinical decision. Actually, it is still a matter of discussion whether early highly active antiretroviral therapy (HAART) during primary HIV infection may influence the dynamics of the viral rebound, in case of therapy interruption, and overall the main disease course. Methods: In this article we use a computational model and clinical data to identify the role of HAART timing on the residual capability to control HIV rebound after treatment suspension. Analyses of clinical data from three groups of patients initiating HAART respectively before seroconversion (very early), during the acute phase (early) and in the chronic phase (late), evidence differences arising from the very early events of the viral infection. Results: The computational model allows a fine grain assessment of the impact of HAART timing on the disease outcome, from acute to chronic HIV-1 infection. Both patients' data and computer simulations reveal that HAART timing may indeed affect the HIV control capability after treatment discontinuation. In particular, we find a median time to viral rebound that is significantly longer in very early than in late patients. Conclusions: A timing threshold is identified, corresponding to approximately three weeks post-infection, after which the capability to control HIV replication is lost. Conversely, HAART initiation occurring within three weeks from the infection could allow to preserve a significant control capability. This time could be related to the global triggering of uncontrolled immune activation, affecting residual immune competence preservation and HIV reservoir establishment. [ABSTRACT FROM AUTHOR]

Two T helper (Th) cell subsets, namely Th1 and Th2 cells, play an important role in inflammatory diseases. The two subsets are thought to counter-regulate each other, and alterations in their balance result in different diseases. This paradigm has been challenged by recent clinical and experimental data. Because of the large number of genes involved in regulating Th1 and Th2 cells, assessment of this paradigm by modeling or experiments is difficult. Novel algorithms based on formal methods now permit the analysis of large gene regulatory networks. By combining these algorithms with in silico knockouts and gene expression microarray data from human T cells, we examined if the results were compatible with a counterregulatory role of Th1 and Th2 cells. We constructed a directed network model of genes regulating Th1 and Th2 cells through text mining and manual curation. We identified four attractors in the network, three of which included genes that corresponded to Th0, Th1 and Th2 cells. The fourth attractor contained a mixture of Th1 and Th2 genes. We found that neither in silico knockouts of the Th1 and Th2 attractor genes nor gene expression microarray data from patients with immunological disorders and healthy subjects supported a counter-regulatory role of Th1 and Th2 cells. By combining network modeling with transcriptomic data analysis and in silico knockouts, we have devised a practical way to help unravel complex regulatory network topology and to increase our understanding of how network actions may differ in health and disease. [ABSTRACT FROM AUTHOR]

We present a new approach to the study of the immune system that combines techniques of systems biology with information provided by data-driven prediction methods. To this end, we have extended an agent-based simulator of the immune response, C-IMMSIM, such that it represents pathogens, as well as lymphocytes receptors, by means of their amino acid sequences and makes use of bioinformatics methods for T and B cell epitope prediction. This is a key step for the simulation of the immune response, because it determines immunogenicity. The binding of the epitope, which is the immunogenic part of an invading pathogen, together with activation and cooperation from T helper cells, is required to trigger an immune response in the affected host. To determine a pathogen’s epitopes, we use existing prediction methods. In addition, we propose a novel method, which uses Miyazawa and Jernigan protein–protein potential measurements, for assessing molecular binding in the context of immune complexes. We benchmark the resulting model by simulating a classical immunization experiment that reproduces the development of immune memory. We also investigate the role of major histocompatibility complex (MHC) haplotype heterozygosity and homozygosity with respect to the influenza virus and show that there is an advantage to heterozygosity. Finally, we investigate the emergence of one or more dominating clones of lymphocytes in the situation of chronic exposure to the same immunogenic molecule and show that high affinity clones proliferate more than any other. These results show that the simulator produces dynamics that are stable and consistent with basic immunological knowledge. We believe that the combination of genomic information and simulation of the dynamics of the immune system, in one single tool, can offer new perspectives for a better understanding of the immune system. [ABSTRACT FROM AUTHOR]

Background: When to initiate antiretroviral therapy in HIV infected patients is a difficult clinical decision. Actually, it is still a matter of discussion whether early highly active antiretroviral therapy (HAART) during primary HIV infection may influence the dynamics of the viral rebound, in case of therapy interruption, and overall the main disease course.Methods: In this article we use a computational model and clinical data to identify the role of HAART timing on the residual capability to control HIV rebound after treatment suspension. Analyses of clinical data from three groups of patients initiating HAART respectively before seroconversion (very early), during the acute phase (early) and in the chronic phase (late), evidence differences arising from the very early events of the viral infection.Results: The computational model allows a fine grain assessment of the impact of HAART timing on the disease outcome, from acute to chronic HIV-1 infection. Both patients' data and computer simulations reveal that HAART timing may indeed affect the HIV control capability after treatment discontinuation. In particular, we find a median time to viral rebound that is significantly longer in very early than in late patients.Conclusions: A timing threshold is identified, corresponding to approximately three weeks post-infection, after which the capability to control HIV replication is lost. Conversely, HAART initiation occurring within three weeks from the infection could allow to preserve a significant control capability. This time could be related to the global triggering of uncontrolled immune activation, affecting residual immune competence preservation and HIV reservoir establishment. [ABSTRACT FROM AUTHOR]