Agent-based models (ABMs) have enabled great advances in the study of tumor development and therapeutic response, allowing researchers to explore the spatiotemporal evolution of the tumor and its ...microenvironment. However, these models face serious drawbacks in the realm of parameterization - ABM parameters are typically set individually based on various data and literature sources, rather than through a rigorous parameter estimation approach. While ABMs can be fit to simple time-course data (such as tumor volume), that type of data loses the spatial information that is a defining feature of ABMs. While tumor images provide spatial information, it is exceedingly difficult to compare tumor images to ABM simulations beyond a qualitative visual comparison. Without a quantitative method of comparing the similarity of tumor images to ABM simulations, a rigorous parameter fitting is not possible. Here, we present a novel approach that applies neural networks to represent both tumor images and ABM simulations as low dimensional points, with the distance between points acting as a quantitative measure of difference between the two. This enables a quantitative comparison of tumor images and ABM simulations, where the distance between simulated and experimental images can be minimized using standard parameter-fitting algorithms. Here, we describe this method and present two examples to demonstrate the application of the approach to estimate parameters for two distinct ABMs. Overall, we provide a novel method to robustly estimate ABM parameters.
Mathematical modeling provides the predictive ability to understand the metabolic reprogramming and complex pathways that mediate cancer cells' proliferation. We present a mathematical model using a ...multiscale, multicellular approach to simulate avascular tumor growth, applied to pancreatic cancer. The model spans three distinct spatial and temporal scales. At the extracellular level, reaction diffusion equations describe nutrient concentrations over a span of seconds. At the cellular level, a lattice-based energy driven stochastic approach describes cellular phenomena including adhesion, proliferation, viability and cell state transitions, occurring on the timescale of hours. At the sub-cellular level, we incorporate a detailed kinetic model of intracellular metabolite dynamics on the timescale of minutes, which enables the cells to uptake and excrete metabolites and use the metabolites to generate energy and building blocks for cell growth. This is a particularly novel aspect of the model. Certain defined criteria for the concentrations of intracellular metabolites lead to cancer cell growth, proliferation or death. Overall, we model the evolution of the tumor in both time and space. Starting with a cluster of tumor cells, the model produces an avascular tumor that quantitatively and qualitatively mimics experimental measurements of multicellular tumor spheroids. Through our model simulations, we can investigate the response of individual intracellular species under a metabolic perturbation and investigate how that response contributes to the response of the tumor as a whole. The predicted response of intracellular metabolites under various targeted strategies are difficult to resolve with experimental techniques. Thus, the model can give novel predictions as to the response of the tumor as a whole, identifies potential therapies to impede tumor growth, and predicts the effects of those therapeutic strategies. In particular, the model provides quantitative insight into the dynamic reprogramming of tumor cells at the intracellular level in response to specific metabolic perturbations. Overall, the model is a useful framework to study targeted metabolic strategies for inhibiting tumor growth.
Within the tumor microenvironment, macrophages exist in an immunosuppressive state, preventing T cells from eliminating the tumor. Due to this, research is focusing on immunotherapies that ...specifically target macrophages in order to reduce their immunosuppressive capabilities and promote T cell function. In this study, we develop an agent-based model consisting of the interactions between macrophages, T cells, and tumor cells to determine how the immune response changes due to three macrophage-based immunotherapeutic strategies: macrophage depletion, recruitment inhibition, and macrophage reeducation. We find that reeducation, which converts the macrophages into an immune-promoting phenotype, is the most effective strategy and that the macrophage recruitment rate and tumor proliferation rate (tumor-specific properties) have large impacts on therapy efficacy. We also employ a novel method of using a neural network to reduce the computational complexity of an intracellular signaling mechanistic model.
Angiogenesis plays an important role in the survival of tissues, as blood vessels provide oxygen and nutrients required by the resident cells. Thus, targeting angiogenesis is a prominent strategy in ...many different settings, including both tissue engineering and cancer treatment. However, not all of the approaches that modulate angiogenesis lead to successful outcomes. Angiogenesis-based therapies primarily target pro-angiogenic factors such as vascular endothelial growth factor-A (VEGF) or fibroblast growth factor (FGF) in isolation, and there is a limited understanding of how these promoters combine together to stimulate angiogenesis. Targeting one pathway could be insufficient, as alternative pathways may compensate, diminishing the overall effect of the treatment strategy.
To gain mechanistic insight and identify novel therapeutic strategies, we have developed a detailed mathematical model to quantitatively characterize the crosstalk of FGF and VEGF intracellular signaling. The model focuses on FGF- and VEGF-induced mitogen-activated protein kinase (MAPK) signaling to promote cell proliferation and the phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt) pathway, which promotes cell survival and migration. We fit the model to published experimental datasets that measure phosphorylated extracellular regulated kinase (pERK) and Akt (pAkt) upon FGF or VEGF stimulation. We validate the model with separate sets of data.
We apply the trained and validated mathematical model to characterize the dynamics of pERK and pAkt in response to the mono- and co-stimulation by FGF and VEGF. The model predicts that for certain ranges of ligand concentrations, the maximum pERK level is more responsive to changes in ligand concentration compared to the maximum pAkt level. Also, the combination of FGF and VEGF indicates a greater effect in increasing the maximum pERK compared to the summation of individual effects, which is not seen for maximum pAkt levels. In addition, our model identifies the influential species and kinetic parameters that specifically modulate the pERK and pAkt responses, which represent potential targets for angiogenesis-based therapies.
Overall, the model predicts the combination effects of FGF and VEGF stimulation on ERK and Akt quantitatively and provides a framework to mechanistically explain experimental results and guide experimental design. Thus, this model can be utilized to study the effects of pro- and anti-angiogenic therapies that particularly target ERK and/or Akt activation upon stimulation with FGF and VEGF. Video Abstract.
Pancreatic β-cells respond to increased extracellular glucose levels by initiating a metabolic shift. That change in metabolism is part of the process of glucose-stimulated insulin secretion and is ...of particular interest in the context of diabetes. However, we do not fully understand how the coordinated changes in metabolic pathways and metabolite products influence insulin secretion. In this work, we apply systems biology approaches to develop a detailed kinetic model of the intracellular central carbon metabolic pathways in pancreatic β-cells upon stimulation with high levels of glucose. The model is calibrated to published metabolomics datasets for the INS1 823/13 cell line, accurately capturing the measured metabolite fold-changes. We first employed the calibrated mechanistic model to estimate the stimulated cell’s fluxome. We then used the predicted network fluxes in a data-driven approach to build a partial least squares regression model. By developing the combined kinetic and data-driven modeling framework, we gain insights into the link between β-cell metabolism and glucose-stimulated insulin secretion. The combined modeling framework was used to predict the effects of common anti-diabetic pharmacological interventions on metabolite levels, flux through the metabolic network, and insulin secretion. Our simulations reveal targets that can be modulated to enhance insulin secretion. The model is a promising tool to contextualize and extend the usefulness of metabolomics data and to predict dynamics and metabolite levels that are difficult to measure
in vitro
. In addition, the modeling framework can be applied to identify, explain, and assess novel and clinically-relevant interventions that may be particularly valuable in diabetes treatment.
Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to ...optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates.
In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux.
Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.
The MiniBooNE experiment at Fermilab reports a total excess of 638.0 ± 52.1 (stat .) ± 122.2 (syst.) electronlike events from a data sample corresponding to 18.75 × 10 20 protons-on-target in ...neutrino mode, which is a 46% increase in the data sample with respect to previously published results and 11.27 × 10 20 protons-on-target in antineutrino mode. The overall significance of the excess, 4.8σ, is limited by systematic uncertainties, assumed to be Gaussian, as the statistical significance of the excess is 12.2σ. The additional statistics allow several studies to address questions on the source of the excess. First, we provide two-dimensional plots in visible energy and the cosine of the angle of the outgoing lepton, which can provide valuable input to models for the event excess. Second, we test whether the excess may arise from photons that enter the detector from external events or photons exiting the detector from π 0 decays in two model independent ways. Beam timing information shows that almost all of the excess is in time with neutrinos that interact in the detector. The radius distribution shows that the excess is distributed throughout the volume, while tighter cuts on the fiducial volume increase the significance of the excess. The data likelihood ratio disfavors models that explain the event excess due to entering or exiting photons.
Autophagy is involved with the turnover of intracellular components and the management of stress responses. Genetic studies in mice have shown that suppression of neuronal autophagy can lead to the ...accumulation of protein aggregates and neurodegeneration. However, no study has shown that increasing autophagic gene expression can be beneficial to an aging nervous system. Here we demonstrate that expression of several autophagy genes is reduced in Drosophila neural tissues as a normal part of aging. The age-dependent suppression of autophagy occurs concomitantly with the accumulation of insoluble ubiquitinated proteins (IUP), a marker of neuronal aging and degeneration. Mutations in the Atg8a gene (autophagy-related 8a) result in reduced lifespan, IUP accumulation and increased sensitivity to oxidative stress. In contrast, enhanced Atg8a expression in older fly brains extends the average adult lifespan by 56% and promotes resistance to oxidative stress and the accumulation of ubiquitinated and oxidized proteins. These data indicate that genetic or age-dependent suppression of autophagy is closely associated with the buildup of cellular damage in neurons and a reduced lifespan, while maintaining the expression of a rate-limiting autophagy gene prevents the age-dependent accumulation of damage in neurons and promotes longevity.
Abstract
Background
Chimeric antigen receptor (CAR)-expressing cells are a powerful modality of adoptive cell therapy against cancer. The potency of signaling events initiated upon antigen binding ...depends on the costimulatory domain within the structure of the CAR. One such costimulatory domain is 4-1BB, which affects cellular response via the NFκB pathway. However, the quantitative aspects of 4-1BB-induced NFκB signaling are not fully understood.
Methods
We developed an ordinary differential equation-based mathematical model representing canonical NFκB signaling activated by CD19scFv-4-1BB. After a global sensitivity analysis on model parameters, we ran Monte Carlo simulations of cell population-wide variability in NFκB signaling and quantified the mutual information between the extracellular signal and different levels of the NFκB signal transduction pathway.
Results
In response to a wide range of antigen concentrations, the magnitude of the transient peak in NFκB nuclear concentration varies significantly, while the timing of this peak is relatively consistent. Global sensitivity analysis showed that the model is robust to variations in parameters, and thus, its quantitative predictions would remain applicable to a broad range of parameter values. The model predicts that overexpressing NEMO and disabling IKKβ deactivation can increase the mutual information between antigen levels and NFκB activation.
Conclusions
Our modeling predictions provide actionable insights to guide CAR development. Particularly, we propose specific manipulations to the NFκB signal transduction pathway that can fine-tune the response of CD19scFv-4-1BB cells to the antigen concentrations they are likely to encounter.
ABSTRACT We describe the operation and performance of the difference imaging pipeline (DiffImg) used to detect transients in deep images from the Dark Energy Survey Supernova program (DES-SN) in its ...first observing season from 2013 August through 2014 February. DES-SN is a search for transients in which ten 3 deg2 fields are repeatedly observed in the g, r, i, z passbands with a cadence of about 1 week. The observing strategy has been optimized to measure high-quality light curves and redshifts for thousands of Type Ia supernovae (SNe Ia) with the goal of measuring dark energy parameters. The essential DiffImg functions are to align each search image to a deep reference image, do a pixel-by-pixel subtraction, and then examine the subtracted image for significant positive detections of point-source objects. The vast majority of detections are subtraction artifacts, but after selection requirements and image filtering with an automated scanning program, there are ∼130 detections per deg2 per observation in each band, of which only ∼25% are artifacts. Of the ∼7500 transients discovered by DES-SN in its first observing season, each requiring a detection on at least two separate nights, Monte Carlo (MC) simulations predict that 27% are expected to be SNe Ia or core-collapse SNe. Another ∼30% of the transients are artifacts in which a small number of observations satisfy the selection criteria for a single-epoch detection. Spectroscopic analysis shows that most of the remaining transients are AGNs and variable stars. Fake SNe Ia are overlaid onto the images to rigorously evaluate detection efficiencies and to understand the DiffImg performance. The DiffImg efficiency measured with fake SNe agrees well with expectations from a MC simulation that uses analytical calculations of the fluxes and their uncertainties. In our 8 "shallow" fields with single-epoch 50% completeness depth ∼23.5, the SN Ia efficiency falls to 1/2 at redshift z 0.7; in our 2 "deep" fields with mag-depth ∼24.5, the efficiency falls to 1/2 at z 1.1. A remaining performance issue is that the measured fluxes have additional scatter (beyond Poisson fluctuations) that increases with the host galaxy surface brightness at the transient location. This bright-galaxy issue has minimal impact on the SNe Ia program, but it may lower the efficiency for finding fainter transients on bright galaxies.