COVID-19 related research projects from the CSB group

The novel coronavirus (SARS-CoV-2) that has emerged from Wuhan, China in December 2019 has spread to almost all countries in the World causing a dramatic number of deaths. The current absence of antiviral treatment against the SARS-CoV-2 urges the scientific community to accelerate the drug discovery research process.

One way to identify potential treatments and to be able to administer it swiftly is to focus on drug repurposing studies, i.e. to investigate the SARS-CoV-2 antiviral potential of drugs that have already been approved for human use.

Screening of approved drugs against the SARS-CoV-2 protease

Proteins that are crucial for the survival and replication of the virus are the most attractive targets for such studies. Here we have focused on the SARS-CoV-2 main protease (3CLpro) that plays an essential role in the virus replication process by screening ~2000 approved drugs (and 6 experimental ones - these can be identified by searching the tables using the term ‘Investigational’) against this particular protein. We used for that HADDOCK2.4, following a pharmacophore-based and a shape-based strategy adapted from our successful participation to the D3R Grand Challanges and described in the following paper:

More details of the protocol are provided below.

The results of our effort can be seen in the tables and graphs below. The top part presents the results based on top cluster rankings, while the bottom part is based on single structure rankings. The scores (in arbitrary units) correspond to the HADDOCK score calculated as:

HADDOCKscore =  1.0 Evdw + 0.1 Eelec + 1.0 Edesol


where:

• Evdw is the van der Waals intermolecular energy
• Eelec is the electrostatic intermolecular energy
• Edesol is an empirical desolvation energy term.

The plots show one data point per compound. In total there are 2020 compounds, grouped into one of four categories ‘Protease Inhibitors’, ‘Antivirals’, ‘Antiinfectives’ and ‘Other’. The first is made up compounds that are known to inhibit proteases, the second antiviral medications, the third general antiinfectives and the last everything else. The ‘NA’ group corresponds to compounds that have no specific associations - these tend to be things like dietary supplements, aminoacid residues, etc… These are not shown in the plots but are listed in the tables.

The compounds were obtained from Drugbank and pre-processed using the OpenEye Omega software (Hawkins et al. J. Chem. Inf. Model. 50 572-584 (2010)). In addition to these compounds we also obtained some of their active metaqbolites from PubChem. These can be recognised by the prefix CID. For the protease we used PDB entry 6Y2F (Zhang et al. Science, 2020).

Relevance

In total, 7 molecules of the HADDOCK Top 100 drugs (one ranked #2) are under clinical trials in Europe as reported in the European Union Drug Regulating Authorities Clinical Trials Database (EudraCT) https://www.clinicaltrialsregister.eu/ctr-search/search?query=covid-19:

• 1 kinase targeting antitumor drug: Imatinib
• 2 antiviral drugs: Lopinavir, Remdesivir
• 3 antiinfective drugs: Azithromycin, ceftriaxone, erythromycin
• 1 angiotensin II receptor antagonist: Telmisartan

Acknowledgments

Our group efforts have been led in major part by Dr. Panagiotis Koukos and Dr. Manon Réau. Dr. Ed Moret from the Department of Pharmaceutical Sciences at Utrecht University has provided valuable expertise on the compounds and their relevance.

Screening results

The plots are zoomable and clickable and groups can be hidden by clicking on them in the legend. Hovering over points in the plot will reveal the name of the compounds.

The tables are sorted based on the HADDOCK score. They are also dynamic, sortable and searchable.

Screening of approved drugs against the RNA dependent RNA polymerase

In addition the screening against the main protease we docked the virtual libray against the RNA dependent RNA polymerase (RdRp) as well. We used a truncated from of PDB entry 7BV2 for the docking and the same compound conformers as for the Mpro screening.

Screening of approved drugs against the ACE2 receptor

In addition the screening against the main protease we docked the virtual libray against the RNA dependent RNA polymerase (RdRp) as well. We used a truncated from of PDB entry 1R4L for the docking and the same compound conformers as for the Mpro screening.

Protocol

Virtual library

We considered approved drugs from the DrugBank dataset with molecular weight < 750 g/mol and > 5 heavy atoms. Some of these drugs act as prodrugs, i.e. they are inactive compounds that are metabolized into their active form after administration. When the data was available we collected the active form of the documented prodrugs, referred here as « Active Metabolites », from the PubChem database. We ended up with 1977 drugs and 50 active metabolites and we generated up to 500 conformers per compound using OpenEye.

Shape and pharmacophore-based protocols

Template compounds

The rational behind HADDOCK docking is to make use of experimental information to guide the docking. Herein, we took advantage of the large amount of high quality holo structures of the SARS-CoV-2 3CLpro and related proteins (> 90% identity) published in the Protein Data Bank. Among those crystallographic data, 66 molecules are non-covalent and covalent active-site fragments from the large XChem crystallographic fragment screen against 3CLpro performed by Diamond. In total, we collected 92 molecules targeting the 3CLpro(-related) binding site.

We identified one crystallographic template to be used for the docking of every target compound in the virtual library. For the selection, we calculated the Tanimoto coefficient computed over the Maximum Common Substructure (MCS) between target and template compounds and selected the template with the highest Tanimoto coefficient. After superimposing all templates we transformed the heavy atoms of their compounds into shape beads and defined restraints between the beads and heavy atoms of the target compound to guide it in the binding pocket.