Anesth Analg 2009; 108:168-175
© 2009 International Anesthesia Research Society
doi: 10.1213/ane.0b013e31818de158
ANESTHETIC PHARMACOLOGY
A Comparison of the Molecular Bases for N-Methyl-d-Aspartate-Receptor Inhibition Versus Immobilizing Activities of Volatile Aromatic Anesthetics
Jason C. Sewell, PhD*,
Douglas E. Raines, MD ,
Edmond I. Eger, II, MD ,
Michael J. Laster, DVM , and
John W. Sear, PhD, FFARCS*
From the *Nuffield Department of Anaesthetics, University of Oxford, John Radcliffe Hospital, Headington, UK; Department of Anesthesia and Critical Care, Massachusetts General Hospital, Boston, Massachusetts; Department of Anesthesia, Harvard Medical School, Boston, Massachusetts; and Department of Anesthesia and Perioperative Care, University of California, San Francisco, California.
Address correspondence to John W. Sear, PhD, Nuffield Department of Anesthetics, University of Oxford, The John Radcliffe Hospital, Headington, Oxford OX3 9DU, UK. Address e-mail to john.sear{at}nda.ox.ac.uk.
BACKGROUND: Aromatic anesthetics exhibit a wide range of N-methyl-d-aspartate (NMDA) receptor inhibitory potencies and immobilizing activities. We sought to characterize the molecular basis of NMDA receptor inhibition using comparative molecular field analysis (CoMFA), and compare the results to those from an equivalent model for immobilizing activity.
METHODS: Published potency data for 14 compounds were supplemented with new values for 2 additional agents. The anesthetics were divided into a training set (n = 12) used to formulate the activity models and a test set (n = 4) used to independently assess the models predictive capability. The anesthetic structures were geometry optimized using ab initio quantum mechanics and aligned by field-fit minimization to provide the best correlation between the steric and electrostatic fields of the molecules and one or more lead structures. Orientations that yielded CoMFA models with the greatest predictive capability (assessed by leave-one-out cross-validation) were retained.
RESULTS: The final CoMFA model for the inhibition of NR1/NR2B NMDA receptors explained 99.3% of the variance in the observed activities of the 12 training set agents (F2,9 = 661.5, P < 0.0001). The model effectively predicted inhibitory potency for the training set (cross-validated r2CV = 0.944) and 4 excluded test set compounds (predictive r2Pred = 0.966). The equivalent model for immobility in response to noxious stimuli explained 98.0% of the variance in the observed activities for the training set (F2,9 = 219.2, P < 0.0001) and exhibited adequate predictive capability for both the training set (r2CV = 0.872) and test set (r2Pred = 0.926) agents. Comparison of pharmacophoric maps showed that several key steric and electrostatic regions were common to both activity models, but differences were observed in the relative importance of these key regions with respect to the two aspects of anesthetic activity.
CONCLUSIONS: The similarities in the pharmacophoric maps are consistent with NMDA receptors contributing part of the immobilizing activity of volatile aromatic anesthetics.
|