Vertical Integration of Pharmacogenetics in Population PK/PD Modeling: A Novel Information Theoretic Method
Ontology highlight
ABSTRACT:
SUBMITTER: Knights J
PROVIDER: S-EPMC3600754 | biostudies-literature | 2013 Feb
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:To develop a model synthesis method for parsimoniously modeling gene-environmental interactions (GEI) associated with clinical outcomes and phenotypes. The AMBROSIA model synthesis approach utilizes the k-way interaction information (KWII), an information-theoretic metric capable of identifying variable combinations associated with GEI. For model synthesis, AMBROSIA considers relevance of combinations to the phenotype, it precludes entry of combinations with redundant information, and penalizes for unjustifiable complexity; each step is KWII based. The performance and power of AMBROSIA were evaluated with simulations and Genetic Association Workshop 15 (GAW15) data sets of rheumatoid arthritis (RA). AMBROSIA identified parsimonious models in data sets containing multiple interactions with linkage disequilibrium present. For the GAW15 data set containing 9187 single-nucleotide polymorphisms, the parsimonious AMBROSIA model identified nine RA-associated combinations with power >90%. AMBROSIA was compared with multifactor dimensionality reduction across several diverse models and had satisfactory power. Software source code is available from http://www.cse.buffalo.edu/DBGROUP/bioinformatics/resources.html. AMBROSIA is a promising method for GEI model synthesis.
| S-EPMC3182499 | biostudies-literature
Project description:Studying a gene's regulatory mechanisms is a tedious process that involves identification of candidate regulators by transcription factor (TF) knockout or over-expression experiments, delineation of enhancers by reporter assays, and demonstration of direct TF influence by site mutagenesis, among other approaches. Such experiments are often chosen based on the biologist's intuition, from several testable hypotheses. We pursue the goal of making this process systematic by using ideas from information theory to reason about experiments in gene regulation, in the hope of ultimately enabling rigorous experiment design strategies. For this, we make use of a state-of-the-art mathematical model of gene expression, which provides a way to formalize our current knowledge of cis- as well as trans- regulatory mechanisms of a gene. Ambiguities in such knowledge can be expressed as uncertainties in the model, which we capture formally by building an ensemble of plausible models that fit the existing data and defining a probability distribution over the ensemble. We then characterize the impact of a new experiment on our understanding of the gene's regulation based on how the ensemble of plausible models and its probability distribution changes when challenged with results from that experiment. This allows us to assess the 'value' of the experiment retroactively as the reduction in entropy of the distribution (information gain) resulting from the experiment's results. We fully formalize this novel approach to reasoning about gene regulation experiments and use it to evaluate a variety of perturbation experiments on two developmental genes of D. melanogaster. We also provide objective and 'biologist-friendly' descriptions of the information gained from each such experiment. The rigorously defined information theoretic approaches presented here can be used in the future to formulate systematic strategies for experiment design pertaining to studies of gene regulatory mechanisms.
| S-EPMC6175532 | biostudies-literature
Project description:Tumors are the result of accumulated genomic alterations that cooperate synergistically to produce uncontrollable cell growth. Although identifying recurrent alterations among large collections of tumors provides a way to pinpoint genes that endow a selective advantage in oncogenesis and progression, it fails to address the genetic interactions behind this selection process. A non-random pattern of co-mutated genes is evidence for selective forces acting on tumor cells that harbor combinations of these genetic alterations. Although existing methods have successfully identified mutually exclusive gene sets, no current method can systematically discover more general genetic relationships. We develop Genomic Alteration Modules using Total Correlation (GAMToC), an information theoretic framework that integrates copy number and mutation data to identify gene modules with any non-random pattern of joint alteration. Additionally, we present the Seed-GAMToC procedure, which uncovers the mutational context of any putative cancer gene. The software is publicly available. Applied to glioblastoma multiforme samples, GAMToC results show distinct subsets of co-occurring mutations, suggesting distinct mutational routes to cancer and providing new insight into mutations associated with proneural, proneural/G-CIMP, and classical types of the disease. The results recapitulate known relationships such as mutual exclusive mutations, place these alterations in the context of other mutations, and find more complex relationships such as conditional mutual exclusivity.
| S-EPMC4565985 | biostudies-literature
Project description:A phase 1 clinical trial in healthy male volunteers was conducted with a somatostatin-dopamine chimera (BIM23B065), from which information could be obtained on the concentration-effect relationship of the inhibition of pulsatile endogenous growth hormone and prolactin secretion. Endogenous growth hormone profiles were analyzed using a two-step deconvolution-analysis-informed population pharmacodynamic modeling approach, which was developed for the analyses of pulsatile profiles. Prolactin concentrations were modelled using a population pool model with a circadian component on the prolactin release. During treatment with BIM23B065, growth hormone secretion was significantly reduced (maximal effect [EMAX]?=?- 64.8%) with significant reductions in the pulse frequency in two out of three multiple ascending dose cohorts. A circadian component in prolactin secretion was identified, modelled using a combination of two cosine functions with 24 h and 12 h periods. Dosing of BIM23B065 strongly inhibited (EMAX?=?- 91%) the prolactin release and demonstrated further reduction of prolactin secretion after multiple days of dosing. This study quantified the concentration-effect relationship of BIM23B065 on the release of two pituitary hormones, providing proof of pharmacology of the chimeric actions of BIM23B065.
| S-EPMC7289785 | biostudies-literature
Project description:Apramycin represents a subclass of aminoglycoside antibiotics that has been shown to evade almost all mechanisms of clinically relevant aminoglycoside resistance. Model-informed drug development may facilitate its transition from preclinical to clinical phase. This study explored the potential of pharmacokinetic/pharmacodynamic (PK/PD) modeling to maximize the use of in vitro time-kill and in vivo preclinical data for prediction of a human efficacious dose (HED) for apramycin. PK model parameters of apramycin from four different species (mouse, rat, guinea pig, and dog) were allometrically scaled to humans. A semimechanistic PK/PD model was developed from the rich in vitro data on four Escherichia coli strains and subsequently the sparse in vivo efficacy data on the same strains were integrated. An efficacious human dose was predicted from the PK/PD model and compared with the classical PK/PD index methodology and the aminoglycoside dose similarity. One-compartment models described the PK data and human values for clearance and volume of distribution were predicted to 7.07 L/hour and 26.8 L, respectively. The required fAUC/MIC (area under the unbound drug concentration-time curve over MIC ratio) targets for stasis and 1-log kill in the thigh model were 34.5 and 76.2, respectively. The developed PK/PD model predicted the efficacy data well with strain-specific differences in susceptibility, maximum bacterial load, and resistance development. All three dose prediction approaches supported an apramycin daily dose of 30 mg/kg for a typical adult patient. The results indicate that the mechanistic PK/PD modeling approach can be suitable for HED prediction and serves to efficiently integrate all available efficacy data with potential to improve predictive capacity.
| S-EPMC8048880 | biostudies-literature
Project description:Pasteurella multocida can invade and translocate through endothelial cells and result in vascular-system infection, which can cause severe economic losses in the poultry industry. Antibacterial therapy (especially florfenicol) plays an important part in controlling P. multocida infection. To preserve the effect of florfenicol, in vivo pharmacokinetic/pharmacodynamic (PK/PD) modeling of florfenicol against three P. multocida strains in duck was established. Then, the efficacy of the currently marketed dose, a rational dosage regimen for populations, and the PK/PD cutoff were predicted through Monte Carlo simulations (MCSs). The area under the concentration-time curve from 0 to 24 h/minimum inhibitory concentration (AUC0-24 h /MIC) was the optimal PK/PD parameter. The PK/PD surrogate values of florfenicol against P. multocida were similar using different organs as the PD target, but varied in different strains. For the florfenicol-sensitive strain 0825Y1, when the AUC0-24 h /MIC reached 117.54 and 108.19, florfenicol showed a bactericidal effect in the liver and lung, respectively. For the florfenicol-sensitive strain 0901J1, the corresponding value was 78.39 and 54.30, respectively. For the florfenicol-resistant strain JY160110, florfenicol could attain a maximum effect of 1 - log10 reduction in bacteria in the liver and lung when the AUC0-24 h /MIC reached 2.03 and 2.06, respectively. The PK/PD-based prediction for the population dose indicated a poor effect for the low end of the currently marketed dose (40 mg/kg body weight per day), but a robust effect for the high end of the currently marketed dose (60 mg/kg body weight per day) with a target attainment rate of 92.79% and 81.44% against P. multocida in mainland China and worldwide, respectively. The recommended dose optimized by MCSs was 52 mg/kg body weight in mainland China. The PK/PD cutoff of florfenicol against P. multocida at the low end and high end of the current daily dose (40 and 60 mg/kg body weight) and predicted daily dose in mainland China (52 mg/kg body weight) was 0.25, 4, and 0.5 μg/ml, respectively. These results suggested that more than one strain should be involved for PK/PD modeling and contributed to rational use of florfenicol in populations. We also provided fundamental data for determination of florfenicol breakpoints in poultry.
| S-EPMC7829356 | biostudies-literature
Project description:Influenza pandemics and seasonal outbreaks have shown the potential of Influenza A virus (IAV) to enhance susceptibility to a secondary infection with the bacterial pathogen Streptococcus pneumoniae (Sp). The high morbidity and mortality rate revealed the poor efficacy of antiviral drugs and vaccines to fight IAV infections. Currently, the most effective treatment for IAV is by antiviral neuraminidase inhibitors. Among them, the most frequently stockpiled is Oseltamivir which reduces viral release and transmission. However, effectiveness of Oseltamivir is compromised by the emergence of resistant IAV strains and secondary bacterial infections. To date, little attention has been given to evaluate how Oseltamivir treatment strategies alter Influenza viral infection in presence of Sp coinfection and a resistant IAV strain emergence. In this paper we investigate the efficacy of current approved Oseltamivir treatment regimens using a computational approach. Our numerical results suggest that the curative regimen (75 mg) may yield 47% of antiviral efficacy and 9% of antibacterial efficacy. An increment in dose to 150 mg (pandemic regimen) may increase the antiviral efficacy to 49% and the antibacterial efficacy to 16%. The choice to decrease the intake frequency to once per day is not recommended due to a significant reduction in both antiviral and antibacterial efficacy. We also observe that the treatment duration of 10 days may not provide a clear improvement on the antiviral and antibacterial efficacy compared to 5 days. All together, our in silico study reveals the success and pitfalls of Oseltamivir treatment strategies within IAV-Sp coinfection and calls for testing the validity in clinical trials.
| S-EPMC4906052 | biostudies-literature
Project description:BACKGROUND:DNA methylation is a stable form of epigenetic memory used by cells to control gene expression. Whole genome bisulfite sequencing (WGBS) has emerged as a gold-standard experimental technique for studying DNA methylation by producing high resolution genome-wide methylation profiles. Statistical modeling and analysis is employed to computationally extract and quantify information from these profiles in an effort to identify regions of the genome that demonstrate crucial or aberrant epigenetic behavior. However, the performance of most currently available methods for methylation analysis is hampered by their inability to directly account for statistical dependencies between neighboring methylation sites, thus ignoring significant information available in WGBS reads. RESULTS:We present a powerful information-theoretic approach for genome-wide modeling and analysis of WGBS data based on the 1D Ising model of statistical physics. This approach takes into account correlations in methylation by utilizing a joint probability model that encapsulates all information available in WGBS methylation reads and produces accurate results even when applied on single WGBS samples with low coverage. Using the Shannon entropy, our approach provides a rigorous quantification of methylation stochasticity in individual WGBS samples genome-wide. Furthermore, it utilizes the Jensen-Shannon distance to evaluate differences in methylation distributions between a test and a reference sample. Differential performance assessment using simulated and real human lung normal/cancer data demonstrate a clear superiority of our approach over DSS, a recently proposed method for WGBS data analysis. Critically, these results demonstrate that marginal methods become statistically invalid when correlations are present in the data. CONCLUSIONS:This contribution demonstrates clear benefits and the necessity of modeling joint probability distributions of methylation using the 1D Ising model of statistical physics and of quantifying methylation stochasticity using concepts from information theory. By employing this methodology, substantial improvement of DNA methylation analysis can be achieved by effectively taking into account the massive amount of statistical information available in WGBS data, which is largely ignored by existing methods.
| S-EPMC5842653 | biostudies-literature
Project description:Akaike's information theoretic criterion for model discrimination (AIC) is often stated to "overfit", i.e., it selects models with a higher dimension than the dimension of the model that generated the data. However, with experimental pharmacokinetic data it may not be possible to identify the correct model, because of the complexity of the processes governing drug disposition. Instead of trying to find the correct model, a more useful objective might be to minimize the prediction error of drug concentrations in subjects with unknown disposition characteristics. In that case, the AIC might be the selection criterion of choice. We performed Monte Carlo simulations using a model of pharmacokinetic data (a power function of time) with the property that fits with common multi-exponential models can never be perfect - thus resembling the situation with real data. Prespecified models were fitted to simulated data sets, and AIC and AIC c (the criterion with a correction for small sample sizes) values were calculated and averaged. The average predictive performances of the models, quantified using simulated validation sets, were compared to the means of the AICs. The data for fits and validation consisted of 11 concentration measurements each obtained in 5 individuals, with three degrees of interindividual variability in the pharmacokinetic volume of distribution. Mean AIC c corresponded very well, and better than mean AIC, with mean predictive performance. With increasing interindividual variability, there was a trend towards larger optimal models, but with respect to both lowest AIC c and best predictive performance. Furthermore, it was observed that the mean square prediction error itself became less suitable as a validation criterion, and that a predictive performance measure should incorporate interindividual variability. This simulation study showed that, at least in a relatively simple mixed effects modelling context with a set of prespecified models, minimal mean AIC c corresponded to best predictive performance even in the presence of relatively large interindividual variability.
| S-EPMC4670010 | biostudies-literature
Project description:Pharmacological administration of fibroblast growth factor 21 (FGF21) improves metabolic profile in preclinical species and humans. FGF21 exerts its metabolic effects through formation of beta-klotho (KLB)/FGF receptor 1c FGFR1c complex and subsequent signaling. Data from various in vitro systems demonstrate the intact C- and N-terminus of FGF21 is required for binding with KLB, and interaction with FGFR1c, respectively. However the relative roles of the termini for in vivo pharmacological effects are unclear. Here we report PF-05231023, a long-acting FGF21 analogue which is unique in that the half-life and subcutaneous (s.c.) bioavailability of the intact C-terminus are significantly different from those of the intact N-terminus (2 vs. 22 hr for half-life and 4~7 vs. ~50% SC bioavailability). Therefore, this molecule serves as a valuable tool to evaluate the relative roles of intact C-terminus vs. N-terminus in in vivo pharmacology studies in preclinical species. We determined the effects of PF-05231023 administration on body weight (BW) loss and glucose reduction during an oral glucose tolerance test (OGTT) following SC and intravenous (i.v.) administration in diet-induced obese (DIO) and leptin-deficient obese (ob/ob) mice, respectively. Our data show that the intact N-terminus of FGF21 in PF-05231023 appears to be sufficient to drive glucose lowering during OGTT and sustain BW loss in DIOs. Further, PK/PD modeling suggests that while the intact FGF21 C-terminus is not strictly required for glucose lowering during OGTT in ob/ob mice or for BW reduction in DIO mice, the higher potency conferred by intact C-terminus contributes to a rapid initiation of pharmacodynamic effects immediately following dosing. These results provide additional insight into the strategy of developing stabilized versions of FGF21 analogs to harness the full spectrum of its metabolic benefits.
| S-EPMC4366384 | biostudies-literature