Overview of the multispecies ovary tissue histology electronic repository (MOTHER).

The Multispecies Ovary Tissue Histology Electronic Repository (MOTHER) is a publicly accessible repository of ovary histology images. MOTHER includes hundreds of images from nonhuman primates, as well as ovary histology images from an expanding range of other species. Along with an image, MOTHER provides metadata about the image, and for selected species, follicle identification annotations. Ongoing work includes assisting scientists with contributing their histology images, creation of manual and automated (via machine learning) processing pipelines to identify and count ovarian follicles in different stages of development, and the incorporation of that data into the MOTHER database (MOTHER-DB). MOTHER will be a critical data repository storing and disseminating high-value histology images that are essential for research into ovarian function, fertility, and intra-species variability.

Morphological basis of the lung adenocarcinoma subtypes.

Lung adenocarcinoma (LUAD), which accounts for a large proportion of lung cancers, is divided into five major subtypes based on histologic characteristics. The clinical characteristics, prognosis, and responses to treatments vary among subtypes. Here, we demonstrate that the variations of cell-cell contact energy result in the LUAD subtype-specific morphogenesis. We reproduced the morphologies of the papillary LUAD subtypes with the cellular Potts Model (CPM). Simulations and experimental validations revealed modifications of cell-cell contact energy changed the morphology from a papillary-like structure to micropapillary or solid subtype-like structures. Remarkably, differential gene expression analysis revealed subtype-specific expressions of genes relating to cell adhesion. Knockdown experiments of the micropapillary upregulated ITGA11 gene resulted in the morphological changes of the spheroids produced from an LUAD cell line PC9. This work shows the consequences of gene mutations and gene expressions on patient prognosis through differences in tissue composing physical forces among LUAD subtypes.

Toward a common standard for data and specimen provenance in life sciences.

Open and practical exchange, dissemination, and reuse of specimens and data have become a fundamental requirement for life sciences research. The quality of the data obtained and thus the findings and knowledge derived is thus significantly influenced by the quality of the samples, the experimental methods, and the data analysis. Therefore, a comprehensive and precise documentation of the pre-analytical conditions, the analytical procedures, and the data processing are essential to be able to assess the validity of the research results. With the increasing importance of the exchange, reuse, and sharing of data and samples, procedures are required that enable cross-organizational documentation, traceability, and non-repudiation. At present, this information on the provenance of samples and data is mostly either sparse, incomplete, or incoherent. Since there is no uniform framework, this information is usually only provided within the organization and not interoperably. At the same time, the collection and sharing of biological and environmental specimens increasingly require definition and documentation of benefit sharing and compliance to regulatory requirements rather than consideration of pure scientific needs. In this publication, we present an ongoing standardization effort to provide trustworthy machine-actionable documentation of the data lineage and specimens. We would like to invite experts from the biotechnology and biomedical fields to further contribute to the standard.

Tissue Forge: Interactive biological and biophysics simulation environment.

Tissue Forge is an open-source interactive environment for particle-based physics, chemistry and biology modeling and simulation. Tissue Forge allows users to create, simulate and explore models and virtual experiments based on soft condensed matter physics at multiple scales, from the molecular to the multicellular, using a simple, consistent interface. While Tissue Forge is designed to simplify solving problems in complex subcellular, cellular and tissue biophysics, it supports applications ranging from classic molecular dynamics to agent-based multicellular systems with dynamic populations. Tissue Forge users can build and interact with models and simulations in real-time and change simulation details during execution, or execute simulations off-screen and/or remotely in high-performance computing environments. Tissue Forge provides a growing library of built-in model components along with support for user-specified models during the development and application of custom, agent-based models. Tissue Forge includes an extensive Python API for model and simulation specification via Python scripts, an IPython console and a Jupyter Notebook, as well as C and C++ APIs for integrated applications with other software tools. Tissue Forge supports installations on 64-bit Windows, Linux and MacOS systems and is available for local installation via conda.

Multiscale Model of Antiviral Timing, Potency, and Heterogeneity Effects on an Epithelial Tissue Patch Infected by SARS-CoV-2.

