The most common approach to in vivo pharmacokinetic and pharmacodynamic modeling involves sequential analysis of the plasma concentration versus time and then response versus time data, such that the plasma kinetic model provides an independent function, driving the dynamics. However, response versus time data, even in the absence of measured drug concentrations, inherently contain useful information about the turnover characteristics of response (turnover rate, half-life of response), the drug's biophase kinetics (F, half-life) as well as the pharmacodynamic characteristics (potency, intrinsic activity). Previous analyses have assumed linear kinetics, linear dynamics, no time lag between kinetics and dynamics (single-valued response), and time constant parameters. However, this report demonstrates that the drug effect can be indirect (antinociception, cortisol/adrenocorticotropin (ACTH), body temperature), display nonlinear kinetics, display feedback mechanisms (nonstationarity, cortiso/ACTH) and exhibit hysteresis with the drug levels in the biophase (antinociception, body temperature). It is also demonstrated that crucial determinants of the success of modeling dose-response-time data are the dose selection, multiple dosing, and to some extent different input rates and routes. This report exemplifies the possibility of assigning kinetic forcing functions in pharmacodynamic modeling in both preclinical and clinical studies for the purpose of characterizing (discrimination between turnover and drug-specific parameters) response data and optimizing subsequent clinical protocols, and for identification of inter-individual differences.