In silico local QSAR modeling of bioconcentration factor of organophosphate pesticides

The persistent and accumulative nature of the pesticide of indiscriminate use emerged as ecotoxicological hazards. The bioconcentration factor (BCF) is one of the key elements for environmental assessments of the aquatic compartment. Limitations of prediction accuracy of global model facilitate the use of local predictive models in toxicity modeling of emerging compounds. The BCF data of diverse organophosphate (n = 55) was collected from the Pesticide Properties Database and used as a model data set in the present study to explore physicochemical properties and structural alert concerning BCF. The structures were downloaded from Pubchem, ChemSpider database. Two splitting techniques (biological sorting and structure-based) were used to divide the whole dataset into training and test set compounds. The QSAR study was carried out with two-dimensional descriptors (2D) calculated from PaDEL by applying genetic algorithm (GA) as chemometric tools using QSARINS software. The models were statistically robust enough both internally as well as externally (Q²: 0.709-0.722, Q² _Ext: 0.717-0.903, CCC: 0.857-0.880). Overall molecular mass, presence of fused, and heterocyclic ring with electron-withdrawing groups affect the BCF value. The developed models reflected extended applicability domain (AD) and reliable predictions than the reported models for the studied chemical class. Finally, predictions of unknown organophosphate pesticides and the toxic nature of unknown organophosphate pesticides were commented on. These findings may be useful for the scientific community in prioritizing high potential pesticides of organophosphate class.