We extend our established agent-based multiscale computational model of infection of lung tissue by SARS-CoV-2 to include pharmacokinetic and pharmacodynamic models of remdesivir. We model remdesivir treatment for COVID-19; however, our methods are general to other viral infections and antiviral therapies. We investigate the effects of drug potency, drug dosing frequency, treatment initiation delay, antiviral half-life, and variability in cellular uptake and metabolism of remdesivir and its active metabolite on treatment outcomes in a simulated patch of infected epithelial tissue. Non-spatial deterministic population models which treat all cells of a given class as identical can clarify how treatment dosage and timing influence treatment efficacy. However, they do not reveal how cell-to-cell variability affects treatment outcomes. Our simulations suggest that for a given treatment regime, including cell-to-cell variation in drug uptake, permeability and metabolism increase the likelihood of uncontrolled infection as the cells with the lowest internal levels of antiviral act as super-spreaders within the tissue. The model predicts substantial variability in infection outcomes between similar tissue patches for different treatment options. In models with cellular metabolic variability, antiviral doses have to be increased significantly (>50% depending on simulation parameters) to achieve the same treatment results as with the homogeneous cellular metabolism.

Computational modelling of nephron progenitor cell movement and aggregation during kidney organogenesis.

During early kidney organogenesis, nephron progenitor (NP) cells move from the tip to the corner region of the ureteric bud (UB) branches in order to form the pretubular aggregate, the early structure giving rise to nephron formation. NP cells derive from metanephric mesenchymal cells and physically interact with them during the movement. Chemotaxis and cell-cell adhesion differences are believed to drive the cell patterning during this critical period of organogenesis. However, the effect of these forces to the cell patterns and their respective movements are known in limited details. We applied a Cellular Potts Model to explore how these forces and organizations contribute to directed cell movement and aggregation. Model parameters were estimated based on fitting to experimental data obtained in ex vivo kidney explant and dissociation-reaggregation organoid culture studies. Our simulations indicated that optimal enrichment and aggregation of NP cells in the UB corner niche requires chemoattractant secretion from both the UB epithelial cells and the NP cells themselves, as well as differences in cell-cell adhesion energies. Furthermore, NP cells were observed, both experimentally and by modelling, to move at higher speed in the UB corner as compared to the tip region where they originated. The existence of different cell speed domains along the UB was confirmed using self-organizing map analysis. In summary, we saw faster NP cell movements near aggregation. The applicability of Cellular Potts Model approach to simulate cell movement and patterning was found to be good during for this early nephrogenesis process. Further refinement of the model should allow us to recapitulate the effects of developmental changes of cell phenotypes and molecular crosstalk during further organ development.

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations.

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

A computational model of liver tissue damage and repair.

Drug induced liver injury (DILI) and cell death can result from oxidative stress in hepatocytes. An initial pattern of centrilobular damage in the APAP model of DILI is amplified by communication from stressed cells and immune system activation. While hepatocyte proliferation counters cell loss, high doses are still lethal to the tissue. To understand the progression of disease from the initial damage to tissue recovery or death, we computationally model the competing biological processes of hepatocyte proliferation, necrosis and injury propagation. We parametrize timescales of proliferation (α), conversion of healthy to stressed cells (ß) and further sensitization of stressed cells towards necrotic pathways (γ) and model them on a Cellular Automaton (CA) based grid of lattice sites. 1D simulations show that a small α/ß (fast proliferation), combined with a large γ/ß (slow death) have the lowest probabilities of tissue survival. At large α/ß, tissue fate can be described by a critical γ/ß* ratio alone; this value is dependent on the initial amount of damage and proportional to the tissue size N. Additionally, the 1D model predicts a minimum healthy population size below which damage is irreversible. Finally, we compare 1D and 2D phase spaces and discuss outcomes of bistability where either survival or death is possible, and of coexistence where simulated tissue never completely recovers or dies but persists as a mixture of healthy, stressed and necrotic cells. In conclusion, our model sheds light on the evolution of tissue damage or recovery and predicts potential for divergent fates given different rates of proliferation, necrosis, and injury propagation.

