Dataset Information

Importance of Stability Analysis When Using Nonlinear Semimechanistic Models to Describe Drug-Induced Hematotoxicity.

ABSTRACT: 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.

SUBMITTER: Fornari C

PROVIDER: S-EPMC7499189 | biostudies-literature | 2020 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Importance of Stability Analysis When Using Nonlinear Semimechanistic Models to Describe Drug-Induced Hematotoxicity.

Fornari Chiara C Pin Carmen C Yates James W T JWT Mettetal Jerome T JT Collins Teresa A TA

CPT: pharmacometrics & systems pharmacology 20200708 9

Stability analysis, often overlooked in pharmacometrics, is essential to explore dynamical systems. The model developed by Friberg et al.<sup>1</sup> 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 instabi ...[more]

PMID: 32453487

Similar Datasets

Project description:Taking parameter uncertainty into account is key to make drug development decisions such as testing whether trial endpoints meet defined criteria. Currently used methods for assessing parameter uncertainty in NLMEM have limitations, and there is a lack of diagnostics for when these limitations occur. In this work, a method based on sampling importance resampling (SIR) is proposed, which has the advantage of being free of distributional assumptions and does not require repeated parameter estimation. To perform SIR, a high number of parameter vectors are simulated from a given proposal uncertainty distribution. Their likelihood given the true uncertainty is then approximated by the ratio between the likelihood of the data given each vector and the likelihood of each vector given the proposal distribution, called the importance ratio. Non-parametric uncertainty distributions are obtained by resampling parameter vectors according to probabilities proportional to their importance ratios. Two simulation examples and three real data examples were used to define how SIR should be performed with NLMEM and to investigate the performance of the method. The simulation examples showed that SIR was able to recover the true parameter uncertainty. The real data examples showed that parameter 95 % confidence intervals (CI) obtained with SIR, the covariance matrix, bootstrap and log-likelihood profiling were generally in agreement when 95 % CI were symmetric. For parameters showing asymmetric 95 % CI, SIR 95 % CI provided a close agreement with log-likelihood profiling but often differed from bootstrap 95 % CI which had been shown to be suboptimal for the chosen examples. This work also provides guidance towards the SIR workflow, i.e.,which proposal distribution to choose and how many parameter vectors to sample when performing SIR, using diagnostics developed for this purpose. SIR is a promising approach for assessing parameter uncertainty as it is applicable in many situations where other methods for assessing parameter uncertainty fail, such as in the presence of small datasets, highly nonlinear models or meta-analysis.

Project description:Levofloxacin (LEV) is a broad-spectrum fluoroquinolone used to treat pneumonia, urinary tract infections, chronic bacterial bronchitis, and prostatitis. Efflux transporters, primarily P-glycoprotein (P-gp), are involved in LEV's tissue penetration. In the present work, LEV free lung and prostate interstitial space fluid (ISF) concentrations were evaluated by microdialysis in Wistar rats after intravenous (i.v.) and intratracheal (i.t.) administration (7 mg/kg of body weight) with and without coadministration of the P-gp inhibitor tariquidar (TAR; 15 mg/kg administered i.v.). Plasma and tissue concentration/time profiles were evaluated by noncompartmental analysis (NCA) and population pharmacokinetics (popPK) analysis. The NCA showed significant differences in bioavailability (F) for the control group (0.4) and the TAR group (0.86) after i.t. administration. A four-compartment model simultaneously characterized total plasma and free lung (compartment 2) and prostate (compartment 3) ISF concentrations. Statistically significant differences in lung and prostate average ISF concentrations and levels of kidney active secretion in the TAR group from those measured for the control group (LEV alone) were observed. The estimated population means were as follows: volume of the central compartment (V1), 0.321 liters; total plasma clearance (CL), 0.220 liters/h; TAR plasma clearance (CLTAR), 0.180 liters/h. The intercompartmental distribution rate constants (K values) were as follows: K12, 8.826 h(-1); K21, 7.271 h(-1); K13, 0.047 h(-1); K31, 7.738 h(-1); K14, 0.908 h(-1); K41, 0.409 h(-1); K21 lung TAR (K21LTAR), 8.883 h(-1); K31 prostate TAR (K31PTAR), 4.377 h(-1). The presence of P-gp considerably impacted the active renal secretion of LEV but had only a minor impact on the efflux from the lung following intratracheal dosing. Our results strongly support the idea of a role of efflux transporters other than P-gp contributing to LEV's tissue penetration into the prostrate.

