Novel 1D QSAR approach that allows analysis of non-additive effects of molecular fragments on toxicity has been proposed. Twenty-eight nitroaromatic compounds including some well-known explosives have been chosen for this study. The 50% lethal dose concentration for rats (LD50) was used as the estimation of toxicity in vivo to develop 1D QSAR models on the framework of Simplex representation of molecular structure. The results of 1D QSAR analysis show that even the information about the composition of molecules provides the main trends of toxicity changes. The necessity of consideration of substituents' mutual impacts for the development of adequate QSAR models of nitroaromatics' toxicity was demonstrated. Statistic characteristics for all the developed partial least squares QSAR models, except the additive ones are quite satisfactory (R2=0.81-0.92; Q2=0.64-0.83; R2 test=0.84-0.87). A successful performance of such models is due to their non-additivity i.e. possibility of taking into account the mutual influence of substituents in benzene ring which plays the governing role for toxicity change and could be mediated through the different C-H fragments of the ring. The correspondence between observed and predicted by these models toxicity values is good. This allowing combine advantages of such approaches and develop adequate consensus model that can be used as a toxicity virtual screening tool.