Therapeutic drug monitoring (TDM) is gaining importance for improving the success of antiretroviral treatment in human immunodeficiency virus-infected patients. However, enfuvirtide (ENF) concentrations are not regularly determined. The objective of this work was to study the pharmacokinetics (PK) of ENF in patients treated in routine clinical settings, to develop a population PK model describing the concentration-time profile, and to establish PK reference values. A liquid chromatography-tandem mass spectrometry method was developed and applied to serum samples submitted for TDM. A two-compartment model with linear absorption and elimination was fitted to 329 concentrations from 131 patients. The PK model was used for simulations resulting in percentile curves for ENF levels for the full dosing interval. The model predicted that a median concentration of 1,968 ng/ml would be reached 12 h after administration of 90 mg of ENF, and 23% and 58% of patients are expected to have concentrations below 1,000 ng/ml and 2,200 ng/ml, respectively. Both values have been proposed as cutoffs for virological efficacy. The median maximum concentration of drug in serum (Cmax) of 3,943 ng/ml, predicted for 3 h after drug administration, is lower than the Cmax reported previously. We found an enormous interpatient variability at every time point, with concentration spectrums covering >1 log and 52% and 123% interindividual variabilities in the typical clearance and volume of distribution, respectively, in contrast to preexisting PK data. In summary, ENF levels are lower and more variable than expected. Many patients may achieve insufficient concentrations. Further covariate analysis in the population PK model might help to identify factors influencing the variability in ENF concentrations.