Refine
Language
- English (14) (remove)
Is part of the Bibliography
- yes (14)
Keywords
- APC concentration gradient (2)
- Dickkopf diffusion and feedback regulation (2)
- Induced neutropenia (2)
- Model (2)
- Myelosuppression (2)
- Prediction (2)
- Al-26 (1)
- Aluminium (1)
- Continuous renal replacement therapy (1)
- Intensive care (1)
Background
Cells are able to communicate and coordinate their function within tissues via secreted factors. Aberrant secretion by cancer cells can modulate this intercellular communication, in particular in highly organised tissues such as the liver. Hepatocytes, the major cell type of the liver, secrete Dickkopf (Dkk), which inhibits Wnt/ β-catenin signalling in an autocrine and paracrine manner. Consequently, Dkk modulates the expression of Wnt/ β-catenin target genes. We present a mathematical model that describes the autocrine and paracrine regulation of hepatic gene expression by Dkk under wild-type conditions as well as in the presence of mutant cells.
Results
Our spatial model describes the competition of Dkk and Wnt at receptor level, intra-cellular Wnt/ β-catenin signalling, and the regulation of target gene expression for 21 individual hepatocytes. Autocrine and paracrine regulation is mediated through a feedback mechanism via Dkk and Dkk diffusion along the porto-central axis. Along this axis an APC concentration gradient is modelled as experimentally detected in liver. Simulations of mutant cells demonstrate that already a single mutant cell increases overall Dkk concentration. The influence of the mutant cell on gene expression of surrounding wild-type hepatocytes is limited in magnitude and restricted to hepatocytes in close proximity. To explore the underlying molecular mechanisms, we perform a comprehensive analysis of the model parameters such as diffusion coefficient, mutation strength and feedback strength.
Conclusions
Our simulations show that Dkk concentration is elevated in the presence of a mutant cell. However, the impact of these elevated Dkk levels on wild-type hepatocytes is confined in space and magnitude. The combination of inter- and intracellular processes, such as Dkk feedback, diffusion and Wnt/ β-catenin signal transduction, allow wild-type hepatocytes to largely maintain their gene expression.
Background: Cells are able to communicate and coordinate their function within tissues via secreted factors. Aberrant secretion by cancer cells can modulate this intercellular communication, in particular in highly organised tissues such as the liver. Hepatocytes, the major cell type of the liver, secrete Dickkopf (Dkk), which inhibits Wnt/beta-catenin signalling in an autocrine and paracrine manner. Consequently, Dkk modulates the expression of Wnt/beta-catenin target genes. We present a mathematical model that describes the autocrine and paracrine regulation of hepatic gene expression by Dkk under wild-type conditions as well as in the presence of mutant cells. Results: Our spatial model describes the competition of Dkk and Wnt at receptor level, intra-cellular Wnt/beta-catenin signalling, and the regulation of target gene expression for 21 individual hepatocytes. Autocrine and paracrine regulation is mediated through a feedback mechanism via Dkk and Dkk diffusion along the porto-central axis. Along this axis an APC concentration gradient is modelled as experimentally detected in liver. Simulations of mutant cells demonstrate that already a single mutant cell increases overall Dkk concentration. The influence of the mutant cell on gene expression of surrounding wild-type hepatocytes is limited in magnitude and restricted to hepatocytes in close proximity. To explore the underlying molecular mechanisms, we perform a comprehensive analysis of the model parameters such as diffusion coefficient, mutation strength and feedback strength. Conclusions: Our simulations show that Dkk concentration is elevated in the presence of a mutant cell. However, the impact of these elevated Dkk levels on wild-type hepatocytes is confined in space and magnitude. The combination of inter-and intracellular processes, such as Dkk feedback, diffusion and Wnt/beta-catenin signal transduction, allow wild-type hepatocytes to largely maintain their gene expression.
Background Evolution of metastatic melanoma (MM) under B-RAF inhibitors (BRAFi) is unpredictable, but anticipation is crucial for therapeutic decision. Kinetics changes in metastatic growth are driven by molecular and immune events, and thus we hypothesized that they convey relevant information for decision making. Patients and methods We used a retrospective cohort of 37 MM patients treated by BRAFi only with at least 2 close CT-scans available before BRAFi, as a model to study kinetics of metastatic growth before, under and after BRAFi. All metastases (mets) were individually measured at each CT-scan. From these measurements, different measures of growth kinetics of each met and total tumor volume were computed at different time points. A historical cohort permitted to build a reference model for the expected spontaneous disease kinetics without BRAFi. All variables were included in Cox and multistate regression models for survival, to select best candidates for predicting overall survival. Results Before starting BRAFi, fast kinetics and moreover a wide range of kinetics (fast and slow growing mets in a same patient) were pejorative markers. At the first assessment after BRAFi introduction, high heterogeneity of kinetics predicted short survival, and added independent information over RECIST progression in multivariate analysis. Metastatic growth rates after BRAFi discontinuation was usually not faster than before BRAFi introduction, but they were often more heterogeneous than before. Conclusions Monitoring kinetics of different mets before and under BRAFi by repeated CT-scan provides information for predictive mathematical modelling. Disease kinetics deserves more interest
Background: Severe bacterial infections remain a major challenge in intensive care units because of their high prevalence and mortality. Adequate antibiotic exposure has been associated with clinical success in critically ill patients. The objective of this study was to investigate the target attainment of standard meropenem dosing in a heterogeneous critically ill population, to quantify the impact of the full renal function spectrum on meropenem exposure and target attainment, and ultimately to translate the findings into a tool for practical application. Methods: A prospective observational single-centre study was performed with critically ill patients with severe infections receiving standard dosing of meropenem. Serial blood samples were drawn over 4 study days to determine meropenem serum concentrations. Renal function was assessed by creatinine clearance according to the Cockcroft and Gault equation (CLCRCG). Variability in meropenem serum concentrations was quantified at the middle and end of each monitored dosing interval. The attainment of two pharmacokinetic/pharmacodynamic targets (100% T->MIC, 50% T->4xMIC) was evaluated for minimum inhibitory concentration (MIC) values of 2 mg/L and 8 mg/L and standard meropenem dosing (1000 mg, 30-minute infusion, every 8 h). Furthermore, we assessed the impact of CLCRCG on meropenem concentrations and target attainment and developed a tool for risk assessment of target non-attainment. Results: Large inter-and intra-patient variability in meropenem concentrations was observed in the critically ill population (n = 48). Attainment of the target 100% T->MIC was merely 48.4% and 20.6%, given MIC values of 2 mg/L and 8 mg/L, respectively, and similar for the target 50% T->4xMIC. A hyperbolic relationship between CLCRCG (25-255 ml/minute) and meropenem serum concentrations at the end of the dosing interval (C-8h) was derived. For infections with pathogens of MIC 2 mg/L, mild renal impairment up to augmented renal function was identified as a risk factor for target non-attainment (for MIC 8 mg/L, additionally, moderate renal impairment). Conclusions: The investigated standard meropenem dosing regimen appeared to result in insufficient meropenem exposure in a considerable fraction of critically ill patients. An easy-and free-to-use tool (the MeroRisk Calculator) for assessing the risk of target non-attainment for a given renal function and MIC value was developed.
Although the detection of metastases radically changes prognosis of and treatment decisions for a cancer patient, clinically undetectable micrometastases hamper a consistent classification into localized or metastatic disease. This chapter discusses mathematical modeling efforts that could help to estimate the metastatic risk in such a situation. We focus on two approaches: (1) a stochastic framework describing metastatic emission events at random times, formalized via Poisson processes, and (2) a deterministic framework describing the micrometastatic state through a size-structured density function in a partial differential equation model. Three aspects are addressed in this chapter. First, a motivation for the Poisson process framework is presented and modeling hypotheses and mechanisms are introduced. Second, we extend the Poisson model to account for secondary metastatic emission. Third, we highlight an inherent crosslink between the stochastic and deterministic frameworks and discuss its implications. For increased accessibility the chapter is split into an informal presentation of the results using a minimum of mathematical formalism and a rigorous mathematical treatment for more theoretically interested readers.
Background and objective Optimisation of hydrocortisone replacement therapy in children is challenging as there is currently no licensed formulation and dose in Europe for children under 6 years of age. In addition, hydrocortisone has non-linear pharmacokinetics caused by saturable plasma protein binding. A paediatric hydrocortisone formulation, Infacort (R) oral hydrocortisone granules with taste masking, has therefore been developed. The objective of this study was to establish a population pharmacokinetic model based on studies in healthy adult volunteers to predict hydrocortisone exposure in paediatric patients with adrenal insufficiency. Methods Cortisol and binding protein concentrations were evaluated in the absence and presence of dexamethasone in healthy volunteers (n = 30). Dexamethasone was used to suppress endogenous cortisol concentrations prior to and after single doses of 0.5, 2, 5 and 10 mg of Infacort (R) or 20 mg of Infacort (R)/hydrocortisone tablet/hydrocortisone intravenously. A plasma protein binding model was established using unbound and total cortisol concentrations, and sequentially integrated into the pharmacokinetic model. Results Both specific (non-linear) and non-specific (linear) protein binding were included in the cortisol binding model. A two-compartment disposition model with saturable absorption and constant endogenous cortisol baseline (Baseline (cort),15.5 nmol/L) described the data accurately. The predicted cortisol exposure for a given dose varied considerably within a small body weight range in individuals weighing < 20 kg. Conclusions Our semi-mechanistic population pharmacokinetic model for hydrocortisone captures the complex pharmacokinetics of hydrocortisone in a simplified but comprehensive framework. The predicted cortisol exposure indicated the importance of defining an accurate hydrocortisone dose to mimic physiological concentrations for neonates and infants weighing < 20 kg.
Development of a tool to identify intensive care patients at risk of meropenem therapy failure
(2018)
Context Optimization of hydrocortisone replacement therapy is important to prevent under- and over dosing. Hydrocortisone pharmacokinetics is complex as circulating cortisol is protein bound mainly to corticosteroid-binding globulin (CBG) that has a circadian rhythm. Objective A detailed analysis of the CBG circadian rhythm and its impact on cortisol exposure after hydrocortisone administration. Design and Methods CBG was measured over 24 hours in 14 healthy individuals and, employing a modelling and simulation approach using a semi-mechanistic hydrocortisone pharmacokinetic model, we evaluated the impact on cortisol exposure (area under concentration-time curve and maximum concentration of total cortisol) of hydrocortisone administration at different clock times and of the changing CBG concentrations. Results The circadian rhythm of CBG was well described with two cosine terms added to the baseline of CBG: baseline CBG was 21.8 mu g/mL and interindividual variability 11.9%; the amplitude for the 24 and 12 hours cosine functions were relatively small (24 hours: 5.53%, 12 hours: 2.87%) and highest and lowest CBG were measured at 18:00 and 02:00, respectively. In simulations, the lowest cortisol exposure was observed after administration of hydrocortisone at 23:00-02:00, whereas the highest was observed at 15:00-18:00. The differences between the highest and lowest exposure were minor (<= 12.2%), also regarding the free cortisol concentration and free fraction (<= 11.7%). Conclusions Corticosteroid-binding globulin has a circadian rhythm but the difference in cortisol exposure is <= 12.2% between times of highest and lowest CBG concentrations; therefore, hydrocortisone dose adjustment based on time of dosing to adjust for the CBG concentrations is unlikely to be of clinical benefit.
Author summary <br /> The use of orally inhaled drugs for treating lung diseases is appealing since they have the potential for lung selectivity, i.e. high exposure at the site of action -the lung- without excessive side effects. However, the degree of lung selectivity depends on a large number of factors, including physiochemical properties of drug molecules, patient disease state, and inhalation devices. To predict the impact of these factors on drug exposure and thereby to understand the characteristics of an optimal drug for inhalation, we develop a predictive mathematical framework (a "pharmacokinetic model"). In contrast to previous approaches, our model allows combining knowledge from different sources appropriately and its predictions were able to adequately predict different sets of clinical data. Finally, we compare the impact of different factors and find that the most important factors are the size of the inhaled particles, the affinity of the drug to the lung tissue, as well as the rate of drug dissolution in the lung. In contrast to the common belief, the solubility of a drug in the lining fluids is not found to be relevant. These findings are important to understand how inhaled drugs should be designed to achieve best treatment results in patients. <br /> The fate of orally inhaled drugs is determined by pulmonary pharmacokinetic processes such as particle deposition, pulmonary drug dissolution, and mucociliary clearance. Even though each single process has been systematically investigated, a quantitative understanding on the interaction of processes remains limited and therefore identifying optimal drug and formulation characteristics for orally inhaled drugs is still challenging. To investigate this complex interplay, the pulmonary processes can be integrated into mathematical models. However, existing modeling attempts considerably simplify these processes or are not systematically evaluated against (clinical) data. In this work, we developed a mathematical framework based on physiologically-structured population equations to integrate all relevant pulmonary processes mechanistically. A tailored numerical resolution strategy was chosen and the mechanistic model was evaluated systematically against data from different clinical studies. Without adapting the mechanistic model or estimating kinetic parameters based on individual study data, the developed model was able to predict simultaneously (i) lung retention profiles of inhaled insoluble particles, (ii) particle size-dependent pharmacokinetics of inhaled monodisperse particles, (iii) pharmacokinetic differences between inhaled fluticasone propionate and budesonide, as well as (iv) pharmacokinetic differences between healthy volunteers and asthmatic patients. Finally, to identify the most impactful optimization criteria for orally inhaled drugs, the developed mechanistic model was applied to investigate the impact of input parameters on both the pulmonary and systemic exposure. Interestingly, the solubility of the inhaled drug did not have any relevant impact on the local and systemic pharmacokinetics. Instead, the pulmonary dissolution rate, the particle size, the tissue affinity, and the systemic clearance were the most impactful potential optimization parameters. In the future, the developed prediction framework should be considered a powerful tool for identifying optimal drug and formulation characteristics.
Bayesian Data Assimilation to Support Informed Decision Making in Individualized Chemotherapy
(2020)
An essential component of therapeutic drug/biomarker monitoring (TDM) is to combine patient data with prior knowledge for model-based predictions of therapy outcomes. Current Bayesian forecasting tools typically rely only on the most probable model parameters (maximum a posteriori (MAP) estimate). This MAP-based approach, however, does neither necessarily predict the most probable outcome nor does it quantify the risks of treatment inefficacy or toxicity. Bayesian data assimilation (DA) methods overcome these limitations by providing a comprehensive uncertainty quantification. We compare DA methods with MAP-based approaches and show how probabilistic statements about key markers related to chemotherapy-induced neutropenia can be leveraged for more informative decision support in individualized chemotherapy. Sequential Bayesian DA proved to be most computationally efficient for handling interoccasion variability and integrating TDM data. For new digital monitoring devices enabling more frequent data collection, these features will be of critical importance to improve patient care decisions in various therapeutic areas.