Project description:BackgroundDrug-induced liver injury (DILI) is a critical issue in drug development because DILI causes failures in clinical trials and the withdrawal of approved drugs from the market. There have been many attempts to predict the risk of DILI based on in vivo and in silico identification of hepatotoxic compounds. In the current study, we propose the in silico prediction model predicting DILI using weighted molecular fingerprints.ResultsIn this study, we used 881 bits of molecular fingerprint and used as features describing presence or absence of each substructure of compounds. Then, the Bayesian probability of each substructure was calculated and labeled (positive or negative for DILI), and a weighted fingerprint was determined from the ratio of DILI-positive to DILI-negative probability values. Using weighted fingerprint features, the prediction models were trained and evaluated with the Random Forest (RF) and Support Vector Machine (SVM) algorithms. The constructed models yielded accuracies of 73.8% and 72.6%, AUCs of 0.791 and 0.768 in cross-validation. In independent tests, models achieved accuracies of 60.1% and 61.1% for RF and SVM, respectively. The results validated that weighted features helped increase overall performance of prediction models. The constructed models were further applied to the prediction of natural compounds in herbs to identify DILI potential, and 13,996 unique herbal compounds were predicted as DILI-positive with the SVM model.ConclusionsThe prediction models with weighted features increased the performance compared to non-weighted models. Moreover, we predicted the DILI potential of herbs with the best performed model, and the prediction results suggest that many herbal compounds could have potential to be DILI. We can thus infer that taking natural products without detailed references about the relevant pathways may be dangerous. Considering the frequency of use of compounds in natural herbs and their increased application in drug development, DILI labeling would be very important.

Project description:BACKGROUND: The reconstruction of gene regulatory networks from time series gene expression data is one of the most difficult problems in systems biology. This is due to several reasons, among them the combinatorial explosion of possible network topologies, limited information content of the experimental data with high levels of noise, and the complexity of gene regulation at the transcriptional, translational and post-translational levels. At the same time, quantitative, dynamic models, ideally with probability distributions over model topologies and parameters, are highly desirable. RESULTS: We present a novel approach to infer such models from data, based on nonlinear differential equations, which we embed into a stochastic Bayesian framework. We thus address both the stochasticity of experimental data and the need for quantitative dynamic models. Furthermore, the Bayesian framework allows it to easily integrate prior knowledge into the inference process. Using stochastic sampling from the Bayes' posterior distribution, our approach can infer different likely network topologies and model parameters along with their respective probabilities from given data. We evaluate our approach on simulated data and the challenge #3 data from the DREAM 2 initiative. On the simulated data, we study effects of different levels of noise and dataset sizes. Results on real data show that the dynamics and main regulatory interactions are correctly reconstructed. CONCLUSIONS: Our approach combines dynamic modeling using differential equations with a stochastic learning framework, thus bridging the gap between biophysical modeling and stochastic inference approaches. Results show that the method can reap the advantages of both worlds, and allows the reconstruction of biophysically accurate dynamic models from noisy data. In addition, the stochastic learning framework used permits the computation of probability distributions over models and model parameters, which holds interesting prospects for experimental design purposes.

Project description:Although the mechanisms underlying benzene-induced toxicity and leukemogenicity are not yet fully understood, they are likely to be complicated by various pathways, including those of metabolism, growth factor regulation, oxidative stress, DNA damage, cell cycle regulation, and programmed cell death. With this as a background, we performed cDNA microarray analyses on mouse bone marrow tissue during and after a 2-week benzene exposure by inhalation. Our goal was to clarify the mechanisms underlying the hematotoxicity and leukemogenicity induced by benzene at the level of altered multigene expression. Because a few researchers have postulated that the cell cycle regulation mediated by p53 is a critical event for benzene-induced hematotoxicity, the present study was carried out using p53-knockout (KO) mice and C57BL/6 mice. On the basis of the results of large-scale gene expression studies, we conclude the following: (a) Benzene induces DNA damage in cells at any phase of the cell cycle through myeloperoxidase and in the redox cycle, resulting in p53 expression through Raf-1 and cyclin D-interacting myb-like protein 1. (b) For G1/S cell cycle arrest, the p53-mediated pathway through p21 is involved, as well as the pRb gene-mediated pathway. (c) Alteration of cyclin G1 and Wee-1 kinase genes may be related to the G2/M arrest induced by benzene exposure. (d) DNA repair genes such as Rad50 and Rad51 are markedly downregulated in p53-KO mice. (e) p53-mediated caspase 11 activation, aside from p53-mediated Bax gene induction, may be an important pathway for cellular apoptosis after benzene exposure. Our results strongly suggest that the dysfunction of the p53 gene, possibly caused by strong and repeated genetic and epigenetic effects of benzene on candidate leukemia cells, may induce fatal problems such as those of cell cycle checkpoint, apoptosis, and the DNA repair system, finally resulting in hemopoietic malignancies. Our cDNA microarray data provide valuable information for future investigations of the mechanisms underlying the toxicity and leukemogenicity of benzene.

Dataset Information

Importance of Stability Analysis When Using Nonlinear Semimechanistic Models to Describe Drug-Induced Hematotoxicity.

Publications

Importance of Stability Analysis When Using Nonlinear Semimechanistic Models to Describe Drug-Induced Hematotoxicity.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets