Background and Aim: Tacrolimus (TAC) is a first-line immunosuppressant for the treatment of refractory nephrotic syndrome (RNS), but the pharmacokinetics of TAC varies widely among individuals, and there is still no accurate model to predict the pharmacokinetics of TAC in RNS. Therefore, this study aimed to combine population pharmacokinetic (PPK) model and machine learning algorithms to develop a simple and accurate prediction model for TAC.Methods: 139 children with RNS from August 2013 to December 2018 were included, and blood samples of TAC trough and partial peak concentrations were collected. The blood concentration of TAC was determined by enzyme immunoassay; CYP3A5 was genotyped by polymerase chain reaction-restriction fragment length polymorphism method; MYH9, LAMB2, ACTN4 and other genotypes were determined by MALDI-TOF MS method; PPK model was established by nonlinear mixed-effects method. Based on this, six machine learning algorithms, including eXtreme Gradient Boosting (XGBoost), Random Forest (RF), Extra-Trees, Gradient Boosting Decision Tree (GBDT), Adaptive boosting (AdaBoost) and Lasso, were used to establish the machine learning model of TAC clearance.Results: A one-compartment model of first-order absorption and elimination adequately described the pharmacokinetics of TAC. Age, co-administration of Wuzhi capsules, CYP3A5 *3/*3 genotype and CTLA4 rs4553808 genotype were significantly affecting the clearance of TAC. Among the six machine learning models, the Lasso algorithm model performed the best (R2 = 0.42).Conclusion: For the first time, a clearance prediction model of TAC in pediatric patients with RNS was established using PPK combined with machine learning, by which the individual clearance of TAC can be predicted more accurately, and the initial dose of administration can be optimized to achieve the goal of individualized treatment.