» Articles » PMID: 33326239

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

Abstract

We present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

Citing Articles

Density functional theory and molecular dynamics simulation-based bioprospection of essential oil metabolites against protein tyrosine phosphatase 1B for interventive antidiabetic therapy.

Adedirin O, Abdulsalam R, Nasir-Naeem K, Oke A, Jubril A, Sabiu S Heliyon. 2025; 11(3):e42239.

PMID: 39991230 PMC: 11847251. DOI: 10.1016/j.heliyon.2025.e42239.


Toward the Prediction of Binding Events in Very Flexible, Allosteric, Multidomain Proteins.

Basciu A, Athar M, Kurt H, Neville C, Malloci G, Muredda F J Chem Inf Model. 2025; 65(4):2052-2065.

PMID: 39907634 PMC: 11863385. DOI: 10.1021/acs.jcim.4c01810.


Multi-scale computational modeling to identify novel chemical scaffolds as trehalose-6-phosphate phosphatase inhibitors to combat .

Noor S, Aljasir M, Bashir M, Khan K, Ahmad S, Abideen S In Silico Pharmacol. 2025; 13(1):21.

PMID: 39901924 PMC: 11787118. DOI: 10.1007/s40203-025-00309-5.


Innovative applications of artificial intelligence in zoonotic disease management.

Guo W, Lv C, Guo M, Zhao Q, Yin X, Zhang L Sci One Health. 2024; 2:100045.

PMID: 39077042 PMC: 11262289. DOI: 10.1016/j.soh.2023.100045.


Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms.

Ngo V Biophys J. 2024; 123(16):2623-2637.

PMID: 38946141 PMC: 11365224. DOI: 10.1016/j.bpj.2024.06.028.


References
1.
Wang F, Chen C, Tan W, Yang K, Yang H . Structure of Main Protease from Human Coronavirus NL63: Insights for Wide Spectrum Anti-Coronavirus Drug Design. Sci Rep. 2016; 6:22677. PMC: 4780191. DOI: 10.1038/srep22677. View

2.
Ton A, Gentile F, Hsing M, Ban F, Cherkasov A . Rapid Identification of Potential Inhibitors of SARS-CoV-2 Main Protease by Deep Docking of 1.3 Billion Compounds. Mol Inform. 2020; 39(8):e2000028. PMC: 7228259. DOI: 10.1002/minf.202000028. View

3.
Barca G, Bertoni C, Carrington L, Datta D, De Silva N, Deustua J . Recent developments in the general atomic and molecular electronic structure system. J Chem Phys. 2020; 152(15):154102. DOI: 10.1063/5.0005188. View

4.
Kim Y, Jedrzejczak R, Maltseva N, Wilamowski M, Endres M, Godzik A . Crystal structure of Nsp15 endoribonuclease NendoU from SARS-CoV-2. Protein Sci. 2020; 29(7):1596-1605. PMC: 7264519. DOI: 10.1002/pro.3873. View

5.
Hagemeijer M, Rottier P, de Haan C . Biogenesis and dynamics of the coronavirus replicative structures. Viruses. 2012; 4(11):3245-69. PMC: 3509692. DOI: 10.3390/v4113245. View