{"title":"Steady-State Drug Exposure of Repeated IV Bolus Administration for a One Compartment Pharmacokinetic Model with Sigmoidal Hill Elimination.","authors":"Meizhu Cao, Xiaotian Wu, Jun Li","doi":"10.1007/s11538-024-01375-0","DOIUrl":"10.1007/s11538-024-01375-0","url":null,"abstract":"<p><p>Drugs exhibiting nonlinear pharmacokinetics hold significant importance in drug research and development. However, evaluating drug exposure accurately is challenging with the current formulae established for linear pharmacokinetics. This article aims to investigate the steady-state drug exposure for a one-compartment pharmacokinetic (PK) model with sigmoidal Hill elimination, focusing on three key topics: the comparison between steady-state drug exposure of repeated intravenous (IV) bolus ( <math><msub><mtext>AUC</mtext> <mrow><mi>ss</mi></mrow> </msub> </math> ) and total drug exposure after a single IV bolus ( <math><msub><mtext>AUC</mtext> <mrow><mn>0</mn> <mo>-</mo> <mi>∞</mi></mrow> </msub> </math> ); the evolution of steady-state drug concentration with varying dosing frequencies; and the control of drug pharmacokinetics in multiple-dose therapeutic scenarios. For the first topic, we established conditions for the existence of <math><msub><mtext>AUC</mtext> <mrow><mi>ss</mi></mrow> </msub> </math> , derived an explicit formula for its calculation, and compared it with <math><msub><mtext>AUC</mtext> <mrow><mn>0</mn> <mo>-</mo> <mi>∞</mi></mrow> </msub> </math> . For the second, we identified the trending properties of steady-state average and trough concentrations concerning dosing frequency. For the third, we developed formulae to compute dose and dosing time for both regular and irregular dosing scenarios. As an example, our findings were applied to a real drug model of progesterone used in lactating dairy cows. In conclusion, these results provide a theoretical foundation for designing rational dosage regimens and conducting therapeutic trials.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"143"},"PeriodicalIF":2.0,"publicationDate":"2024-11-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142564011","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}
{"title":"Accumulation of Oncogenic Mutations During Progression from Healthy Tissue to Cancer.","authors":"Ruibo Zhang, Ivana Bozic","doi":"10.1007/s11538-024-01372-3","DOIUrl":"10.1007/s11538-024-01372-3","url":null,"abstract":"<p><p>Cancers are typically fueled by sequential accumulation of driver mutations in a previously healthy cell. Some of these mutations, such as inactivation of the first copy of a tumor suppressor gene, can be neutral, and some, like those resulting in activation of oncogenes, may provide cells with a selective growth advantage. We study a multi-type branching process that starts with healthy tissue in homeostasis and models accumulation of neutral and advantageous mutations on the way to cancer. We provide results regarding the sizes of premalignant populations and the waiting times to the first cell with a particular combination of mutations, including the waiting time to malignancy. Finally, we apply our results to two specific biological settings: initiation of colorectal cancer and age incidence of chronic myeloid leukemia. Our model allows for any order of neutral and advantageous mutations and can be applied to other evolutionary settings.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"142"},"PeriodicalIF":4.6,"publicationDate":"2024-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11522190/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142543859","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":"Bounding Seed Loss from Isolated Habitat Patches.","authors":"Benjamin Hafner, Katherine Meyer","doi":"10.1007/s11538-024-01367-0","DOIUrl":"10.1007/s11538-024-01367-0","url":null,"abstract":"<p><p>Dispersal of propagules (seeds, spores) from a geographically isolated habitat into an uninhabitable matrix can play a decisive role in driving population dynamics. ODE and integrodifference models of these dynamics commonly feature a \"dispersal success\" parameter representing the average proportion of dispersing propagules that remain in viable habitat. While dispersal success can be estimated by empirical measurements or by integration of dispersal kernels, one may lack resources for fieldwork or details on dispersal kernels for numerical computation. Here we derive simple upper bounds on the proportion of propagule loss-the complement of dispersal success-that require only habitat area, habitat perimeter, and the mean dispersal distance of a propagule. Using vector calculus in a probabilistic framework, we rigorously prove bounds for the cases of both symmetric and asymmetric dispersal. We compare the bounds to simulations of integral models for the population of Asclepias syriaca (common milkweed) at McKnight Prairie-a 14 hectare reserve surrounded by agricultural fields in Goodhue County, Minnesota-and identify conditions under which the bounds closely estimate propagule loss.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"141"},"PeriodicalIF":4.6,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11519115/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142521065","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":"Evolution of dispersal in river networks.","authors":"Olga Vasilyeva, Dylan Smith, Frithjof Lutscher","doi":"10.1007/s11538-024-01370-5","DOIUrl":"10.1007/s11538-024-01370-5","url":null,"abstract":"<p><p>Evolution of dispersal is a fascinating topic at the intersection of ecology and evolutionary dynamics that has generated many challenging problems in the analysis of reaction-diffusion equations. Early results indicated that lower random diffusion rates are generally beneficial. However, in riverine environments with downstream drift, high diffusion may be optimal, depending on downstream boundary conditions. Most of these results were obtained from modeling a single river reach, yet many rivers form intricate tree-shaped networks. We study the evolution of dispersal on a metric graph representing the simplest such possible network: two upstream segments joining to form one downstream segment. We first show that the shape of the positive steady state of a single population depends crucially on the geometry of the network, here considered as the relative length of the three segments. We then study the evolution of dispersal by considering the possibility of \"invasion\" of a second type (invader) at the steady state of the first type (resident). We show that the geometry of the network determines whether higher or intermediate dispersal is favored.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"140"},"PeriodicalIF":2.0,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142521066","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}
Louis V Kunz, Jesús J Bosque, Mohammad Nikmaneshi, Ibrahim Chamseddine, Lance L Munn, Jan Schuemann, Harald Paganetti, Alejandro Bertolet
{"title":"AMBER: A Modular Model for Tumor Growth, Vasculature and Radiation Response.","authors":"Louis V Kunz, Jesús J Bosque, Mohammad Nikmaneshi, Ibrahim Chamseddine, Lance L Munn, Jan Schuemann, Harald Paganetti, Alejandro Bertolet","doi":"10.1007/s11538-024-01371-4","DOIUrl":"10.1007/s11538-024-01371-4","url":null,"abstract":"<p><p>Computational models of tumor growth are valuable for simulating the dynamics of cancer progression and treatment responses. In particular, agent-based models (ABMs) tracking individual agents and their interactions are useful for their flexibility and ability to model complex behaviors. However, ABMs have often been confined to small domains or, when scaled up, have neglected crucial aspects like vasculature. Additionally, the integration into tumor ABMs of precise radiation dose calculations using gold-standard Monte Carlo (MC) methods, crucial in contemporary radiotherapy, has been lacking. Here, we introduce AMBER, an Agent-based fraMework for radioBiological Effects in Radiotherapy that computationally models tumor growth and radiation responses. AMBER is based on a voxelized geometry, enabling realistic simulations at relevant pre-clinical scales by tracking temporally discrete states stepwise. Its hybrid approach, combining traditional ABM techniques with continuous spatiotemporal fields of key microenvironmental factors such as oxygen and vascular endothelial growth factor, facilitates the generation of realistic tortuous vascular trees. Moreover, AMBER is integrated with TOPAS, an MC-based particle transport algorithm that simulates heterogeneous radiation doses. The impact of radiation on tumor dynamics considers the microenvironmental factors that alter radiosensitivity, such as oxygen availability, providing a full coupling between the biological and physical aspects. Our results show that simulations with AMBER yield accurate tumor evolution and radiation treatment outcomes, consistent with established volumetric growth laws and radiobiological understanding. Thus, AMBER emerges as a promising tool for replicating essential features of tumor growth and radiation response, offering a modular design for future expansions to incorporate specific biological traits.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"139"},"PeriodicalIF":2.0,"publicationDate":"2024-10-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142495576","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}
{"title":"Host-Parasitoid Systems are Vulnerable to Extinction via P-Tipping: Forest Tent Caterpillar as an Example.","authors":"Bryce F Dyck, Rebecca Tyson","doi":"10.1007/s11538-024-01358-1","DOIUrl":"10.1007/s11538-024-01358-1","url":null,"abstract":"<p><p>Continuous-time predator-prey models admit limit cycle solutions that are vulnerable to the phenomenon of phase-sensitive tipping (P-tipping): The predator-prey system can tip to extinction following a rapid change in a key model parameter, even if the limit cycle remains a stable attractor. In this paper, we investigate the existence of P-tipping in an analogous discrete-time system: a host-parasitoid system, using the economically damaging forest tent caterpillar as our motivating example. We take the intrinsic growth rate of the consumer as our key parameter, allowing it to vary with environmental conditions in ways consistent with the predictions of global warming. We find that the discrete-time system does admit P-tipping, and that the discrete-time P-tipping phenomenon shares characteristics with the continuous-time one: Both require an Allee effect on the resource population, occur in small subsets of the phase plane, and exhibit stochastic resonance as a function of the autocorrelation in the environmental variability. In contrast, the discrete-time P-tipping phenomenon occurs when the environmental conditions switch from low to high productivity, can occur even if the magnitude of the switch is relatively small, and can occur from multiple disjoint regions in the phase plane.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"138"},"PeriodicalIF":2.0,"publicationDate":"2024-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142458595","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}
{"title":"Wolbachia Invasion in Mosquitoes with Incomplete CI, Imperfect Maternal Transmission and Maturation Delay.","authors":"Xiaoke Ma, Ying Su","doi":"10.1007/s11538-024-01363-4","DOIUrl":"10.1007/s11538-024-01363-4","url":null,"abstract":"<p><p>The mechanism of cytoplasmic incompatibility (CI) is important in the study of Wolbachia invasion in wild mosquitoes. Su et al. (Bull Math Biol 84(9):95, 2022) proposed a delay differential equation model by relating the CI effect to maturation delay. In this paper, we investigate the dynamics of this model by allowing the same density-dependent death rate and distinct density-independent death rates. Through analyzing the existence and stability of equilibria, we obtain the parameter conditions for Wolbachia successful invasion if the maternal transmission is perfect. While if the maternal transmission is imperfect, we give the ranges of parameters to ensure failure invasion, successful invasion and partially suppressing, respectively. Meanwhile, numerical simulations indicate that the system may exhibit monostable and bistable dynamics when parameters vary. Particularly, in the bistable situation an unstable separatrix, like a line, exists when choosing constant functions as initial values; and the maturation delay affects this separatrix in an interesting way.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 12","pages":"137"},"PeriodicalIF":2.0,"publicationDate":"2024-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142458604","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}
{"title":"Fungal Parasite Transmission in a Planktonic Ecosystem Under Light and Nutrient Constraints.","authors":"Yawen Yan, Juping Ji, Hao Wang","doi":"10.1007/s11538-024-01365-2","DOIUrl":"https://doi.org/10.1007/s11538-024-01365-2","url":null,"abstract":"<p><p>The two main components of the planktonic ecosystem are phytoplankton and zooplankton. Fungal parasites can infect zooplankton and spread between them. In this paper, we construct a dynamic model to describe the spread of fungal parasites among zooplankton. Basic reproduction number for fungal parasite transmission among zooplankton are rigorously derived. The dynamics of this system are analyzed including dissipativity and equilibria. We further explore the effects of ecological factors on population dynamics and the relationship between fungal parasite transmission and phytoplankton blooms. Interestingly, our theoretical and numerical results indicate that a low-light or oligotrophic aquatic environment is helpful in mitigating the transmission of fungal parasites. We also show that fungal parasites on zooplankton can increase phytoplankton biomass and induce blooms.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 11","pages":"136"},"PeriodicalIF":2.0,"publicationDate":"2024-10-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142458594","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}
Stefano Pasetto, Isha Harshe, Renee Brady-Nicholls, Robert A Gatenby, Heiko Enderling
{"title":"Harnessing Flex Point Symmetry to Estimate Logistic Tumor Population Growth.","authors":"Stefano Pasetto, Isha Harshe, Renee Brady-Nicholls, Robert A Gatenby, Heiko Enderling","doi":"10.1007/s11538-024-01361-6","DOIUrl":"10.1007/s11538-024-01361-6","url":null,"abstract":"<p><p>The observed time evolution of a population is well approximated by a logistic growth function in many research fields, including oncology, ecology, chemistry, demography, economy, linguistics, and artificial neural networks. Initial growth is exponential, then decelerates as the population approaches its limit size, i.e., the carrying capacity. In mathematical oncology, the tumor carrying capacity has been postulated to be dynamically evolving as the tumor overcomes several evolutionary bottlenecks and, thus, to be patient specific. As the relative tumor-over-carrying capacity ratio may be predictive and prognostic for tumor growth and treatment response dynamics, it is paramount to estimate it from limited clinical data. We show that exploiting the logistic function's rotation symmetry can help estimate the population's growth rate and carry capacity from fewer data points than conventional regression approaches. We test this novel approach against published pan-cancer animal and human breast cancer data, achieving a 30% to 40% reduction in the time at which subsequent data collection is necessary to estimate the logistic growth rate and carrying capacity correctly. These results could improve tumor dynamics forecasting and augment the clinical decision-making process.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 11","pages":"135"},"PeriodicalIF":2.0,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142388257","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}
{"title":"Modelling Mucus Clearance in Sinuses: Thin-Film Flow Inside a Fluid-Producing Cavity Lined with an Active Surface.","authors":"Nikhil Desai, Eric Lauga","doi":"10.1007/s11538-024-01360-7","DOIUrl":"10.1007/s11538-024-01360-7","url":null,"abstract":"<p><p>The paranasal sinuses are a group of hollow spaces within the human skull, surrounding the nose. They are lined with an epithelium that contains mucus-producing cells and tiny hairlike active appendages called cilia. The cilia beat constantly to sweep mucus out of the sinus into the nasal cavity, thus maintaining a clean mucus layer within the sinuses. This process, called mucociliary clearance, is essential for a healthy nasal environment and disruption in mucus clearance leads to diseases such as chronic rhinosinusitis, specifically in the maxillary sinuses, which are the largest of the paranasal sinuses. We present here a continuum mathematical model of mucociliary clearance inside the human maxillary sinus. Using a combination of analysis and computations, we study the flow of a thin fluid film inside a fluid-producing cavity lined with an active surface: fluid is continuously produced by a wall-normal flux in the cavity and then is swept out, against gravity, due to an effective tangential flow induced by the cilia. We show that a steady layer of mucus develops over the cavity surface only when the rate of ciliary clearance exceeds a threshold, which itself depends on the rate of mucus production. We then use a scaling analysis, which highlights the competition between gravitational retention and cilia-driven drainage of mucus, to rationalise our computational results. We discuss the biological relevance of our findings, noting that measurements of mucus production and clearance rates in healthy sinuses fall within our predicted regime of steady-state mucus layer development.</p>","PeriodicalId":9372,"journal":{"name":"Bulletin of Mathematical Biology","volume":"86 11","pages":"134"},"PeriodicalIF":2.0,"publicationDate":"2024-10-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11455677/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142379124","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}