This research article by Dr. Willian E.A. de Lima et al. is published in Letters in Drug Design & Discovery Journal Volume 13, Issue 5, 2016
The study of the mechanism of reactivation of acetylcholinesterase (AChE) inhibited by organophosphates (OP) is still a challenge for theoretical chemists since mechanistic studies, involve electronic transfer and breaking and formation of chemical bonds. The solution for this is the Quantum Mechanics/Molecular Mechanics (QM/MM) approach. The QM/MM approach has been implemented with docking, Monte Carlo and molecular dynamics methods. As many configurations are generated in the MM step, the number of QM calculations required in a hybrid QM/MM simulation is too high. In order to circumvent this situation, in this work, an occupancy frequency approach based on chemometric techniques is proposed to select promising representative configurations for reaction mechanism calculations. In fact, this selection method focuses on representative configurations of the whole system that could mimetize the average structure. Thus, the occupancy frequency calculations seem to be a powerful tool in multiscale simulations for sequential QM/MM calculations. It is necessary to emphasize that the investigation about the dynamic contribution ought to be taken into account to select an adequate group of statistically representative conformations. In this way, the configurations produced during the simulation may be radically decreased with no loss of statistical data. We intensively feel that this theoretical investigation may be useful for the design and choice of new acetylcholinesterase reactivators, because the flexibility and stereoselectivity of the reactivation process of AChE by oximes is one of the most serious issues. The occupancy frequency approach can also be extended to other examples, mainly when there is no information about the bioactive conformation from experimental results. Theoretical methods have shown to be a useful tool in the mechanistic study inside enzymes active site.