The hypothesis that acetylcholinesterase (AChE) inhibition is the mechanism of toxicity of organophosphorus (OP) compounds was examined by mathematically modeling the in vivo lethal effects of OP compounds and determining the amount of variation in OP toxicity that is explained by AChE inhibition. Mortality dose-response curves for several OP compounds (i.e., VX, soman, cyclosarin, sarin, tabun, diisopropylfluorophosphate and paraoxon) exhibited steep probit slopes (> 9.6) in guinea pigs. Steep probit slopes were also observed when the mortality dose-response curves for soman were examined in mice, rats, rabbits and non-human primates. The consistently steep probit slopes of the dose-response curves for highly toxic OP compounds suggested that these compounds have a single specific mechanism of toxicity regardless of the OP compound or the species in which it was tested. Regression analysis indicated that 93% of the 3,280-fold variation in the median lethal doses (i.e., LD(50)) of OP compounds in rats was explained by the variation in their in vitro rate constants for inhibition of AChE. Conversely, 91% of the 23-fold variation in the ability of the oximes pralidoxime and obidoxime to protect against the toxicity of OP compounds in guinea pigs was explained by the variation in the in vitro ability of oximes to reactivate OP-inhibited AChE. The best explanation for this variety of observations was that the primary mechanism of in vivo toxicity for highly toxic OP compounds is the inhibition of AChE, and the residual unexplained variation in OP toxicity that might be explained by other mechanisms represents < 10% of the total variation in OP toxicity.