1
|
Ladib M, Ouhinou A, Yakubu AA. Mathematical modeling of contact tracing and stability analysis to inform its impact on disease outbreaks; an application to COVID-19. Infect Dis Model 2024; 9:329-353. [PMID: 38371875 PMCID: PMC10867662 DOI: 10.1016/j.idm.2024.01.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2023] [Revised: 01/29/2024] [Accepted: 01/29/2024] [Indexed: 02/20/2024] Open
Abstract
We develop a mathematical model to investigate the effect of contact tracing on containing epidemic outbreaks and slowing down the spread of transmissible diseases. We propose a discrete-time epidemic model structured by disease-age which includes general features of contact tracing. The model is fitted to data reported for the early spread of COVID-19 in South Korea, Brazil, and Venezuela. The calibrated values for the contact tracing parameters reflect the order pattern observed in its performance intensity within the three countries. Using the fitted values, we estimate the effective reproduction number R e and investigate its responses to varied control scenarios of contact tracing. Alongside the positivity of solutions, and a stability analysis of the disease-free equilibrium are provided.
Collapse
Affiliation(s)
- Mohamed Ladib
- University of Sultan Moulay Slimane, Faculty of Sciences and Techniques, Team of Mathematics and Interactions, Béni-Mellal, Morocco
| | - Aziz Ouhinou
- University of Sultan Moulay Slimane, Faculty of Sciences and Techniques, Team of Mathematics and Interactions, Béni-Mellal, Morocco
| | | |
Collapse
|
2
|
Dumont Y, Yatat-Djeumen IV. About contamination by sterile females and residual male fertility on the effectiveness of the sterile insect technique. Impact on disease vector control and disease control. Math Biosci 2024; 370:109165. [PMID: 38387836 DOI: 10.1016/j.mbs.2024.109165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 11/19/2023] [Accepted: 02/13/2024] [Indexed: 02/24/2024]
Abstract
The sterile insect technique (SIT) is a technique to control pests and vectors of diseases by releasing mainly sterile males. Several challenges need to be solved before large-scale field application in order to guarantee its success. In this paper we intend to focus on two important issues: residual fertility in released (sterile) males and contamination of each release by sterile females. Indeed, sterile males are never 100% sterile, that is there is always a small proportion, ɛ, of fertile males (sperm of) within the sterile males population. Among the sterile insects that are released, a certain proportion, ϵF, of them are sterile females due to imperfect mechanical sex-separation technique. This can be particularly problematic when arthropod viruses are circulating, because mosquito females, even sterile, are vectors of diseases. Various upper bound values are given in the entomological literature for ϵF and ɛ without clear explanations. In this work, we aim to show that these values are related to the biological parameters of the targeted vector, the sterile insects release rate, and the epidemiological parameters of a vector-borne disease, like Dengue. We extend results studied separately in Aronna and Dumont (2020), Dumont and Yatat-Djeumen (2022). To study the impact of both issues, we develop and study a SIT-entomological-epidemiological mathematical model, with application to Dengue. Qualitative analysis of the model is carried out to highlight threshold values that shape the overall dynamics of the system. We show that vector elimination is possible only when Nɛ<1, where N is the basic-offspring number related to the targeted wild population. To ensure the success of SIT control, we recommend that the issue of residual fertility be addressed as a priority and then that contamination by sterile females be minimized with each release.
Collapse
Affiliation(s)
- Y Dumont
- CIRAD, Umr AMAP, Pôle de Protection des Plantes, F-97410 Saint Pierre, France; AMAP, Univ Montpellier, CIRAD, CNRS, INRA, IRD, Montpellier, France; University of Pretoria, Department of Mathematics and Applied Mathematics, Pretoria, South Africa.
| | - I V Yatat-Djeumen
- CIRAD, Umr AMAP, F-34398 Montpellier, France; University of Yaoundé I, National Advanced School of Engineering of Yaoundé, Department of Mathematics and Physics, Yaoundé, Cameroon; UMI 209 IRD/UPMC UMMISCO, Bondy, France
| |
Collapse
|
3
|
Singh S, Bianchi L, Korganbayev S, Namakshenas P, Melnik R, Saccomandi P. Non-Fourier Bioheat Transfer Analysis in Brain Tissue During Interstitial Laser Ablation: Analysis of Multiple Influential Factors. Ann Biomed Eng 2024; 52:967-981. [PMID: 38236341 DOI: 10.1007/s10439-023-03433-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Accepted: 12/22/2023] [Indexed: 01/19/2024]
Abstract
This work presents the dual-phase lag-based non-Fourier bioheat transfer model of brain tissue subjected to interstitial laser ablation. The finite element method has been utilized to predict the brain tissue's temperature distributions and ablation volumes. A sensitivity analysis has been conducted to quantify the effect of variations in the input laser power, treatment time, laser fiber diameter, laser wavelength, and non-Fourier phase lags. Notably, in this work, the temperature-dependent thermal properties of brain tissue have been considered. The developed model has been validated by comparing the temperature obtained from the numerical and ex vivo brain tissue during interstitial laser ablation. The ex vivo brain model has been further extended to in vivo settings by incorporating the blood perfusion effects. The results of the systematic analysis highlight the importance of considering temperature-dependent thermal properties of the brain tissue, non-Fourier behavior, and microvascular perfusion effects in the computational models for accurate predictions of the treatment outcomes during interstitial laser ablation, thereby minimizing the damage to surrounding healthy tissue. The developed model and parametric analysis reported in this study would assist in a more accurate and precise prediction of the temperature distribution, thus allowing to optimize the thermal dosage during laser therapy in the brain.
Collapse
Affiliation(s)
- Sundeep Singh
- Faculty of Sustainable Design Engineering, University of Prince Edward Island, Charlottetown, PE, C1A 4P3, Canada
| | - Leonardo Bianchi
- Department of Mechanical Engineering, Politecnico di Milano, 20156, Milan, Italy
| | - Sanzhar Korganbayev
- Department of Mechanical Engineering, Politecnico di Milano, 20156, Milan, Italy
| | - Pouya Namakshenas
- Department of Mechanical Engineering, Politecnico di Milano, 20156, Milan, Italy
| | - Roderick Melnik
- MS2Discovery Interdisciplinary Research Institute, Wilfrid Laurier University, Waterloo, ON, Canada
| | - Paola Saccomandi
- Department of Mechanical Engineering, Politecnico di Milano, 20156, Milan, Italy.
| |
Collapse
|
4
|
Kuznetsov AV. Lewy body radius growth: The hypothesis of the cube root of time dependency. J Theor Biol 2024; 581:111734. [PMID: 38246486 DOI: 10.1016/j.jtbi.2024.111734] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 01/15/2024] [Indexed: 01/23/2024]
Abstract
This paper presents a model for the growth of Lewy bodies (LBs), which are pathological hallmarks of Parkinson's disease (PD). The model simulates the growth of classical LBs, consisting of a core and a halo. The core is assumed to comprise lipid membrane fragments and damaged organelles, while the halo consists of radiating alpha-synuclein (α-syn) fibrils. The Finke-Watzky model is employed to simulate the aggregation of lipid fragments and α-syn monomers. Analytical and numerical exploration of the governing equations yielded approximate solutions applicable for larger times. The application of these approximate solutions to simulate LB radius growth led to the discovery of the cube root hypothesis, which posits that the LB radius is proportional to the cube root of its growth time. Sensitivity analysis revealed that the LB radius is unaffected by the kinetic rates of nucleation and autocatalytic growth, with growth primarily regulated by the production rates of lipid membrane fragments and α-syn monomers. The model indicates that the formation of large LBs associated with PD is dependent on the malfunction of the machinery responsible for the degradation of lipid membrane fragments, α-syn monomers, and their aggregates.
Collapse
Affiliation(s)
- Andrey V Kuznetsov
- Department of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC 27695-7910, USA.
| |
Collapse
|
5
|
Karametos I, Vasiliadou IA, Papaevangelou V, Sultana MY, Tekerlekopoulou AG, Vayenas DV, Akratos CS. Mathematical modeling of constructed wetlands for hexavalent chromium removal. Sci Total Environ 2024; 917:170088. [PMID: 38220001 DOI: 10.1016/j.scitotenv.2024.170088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 12/24/2023] [Accepted: 01/09/2024] [Indexed: 01/16/2024]
Abstract
Constructed wetlands (CWs) have been extensively used in Cr(VI) removal and have proven their ability to achieve high removal efficiencies. Although, numerous studies have been published in the past years presenting experimental results of CWs treating wastewater with Cr(VI) concentrations, a mathematical modeling describing the processes for Cr(VI) removal in CWs is lacking. In this work a mathematical model was developed, able to accurately describe the main mechanisms and reactions (i.e. biological reduction, plant biomass uptake-sorption) which are responsible for Cr(VI) removal in a wetland system. The model was calibrated and validated using data from a previously reported experimental study of horizontal subsurface CWs. Mathematical simulation indicates that in an unplanted wetland Cr(VI) was mainly removed through the diffusion/reduction process inside biofilm, attached on the porous media, while in the planted unit Cr(VI) was mainly removed through the sorption process to the root system of the plants. The developed model's simulations showed high correlation between predicted and experimental data, indicating that the proposed model can be used to design and predict full scale constructed wetland process for Cr(VI) removal.
Collapse
Affiliation(s)
- Ioannis Karametos
- Department of Civil Engineering, Democritus University of Thrace, Xanthi 67100, Greece
| | - Ioanna A Vasiliadou
- Department of Chemical Engineering, Western University of Macedonia, Kozani 50100, Greece.
| | | | - Mar-Yam Sultana
- Department of Sustainable Agriculture, University of Patras, 2 G. Seferi Str., Agrinio 30100, Greece
| | | | - Dimitrios V Vayenas
- Department of Chemical Engineering, University of Patras, Rio, GR-26504 Patras, Greece; Institute of Chemical Engineering Sciences (FORTH/ICE-HT), Stadiou Str., Platani, 26504 Patras, Greece
| | - Christos S Akratos
- Department of Civil Engineering, Democritus University of Thrace, Xanthi 67100, Greece
| |
Collapse
|
6
|
Yan PF, Dong S, Woodcock MJ, Manz KE, Garza-Rubalcava U, Abriola LM, Pennell KD, Cápiro NL. Biotransformation of 6:2 fluorotelomer sulfonate and microbial community dynamics in water-saturated one-dimensional flow-through columns. Water Res 2024; 252:121146. [PMID: 38306753 DOI: 10.1016/j.watres.2024.121146] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 01/10/2024] [Accepted: 01/14/2024] [Indexed: 02/04/2024]
Abstract
Nearly all per- and polyfluoroalkyl substances (PFAS) biotransformation studies reported to date have been limited to laboratory-scale batch reactors. The fate and transport of PFAS in systems that more closely represent field conditions, i.e., in saturated porous media under flowing conditions, remain largely unexplored. This study investigated the biotransformation of 6:2 fluorotelomer sulfonate (6:2 FTS), a representative PFAS of widespread environmental occurrence, in one-dimensional water-saturated flow-through columns packed with soil obtained from a PFAS-contaminated site. The 305-day column experiments demonstrated that 6:2 FTS biotransformation was rate-limited, where a decrease in pore-water velocity from 3.7 to 2.4 cm/day, resulted in a 21.7-26.1 % decrease in effluent concentrations of 6:2 FTS and higher yields (1.0-1.4 mol% vs. 0.3 mol%) of late-stage biotransformation products (C4C7 perfluoroalkyl carboxylates). Flow interruptions (2 and 7 days) were found to enhance 6:2 FTS biotransformation during the 6-7 pore volumes following flow resumption. Model-fitted 6:2 FTS column biotransformation rates (0.039-0.041 cmw3/gs/d) were ∼3.5 times smaller than those observed in microcosms (0.137 cmw3/gs/d). Additionally, during column experiments, planktonic microbial communities remained relatively stable, whereas the composition of the attached microbial communities shifted along the flow path, which may have been attributed to oxygen availability and the toxicity of 6:2 FTS and associated biotransformation products. Genus Pseudomonas dominated in planktonic microbial communities, while in the attached microbial communities, Rhodococcus decreased and Pelotomaculum increased along the flow path, suggesting their potential involvement in early- and late-stage 6:2 FTS biotransformation, respectively. Overall, this study highlights the importance of incorporating realistic environmental conditions into experimental systems to obtain a more representative assessment of in-situ PFAS biotransformation.
Collapse
Affiliation(s)
- Peng-Fei Yan
- Department of Biological and Environmental Engineering, Cornell University, Ithaca, NY, United States.
| | - Sheng Dong
- Department of Biological and Environmental Engineering, Cornell University, Ithaca, NY, United States
| | | | - Katherine E Manz
- School of Public Health, University of Michigan, Ann Arbor, MI, United States
| | | | - Linda M Abriola
- School of Engineering, Brown University, Providence, RI, United States
| | - Kurt D Pennell
- School of Engineering, Brown University, Providence, RI, United States
| | - Natalie L Cápiro
- Department of Biological and Environmental Engineering, Cornell University, Ithaca, NY, United States.
| |
Collapse
|
7
|
Henriot P, Buelow E, Petit F, Ploy MC, Dagot C, Opatowski L. Modeling the impact of urban and hospital eco-exposomes on antibiotic-resistance dynamics in wastewaters. Sci Total Environ 2024; 924:171643. [PMID: 38471588 DOI: 10.1016/j.scitotenv.2024.171643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Revised: 01/10/2024] [Accepted: 03/09/2024] [Indexed: 03/14/2024]
Abstract
The emergence and selection of antibiotic resistance is a major public health problem worldwide. The presence of antibiotic-resistant bacteria (ARBs) in natural and anthropogenic environments threatens the sustainability of efforts to reduce resistance in human and animal populations. Here, we use mathematical modeling of the selective effect of antibiotics and contaminants on the dynamics of bacterial resistance in water to analyze longitudinal spatio-temporal data collected in hospital and urban wastewater between 2012 and 2015. Samples were collected monthly during the study period at four different sites in Haute-Savoie, France: hospital and urban wastewater, before and after water treatment plants. Three different categories of exposure variables were collected simultaneously: 1) heavy metals, 2) antibiotics and 3) surfactants for a total of 13 drugs/molecules; in parallel to the normalized abundance of 88 individual genes and mobile genetic elements, mostly conferring resistance to antibiotics. A simple hypothesis-driven model describing weekly antibiotic resistance gene (ARG) dynamics was proposed to fit the available data, assuming that normalized gene abundance is proportional to antibiotic resistant bacteria (ARB) populations in water. The detected compounds were found to influence the dynamics of 17 genes found at multiple sites. While mercury and vancomycin were associated with increased ARG and affected the dynamics of 10 and 12 identified genes respectively, surfactants antagonistically affected the dynamics of three genes. The models proposed here make it possible to analyze the relationship between the persistence of resistance genes in the aquatic environment and specific compounds associated with human activities from longitudinal data. Our analysis of French data over 2012-2015 identified mercury and vancomycin as co-selectors for some ARGs.
Collapse
Affiliation(s)
- Paul Henriot
- Epidemiology and Modeling of bacterial Evasion to Antibacterials Unit (EMEA), Institut Paris, France; MESuRS Laboratory, Conservatoire National des Arts et Métiers Paris, France; Université Paris-Saclay, UVSQ, Inserm, CESP, Anti-Infective Evasion and Pharmacoepidemiology Team, Montigny-le-Bretonneux, France.
| | - Elena Buelow
- Université Limoges, INSERM, CHU Limoges, RESINFIT, U1092 Limoges, France; Univ. Grenoble Alpes, CNRS, UMR 5525, VetAgro Sup, Grenoble INP, TIMC, 38000 Grenoble, France
| | - Fabienne Petit
- UNIROUEN, UNICAEN, CNRS, M2C, Normandie Université, Rouen, France; Sorbonne Université, CNRS, EPHE, PSL, UMR METIS, Paris, France
| | - Marie-Cécile Ploy
- Université Limoges, INSERM, CHU Limoges, RESINFIT, U1092 Limoges, France
| | - Christophe Dagot
- Université Limoges, INSERM, CHU Limoges, RESINFIT, U1092 Limoges, France
| | - Lulla Opatowski
- Epidemiology and Modeling of bacterial Evasion to Antibacterials Unit (EMEA), Institut Paris, France; Université Paris-Saclay, UVSQ, Inserm, CESP, Anti-Infective Evasion and Pharmacoepidemiology Team, Montigny-le-Bretonneux, France
| |
Collapse
|
8
|
Plyusnina TY, Khruschev SS, Degtereva NS, Voronova EN, Volgusheva AA, Riznichenko GY, Rubin AB. Three-state mathematical model for the assessment of DCMU-treated photosystem II heterogeneity. Photosynth Res 2024:10.1007/s11120-024-01077-7. [PMID: 38466456 DOI: 10.1007/s11120-024-01077-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2023] [Accepted: 01/15/2024] [Indexed: 03/13/2024]
Abstract
Photosystem II (PSII) is one of the main pigment-protein complexes of photosynthesis which is highly sensitive to unfavorable environmental factors. The heterogeneity of PSII properties is essential for the resistance of autotrophic organisms to stress factors. Assessment of the PSII heterogeneity may be used in environmental monitoring for on-line detection of contamination of the environment. We propose an approach to assess PSII oxygen-evolving complex and light-harvesting antenna heterogeneity that is based on mathematical modeling of the shape of chlorophyll a fluorescence rise of 3-(3,4-dichlorophenyl)-1,1-dimethylurea-treated samples. The hierarchy of characteristic times of the processes considered in the model makes it possible to reduce the model to a system of three ordinary differential equations. The analytic solution of the reduced three-state model is expressed as a sum of two exponential functions, and it exactly reproduces the solution of the complete system within the time range from microseconds to hundreds of milliseconds. The combination of several such models for reaction centers with different properties made it possible to use it as an instrument to study PSII heterogeneity. PSII heterogeneity was studied for Chlamydomonas at different intensities of actinic light, for Scenedesmus under short-term heating, and for Chlorella grown in nitrate-enriched and nitrate-depleted media.
Collapse
Affiliation(s)
- Tatiana Yu Plyusnina
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia.
| | - Sergei S Khruschev
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
| | - Natalia S Degtereva
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
| | - Elena N Voronova
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
| | - Alena A Volgusheva
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
| | - Galina Yu Riznichenko
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
| | - Andrew B Rubin
- Department of Biophysics, Faculty of Biology, Lomonosov Moscow State University, Moscow, 119234, Russia
| |
Collapse
|
9
|
Sewer A, Talikka M, Calvino-Martin F, Luettich K, Iskandar A. Quantitative modeling of in vitro data using an adverse outcome pathway for the risk assessment of decreased lung function in humans. Toxicol Lett 2024; 393:107-113. [PMID: 38350531 DOI: 10.1016/j.toxlet.2024.02.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 01/18/2024] [Accepted: 02/07/2024] [Indexed: 02/15/2024]
Abstract
In the absence of epidemiological data, there is a need to develop computational models that convert in vitro findings to human disease risk predictions following toxicant exposure. In such efforts, in vitro data can be evaluated in the context of adverse outcome pathways (AOPs) that organize mechanistic knowledge based on empirical evidence into a sequence of molecular-, cellular-, tissue-, and organ-level key events that precede an adverse outcome (AO). Here we combined data from advanced in vitro organotypic airway models exposed to combustible cigarette (CC) smoke or Tobacco Heating System (THS) aerosol with an AOP for increased oxidative stress leads to decreased lung function. The mathematical modeling predicted reduced risk of decreased ciliary beating frequency (CBF) based on oxidative stress measurements and reduced risk of decreased mucociliary clearance (MCC) based on CBF measurements in THS aerosol- compared with CC smoke-exposed cultures. To extend the predictions to the AO of decreased lung function, we leveraged human MCC data from current smokers, nonsmokers, former smokers, and users of heated tobacco products. This approach provided a plausible prediction of diminished reduction in lung function in response to THS use compared with continued smoking. The current approach may also present a basis for an integrated approach to testing and assessment of tobacco products for future regulatory decision-making.
Collapse
Affiliation(s)
- Alain Sewer
- PMI R&D, Philip Morris Products S.A., Neuchâtel, Switzerland.
| | - Marja Talikka
- PMI R&D, Philip Morris Products S.A., Neuchâtel, Switzerland
| | | | - Karsta Luettich
- PMI R&D, Philip Morris Products S.A., Neuchâtel, Switzerland
| | - Anita Iskandar
- PMI R&D, Philip Morris Products S.A., Neuchâtel, Switzerland
| |
Collapse
|
10
|
Enderling H. Fractal calculus in tumor growth simulations: The proof is in the pudding. Biosystems 2024; 237:105141. [PMID: 38355079 DOI: 10.1016/j.biosystems.2024.105141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Revised: 02/05/2024] [Accepted: 02/06/2024] [Indexed: 02/16/2024]
Abstract
Mathematical modeling in oncology has a long history. Recently, mathematical models and their predictions have made inroads into prospective clinical trials with encouraging results. The goal of many such modeling efforts is to make predictions, either to clinician's choice therapy or into "optimal" therapy - often for individual patients. The mathematical oncology community rightfully puts great hope into predictive modeling and mechanistic digital twins - but with this great opportunity comes great responsibility. Mathematical models need to be rigorously calibrated and validated, and their predictive performance ascertained, before conclusions about predictions into the unknown can be drawn. The recent article "Modeling tumor growth using fractal calculus: Insights into tumor dynamics" (Golmankhaneh et al., 2023), applied fractal calculus to tumor growth data. In this short commentary, I raise concerns about the study design and interpretation. In its current form, this study is poised to put cancer patients at risk if interpreted as concluded by the authors.
Collapse
Affiliation(s)
- Heiko Enderling
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA; Institute for Data Science in Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA.
| |
Collapse
|
11
|
Fauzi IS, Nuraini N, Sari AM, Wardani IB, Taurustiati D, Simanullang PM, Lestari BW. Assessing the impact of booster vaccination on diphtheria transmission: Mathematical modeling and risk zone mapping. Infect Dis Model 2024; 9:245-262. [PMID: 38312350 PMCID: PMC10837633 DOI: 10.1016/j.idm.2024.01.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2023] [Revised: 12/23/2023] [Accepted: 01/11/2024] [Indexed: 02/06/2024] Open
Abstract
The COVID-19 pandemic caused significant disruptions in the healthcare system, affecting vaccinations and the management of diphtheria cases. As a consequence of these disruptions, numerous countries have experienced a resurgence or an increase in diphtheria cases. West Java province in Indonesia is identified as one of the high-risk areas for diphtheria, experiencing an upward trend in cases from 2021 to 2023. To analyze the situation, we developed an SIR model, which integrated DPT and booster vaccinations to determine the basic reproduction number, an essential parameter for infectious diseases. Through spatial analysis of geo-referenced data, we identified hotspots and explained diffusion in diphtheria case clusters. The calculation of R0 resulted in an R0 = 1.17, indicating the potential for a diphtheria outbreak in West Java. To control the increasing cases, one possible approach is to raise the booster vaccination coverage from the current 64.84% to 75.15%, as suggested by simulation results. Furthermore, the spatial analysis revealed that hot spot clusters were present in the western, central, and southern regions, posing a high risk not only in densely populated areas but also in rural regions. The diffusion pattern of diphtheria clusters displayed an expansion-contagious pattern. Understanding the rising trend of diphtheria cases and their geographic distribution can offer crucial insights for government and health authorities to manage the number of diphtheria cases and make informed decisions regarding the best prevention and intervention strategies.
Collapse
Affiliation(s)
| | - Nuning Nuraini
- Department of Mathematics, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Bandung, Indonesia
- Center for Mathematical Modeling and Simulation, Institut Teknologi Bandung, Bandung, Indonesia
| | - Ade Maya Sari
- Department of Mathematics, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Bandung, Indonesia
| | - Imaniah Bazlina Wardani
- Study Program of Biology Education, Faculty of Education and Teacher Training, UIN Kiai Haji Achmad Siddiq Jember, Jember, Indonesia
| | | | | | - Bony Wiem Lestari
- Department of Public Health, Faculty of Medicine, Universitas Padjadjaran, Bandung, Indonesia
- Department of Internal Medicine, Radboud Institute for Health Sciences, Radboud University Medical Centre, Nijmegen, the Netherlands
| |
Collapse
|
12
|
Bagkur C, Amilo D, Kaymakamzade B. A fractional-order model for nosocomial infection caused by Pseudomonas aeruginosa in Northern Cyprus. Comput Biol Med 2024; 171:108094. [PMID: 38335823 DOI: 10.1016/j.compbiomed.2024.108094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Revised: 01/09/2024] [Accepted: 01/31/2024] [Indexed: 02/12/2024]
Abstract
Pseudomonas aeruginosa, a resilient gram-negative bacterium, poses a persistent threat as a leading cause of nosocomial infections, particularly in resource-constrained regions. Despite existing treatment and control measures, the bacterium continues to challenge healthcare systems, especially in developing nations. This paper introduces a fractional-order model to elucidate the dynamic behavior of nosocomial infections caused by P. aeruginosa and to compare the efficacy of carbapenems and aminoglycosides in treatment. The model's existence and uniqueness are established, and both global and local stability are confirmed. The effective reproduction number is computed, revealing an epidemic potential with a value of 1.02 in Northern Cyprus. Utilizing real-life data from a university hospital and employing numerical simulations, our results indicate that patients exhibit higher sensitivity and lower resistance to aminoglycoside treatment compared to carbapenems. Aminoglycosides consistently outperform carbapenems across key metrics, including the reduction of susceptible population, infection numbers, treatment efficacy, total infected population, hospital occupancy, and effective reproduction number. The fractional-order approach emerges as a suitable and insightful tool for studying the transmission dynamics of the disease and assessing treatment effectiveness. This research provides a robust foundation for refining treatment strategies against P. aeruginosa infections, contributing valuable insights for healthcare practitioners and policymakers alike.
Collapse
Affiliation(s)
- Cemile Bagkur
- DESAM Research Institute Near East University, Nicosia, 99010, Cyprus.
| | - David Amilo
- Mathematics Research Center, Near East University, Nicosia, 99010, Cyprus; Department of Mathematics, Near East University, Nicosia, 99010, Cyprus.
| | - Bilgen Kaymakamzade
- Mathematics Research Center, Near East University, Nicosia, 99010, Cyprus; Department of Mathematics, Near East University, Nicosia, 99010, Cyprus.
| |
Collapse
|
13
|
Olejarz JW, Roster KIO, Kissler SM, Lipsitch M, Grad YH. Optimal environmental testing frequency for outbreak surveillance. Epidemics 2024; 46:100750. [PMID: 38394927 DOI: 10.1016/j.epidem.2024.100750] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 01/12/2024] [Accepted: 02/14/2024] [Indexed: 02/25/2024] Open
Abstract
Public health surveillance for pathogens presents an optimization problem: we require enough sampling to identify intervention-triggering shifts in pathogen epidemiology, such as new introductions or sudden increases in prevalence, but not so much that costs due to surveillance itself outweigh those from pathogen-associated illness. To determine this optimal sampling frequency, we developed a general mathematical model for the introduction of a new pathogen that, once introduced, increases in prevalence exponentially. Given the relative cost of infection vs. sampling, we derived equations for the expected combined cost per unit time of disease burden and surveillance for a specified sampling frequency, and thus the sampling frequency for which the expected total cost per unit time is lowest.
Collapse
Affiliation(s)
- Jason W Olejarz
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA.
| | - Kirstin I Oliveira Roster
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Stephen M Kissler
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Department of Computer Science, University of Colorado Boulder, Boulder, CO 80309, USA
| | - Marc Lipsitch
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| | - Yonatan H Grad
- Department of Immunology and Infectious Diseases, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA; Center for Communicable Disease Dynamics, Harvard T. H. Chan School of Public Health, Boston, MA 02115, USA
| |
Collapse
|
14
|
Omeje CO, Salau AO, Eya CU. Dynamics analysis of permanent magnet synchronous motor speed control with enhanced state feedback controller using a linear quadratic regulator. Heliyon 2024; 10:e26018. [PMID: 38379983 PMCID: PMC10877355 DOI: 10.1016/j.heliyon.2024.e26018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 02/06/2024] [Accepted: 02/06/2024] [Indexed: 02/22/2024] Open
Abstract
This paper presents a robust speed regulation and control of a permanent magnet synchronous motor (PMSM). A linear quadratic regulator (LQR) based state feedback controller was developed to achieve a successful suppression of periodic disturbance of speed and torque. Sliding mode observer in conjunction with the disturbance observer was deployed in the control of motor speed. Simulations were carried out based on two compared controllers such as the state feedback controller and the conventional proportional-integral-derivative (PID) controller to attenuate the noisy effects of the external disturbance. A comparative analysis of results showed that a robust as well as an improved speed and torque dynamic performance was achieved with the state feedback (SFC) controller. A reduced periodic disturbance with percentage steady state error values of 24.17% and 23.51% was obtained with the SFC controller as compared to 38.0% and 38.37% obtained using a PID controller. The Eigen values obtained from the derived state feedback matrix (K) based on Ackerman's rule proved that the entire system operation is controllable and the performance index is marginally stable. All simulations were performed using MATLAB/SIMULINK version 2021.
Collapse
Affiliation(s)
- Crescent Onyebuchi Omeje
- Department of Electrical/Electronic Engineering, University of Port Harcourt, Rivers State, Nigeria
| | - Ayodeji Olalekan Salau
- Department of Electrical//Electronics and Computer Engineering, Afe Babalola University, Ado-Ekiti, Nigeria
- Saveetha School of Engineering, Saveetha Institute of Medical and Technical Sciences, Chennai, Tamil Nadu, India
| | - Candidus Ugwoke Eya
- Department of Electrical Engineering, University of Nigeria, Nsukka Enugu State, Nigeria
| |
Collapse
|
15
|
Una R, Glimm T. A Cellular Potts Model of the interplay of synchronization and aggregation. PeerJ 2024; 12:e16974. [PMID: 38435996 PMCID: PMC10909357 DOI: 10.7717/peerj.16974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 01/29/2024] [Indexed: 03/05/2024] Open
Abstract
We investigate the behavior of systems of cells with intracellular molecular oscillators ("clocks") where cell-cell adhesion is mediated by differences in clock phase between neighbors. This is motivated by phenomena in developmental biology and in aggregative multicellularity of unicellular organisms. In such systems, aggregation co-occurs with clock synchronization. To account for the effects of spatially extended cells, we use the Cellular Potts Model (CPM), a lattice agent-based model. We find four distinct possible phases: global synchronization, local synchronization, incoherence, and anti-synchronization (checkerboard patterns). We characterize these phases via order parameters. In the case of global synchrony, the speed of synchronization depends on the adhesive effects of the clocks. Synchronization happens fastest when cells in opposite phases adhere the strongest ("opposites attract"). When cells of the same clock phase adhere the strongest ("like attracts like"), synchronization is slower. Surprisingly, the slowest synchronization happens in the diffusive mixing case, where cell-cell adhesion is independent of clock phase. We briefly discuss potential applications of the model, such as pattern formation in the auditory sensory epithelium.
Collapse
Affiliation(s)
- Rose Una
- Department of Mathematics, Western Washington University, Bellingham, WA, United States of America
| | - Tilmann Glimm
- Department of Mathematics, Western Washington University, Bellingham, WA, United States of America
| |
Collapse
|
16
|
Brhane KW, Ahmad AG, Hina H, Emadifar H. Mathematical modeling of cholera dynamics with intrinsic growth considering constant interventions. Sci Rep 2024; 14:4616. [PMID: 38409239 PMCID: PMC10897316 DOI: 10.1038/s41598-024-55240-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Accepted: 02/21/2024] [Indexed: 02/28/2024] Open
Abstract
A mathematical model that describes the dynamics of bacterium vibrio cholera within a fixed population considering intrinsic bacteria growth, therapeutic treatment, sanitation and vaccination rates is developed. The developed mathematical model is validated against real cholera data. A sensitivity analysis of some of the model parameters is also conducted. The intervention rates are found to be very important parameters in reducing the values of the basic reproduction number. The existence and stability of equilibrium solutions to the mathematical model are also carried out using analytical methods. The effect of some model parameters on the stability of equilibrium solutions, number of infected individuals, number of susceptible individuals and bacteria density is rigorously analyzed. One very important finding of this research work is that keeping the vaccination rate fixed and varying the treatment and sanitation rates provide a rapid decline of infection. The fourth order Runge-Kutta numerical scheme is implemented in MATLAB to generate the numerical solutions.
Collapse
Affiliation(s)
| | - Abdulaziz Garba Ahmad
- Department of Applied Mathematics, Federal University of Technology, Babura, Jigawa State, Nigeria
| | - Hina Hina
- Department of Mathematics and Statistics, Women University Swabi, Swabi, KP, Pakistan
| | - Homan Emadifar
- Department of Mathematics, Saveetha School of Engineering, Saveetha Institute of Medical and Technical Sciences, Saveetha University, Chennai, 602 105, Tamil Nadu, India.
- MEU Research Unit, Middle East University, Amman, Jordan.
- Department of Mathematics, Hamedan Branch, Islamic Azad University of Hamedan, Hamadan, Iran.
| |
Collapse
|
17
|
Ardakani MH, Rezadoost H, Norouzi HR. Sequential purification of cannabidiol by two-dimensional liquid chromatography combined with modeling and simulation of elution profiles. J Chromatogr A 2024; 1717:464702. [PMID: 38310701 DOI: 10.1016/j.chroma.2024.464702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Revised: 12/15/2023] [Accepted: 01/30/2024] [Indexed: 02/06/2024]
Abstract
Cannabidiol (CBD) has garnered significant attention for its neuroprotective properties, and research on its therapeutic effects has increased dramatically in recent years. However, the systematic purification of CBD through scalable processes has remained bottleneck due to the structural similarities of the cannabinoids. Although preparative chromatography is considered as a potential solution, it is usually time-consuming and expensive. Therefore, the development of scalable strategy via fast and accurate optimization approach is crucial. The present study aimed to develop a sequential process for the scalable purification of CBD through an eco-friendly ethanolic extraction using ultrasonic assisted extraction, decarboxylation of cannabidiolic acid optimized by response surface methodology, followed by the development of off-line two-dimensional semi-preparative chromatography, boosted with stacked injection overloading. In the first dimension, a column packed with macroporous resin allows to enrich the target substance and then, the behavior of resin column for scale-up procedure were predicted and optimized by developed mathematical model. A C18 column was used in the second dimension. The CBD purity and recovery obtained were 94.3 and 82.1 %, respectively. A robust and reliable method was employed for CBD enrichment/purification, which can be generalized to other bioactive compounds in complex matrices.
Collapse
Affiliation(s)
- Mohammad Hooshyari Ardakani
- Department of Phytochemistry, Medicinal Plants and Drugs Research Institute, Shahid Beheshti University, G.C., Evin, Tehran, Iran
| | - Hassan Rezadoost
- Department of Phytochemistry, Medicinal Plants and Drugs Research Institute, Shahid Beheshti University, G.C., Evin, Tehran, Iran.
| | - Hamid Reza Norouzi
- Center of Engineering and Multiscale Modeling of Fluid Flow (CEMF), Department of Chemical Engineering, Amirkabir University of Technology (Tehran Poly Technique), Tehran, Iran
| |
Collapse
|
18
|
Go N, Arsène S, Faddeenkov I, Galland T, Martis B S, Lefaudeux D, Wang Y, Etheve L, Jacob E, Monteiro C, Bosley J, Sansone C, Pasquali C, Lehr L, Kulesza A. A quantitative systems pharmacology workflow towards optimal design and biomarker stratification of atopic dermatitis clinical trials. J Allergy Clin Immunol 2024:S0091-6749(24)00158-1. [PMID: 38369029 DOI: 10.1016/j.jaci.2023.12.031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 11/03/2023] [Accepted: 12/22/2023] [Indexed: 02/20/2024]
Abstract
BACKGROUND The development of atopic dermatitis (AD) drugs is confronted by many disease phenotypes and trial design options, which are hard to explore experimentally. OBJECTIVE Optimize AD trial design using simulations. METHODS We constructed a quantitative systems pharmacology (QSP) model of AD and standard of care (SoC) treatments and generated a phenotypically diverse virtual population whose parameter distribution is a) derived from known relationships between AD biomarkers and disease severity and b) calibrated using disease severity evolution under SoC regimens. RESULTS We applied this workflow to the immunomodulator OM-85, currently being investigated for its potential use in AD, and calibrated investigational treatment model with the efficacy profile of an existing trial (thereby enriching it with plausible marker levels and dynamics). We assessed the sensitivity of trial outcomes to trial protocol and found that for this particular example, a) the choice of endpoint is more important than the choice of dosing-regimen and b) patient selection by model-based responder enrichment could increase the expected effect size. A global sensitivity analysis revealed that only a limited subset of baseline biomarkers is needed to predict the drug response of the full virtual population. CONCLUSION This AD QSP workflow built around knowledge of marker-severity relationships as well as SoC efficacy can be tailored to specific development cases so as to optimize several trial protocol parameters and biomarker-stratification and therefore holds promise to become a powerful model-informed AD drug development and personalized medicine tool.
Collapse
|
19
|
Tada S, Yamazaki Y, Yamamoto K, Fujii K, Yamada TG, Hiroi NF, Kimura A, Funahashi A. Switching from weak to strong cortical attachment of microtubules accounts for the transition from nuclear centration to spindle elongation in metazoans. Heliyon 2024; 10:e25494. [PMID: 38356608 PMCID: PMC10865266 DOI: 10.1016/j.heliyon.2024.e25494] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2023] [Revised: 01/06/2024] [Accepted: 01/29/2024] [Indexed: 02/16/2024] Open
Abstract
The centrosome is a major microtubule organizing center in animal cells. The position of the centrosomes inside the cell is important for cell functions such as cell cycle, and thus should be tightly regulated. Theoretical models based on the forces generated along the microtubules have been proposed to account for the dynamic movements of the centrosomes during the cell cycle. These models, however, often adopted inconsistent assumptions to explain distinct but successive movements, thus preventing a unified model for centrosome positioning. For the centration of the centrosomes, weak attachment of the astral microtubules to the cell cortex was assumed. In contrast, for the separation of the centrosomes during spindle elongation, strong attachment was assumed. Here, we mathematically analyzed these processes at steady state and found that the different assumptions are proper for each process. We experimentally validated our conclusion using nematode and sea urchin embryos by manipulating their shapes. Our results suggest the existence of a molecular mechanism that converts the cortical attachment from weak to strong during the transition from centrosome centration to spindle elongation.
Collapse
Affiliation(s)
- Shohei Tada
- Center for Biosciences and Informatics, Graduate School of Fundamental Science and Technology, Keio University, Yokohama, Kanagawa, 223-8522, Japan
| | - Yoshitaka Yamazaki
- Center for Biosciences and Informatics, Graduate School of Fundamental Science and Technology, Keio University, Yokohama, Kanagawa, 223-8522, Japan
| | - Kazunori Yamamoto
- Cell Architecture Laboratory, Department of Chromosome Science, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
- Genetics Program, The Graduate University for Advanced Studies, SOKENDAI, Mishima, Shizuoka 411-8540, Japan
- Faculty of Applied Bioscience, Kanagawa Institute of Technology, Atsugi, Kanagawa, 243-0292, Japan
- Division of Developmental Physiology, Institute for Genetic Medicine, Hokkaido University, Sapporo, Hokkaido, 060-0815, Japan
| | - Ken Fujii
- Cell Architecture Laboratory, Department of Chromosome Science, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
- Genetics Program, The Graduate University for Advanced Studies, SOKENDAI, Mishima, Shizuoka 411-8540, Japan
| | - Takahiro G. Yamada
- Center for Biosciences and Informatics, Graduate School of Fundamental Science and Technology, Keio University, Yokohama, Kanagawa, 223-8522, Japan
- Department of Biosciences and Informatics, Keio University, Yokohama, Kanagawa, 223-8522, Japan
| | - Noriko F. Hiroi
- School of Medicine, Keio University, Shinjuku-ward, Tokyo, 160-8582, Japan
- Faculty of Creative Engineering, Kanagawa Institute of Technology, Atsugi, Kanagawa, 243-0292, Japan
| | - Akatsuki Kimura
- Cell Architecture Laboratory, Department of Chromosome Science, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
- Genetics Program, The Graduate University for Advanced Studies, SOKENDAI, Mishima, Shizuoka 411-8540, Japan
- Center for Data Assimilation Research and Applications, Joint Support-Center for Data Science Research, Research Organization of Information and Systems (ROIS), Tachikawa, 190-8562, Japan
| | - Akira Funahashi
- Center for Biosciences and Informatics, Graduate School of Fundamental Science and Technology, Keio University, Yokohama, Kanagawa, 223-8522, Japan
- Department of Biosciences and Informatics, Keio University, Yokohama, Kanagawa, 223-8522, Japan
| |
Collapse
|
20
|
Bhat R, Khanday MA. Traveling wave solution and the stability of critical points of an enzyme-inhibitor system under diffusion effects: with special reference to dimer molecule. Comput Methods Biomech Biomed Engin 2024:1-7. [PMID: 38321971 DOI: 10.1080/10255842.2024.2311321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Accepted: 01/21/2024] [Indexed: 02/08/2024]
Abstract
Enzymes are absolutely essential biological catalysts in human body that catalyze all cellular processes in physiological network. However, there are certain low molecular weight chemical compounds known as inhibitors, that reduce or completely inhibit the enzyme catalytic activity. Mathematical modeling plays a key role in the control and stability of metabolic enzyme inhibition. Enzyme stability is an important issue for protein engineers, because of its great importance and impact on optimal utility of material in biological tissues. In this outlook, we have first determined the existence of traveling wave solution for the enzyme-inhibitor system and then emphasized the stability of critical points that arise in the reactions. The study of traveling wave solution of an enzyme-inhibitor system with reaction diffusion equations involve quite complex mathematical analysis. The results obtained in this model indicate that the traveling wave solution may give a well explained method for improving enzyme kinetic stability. The present study will be helpful in understanding the stability of critical points of an enzyme-inhibitor system to give an idea about the inhibition of less stable enzymes. Moreover, the role of diffusion on the enzyme activity has been exhaustively discussed using mathematical tools related to eigen values and eigen function analysis.
Collapse
Affiliation(s)
- Roohi Bhat
- Department of Mathematics, University of Kashmir, Srinagar, India
- Department of Mathematics, School of Chemical Engineering and Physical Sciences, Lovely Professional University, Phagwara, Punjab
| | - M A Khanday
- Department of Mathematics, University of Kashmir, Srinagar, India
| |
Collapse
|
21
|
Jannat MAH, Park SH, Hwang S. Modeling interactions of Clostridium cadaveris and Clostridium sporogenes in anaerobic acidogenesis of glucose and peptone. Bioresour Technol 2024; 393:130099. [PMID: 38013037 DOI: 10.1016/j.biortech.2023.130099] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 11/22/2023] [Accepted: 11/22/2023] [Indexed: 11/29/2023]
Abstract
This study focuses on developing a mathematical model to assess interaction among acidogenic bacteria during the anaerobic degradation of two substrates. Clostridium cadaveris and Clostridium sporogenes were cultured in various combinations with glucose and peptone. Parameter estimates are given for both conventional Monod parameters from single substrate-single species cultures and sum kinetics with interaction parameters obtained from dual substrate-single species cultures. The presence of multiple substrates led to both inhibitory and enhancing effects on biodegradation rates for dual substrates compared to single substrate cultures. A new model of interspecies interaction was developed within the framework of Lotka-Volterra incorporating substrate interaction parameters, with a focus on accuracy, realism, simplicity, and biological significance. The model demonstrated competitive interaction for resource sharing and the additional non-linearity parameter eliminated the constraint of the linear relationship between growth rate and population density.
Collapse
Affiliation(s)
- Md Abu Hanifa Jannat
- Division of Environmental Science and Engineering, Pohang University of Science and Technology, 77 Cheongam-Ro, Nam-Gu, Pohang, Gyeongbuk, South Korea.
| | - Sang Hyeok Park
- Division of Environmental Science and Engineering, Pohang University of Science and Technology, 77 Cheongam-Ro, Nam-Gu, Pohang, Gyeongbuk, South Korea.
| | - Seokhwan Hwang
- Division of Environmental Science and Engineering, Pohang University of Science and Technology, 77 Cheongam-Ro, Nam-Gu, Pohang, Gyeongbuk, South Korea; Yonsei University Institute for Convergence Research and Education in Advanced Technology (I-CREATE), 85, Songdogwahak-ro, Yeonsu-gu, Incheon 21983, Republic of Korea.
| |
Collapse
|
22
|
Jayathungage Don TD, Safaei S, Maso Talou GD, Russell PS, Phillips ARJ, Reynolds HM. Computational fluid dynamic modeling of the lymphatic system: a review of existing models and future directions. Biomech Model Mechanobiol 2024; 23:3-22. [PMID: 37902894 PMCID: PMC10901951 DOI: 10.1007/s10237-023-01780-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 10/02/2023] [Indexed: 11/01/2023]
Abstract
Historically, research into the lymphatic system has been overlooked due to both a lack of knowledge and limited recognition of its importance. In the last decade however, lymphatic research has gained substantial momentum and has included the development of a variety of computational models to aid understanding of this complex system. This article reviews existing computational fluid dynamic models of the lymphatics covering each structural component including the initial lymphatics, pre-collecting and collecting vessels, and lymph nodes. This is followed by a summary of limitations and gaps in existing computational models and reasons that development in this field has been hindered to date. Over the next decade, efforts to further characterize lymphatic anatomy and physiology are anticipated to provide key data to further inform and validate lymphatic fluid dynamic models. Development of more comprehensive multiscale- and multi-physics computational models has the potential to significantly enhance the understanding of lymphatic function in both health and disease.
Collapse
Affiliation(s)
| | - Soroush Safaei
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Gonzalo D Maso Talou
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Peter S Russell
- School of Biological Sciences, The University of Auckland, Auckland, New Zealand
- Surgical and Translational Research Centre, Department of Surgery, Faculty of Medical and Health Sciences, The University of Auckland, Auckland, New Zealand
| | - Anthony R J Phillips
- School of Biological Sciences, The University of Auckland, Auckland, New Zealand
- Surgical and Translational Research Centre, Department of Surgery, Faculty of Medical and Health Sciences, The University of Auckland, Auckland, New Zealand
| | - Hayley M Reynolds
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.
| |
Collapse
|
23
|
Simões A, Veiga F, Vitorino C. Question-based review for pharmaceutical development: An enhanced quality approach. Eur J Pharm Biopharm 2024; 195:114174. [PMID: 38160986 DOI: 10.1016/j.ejpb.2023.114174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Revised: 12/21/2023] [Accepted: 12/26/2023] [Indexed: 01/03/2024]
Abstract
Over the last years, the pharmaceutical industry has faced real challenges regarding quality assurance. In this context, the establishment of more holistic approaches to the pharmaceutical development has been encouraged. The emergence of the Quality by Design (QbD) paradigm as systematic, scientific and risk-based methodology introduced a new concept of pharmaceutical quality. In essence, QbD can be interpreted as a strategy to maximize time and cost savings. An in-depth understanding of the formulation and manufacturing process is demanded to optimize the safety, efficacy and quality of a drug product at all stages of development. This innovative approach streamlines the pharmaceutical Research and Development (R&D) process, provides greater manufacturing flexibility and reduces regulatory burden. To assist in QbD implementation, International Conference on Harmonisation (ICH), U.S. Food and Drug Administration (FDA) and European Medicines Agency (EMA) organized and launched QbD principles in their guidance for industry, identifying key concepts and tools to design and develop a high-quality drug product. Despite the undeniable advantages of the QbD approach, and the widespread information on QbD regulatory expectations, its full implementation in the pharmaceutical field is still limited. The present review aims to establish a crosswise overview on the current application status of QbD within the framework of the ICH guidelines (ICH Q8(R2) - Q14 and ICH Q2(R2)). Moreover, it outlines the way information gathered from the QbD methodology is being harmonized in Marketing Authorization Applications (MAAs) for European market approval. This work also highlights the challenges that hinder the deployment of the QbD strategy as a standard practice.
Collapse
Affiliation(s)
- Ana Simões
- Faculty of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal; Associated Laboratory for Green Chemistry of the Network of Chemistry and Technology (LAQV/REQUIMTE), Group of Pharmaceutical Technology, Faculty of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal
| | - Francisco Veiga
- Faculty of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal; Associated Laboratory for Green Chemistry of the Network of Chemistry and Technology (LAQV/REQUIMTE), Group of Pharmaceutical Technology, Faculty of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal
| | - Carla Vitorino
- Faculty of Pharmacy, University of Coimbra, 3000-548 Coimbra, Portugal; Coimbra Chemistry Centre, Institute of Molecular Sciences - IMS, Department of Chemistry, University of Coimbra 3004-535 Coimbra, Portugal.
| |
Collapse
|
24
|
Hatami T, Jarles Santos de Araújo E, Luiz Baião Dias A, Helena Innocentini Mei L, Martínez J. Mechanism of multicyclic β-carotene impregnation into corn starch aerogels via supercritical CO 2 with mathematical modeling. Food Res Int 2024; 178:114002. [PMID: 38309888 DOI: 10.1016/j.foodres.2024.114002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 01/03/2024] [Accepted: 01/08/2024] [Indexed: 02/05/2024]
Abstract
β-carotene, a natural dye renowned for its antioxidant and provitamin A activities, is hindered from direct use in food and drug products due to its susceptibility to oxidation, easy isomerization under light, heat, or acidic conditions, as well as its low water solubility and oral bioavailability. In this study, we addressed these challenges by loading β-carotene into corn starch aerogels via supercritical carbon dioxide (sc-CO2) and assessed its loading contributions through adsorption during contact time and precipitation during depressurization. The loading process was studied under two cycles at pressure of 30 MPa, temperature of 40 °C, depressurization rate of 0.4 MPa/min, and co-solvent (ethanol) mass percentage of 1.2 %. Experiments found adsorption minimally contributed to impregnation, while precipitation became the primary loading mechanism. The subsequent work focused on a mathematical model describing β-carotene loading into corn starch aerogels via precipitation, using the law of conservation of mass and classical nucleation theory. The model shows that using pure CO2 results in a loading efficiency of 0.10 mg β-carotene/g aerogel, while with CO2 and 1.2 % ethanol as the co-solvent, the loading efficiency increases threefold to 0.30 mg β-carotene/g aerogel.
Collapse
Affiliation(s)
- Tahmasb Hatami
- Departamento de Engenharia de Materiais e Bioprocessos, Faculdade de Engenharia Química, Universidade Estadual de Campinas (UNICAMP), Av. Albert Einstein 500, 13083-852, Campinas, SP, Brasil.
| | - Erick Jarles Santos de Araújo
- Laboratório de Alta Pressão em Engenharia de Alimentos (LAPEA), Faculdade de Engenharia de Alimentos (FEA), Universidade Estadual de Campinas (UNICAMP), 13083-862 Campinas, Brasil
| | - Arthur Luiz Baião Dias
- Laboratório de Alta Pressão em Engenharia de Alimentos (LAPEA), Faculdade de Engenharia de Alimentos (FEA), Universidade Estadual de Campinas (UNICAMP), 13083-862 Campinas, Brasil
| | - Lucia Helena Innocentini Mei
- Departamento de Engenharia de Materiais e Bioprocessos, Faculdade de Engenharia Química, Universidade Estadual de Campinas (UNICAMP), Av. Albert Einstein 500, 13083-852, Campinas, SP, Brasil
| | - Julian Martínez
- Laboratório de Alta Pressão em Engenharia de Alimentos (LAPEA), Faculdade de Engenharia de Alimentos (FEA), Universidade Estadual de Campinas (UNICAMP), 13083-862 Campinas, Brasil
| |
Collapse
|
25
|
Almutairi D. Mathematical modelling and heat transfer observations for Jeffrey nanofluid with applications of extended Fourier theory and temperature dependent thermal conductivity. Heliyon 2024; 10:e24353. [PMID: 38298636 PMCID: PMC10828649 DOI: 10.1016/j.heliyon.2024.e24353] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Revised: 11/01/2023] [Accepted: 01/08/2024] [Indexed: 02/02/2024] Open
Abstract
The suspension of non-Newtonian materials with nanoparticles is important to enhance the thermal phenomenon in various engineering and industrial processes. The versatile research in nanomaterials provide different applications in thermal processes, heat exchangers, thermoelectric devices, HVAC systems, energy processes etc. Following to such novel motivations in mind, current research endorsed the enhancement in heat transfer due to suspension of Jeffrey nanofluid comprising the variable thermal conductivity. The cause of flow is associated to two disks attaining fixed distance. The modified developed relations for Fourier's hypothesis are utilized to model the problem. The flow problem is modeled with appliance of fundamental novel laws. By applying suitable transformations, corresponding differential equations are renovated into dimensionless forms which are solved with applications of analytic homotopic algorithm. The behavior of temperature and velocity due to various parameters is discussed. The numerical calculations have been done for wall shear force and Nusselt number. The results show that the velocity profile boosted due to variation of stretching ratio constant. The enhancement in heat transfer is observed due to Reynolds number. Moreover, the increasing observations for wall shear force in upper and lower disk surfaces are obtained against larger material parameter. The simulated results may find applications in improving heat transfer phenomenon, manufacturing systems, recovery processes, cooling systems, chemical phenomenon, fuel cells etc.
Collapse
Affiliation(s)
- D.K. Almutairi
- Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Al-Majmaah 11952, Saudi Arabia
| |
Collapse
|
26
|
Arabameri A, Arab S. Understanding the Interplay of CAR-NK Cells and Triple-Negative Breast Cancer: Insights from Computational Modeling. Bull Math Biol 2024; 86:20. [PMID: 38240892 DOI: 10.1007/s11538-023-01247-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2023] [Accepted: 12/14/2023] [Indexed: 01/23/2024]
Abstract
Chimeric antigen receptor (CAR)-engineered natural killer (NK) cells have recently emerged as a promising and safe alternative to CAR-T cells for targeting solid tumors. In the case of triple-negative breast cancer (TNBC), traditional cancer treatments and common immunotherapies have shown limited effectiveness. However, CAR-NK cells have been successfully employed to target epidermal growth factor receptor (EGFR) on TNBC cells, thereby enhancing the efficacy of immunotherapy. The effectiveness of CAR-NK-based immunotherapy is influenced by various factors, including the vaccination dose, vaccination pattern, and tumor immunosuppressive factors in the microenvironment. To gain insights into the dynamics and effects of CAR-NK-based immunotherapy, we propose a computational model based on experimental data and immunological theories. This model integrates an individual-based model that describes the interplay between the tumor and the immune system, along with an ordinary differential equation model that captures the variation of inflammatory cytokines. Computational results obtained from the proposed model shed light on the conditions necessary for initiating an effective anti-tumor response. Furthermore, global sensitivity analysis highlights the issue of low persistence of CAR-NK cells in vivo, which poses a significant challenge for the successful clinical application of these cells. Leveraging the model, we identify the optimal vaccination time, vaccination dose, and time interval between injections for maximizing therapeutic outcomes.
Collapse
Affiliation(s)
- Abazar Arabameri
- Department of Electrical Engineering, University of Zanjan, Zanjan, Iran.
| | - Samaneh Arab
- Department of Tissue Engineering and Applied Cell Sciences, School of Medicine, Semnan University of Medical Sciences, Semnan, Iran
| |
Collapse
|
27
|
Meyer L, Ribot M, Yvinec R. A Lifshitz-Slyozov type model for adipocyte size dynamics: limit from Becker-Döring system and numerical simulation. J Math Biol 2024; 88:16. [PMID: 38231273 DOI: 10.1007/s00285-023-02036-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 12/04/2023] [Accepted: 12/10/2023] [Indexed: 01/18/2024]
Abstract
Biological data show that the size distribution of adipose cells follows a bimodal distribution. In this work, we introduce a Lifshitz-Slyozov type model, based on a transport partial differential equation, for the dynamics of the size distribution of adipose cells. We prove a new convergence result from the related Becker-Döring model, a system composed of several ordinary differential equations, toward mild solutions of the Lifshitz-Slyozov model using distribution tail techniques. Then, this result allows us to propose a new advective-diffusive model, the second-order diffusive Lifshitz-Slyozov model, which is expected to better fit the experimental data. Numerical simulations of the solutions to the diffusive Lifshitz-Slyozov model are performed using a well-balanced scheme and compared to solutions to the transport model. Those simulations show that both bimodal and unimodal profiles can be reached asymptotically depending on several parameters. We put in evidence that the asymptotic profile for the second-order system does not depend on initial conditions, unlike for the transport Lifshitz-Slyozov model.
Collapse
Affiliation(s)
- Léo Meyer
- Institut Denis Poisson, Université d'Orléans, Orléans, France.
| | - Magali Ribot
- Institut Denis Poisson, Université d'Orléans, Orléans, France
| | - Romain Yvinec
- Inria, Centre Inria de Saclay, Université Paris-Saclay, Palaiseau, France
- PRC, INRAE, CNRS, Université de Tours, Nouzilly, France
| |
Collapse
|
28
|
Zhao H, Zhan J, Yun H, Mu H, Zhang D, An L, Yao H, Zhang G. Deciphering the intricate dielectric relaxation processes of cellulose paper: Extraction of distribution of relaxation time and analysis of degradation characteristics. Carbohydr Polym 2024; 324:121497. [PMID: 37985048 DOI: 10.1016/j.carbpol.2023.121497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 10/11/2023] [Accepted: 10/13/2023] [Indexed: 11/22/2023]
Abstract
Cellulose material is a dielectric with intricate microscopic relaxation processes due to its complex structure. However, conventional models and curve fitting methods used for tracing and analyzing these processes often fail to capture crucial dielectric information. This paper aimed to extract the Distribution of Relaxation Time (DRT), the most fundamental and effective dielectric information providing the time scale and relative contribution of all microscopic relaxation processes. First, a distributed extended Debye model with infinite branches was constructed based on the microscopic nature of dielectric relaxation. Then, an implicit equation of the DRT function was established, inspired by the mathematical principles of infinite subdivision and summation. To obtain the numeral solution of the DRT function, a regularization method was proposed and validated. Finally, the approach was applied to cellulose insulating paper with varying degradation degrees. The relaxation process with a long time constant played a significant role, and variations during the degradation process were attributed to reduced activation energy. With clear physical interpretation and robust mathematical foundation, our method sheds light on the intricate dielectric relaxation processes in cellulose. This not only enhances the theoretical understanding and practical application of cellulose materials but also provides valuable insights for the analysis and application of other materials.
Collapse
Affiliation(s)
- Haoxiang Zhao
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an, China; Department of Mathematics and Computer Science, Eindhoven University of Technology, Eindhoven, the Netherlands.
| | - Jiangyang Zhan
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an, China.
| | - Hao Yun
- China Nuclear Power Operation Technology Corporation., Ltd, Wuhan, China.
| | - Haibao Mu
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an, China.
| | - Daning Zhang
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an, China.
| | - Lixuan An
- Department of Data Analysis and Mathematical Modelling, Ghent University, Ghent, Belgium.
| | - Huanmin Yao
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an, China
| | - Guanjun Zhang
- State Key Laboratory of Electrical Insulation and Power Equipment, School of Electrical Engineering, Xi'an Jiaotong University, Xi'an, China.
| |
Collapse
|
29
|
Bai J, Wang J. Modeling long COVID dynamics: Impact of underlying health conditions. J Theor Biol 2024; 576:111669. [PMID: 37977479 PMCID: PMC10754059 DOI: 10.1016/j.jtbi.2023.111669] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 10/03/2023] [Accepted: 11/13/2023] [Indexed: 11/19/2023]
Abstract
We propose a new mathematical model to investigate the population dynamics of long COVID, with a focus on the impact of chronic health conditions. Our model connects long COVID with the transmission of COVID-19 so as to accurately predict the prevalence of long COVID from the progression of the infection in the host population. The model additionally incorporates the effects of COVID-19 vaccination. We implement the model with data from both the US and the UK to demonstrate the real-world applications of this modeling framework.
Collapse
Affiliation(s)
- Jie Bai
- School of Mathematics and Statistics, Liaoning University, Shenyang 110036, China.
| | - Jin Wang
- Department of Mathematics, University of Tennessee at Chattanooga, Chattanooga TN 37403, USA.
| |
Collapse
|
30
|
Mohsin N, Enderling H, Brady-Nicholls R, Zahid MU. Simulating tumor volume dynamics in response to radiotherapy: Implications of model selection. J Theor Biol 2024; 576:111656. [PMID: 37952611 DOI: 10.1016/j.jtbi.2023.111656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2022] [Revised: 10/27/2023] [Accepted: 11/05/2023] [Indexed: 11/14/2023]
Abstract
From the beginning of the usage of radiotherapy (RT) for cancer treatment, mathematical modeling has been integral to understanding radiobiology and for designing treatment approaches and schedules. There has been extensive modeling of response to RT with the inclusion of various degrees of biological complexity. In this study, we compare three models of tumor volume dynamics: (1) exponential growth with RT directly reducing tumor volume, (2) logistic growth with direct tumor volume reduction, and (3) logistic growth with RT reducing the tumor carrying capacity with the objective of understanding the implications of model selection and informing the process of model calibration and parameterization. For all three models, we: examined the rates of change in tumor volume during and RT treatment course; performed parameter sensitivity and identifiability analyses; and investigated the impact of the parameter sensitivity on the tumor volume trajectories. In examining the tumor volume dynamics trends, we coined a new metric - the point of maximum reduction of tumor volume (MRV) - to quantify the magnitude and timing of the expected largest impact of RT during a treatment course. We found distinct timing differences in MRV, dependent on model selection. The parameter identifiability and sensitivity analyses revealed the interdependence of the different model parameters and that it is only possible to independently identify tumor growth and radiation response parameters if the underlying tumor growth rate is sufficiently large. Ultimately, the results of these analyses help us to better understand the implications of model selection while simultaneously generating falsifiable hypotheses about MRV timing that can be tested on longitudinal measurements of tumor volume from pre-clinical or clinical data with high acquisition frequency. Although, our study only compares three particular models, the results demonstrate that caution is necessary in selecting models of response to RT, given the artifacts imposed by each model.
Collapse
Affiliation(s)
- Nuverah Mohsin
- Dr. Kiran C. Patel College of Allopathic Medicine, Nova Southeastern University, Fort Lauderdale, FL, United States
| | - Heiko Enderling
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA; Institute for Data Science in Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - Renee Brady-Nicholls
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, United States.
| | - Mohammad U Zahid
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, United States.
| |
Collapse
|
31
|
Montévil M. Mathematical Modeling in the Study of Organisms and Their Parts. Methods Mol Biol 2024; 2745:105-119. [PMID: 38060182 DOI: 10.1007/978-1-0716-3577-3_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2023]
Abstract
Mathematical modeling is a very powerful tool to understand natural phenomena. Such a tool carries its own assumptions and should always be used critically. In this chapter we highlight the key ingredients and steps of modeling and focus on their biological interpretation. Particularly, we discuss the role of theoretical principles in writing models. We also highlight the meaning and interpretation of equations. The main aim of this chapter is to facilitate the interaction between biologists and mathematical modelers. We focus on the case of cell proliferation and motility in the context of multicellular organisms.
Collapse
Affiliation(s)
- Maël Montévil
- Centre Cavaillès, République des savoirs UAR 3608, ÉNS-PSL and CNRS, Paris, France
| |
Collapse
|
32
|
Gouw ZAR, Jeong J, Rimner A, Lee NY, Jackson A, Fu A, Sonke JJ, Deasy JO. "Primer shot" fractionation with an early treatment break is theoretically superior to consecutive weekday fractionation schemes for early-stage non-small cell lung cancer. Radiother Oncol 2024; 190:110006. [PMID: 37972733 DOI: 10.1016/j.radonc.2023.110006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 10/14/2023] [Accepted: 11/02/2023] [Indexed: 11/19/2023]
Abstract
PURPOSE Radiotherapy is traditionally given in equally spaced weekday fractions. We hypothesize that heterogeneous interfraction intervals can increase radiosensitivity via reoxygenation. Through modeling, we investigate whether this minimizes local failures and toxicity for early-stage non-small cell lung cancer (NSCLC). METHODS Previously, a tumor dose-response model based on resource competition and cell-cycle-dependent radiosensitivity accurately predicted local failure rates for early-stage NSCLC cohorts. Here, the model mathematically determined non-uniform inter-fraction intervals minimizing local failures at similar normal tissue toxicity risk, i.e., iso-BED3 (iso-NTCP) for fractionation schemes 18Gyx3, 12Gyx4, 10Gyx5, 7.5Gyx8, 5Gyx12, 4Gyx15. Next, we used these optimized schedules to reduce toxicity risk (BED3) while maintaining stable local failures (TCP). RESULTS Optimal schedules consistently favored a "primer shot" fraction followed by a 2-week break, allowing tumor reoxygenation. Increasing or decreasing the assumed baseline hypoxia extended or shortened this optimal break by up to one week. Fraction sizes of 7.5 Gy and up required a single primer shot, while smaller fractions needed one or two extra fractions for full reoxygenation. The optimized schedules, versus consecutive weekday fractionation, predicted absolute LF reductions of 4.6%-7.4%, except for the already optimal LF rate seen for 18Gyx3. Primer shot schedules could also reduce BED3 at iso-TCP with the biggest improvements for the shortest schedules (94.6Gy reduction for 18Gyx3). CONCLUSION A validated simulation model clearly supports non-standard "primer shot" fractionation, reducing the impact of hypoxia-induced radioresistance. A limitation of this study is that primer-shot fractionation is outside prior clinical experience and therefore will require clinical studies for definitive testing.
Collapse
Affiliation(s)
- Z A R Gouw
- Memorial Sloan Kettering Cancer Center, Department of Medical Physics, New York, NY, USA; The Netherlands Cancer Institute, Amsterdam, Department of Radiation Oncology, the Netherlands.
| | - J Jeong
- Memorial Sloan Kettering Cancer Center, Department of Medical Physics, New York, NY, USA
| | - A Rimner
- Memorial Sloan Kettering Cancer Center, Department of Radiation Oncology, New York, NY, USA
| | - N Y Lee
- Memorial Sloan Kettering Cancer Center, Department of Radiation Oncology, New York, NY, USA
| | - A Jackson
- Memorial Sloan Kettering Cancer Center, Department of Medical Physics, New York, NY, USA
| | - A Fu
- Memorial Sloan Kettering Cancer Center, Department of Medical Physics, New York, NY, USA
| | - J-J Sonke
- The Netherlands Cancer Institute, Amsterdam, Department of Radiation Oncology, the Netherlands
| | - J O Deasy
- Memorial Sloan Kettering Cancer Center, Department of Medical Physics, New York, NY, USA
| |
Collapse
|
33
|
Yonekura T, Sugiyama M. A new mathematical model of phyllotaxis to solve the genuine puzzle spiromonostichy. J Plant Res 2024; 137:143-155. [PMID: 37833503 PMCID: PMC10764405 DOI: 10.1007/s10265-023-01503-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Accepted: 10/06/2023] [Indexed: 10/15/2023]
Abstract
Arrangement of plant leaves around the stem, termed phyllotaxis, exhibits beautiful and mysterious regularities and has been one of the most attractive subjects of biological pattern formation. After the long history of studies on phyllotaxis, it is now widely accepted that the inhibitory effect of existing leaf primordia on new primordium formation determines phyllotactic patterning. However, costoid phyllotaxis unique to Costaceae of Zingiberales, displaying spiromonostichy characterized by a steep spiral with a small divergence angle, seems to disagree with the inhibitory effect-based mechanism and has remained as a "genuine puzzle". We developed a new mathematical model, hypothesizing that each leaf primordium exerts not only the inhibitory effect but also some inductive effect. Computer simulations with the new model successfully generated a spiromonostichous pattern when these two effects met a certain relationship. The obtained spiromonostichy matched the real costoid phyllotaxis observed with Costus megalobractea, particularly for the decrease of the divergence angle associated with the enlargement of the shoot apical meristem. The new model was also shown to be able to produce a one-sided distichous pattern that is seen in phyllotaxis of a few plants of Zingiberales and has never been addressed in the previous model studies. These results implicated inductive as well as inhibitory mechanisms in phyllotactic patterning, at least in Zingiberales.
Collapse
Affiliation(s)
- Takaaki Yonekura
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo, 113-0033, Japan.
| | - Munetaka Sugiyama
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo, 113-0033, Japan.
| |
Collapse
|
34
|
Mody H, Nair S, Rump A, Vaidya TR, Garrett TJ, Lesko L, Ait-Oudhia S. Identification of Novel and Early Biomarkers for Cisplatin-induced Nephrotoxicity and the Nephroprotective Role of Cimetidine using a Pharmacometabolomic-based Approach Coupled with In Vitro Toxicodynamic Modeling and Simulation. J Pharm Sci 2024; 113:268-277. [PMID: 37992870 DOI: 10.1016/j.xphs.2023.11.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Revised: 11/16/2023] [Accepted: 11/16/2023] [Indexed: 11/24/2023]
Abstract
Cisplatin is widely used for the treatment of various types of cancer. However, cisplatin-induced nephrotoxicity (CIN) is frequently observed in patients receiving cisplatin therapy which poses a challenge in its clinical utility. Currently used clinical biomarkers for CIN are not adequate for early detection of nephrotoxicity, hence there is a need to identify potential early biomarkers in predicting CIN. In the current study, a combination of in vitro toxicodynamic (TD) modeling and untargeted global metabolomics approach was used to identify novel potential metabolite biomarkers for early detection of CIN. In addition, we investigated the protective role of cimetidine (CIM), an inhibitor of the organic cation transporter 2 (OCT2), in suppressing CIN. We first characterized the time-course of nephrotoxic effects of cisplatin (CIS) and the protective effects of CIM in a human pseudo-immortalized renal proximal tubule epithelial cell line (RPTEC), SA7K cell line. Secondly, we used a mathematical cell-level, in vitro TD modeling approach to quantitatively characterize the time-course effects of CIS and CIM as single agents and combination in SA7K cells. Based on the experimental and modeling results, we selected relevant concentrations of CIS and CIM for our metabolomics study. With the help of PCA (Principal Component Analysis) and PLS-DA (Projection to Latent Structure - Discriminate Analysis) analyses, we confirmed global metabolome changes for different groups (CIS, CIM, CIS+CIM vs control) in SA7K cells. Based on the criterion of a p-value ≤ 0.05 and a fold change ≥ 2 or ≤ 0.5, we identified 20 top metabolites that were significantly changed during the early phase i.e. within first 12 h of CIS treatment. Finally, pathway analysis was conducted that revealed the key metabolic pathways that were most impacted in CIN.
Collapse
Affiliation(s)
- Hardik Mody
- Center for Pharmacometrics and Systems Pharmacology, Department of Pharmaceutics, College of Pharmacy, University of Florida, FL, USA
| | - Sreenath Nair
- Pharmaceutical Sciences Department, St. Jude Children's Research Hospital, Memphis, TN, USA
| | - Adrian Rump
- Center for Pharmacometrics and Systems Pharmacology, Department of Pharmaceutics, College of Pharmacy, University of Florida, FL, USA
| | - Tanaya R Vaidya
- Center for Pharmacometrics and Systems Pharmacology, Department of Pharmaceutics, College of Pharmacy, University of Florida, FL, USA
| | - Timothy J Garrett
- Southeast Center for Integrated Metabolomics, University of Florida, Gainesville, FL, USA
| | - Lawrence Lesko
- Center for Pharmacometrics and Systems Pharmacology, Department of Pharmaceutics, College of Pharmacy, University of Florida, FL, USA
| | - Sihem Ait-Oudhia
- Quantitative Pharmacology and Pharmacometrics (QP2), Merck & Co., Inc, Rahway, NJ, USA.
| |
Collapse
|
35
|
Tonnis K, Jaworska J, Kasting GB. Modeling the percutaneous absorption of solvent-deposited solids over a wide dose range: II. Weak electrolytes. J Control Release 2024; 365:435-447. [PMID: 37996054 DOI: 10.1016/j.jconrel.2023.11.038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Revised: 11/13/2023] [Accepted: 11/19/2023] [Indexed: 11/25/2023]
Abstract
Dermal absorption of weak electrolytes applied to skin from pharmaceutical and cosmetic compositions is an important consideration for both their efficacy and skin safety. We developed a mechanistic, physics-based framework that simulates this process for leave on applications following solvent deposition. We incorporated this framework into our finite dose computational skin permeation model previously tested with nonelectrolytes to generate quantitative predictions for weak electrolytes. To test the model, we analyzed experimental data from an in vitro human skin permeation study of a weak acid (benzoic acid) and a weak base (propranolol) and their sodium and hydrochloride salts from simple, ethanol/water vehicles as a function of dose and ionization state. Key factors controlling absorption are the pH and buffer capacity of the dose solution, the dissolution rate of precipitated solids into a lipid boundary layer and the rate of conversion of the deposited solid to its conjugate form as the nonionized component permeates and (sometimes) evaporates from the skin surface. The resulting framework not only describes the current test data but has the potential to predict the absorption of other weak electrolytes following topical application.
Collapse
Affiliation(s)
- Kevin Tonnis
- College of Engineering and Applied Science, The University of Cincinnati, Cincinnati, OH 45221, USA
| | - Joanna Jaworska
- The Procter & Gamble Company, Data and Modeling Sciences, Brussels Innovation Center, Belgium
| | - Gerald B Kasting
- The James L. Winkle College of Pharmacy, The University of Cincinnati, Cincinnati, OH 45267-0514, USA.
| |
Collapse
|
36
|
Li Y, Lu L, Androulakis IP. The Physiological and Pharmacological Significance of the Circadian Timing of the HPA Axis: A Mathematical Modeling Approach. J Pharm Sci 2024; 113:33-46. [PMID: 37597751 PMCID: PMC10840710 DOI: 10.1016/j.xphs.2023.08.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 08/02/2023] [Accepted: 08/02/2023] [Indexed: 08/21/2023]
Abstract
As a potent endogenous regulator of homeostasis, the circadian time-keeping system synchronizes internal physiology to periodic changes in the external environment to enhance survival. Adapting endogenous rhythms to the external time is accomplished hierarchically with the central pacemaker located in the suprachiasmatic nucleus (SCN) signaling the hypothalamus-pituitary-adrenal (HPA) axis to release hormones, notably cortisol, which help maintain the body's circadian rhythm. Given the essential role of HPA-releasing hormones in regulating physiological functions, including immune response, cell cycle, and energy metabolism, their daily variation is critical for the proper function of the circadian timing system. In this review, we focus on cortisol and key fundamental properties of the HPA axis and highlight their importance in controlling circadian dynamics. We demonstrate how systems-driven, mathematical modeling of the HPA axis complements experimental findings, enhances our understanding of complex physiological systems, helps predict potential mechanisms of action, and elucidates the consequences of circadian disruption. Finally, we outline the implications of circadian regulation in the context of personalized chronotherapy. Focusing on the chrono-pharmacology of synthetic glucocorticoids, we review the challenges and opportunities associated with moving toward personalized therapies that capitalize on circadian rhythms.
Collapse
Affiliation(s)
- Yannuo Li
- Chemical & Biochemical Engineering Department, Piscataway, NJ 08854, USA
| | - Lingjun Lu
- Chemical & Biochemical Engineering Department, Piscataway, NJ 08854, USA
| | - Ioannis P Androulakis
- Chemical & Biochemical Engineering Department, Piscataway, NJ 08854, USA; Biomedical Engineering Department, Rutgers University, Piscataway, NJ 08540, USA.
| |
Collapse
|
37
|
Petrella JR, Jiang J, Sreeram K, Dalziel S, Doraiswamy PM, Hao W. Personalized Computational Causal Modeling of the Alzheimer Disease Biomarker Cascade. J Prev Alzheimers Dis 2024; 11:435-444. [PMID: 38374750 DOI: 10.14283/jpad.2023.134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
BACKGROUND Mathematical models of complex diseases, such as Alzheimer's disease, have the potential to play a significant role in personalized medicine. Specifically, models can be personalized by fitting parameters with individual data for the purpose of discovering primary underlying disease drivers, predicting natural history, and assessing the effects of theoretical interventions. Previous work in causal/mechanistic modeling of Alzheimer's Disease progression has modeled the disease at the cellular level and on a short time scale, such as minutes to hours. No previous studies have addressed mechanistic modeling on a personalized level using clinically validated biomarkers in individual subjects. OBJECTIVES This study aimed to investigate the feasibility of personalizing a causal model of Alzheimer's Disease progression using longitudinal biomarker data. DESIGN/SETTING/PARTICIPANTS/MEASUREMENTS We chose the Alzheimer Disease Biomarker Cascade model, a widely-referenced hypothetical model of Alzheimer's Disease based on the amyloid cascade hypothesis, which we had previously implemented mathematically as a mechanistic model. We used available longitudinal demographic and serial biomarker data in over 800 subjects across the cognitive spectrum from the Alzheimer's Disease Neuroimaging Initiative. The data included participants that were cognitively normal, had mild cognitive impairment, or were diagnosed with dementia (probable Alzheimer's Disease). The model consisted of a sparse system of differential equations involving four measurable biomarkers based on cerebrospinal fluid proteins, imaging, and cognitive testing data. RESULTS Personalization of the Alzheimer Disease Biomarker Cascade model with individual serial biomarker data yielded fourteen personalized parameters in each subject reflecting physiologically meaningful characteristics. These included growth rates, latency values, and carrying capacities of the various biomarkers, most of which demonstrated significant differences across clinical diagnostic groups. The model fits to training data across the entire cohort had a root mean squared error (RMSE) of 0.09 (SD 0.081) on a variable scale between zero and one, and were robust, with over 90% of subjects showing an RMSE of < 0.2. Similarly, in a subset of subjects with data on all four biomarkers in at least one test set, performance was high on the test sets, with a mean RMSE of 0.15 (SD 0.117), with 80% of subjects demonstrating an RMSE < 0.2 in the estimation of future biomarker points. Cluster analysis of parameters revealed two distinct endophenotypic groups, with distinct biomarker profiles and disease trajectories. CONCLUSION Results support the feasibility of personalizing mechanistic models based on individual biomarker trajectories and suggest that this approach may be useful for reclassifying subjects on the Alzheimer's clinical spectrum. This computational modeling approach is not limited to the Alzheimer Disease Biomarker Cascade hypothesis, and can be applied to any mechanistic hypothesis of disease progression in the Alzheimer's field that can be monitored with biomarkers. Thus, it offers a computational platform to compare and validate various disease hypotheses, personalize individual biomarker trajectories and predict individual response to theoretical prevention and therapeutic intervention strategies.
Collapse
Affiliation(s)
- J R Petrella
- Jeffrey R. Petrella, Department of Radiology, Duke University School of Medicine, DUMC - Box 3808 , 27710-3808, NC, USA
| | | | | | | | | | | |
Collapse
|
38
|
Spoerri L, Beaumont KA, Anfosso A, Murphy RJ, Browning AP, Gunasingh G, Haass NK. Real-Time Cell Cycle Imaging in a 3D Cell Culture Model of Melanoma, Quantitative Analysis, Optical Clearing, and Mathematical Modeling. Methods Mol Biol 2024; 2764:291-310. [PMID: 38393602 DOI: 10.1007/978-1-0716-3674-9_19] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2024]
Abstract
Aberrant cell cycle progression is a hallmark of solid tumors. Therefore, cell cycle analysis is an invaluable technique to study cancer cell biology. However, cell cycle progression has been most commonly assessed by methods that are limited to temporal snapshots or that lack spatial information. In this chapter, we describe a technique that allows spatiotemporal real-time tracking of cell cycle progression of individual cells in a multicellular context. The power of this system lies in the use of 3D melanoma spheroids generated from melanoma cells engineered with the fluorescent ubiquitination-based cell cycle indicator (FUCCI). This technique, combined with mathematical modeling, allows us to gain further and more detailed insight into several relevant aspects of solid cancer cell biology, such as tumor growth, proliferation, invasion, and drug sensitivity.
Collapse
Affiliation(s)
- Loredana Spoerri
- Frazer Institute, The University of Queensland, Brisbane, QLD, Australia
| | - Kimberley A Beaumont
- The Centenary Institute, Sydney, NSW, Australia
- Uniquest, The University of Queensland, Brisbane, QLD, Australia
| | | | - Ryan J Murphy
- Mathematical Sciences, Queensland University of Technology, Brisbane, QLD, Australia
| | - Alexander P Browning
- Mathematical Sciences, Queensland University of Technology, Brisbane, QLD, Australia
| | - Gency Gunasingh
- Frazer Institute, The University of Queensland, Brisbane, QLD, Australia
| | - Nikolas K Haass
- Frazer Institute, The University of Queensland, Brisbane, QLD, Australia.
- The Centenary Institute, Sydney, NSW, Australia.
| |
Collapse
|
39
|
Chew YH, Marucci L. Mechanistic Model-Driven Biodesign in Mammalian Synthetic Biology. Methods Mol Biol 2024; 2774:71-84. [PMID: 38441759 DOI: 10.1007/978-1-0716-3718-0_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Mathematical modeling plays a vital role in mammalian synthetic biology by providing a framework to design and optimize design circuits and engineered bioprocesses, predict their behavior, and guide experimental design. Here, we review recent models used in the literature, considering mathematical frameworks at the molecular, cellular, and system levels. We report key challenges in the field and discuss opportunities for genome-scale models, machine learning, and cybergenetics to expand the capabilities of model-driven mammalian cell biodesign.
Collapse
Affiliation(s)
- Yin Hoon Chew
- School of Mathematics, University of Birmingham, Birmingham, UK
| | - Lucia Marucci
- Department of Engineering Mathematics, University of Bristol, Bristol, UK.
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, UK.
| |
Collapse
|
40
|
Seo S, Lee JW. Applications of Big Data and AI-Driven Technologies in CADD (Computer-Aided Drug Design). Methods Mol Biol 2024; 2714:295-305. [PMID: 37676605 DOI: 10.1007/978-1-0716-3441-7_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
In the field of computer-aided drug design (CADD), there has been dramatic progress in the development of big data and AI-driven methodologies. The expensive and time-consuming process of drug design is related to biomedical complexity. CADD can be used to apply effective and efficient strategies to overcome obstacles in the field of drug design in order to properly design and develop a new medicine. To prepare the raw data for consistent and repeatable applications of big data and AI methodologies, data pre-processing methods are introduced. Big data and AI technologies can be used to develop drugs in areas including predicting absorption, distribution, metabolism, excretion, and toxicity properties as well as finding binding sites in target proteins and conducting structure-based virtual screenings. The accurate and thorough analysis of large amounts of biomedical data as well as the design of prediction models in the area of drug design is made possible by data pre-processing and applications of big data and AI skills. In the biomedical big data era, knowledge on the biological, chemical, or pharmacological structures of biomedical entities relevant to drug design should be analyzed with significant big data and AI approaches.
Collapse
Affiliation(s)
- Seongmin Seo
- Department of Mechanical Engineering, Ulsan National Institute of Science and Technology, Ulsan, Republic of Korea
| | - Jai Woo Lee
- Department of Big Data Science, College of Public Policy, Korea University, Sejong, Republic of Korea.
| |
Collapse
|
41
|
Jalili-Jahani N, Abdollahi H, Rabbani F. Spectrophotometric thermodynamic analysis of methylene blue aggregation in water by chemometrics; what comes in the presence of different cyclodextrins? Spectrochim Acta A Mol Biomol Spectrosc 2023; 303:123126. [PMID: 37506453 DOI: 10.1016/j.saa.2023.123126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/28/2022] [Revised: 06/28/2023] [Accepted: 07/08/2023] [Indexed: 07/30/2023]
Abstract
Spectrophotometry has been utilized to characterize the thermodynamic/dynamic properties of self-aggregation of methylene blue (MB) in water, particularly while interacting with a modulator like different cyclodextrins (α-, β-, hydroxypropyl-β- (HP-β-), and γ-CDs). These systems comprise many interactions that make such chemical systems sophisticated. We developed a mathematical modeling-fitting analysis for the simultaneous quantitative analysis of thermodynamic parameters of chemical reactions, relying on the fitting algorithm. Through analyzing simulated photometric titration data, we demonstrate the simultaneous determination of thermodynamic parameters of the different guest/host interactions. This first has brought the need for the calculation of the visible-light absorption spectrum and the thermodynamic parameters for the pure dimerization system. Therefore, the multiwavelength spectral-mole ratio data of aqueous solutions of MB over a concentration range of 2.5 × 10-5 to 4.5 × 10-5 M while temperature is changing; or being titrated with CDs solutions at various temperatures were collected, augmented, and then have been fed to solid mathematical routines to determine the potential existence of dimeric aggregates. The results of thermodynamics indicated that the positions of the monomer/dimer equilibria do not alter by the presence of α-CD. The apparent dimerization was suppressed upon addition of β- or HP-β-CDs, while the addition of γ-CD enhanced the dimerization.
Collapse
Affiliation(s)
- Nasser Jalili-Jahani
- Green Land Shiraz Eksir Chemical and Agricultural Industries Company, Shiraz 7137753451, Iran; Department of Chemistry, Idaho State University, Pocatello, ID 83209, United States
| | - Hamid Abdollahi
- Institute for Advanced Studies in Basic Sciences (IASBS), Department of Chemistry, Zanjan 4513766731, Iran.
| | - Faride Rabbani
- Institute for Advanced Studies in Basic Sciences (IASBS), Department of Chemistry, Zanjan 4513766731, Iran; Monitoring the Human Hygience Condition and Standard Qeshm, Qeshm Island, 7951919697, Iran.
| |
Collapse
|
42
|
Shekhtman L, Duehren S, Etzion O, Cotler SJ, Dahari H. Hepatitis D Virus and HBsAg Dynamics in the era of new Antiviral Treatments. Curr Gastroenterol Rep 2023; 25:401-412. [PMID: 37819559 PMCID: PMC10842234 DOI: 10.1007/s11894-023-00901-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/19/2023] [Indexed: 10/13/2023]
Abstract
PURPOSE OF REVIEW Hepatitis D virus (HDV) infection is the most severe form of chronic viral hepatitis, with no FDA-approved therapy. Progress in the development of effective HDV treatments is accelerating. This review highlights how mathematical modeling is improving understanding of HDV-HBsAg-host dynamics during antiviral therapy and generating insights into the efficacy and modes of action (MOA) of new antiviral agents. RECENT FINDINGS Clinical trials with pegylated-interferon-λ, bulevertide, nucleic acid polymers, and/or lonafarnib against various steps of the HDV-life cycle have revealed new viral-kinetic patterns that were not observed under standard treatment with pegylated-interferon-α. Modeling indicated that the half-lives of circulating HDV and HBsAg are ~ 1.7 d and ~ 1.3 d, respectively, estimated the relative response of HDV and HBsAg during different antiviral therapies, and provided insights into the efficacy and MOA of drugs in development for treating HDV, which can inform response-guided therapy to individualize treatment duration. Mathematical modeling of HDV and HBsAg kinetics provides a window into the HDV virus lifecycle, HDV-HBsAg-host dynamics during antiviral therapy, and the MOA of new drugs for HDV.
Collapse
Affiliation(s)
- Louis Shekhtman
- The Program for Experimental & Theoretical Modeling, Division of Hepatology, Department of Medicine, Stritch School of Medicine, Loyola University Medical Center, Maywood, IL, USA
- Department of Information Science, Bar-Ilan University, Ramat Gan, Israel
| | - Sarah Duehren
- The Program for Experimental & Theoretical Modeling, Division of Hepatology, Department of Medicine, Stritch School of Medicine, Loyola University Medical Center, Maywood, IL, USA
| | - Ohad Etzion
- Department of Gastroenterology and Liver Diseases, Soroka University Medical Center, Beer-Sheva, Israel
| | - Scott J Cotler
- The Program for Experimental & Theoretical Modeling, Division of Hepatology, Department of Medicine, Stritch School of Medicine, Loyola University Medical Center, Maywood, IL, USA
| | - Harel Dahari
- The Program for Experimental & Theoretical Modeling, Division of Hepatology, Department of Medicine, Stritch School of Medicine, Loyola University Medical Center, Maywood, IL, USA.
| |
Collapse
|
43
|
Savagar B, Jones BA, Arnold M, Walker M, Fournié G. Modelling flock heterogeneity in the transmission of peste des petits ruminants virus and its impact on the effectiveness of vaccination for eradication. Epidemics 2023; 45:100725. [PMID: 37935076 DOI: 10.1016/j.epidem.2023.100725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Revised: 09/29/2023] [Accepted: 10/30/2023] [Indexed: 11/09/2023] Open
Abstract
Peste des petits ruminants (PPR) is an acute infectious disease of small ruminants targeted for global eradication by 2030. The Global Strategy for Control and Eradication (GSCE) recommends mass vaccination targeting 70% coverage of small ruminant populations in PPR-endemic regions. These small ruminant populations are diverse with heterogeneous mixing patterns that may influence PPR virus (PPRV) transmission dynamics. This paper evaluates the impact of heterogeneous mixing on (i) PPRV transmission and (ii) the likelihood of different vaccination strategies achieving PPRV elimination, including the GSCE recommended strategy. We develop models simulating heterogeneous transmission between hosts, including a metapopulation model of PPRV transmission between villages in lowland Ethiopia fitted to serological data. Our results demonstrate that although heterogeneous mixing of small ruminant populations increases the instability of PPRV transmission-increasing the chance of fadeout in the absence of intervention-a vaccination coverage of 70% may be insufficient to achieve elimination if high-risk populations are not targeted. Transmission may persist despite very high vaccination coverage (>90% small ruminants) if vaccination is biased towards more accessible but lower-risk populations such as sedentary small ruminant flocks. These results highlight the importance of characterizing small ruminant mobility patterns and identifying high-risk populations for vaccination and support a move towards targeted, risk-based vaccination programmes in the next phase of the PPRV eradication programme. Our modelling approach also illustrates a general framework for incorporating heterogeneous mixing patterns into models of directly transmitted infectious diseases where detailed contact data are limited. This study improves understanding of PPRV transmission and elimination in heterogeneous small ruminant populations and should be used to inform and optimize the design of PPRV vaccination programmes.
Collapse
Affiliation(s)
- Bethan Savagar
- Veterinary Epidemiology, Economics and Public Health Group, WOAH Collaborating Centre for Risk Analysis and Modelling, Department of Pathobiology and Population Sciences, The Royal Veterinary College, London, UK.
| | - Bryony A Jones
- Department of Epidemiological Sciences, WOAH Collaborating Centre in Risk Analysis and Modelling, Animal and Plant Health Agency (APHA), Addlestone, Surrey, UK
| | - Mark Arnold
- Department of Epidemiological Sciences, WOAH Collaborating Centre in Risk Analysis and Modelling, Animal and Plant Health Agency (APHA), Addlestone, Surrey, UK
| | - Martin Walker
- Veterinary Epidemiology, Economics and Public Health Group, WOAH Collaborating Centre for Risk Analysis and Modelling, Department of Pathobiology and Population Sciences, The Royal Veterinary College, London, UK; London Centre for Neglected Tropical Disease Research, Department of Infectious Disease Epidemiology, Imperial College London, UK
| | - Guillaume Fournié
- Veterinary Epidemiology, Economics and Public Health Group, WOAH Collaborating Centre for Risk Analysis and Modelling, Department of Pathobiology and Population Sciences, The Royal Veterinary College, London, UK; Université de Lyon, INRAE, VetAgro Sup, UMR EPIA, Marcy l'Etoile, France; Université Clermont Auvergne, INRAE, VetAgro Sup, UMR EPIA, Saint Genes Champanelle, France
| |
Collapse
|
44
|
Souza VC, Remus A, Batonon-Alavo DI, Rouffineau F, Mercier Y, Pomar C, Kebreab E. Systematic review and meta-analysis of the methionine utilization efficiency in piglets receiving different methionine sources. Animal 2023; 17 Suppl 5:100894. [PMID: 37482458 DOI: 10.1016/j.animal.2023.100894] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 06/13/2023] [Accepted: 06/15/2023] [Indexed: 07/25/2023] Open
Abstract
Methionine (Met) is an essential amino acid that can be supplied in different chemical forms: DL-Met, L-Met, and OH-Met. This study aimed (i) to model and compare the utilization efficiency of Met for protein deposition (PD) from all sources and (ii) to determine the efficacy and efficiency of these three free Met sources in average daily gain (ADG) of post-weaning pigs fed at or below the Met requirement. A systematic review of the literature resulted in 1 898 papers being screened for title and abstract, with 24 papers meeting the inclusion criteria. The resulting database containing 208 treatment means was used. Prior to model development, the standardized ileal digestible (SID) Met requirements in percentage in the diet were determined using initial and final BW according to the NRC (2012). Data from piglets fed above the SID Met requirements were excluded from the database prior to statistical analysis. Linear mixed-effects regression models predicting ADG as a function of free Met source and SID methionine intake (Meti) or methionine + cysteine intake (Met + cysi) were used to evaluate the efficacy and efficiency of free Met source for weight gain. Moreover, Met retention was modeled assuming that 16% of ADG is deposited as PD, and that Met accounts for 2% of PD. Met utilization efficiency was calculated as Meti after maintenance divided by Met retained in PD. Met utilization efficiency was 77% for the basal diet, decreased (P < 0.01) as Meti increased, and was equal among the three free Met sources. The mixed-effects models showed no difference in ADG for three free Met sources evaluated (P > 0.05). However, the efficacy (ADG per unit of SID Meti) of free Met sources for weight gain differed between piglets fed L and DL-Met (P < 0.05), while there was no difference (P > 0.05) between piglets fed DL and OH-Met or OH and L-Met. On average, piglets fed L-Met gained 40.3 g/d more weight per unit of increase in SID Meti than those fed DL-Met (model 4; P = 0.05). The efficacy of free Met sources for ADG was also compared using SID Met + cysi as covariable. Piglets fed L- (+11.7 g/d; P = 0.02) or OH-Met (+11.5 g/d; P = 0.04) gained more weight per gram of SID Met + cysi compared to those fed DL-Met. In conclusion, although the efficacy of DL- and L-Met for ADG differed, the efficiency for PD of L-, DL-, and OH-Met were not different in piglets fed at or below Meti requirement.
Collapse
Affiliation(s)
- V C Souza
- Department of Animal Science, University of California, Davis, CA 95616, USA
| | - A Remus
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, Québec J1M 0C8, Canada
| | | | | | - Y Mercier
- Adisseo France SAS, Malicorne F-03630, France
| | - C Pomar
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, Québec J1M 0C8, Canada
| | - E Kebreab
- Department of Animal Science, University of California, Davis, CA 95616, USA.
| |
Collapse
|
45
|
Casson CL, John SA, Ferrall-Fairbanks MC. Mathematical modeling of cardio-oncology: Modeling the systemic effects of cancer therapeutics on the cardiovascular system. Semin Cancer Biol 2023; 97:30-41. [PMID: 37979714 DOI: 10.1016/j.semcancer.2023.11.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Revised: 08/25/2023] [Accepted: 11/12/2023] [Indexed: 11/20/2023]
Abstract
Cardiotoxicity is a common side-effect of many cancer therapeutics; however, to-date there has been very little push to understand the mechanisms underlying this group of pathologies. This has led to the emergence of cardio-oncology, a field of medicine focused on understanding the effects of cancer and its treatment on the human heart. Here, we describe how mechanistic modeling approaches have been applied to study open questions in the cardiovascular system and how these approaches are being increasingly applied to advance knowledge of the underlying effects of cancer treatments on the human heart. A variety of mechanistic, mathematical modeling techniques have been applied to explore the link between common cancer treatments, such as chemotherapy, radiation, targeted therapy, and immunotherapy, and cardiotoxicity, nevertheless there is limited coverage in the different types of cardiac dysfunction that may be associated with these treatments. Moreover, cardiac modeling has a rich heritage of mathematical modeling and is well suited for the further development of novel approaches for understanding the cardiotoxicities associated with cancer therapeutics. There are many opportunities to combine mechanistic, bottom-up approaches with data-driven, top-down approaches to improve personalized, precision oncology to better understand, and ultimately mitigate, cardiac dysfunction in cancer patients.
Collapse
Affiliation(s)
- Camara L Casson
- J. Crayton Pruitt Family Department of Biomedical Engineering, University of Florida, Gainesville, FL 32611, USA
| | - Sofia A John
- Department of Statistics, University of Florida, Gainesville, FL 32611, USA
| | - Meghan C Ferrall-Fairbanks
- J. Crayton Pruitt Family Department of Biomedical Engineering, University of Florida, Gainesville, FL 32611, USA; University of Florida Health Cancer Center, University of Florida, Gainesville, FL 32611, USA.
| |
Collapse
|
46
|
Vynnycky E, Knapp JK, Papadopoulos T, Cutts FT, Hachiya M, Miyano S, Reef SE. Estimates of the global burden of Congenital Rubella Syndrome, 1996-2019. Int J Infect Dis 2023; 137:149-156. [PMID: 37690575 PMCID: PMC10689248 DOI: 10.1016/j.ijid.2023.09.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Revised: 08/18/2023] [Accepted: 09/05/2023] [Indexed: 09/12/2023] Open
Abstract
OBJECTIVES Many countries introduced rubella-containing vaccination (RCV) after 2011, following changes in recommended World Health Organization (WHO) vaccination strategies and external support. We evaluated the impact of these introductions. METHODS We estimated the country-specific, region-specific, and global Congenital Rubella Syndrome (CRS) incidence during 1996-2019 using mathematical modeling, including routine and campaign vaccination coverage and seroprevalence data. RESULTS In 2019, WHO African and Eastern Mediterranean regions had the highest estimated CRS incidence (64 [95% confidence intervals (CI): 24-123] and 27 [95% CI: 4-67] per 100,000 live births respectively), where nearly half of births occur in countries that have introduced RCV. Other regions, where >95% of births occurred in countries that had introduced RCV, had a low estimated CRS incidence (<1 [95% CI: <1 to 8] and <1 [95% CI: <1 to 12] per 100,000 live births in South-East Asia [SEAR] and the Western Pacific [WPR] respectively, and similarly in Europe and the Americas). The estimated number of CRS births globally declined by approximately two-thirds during 2010-2019, from 100,000 (95% CI: 54,000-166,000) to 32,000 (95% CI: 13,000-60,000), representing a 73% reduction since 1996, largely following RCV introductions in WPR and SEAR, where the greatest reductions occurred. CONCLUSIONS Further reductions can occur by introducing RCV in remaining countries and maintaining high RCV coverage.
Collapse
Affiliation(s)
- Emilia Vynnycky
- Statistics Modelling and Economics Department, United Kingdom Health Security Agency, London, UK; Department of Infectious Disease Epidemiology, London School of Hygiene & Tropical Medicine, London, UK; TB Modelling Group and Centre for Mathematical Modelling of Infectious Diseases, London School of Hygiene & Tropical Medicine, London, UK.
| | - Jennifer K Knapp
- Global Immunization Division, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
| | - Timos Papadopoulos
- Statistics Modelling and Economics Department, United Kingdom Health Security Agency, London, UK
| | - Felicity T Cutts
- Department of Infectious Disease Epidemiology, London School of Hygiene & Tropical Medicine, London, UK
| | - Masahiko Hachiya
- Bureau of International Health Cooperation, National Center for Global Health and Medicine, Toyama, Shinjuku-ku, Tokyo, Japan
| | - Shinsuke Miyano
- Bureau of International Health Cooperation, National Center for Global Health and Medicine, Toyama, Shinjuku-ku, Tokyo, Japan
| | - Susan E Reef
- Global Immunization Division, Centers for Disease Control and Prevention, Atlanta, Georgia, USA
| |
Collapse
|
47
|
Andrade DLLS, Pintarelli GB, Rosa JV, Paro IB, Pagano PJT, Silva JCN, Suzuki DOH. Musa acuminata as electroporation model. Bioelectrochemistry 2023; 154:108549. [PMID: 37639773 DOI: 10.1016/j.bioelechem.2023.108549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2023] [Revised: 08/17/2023] [Accepted: 08/19/2023] [Indexed: 08/31/2023]
Abstract
Electrochemotherapy (ECT) and Irreversible electroporation (IRE) are cancer treatments based on electric field distribution in tissues. Solanum tuberosum (potato tissue) phantom is known to mimic changes in the electrical conductivity that occur in animal tissues during electroporation (EP). Electric field distribution is assessed through enzymatic staining. However, the 24-h wait for this assessment could slow agile response scenarios. We developed and validated the Musa acuminata (cavendish banana) conductivity model, which quickly evaluates EP by tissue staining. We investigated the frequency response of the tissue using impedance spectroscopy analysis, conductivity changes, and enzymatic staining. We optimized three usual EP models: adapted Gompertz, smoothed Heaviside, and the sigmoid or logistic function. We found dielectric parameters in banana tissue similar to those in potato (electrical conductivity of 0.035 S/m and relative permittivity of 4.1×104). The coefficients of determination R2 were 99.94% (Gompertz), 99.85% (Heaviside), and 99.58% (sigmoid). The sigmoid and Heaviside functions described the calibration and validation electric currents with 95% confidence. We observed the electroporated areas in bananas 3h30m after EP. Staining was significant after 450 V/cm. The conductivity model of Musa acuminata suits treatment planning, hardware development, and training scenarios. Banana phantom supports the 3Rs practice and is a reliable alternative for potato in EP studies.
Collapse
Affiliation(s)
- Daniella L L S Andrade
- Institute of Biomedical Engineering, Department of Electrical and Electronics Engineering, Federal University of Santa Catarina, Florianópolis, SC, Brazil
| | - Guilherme B Pintarelli
- Department of Control and Automation Engineering, Federal University of Santa Catarina, Blumenau, SC, Brazil
| | - Juliana V Rosa
- Institute of Biomedical Engineering, Department of Electrical and Electronics Engineering, Federal University of Santa Catarina, Florianópolis, SC, Brazil
| | - Isabela B Paro
- Institute of Biomedical Engineering, Department of Electrical and Electronics Engineering, Federal University of Santa Catarina, Florianópolis, SC, Brazil
| | - Pedro J T Pagano
- Institute of Biomedical Engineering, Department of Electrical and Electronics Engineering, Federal University of Santa Catarina, Florianópolis, SC, Brazil
| | - Julia C N Silva
- Institute of Biomedical Engineering, Department of Electrical and Electronics Engineering, Federal University of Santa Catarina, Florianópolis, SC, Brazil
| | - Daniela O H Suzuki
- Institute of Biomedical Engineering, Department of Electrical and Electronics Engineering, Federal University of Santa Catarina, Florianópolis, SC, Brazil
| |
Collapse
|
48
|
Kiti MC, Aguolu OG, Zelaya A, Chen HY, Ahmed N, Batross J, Liu CY, Nelson KN, Jenness SM, Melegaro A, Ahmed F, Malik F, Omer SB, Lopman BA. Changing social contact patterns among US workers during the COVID-19 pandemic: April 2020 to December 2021. Epidemics 2023; 45:100727. [PMID: 37948925 PMCID: PMC10730080 DOI: 10.1016/j.epidem.2023.100727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Revised: 10/21/2023] [Accepted: 11/03/2023] [Indexed: 11/12/2023] Open
Abstract
Non-pharmaceutical interventions minimize social contacts, hence the spread of respiratory pathogens such as influenza and SARS-CoV-2. Globally, there is a paucity of social contact data from the workforce. In this study, we quantified two-day contact patterns among USA employees. Contacts were defined as face-to-face conversations, involving physical touch or proximity to another individual and were collected using electronic self-kept diaries. Data were collected over 4 rounds from 2020 to 2021 during the COVID-19 pandemic. Mean (standard deviation) contacts reported by 1456 participants were 2.5 (2.5), 8.2 (7.1), 9.2 (7.1) and 10.1 (9.5) across round 1 (April-June 2020), 2 (November 2020-January 2021), 3 (June-August 2021), and 4 (November-December 2021), respectively. Between round 1 and 2, we report a 3-fold increase in the mean number of contacts reported per participant with no major increases from round 2-4. We then modeled SARS-CoV-2 transmission at home, work, and community settings. The model revealed reduced relative transmission in all settings in round 1. Subsequently, transmission increased at home and in the community but remained exceptionally low in work settings. To accurately parameterize models of infection transmission and control, we need empirical social contact data that capture human mixing behavior across time.
Collapse
Affiliation(s)
- Moses C Kiti
- Rollins School of Public Health, Emory University, GA, USA.
| | - Obianuju G Aguolu
- Yale Institute for Global Health, Yale University, CT, USA; Yale School of Medicine, Yale University, CT, USA
| | - Alana Zelaya
- Rollins School of Public Health, Emory University, GA, USA
| | - Holin Y Chen
- Rollins School of Public Health, Emory University, GA, USA
| | - Noureen Ahmed
- Yale Institute for Global Health, Yale University, CT, USA
| | | | - Carol Y Liu
- Rollins School of Public Health, Emory University, GA, USA
| | | | | | - Alessia Melegaro
- DONDENA Centre for Research in Social Dynamics and Public Policy, Bocconi University, Italy
| | - Faruque Ahmed
- Division of Global Migration and Quarantine, Centers for Disease Control and Prevention, Atlanta, GA, USA
| | - Fauzia Malik
- Yale Institute for Global Health, Yale University, CT, USA
| | - Saad B Omer
- Yale Institute for Global Health, Yale University, CT, USA; Yale School of Medicine, Yale University, CT, USA
| | - Ben A Lopman
- Rollins School of Public Health, Emory University, GA, USA
| |
Collapse
|
49
|
Saldaña F, Kebir A, Camacho-Gutiérrez JA, Aguiar M. Optimal vaccination strategies for a heterogeneous population using multiple objectives: The case of L 1- and L 2-formulations. Math Biosci 2023; 366:109103. [PMID: 37918477 DOI: 10.1016/j.mbs.2023.109103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Revised: 09/29/2023] [Accepted: 10/28/2023] [Indexed: 11/04/2023]
Abstract
The choice of the objective functional in optimization problems coming from biomedical and epidemiological applications plays a key role in optimal control outcomes. In this study, we investigate the role of the objective functional on the structure of the optimal control solution for an epidemic model for sexually transmitted infections that includes a core group with higher sexual activity levels than the rest of the population. An optimal control problem is formulated to find a targeted vaccination program able to control the spread of the infection with minimum vaccine deployment. Both L1- and L2-objectives are considered as an attempt to explore the trade-offs between control dynamics and the functional form characterizing optimality. The results show that the optimal vaccination policies for both the L1- and the L2-formulation share one important qualitative property, that is, immunization of the core group should be prioritized by policymakers to achieve a fast reduction of the epidemic. However, quantitative aspects of this result can be significantly affected depending on the choice of the control weights between formulations. Overall, the results suggest that with appropriate weight constants, the optimal control outcomes are reasonably robust with respect to the L1- or L2-formulation. This is particularly true when the monetary cost of the control policy is substantially lower than the cost associated with the disease burden. Under these conditions, even if the L1-formulation is more realistic from a modeling perspective, the L2-formulation can be used as an approximation and yield qualitatively comparable outcomes.
Collapse
Affiliation(s)
| | - Amira Kebir
- Basque Center for Applied Mathematics (BCAM), Bilbao, Spain; IPEIT, Tunis University, Tunis, Tunisia; BIMS-IPT, Tunis El Manar University, Tunis, Tunisia
| | | | - Maíra Aguiar
- Basque Center for Applied Mathematics (BCAM), Bilbao, Spain; Ikerbasque, Basque Foundation for Science, Bilbao, Spain; Dipartimento di Matematica, Università degli Studi di Trento, Trento, Italy
| |
Collapse
|
50
|
Lucey PM, Rossow HA. Describing the distribution type of DM intake for dairy cow pens based on pen characteristics. Animal 2023; 17 Suppl 5:100888. [PMID: 37451902 DOI: 10.1016/j.animal.2023.100888] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Revised: 06/11/2023] [Accepted: 06/12/2023] [Indexed: 07/18/2023] Open
Abstract
In practice, cows are fed by pen, but a diet is formulated to the nutrient requirements of a single cow. If the DM intake (DMI) of a pen were equal for all cows, this approach would have no error, but cows are grouped into pens on pregnancy and other management factors creating a distribution of DMI. The goal of precision feeding is to meet the requirements of individual animals to increase efficiency and reduce environmental impact but is not achieved when a group is fed as if the individuals have uniform requirements and the DMI distribution is not normal. The hypothesis of this work is that the DMI of cow pens are not normally distributed and the total DMI from the best-fit distribution shape for a cow pen will have lower percentage error to the observed DMI than a prediction of a single DMI that is fed at a uniform level and assumes a normal distribution. Our objective was to describe the distribution shape of DMI by week of lactation, and for different pen types. Pens were generated by randomly assorting cows by the week of lactation from a database into different categories of pen for size and lactation period. These pens were fitted to the best distribution type, and its parameters were used to randomly generate distribution plots that predict the total DMI for each pen. A second predictive model estimated the DMI of each pen using an empirical equation of DMI that was multiplied by the number of cows in the pen to represent feeding of a uniform DMI quantity. The percentage error for the distribution shape model was significantly lower than the empirical model with pen errors being less than 1%. The beta distribution type was the most common distribution to best represent the data of pen DMI. Describing the distribution and using it to predict a total pen DMI provides accurate estimates of feed quantity for a group. Reducing error by using the distribution of DMI for feed formulation, instead of the nutrient requirements of an individual animal can provide a precision nutrition approach to group feeding.
Collapse
Affiliation(s)
- P M Lucey
- Veterinary Medicine Teaching and Research Center, School of Veterinary Medicine, University of California Davis, Tulare 93274, USA
| | - H A Rossow
- Veterinary Medicine Teaching and Research Center, School of Veterinary Medicine, University of California Davis, Tulare 93274, USA.
| |
Collapse
|