A classification consensus ensemble multitarget neural network model of the dependence of the anxiolytic activity of chemical compounds on the energy of their docking in 17 biotargets was developed. The training set included compounds thathadalready been tested for anxiolytic activity and were structurally similar to the 15 studied nitrogen-containing heterocyclic chemotypes. Seventeen biotargets relevant to anxiolytic activity were selected, taking into account the possible effect on them of the derivatives of these chemotypes. The generated model consistedof three ensembles of artificial neural networks for predicting three levels of anxiolytic activity, with sevenneural networks in each ensemble. A sensitive analysis of neurons in an ensemble of neural networks for a high level of activity made it possible to identify four biotargets ADRA1B, ADRA2A, AGTR1, and NMDA-Glut, which were the most significant for the manifestation of the anxiolytic effect. For these four key biotargets for 2,3,4,5-tetrahydro-11H-[1,3]diazepino[1,2-a]benzimidazole and [1,2,4]triazolo[3,4-a][2,3]benzodiazepine derivatives, eight monotarget pharmacophores of high anxiolytic activity were built. Superposition of monotarget pharmacophores built two multitarget pharmacophores of high anxiolytic activity, reflecting the universal features of interaction 2,3,4,5-tetrahydro-11H-[1,3]diazepino[1,2-a]benzimidazole and [1,2,4]triazolo[3,4-a][2,3]benzodiazepine derivatives with the most significant biotargets ADRA1B, ADRA2A, AGTR1, and NMDA-Glut.