Efficient inference for stochastic differential equation mixed-effects models using correlated particle pseudo-marginal algorithms

Dynamic learning of individual-level suicidal ideation trajectories to enhance mental health care

There has recently been an increase in ongoing patient-report routine outcome monitoring for individuals within clinical care, which has corresponded to increased longitudinal information about an individual. However, many models that are aimed at clinical practice have difficulty fully incorporating this information. This is in part due to the difficulty in dealing with the irregularly time-spaced observations that are common in clinical data. Consequently, we built individual-level continuous-time trajectory models of suicidal ideation for a clinical population (N = 585) with data collected via a digital platform. We demonstrate how such models predict an individual's level and variability of future suicide ideation, with implications for the frequency that individuals may need to be observed. These individual-level predictions provide a more personalised understanding than other predictive methods and have implications for enhanced measurement-based care.

Stochastic Differential Equation-based Mixed Effects Model of the Fluid Volume in the Fasted Stomach in Healthy Adult Human

The rate and extent of drug dissolution and absorption from a solid oral dosage form depend largely on the fluid volume along the gastrointestinal tract. Hence, a model built upon the gastric fluid volume profiles can help to predict drug dissolution and subsequent absorption. To capture the great inter- and intra-individual variability (IAV) of the gastric fluid volume in fasted human, a stochastic differential equation (SDE)-based mixed effects model was developed and compared with the ordinary differential equation (ODE)-based model. Twelve fasted healthy adult subjects were enrolled and had their gastric fluid volume measured before and after consumption of 240 mL of water at pre-determined intervals for up to 2 hours post ingestion. The SDE- and ODE-based mixed effects models were implemented and compared using extended Kalman filter algorithm via NONMEM. The SDE approach greatly improved the goodness of fit compared with the ODE counterpart. The proportional and additive measurement error of the final SDE model decreased from 14.4 to 4.10% and from 17.6 to 4.74 mL, respectively. The SDE-based mixed effects model successfully characterized the gastric volume profiles in the fasted healthy subjects, and provided a robust approximation of the physiological parameters in the very dynamic system. The remarkable IAV could be further separated into system dynamics terms and measurement error terms in the SDE model instead of only empirically attributing IAV to measurement errors in the traditional ODE method. The system dynamics were best captured by the random fluctuations of gastric emptying coefficient Kge.

Identifiability analysis for models of the translation kinetics after mRNA transfection

Mechanistic models are a powerful tool to gain insights into biological processes. The parameters of such models, e.g. kinetic rate constants, usually cannot be measured directly but need to be inferred from experimental data. In this article, we study dynamical models of the translation kinetics after mRNA transfection and analyze their parameter identifiability. That is, whether parameters can be uniquely determined from perfect or realistic data in theory and practice. Previous studies have considered ordinary differential equation (ODE) models of the process, and here we formulate a stochastic differential equation (SDE) model. For both model types, we consider structural identifiability based on the model equations and practical identifiability based on simulated as well as experimental data and find that the SDE model provides better parameter identifiability than the ODE model. Moreover, our analysis shows that even for those parameters of the ODE model that are considered to be identifiable, the obtained estimates are sometimes unreliable. Overall, our study clearly demonstrates the relevance of considering different modeling approaches and that stochastic models can provide more reliable and informative results.

Scalable and flexible inference framework for stochastic dynamic single-cell models

Understanding the inherited nature of how biological processes dynamically change over time and exhibit intra- and inter-individual variability, due to the different responses to environmental stimuli and when interacting with other processes, has been a major focus of systems biology. The rise of single-cell fluorescent microscopy has enabled the study of those phenomena. The analysis of single-cell data with mechanistic models offers an invaluable tool to describe dynamic cellular processes and to rationalise cell-to-cell variability within the population. However, extracting mechanistic information from single-cell data has proven difficult. This requires statistical methods to infer unknown model parameters from dynamic, multi-individual data accounting for heterogeneity caused by both intrinsic (e.g. variations in chemical reactions) and extrinsic (e.g. variability in protein concentrations) noise. Although several inference methods exist, the availability of efficient, general and accessible methods that facilitate modelling of single-cell data, remains lacking. Here we present a scalable and flexible framework for Bayesian inference in state-space mixed-effects single-cell models with stochastic dynamic. Our approach infers model parameters when intrinsic noise is modelled by either exact or approximate stochastic simulators, and when extrinsic noise is modelled by either time-varying, or time-constant parameters that vary between cells. We demonstrate the relevance of our approach by studying how cell-to-cell variation in carbon source utilisation affects heterogeneity in the budding yeast Saccharomyces cerevisiae SNF1 nutrient sensing pathway. We identify hexokinase activity as a source of extrinsic noise and deduce that sugar availability dictates cell-to-cell variability.

Modelling of glucose repression signalling in yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae has a sophisticated signalling system that plays a crucial role in cellular adaptation to changing environments. The SNF1 pathway regulates energy homeostasis upon glucose derepression; hence, it plays an important role in various processes, such as metabolism, cell cycle and autophagy. To unravel its behaviour, SNF1 signalling has been extensively studied. However, the pathway components are strongly interconnected and inconstant; therefore, elucidating its dynamic behaviour based on experimental data only is challenging. To tackle this complexity, systems biology approaches have been successfully employed. This review summarizes the progress, advantages and disadvantages of the available mathematical modelling frameworks covering Boolean, dynamic kinetic, single-cell models, which have been used to study processes and phenomena ranging from crosstalks to sources of cell-to-cell variability in the context of SNF1 signalling. Based on the lessons from existing models, we further discuss how to develop a consensus dynamic mechanistic model of the entire SNF1 pathway that can provide novel insights into the dynamics of nutrient signalling.

Symmetric and Asymmetric Diffusions through Age-Varying Mixed-Species Stand Parameters

(1) Background: This paper deals with unevenly aged, whole-stand models from mixed-effect parameters diffusion processes and Voronoi diagram points of view and concentrates on the mixed-species stands in Lithuania. We focus on the Voronoi diagram of potentially available areas to tree positions as the measure of the competition effect of individual trees and the tree diameter at breast height to relate their evolution through time. (2) Methods: We consider a bivariate hybrid mixed-effect parameters stochastic differential equation for the parameterization of the diameter and available polygon area at age to ensure a proper description of the link between them during the age (time) span of a forest stand. In this study, the Voronoi diagram was used as a mathematical tool for the quantitative characterization of inter-tree competition. (3) Results: The newly derived model considers bivariate correlated observations, tree diameter, and polygon area arising from a particular stand and enables defining equations for calculating diameter, polygon-area, and stand-density predictions and forecasts. (4) Conclusions: From a statistical point of view, the newly developed models produced acceptable statistical measures of predictions and forecasts. All the results were implemented in the Maple computer algebra system.