Aim To develop a method to estimate population pharmacokinetic parameters with the limited sampling time points provided clinically during therapeutic drug monitoring. Methods Various simulations were attempted using a one-compartment open model with the first order absorption to determine PK parameter estimates with different sampling strategies as a validation of the method. The estimated parameters were further verified by comparing to the observed values. Results The samples collected at the single time point close to the non-informative sampling time point designed by this method led to bias and inaccurate parameter estimations. Furthermore, the relationship between the estimated non-informative sampling time points and the values of the parameter was examined. The non-informative sampling time points have been developed under some typical occasions and the results were plotted to show the tendency. As a result, one non-informative time point was demonstrated to be appropriate for clearance and two for both volume of distribution and constant of absorption in the present study. It was found that the estimates of the non-informative sampling time points developed in the method increase with increases of volume of distribution and the decrease of clearance and constant of absorption. Conclusion A rational sampling strategy during therapeutic drug monitoring can be established using the method present in the study.
