Development of simple QSPR models for the impact sensitivity of nitramines

A General Model of Impact Sensitivity for Nitrogen-Rich Energetic Materials: A Combined Incremental Theory and Genetic Function Approximation Study

In this paper, we report the first example of impact sensitivity prediction based on the genetic function approximation (GFA) as a regression method. The prediction is applicable for a wide variety of chemical families, which include nitro compounds, peroxides, nitrogen-rich salts, heterocycles, etc. Within this work, we have obtained 7 empirical models (with 27-32 basis functions), which all provide 0.80≤R2≤0.83 and 7.2 J≤RMSE≤7.8 J (for 450 training set compounds) and 0.64≤R2≤0.70 and 11.2 J≤RMSE≤12.4 J (for 170 test set compounds). The models were developed using Friedman Lack-of-Fit as a scoring function, which allows avoiding an overfitting. All the models have simple descriptors as basis functions and include linear splines. Furthermore, the applied descriptors do not require expensive calculation procedures, namely, non-empirical quantum-chemical calculations, complex iterative procedures, real space electron density analysis, etc. Most descriptors are based on structural and topological analysis and a part of them require very cheap semi-empirical PM6 calculations. The prediction takes a few minutes as an average, and most of the time is for the structure preparation and manual calculation of the descriptor "Increment", which is based on our recent incremental theory.

The pursuit of accurate predictive models of the bioactivity of small molecules

Property prediction is a key interest in chemistry. For several decades there has been a continued and incremental development of mathematical models to predict properties. As more data is generated and accumulated, there seems to be more areas of opportunity to develop models with increased accuracy. The same is true if one considers the large developments in machine and deep learning models. However, along with the same areas of opportunity and development, issues and challenges remain and, with more data, new challenges emerge such as the quality and quantity and reliability of the data, and model reproducibility. Herein, we discuss the status of the accuracy of predictive models and present the authors' perspective of the direction of the field, emphasizing on good practices. We focus on predictive models of bioactive properties of small molecules relevant for drug discovery, agrochemical, food chemistry, natural product research, and related fields.

A novel model for prediction of stability constants of the thiosemicarbazone ligands with different types of toxic heavy metal ions using structural parameters and multivariate linear regression method

A novel model is presented for reliable estimation of the stability constants of the thiosemicarbazone ligands with different types of toxic heavy metal ions (log β₁₁) in an aqueous solution, which has wide usage in environmental safety and ecotoxicology applications. The biggest reported data of log β₁₁ for 120 metalthiosemicarbazone complexes are used for deriving and testing the novel model. In contrast to available methods where they need the two-dimensional (2D) and three-dimensional (3D) complex molecular descriptors as well as expert users and computer codes, the novel correlation uses four additive and two non-additive structural parameters of thiosemicarbazone ligands. The calculated results of the novel correlation are compared with the outputs of the genetic algorithm with multivariate linear regression method (GA-MLR) as one of the best existing methods, which requires seven complex descriptors. The estimated results for 78 of training as well as 42 of two different test sets were established by external and internal validations. The values of statistical parameters comprising average deviation, average absolute deviation, average absolute relative deviation, absolute maximum deviation, and the coefficient of determination for 73 data of training set of New model/GA-MLR are 0.04/ - 0.25, 1.06/1.31, 14.4/18.7, 3.18/7.92, and 0.830/0.652, respectively. Thus, the predicted results of the new model are worthy as compared to the complex GA-MLR model. Moreover, assessments of various statistical parameters confirm that the new model provides great reliability, goodness-of-fit, accuracy, and precision.

Chemoinformatics for the Safety of Energetic and Reactive Materials at Ineris

