Identification of a potential SARS-CoV2 inhibitor via molecular dynamics simulations and amino acid decomposition analysis

Nima Razzaghi-Asl,Ahmad Ebadi,Sara Shahabipour,Danial Gholamin
DOI: https://doi.org/10.1080/07391102.2020.1797536
2020-07-24
Journal of Biomolecular Structure and Dynamics
Abstract:Considering lack of validated therapeutic drugs or vaccines against contagious SARS-CoV2, various efforts have been focused on repurposing of existing drugs or identifying new agents. In an attempt to identify new and potential SARS-CoV2 inhibitors targeting specific enzyme of the pathogen, a few induced fit models of SARS-CoV2 main protease (Mpro) including <i>N</i>-aryl amide and aryl sulfonamide based fragments were subjected to a multi-step <i>in silico</i> strategy. Sub-structure query of co-crystallographic fragments provided numerous ZINC15 driven commercially available compounds that entered molecular docking stage to find binding interactions/modes inside Mpro active site. Docking results were reevaluated through time dependent stability of top-ranked ligand-protease complexes by molecular dynamics (MD) simulations within 50 ns. Relative contribution of interacted residues in binding to the most probable binding pose was estimated through amino acid decomposition analysis in B3LYP level of theory with Def2-TZVPP split basis set. In confirmation of docking results, MD simulations revealed less perceptible torsional distortions (more stable binding mode) in binding of ZINC_252512772 (ΔG<sub>b</sub> −9.18 kcal/mol) into Mpro active site. H-bond interactions and hydrophobic contacts were determinant forces in binding interactions of <i>in silico</i> hit. Quantum chemical calculations confirmed MD results and proved the pivotal role of a conserved residue (Glu166) in making permanent hydrogen bond (98% of MD simulations time) with ZINC_252512772. Drug-like physicochemical properties as well as desirable target binding interactions nominated ZINC_252512772 as a desirable <i>in silico</i> hit for further development toward SARS-CoV2 inhibitors. A few <i>N</i>-aryl amide/aryl sulfonamide based fragments were subjected to a multi-step <i>in silico</i> strategy to afford potential SARS-CoV2 Mpro inhibitors.MD simulations revealed less perceptible torsional distortions (more stable binding mode) in binding of ZINC_252512772 (ΔG<sub>b</sub> -9.18 kcal/mol) into Mpro active site.H-bond interactions and hydrophobic contacts were determinant forces in binding interactions of <i>in silico</i> hit.Quantum chemical calculations confirmed MD results and proved pivotal role of a conserved residue (Glu166) in making permanent hydrogen bond (98% of MD simulations time) with ZINC_252512772.Communicated by Ramaswamy H. Sarma
What problem does this paper attempt to address?