Physiological-based pharmacokinetic modeling trends in pharmaceutical drug development over the last 20-years; in-depth analysis of applications, organizations, and platforms.

We assess the advancement of physiologically based pharmacokinetic (PBPK) modeling and simulation (M&S) over the last 20 years (start of 2000 to end of 2019) focusing on the trends in each decade with the relative contributions from different organizations, areas of applications, and software tools used. Unlike many of the previous publications which focused on regulatory applications, our analysis is based on PBPK publications in peer-reviewed journals based on a large sample (>700 original articles). We estimated a rate of growth for PBPK (>40 fold/20 years) that was much steeper than the general pharmacokinetic modeling (<3 fold/20 years) or overall scientific publications (∼3 fold/20 years). The analyses demonstrated that contrary to commonly held belief, commercial specialized PBPK platforms with graphical-user interface were a much more popular choice than open-source alternatives even within academic organizations. These platforms constituted 81% of the whole set of the sample we assessed. The major PBPK applications (top 3) were associated with the study design, predicting formulation effects, and metabolic drug-drug interactions, while studying the fate of drugs in special populations, predicting kinetics in early drug development, and investigating transporter drug interactions have increased proportionally over the last decade. The proportions of application areas based on published research were distinctively different from those shown previously for the regulatory submissions and impact on labels. This may demonstrate the lag time between the research applications versus verified usage within the regulatory framework. The report showed the trend of overall PBPK publications in pharmacology drug development from the past 2 decades stratified by the organizations involved, software used, and area of applications. The analysis showed a more rapid increase in PBPK than that of the pharmacokinetic space itself with an equal contribution from academia and industry. By establishing and recording the journey of PBPK modeling in the past and looking at its current status, the analysis can be used for devising plans based on the anticipated trajectory of future regulatory applications.

Development of a whole body physiologically based model to characterise the pharmacokinetics of benzodiazepines. 1: Estimation of rat tissue-plasma partition ratios.

Three methods for estimation of the equilibrium tissue-to-plasma partition ratios (Kp values) in the presence of tissue concentration time data have been investigated. These are the area method, the open loop (tissue specific) method and the whole body model(closed loop) method, each with different model assumptions. Additionally, multiple imputations, a technique for dealing with deficiencies in data sets (i.e., missing tissues) is used. The estimated Kp values by the three methods have been compared and the limitations and advantages of each approach drawn. The area method, which is essentially model free, gives only a crude estimate of Kp without making any statement of its uncertainty; whereas both the open and closed loop methods provide an estimate of this. The closed loop method, where the most assumptions are made, is the approach that gives the best overall estimates of Kp, which was confirmed by comparing the predicted concentration-time profiles with experimental data. Although the estimates from the closed loop method, as well as the other two methods, are conditioned on the data, they are the most reliable for both propagating parameter variability and uncertainty through a whole body physiologically based model, as well as for extrapolation to human. A series of benzodiazepines, namely alprazolam, chlordiazepoxide, clobazam, diazepam, flunitrazepam, midazolam and triazolam in rat is used as a case study in the current investigation.