Browse
Submit Data
Databases
API
Help

Dataset Information

17 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Well-tempered MCMC simulations for population pharmacokinetic models.

ABSTRACT: A full Bayesian statistical treatment of complex pharmacokinetic or pharmacodynamic models, in particular in a population context, gives access to powerful inference, including on model structure. Markov Chain Monte Carlo (MCMC) samplers are typically used to estimate the joint posterior parameter distribution of interest. Among MCMC samplers, the simulated tempering algorithm (TMCMC) has a number of advantages: it can sample from sharp multi-modal posteriors; it provides insight into identifiability issues useful for model simplification; it can be used to compute accurate Bayes factors for model choice; the simulated Markov chains mix quickly and have assured convergence in certain conditions. The main challenge when implementing this approach is to find an adequate scale of auxiliary inverse temperatures (perks) and associated scaling constants. We solved that problem by adaptive stochastic optimization and describe our implementation of TMCMC sampling in the GNU MCSim software. Once a grid of perks is obtained, it is easy to perform posterior-tempered MCMC sampling or likelihood-tempered MCMC (thermodynamic integration, which bridges the joint prior and the posterior parameter distributions, with assured convergence of a single sampling chain). We compare TMCMC to other samplers and demonstrate its efficient sampling of multi-modal posteriors and calculation of Bayes factors in two stylized case-studies and two realistic population pharmacokinetic inference problems, one of them involving a large PBPK model.

SUBMITTER: Bois FY

PROVIDER: S-EPMC8082542 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Linking well-tempered metadynamics simulations with experiments.

Project description:Linking experiments with the atomistic resolution provided by molecular dynamics simulations can shed light on the structure and dynamics of protein-disordered states. The sampling limitations of classical molecular dynamics can be overcome using metadynamics, which is based on the introduction of a history-dependent bias on a small number of suitably chosen collective variables. Even if such bias distorts the probability distribution of the other degrees of freedom, the equilibrium Boltzmann distribution can be reconstructed using a recently developed reweighting algorithm. Quantitative comparison with experimental data is thus possible. Here we show the potential of this combined approach by characterizing the conformational ensemble explored by a 13-residue helix-forming peptide by means of a well-tempered metadynamics/parallel tempering approach and comparing the reconstructed nuclear magnetic resonance scalar couplings with experimental data.

| S-EPMC2862196 | biostudies-literature

Computational screening combined with well-tempered metadynamics simulations identifies potential TMPRSS2 inhibitors.

Project description:Type-II transmembrane serine proteases are effective pharmacological targets for host defence against viral entry and in certain cancer cell progressions. These serine proteases cleave viral spike proteins to expose the fusion peptide for cell entry, which is essential to the life cycle of the virus. TMPRSS2 inhibitors can also fight against respiratory viruses that employ them for cell entry. Our study combining virtual screening, all-atom molecular dynamics, and well-tempered metadynamics simulation identifies vicenin-2, neohesperidin, naringin, and rhoifolin as promising TMPRSS2 antagonists. The binding energies obtained are - 16.3, - 15.4, - 13.6, and - 13.8 kcal/mol for vicenin-2, neohesperidin, naringin, and rhoifolin respectively. The RMSD, RMSF, PCA, DCCM, and binding free energy profiles also correlate with the stable binding of these ligands at the active site of TMPRSS2. The study reveals that these molecules could be promising lead molecules for combating future outbreaks of coronavirus and other respiratory viruses.

| S-EPMC11246518 | biostudies-literature

Activation of Abl1 Kinase Explored Using Well-Tempered Metadynamics Simulations on an Essential Dynamics Sampled Path.

