INTRODUCTION

In 2002, the global community was faced with the SARS-CoV coronavirus, with caused an epidemic of atypical pneumonia (severe acute respiratory syndrome). The first cases of infection with the SARS-CoV coronavirus were detected in Southern China, before the epidemic spread to 29 countries, as a result of which more than 8000 people were infected, while 916 persons died [1]. In 2012, the second epidemic, caused by the spread of a coronavirus, in this case the MERS-CoV virus, which is characterized by higher lethality, was recorded in Saudi Arabia. By the end of August 2015, 574 persons died among the total number of infected people (1511 patients); i.e., the lethality from the MERS-CoV coronavirus was over 37% in contrast to SARS-CoV, for which the lethality was estimated at 10%.

A number of experts forecast the emergence of a threat for the human race, which will be caused by a new type of coronavirus [3]. The authors of this study proposed that the new coronavirus, circulating in the Chinese populations of horseshoe bats, will bind to the human angiotensin converting enzyme II (ACE2) followed by the efficient replication in the cell of the respiratory system.

Unfortunately, the prognoses of the above-mentioned experts turned out to be right, and the Chinese public health authorities recorded the first case of infection with the new SARS-CoV-2 coronavirus on December 8, 2019 [4]. The level of the lethality from the new SARS-CoV-2 coronavirus is estimated to be lower (about 7%) than for SARS-CoV and MERS-CoV; however, it was shown that only 48 days are required for infection of the first 1000 patients with the SARS-CoV-2 coronavirus, while 130 days are required for infection with SARS-CoV, and two-and-a-half years for MERS-CoV [5].

The pandemic caused by the new SARS-CoV-2 coronavirus represents a serious medical and socioeconomic problem for all mankind.

The drug Favipiravir recommended in the Russian Federation for the treatment of COVID-19 [6], according to a number of researchers, has a teratogenic effect [7]. Thus, the search and development of highly effective and safe drugs that can stop the spread of the COVID-19 pandemic is an urgent issue.

Chemoinformatics methods can provide significant assistance in reducing the time and financial costs in repositioning and developing new drugs [8–10]. Since the emergence of the new SARS-CoV-2 coronavirus, a number of studies using the methods of molecular docking, molecular dynamics, and pharmacophore analysis have been conducted in this direction [11–22].

Currently, the molecular structure of potential SARS-CoV-2 inhibitors is being considered from different points of view. The resulting diverse information is of significant interest to the world scientific and medical community. In the studies conducted, the most significant target for developing drugs is the main protease of the virus Mpro, also known as 3-chymotrypsin-like protease (3CLpro), which plays a key role in replication of coronaviruses. It was found that this enzyme, being evolutionarily conservative, did not undergo significant changes in contrast to the main Mpro protease of the SARS-CoV virus, which caused an outbreak of acute respiratory syndrome in 2002–2003 [23–27]. Previous studies also indicate the conservatism of the Mpro sequences and spatial Mpro structures of different types of coronaviruses [28]. At the same time, close homologues of this enzyme have not been identified in the human body, which has a positive effect on the specificity and a decrease in the number of potential side effects of inhibitors of the main Mpro protease [29].

In the study [25], a virtual screening of the DrugBank library of chemical compounds was performed based on the similarity of the binding sites of the main Mpro protease of SARS-CoV and SARS-CoV-2 viruses using the method of molecular docking [30]. As a result, a list of ten potential inhibitors of the main Mpro protease was proposed, which, according to the authors [25], are the most promising for combating SARS-CoV-2.

In the study [31], compounds included in the list of drugs of traditional Chinese medicine were initially selected. For these compounds, an assessment of their pharmacokinetic characteristics such as adsorption, distribution, metabolism, and excretion is given. The most promising compounds were studied using molecular docking. The next step was the selection of medicinal herbs that contain at least two compounds proposed in the course of molecular docking. As a result of the study, the authors of [31] identified 26 medicinal herbs of Chinese medicine that are potentially promising for the treatment of the COVID-19 disease caused by the SARS-CoV-2 coronavirus.

