Explore the vast range of titles available.

Many multicellular systems problems can only be understood by studying how cells move, grow, divide, interact, and die. Tissue-scale dynamics emerge from systems of many interacting cells as they respond to and influence their microenvironment. The ideal \"virtual laboratory\" for such multicellular systems simulates both the biochemical microenvironment (the \"stage\") and many mechanically and biochemically interacting cells (the \"players\" upon the stage). PhysiCell-physics-based multicellular simulator-is an open source agent-based simulator that provides both the stage and the players for studying many interacting cells in dynamic tissue microenvironments. It builds upon a multi-substrate biotransport solver to link cell phenotype to multiple diffusing substrates and signaling factors. It includes biologically-driven sub-models for cell cycling, apoptosis, necrosis, solid and fluid volume changes, mechanics, and motility \"out of the box.\" The C++ code has minimal dependencies, making it simple to maintain and deploy across platforms. PhysiCell has been parallelized with OpenMP, and its performance scales linearly with the number of cells. Simulations up to 105-106 cells are feasible on quad-core desktop workstations; larger simulations are attainable on single HPC compute nodes. We demonstrate PhysiCell by simulating the impact of necrotic core biomechanics, 3-D geometry, and stochasticity on the dynamics of hanging drop tumor spheroids and ductal carcinoma in situ (DCIS) of the breast. We demonstrate stochastic motility, chemical and contact-based interaction of multiple cell types, and the extensibility of PhysiCell with examples in synthetic multicellular systems (a \"cellular cargo delivery\" system, with application to anti-cancer treatments), cancer heterogeneity, and cancer immunology. PhysiCell is a powerful multicellular systems simulator that will be continually improved with new capabilities and performance improvements. It also represents a significant independent code base for replicating results from other simulation platforms. The PhysiCell source code, examples, documentation, and support are available under the BSD license at http://PhysiCell.MathCancer.org and http://PhysiCell.sf.net.

Water-to-air carbon dioxide fluxes from tropical lakes and reservoirs (artificial lakes) may be an important but understudied component of global carbon fluxes. Here, we investigate the seasonal dissolved carbon dioxide (CO2) dynamics in a lake and a reservoir on a tropical volcanic island (Bali, Indonesia). Observations were performed over four seasonal surveys in Bali's largest natural lake (Lake Batur) and largest reservoir (Palasari Reservoir). Average CO2 partial pressures in the natural lake and reservoir were 263.7±12.2 μatm and 785.0±283.6 μatm respectively, with the highest area-weighted partial pressures in the wet season for both systems. The strong correlations between seasonal mean values of dissolved oxygen (DO) and pCO2 in the natural lake (r2 = 0.92) suggest that surface water metabolism was an important driver of CO2 dynamics in this deep system. Radon (222Rn, a natural groundwater discharge tracer) explained up to 77% of the variability in pCO2 in the shallow reservoir, suggesting that groundwater seepage was the major CO2 driver in the reservoir. Overall, the natural lake was a sink of atmospheric CO2 (average fluxes of -2.8 mmol m-2 d-1) while the reservoir was a source of CO2 to the atmosphere (average fluxes of 7.3 mmol m-2 d-1). Reservoirs are replacing river valleys and terrestrial ecosystems, particularly throughout developing tropical regions. While the net effect of this conversion on atmospheric CO2 fluxes remains to be resolved, we speculate that reservoir construction will partially offset the CO2 sink provided by deep, volcanic, natural lakes and terrestrial environments.

Mechanistic learning refers to the synergistic combination of mechanistic mathematical modeling and data-driven machine or deep learning. This emerging field finds increasing applications in (mathematical) oncology. This review aims to capture the current state of the field and provides a perspective on how mechanistic learning may progress in the oncology domain. We highlight the synergistic potential of mechanistic learning and point out similarities and differences between purely data-driven and mechanistic approaches concerning model complexity, data requirements, outputs generated, and interpretability of the algorithms and their results. Four categories of mechanistic learning (sequential, parallel, extrinsic, intrinsic) of mechanistic learning are presented with specific examples. We discuss a range of techniques including physics-informed neural networks, surrogate model learning, and digital twins. Example applications address complex problems predominantly from the domain of oncology research such as longitudinal tumor response predictions or time-to-event modeling. As the field of mechanistic learning advances, we aim for this review and proposed categorization framework to foster additional collaboration between the data- and knowledge-driven modeling fields. Further collaboration will help address difficult issues in oncology such as limited data availability, requirements of model transparency, and complex input data which are embraced in a mechanistic learning framework