Project description:Well-tempered metadynamics (wT-metaD) simulations using path collective variables (CVs) have been successfully applied in recent years to explore conformational transitions in protein kinases and other biomolecular systems. While this methodology has the advantage of describing the transitions with a limited number of predefined path CVs, it requires as an input a reference path connecting the initial and target states of the system. It is desirable to automate the path generation using approaches that do not rely on the choice of geometric CVs to describe the transition of interest. To this end, we developed an approach that couples essential dynamics sampling with wT-metaD simulations. We used this newly developed procedure to explore the activation mechanism of Abl1 kinase and compute the associated free energy barriers. Through these simulations, we identified a three-step mechanism for the activation that involved two metastable intermediates that possessed a partially open activation loop and differed primarily in the "in" or "out" conformation of the aspartate residue of the DFG motif. One of these states is similar to a conformation that was detected in previous spectroscopic studies of Abl1 kinase, albeit its mechanistic role in the activation was hitherto not well understood. The present study establishes its intermediary role in the activation and predicts a rate-determining free energy barrier of 13.8 kcal/mol that is in good agreement with previous experimental and computational estimates. Overall, our study demonstrates the usability of essential dynamics sampling as a path CV in wT-metaD to conveniently study conformational transitions and accurately calculate the associated barriers.

| S-EPMC8582261 | biostudies-literature

Towards the Well-Tempered Chloroplast DNA Sequences.

Project description:With the development of next-generation sequencing technology and bioinformatics tools, the process of assembling DNA sequences has become cheaper and easier, especially in the case of much shorter organelle genomes. The number of available DNA sequences of complete chloroplast genomes in public genetic databases is constantly increasing and the data are widely used in plant phylogenetic and biotechnological research. In this work, we investigated possible inconsistencies in the stored form of publicly available chloroplast genome sequence data. The impact of these inconsistencies on the results of the phylogenetic analysis was investigated and the bioinformatic solution to identify and correct inconsistencies was implemented. The whole procedure was demonstrated using five plant families (Apiaceae, Asteraceae, Campanulaceae, Lamiaceae and Rosaceae) as examples.

| S-EPMC8309291 | biostudies-literature

Well-Tempered Metadynamics Simulations Predict the Structural and Dynamic Properties of a Chiral 24-Atom Macrocycle in Solution.

Project description:Inspired by therapeutic potential, the molecular engineering of macrocycles is garnering increased interest. Exercising control with design, however, is challenging due to the dynamic behavior that these molecules must demonstrate in order to be bioactive. Herein, the value of metadynamics simulations is demonstrated: the free-energy surfaces calculated reveal folded and flattened accessible conformations of a 24-atom macrocycle separated by barriers of ∼6 kT under experimentally relevant conditions. Simulations reveal that the dominant conformer is folded-an observation consistent with a solid-state structure determined by X-ray crystallography and a network of rOes established by 1H NMR. Simulations suggest that the macrocycle exists as a rapidly interconverting pair of enantiomeric, folded structures. Experimentally, 1H NMR shows a single species at room temperature. However, at lower temperature, the interconversion rate between these enantiomers becomes markedly slower, resulting in the decoalescence of enantiotopic methylene protons into diastereotopic, distinguishable resonances due to the persistence of conformational chirality. The emergence of conformational chirality provides critical experimental support for the simulations, revealing the dynamic nature of the scaffold-a trait deemed critical for oral bioactivity.

| S-EPMC9434777 | biostudies-literature

MCMC can detect nonidentifiable models.

Project description:Continuous-time Markov models have been considered the best representation for the stochastic dynamics of ion channels for more than thirty years. For most single-channel data sets, several open and closed states are required for accurately representing the dynamics. However, each data point only shows if the channel is open or closed but not in which state it is. Consequently, some model structures are inherently overparameterized and therefore, in principle, unsuitable for representing any data--those models are called "nonidentifiable". As of this writing, it seems to be poorly understood which continuous-time Markov models are identifiable and which are not, therefore the unconscious use of a nonidentifiable model is a considerable concern. To address this problem, an improved variant of a recently published Markov-chain Monte Carlo method is presented. The algorithm is tested using test data as well as experimental data. We demonstrate that, opposed to a widely used maximum-likelihood estimator, it gives clear warning signs when a nonidentifiable model is used for fitting. Furthermore, for test data that was generated from a nonidentifiable model, the Markov-chain Monte Carlo results recover much more information from the data than maximum-likelihood estimation.

