Pharmacometrics and Machine Learning Partner to Advance Clinical Data Analysis

Clinical pharmacology is a multidisciplinary data sciences field that utilizes mathematical and statistical methods to generate maximal knowledge from data. Pharmacometrics (PMX) is a well‐recognized tool to characterize disease progression, pharmacokinetics, and risk factors. Because the amount of data produced keeps growing with increasing pace, the computational effort necessary for PMX models is also increasing. Additionally, computationally efficient methods, such as machine learning (ML) are becoming increasingly important in medicine. However, ML is currently not an integrated part of PMX, for various reasons. The goals of this article are to (i) provide an introduction to ML classification methods, (ii) provide examples for a ML classification analysis to identify covariates based on specific research questions, (iii) examine a clinically relevant example to investigate possible relationships of ML and PMX, and (iv) present a summary of ML and PMX tasks to develop clinical decision support tools.ISSN:0009-9236ISSN:1532-653

Leveraging Predictive Pharmacometrics-Based Algorithms to Enhance Perinatal Care—Application to Neonatal Jaundice

The field of medicine is undergoing a fundamental change, transforming towards a modern data-driven patient-oriented approach. This paradigm shift also affects perinatal medicine as predictive algorithms and artificial intelligence are applied to enhance and individualize maternal, neonatal and perinatal care. Here, we introduce a pharmacometrics-based mathematical-statistical computer program (PMX-based algorithm) focusing on hyperbilirubinemia, a medical condition affecting half of all newborns. Independent datasets from two different centers consisting of total serum bilirubin measurements were utilized for model development (342 neonates, 1,478 bilirubin measurements) and validation (1,101 neonates, 3,081 bilirubin measurements), respectively. The mathematical-statistical structure of the PMX-based algorithm is a differential equation in the context of non-linear mixed effects modeling, together with Empirical Bayesian Estimation to predict bilirubin kinetics for a new patient. Several clinically relevant prediction scenarios were validated, i.e., prediction up to 24 h based on one bilirubin measurement, and prediction up to 48 h based on two bilirubin measurements. The PMX-based algorithm can be applied in two different clinical scenarios. First, bilirubin kinetics can be predicted up to 24 h based on one single bilirubin measurement with a median relative (absolute) prediction difference of 8.5% (median absolute prediction difference 17.4 mu mol/l), and sensitivity and specificity of 95.7 and 96.3%, respectively. Second, bilirubin kinetics can be predicted up to 48 h based on two bilirubin measurements with a median relative (absolute) prediction difference of 9.2% (median absolute prediction difference 21.5 mu mol/l), and sensitivity and specificity of 93.0 and 92.1%, respectively. In contrast to currently available nomogram-based static bilirubin stratification, the PMX-based algorithm presented here is a dynamic approach predicting individual bilirubin kinetics up to 48 h, an intelligent, predictive algorithm that can be incorporated in a clinical decision support tool. Such clinical decision support tools have the potential to benefit perinatal medicine facilitating personalized care of mothers and their born and unborn infants

Contribution of machine learning to tumor growth inhibition modeling for hepatocellular carcinoma patients under Roblitinib (FGF401) drug treatment.

Machine learning (ML) opens new perspectives in identifying predictive factors of efficacy among a large number of patients' characteristics in oncology studies. The objective of this work was to combine ML with population pharmacokinetic/pharmacodynamic (PK/PD) modeling of tumor growth inhibition to understand the sources of variability between patients and therefore improve model predictions to support drug development decisions. Data from 127 patients with hepatocellular carcinoma enrolled in a phase I/II study evaluating once-daily oral doses of the fibroblast growth factor receptor FGFR4 kinase inhibitor, Roblitinib (FGF401), were used. Roblitinib  PKs was best described by a two-compartment model with a delayed zero-order absorption and linear elimination. Clinical efficacy using the longitudinal sum of the longest lesion diameter data was described with a population PK/PD model of tumor growth inhibition including resistance to treatment. ML, applying elastic net modeling of time to progression data, was associated with cross-validation, and allowed to derive a composite predictive risk score from a set of 75 patients' baseline characteristics. The two approaches were combined by testing the inclusion of the continuous risk score as a covariate on PD model parameters. The score was found as a significant covariate on the resistance parameter and resulted in 19% reduction of its variability, and 32% variability reduction on the average dose for stasis. The final PK/PD model was used to simulate effect of patients' characteristics on tumor growth inhibition profiles. The proposed methodology can be used to support drug development decisions, especially when large interpatient variability is observed

Population Pharmacokinetic Modeling in the Presence of Missing Time-Dependent Covariates: Impact of Body Weight on Pharmacokinetics of Paracetamol in Neonates

Body weight is the primary covariate in pharmacokinetics of many drugs and dramatically changes during the first weeks of life of neonates. The objective of this study is to determine if missing body weights in preterm and term neonates affect estimates of model parameters and which methods can be used to improve performance of a population pharmacokinetic model of paracetamol. Data for our analysis were obtained from previously published studies on the pharmacokinetics of intravenous paracetamol in neonates. We adopted a population model of body weight change in neonates to implement three previously introduced methods of handling missing covariates based on data imputation, likelihood function modification, and full random effects modeling. All models were implemented in NONMEM 7.4, and population parameters were estimated using the FOCE method. Our major finding was that missing body weights minimally affect population estimates of pharmacokinetic parameters but do affect the covariate relationship parameters, particularly the one describing dependence of clearance on body weight. None of the tested methods changed estimates of between-subject variability nor impacted the predictive performance of the model. Our analysis shows that a modeling approach towards handling missing covariates allows borrowing information gathered in various studies as long as they target the same population. This approach is particularly useful for handling time-dependent missing covariates

