This study developed a population pharmacokinetic model for sodium tanshinone IIA sulfonate (STS) in healthy volunteers and coronary heart disease (CHD) patients in order to identify significant covariates for the pharmacokinetics of STS. Blood samples were obtained by intense sampling approach from 10 healthy volunteers and sparse sampling from 25 CHD patients, and a population pharmacokinetic analysis was performed by nonlinear mixed-effect modeling. The final model was evaluated by bootstrap and visual predictive check. A total of 230 plasma concentrations were included, 137 from healthy volunteers and 93 from CHD patients. It was a two-compartment model with first-order elimination. The typical value of the apparent clearance (CL) of STS in CHD patients with total bilirubin (TBIL) level of 10 µmol·L –1 was 48.7 L·h –1 with inter individual variability of 27.4%, whereas that in healthy volunteers with the same TBIL level was 63.1 L·h –1. Residual variability was described by a proportional error model and estimated at 5.2%. The CL of STS in CHD patients was lower than that in healthy volunteers and decreased when TBIL levels increased. The bootstrap and visual predictive check confirmed the stability and validity of the final model. These results suggested that STS dosage adjustment might be considered based on TBIL levels in CHD patients.