1
|
Alamoudi E, Reck F, Bundgaard N, Graw F, Brusch L, Hasenauer J, Schälte Y. A wall-time minimizing parallelization strategy for approximate Bayesian computation. PLoS One 2024; 19:e0294015. [PMID: 38386671 PMCID: PMC10883530 DOI: 10.1371/journal.pone.0294015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Accepted: 10/24/2023] [Indexed: 02/24/2024] Open
Abstract
Approximate Bayesian Computation (ABC) is a widely applicable and popular approach to estimating unknown parameters of mechanistic models. As ABC analyses are computationally expensive, parallelization on high-performance infrastructure is often necessary. However, the existing parallelization strategies leave computing resources unused at times and thus do not optimally leverage them yet. We present look-ahead scheduling, a wall-time minimizing parallelization strategy for ABC Sequential Monte Carlo algorithms, which avoids idle times of computing units by preemptive sampling of subsequent generations. This allows to utilize all available resources. The strategy can be integrated with e.g. adaptive distance function and summary statistic selection schemes, which is essential in practice. Our key contribution is the theoretical assessment of the strategy of preemptive sampling and the proof of unbiasedness. Complementary, we provide an implementation and evaluate the strategy on different problems and numbers of parallel cores, showing speed-ups of typically 10-20% and up to 50% compared to the best established approach, with some variability. Thus, the proposed strategy allows to improve the cost and run-time efficiency of ABC methods on high-performance infrastructure.
Collapse
Affiliation(s)
- Emad Alamoudi
- Life and Medical Sciences Institute, University of Bonn, Bonn, Germany
| | - Felipe Reck
- Life and Medical Sciences Institute, University of Bonn, Bonn, Germany
| | - Nils Bundgaard
- BioQuant—Center for Quantitative Biology, Heidelberg University, Heidelberg, Germany
| | - Frederik Graw
- BioQuant—Center for Quantitative Biology, Heidelberg University, Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
- Department of Medicine 5, Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany
| | - Lutz Brusch
- Center of Information Services and High Performance Computing (ZIH), Technische Universität Dresden, Dresden, Germany
| | - Jan Hasenauer
- Life and Medical Sciences Institute, University of Bonn, Bonn, Germany
- Helmholtz Zentrum München, Institute of Computational Biology, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| | - Yannik Schälte
- Life and Medical Sciences Institute, University of Bonn, Bonn, Germany
- Helmholtz Zentrum München, Institute of Computational Biology, Neuherberg, Germany
- Center for Mathematics, Technische Universität München, Garching, Germany
| |
Collapse
|
2
|
Alamoudi E, Schälte Y, Müller R, Starruß J, Bundgaard N, Graw F, Brusch L, Hasenauer J. FitMultiCell: simulating and parameterizing computational models of multi-scale and multi-cellular processes. Bioinformatics 2023; 39:btad674. [PMID: 37947308 PMCID: PMC10666203 DOI: 10.1093/bioinformatics/btad674] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Revised: 10/25/2023] [Accepted: 11/07/2023] [Indexed: 11/12/2023] Open
Abstract
MOTIVATION Biological tissues are dynamic and highly organized. Multi-scale models are helpful tools to analyse and understand the processes determining tissue dynamics. These models usually depend on parameters that need to be inferred from experimental data to achieve a quantitative understanding, to predict the response to perturbations, and to evaluate competing hypotheses. However, even advanced inference approaches such as approximate Bayesian computation (ABC) are difficult to apply due to the computational complexity of the simulation of multi-scale models. Thus, there is a need for a scalable pipeline for modeling, simulating, and parameterizing multi-scale models of multi-cellular processes. RESULTS Here, we present FitMultiCell, a computationally efficient and user-friendly open-source pipeline that can handle the full workflow of modeling, simulating, and parameterizing for multi-scale models of multi-cellular processes. The pipeline is modular and integrates the modeling and simulation tool Morpheus and the statistical inference tool pyABC. The easy integration of high-performance infrastructure allows to scale to computationally expensive problems. The introduction of a novel standard for the formulation of parameter inference problems for multi-scale models additionally ensures reproducibility and reusability. By applying the pipeline to multiple biological problems, we demonstrate its broad applicability, which will benefit in particular image-based systems biology. AVAILABILITY AND IMPLEMENTATION FitMultiCell is available open-source at https://gitlab.com/fitmulticell/fit.
Collapse
Affiliation(s)
- Emad Alamoudi
- Life and Medical Sciences Institute, University of Bonn, Bonn 53113, Germany
| | - Yannik Schälte
- Life and Medical Sciences Institute, University of Bonn, Bonn 53113, Germany
- Institute of Computational Biology, Helmholtz Zentrum München—German Research Center for Environmental Health, Neuherberg 85764, Germany
- Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, Technische Universität München, Garching 85748, Germany
| | - Robert Müller
- Center of Information Services and High Performance Computing (ZIH), Technische Universität Dresden, Dresden 01062, Germany
| | - Jörn Starruß
- Center of Information Services and High Performance Computing (ZIH), Technische Universität Dresden, Dresden 01062, Germany
| | - Nils Bundgaard
- BioQuant—Center for Quantitative Biology, Heidelberg University, Heidelberg 69120, Germany
| | - Frederik Graw
- BioQuant—Center for Quantitative Biology, Heidelberg University, Heidelberg 69120, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg 69120, Germany
- Department of Medicine 5, Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Lutz Brusch
- Center of Information Services and High Performance Computing (ZIH), Technische Universität Dresden, Dresden 01062, Germany
| | - Jan Hasenauer
- Life and Medical Sciences Institute, University of Bonn, Bonn 53113, Germany
- Institute of Computational Biology, Helmholtz Zentrum München—German Research Center for Environmental Health, Neuherberg 85764, Germany
- Center for Mathematics, Chair of Mathematical Modeling of Biological Systems, Technische Universität München, Garching 85748, Germany
| |
Collapse
|
3
|
Raach B, Bundgaard N, Haase MJ, Starruß J, Sotillo R, Stanifer ML, Graw F. Influence of cell type specific infectivity and tissue composition on SARS-CoV-2 infection dynamics within human airway epithelium. PLoS Comput Biol 2023; 19:e1011356. [PMID: 37566610 PMCID: PMC10446191 DOI: 10.1371/journal.pcbi.1011356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2023] [Revised: 08/23/2023] [Accepted: 07/13/2023] [Indexed: 08/13/2023] Open
Abstract
Human airway epithelium (HAE) represents the primary site of viral infection for SARS-CoV-2. Comprising different cell populations, a lot of research has been aimed at deciphering the major cell types and infection dynamics that determine disease progression and severity. However, the cell type-specific replication kinetics, as well as the contribution of cellular composition of the respiratory epithelium to infection and pathology are still not fully understood. Although experimental advances, including Air-liquid interface (ALI) cultures of reconstituted pseudostratified HAE, as well as lung organoid systems, allow the observation of infection dynamics under physiological conditions in unprecedented level of detail, disentangling and quantifying the contribution of individual processes and cells to these dynamics remains challenging. Here, we present how a combination of experimental data and mathematical modelling can be used to infer and address the influence of cell type specific infectivity and tissue composition on SARS-CoV-2 infection dynamics. Using a stepwise approach that integrates various experimental data on HAE culture systems with regard to tissue differentiation and infection dynamics, we develop an individual cell-based model that enables investigation of infection and regeneration dynamics within pseudostratified HAE. In addition, we present a novel method to quantify tissue integrity based on image data related to the standard measures of transepithelial electrical resistance measurements. Our analysis provides a first aim of quantitatively assessing cell type specific infection kinetics and shows how tissue composition and changes in regeneration capacity, as e.g. in smokers, can influence disease progression and pathology. Furthermore, we identified key measurements that still need to be assessed in order to improve inference of cell type specific infection kinetics and disease progression. Our approach provides a method that, in combination with additional experimental data, can be used to disentangle the complex dynamics of viral infection and immunity within human airway epithelial culture systems.
Collapse
Affiliation(s)
- Benjamin Raach
- BioQuant-Center for Quantitative Biology, Heidelberg University, Heidelberg, Germany
| | - Nils Bundgaard
- BioQuant-Center for Quantitative Biology, Heidelberg University, Heidelberg, Germany
| | - Marika J. Haase
- BioQuant-Center for Quantitative Biology, Heidelberg University, Heidelberg, Germany
| | - Jörn Starruß
- Center for Information Services and High Performance Computing, TU Dresden, Dresden, Germany
| | - Rocio Sotillo
- Division of Molecular Thoracic Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Megan L. Stanifer
- Department of Infectious Diseases, Molecular Virology, University Hospital Heidelberg, Heidelberg, Germany
- University of Florida, College of Medicine, Dept. of Molecular Genetics and Microbiology, Gainesville, Florida, United States of America
| | - Frederik Graw
- BioQuant-Center for Quantitative Biology, Heidelberg University, Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
- Friedrich-Alexander-Universität Erlangen-Nürnberg, Department of Medicine 5, Erlangen, Germany
| |
Collapse
|
4
|
Haun A, Fain B, Dobrovolny HM. Effect of cellular regeneration and viral transmission mode on viral spread. J Theor Biol 2023; 558:111370. [PMID: 36460057 DOI: 10.1016/j.jtbi.2022.111370] [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/21/2022] [Revised: 11/03/2022] [Accepted: 11/18/2022] [Indexed: 12/03/2022]
Abstract
Illness negatively affects all aspects of life and one major cause of illness is viral infections. Some viral infections can last for weeks; others, like influenza (the flu), can resolve quickly. During infections, uninfected cells can replicate in order to replenish the cells that have died due to the virus. Many viral models, especially those for short-lived infections like influenza, tend to ignore cellular regeneration since many think that uncomplicated influenza resolves much faster than cells regenerate. This research accounts for cellular regeneration, using an agent-based framework, and varies the regeneration rate in order to understand how cell regeneration affects viral infection dynamics under assumptions of different modes of transmission. We find that although the general trends in peak viral load, time of viral peak, and chronic viral load as regeneration rate changes are the same for cell-free or cell-to-cell transmission, the changes are more extreme for cell-to-cell transmission due to limited access of infected cells to newly generated cells.
Collapse
Affiliation(s)
- Asher Haun
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, United States of America
| | - Baylor Fain
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, United States of America
| | - Hana M Dobrovolny
- Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, United States of America.
| |
Collapse
|
5
|
Tiwari S, Yadav N, Singh S. COVID-19 and Parkinson's Disease: Possible Links in Pathology and Therapeutics. Neurotox Res 2022. [PMID: 35829997 DOI: 10.1007/s12640-022-00540-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 05/30/2022] [Accepted: 06/28/2022] [Indexed: 02/07/2023]
Abstract
The outbreak of SARs-CoV-2 with emerging new variants is leading to global health crisis and has brought a major concern for patients with comorbidities. Parkinson's disease (PD) is a motor neurodegenerative disease involving various metabolic and psychological ailments along with the common occurrence of hyposmia as observed in COVID-19 patients. In addition, the observed surplus inflammatory responses in both diseases are also alarming. Alongside, angiotensin-converting enzyme 2 (ACE2) receptor, essentially required by SARS-CoV-2 to enter the cell and dopamine decarboxylase (DDC), required for dopamine synthesis is known to co-regulate in the non-neuronal cells. Taken together, these conditions suggested the probable reciprocal pathological relation between COVID-19 and PD and also suggested that during comorbidities, the disease diagnosis and therapeutics are critical and may engender severe health complications. In this review, we discuss various events and mechanisms which may have implications for the exacerbation of PD conditions and must be taken into account during the treatment of patients.
Collapse
|
6
|
Miller J, Burch-Smith TM, Ganusov VV. Mathematical Modeling Suggests Cooperation of Plant-Infecting Viruses. Viruses 2022; 14:741. [PMID: 35458472 PMCID: PMC9029262 DOI: 10.3390/v14040741] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Revised: 03/12/2022] [Accepted: 03/25/2022] [Indexed: 02/05/2023] Open
Abstract
Viruses are major pathogens of agricultural crops. Viral infections often start after the virus enters the outer layer of a tissue, and many successful viruses, after local replication in the infected tissue, are able to spread systemically. Quantitative details of virus dynamics in plants, however, are poorly understood, in part, because of the lack of experimental methods which allow the accurate measurement of the degree of infection in individual plant tissues. Recently, a group of researchers followed the kinetics of infection of individual cells in leaves of Nicotiana tabacum plants using Tobacco etch virus (TEV) expressing either Venus or blue fluorescent protein (BFP). Assuming that viral spread occurs from lower to upper leaves, the authors fitted a simple mathematical model to the frequency of cellular infection by the two viral variants found using flow cytometry. While the original model could accurately describe the kinetics of viral spread locally and systemically, we found that many alternative versions of the model, for example, if viral spread starts at upper leaves and progresses to lower leaves or when virus dissemination is stopped due to an immune response, fit the data with reasonable quality, and yet with different parameter estimates. These results strongly suggest that experimental measurements of the virus infection in individual leaves may not be sufficient to identify the pathways of viral dissemination between different leaves and reasons for viral control. We propose experiments that may allow discrimination between the alternatives. By analyzing the kinetics of coinfection of individual cells by Venus and BFP strains of TEV we found a strong deviation from the random infection model, suggesting cooperation between the two strains when infecting plant cells. Importantly, we showed that many mathematical models on the kinetics of coinfection of cells with two strains could not adequately describe the data, and the best fit model needed to assume (i) different susceptibility of uninfected cells to infection by two viruses locally in the leaf vs. systemically from other leaves, and (ii) decrease in the infection rate depending on the fraction of uninfected cells which could be due to a systemic immune response. Our results thus demonstrate the difficulty in reaching definite conclusions from extensive and yet limited experimental data and provide evidence of potential cooperation between different viral variants infecting individual cells in plants.
Collapse
Affiliation(s)
- Joshua Miller
- Department of Mathematics, University of Tennessee, Knoxville, TN 37996, USA;
| | | | - Vitaly V. Ganusov
- Department of Mathematics, University of Tennessee, Knoxville, TN 37996, USA;
- Department of Microbiology, University of Tennessee, Knoxville, TN 37996, USA
| |
Collapse
|