Docking and other computing tools in drug design against SARS-CoV-2

The use of computer simulation methods has become an indispensable component in identifying drugs against the SARS-CoV-2 coronavirus. There is a huge body of literature on application of molecular modelling to predict inhibitors against target proteins of SARS-CoV-2. To keep our review clear and readable, we limited ourselves primarily to works that use computational methods to find inhibitors and test the predicted compounds experimentally either in target protein assays or in cell culture with live SARS-CoV-2. Some works containing results of experimental discovery of corresponding inhibitors without using computer modelling are included as examples of a success. Also, some computational works without experimental confirmations are also included if they attract our attention either by simulation methods or by databases used. This review collects studies that use various molecular modelling methods: docking, molecular dynamics, quantum mechanics, machine learning, and others. Most of these studies are based on docking, and other methods are used mainly for post-processing to select the best compounds among those found through docking. Simulation methods are presented concisely, information is also provided on databases of organic compounds that can be useful for virtual screening, and the review itself is structured in accordance with coronavirus target proteins.

[Development of antiviral drugs based on inhibitors of the SARS-COV-2 main protease]

Docking and quantum-chemical methods have been used for screening of drug-like compounds from the own database of the Voronezh State University to find inhibitors the SARS-CoV-2 main protease, an important enzyme of the coronavirus responsible for the COVID-19 pandemic. Using the SOL program more than 42000 3D molecular structures were docked into the active site of the main protease, and more than 1000 ligands with most negative values of the SOL score were selected for further processing. For all these top ligands, the protein-ligand binding enthalpy has been calculated using the PM7 semiempirical quantum-chemical method with the COSMO implicit solvent model. 20 ligands with the most negative SOL scores and the most negative binding enthalpies have been selected for further experimental testing. The latter has been made by measurements of the inhibitory activity against the main protease and suppression of SARS-CoV-2 replication in a cell culture. The inhibitory activity \of the compounds was determined using a synthetic fluorescently labeled peptide substrate including the proteolysis site of the main protease. The antiviral activity was tested against SARS-CoV-2 virus in the Vero cell culture. Eight compounds showed inhibitory activity against the main protease of SARS-CoV-2 in the submicromolar and micromolar ranges of the IC50 values. Three compounds suppressed coronavirus replication in the cell culture at the micromolar range of EC50 values and had low cytotoxicity. The found chemically diverse inhibitors can be used for optimization in order to obtain a leader compound, the basis of new direct-acting antiviral drugs against the SARS-CoV-2 coronavirus.

Supercomputing, Docking and Quantum Mechanics in Quest for Inhibitors of Papain-like Protease of SARS-CoV-2

Lomonosov-2 supercomputer is used to search for new organic compounds that can suppress the replication of the SARS-CoV-2 coronavirus. The latter is responsible for the COVID-19 pandemic. Docking and a quantum-chemical semiempirical atomistic modeling method are used to find inhibitors of the SARS-CoV-2 papain-like protease, which is one of the key coronavirus enzymes responsible for its replication. The atomistic model of the papain-like protease of this coronavirus is based on the high-resolution structure deposited in the Protein Data Bank. The SOL docking program has been used for virtual screening of more than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$40000$$\end{document} low molecular weight molecules (ligands). Ligands with the highest protein-ligand binding energy, selected using the docking results, were subjected to quantum-chemical calculations. The latters are performed by the PM7 semiempirical method with the COSMO implicit solvent model using the MOPAC program. The enthalpy of protein-ligand binding is calculated for the best position of the ligand in the protein. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$19$$\end{document} ligands were selected for experimental in vitro testing as candidates for papain-like protease inhibitors base on docking and quantum-chemical results. In case of experimental confirmation, these compounds may become the basis for direct-acting antiviral drugs for the SARS-CoV-2 coronavirus.

Combined Docking with Classical Force Field and Quantum Chemical Semiempirical Method PM7

Results of the combined use of the classical force field and the recent quantum chemical PM7 method for docking are presented. Initially the gridless docking of a flexible low molecular weight ligand into the rigid target protein is performed with the energy function calculated in the MMFF94 force field with implicit water solvent in the PCM model. Among several hundred thousand local minima, which are found in the docking procedure, about eight thousand lowest energy minima are chosen and then energies of these minima are recalculated with the recent quantum chemical semiempirical PM7 method. This procedure is applied to 16 test complexes with different proteins and ligands. For almost all test complexes such energy recalculation results in the global energy minimum configuration corresponding to the ligand pose near the native ligand position in the crystalized protein-ligand complex. A significant improvement of the ligand positioning accuracy comparing with MMFF94 energy calculations is demonstrated.