Mitochondrial depolarization and repolarization in the early stages of acetaminophen hepatotoxicity in mice.

Mitochondrial injury and depolarization are primary events in acetaminophen hepatotoxicity. Previous studies have shown that restoration of mitochondrial function in surviving hepatocytes, which is critical to recovery, is at least partially accomplished via biogenesis of new mitochondria. However, other studies indicate that mitochondria also have the potential to spontaneously repolarize. Although repolarization was previously observed only at a sub-hepatotoxic dose of acetaminophen, we postulated that mitochondrial repolarization in hepatocytes outside the centrilobular regions of necrosis might contribute to recovery of mitochondrial function following acetaminophen-induced injury. Our studies utilized longitudinal intravital microscopy of millimeter-scale regions of the mouse liver to characterize the spatio-temporal relationship between mitochondrial polarization and necrosis early in acetaminophen-induced liver injury. Treatment of male C57BL/6J mice with a single intraperitoneal 250 mg/kg dose of acetaminophen resulted in hepatotoxicity that was apparent histologically within 2 h of treatment, leading to 20 and 60-fold increases in serum aspartate aminotransferase and alanine aminotransferase, respectively, within 6 h. Intravital microscopy of the livers of mice injected with rhodamine123, TexasRed-dextran, propidium iodide and Hoechst 33342 detected centrilobular foci of necrosis within extended regions of mitochondrial depolarization within 2 h of acetaminophen treatment. Although regions of necrosis were more apparent 6 h after acetaminophen treatment, the vast majority of hepatocytes with depolarized mitochondria did not progress to necrosis, but rather recovered mitochondrial polarization within 6 h. Recovery of mitochondrial function following acetaminophen hepatotoxicity thus involves not only biogenesis of new mitochondria, but also repolarization of existing mitochondria. These studies also revealed a spatial distribution of necrosis and mitochondrial depolarization whose single-cell granularity is inconsistent with the hypothesis that communication between neighboring cells plays an important role in the propagation of necrosis during the early stages of APAP hepatotoxicity. Small islands of healthy, intact cells were frequently found surrounded by necrotic cells, and small islands of necrotic cells were frequently found surrounded by healthy, intact cells. Time-series studies demonstrated that these "islands", consisting in some cases of single cells, are persistent; over a period of hours, injury does not spread from individual necrotic cells to their neighbors.

Spatial Temporal Analysis of Fieldwise Flow in Microvasculature.

Changes in blood flow velocity and distribution are vital in maintaining tissue and organ perfusion in response to varying cellular needs. Further, appearance of defects in microcirculation can be a primary indicator in the development of multiple pathologies. Advances in optical imaging have made intravital microscopy (IVM) a practical approach, permitting imaging at the cellular and subcellular level in live animals at high-speed over time. Yet, despite the importance of maintaining adequate tissue perfusion, spatial and temporal variability in capillary flow is seldom documented. In the standard approach, a small number of capillary segments are chosen for imaging over a limited time. To comprehensively quantify capillary flow in an unbiased way we developed Spatial Temporal Analysis of Fieldwise Flow (STAFF), a macro for FIJI open-source image analysis software. Using high-speed image sequences of full fields of blood flow within capillaries, STAFF produces images that represent motion over time called kymographs for every time interval for every vascular segment. From the kymographs STAFF calculates velocities from the distance that red blood cells move over time, and outputs the velocity data as a sequence of color-coded spatial maps for visualization and tabular output for quantitative analyses. In normal mouse livers, STAFF analyses quantified profound differences in flow velocity between pericentral and periportal regions within lobules. Even more unexpected are the differences in flow velocity seen between sinusoids that are side by side and fluctuations seen within individual vascular segments over seconds. STAFF is a powerful new tool capable of providing novel insights by enabling measurement of the complex spatiotemporal dynamics of capillary flow.

The 2019 mathematical oncology roadmap.

Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology-defined here simply as the use of mathematics in cancer research-complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between mathematics and cancer.

A simple automated method for continuous fieldwise measurement of microvascular hemodynamics.

