Abstract The objective of the study was to develop population pharmacokinetic model of topiramate (TPM) using nonlinear mixed effects modelling approach. Data were collected from 78 adult epileptic patients on mono- or co-therapy of TPM and other antiepileptic drugs, such as carbamazepine (CBZ), valproic acid, lamotrigine, levetiracetam, phenobarbital and pregabalin. Steady-state TPM concentrations were determined in blood samples by high performance liquid chromatography with fluorescence detection. A one-compartment model with first order absorption and elimination was used to fit the concentration-time TPM data. Volume of distribution of TPM was estimated at 0.575l/kg. The influence of demographic, biochemical parameters and therapy characteristics of the patients on oral clearance (CL/F) was evaluated. Daily carbamazepine dose (DCBZ) and renal function estimated by Modification of diet in renal disease (MDRD) formula significantly (p<0.001) influenced CL/F and were included in the final model: CL/F·(l/h)=1.53(l/h)·[1+0.476·DCBZ(mg/day)/1000(mg/day)]·EXP[0.00476·[MDRD(ml/min)-95.72(ml/min)]].. Increase of CL/F with DCBZ and MDRD was best described by linear and exponential models. Mean TPM CL/F during CBZ co-therapy was 2.46l/h, which is higher for 60.8% than in patients not co-treated with CBZ. Evaluation by bootstrapping showed that the final model was stable. The predictive performance was evaluated by adequate plots and indicated satisfactory precision. This model allows individualisation of TPM dosing in routine patient care, especially useful for patients on different CBZ dosing regimen.