Modeling BK Virus Infection in Renal Transplant Recipients.

Kidney transplant recipients require a lifelong protocol of immunosuppressive therapy to prevent graft rejection. However, these same medications leave them susceptible to opportunistic infections. One pathogen of particular concern is human polyomavirus 1, also known as BK virus (BKPyV).

Evaluating the impact of anatomical and physiological variability on human equivalent doses using PBPK models.

Addressing human anatomical and physiological variability is a crucial component of human health risk assessment of chemicals. Experts have recommended probabilistic chemical risk assessment paradigms in which distributional adjustment factors are used to account for various sources of uncertainty and variability, including variability in the pharmacokinetic behavior of a given substance in different humans. In practice, convenient assumptions about the distribution forms of adjustment factors and human equivalent doses (HEDs) are often used.

Topological and geometric analysis of cell states in single-cell transcriptomic data.

Single-cell RNA sequencing (scRNA-seq) enables dissecting cellular heterogeneity in tissues, resulting in numerous biological discoveries. Various computational methods have been devised to delineate cell types by clustering scRNA-seq data, where clusters are often annotated using prior knowledge of marker genes. In addition to identifying pure cell types, several methods have been developed to identify cells undergoing state transitions, which often rely on prior clustering results.

Quantifying collective motion patterns in mesenchymal cell populations using topological data analysis and agent-based modeling.

Fibroblasts in a confluent monolayer are known to adopt elongated morphologies in which cells are oriented parallel to their neighbors. We collected and analyzed new microscopy movies to show that confluent fibroblasts are motile and that neighboring cells often move in anti-parallel directions in a collective motion phenomenon we refer to as "fluidization" of the cell population. We used machine learning to perform cell tracking for each movie and then leveraged topological data analysis (TDA) to show that time-varying point-clouds generated by the tracks contain significant topological information content that is driven by fluidization, i.

Parameter estimation and identifiability analysis for a bivalent analyte model of monoclonal antibody-antigen binding.

Surface plasmon resonance (SPR) is an extensively used technique to characterize antigen-antibody interactions. Affinity measurements by SPR typically involve testing the binding of antigen in solution to monoclonal antibodies (mAbs) immobilized on a chip and fitting the kinetics data using 1:1 Langmuir binding model to derive rate constants. However, when it is necessary to immobilize antigens instead of the mAbs, a bivalent analyte (1:2) binding model is required for kinetics analysis.

Estimation of Parameter Distributions for Reaction-Diffusion Equations with Competition using Aggregate Spatiotemporal Data.

Reaction-diffusion equations have been used to model a wide range of biological phenomenon related to population spread and proliferation from ecology to cancer. It is commonly assumed that individuals in a population have homogeneous diffusion and growth rates; however, this assumption can be inaccurate when the population is intrinsically divided into many distinct subpopulations that compete with each other. In previous work, the task of inferring the degree of phenotypic heterogeneity between subpopulations from total population density has been performed within a framework that combines parameter distribution estimation with reaction-diffusion models.

Screening cell-cell communication in spatial transcriptomics via collective optimal transport.

Spatial transcriptomic technologies and spatially annotated single-cell RNA sequencing datasets provide unprecedented opportunities to dissect cell-cell communication (CCC). However, incorporation of the spatial information and complex biochemical processes required in the reconstruction of CCC remains a major challenge. Here, we present COMMOT (COMMunication analysis by Optimal Transport) to infer CCC in spatial transcriptomics, which accounts for the competition between different ligand and receptor species as well as spatial distances between cells.

A Comparison of Mathematical and Statistical Modeling with Longitudinal Data: An Application to Ecological Momentary Assessment of Behavior Change in Individuals with Alcohol Use Disorder.

Ecological momentary assessment (EMA) has been broadly used to collect real-time longitudinal data in behavioral research. Several analytic methods have been applied to EMA data to understand the changes of motivation, behavior, and emotions on a daily or within-day basis. One challenge when utilizing those methods on intensive datasets in the behavioral field is to understand when and why the methods are appropriate to investigate particular research questions.

Beyond body size-new traits for new heights in trait-based modelling of predator-prey dynamics.

Food webs map feeding interactions among species, providing a valuable tool for understanding and predicting community dynamics. Using species' body sizes is a promising avenue for parameterizing food-web models, but such approaches have not yet been able to fully recover observed community dynamics. Such discrepancies suggest that traits other than body size also play important roles.

Simultaneous Measurement of Striatal Dopamine and Hydrogen Peroxide Transients Associated with L-DOPA Induced Rotation in Hemiparkinsonian Rats.

Parkinson's disease (PD) is a neurodegenerative disorder commonly treated with levodopa (L-DOPA), which eventually induces abnormal involuntary movements (AIMs). The neurochemical contributors to these dyskinesias are unknown; however, several lines of evidence indicate an interplay of dopamine (DA) and oxidative stress. Here, DA and hydrogen peroxide (HO) were simultaneously monitored at discrete recording sites in the dorsal striata of hemiparkinsonian rats using fast-scan cyclic voltammetry.

Modeling the synergistic properties of drugs in hormonal treatment for prostate cancer.

Prostate cancer is one of the most prevalent cancers in men, with increasing incidence worldwide. This public health concern has inspired considerable effort to study various aspects of prostate cancer treatment using dynamical models, especially in clinical settings. The standard of care for metastatic prostate cancer is hormonal therapy, which reduces the production of androgen that fuels the growth of prostate tumor cells prior to treatment resistance.