Microvascular perfusion dynamics are vital to physiological function and are frequently dysregulated in injury and disease. Typically studies measure microvascular flow in a few selected vascular segments over limited time, failing to capture spatial and temporal variability. To quantify microvascular flow in a more complete and unbiased way we developed STAFF (Spatial Temporal Analysis of Fieldwise Flow), a macro for FIJI open-source image analysis software. Using high-speed microvascular flow movies, STAFF generates kymographs for every time interval for every vascular segment, calculates flow velocities from red blood cell shadow angles, and outputs the data as color-coded velocity map movies and spreadsheets. In untreated mice, analyses demonstrated profound variation even between adjacent sinusoids over seconds. In acetaminophen-treated mice we detected flow reduction localized to pericentral regions. STAFF is a powerful new tool capable of providing novel insights by enabling measurement of the complex spatiotemporal dynamics of microvascular flow.

Modeling of xenobiotic transport and metabolism in virtual hepatic lobule models.

Computational models of normal liver function and xenobiotic induced liver damage are increasingly being used to interpret in vitro and in vivo data and as an approach to the de novo prediction of the liver's response to xenobiotics. The microdosimetry (dose at the level of individual cells) of xenobiotics vary spatially within the liver because of both compound-independent and compound-dependent factors. In this paper, we build model liver lobules to investigate the interplay between vascular structure, blood flow and cellular transport that lead to regional variations in microdosimetry. We then compared simulation results obtained using this complex spatial model with a simpler linear pipe model of a sinusoid and a very simple single box model. We found that variations in diffusive transport, transporter-mediated transport and metabolism, coupled with complex liver sinusoid architecture and blood flow distribution, led to three essential patterns of xenobiotic exposure within the virtual liver lobule: (1) lobular-wise uniform, (2) radially varying and (3) both radially and azimuthally varying. We propose to use these essential patterns of exposure as a reference for selection of model representations when a computational study involves modeling detailed hepatic responses to xenobiotics.

A Liver-Centric Multiscale Modeling Framework for Xenobiotics.

We describe a multi-scale, liver-centric in silico modeling framework for acetaminophen pharmacology and metabolism. We focus on a computational model to characterize whole body uptake and clearance, liver transport and phase I and phase II metabolism. We do this by incorporating sub-models that span three scales; Physiologically Based Pharmacokinetic (PBPK) modeling of acetaminophen uptake and distribution at the whole body level, cell and blood flow modeling at the tissue/organ level and metabolism at the sub-cellular level. We have used standard modeling modalities at each of the three scales. In particular, we have used the Systems Biology Markup Language (SBML) to create both the whole-body and sub-cellular scales. Our modeling approach allows us to run the individual sub-models separately and allows us to easily exchange models at a particular scale without the need to extensively rework the sub-models at other scales. In addition, the use of SBML greatly facilitates the inclusion of biological annotations directly in the model code. The model was calibrated using human in vivo data for acetaminophen and its sulfate and glucuronate metabolites. We then carried out extensive parameter sensitivity studies including the pairwise interaction of parameters. We also simulated population variation of exposure and sensitivity to acetaminophen. Our modeling framework can be extended to the prediction of liver toxicity following acetaminophen overdose, or used as a general purpose pharmacokinetic model for xenobiotics.

Multiscale Modeling in the Clinic: Drug Design and Development.

A wide range of length and time scales are relevant to pharmacology, especially in drug development, drug design and drug delivery. Therefore, multiscale computational modeling and simulation methods and paradigms that advance the linkage of phenomena occurring at these multiple scales have become increasingly important. Multiscale approaches present in silico opportunities to advance laboratory research to bedside clinical applications in pharmaceuticals research. This is achievable through the capability of modeling to reveal phenomena occurring across multiple spatial and temporal scales, which are not otherwise readily accessible to experimentation. The resultant models, when validated, are capable of making testable predictions to guide drug design and delivery. In this review we describe the goals, methods, and opportunities of multiscale modeling in drug design and development. We demonstrate the impact of multiple scales of modeling in this field. We indicate the common mathematical and computational techniques employed for multiscale modeling approaches used in pharmacometric and systems pharmacology models in drug development and present several examples illustrating the current state-of-the-art models for (1) excitable systems and applications in cardiac disease; (2) stem cell driven complex biosystems; (3) nanoparticle delivery, with applications to angiogenesis and cancer therapy; (4) host-pathogen interactions and their use in metabolic disorders, inflammation and sepsis; and (5) computer-aided design of nanomedical systems. We conclude with a focus on barriers to successful clinical translation of drug development, drug design and drug delivery multiscale models.