The publication [23], which presents the results of the consensus in vitro and in silico screening, deserves special attention. The authors studied a database of more than 10 000 compounds, for which the binding to Mpro of the SARS-CoV-2 coronavirus was experimentally measured by the method of fluorescent resonance energy transfer. The most promising seven compounds were further investigated for their ability to prevent the infection of cells with the SARS-CoV-2 virus. The in silico screening of these 10 000 compounds was performed using the Glide v8.2 and Maestro software (Schrödinger). According to the consensus results, the most promising were ebselen (2-phenyl-1,2-benzoselenazol-3-one, CAS no. 60940-34-3), carmofur (1-hexylcarbamoyl-5-fluorouracil, CAS no. 61422-45-5), compound TDZD-8 (2-methyl-4-(phenylmethyl)-1,2-athiadiazolidine-3,5-dione, CAS no. 327036-89-5), and N3 peptidomimetic, the previously proposed SARS-CoV inhibitor, which covalently binds to Mpro according to the Michael reaction.

The authors of [32] conducted a virtual screening of 1.3 billion molecules in order to identify the most active inhibitors of the main Mpro protease of the SARS-CoV-2 virus. In this study, the recently developed Deep Docking algorithm, which integrated classical docking and the methodological foundations for constructing QSAR (Quantitative Structure–Activity Relationship) models, was used; this allowed them to increase the screening performance compared with traditional docking methods. The authors note that the scoring functions are determined by the used docking methods, and the QSAR models were used to optimize the virtual screening. Based on the results of the virtual screening, a hit list of 1000 compounds was proposed, which is available for free download at https://drive.google.com/drive/ folders/1xgA8ScPRqIunxEAXFrUEkavS7y3tLIMN.

Using deep learning the authors performed the study [33] in which models that describe the structure of compounds using character strings composed according to the rules of SMILES [34] were developed. The principle of the used prediction method is based on a technology called natural language processing [35], used in the analysis of human speech by a computer, but in this case the language is a string of characters written according to the SMILES rules and the sequence of the target protein. In order to identify and analyze regularities, convolutional neural networks were used [36]. The study predicted the activity of inhibitors of the main Mpro protease, RNA replicase (RNA-dependent RNA-polymerase, RdRP), helicase, and a number of other enzymes of the SARS-CoV-2 virus. Molecular docking (in particular, the AutoDock Vina v.1.1.2 program) was used for the comparative study. As a result, the authors identified three drugs against HIV (ritonavir, atazanavir, efavirenz), as well as the antiviral agent ganciclovir.

The TMPRSS2 protease (Transmembrane protease, serine 2, membrane-bound serine protease) serves as another target for the fight against the coronavirus; its inhibitors can prevent the entry of the virion into the cell [37]. However, the number of studies devoted to the computer modeling of TMPRSS2 protease inhibitors is significantly less than that of the main Mpro protease. Thus, we can note the publication [38], in which a virtual screening of a database containing more than 30 000 natural compounds was performed using molecular docking and pharmacophore analysis. For the selected 12 compounds, the authors evaluated adsorption, distribution, metabolism, elimination, and toxicity. The small molecule compound geniposide (CAS no. 24512-63-8) turned out to be the most promising.

Thus, the majority of works was carried out using the method of molecular docking, which, like any method of research, has limitations. In particular, the methodological difficulties of docking are associated with considering conformations of the ligand, the choice of methods for constructing the scoring function, and the flexibility of receptors. The main methods of molecular docking and their inherent limitations are described in detail in the reviews [39–41].

QSAR is an alternative method of computer drug development, which has successfully solved various problems [42]. The literature contains information on the similarity of the binding sites of the main Mpro protease of SARS-CoV and SARS-CoV-2 viruses, which were confirmed in the course of independent studies [23–27]. In relation to this, we assumed that the assessment of potential SARS-CoV-2 inhibitors during drug development can be realized using QSAR models of SARS-CoV inhibitors.

The authors of [43] developed a QSAR model of SARS-CoV inhibitors using 3D-QSAR methods (CoMFA, CoMSIA), the limiting feature of which is the ambiguity of the three-dimensional alignment of the structures of the studied compounds [44]. In another study [45], 33 QSAR models of SARS-CoV inhibitors were developed, but the authors did not provide indicators of their predictive ability, assessed using compounds of the test sample.

In accordance with the fifth principle of QSAR modeling, developed by the OECD expert group [46], the interpretation of the obtained models is desirable. In the reviewed publications [43, 45], there is no structural interpretation of the QSAR models, which does not allow conducting molecular design and limits the use of simulation results for studying the mechanisms of biological reactions [47].

