{"title":"Correction to: Combining Mechanisms of Growth Arrest in Solid Tumours: A Mathematical Investigation.","authors":"Chloé Colson, Helen M Byrne, Philip K Maini","doi":"10.1007/s11538-024-01382-1","DOIUrl":"10.1007/s11538-024-01382-1","url":null,"abstract":"","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"42"},"PeriodicalIF":2.0,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11821782/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143398050","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
John Metzcar, Ben S Duggan, Brandon Fischer, Matthew Murphy, Randy Heiland, Paul Macklin
{"title":"A Simple Framework for Agent-Based Modeling with Extracellular Matrix.","authors":"John Metzcar, Ben S Duggan, Brandon Fischer, Matthew Murphy, Randy Heiland, Paul Macklin","doi":"10.1007/s11538-024-01408-8","DOIUrl":"10.1007/s11538-024-01408-8","url":null,"abstract":"<p><p>Extracellular matrix (ECM) is a key component of the cellular microenvironment and critical in multiple disease and developmental processes. Representing ECM and cell-ECM interactions is a challenging multiscale problem as they span molecular-level details to tissue-level dynamics. While several computational frameworks exist for ECM modeling, they often focus on very detailed modeling of individual ECM fibers or represent only a single aspect of the ECM. Using the PhysiCell agent-based modeling platform, we developed a framework of intermediate detail with the ability to capture bidirectional cell-ECM interactions. We represent a small region of ECM, an ECM element, with three variables describing its local microstructure: anisotropy, density, and overall fiber orientation. To spatially model the ECM, we use an array of ECM elements. Cells remodel local ECM microstructure and in turn, local microstructure impacts cellular motility. We demonstrate the utility of this framework and reusability of its core cell-ECM interaction model through examples in cellular invasion, wound healing, basement membrane degradation, and leader-follower collective migration. Despite the relative simplicity of the framework, it is able to capture a broad range of cell-ECM interactions of interest to the modeling community. Furthermore, variables representing the ECM microstructure are accessible through simple programming interfaces. This allows them to impact cell behaviors, such as proliferation and death, without requiring custom code for each interaction, particularly through PhysiCell's modeling grammar, enabling rapid modeling of a diverse range of cell-matrix biology. We make this framework available as a free and open source software package at https://github.com/PhysiCell-Models/collective-invasion .</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"43"},"PeriodicalIF":2.0,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11821717/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143397890","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Simplifying and Characterizing DAGs and Phylogenetic Networks via Least Common Ancestor Constraints.","authors":"Anna Lindeberg, Marc Hellmuth","doi":"10.1007/s11538-025-01419-z","DOIUrl":"10.1007/s11538-025-01419-z","url":null,"abstract":"<p><p>Rooted phylogenetic networks, or more generally, directed acyclic graphs (DAGs), are widely used to model species or gene relationships that traditional rooted trees cannot fully capture, especially in the presence of reticulate processes or horizontal gene transfers. Such networks or DAGs are typically inferred from observable data (e.g., genomic sequences of extant species), providing only an estimate of the true evolutionary history. However, these inferred DAGs are often complex and difficult to interpret. In particular, many contain vertices that do not serve as least common ancestors (LCAs) for any subset of the underlying genes or species, thus may lack direct support from the observable data. In contrast, LCA vertices are witnessed by historical traces justifying their existence and thus represent ancestral states substantiated by the data. To reduce unnecessary complexity and eliminate unsupported vertices, we aim to simplify a DAG to retain only LCA vertices while preserving essential evolutionary information. In this paper, we characterize <math><mtext>LCA</mtext></math> -relevant and <math><mtext>lca</mtext></math> -relevant DAGs, defined as those in which every vertex serves as an LCA (or unique LCA) for some subset of taxa. We introduce methods to identify LCAs in DAGs and efficiently transform any DAG into an <math><mtext>LCA</mtext></math> -relevant or <math><mtext>lca</mtext></math> -relevant one while preserving key structural properties of the original DAG or network. This transformation is achieved using a simple operator \" <math><mo>⊖</mo></math> \" that mimics vertex suppression.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"44"},"PeriodicalIF":2.0,"publicationDate":"2025-02-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11821803/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143398056","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Effects of Non-linear Interaction Between Oxygen and Lactate on Solid Tumor Growth Under Cyclic Hypoxia.","authors":"Gopinath Sadhu, D C Dalal","doi":"10.1007/s11538-025-01420-6","DOIUrl":"10.1007/s11538-025-01420-6","url":null,"abstract":"<p><p>Oxygen is a crucial element for cellular respiration. Based on oxygen concentration, tumor regions can be categorized as normoxic, hypoxic, and necrotic. Hypoxic tumor cells switch their metabolism from aerobic glycolysis to anaerobic glycolysis. As a result, lactate is produced in hypoxic regions and is used as an alternative metabolic fuel by normoxic tumor cells. The consumption of lactate and oxygen by tumor cells does not follow a linear pattern. Scientific studies suggest that oxygen consumption and lactate production are non-linear phenomena. In this study, we propose a two-dimensional mathematical model to investigate lactate dynamics in avascular tumors with various initial shapes, such as circular, elliptical, and petal, and to explore its growth patterns in the context of non-linear interactions between oxygen and lactate. In certain human tumors, particularly in kidney, skin, and liver, multiple tumors may emerge within a tissue domain simultaneously. We also examine how the growth patterns of multiple tumors evolve within a shared domain. Cyclic hypoxia, commonly observed in solid tumors, results from oxygen fluctuations over time at the tumor site. Additionally, we analyze lactate dynamics and tumor growth patterns in environments with cyclic hypoxia. In order to simulate the proposed model, we use finite element based COMSOL Multiphysics 6.0 interface. The simulated results show excellent agreement with experimental data. Our findings reveal that the initial tumor shape significantly influences the lactate distribution and the tumor's internal structure. Furthermore, the simulations indicate that multiple tumors eventually merge into a single tumor. We also observe that cyclic hypoxia with short periodicity increases tumor volume.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"41"},"PeriodicalIF":2.0,"publicationDate":"2025-02-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143398055","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Soukaina Sabir, Odelaisy León-Triana, Sergio Serrano, Roberto Barrio, Victor M Pérez-García
{"title":"Mathematical Model of CAR T-Cell Therapy for a B-Cell Lymphoma Lymph Node.","authors":"Soukaina Sabir, Odelaisy León-Triana, Sergio Serrano, Roberto Barrio, Victor M Pérez-García","doi":"10.1007/s11538-025-01417-1","DOIUrl":"10.1007/s11538-025-01417-1","url":null,"abstract":"<p><p>CAR T-cell therapies have demonstrated significant success in treating B-cell leukemia in children and young adults. However, their effectiveness in treating B-cell lymphomas has been limited in comparison to leukemia. In this paper we present a mathematical model that elucidates the dynamics of diffuse large B-cell lymphoma and CAR T-cells in a lymph node. The mathematical model aids in understanding the complex interplay between the cell populations involved and proposes ways to identify potential underlying dynamical causes of treatment failure. We also study the phenomenon of immunosuppression induced by tumor cells and theoretically demonstrate its impact on cell dynamics. Through the examination of various response scenarios, we underscore the significance of product characteristics in treatment outcomes.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"40"},"PeriodicalIF":2.0,"publicationDate":"2025-02-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11805830/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143363717","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Employing Observability Rank Conditions for Taking into Account Experimental Information a priori.","authors":"Alejandro F Villaverde","doi":"10.1007/s11538-025-01415-3","DOIUrl":"10.1007/s11538-025-01415-3","url":null,"abstract":"<p><p>The concept of identifiability describes the possibility of inferring the parameters of a dynamic model by observing its output. It is common and useful to distinguish between structural and practical identifiability. The former property is fully determined by the model equations, while the latter is also influenced by the characteristics of the available experimental data. Structural identifiability can be determined by means of symbolic computations, which may be performed before collecting experimental data, and are hence sometimes called a priori analyses. Practical identifiability is typically assessed numerically, with methods that require simulations-and often also optimization-and are applied a posteriori. An approach to study structural local identifiability is to consider it as a particular case of observability, which is the possibility of inferring the internal state of a system from its output. Thus, both properties can be analysed jointly, by building a generalized observability matrix and computing its rank. The aim of this paper is to investigate to which extent such observability-based methods can also inform about practical aspects related with the experimental setup, which are usually not approached in this way. To this end, we explore a number of possible extensions of the rank tests, and discuss the purposes for which they can be informative as well as others for which they cannot.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"39"},"PeriodicalIF":2.0,"publicationDate":"2025-02-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11802610/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143254620","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Victoria Chebotaeva, Anish Srinivasan, Paula A Vasquez
{"title":"Differentiating Contact with Symptomatic and Asymptomatic Infectious Individuals in a SEIR Epidemic Model.","authors":"Victoria Chebotaeva, Anish Srinivasan, Paula A Vasquez","doi":"10.1007/s11538-025-01416-2","DOIUrl":"10.1007/s11538-025-01416-2","url":null,"abstract":"<p><p>This manuscript introduces a new Erlang-distributed SEIR model. The model incorporates asymptomatic spread through a subdivided exposed class, distinguishing between asymptomatic ( <math><msub><mtext>E</mtext> <mi>a</mi></msub> </math> ) and symptomatic ( <math><msub><mtext>E</mtext> <mi>s</mi></msub> </math> ) cases. The model identifies two key parameters: relative infectiousness, <math><msub><mi>β</mi> <mrow><mi>SA</mi></mrow> </msub> </math> , and the percentage of people who become asymptomatic after being infected by a symptomatic individual, <math><mi>κ</mi></math> . Lower values of these parameters reduce the peak magnitude and duration of the infectious period, highlighting the importance of isolation measures. Additionally, the model underscores the need for strategies addressing both symptomatic and asymptomatic transmissions.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"38"},"PeriodicalIF":2.0,"publicationDate":"2025-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11794362/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143188115","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Matthew C Nitschke, Viney Kumar, Katrina E Milliner, Kristen Hawkes, Peter S Kim
{"title":"Evolution of Human Pair Bonds as a Consequence of Male-Biased Mating Sex Ratios?","authors":"Matthew C Nitschke, Viney Kumar, Katrina E Milliner, Kristen Hawkes, Peter S Kim","doi":"10.1007/s11538-025-01414-4","DOIUrl":"10.1007/s11538-025-01414-4","url":null,"abstract":"<p><p>Compared to our closest primate relatives, human life history involves greater longevity, which includes a distinctive postmenopausal life stage. Given mammalian reproductive physiology in which females build a finite stock of cells that can become oocytes early in life, which then continuously deplete mostly through cell death while males produce new sperm throughout adulthood, the postmenopausal stage makes the sex ratio in the fertile pool, called the adult sex ratio (ASR), male biased. Additionally, this affects a more fine-grained ratio, the operational sex ratio (OSR), defined as the ratio of males to females currently able to conceive. Here, we construct an ODE model in which males compete for paternities using either a multiple-mating or mate-guarding strategy. Our focus is on investigating the differences of strategy choice between populations with varying life histories, which include a distinct post-fertile stage for adult females. By simulating the system, we determine the dominant strategy and its dependence on various parameter combinations. Our results show that an increase in OSR and ASR correlates well with a change in the dominant strategy from multiple mating to guarding.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"37"},"PeriodicalIF":2.0,"publicationDate":"2025-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11782467/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143064006","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Adaptive Immunity Determines the Cancer Treatment Outcome of Oncolytic Virus and Anti-PD-1.","authors":"Kang-Ling Liao, Kenton D Watt","doi":"10.1007/s11538-025-01413-5","DOIUrl":"10.1007/s11538-025-01413-5","url":null,"abstract":"<p><p>The immune checkpoint inhibitor, anti-programmed death protein-1 (anti-PD-1), enhances adaptive immunity to kill tumor cells, and the oncolytic virus (OV) triggers innate immunity to clear the infected tumor cells. We create a mathematical model to investigate how the interaction between adaptive and innate immunities under OV and anti-PD-1 affects tumor reduction. For different immunity strength, we create the corresponding virtual baseline patients and cohort patients to decipher the major factors determining the treatment outcome. Global sensitivity analysis indicates that adaptive immunity has more control on the treatment outcome than innate immunity, and whether anti-PD-1 cancels out the OV treatment efficacy depends on the OV dosage and the balance between clearance of infected tumor cells and OV by T cells. The optimal OV infection rate and dosage suggest that OV treatment is more sensitive to adaptive immunity than innate immunity. Our model prediction also indicates that tumor reduction is more sensitive to anti-PD-1 efficacy as adaptive immunity becomes stronger, and anti-PD-1 trends to cancel out the OV treatment efficacy as innate immunity becomes stronger. Based on these results, the recommended treatment protocol for patients with different immunity strength can be determined.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 3","pages":"36"},"PeriodicalIF":2.0,"publicationDate":"2025-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143058219","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
George Atkinson, Yaron Ben-Ami, Philip Maini, Joe Pitt-Francis, Helen Byrne
{"title":"Key Structural Features of Microvascular Networks Leading to the Formation of Multiple Equilibria.","authors":"George Atkinson, Yaron Ben-Ami, Philip Maini, Joe Pitt-Francis, Helen Byrne","doi":"10.1007/s11538-024-01404-y","DOIUrl":"10.1007/s11538-024-01404-y","url":null,"abstract":"<p><p>We analyse mathematical models of blood flow in two simple vascular networks in order to identify structural features that lead to the formation of multiple equilibria. Our models are based on existing rules for blood rheology and haematocrit splitting. By performing bifurcation analysis on these simple network flow models, we identify a link between the changing flow direction in key vessels and the existence of multiple equilibria. We refer to these key vessels as redundant vessels, and relate the maximum number of equilibria with the number of redundant vessels. We vary geometric parameters of the two networks, such as vessel length ratios and vessel diameters, to demonstrate that equilibria are uniquely defined by the flow in the redundant vessels. Equilibria typically emerge in sets of three, each having a different flow characteristic in one of the network's redundant vessels. For one of the three equilibria, the flow within the relevant redundant vessel will be smaller in magnitude than the other two and the redundant vessel will contain few Red Blood Cells (RBCs), if any. For the other two equilibria, the redundant vessel contains RBCs and significant flow in the two available directions. These structural features of networks provide a useful geometric property when studying the equilibria of blood flow in microvascular networks.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"87 2","pages":"30"},"PeriodicalIF":2.0,"publicationDate":"2025-01-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11757897/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143022274","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}