The characterization of physical hazards of substances is a key information to manage the risks associated to their use, storage and transport. With decades of work in this area, Ineris develops and implements cutting-edge experimental facilities allowing such characterizations at different scales and under various conditions to study all of the dreaded accident scenarios. This review presents the efforts engaged by Ineris more recently in the field of chemoinformatics to develop and use new predictive methods for the anticipation and management of industrials risks associated to energetic and reactive materials as a complement to experiments. An overview of the methods used for the development of Quantitative Structure-Property Relationships for physical hazards are presented and discussed regarding the specificities associated to this class of properties. A review of models developed at Ineris is also provided from the first tentative models on the explosivity of nitro compounds to the successful application to the flammability of organic mixtures. Then, a discussion is proposed on the use of QSPR models. Good practices for robust use for QSPR models are recalled with specific comments related to physical hazards, notably for regulatory purpose. Dissemination and training efforts engaged by Ineris are also presented. The potential offered by these predictive methods in terms of in silico design and for the development of new intrinsically safer technologies in safety-by-design strategies is finally discussed. At last, challenges and perspectives to extend the application of chemoinformatics in the field of safety and in particular for the physical hazards of energetic and reactive substances are proposed.

A unified model of impact sensitivity of metal azides

A comprehensive theoretical study of 18 metal azides is reported.

Computational study of the structure and properties of bicyclo[3.1.1]heptane derivatives for new high-energy density compounds with low impact sensitivity

To design new high-energy density compounds (HEDCs), a series of new bicyclo[2.2.1]heptane derivatives containing an aza nitrogen atom and nitro substituent were designed and studied theoretically. The density, heat of sublimation and impact sensitivity were estimated by electrostatic potential analysis of the molecular surface. Based on the designed isodesmic reaction, and the reliable heat of formation (HOF) of the reference compounds, HOFs were calculated and compared at B3LYP/6-311G(d,p) and B3P86/6-311G(d,p), respectively. The detonation performances, bond dissociation energies (BDE) and impact sensitivity were calculated to evaluate the designed compounds. The calculated results show that the number of aza nitrogen atoms and NO₂ groups are two important factors for improving HOF, density and detonation properties. Thermal stability generally decreases with increasing nitro groups. And the N-NO₂ bond is the trigger bond for all designed compounds except B8, whose trigger bond is C-NO₂. Importantly, the BDE values are between 86.95 and 179.71 kJ mol^-1 and meet the requirement for HEDCs. Detonation velocity and detonation pressure were found to be 5.77-9.65 km s^-1 and 12.30-43.64 GPa, respectively. After comprehensive consideration of thermal stability, impact sensitivity and detonation properties, A7, A8, B8, C8, D7, E7, F7 and G6 may be considered as potential HEDCs. Especially, A8, B8, C8, and D7 have better detonation properties than the famous caged nitramine CL-20 (D = 9.40 km/s, P = 42.00GPa). Besides, all the designed potential HEDCs have reasonable impact sensitivity. Graphical abstract New high-energy density compounds (HEDCs) with low impact sensitivity (A8, B8, C8 and D7 have better detonation properties than CL-20).

Development of Simple QSPR Models for the Prediction of the Heat of Decomposition of Organic Peroxides

Quantitative structure-property relationships represent alternative method to experiments to access the estimation of physico-chemical properties of chemicals for screening purpose at R&D level but also to gather missing data in regulatory context. In particular, such predictions were encouraged by the REACH regulation for the collection of data, provided that they are developed respecting the rigorous principles of validation proposed by OECD. In this context, a series of organic peroxides, unstable chemicals which can easily decompose and may lead to explosion, were investigated to develop simple QSPR models that can be used in a regulatory framework. Only constitutional and topological descriptors were employed to achieve QSPR models predicting the heat of decomposition, which could be used without any time consuming preliminary structure calculations at quantum chemical level. To validate the models, the original experimental dataset was divided into a training and a validation set according to two methods of partitioning, one based on the property value and the other based on the structure of the molecules by the mean of PCA. Four QSPR models were developed upon the type of descriptors and the methods of partitioning. The 2 models issuing from the PCA based method were highlighted as they presented good predictive power and they are easier to apply than our previous quantum chemical based model, since they do not need any preliminary calculations.