In silico approach for the development of novel antiviral compounds based on SARS-COV-2 protease inhibition

The COVID-19 pandemic emerged in 2019, bringing with it the need for greater stores of effective antiviral drugs. This paper deals with the conformation-independent, QSAR model, developed by employing the Monte Carlo optimization method, as well as molecular graphs and the SMILES notation-based descriptors for the purpose of modeling the SARS-CoV-3CLpro enzyme inhibition. The main purpose was developing a reproducible model involving easy interpretation, utilized for a quick prediction of the inhibitory activity of SAR-CoV-3CLpro. The following statistical parameters were present in the best-developed QSAR model: (training set) R(2) = 0.9314, Q(2) = 0.9271; (test set) R(2) = 0.9243, Q(2) = 0.8986. Molecular fragments, defined as SMILES notation descriptors, that have a positive and negative impact on 3CLpro inhibition were identified on the basis of the results obtained for structural indicators, and were applied to the computer-aided design of five new compounds with (4-methoxyphenyl)[2-(methylsulfanyl)-6,7-dihydro-1H-[1,4]dioxino[2,3-f]benzimidazol-1-yl]methanone as a template molecule. Molecular docking studies were used to examine the potential inhibition effect of designed molecules on SARS-CoV-3CLpro enzyme inhibition and obtained results have high correlation with the QSAR modeling results. In addition, the interactions between the designed molecules and amino acids from the 3CLpro active site were determined, and the energies they yield were calculated. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11696-022-02170-8.