Population models for the pharmacokinetic-pharmacodynamic relationship for cetrorelix (CET), a luteinising hormone-releasing hormone (LH-RH) antagonist, and the pharmacodynamic response on testosterone production were investigated in rats and dogs. The plasma concentrations of CET and testosterone were determined after intravenous and subcutaneous injections. The population PK/PD-models were developed using P-PHARM software. Absolute bioavailability of cetrorelix was 100% in rats and 97% in dogs. In rats, the pharmacokinetics was explained by a two-compartment model with saturable absorption, while a three-compartment model was used in dogs. Testosterone suppression in both species was described by a sigmoid E(max) model with maximum effect (E(max)) considered as total hormonal suppression. The duration of testosterone suppression in rats was longer at higher doses. The population elimination half-lives after iv-dose were 3.0 h in rats and 9.3 h in dogs. Population mean estimates of IC50 were 1.39 and 1.24 ng/ml in rats and dogs, respectively. A population pharmacokinetic model was developed to explain the dissolution rate limited absorption from the injection site. The suppression of testosterone could be described by an indirect inhibitory sigmoid E(max) model. In both species 1-2 ng/ml CET in plasma was necessary to suppress testosterone production.