Recently, a study was published [48] in which acceptable QSAR models of inhibitors of the main Mpro protease of the SARS-CoV virus, which were developed using PaDEL and Dragon descriptors, and the method of multiple linear regression (MLR) were proposed in the search for effective drugs against COVID-19. Using the developed QSAR models, the authors conducted a virtual screening of more than 50 000 different compounds in order to identify the most active inhibitors of the main Mpro protease of the virus. Based on the proposed regression equations, namely, the contributions of some significant descriptors, the authors of [48] analyzed the effect of the structural features of the studied compounds on a change in the inhibitory activity.

In the study [49], adequate QSAR models of SARS-CoV inhibitors, in the course of the structural interpretation of which molecular fragments that decrease and increase this type of activity were identified, were developed for 54 peptidomimetics. To construct the models, the authors of [49] also used the MLR method.

The regression method of the MLR data analysis applied in [48, 49] can give adequate results only in the presence of a linear relationship between the structure and activity [50]. One of the ways to overcome this disadvantage can be the use of nonparametric methods, in particular, various methods of machine learning (ML).

The study [27], in the course of which high-quality classification-based QSAR models of SARS-CoV inhibitors were developed, deserves special attention. The reliability of the constructed QSAR models was confirmed by subsequent experimental studies, as a result of which some compounds demonstrated high activity and were recommended for further study. In parallel with the QSAR analysis, the authors of [27] carried out a study using methods of molecular docking, while the revealed unacceptably low level of predictive ability did not allow considering the results of studying SARS-CoV inhibitors by molecular docking methods. Undoubtedly, the study [27] is very successful, but it lacks the structural interpretation of QSAR models.

The present study consisted of the following stages:

(1) construction of QSAR models of the main Mpro protease inhibitors of SARS-CoV;

(2) performing a virtual screening of the most promising compounds of potential drugs for the treatment of COVID-19;

(3) structural interpretation of QSAR models and rational molecular design of the main Mpro protease inhibitors.

EXPERIMENTAL

The well-known CHEMBL database (ID: CHEMBL3927) [51] served as the source for the sample formation for QSAR modeling. Inorganic compounds, polymers, mixtures, and compounds in the salt form were removed from the obtained sample. The final set of inhibitors of the main Mpro protease of SARS-CoV contained 65 compounds.

The experimental values of the activity of inhibitors of the SARS-CoV main Mpro protease, expressed in terms of the half-maximal inhibitory concentration (IC50, nM or μM), which were given in the primary sources, were converted (1) into the negative common logarithm of the pIC50 value, which is generally accepted in QSAR studies and is used in cases when a linear increase in the concentration causes an exponential increase in the effect:

$${\text{pI}}{{{\text{C}}}_{{50}}} = - {\text{log}}\left( {{\text{I}}{{{\text{C}}}_{{50}}}} \right).$$
(1)

The exported experimental values and structures of compounds are given in the Appendix A (Table A1).

To describe the molecular structure, 2D simplex descriptors, which are calculated in the simplex representation of the molecular structure (SIRMS), were used [52]. Within the SIRMS, a molecule is considered as a system of various simplexes, tetraatomic molecular fragments of a fixed structure (Fig. 1).

Fig. 1.
figure 1

An example of 2D generation of simplex descriptors for alanine at the 2D level using differentiation of atoms by their partial charges

Full size image

The descriptor in this case is the number of simplices of a certain type. At the 2D level, atoms (vertices of a simplex) are differentiated not only by the nature “label” of the atom but also considering different physicochemical properties (partial charge on the atom, lipophilicity, refraction, and the ability to act as a donor or acceptor of hydrogen during the formation of a hydrogen bond).

Structural interpretation was performed in accordance with the approach [53], in which the contribution of the studied fragment (C) was calculated by the difference between the calculated values of the activity for the parent structure (A) and the hypothetical structure (B) obtained by removing the studied fragment (C) from the parent structure (A) (Fig. 2).

Fig. 2.
figure 2

The principle of structural interpretation used. W(C) is the contribution of the fragment (C); X(A) is the predicted activity value of the parent structure (A); X (B) is the predicted activity value for a hypothetical structure (B)

Full size image

When simulating using simplex descriptors, we used the Scikit-learn package [54] for the Python programming language, which implements the methods of partial least squares (PLS), random forest (RF), support vector machine (SVM), and gradient boosting method (GBM).