The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-β). The TGF-β-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-β sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-β, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-β, and the duration of TGF-β sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-β sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-β, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-β sources changed the fibrotic patterns; in the presence of persistent TGF-β sources, partial removal of TGF-β from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-β sources may find applications in the current clinical trials targeting TGF-β for the resolution of lung fibrosis.

Background Cancer is a complex, multiscale dynamical system, with interactions between tumor cells and non-cancerous host systems. Therapies act on this combined cancer-host system, sometimes with unexpected results. Systematic investigation of mechanistic computational models can augment traditional laboratory and clinical studies, helping identify the factors driving a treatment’s success or failure. However, given the uncertainties regarding the underlying biology, these multiscale computational models can take many potential forms, in addition to encompassing high-dimensional parameter spaces. Therefore, the exploration of these models is computationally challenging. We propose that integrating two existing technologies—one to aid the construction of multiscale agent-based models, the other developed to enhance model exploration and optimization—can provide a computational means for high-throughput hypothesis testing, and eventually, optimization. Results In this paper, we introduce a high throughput computing (HTC) framework that integrates a mechanistic 3-D multicellular simulator (PhysiCell) with an extreme-scale model exploration platform (EMEWS) to investigate high-dimensional parameter spaces. We show early results in applying PhysiCell-EMEWS to 3-D cancer immunotherapy and show insights on therapeutic failure. We describe a generalized PhysiCell-EMEWS workflow for high-throughput cancer hypothesis testing, where hundreds or thousands of mechanistic simulations are compared against data-driven error metrics to perform hypothesis optimization. Conclusions While key notational and computational challenges remain, mechanistic agent-based models and high-throughput model exploration environments can be combined to systematically and rapidly explore key problems in cancer. These high-throughput computational experiments can improve our understanding of the underlying biology, drive future experiments, and ultimately inform clinical practice.

Colorectal cancer and other cancers often metastasize to the liver in later stages of the disease, contributing significantly to patient death. While the biomechanical properties of the liver parenchyma (normal liver tissue) are known to affect tumor cell behavior in primary and metastatic tumors, the role of these properties in driving or inhibiting metastatic inception remains poorly understood, as are the longer-term multicellular dynamics. This study adopts a multi-model approach to study the dynamics of tumor-parenchyma biomechanical interactions during metastatic seeding and growth. We employ a detailed poroviscoelastic model of a liver lobule to study how micrometastases disrupt flow and pressure on short time scales. Results from short-time simulations in detailed single hepatic lobules motivate constitutive relations and biological hypotheses for a minimal agent-based model of metastatic growth in centimeter-scale tissue over months-long time scales. After a parameter space investigation, we find that the balance of basic tumor-parenchyma biomechanical interactions on shorter time scales (adhesion, repulsion, and elastic tissue deformation over minutes) and longer time scales (plastic tissue relaxation over hours) can explain a broad range of behaviors of micrometastases, without the need for complex molecular-scale signaling. These interactions may arrest the growth of micrometastases in a dormant state and prevent newly arriving cancer cells from establishing successful metastatic foci. Moreover, the simulations indicate ways in which dormant tumors could “reawaken” after changes in parenchymal tissue mechanical properties, as may arise during aging or following acute liver illness or injury. We conclude that the proposed modeling approach yields insight into the role of tumor-parenchyma biomechanics in promoting liver metastatic growth, and advances the longer term goal of identifying conditions to clinically arrest and reverse the course of late-stage cancer.

Metastasis is the leading cause of death in patients with cancer, driving considerable scientific and clinical interest in immunosurveillance of micrometastases. We investigated this process by creating a multiscale mathematical model to study the interactions between the immune system and the progression of micrometastases in general epithelial tissue. We analyzed the parameter space of the model using high-throughput computing resources to generate over 100,000 virtual patient trajectories. We demonstrated that the model could recapitulate a wide variety of virtual patient trajectories, including uncontrolled growth, partial response, and complete immune response to tumor growth. We classified the virtual patients and identified key patient parameters with the greatest effect on the simulated immunosurveillance. We highlight the lessons derived from this analysis and their impact on the nascent field of cancer patient digital twins (CPDTs). While CPDTs could enable clinicians to systematically dissect the complexity of cancer in each individual patient and inform treatment choices, our work shows that key challenges remain before we can reach this vision. In particular, we show that there remain considerable uncertainties in immune responses, unreliable patient stratification, and unpredictable personalized treatment. Nonetheless, we also show that in spite of these challenges, patient-specific models suggest strategies to increase control of clinically undetectable micrometastases even without complete parameter certainty.

