Introduction. Molecular modeling methods are very popular in the scientific community at the present time. The value of the pharmacological action depends on the affinity of the substance to the biological target. Molecular docking makes it possible to assess the degree of affinity of the studied compound with the active center of the molecular target. The enzyme cyclooxygenase (COX) plays a key role in the cascade of synthesis of proinflammatory cytokines and, as a consequence, in pain.Aim. To identify the dependence "structure – analgesic activity" using the method of molecular docking for cyclooxygenase type 1 and type 2 in the series of 5-N-arylaminocarbonyl-6-(get)aryl-4-methyl-1,2,3,6-tetrahydropyrimidine-2-thions.Materials and methods. 19 compounds 5-N-arylaminocarbonyl-6-(get)aryl-4-methyl-1,2,3,6-tetrahydropyrimidine-2-thions are objects of research The study of the interaction of tetrahydropyridine derivatives with COX 1 and 2 was carried out by the method of molecular docking by the AutoDock 4 program using scoring functions.Results and discussion. Molecular docking of 5-N-arylaminocarbonyl-6-(get)aryl-4-methyl-1,2,3,6-tetrahydropyrimidine-2-thions (I–XIX) with COX 1 and 2 is described. Quantitative studies of the "structure – analgesic activity" in the studied series of compounds of the dependence of experimental values of analgesic activity (ААex.) on scoring functions (BeCOX1, KiCOX1, BeCOX2, KiCOX2) and physico-chemical descriptors (log Ppred., рКаpred., рКvpred.) were performed. By checking on an independent sample of 5 compounds, equation No. 3 (ААpred. 3 = 32,6215 – 4,4894 × Becox1 + 0,0066 × Kicox1 + 3,6032 × log Ppred. (R = 0,854, F = 9,01, S = 7,73, Q2LOO = 0,53) was found, with a high value of the correlation coefficient of the predicted values ААex. with experimental values (Rpred. = 0,878) and a minimum forecast error (Spred. = 6,74).Conclusion. The "structure-activity" models for predicting analgesic activity in a series of 5-N-arylaminocarbonyl-6-(get)aryl-4-methyl-1,2,3,6-tetrahydropyrimidine-2-thions are obtained. The result of the prediction of biological activity is confirmed by the values of the correlation coefficient (R) obtained when testing models on independent samples.