Development of a Predictive Statistical Pharmacological Model for Local Anesthetic Agent Effects with Bayesian Hierarchical Model Parameter Estimation

Abstract

As an alternative to animal use, computer simulations are useful for predicting pharmacokinetics and cardiovascular activities. For this purpose, we constructed a statistical model to simulate the effects of local anesthetic agents. To train the model, animal experiments were performed on 6-week-old male Hartley guinea pigs. Firstly, the guinea pigs’ backs were shaved, then local anesthetic agents were subcutaneously injected, with subsequent stimulation of the anesthetized site with a needle six times at regular intervals. The number of reactions (score value) was counted. In this statistical model, the probability of reacting to needle stimulation was calculated using the elapsed time, type of local anesthetic agent, and presence or absence of adrenaline. Score values were assumed to follow a binomial distribution at the calculated probability. Parameters were estimated using the Bayesian hierarchical model and Hamiltonian Monte Carlo method. The predicted curves using the estimated parameters fitted well the observed animal values. When score values were predicted using randomly generated parameters, the median of duration was similar between animal experiments and simulations (Procaine: 55 min vs. 50 min, Lidocaine: both 60 min, and Mepivacaine: both 85 min). This approach effectively modeled the effects of local anesthetic agents. It is possible to create the simulator using the parameter values estimated in this study.