Evaluation of SARS-CoV-2 main protease and inhibitor interactions using dihedral angle distributions and radial distribution function

Tika Ram Lamichhane,Madhav Prasad Ghimire
DOI: https://doi.org/10.1016/j.heliyon.2021.e08220
IF: 3.776
2021-10-01
Heliyon
Abstract:In order to evaluate the interactions between a potential drug candidate like inhibitor N3 and the residues in substrate binding site of SARS-CoV-2 main protease ( <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> ), we used molecular docking and dynamics simulations. The structural features describing the degrees of folding states of <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> formed by beta-barrels and alpha-helices were analyzed by means of root mean square deviation, root mean square fluctuation, radius of gyration, residue velocity, H-bonding, dihedral angle distributions and radial distribution function. All of the residues forming ligand binding domain (LBD) of <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> lie within the allowed region of the dihedral angle distributions as observed from the equilibrating best pose of <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> -N3 system. Sharp peaks of radial distribution function (RDF) for H-bonding atom pairs (about 2 Å radial distance apart) describe the strong interactions between inhibitor and SARS-CoV-2 <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> . During MD simulations, HSE163 has the lowest residue speed offering a sharp RDF peak whereas GLN192 has the highest residue speed resulting a flat RDF peak for the H-bonding atom pairs of <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> -N3 system. Along with negative values of coulombic and Lenard-Jones energies, MM/PBSA free energy of binding contributed by the non-covalent interactions between <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> and N3 has been obtained to be -19.45 ± 3.6 kcal/mol. These physical parameters demonstrate the binding nature of an inhibitor in <math> <msup><mrow><mi>M</mi></mrow> <mrow><mi>pro</mi></mrow> </msup> </math> -LBD. This study will be helpful in evaluating the drug candidates which are expected to inhibit the SARS-CoV-2 structural proteins.
What problem does this paper attempt to address?