Formalizing Knowledge in Multi-Scale Agent-Based Simulations.

Multi-scale, agent-based simulations of cellular and tissue biology are increasingly common. These simulations combine and integrate a range of components from different domains. Simulations continuously create, destroy and reorganize constituent elements causing their interactions to dynamically change. For example, the multi-cellular tissue development process coordinates molecular, cellular and tissue scale objects with biochemical, biomechanical, spatial and behavioral processes to form a dynamic network. Different domain specific languages can describe these components in isolation, but cannot describe their interactions. No current programming language is designed to represent in human readable and reusable form the domain specific knowledge contained in these components and interactions. We present a new hybrid programming language paradigm that naturally expresses the complex multi-scale objects and dynamic interactions in a unified way and allows domain knowledge to be captured, searched, formalized, extracted and reused.

Toxicokinetic Triage for Environmental Chemicals.

Toxicokinetic (TK) models link administered doses to plasma, blood, and tissue concentrations. High-throughput TK (HTTK) performs in vitro to in vivo extrapolation to predict TK from rapid in vitro measurements and chemical structure-based properties. A significant toxicological application of HTTK has been "reverse dosimetry," in which bioactive concentrations from in vitro screening studies are converted into in vivo doses (mg/kg BW/day). These doses are predicted to produce steady-state plasma concentrations that are equivalent to in vitro bioactive concentrations. In this study, we evaluate the impact of the approximations and assumptions necessary for reverse dosimetry and develop methods to determine whether HTTK tools are appropriate or may lead to false conclusions for a particular chemical. Based on literature in vivo data for 87 chemicals, we identified specific properties (eg, in vitro HTTK data, physico-chemical descriptors, and predicted transporter affinities) that correlate with poor HTTK predictive ability. For 271 chemicals we developed a generic HT physiologically based TK (HTPBTK) model that predicts non-steady-state chemical concentration time-courses for a variety of exposure scenarios. We used this HTPBTK model to find that assumptions previously used for reverse dosimetry are usually appropriate, except most notably for highly bioaccumulative compounds. For the thousands of man-made chemicals in the environment that currently have no TK data, we propose a 4-element framework for chemical TK triage that can group chemicals into 7 different categories associated with varying levels of confidence in HTTK predictions. For 349 chemicals with literature HTTK data, we differentiated those chemicals for which HTTK approaches are likely to be sufficient, from those that may require additional data.

The cell behavior ontology: describing the intrinsic biological behaviors of real and model cells seen as active agents.

MOTIVATION: Currently, there are no ontologies capable of describing both the spatial organization of groups of cells and the behaviors of those cells. The lack of a formalized method for describing the spatiality and intrinsic biological behaviors of cells makes it difficult to adequately describe cells, tissues and organs as spatial objects in living tissues, in vitro assays and in computational models of tissues. RESULTS: We have developed an OWL-2 ontology to describe the intrinsic physical and biological characteristics of cells and tissues. The Cell Behavior Ontology (CBO) provides a basis for describing the spatial and observable behaviors of cells and extracellular components suitable for describing in vivo, in vitro and in silico multicell systems. Using the CBO, a modeler can create a meta-model of a simulation of a biological model and link that meta-model to experiment or simulation results. Annotation of a multicell model and its computational representation, using the CBO, makes the statement of the underlying biology explicit. The formal representation of such biological abstraction facilitates the validation, falsification, discovery, sharing and reuse of both models and experimental data. AVAILABILITY AND IMPLEMENTATION: The CBO, developed using Protégé 4, is available at http://cbo.biocomplexity.indiana.edu/cbo/ and at BioPortal (http://bioportal.bioontology.org/ontologies/CBO).