Bayesian beta nonlinear models with constrained parameters to describe ruminal degradation kinetics.
Ontology highlight
ABSTRACT:
SUBMITTER: Salmeron D
PROVIDER: S-EPMC9225230 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Stability analysis, often overlooked in pharmacometrics, is essential to explore dynamical systems. The model developed by Friberg et al.1 to describe drug-induced hematotoxicity is widely used to support decisions across drug development, and parameter values are often identified from observed blood counts. We use stability analysis to study the parametric dependence of stable and unstable solutions of several Friberg-type models and highlight the risks associated with system instability in the context of nonlinear mixed effects modeling. We emphasize the consequences of unstable solutions on prediction performance by demonstrating nonbiological system behaviors in a real case study of drug-induced thrombocytopenia. Ultimately, we provide simple criteria for identifying parameters associated with stable solutions of Friberg-type models. For instance, in the original Friberg model, we find that stability depends only on the parameter that governs the feedback from peripheral cells to progenitors and provide the exact range of values that results in stable solutions.
| S-EPMC7499189 | biostudies-literature
Project description:Neural mass models (NMMs) are important for helping us interpret observations of brain dynamics. They provide a means to understand data in terms of mechanisms such as synaptic interactions between excitatory and inhibitory neuronal populations. To interpret data using NMMs we need to quantitatively compare the output of NMMs with data, and thereby find parameter values for which the model can produce the observed dynamics. Mapping dynamics to NMM parameter values in this way has the potential to improve our understanding of the brain in health and disease. Though abstract, NMMs still comprise of many parameters that are difficult to constrain a priori. This makes it challenging to explore the dynamics of NMMs and elucidate regions of parameter space in which their dynamics best approximate data. Existing approaches to overcome this challenge use a combination of linearising models, constraining the values they can take and exploring restricted subspaces by fixing the values of many parameters a priori. As such, we have little knowledge of the extent to which different regions of parameter space of NMMs can yield dynamics that approximate data, how nonlinearities in models can affect parameter mapping or how best to quantify similarities between model output and data. These issues need to be addressed in order to fully understand the potential and limitations of NMMs, and to aid the development of new models of brain dynamics in the future. To begin to overcome these issues, we present a global nonlinear approach to recovering parameters of NMMs from data. We use global optimisation to explore all parameters of nonlinear NMMs simultaneously, in a minimally constrained way. We do this using multi-objective optimisation (multi-objective evolutionary algorithm, MOEA) so that multiple data features can be quantified. In particular, we use the weighted horizontal visibility graph (wHVG), which is a flexible framework for quantifying different aspects of time series, by converting them into networks. We study EEG alpha activity recorded during the eyes closed resting state from 20 healthy individuals and demonstrate that the MOEA performs favourably compared to single objective approaches. The addition of the wHVG objective allows us to better constrain the model output, which leads to the recovered parameter values being restricted to smaller regions of parameter space, thus improving the practical identifiability of the model. We then use the MOEA to study differences in the alpha rhythm observed in EEG recorded from 20 people with epilepsy. We find that a small number of parameters can explain this difference and that, counterintuitively, the mean excitatory synaptic gain parameter is reduced in people with epilepsy compared to control. In addition, we propose that the MOEA could be used to mine for the presence of pathological rhythms, and demonstrate the application of this to epileptiform spike-wave discharges.
| S-EPMC10075456 | biostudies-literature
Project description:Amyloid formation is implicated in a number of human diseases, and is thought to proceed via a nucleation-dependent polymerization mechanism. Experimenters often wish to relate changes in amyloid formation kinetics, for example, in response to small molecules to specific mechanistic steps along this pathway. However, fitting kinetic fibril formation data to a complex model including explicit rate constants results in an ill-posed problem with a vast number of potential solutions. The levels of uncertainty remaining in parameters calculated from these models, arising both from experimental noise and high levels of degeneracy or codependency in parameters, is often unclear. Here, we demonstrate that a combination of explicit mathematical models with an approximate Bayesian computation approach can be used to assign the mechanistic effects of modulators on amyloid fibril formation. We show that even when exact rate constants cannot be extracted, parameters derived from these rate constants can be recovered and used to assign mechanistic effects and their relative magnitudes with a great deal of confidence. Furthermore, approximate Bayesian computation provides a robust method for visualizing uncertainty remaining in the model parameters, regardless of its origin. We apply these methods to the problem of heparin-mediated tau polymerization, which displays complex kinetic behavior not amenable to analysis by more traditional methods. Our analysis indicates that the role of heparin cannot be explained by enhancement of nucleation alone, as has been previously proposed. The methods described here are applicable to a wide range of systems, as models can be easily adapted to account for new reactions and reversibility.
| S-EPMC5355488 | biostudies-literature
Project description:The mechanical behavior of tissues at the macroscale is tightly coupled to cellular activity at the microscale. Dermal wound healing is a prominent example of a complex system in which multiscale mechanics regulate restoration of tissue form and function. In cutaneous wound healing, a fibrin matrix is populated by fibroblasts migrating in from a surrounding tissue made mostly out of collagen. Fibroblasts both respond to mechanical cues, such as fiber alignment and stiffness, as well as exert active stresses needed for wound closure. Here, we develop a multiscale model with a two-way coupling between a microscale cell adhesion model and a macroscale tissue mechanics model. Starting from the well-known model of adhesion kinetics proposed by Bell, we extend the formulation to account for nonlinear mechanics of fibrin and collagen and show how this nonlinear response naturally captures stretch-driven mechanosensing. We then embed the new nonlinear adhesion model into a custom finite element implementation of tissue mechanical equilibrium. Strains and stresses at the tissue level are coupled with the solution of the microscale adhesion model at each integration point of the finite element mesh. In addition, solution of the adhesion model is coupled with the active contractile stress of the cell population. The multiscale model successfully captures the mechanical response of biopolymer fibers and gels, contractile stresses generated by fibroblasts, and stress-strain contours observed during wound healing. We anticipate that this framework will not only increase our understanding of how mechanical cues guide cellular behavior in cutaneous wound healing, but will also be helpful in the study of mechanobiology, growth, and remodeling in other tissues.
| S-EPMC8874030 | biostudies-literature
Project description:Polar cell growth is a process that couples the establishment of cell polarity with growth and is extremely important in the growth, development, and reproduction of eukaryotic organisms, such as pollen tube growth during plant fertilization and neuronal axon growth in animals. Pollen tube growth requires dynamic but polarized distribution and activation of a signaling protein named ROP1 to the plasma membrane via three processes: positive feedback and negative feedback regulation of ROP1 activation and its lateral diffusion along the plasma membrane. In this paper, we introduce a mechanistic integro-differential equation (IDE) along with constrained semiparametric regression to quantitatively describe the interplay among these three processes that lead to the polar distribution of active ROP1 at a steady state. Moreover, we introduce a population variability by a constrained nonlinear mixed model. Our analysis of ROP1 activity distributions from multiple pollen tubes revealed that the equilibrium between the positive and negative feedbacks for pollen tubes with similar shapes are remarkably stable, permitting us to infer an inherent quantitative relationship between the positive and negative feedback loops that defines the tip growth of pollen tubes and the polarity of tip growth.
| S-EPMC9175011 | biostudies-literature
Project description:In large-scale applications of undirected graphical models, such as social networks and biological networks, similar patterns occur frequently and give rise to similar parameters. In this situation, it is beneficial to group the parameters for more efficient learning. We show that even when the grouping is unknown, we can infer these parameter groups during learning via a Bayesian approach. We impose a Dirichlet process prior on the parameters. Posterior inference usually involves calculating intractable terms, and we propose two approximation algorithms, namely a Metropolis-Hastings algorithm with auxiliary variables and a Gibbs sampling algorithm with "stripped" Beta approximation (Gibbs_SBA). Simulations show that both algorithms outperform conventional maximum likelihood estimation (MLE). Gibbs_SBA's performance is close to Gibbs sampling with exact likelihood calculation. Models learned with Gibbs_SBA also generalize better than the models learned by MLE on real-world Senate voting data.
| S-EPMC4232941 | biostudies-literature
Project description:The Kalman filter provides a simple and efficient algorithm to compute the posterior distribution for state-space models where both the latent state and measurement models are linear and gaussian. Extensions to the Kalman filter, including the extended and unscented Kalman filters, incorporate linearizations for models where the observation model p(observation|state) is nonlinear. We argue that in many cases, a model for p(state|observation) proves both easier to learn and more accurate for latent state estimation. Approximating p(state|observation) as gaussian leads to a new filtering algorithm, the discriminative Kalman filter (DKF), which can perform well even when p(observation|state) is highly nonlinear and/or nongaussian. The approximation, motivated by the Bernstein-von Mises theorem, improves as the dimensionality of the observations increases. The DKF has computational complexity similar to the Kalman filter, allowing it in some cases to perform much faster than particle filters with similar precision, while better accounting for nonlinear and nongaussian observation models than Kalman-based extensions. When the observation model must be learned from training data prior to filtering, off-the-shelf nonlinear and nonparametric regression techniques can provide a gaussian model for p(observation|state) that cleanly integrates with the DKF. As part of the BrainGate2 clinical trial, we successfully implemented gaussian process regression with the DKF framework in a brain-computer interface to provide real-time, closed-loop cursor control to a person with a complete spinal cord injury. In this letter, we explore the theory underlying the DKF, exhibit some illustrative examples, and outline potential extensions.
| S-EPMC8259355 | biostudies-literature
Project description:Mixed effects Beta regression models based on Beta distributions have been widely used to analyze longitudinal percentage or proportional data ranging between zero and one. However, Beta distributions are not flexible to extreme outliers or excessive events around tail areas, and they do not account for the presence of the boundary values zeros and ones because these values are not in the support of the Beta distributions. To address these issues, we propose a mixed effects model using Beta rectangular distribution and augment it with the probabilities of zero and one. We conduct extensive simulation studies to assess the performance of mixed effects models based on both the Beta and Beta rectangular distributions under various scenarios. The simulation studies suggest that the regression models based on Beta rectangular distributions improve the accuracy of parameter estimates in the presence of outliers and heavy tails. The proposed models are applied to the motivating Neuroprotection Exploratory Trials in Parkinson's Disease (PD) Long-term Study-1 (LS-1 study, n = 1741), developed by The National Institute of Neurological Disorders and Stroke Exploratory Trials in Parkinson's Disease (NINDS NET-PD) network.
| S-EPMC5064841 | biostudies-literature
Project description:Detecting disease-associated genomic outcomes is one of the key steps in precision medicine research. Cutting-edge high-throughput technologies enable researchers to unbiasedly test if genomic outcomes are associated with disease of interest. However, these technologies also include the challenges associated with the analysis of genome-wide data. Two big challenges are (1) how to reduce the effects of technical noise; and (2) how to handle the curse of dimensionality (i.e., number of variables are way larger than the number of samples). To tackle these challenges, we propose a constrained mixture of Bayesian hierarchical models (MBHM) for detecting disease-associated genomic outcomes for data obtained from paired/matched designs. Paired/matched designs can effectively reduce effects of confounding factors. MBHM does not involve multiple testing, hence does not have the problem of the curse of dimensionality. It also could borrow information across genes so that it can be used for whole genome data with small sample sizes.
| S-EPMC5373614 | biostudies-literature
Project description:MotivationThe selection among substitution models of molecular evolution is fundamental for obtaining accurate phylogenetic inferences. At the protein level, evolutionary analyses are traditionally based on empirical substitution models but these models make unrealistic assumptions and are being surpassed by structurally constrained substitution (SCS) models. The SCS models often consider site-dependent evolution, a process that provides realism but complicates their implementation into likelihood functions that are commonly used for substitution model selection.ResultsWe present a method to perform selection among site-dependent SCS models, also among empirical and site-dependent SCS models, based on the approximate Bayesian computation (ABC) approach and its implementation into the computational framework ProteinModelerABC. The framework implements ABC with and without regression adjustments and includes diverse empirical and site-dependent SCS models of protein evolution. Using extensive simulated data, we found that it provides selection among SCS and empirical models with acceptable accuracy. As illustrative examples, we applied the framework to analyze a variety of protein families observing that SCS models fit them better than the corresponding best-fitting empirical substitution models.Availability and implementationProteinModelerABC is freely available from https://github.com/DavidFerreiro/ProteinModelerABC, can run in parallel and includes a graphical user interface. The framework is distributed with detailed documentation and ready-to-use examples.
| S-EPMC10914458 | biostudies-literature