| S-EPMC3514526 | biostudies-literature

Individualized Absorption Models in Population Pharmacokinetic Analyses.

Project description: Not available

| S-EPMC7306615 | biostudies-literature

MCMC estimation of Markov models for ion channels.

Project description:Ion channels are characterized by inherently stochastic behavior which can be represented by continuous-time Markov models (CTMM). Although methods for collecting data from single ion channels are available, translating a time series of open and closed channels to a CTMM remains a challenge. Bayesian statistics combined with Markov chain Monte Carlo (MCMC) sampling provide means for estimating the rate constants of a CTMM directly from single channel data. In this article, different approaches for the MCMC sampling of Markov models are combined. This method, new to our knowledge, detects overparameterizations and gives more accurate results than existing MCMC methods. It shows similar performance as QuB-MIL, which indicates that it also compares well with maximum likelihood estimators. Data collected from an inositol trisphosphate receptor is used to demonstrate how the best model for a given data set can be found in practice.

| S-EPMC3077709 | biostudies-literature

External Evaluation of Two Fluconazole Infant Population Pharmacokinetic Models.

Project description:Fluconazole is an antifungal agent used for the treatment of invasive candidiasis, a leading cause of morbidity and mortality in premature infants. Population pharmacokinetic (PK) models of fluconazole in infants have been previously published by Wade et al. (Antimicrob Agents Chemother 52:4043-4049, 2008, https://doi.org/10.1128/AAC.00569-08) and Momper et al. (Antimicrob Agents Chemother 60:5539-5545, 2016, https://doi.org/10.1128/AAC.00963-16). Here we report the results of the first external evaluation of the predictive performance of both models. We used patient-level data from both studies to externally evaluate both PK models. The predictive performance of each model was evaluated using the model prediction error (PE), mean prediction error (MPE), mean absolute prediction error (MAPE), prediction-corrected visual predictive check (pcVPC), and normalized prediction distribution errors (NPDE). The values of the parameters of each model were reestimated using both the external and merged data sets. When evaluated with the external data set, the model proposed by Wade et al. showed lower median PE, MPE, and MAPE (0.429 ?g/ml, 41.9%, and 57.6%, respectively) than the model proposed by Momper et al. (2.45 ?g/ml, 188%, and 195%, respectively). The values of the majority of reestimated parameters were within 20% of their respective original parameter values for all model evaluations. Our analysis determined that though both models are robust, the model proposed by Wade et al. had greater accuracy and precision than the model proposed by Momper et al., likely because it was derived from a patient population with a wider age range. This study highlights the importance of the external evaluation of infant population PK models.

| S-EPMC5700313 | biostudies-literature

An approach for identifiability of population pharmacokinetic-pharmacodynamic models.

Project description:Mathematical models are routinely used in clinical pharmacology to study the pharmacokinetic and pharmacodynamic properties of a drug in the body. Identifiability of these models is an important requirement for the success of these clinical studies. Identifiability is classified into two types, structural identifiability related to the structure of the mathematical model and deterministic identifiability which is related to the study design. There are existing approaches for assessment of structural identifiability of fixed-effects models, although their use appears uncommon in the literature. In this study, we develop an informal unified approach for simultaneous assessment of structural and deterministic identifiability for fixed and mixed-effects pharmacokinetic or pharmacokinetic-pharmacodynamic models. This approach uses an information theoretic framework. The method is applied both to simple examples to explore known identifiability properties and to a more complex example to illustrate its utility.CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e49; doi:10.1038/psp.2013.25; advance online publication 19 June 2013.

| S-EPMC3697038 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data