Correctly modeling plant-insect-herbivore-pesticide interactions as aggregate data.

We consider a population dynamics model in investigating data from controlled experiments with aphids in broccoli patches surrounded by different margin types (bare or weedy ground) and three levels of insecticide spray (no, light, or heavy spray). The experimental data is clearly aggregate in nature. In previous efforts [1], the aggregate nature of the data was ignored.

Learning partial differential equations for biological transport models from noisy spatio-temporal data.

We investigate methods for learning partial differential equation (PDE) models from spatio-temporal data under biologically realistic levels and forms of noise. Recent progress in learning PDEs from data have used sparse regression to select candidate terms from a denoised set of data, including approximated partial derivatives. We analyse the performance in using previous methods to denoise data for the task of discovering the governing system of PDEs.

Estimation of probability distributions of parameters using aggregate population data: analysis of a CAR T-cell cancer model.

In this effort we explain fundamental formulations for aggregate data inverse problems requiring estimation of probability distribution parameters. We use as a motivating example a class of CAR T-call cancer models in mice. After ascertaining results on model stability and sensitivity with respect to parameters, we carry out first elementary computations on the question how much data is needed for successful estimation of probability distributions.

Modeling Zika Virus Transmission Dynamics: Parameter Estimates, Disease Characteristics, and Prevention.

Because of limited data, much remains uncertain about parameters related to transmission dynamics of Zika virus (ZIKV). Estimating a large number of parameters from the limited information in data may not provide useful knowledge about the ZIKV. Here, we developed a method that utilizes a mathematical model of ZIKV dynamics and the complex-step derivative approximation technique to identify parameters that can be estimated from the available data.

Population model for the decline of Homalodisca vitripennis (Hemiptera: Cicadellidae) over a ten-year period.

The glassy-winged sharpshooter, Homalodisca vitripennis (Germar), is an invasive pest which presents a major economic threat to grape industries in California, because it spreads a disease-causing bacterium, Xylella fastidiosa. In this note we develop a time and temperature dependent mathematical model to analyze aggregate population data for H. vitripennis from a 10-year study consisting of biweekly monitoring of H.

Analysis of compartments-in-series models of liver metabolism as partial differential equations: the effect of dispersion and number of compartments.

Non-alcoholic fatty liver disease is the most common cause of chronic liver disease. Precipitated by the build up of extra fat in the liver not caused by alcohol, it is still not understood why steatosis occurs where it does in the liver microstructure in non-alcoholic fatty liver disease. It is likely, however, that the location of steatosis is due, at least in part, to metabolic zonation (heterogeneity among liver cells in function and enzyme expression).

Optimal design for dynamical modeling of pest populations.

We apply SE-optimal design methodology to investigate optimal data collection procedures as a first step in investigating information content in ecoinformatics data sets. To illustrate ideas we use a simple phenomenological citrus red mite population model for pest dynamics. First the optimal sampling distributions for a varying number of data points are determined.

The effect of pore size and density on ultrasonic attenuation in porous structures with mono-disperse random pore distribution: A two-dimensional in-silico study.

This work proposes a power law model to describe the attenuation of ultrasonic waves in non-absorbing heterogeneous media with randomly distributed scatterers, mimicking a simplified structure of cortical bone. This paper models the propagation in heterogeneous structures with controlled porosity using a two-dimensional finite-difference time domain numerical simulation in order to measure the frequency dependent attenuation. The paper then fits a phenomenological model to the simulated frequency dependent attenuation by optimizing parameters under an ordinary least squares framework.

Ecosystem function in predator-prey food webs-confronting dynamic models with empirical data.

Most ecosystem functions and related services involve species interactions across trophic levels, for example, pollination and biological pest control. Despite this, our understanding of ecosystem function in multitrophic communities is poor, and research has been limited to either manipulation in small communities or statistical descriptions in larger ones. Recent advances in food web ecology may allow us to overcome the trade-off between mechanistic insight and ecological realism.

Estimating intratumoral heterogeneity from spatiotemporal data.

Glioblastoma multiforme (GBM) is a malignant brain cancer with a tendency to both migrate and proliferate. We propose modeling GBM with heterogeneity in cell phenotypes using a random differential equation version of the reaction-diffusion equation, where the parameters describing diffusion (D) and proliferation ([Formula: see text]) are random variables. We investigate the ability to perform the inverse problem to recover the probability distributions of D and [Formula: see text] using the Prohorov metric, for a variety of probability distribution functions.

From theory to experimental design-Quantifying a trait-based theory of predator-prey dynamics.

Successfully applying theoretical models to natural communities and predicting ecosystem behavior under changing conditions is the backbone of predictive ecology. However, the experiments required to test these models are dictated by practical constraints, and models are often opportunistically validated against data for which they were never intended. Alternatively, we can inform and improve experimental design by an in-depth pre-experimental analysis of the model, generating experiments better targeted at testing the validity of a theory.

Forecasting and Uncertainty Quantification Using a Hybrid of Mechanistic and Non-mechanistic Models for an Age-Structured Population Model.

In this paper, we present a new method for the prediction and uncertainty quantification of data-driven multivariate systems. Traditionally, either mechanistic or non-mechanistic modeling methodologies have been used for prediction; however, it is uncommon for the two to be incorporated together. We compare the forecast accuracy of mechanistic modeling, using Bayesian inference, a non-mechanistic modeling approach based on state space reconstruction, and a novel hybrid methodology composed of the two for an age-structured population data set.