Due to the small number of studied compounds and their structural diversity, a five-fold internal cross-validation (CV) was performed. For this, all compounds of the training sample are randomly divided into five parts. Then, a QSAR model is built (trained) on four pieces of data combined into a training sample, and the rest of the data is used as an external test sample; i.e., the predictive ability of the model is checked on the compounds of this group. This procedure is repeated 5 times; as a result, each of the five portions of the data is sequentially used for testing. Note that the studied compounds are never simultaneously used as a part of both the training and the external test set.

During QSAR modeling, the inclusion of the compounds in the applicability domain (AD) [55] was considerered for test samples, while if the value of at least one descriptor went beyond its minimum or maximum value for the training sample, then the compound of the test sample containing this descriptor was not included in the bounding box. This approach for QSAR modeling using simplex descriptors is implemented in the form of the SPCI software, which is freely available at http://qsar4u.com/pages/sirms_ qsar.php.

In addition, we used the OCHEM Internet platform ((https://ochem.eu) for QSAR analysis. The best modeling results were achieved using a number of descriptors (ALogPS, OEstate, Dragon, CDK, ISIDA Fragments, StructuralAlerts, alvaDesc) and the RF and associative neural networks (ASNN) methods. The consensus model was constructed by averaging the predictions of the best individual models. In this case, the applicability domain was assessed using the concept of the distance to the model (in particular, the CLASS-LAG approach). A brief description of the used methods and descriptors, as well as links to the original works are given in the OCHEM user manual [56].

The OCHEM internet resource implements the method of molecular pairs [57], which also allows us to interpret models constructed on any descriptors.

The assessment of the accuracy and predictive ability of the models proposed in this study and their comparison with other QSAR models was performed based on the following criteria.

1. The coefficient of determination (R2):

$${{R}^{2}} = 1 - \frac{{\sum\limits_{i = 1}^m {{{{\left( {{{y}_{i}} - {{{\hat {y}}}_{i}}} \right)}}^{2}}} }}{{\sum\limits_{i = 1}^m {{{{\left( {{{y}_{i}} - {{y}_{{i{\text{mean}}}}}} \right)}}^{2}}} }},$$
(2)

where \({{\hat {y}}_{i}}\) is the calculated value of the property for the ith molecule, yi is the observed (experimental) value of the property for the ith molecule, m is the number of molecules in the sample, and yi mean is the mean value of the observed property.

2. Root mean square error (RMSE):

$${\text{RMSE}} = \sqrt {\frac{{\sum\limits_{i = 1}^m {{{{\left( {{{y}_{i}} - {{{\hat {y}}}_{i}}} \right)}}^{2}}} }}{m}} .$$
(3)

Considering the fact that the QSAR modeling mentioned above requires the use of a large number (hundreds and thousands) of descriptors, for comparison, regression models containing a small number of HYBOT variables [59] with the addition of a number of topological and fractal characteristics [60] were constructed using the RF method [58]. In this case, AD was assessed using the interval method.

For the virtual screening, we used the DrugBank database [30], represented by FDA-approved drugs. Inorganic compounds, polymers, mixtures, and compounds in the salt form were removed from the exported DrugBank database. The final sample for the virtual screening contained 2087 FDA-approved drugs.

RESULTS AND DISCUSSION

The results of the QSAR modeling are shown in Tables 1–3. The consensus model is freely available at the link given in Table 1. All the constructed models have satisfactory statistical characteristics and possess comparable predictive power.

Table 1.   Statistical characteristics of QSAR models developed using the OCHEM internet resource
Full size table

For the consensus QSAR models obtained (Tables 1, 2), a structural interpretation was performed. When interpreting the consensus model built using simplex descriptors (Table 2), the contributions of molecular fragments to the activity of inhibitors of the SARS-CoV main Mpro protease were determined (Fig. 3).

Table 2.   Statistical characteristics of QSAR models developed using 2D simplex descriptors
Full size table
Fig. 3.
figure 3

Contributions of molecular fragments to the ability of compounds to inhibit the SARS-CoV main Mpro protease. A, the place of the fragment’s attachment to the other part of the molecule; f, the order number of the fragment; M, the number of compounds containing the given fragment; N, the number of detections of the corresponding fragment in the sample.

Full size image
Table 3.   Statistical characteristics of QSAR model developed using HYBOT, topological, and fractal descriptors
Full size table

This set of molecular fragments is formed from the standard functional groups (51 fragments) and six molecular fragments obtained during the automatic fragmentation of compounds of the training set using the SPCI program using the SMART template [#6+0;!$(*=,#[!#6])]!@!=!#[*], which encodes breakable bonds [53]. Only those molecular fragments that were found in three or more compounds were subjected to interpretation, which, from our point of view, allowed us to focus on the fragments that stably affect the inhibitory activity and to avoid, to some extent, the influence of random factors, for example errors in the experimental data or predicted values of the activity and contributions of fragments. The complete list of identified molecular fragments in the form of SMARTS with the calculated average contributions to the activity is given (Appendix A, Table A2).

The interpretation allowed us to quantitatively describe and rank the effect of molecular fragments on the change in the activity of SARS-CoV Mpro inhibitors and detail the molecular environment of the known functional groups, highlighting derivative fragments that increase and decrease contributions to the indicated type of the activity. For example, when detailing pyrimidine, the 2-sulfanylpyrimidin-4-ol molecular fragment (f9 in Fig. 3), which significantly reduces the activity of SARS-CoV Mpro inhibitors, was isolated. In this case, carboxyl derivatives of furan and pyridine (fragments f1 and f2 in Fig. 3), on the contrary, increase the activity of SARS-CoV Mpro inhibitors.

The interpretation was also performed for the consensus model (Table 1), built using the OCHEM internet resource. Table 4 shows the results of the interpretation, according to which the inhibitory activity increases under the substitution of hydrogen atoms with chlorine or methyl group. An increase in the activity of SARS-CoV Mpro inhibitors is also observed when phenyl and n-propyl radicals are replaced by naphthyl radicals. The results of the interpretations for the consensus models described above consistently indicate an increase in the activity of SARS-CoV Mpro inhibitors under the substitution of fragments containing iodine (f5) by the carbamoyl group (f4).

Table 4.   The interpretation results for the consensus model developed using the OCHEM internet resource
Full size table

Considering the trends in the effects of the structure of compounds on the change in the activity revealed during the interpretation, we carried out a rational molecular design and proposed a number of promising agents against COVID-19. In this case, molecular fragments that reduce the activity were replaced by the fragments which increase the activity of Mpro inhibitors according to the interpretation results. As a result, hypothetical compounds (Table 5, substances 2, 4, 6), which possess a significant calculated inhibitory activity and fall into the applicability domain of the consensus QSAR model, developed using OCHEM, were proposed. For example, when the residue of 6-methyl-2-sulfanylpyridine-4-ol (compound 1, Table 5) is substituted with carboxyl derivative of pyridine (compound 2, Table 5), a significant increase in the activity of the SARS-CoV main Mpro protease inhibitors is noted. Also an increase in the activity is characteristic for the substitution of a fragment, containing nitrile (compound 3), with trifluoromethyl (compound 4) or of the 4-(1,3-thiazol-4-yl)pyrimidine-2-thiol residue (compound 5, Table 5) with the above-mentioned carboxyl pyridine derivative (compound 6). It should be noted that during molecular design in these examples the results of the interpretation of the QSAR model built using simplex descriptors were considered, while the prediction of the inhibitor activity was performed using the QSAR model built by the OCHEM internet resource.

Table 5.   The results of molecular design
Full size table

When determining the strategies of synthesis and testing, it is important to evaluate various types of toxicity and lipophilicity in addition to the target property (activity), which are important factors when deciding whether to recommend the use of a compound as an active substance of the drug. For this purpose, the acute toxicity (LD50) after oral administration to rats and the probability of mutagenicity (the Ames test) was assessed for compounds 1–9 using the T.E.S.T. v.4.2. program, developed by experts from the Environmental Protection Agency of the United States [61]. Also, using the swissADME Internet platform of the Swiss Bioinformatics Institute (http://www.swissadme.ch/) [62], the lipophilicity (log Po/w), compliance with Lipinski’s rules [63], the presence of PAINS fragments [64], the synthetic availability on a ten-point scale (0 is the maximum degree of synthetic availability, 10 is the minimum degree of synthetic availability) [65] were assessed for these compounds, which is extremely important for the proposed, but not yet synthesized, compounds. The prediction results are shown in Table 6, from which it can be seen that compounds 2, 4, 6, 8, and 9 proposed in the course of molecular design have comparable synthetic availability in comparison with the synthesized substances 1, 3, 5, and 7. All substances satisfy Lipinski’s rules of five, except for compound 4, and do not contain PAINS fragments. The studied compounds are characterized by a wide range of lipophilicity values, which should be considered when studying pharmacokinetics and choosing dosage forms. According to the calculations performed, compound 7 has a nonzero probability of mutagenicity, which can reduce its attractiveness as a lead compound, even though it has the maximum experimentally measured inhibitory activity (Table 5) among the compounds of the exported sample from the CHEMBL database (ID: CHEMBL3927). When modifying compound 7, a hypothetical compound 9, which does not have the probability of mutagenicity according to the calculated data, while the indicators of LD50 and the inhibitory activity are comparable with the initial compound 7, was proposed. In addition, the modification of compound 1 into compound 2 allowed us not only to increase the inhibitory activity by almost two orders, reaching comparable values with the most active substances in the exported sample but also to reduce the toxicity (LD50) by a factor of almost 2.7. Thus, compound 2 can be recommended for synthesis and further testing.

Table 6.   The assessment of some physicochemical properties, types of toxicity and synthetic availability of compounds studied within molecular design (the structures of compounds are presented in Table 5)
Full size table

Since the synthesis of new compounds and their clinical trials take a long time, the most important means of combating a new, rapidly spreading pandemic is the repositioning of approved drugs that have passed all the necessary clinical studies. In order to identify promising inhibitors of the SARS-CoV-2 main Mpro protease, 2087 FDA-approved drugs were screened. The consensus model built by the OCHEM expert system was used for screening, since it has better statistical characteristics and can be used by all interested persons for the virtual screening of their own sets of compounds. The QSAR model obtained using simplex descriptors was not used due to the peculiarities of the method described above for determining the applicability domain, which severely limits the structural space of the model.

In the course of the virtual screening, ten compounds that are within the applicability domain of the consensus QSAR model developed using OCHEM and have the highest calculated inhibitory activity were proposed (Table 7). Antitumor, antiprotozoal, antifungal, antibacterial, antihypertensive, and anti-inflammatory drugs are among these compounds. Tilbroquinol, Chloroxine, and Chlorquinaldol, which are halogenated quinoline derivatives as are the well-known chloroquine and hydroxychloroquine used in the treatment of COVID-19, are of particular interest [66].

Table 7.   The most promising FDA-approved drugs for inhibiting SARS-CoV-2 replication according to the results of virtual screening
Full size table

According to the data of the virtual screening, the highest inhibitory activity among the FDA-approved drugs is possessed by Pexidartinib, which is an anticancer agent, a tyrosine kinase inhibitor. In the study [12], conducted using molecular docking and molecular dynamics, the antitumor agent Neratinib, which blocks the functioning of receptor tyrosine kinases, was also proposed as a promising inhibitor of the SARS-CoV-2 main Mpro protease. The conclusions of the authors [12] are based on the assumption of a similar binding of this antitumor agent to the cysteine residue in the active centers of the kinase domains of receptor tyrosine kinases and the SARS-CoV-2 main Mpro protease. Another antitumor agent, carmofur, was also isolated as a promising inhibitor of the SARS-CoV-2 main Mpro protease according to the results of high-throughput screening in the above-mentioned study [23]. Recent additional studies using X-ray diffraction analysis [67] describe the mechanism of the inhibition of the SARS-CoV-2 main Mpro protease by carmofur through covalent binding to the cysteine residue Cys145 in the active center. Based on the foregoing, the proposal to repurpose pexidartinib, identified during the virtual screening, for the treatment of COVID-19 seems logical. It should be noted that the confirmation of the effectiveness of the drugs proposed for repurposing in the fight against COVID-19 requires significant additional experimental research. Drugs should be taken only according to the medical prescription by the physician.

Thus, in the course of computational experiments using conceptually different descriptors and machine learning methods, acceptable QSAR models of the main Mpro protease inhibitors were developed.

The structural interpretation of the QSAR models allowed us to reveal the common regularities in the effect of the structure of chemical compounds on their inhibitory activity by isolating molecular fragments and transformations that increase and decrease the activity of SARS-CoV inhibitors. The results of the structural interpretation were used to perform rational molecular design, in the course of which a number of promising compounds for combating COVID-19 were proposed.

The virtual screening of FDA-approved drugs identified ten substances that can be recommended for repurposing as drugs against the new coronavirus infection.

The results of this study can help to reduce financial, time, and labor costs when determining the strategy for the development of new drugs and repositioning existing drugs that are SARS-CoV-2 inhibitors.