Monte Carlo QSAR models for predicting organophosphate inhibition of acetycholinesterase

A series of 278 organophosphate compounds acting as acetylcholinesterase inhibitors has been studied. The Monte Carlo method was used as a tool for building up one-variable quantitative structure–activity relationship (QSAR) models for acetylcholinesterase inhibition activity based on the principle that the target endpoint is treated as a random event. As an activity, bimolecular rate constants were used. The QSAR models were based on optimal descriptors obtained from Simplified Molecular Input-Line Entry System (SMILES) used for the representation of molecular structure. Two modelling approaches were examined: (1) ‘classic’ training-test system where the QSAR model was built with one random split into a training, test and validation set; and (2) the correlation balance based QSAR models were built with two random splits into a sub-training, calibration, test and validation set. The DModX method was used for defining the applicability domain. The obtained results suggest that studied activity can be determined with the application of QSAR models calculated with the Monte Carlo method since the statistical quality of all build models was very good. Finally, structural indicators for the increase and the decrease of the bimolecular rate constant are defined. The possibility of using these results for the computer-aided design of new organophosphate compounds is presented.

[1]  S. Bryant,et al.  Organophosphate poisoning. , 1968, Lancet.

[2]  Andrey A Toropov,et al.  SMILES-based QSAR model for arylpiperazines as high-affinity 5-HT(1A) receptor ligands using CORAL. , 2013, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[3]  S. Wold,et al.  PLS-regression: a basic tool of chemometrics , 2001 .

[4]  Kevin Barraclough,et al.  I and i , 2001, BMJ : British Medical Journal.

[5]  Paul E. Anderson,et al.  Quantitative structure–activity relationships for organophosphates binding to acetylcholinesterase , 2013, Archives of Toxicology.

[6]  K. Roy,et al.  Further exploring rm2 metrics for validation of QSPR models , 2011 .

[7]  L G Sultatos,et al.  Common mechanism of toxicity: a case study of organophosphorus pesticides. , 1998, Toxicological sciences : an official journal of the Society of Toxicology.

[8]  C. Dary,et al.  Physicochemical and Biological Data for the Development of Predictive Organophosphorus Pesticide QSARs and PBPK/PD Models for Human Risk Assessment , 2004, Critical reviews in toxicology.

[9]  Alexandru T. Balaban,et al.  Chemical graphs , 1979 .

[10]  Emilio Benfenati,et al.  SMILES as an alternative to the graph in QSAR modelling of bee toxicity , 2007, Comput. Biol. Chem..

[11]  M. Karelson,et al.  Quantum-Chemical Descriptors in QSAR/QSPR Studies. , 1996, Chemical reviews.

[12]  P. Roy,et al.  Exploring the impact of size of training sets for the development of predictive QSAR models , 2008 .

[13]  J. Jeyaratnam.  Acute pesticide poisoning: a major global health problem. , 1990, World health statistics quarterly. Rapport trimestriel de statistiques sanitaires mondiales.

[14]  Ovidiu Ivanciuc,et al.  QSAR Comparative Study of Wiener Descriptors for Weighted Molecular Graphs , 2000, J. Chem. Inf. Comput. Sci..

[15]  A Goulon,et al.  Predicting activities without computing descriptors: graph machines for QSAR , 2007, SAR and QSAR in environmental research.

[16]  A. Tropsha,et al.  Beware of q2! , 2002, Journal of molecular graphics & modelling.

[17]  Kunal Roy,et al.  Comparative QSARs for antimalarial endochins: Importance of descriptor-thinning and noise reduction prior to feature selection , 2011 .

[18]  Paola Gramatica,et al.  Principles of QSAR models validation: internal and external , 2007 .

[19]  R. S. Thomas,et al.  Approaches to developing alternative and predictive toxicology based on PBPK/PD and QSAR modeling. , 1998, Environmental health perspectives.

[20]  S. E. Kunz,et al.  Insecticides and acaricides: resistance and environmental impact. , 1994, Revue scientifique et technique.

[21]  F. Prado-Prado,et al.  Predicting antimicrobial drugs and targets with the MARCH-INSIDE approach. , 2008, Current topics in medicinal chemistry.

[22]  Jerzy Leszczynski,et al.  QSAR model as a random event: A case of rat toxicity. , 2015, Bioorganic & medicinal chemistry.

[23]  Michael Eddleston,et al.  Management of acute organophosphorus pesticide poisoning , 2008, The Lancet.

[24]  Umakant Sharma,et al.  Insect vectors of Leishmania: distribution, physiology and their control. , 2008, Journal of vector borne diseases.

[25]  Emilio Benfenati,et al.  Interpretation of Quantitative Structure-Property and -Activity Relationships , 2001, J. Chem. Inf. Comput. Sci..

[26]  A. Toropova,et al.  QSAR models for HEPT derivates as NNRTI inhibitors based on Monte Carlo method. , 2014, European journal of medicinal chemistry.

[27]  Chris Morley,et al.  Open Babel: An open chemical toolbox , 2011, J. Cheminformatics.

[28]  Andrey A Toropov,et al.  SMILES‐Based QSAR Models for the Calcium Channel‐Antagonistic Effect of 1,4‐Dihydropyridines , 2013, Archiv der Pharmazie.

[29]  E Benfenati,et al.  Additive SMILES-based optimal descriptors in QSAR modelling bee toxicity: Using rare SMILES attributes to define the applicability domain. , 2008, Bioorganic & medicinal chemistry.

[30]  Ovidiu Ivanciuc,et al.  Chemical graphs, molecular matrices and topological indices in chemoinformatics and quantitative structure-activity relationships. , 2013, Current computer-aided drug design.

[31]  Kunal Roy,et al.  On some aspects of validation of predictive quantitative structure–activity relationship models , 2007, Expert opinion on drug discovery.

[32]  E Benfenati,et al.  CORAL: building up the model for bioconcentration factor and defining it's applicability domain. , 2011, European journal of medicinal chemistry.

[33]  Fukuto Tr.  Mechanism of action of organophosphorus and carbamate insecticides. , 1990 .

[34]  B. Eskenazi,et al.  Exposures of children to organophosphate pesticides and their potential adverse health effects. , 1999, Environmental health perspectives.

[35]  A P Toropova,et al.  QSPR modeling of octanol/water partition coefficient for vitamins by optimal descriptors calculated with SMILES. , 2008, European journal of medicinal chemistry.

[36]  T. R. Fukuto,et al.  Mechanism of action of organophosphorus and carbamate insecticides. , 1990, Environmental health perspectives.

[37]  E Benfenati,et al.  SMILES-based optimal descriptors: QSAR modeling of carcinogenicity by balance of correlations with ideal slopes. , 2010, European journal of medicinal chemistry.

[38]  Kunal Roy,et al.  Some case studies on application of “rm2” metrics for judging quality of quantitative structure–activity relationship predictions: Emphasis on scaling of response data , 2013, J. Comput. Chem..

[39]  N. Munro.  Toxicity of the Organophosphate Chemical Warfare Agents GA, GB, and VX: Implications for Public Protection. , 1994, Environmental health perspectives.

[40]  J. Hemingway,et al.  Resistance management strategies in malaria vector mosquito control. Baseline data for a large‐scale field trial against Anopheles albimanus in Mexico , 1998, Medical and veterinary entomology.

[41]  Kunal Roy,et al.  QSAR Studies of CYP2D6 Inhibitor Aryloxypropanolamines Using 2D and 3D Descriptors , 2009, Chemical biology & drug design.