Quantum mechanics and mixed quantum mechanics/molecular mechanics simulations of model nerve agents with acetylcholinesterase

The accurate modeling of biological processes presents major computational difficulties owing to the inherent complexity of the macromolecular systems of interest. Simulations of biochemical reactivity tend to require highly computationally intensive quantum mechanical methods, but localized chemical effects tend to depend significantly on properties of the extended biological environment - a regime far more readily examined with lower-level classical empirical models. Mixed quantum/classical techniques are gaining in popularity as a means of bridging these competing requirements. Here we present results comparing two quantum mechanics/molecular mechanics implementations (the SIMOMM technique of Gordon et al. as implemented in GAMESS, and the ONIOM technique of Morokuma et al. found in Gaussian 98) as performed on the enzyme acetylcholinesterase and model nerve agents. This work represents part of the initial phase of a DoD HPCMP Challenge project in which we are attempting to reliably characterize the biochemical processes responsible for nerve agent activity and inhibition, thereby allowing predictions on compounds unrelated to those already studied.