Coastal lagoons are dynamic aquatic systems that are susceptible to eutrophication due to long residence times and high inputs of nutrients. We hypothesise that groundwater-derived nutrients make a significant contribution to primary production, eutrophication and carbon cycling in these systems. Here, we report results from seasonal (pre-bloom, bloom and post bloom), nutrient, pCO₂ and ²²²Rn (a natural groundwater tracer) surveys in a coastal lagoon that frequently experiences macroalgae blooms (Avoca Lagoon, Australia). Groundwater inputs of DIN and DIP deliver the nutrient load equivalent to the entire lagoon inventory in 1.6 to 3.5 days and 1.7 and 6 days respectively, indicating groundwater was a major source of inorganic nutrients to the system. Lagoon pCO₂ displayed significant spatial and temporal variability, with the average pCO₂ shifting from 1717 µatm during pre-bloom, to 137 µatm during the bloom and 3056 µatm post bloom. Radon-222 displayed a significant positive relationship with dissolved inorganic nitrogen (DIN) and pCO₂ during the pre-bloom period. This suggests a hydrological dominance of surface water DIN and pCO₂ during the pre-bloom period. During the bloom period, DIN displayed a positive relationship with pCO₂ and negative relationship with dissolved oxygen, indicating a strong biological control over the nutrient pool. Phosphate did not appear to be controlled by either groundwater inputs, or ecosystem metabolism throughout the study. While groundwater discharge stimulated primary production through nutrient inputs (thus reducing surface water CO₂), it also directly delivered significant quantities of CO₂ to surface waters. The net effect of groundwater inputs of nutrients and dissolved carbon into the lagoon was a stimulation of CO₂ fluxes to the atmosphere.

We are rapidly approaching a future in which cancer patient digital twins will reach their potential to predict cancer prevention, diagnosis, and treatment in individual patients. This will be realized based on advances in high performance computing, computational modeling, and an expanding repertoire of observational data across multiple scales and modalities. In 2020, the US National Cancer Institute, and the US Department of Energy, through a trans-disciplinary research community at the intersection of advanced computing and cancer research, initiated team science collaborative projects to explore the development and implementation of predictive Cancer Patient Digital Twins. Several diverse pilot projects were launched to provide key insights into important features of this emerging landscape and to determine the requirements for the development and adoption of cancer patient digital twins. Projects included exploring approaches to using a large cohort of digital twins to perform deep phenotyping and plan treatments at the individual level, prototyping self-learning digital twin platforms, using adaptive digital twin approaches to monitor treatment response and resistance, developing methods to integrate and fuse data and observations across multiple scales, and personalizing treatment based on cancer type. Collectively these efforts have yielded increased insights into the opportunities and challenges facing cancer patient digital twin approaches and helped define a path forward. Given the rapidly growing interest in patient digital twins, this manuscript provides a valuable early progress report of several CPDT pilot projects commenced in common, their overall aims, early progress, lessons learned and future directions that will increasingly involve the broader research community.

Abstract Mangrove ecosystems in Indonesia are facing tremendous threats from anthropogenic activities. The degradation and depletion of mangroves in Indonesia are driven by the rising demand for commodities, including shrimp, rice, palm oil, and other cash crops. These activities ultimately result in negative environmental impacts on mangrove forests, adjacent ecosystems, and the livelihoods of coastal communities. Reliable and accurate information on the losses from mangrove ecosystems is critical to achieving Indonesia’s national climate and blue economy targets. Here, we analyze the effects of human activities on the vulnerability of mangrove forests in Indonesia, addressing the potential loss of blue carbon and other environmental services. We conduct literature review from publications in major bibliographic databases and academic search engines. From this study we found that the impacts of mangrove loss are detrimental to the environment, social and economy. To halt further loss, there is a need for a comprehensive system of mangrove conservation and restoration in Indonesia that could improve local communities’ livelihood. At the national level, Indonesia will need to reform its policies in coastal and marine management, improve information systems, and enable new funding and coordination. Science-based policies and management would lay the groundwork for a sustainable framework for mangrove management and consequently enhance blue carbon governance at the national level.