Population Pharmacokinetic Modeling in the Presence of Missing Time-Dependent Covariates: Impact of Body Weight on Pharmacokinetics of Paracetamol in Neonates

Body weight is the primary covariate in pharmacokinetics of many drugs and dramatically changes during the first weeks of life of neonates. The objective of this study is to determine if missing body weights in preterm and term neonates affect estimates of model parameters and which methods can be used to improve performance of a population pharmacokinetic model of paracetamol. Data for our analysis were obtained from previously published studies on the pharmacokinetics of intravenous paracetamol in neonates. We adopted a population model of body weight change in neonates to implement three previously introduced methods of handling missing covariates based on data imputation, likelihood function modification, and full random effects modeling. All models were implemented in NONMEM 7.4, and population parameters were estimated using the FOCE method. Our major finding was that missing body weights minimally affect population estimates of pharmacokinetic parameters but do affect the covariate relationship parameters, particularly the one describing dependence of clearance on body weight. None of the tested methods changed estimates of between-subject variability nor impacted the predictive performance of the model. Our analysis shows that a modeling approach towards handling missing covariates allows borrowing information gathered in various studies as long as they target the same population. This approach is particularly useful for handling time-dependent missing covariates.status: publishe

Population Pharmacokinetic Modeling in the Presence of Missing Time-Dependent Covariates: Impact of Body Weight on Pharmacokinetics of Paracetamol in Neonates.

Body weight is the primary covariate in pharmacokinetics of many drugs and dramatically changes during the first weeks of life of neonates. The objective of this study is to determine if missing body weights in preterm and term neonates affect estimates of model parameters and which methods can be used to improve performance of a population pharmacokinetic model of paracetamol. Data for our analysis were obtained from previously published studies on the pharmacokinetics of intravenous paracetamol in neonates. We adopted a population model of body weight change in neonates to implement three previously introduced methods of handling missing covariates based on data imputation, likelihood function modification, and full random effects modeling. All models were implemented in NONMEM 7.4, and population parameters were estimated using the FOCE method. Our major finding was that missing body weights minimally affect population estimates of pharmacokinetic parameters but do affect the covariate relationship parameters, particularly the one describing dependence of clearance on body weight. None of the tested methods changed estimates of between-subject variability nor impacted the predictive performance of the model. Our analysis shows that a modeling approach towards handling missing covariates allows borrowing information gathered in various studies as long as they target the same population. This approach is particularly useful for handling time-dependent missing covariates

Population Pharmacokinetic Modeling in the Presence of Missing Time-Dependent Covariates: Impact of Body Weight on Pharmacokinetics of Paracetamol in Neonates

Body weight is the primary covariate in pharmacokinetics of many drugs and dramatically changes during the first weeks of life of neonates. The objective of this study is to determine if missing body weights in preterm and term neonates affect estimates of model parameters and which methods can be used to improve performance of a population pharmacokinetic model of paracetamol. Data for our analysis were obtained from previously published studies on the pharmacokinetics of intravenous paracetamol in neonates. We adopted a population model of body weight change in neonates to implement three previously introduced methods of handling missing covariates based on data imputation, likelihood function modification, and full random effects modeling. All models were implemented in NONMEM 7.4, and population parameters were estimated using the FOCE method. Our major finding was that missing body weights minimally affect population estimates of pharmacokinetic parameters but do affect the covariate relationship parameters, particularly the one describing dependence of clearance on body weight. None of the tested methods changed estimates of between-subject variability nor impacted the predictive performance of the model. Our analysis shows that a modeling approach towards handling missing covariates allows borrowing information gathered in various studies as long as they target the same population. This approach is particularly useful for handling time-dependent missing covariates

Effect of Kidney Function on Drug Kinetics and Dosing in Neonates, Infants, and Children

Neonates, infants, and children differ from adults in many aspects, not just in age, weight, and body composition. Growth, maturation and environmental factors affect drug kinetics, response and dosing in pediatric patients. Almost 80% of drugs have not been studied in children, and dosing of these drugs is derived from adult doses by adjusting for body weight/size. As developmental and maturational changes are complex processes, such simplified methods may result in subtherapeutic effects or adverse events. Kidney function is impaired during the first 2 years of life as a result of normal growth and development. Reduced kidney function during childhood has an impact not only on renal clearance but also on absorption, distribution, metabolism and nonrenal clearance of drugs. 'Omics'-based technologies, such as proteomics and metabolomics, can be leveraged to uncover novel markers for kidney function during normal development, acute kidney injury, and chronic diseases. Pharmacometric modeling and simulation can be applied to simplify the design of pediatric investigations, characterize the effects of kidney function on drug exposure and response, and fine-tune dosing in pediatric patients, especially in those with impaired kidney function. One case study of amikacin dosing in neonates with reduced kidney function is presented. Collaborative efforts between clinicians and scientists in academia, industry, and regulatory agencies are required to evaluate new renal biomarkers, collect and share prospective pharmacokinetic, genetic and clinical data, build integrated pharmacometric models for key drugs, optimize and standardize dosing strategies, develop bedside decision tools, and enhance labels of drugs utilized in neonates, infants, and children