key: cord-291248-0kuc9jv9
authors: Al-Sehemi, Abdullah G.; Pannipara, Mehboobali; Parulekar, Rishikesh S.; Patil, Omkar; Choudhari, Prafulla B.; Bhatia, M. S.; Zubaidha, P. K.; Tamboli, Yasinalli
title: Potential of NO donor furoxan as SARS-CoV-2 main protease (M(pro)) inhibitors: in silico analysis
date: 2020-07-08
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1790038
sha: 
doc_id: 291248
cord_uid: 0kuc9jv9

The sharp spurt in positive cases of novel coronavirus-19 (SARS-CoV-2) worldwide has created a big threat to human. In view to expedite new drug leads for COVID-19, Main Proteases (M(pro)) of novel Coronavirus (SARS‐CoV‐2) has emerged as a crucial target for this virus. Nitric oxide (NO) inhibits the replication cycle of SARS-CoV. Inhalation of nitric oxide is used in the treatment of severe acute respiratory syndrome. Herein, we evaluated the phenyl furoxan, a well-known exogenous NO donor to identify the possible potent inhibitors through in silico studies such as molecular docking as per target analysis for candidates bound to substrate binding pocket of SARS-COV-2 M(pro). Molecular dynamics (MD) simulations of most stable docked complexes (M(pro)-22 and M(pro)-26) helped to confirm the notable conformational stability of these docked complexes under dynamic state. Furthermore, Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations revealed energetic contributions of key residues of M(pro) in binding with potent furoxan derivatives 22, 26. In the present study to validate the molecular docking, MD simulation and MM-PBSA results, crystal structure of M(pro) bound to experimentally known inhibitor X77 was used as control and the obtained results are presented herein. We envisaged that spiro-isoquinolino-piperidine-furoxan moieties can be used as effective ligand for SARS-CoV-2 M(pro) inhibition due to the presence of key isoquinolino-piperidine skeleton with additional NO effect. Communicated by Ramaswamy H. Sarma

Novel Coronavirus (SARS-CoV-2) causing COVID-19 disease was originated from Hubei province of China, and has been exponentially increasing in whole world with United States of America, China, Europe, India and Middle East countries becoming the greatest concern. It is one of three Coronaviruses which have been associated with pneumonia, including severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle-East respiratory syndrome coronavirus (MERS-CoV). The World Health Organization announced that the outbreak of this deadly and fast-spreading infectious pneumonia constitutes a global health emergency (WHO, 2019 -nCoV Situation Report, 2020 . As on 22nd June 2020, the total numbers of cases around the world were recorded to be 9,072,803 with more than 471,178 deaths (worldometers.info/coronavirus).

WHO has been declared COVID-19 as a global pandemic disease but currently there is no effective treatment available and we only depend on preventive and supportive therapies (Cheng, 2019) . Researchers around the globe are desperately working to discover either effective vaccines or antiviral drugs, however these experimental methods for drug discovery are expensive and time-consuming. In silico based screening methods offer easy way, fast and low-cost techniques for novel testable hypotheses to discover potential lead drug candidate (Cheng, 2019) . There are many clinical trials in progress by pharmaceutical companies to afford the potential agents but very limited data with respect to in vitro and in vivo activities of the lead drugs is available in literatures which are in clinical studies for treatment of COVID-19 (Lu, 2020) .

Nitric oxide (NO) is involved in a broad range of processes and act as an important signaling molecule. Thus, the designing of NO-donors/hybrid molecules with NO releasing group may be beneficial for a variety of disorders, such as: cardiovascular, inflammatory, bacterial, fungal, parasitic, ocular diseases and cancer, arthritis, asthma, cerebral ischemia, Parkinson's disease, neurodegenerative diseases, seizures and viral (Serafim et al., 2012) . The well know NO donor, furoxan is 1,2,5-oxadiazole-2-oxide that could release NO in the presence of thiols as well as by enzymatic action in tissue with complex process that probably involves generation more than one redox form of NO (Gasco et al., 2004) .

Nitric oxide (NO) is an antiviral effector of the immune system and inhibits the replication of a variety of viruses (Saura et al., 1999) . NO can inhibit the replication of viruses that encode cysteine proteases, such as members of the Picornavirus family and the Coronavirus family (Mannick, 1995; Reiss & Komatsu, 1998) . Endogenously generated NO can inhibit the replication cycle of SARS-CoV and hence exhibit the desired antiviral effect (Åkerstr€ om et al., 2009) . NO inhibits the replication of SARS-CoV by dual mechanisms. First, NO or its derivatives cause a reduction in the palmitoylation of nascently expressed spike (S) protein which affects the fusion between the S protein and its cognate receptor i.e. angiotensin converting enzyme 2 (ACE-2). Second, NO or its derivatives cause a reduction in viral RNA production in the early steps of viral replication, and this could possibly be due to an effect on one or both of the cysteine proteases encoded in Orf1a of SARS-CoV (Åkerstr€ om et al., 2009) .

In-vitro studies have identified significant inhibition activity of NO against SARS-coronavirus infection by S-Nitroso-Nacetylpenicillamine, a nitric oxide donor (Keyaerts et al., 2004) . In patients with severe acute respiratory syndrome (SARS), inhalation of nitric oxide (NO) has improved arterial oxygenation and enabled the reduction of inspired oxygen therapy and airway pressure . Currently, clinical trial testing of inducible nitric oxide (iNO) is underway in COVID-19 infected patients complicated with Acute respiratory distress syndrome (ARDS) [NCT04306393 and NCT04305457] (clinicaltrials.gov).

Theophylline derivative and pyrimidone derivative were identified as inhibitor of RNA binding agent against 2019-nCoV N protein (N terminal domain) by in-silico studies (Sarma et al., 2020) . Many natural compounds were studied for their properties to block the activity of the angiotensin converting enzyme 2 (ACE-2) as a receptor for SARS-CoV-2 and as a potential therapeutic target of the COVID-19 virus using in silico study (Abdelli et al., 2020) . The extensive research is required by using computer simulations such as (A) Exploration of the energetic binding affinity of SARS-CoV-2 M pro with new inhibitors based on free energy calculations. (B) Investigation of the structural properties, flexibility, conformational changes of the SARS-CoV-2 M pro , and study interaction pattern between virus and membrane, virus and inhibitor. C) Monitoring of the thermodynamics (virus-water interaction) properties of SARS-CoV-2 M pro in the presence/ absence of the antiviral inhibitor (Boopathi et al., 2020) . Thus computational study could be an effective and swift approach on discovering novel uses for existing medications or approved antiviral drugs to inhibit SARS-CoV-2 M pro (Al-Khafaji et al., 2020; Khan et al., 2020a Khan et al., , 2020b Lobo-Galo et al., 2020; Muralidharan et al., 2020; Pant et al., 2020) or SARS-CoV-2 RdRp (Elfiky, 2020a) . The researchers have also studied the potential application of natural origin phytochemicals and bioactive derivatives with an anti-viral properties as inhibitors of SARS-CoV-2 M pro using molecular docking (Aanouz et al., 2020; Das et al., 2020; Elfiky, 2020b; Gyebi et al., 2020; Islam et al., 2020; Joshi et al., 2020; Sinha et al., 2020) . The advent of high resolution experimental structure of the target, the predictions of virtual screening studies and binding energy calculations are generally more accurate. The high resolution experimental structure of the target, the main protease of SARS-CoV-2 (Jin et al., 2020; Zhang et al., 2020) , was utilized in the current study as the target for molecular docking study.

Keeping in view the importance of early screening of NO donor hybrids for effective anti-viral attributes, the present study focus on evaluation of phenyl furoxan as potential drug candidates by using molecular docking, MD simulation and MM-PBSA approach. The predictions of this study will provide information that can be utilized further for development of efficient drugs by testing it in vitro, in vivo and through clinical trials to control SARS-CoV-2.

The crystal structure of the SARS-CoV-2 M pro bound to potent broad-spectrum non-covalent inhibitor X77 (PDB ID: 6W63) was selected for the docking analysis (www.rcsb.org; www.pdb.org). The protein structure was selected on the basis of resolution and method of determination. The selected protein structure was analyzed for the Ramachandran plot using V life MDS 4.6 (www.vlifesciences.com). Selected protein structure of SARS-CoV-2 M pro was refined for removal of native ligand structure X77 keeping its native conformation intact (Pettersen et al., 2004) . The refined protein structure was utilized further for docking analysis.

Molecular structures of the selected furoxan derivatives was developed using 2D molecular builder in V life MDS and further converted in to 3D geometry and optimized via application of merck molecular force field (MMFF) to get optimized ligand structures (vlifesciences.com). These optimized ligand structures were utilized for docking analysis.

The redocking protocol is utilized for the validation of the docking protocol. The redocking was carried out via docking of co crystallized ligand X77 in SARS-CoV-2 M pro . The RMSD was determined between docked ligand X77 and co crystallized ligand X77. The RMSD value less than 1.0 Å was considered to be satisfactory and which will ascertain that the binding mode of the selected furoxan derivative will be similar to that of co crystallized ligand X77.

Docking simulations was performed to predict the possible binding mode of the selected furoxan derivatives. The Grip based docking analysis followed by redocking using AutoDock Vina (Trott & Olson, 2010) which uses Lamarckian Genetic Algorithm (LGA) was performed keeping ligand structures (Table 1) in the flexible conformations. The rotational angle between the two conformations was kept at 10 and total number of rotation to 30 in grip based docking.

The best docking pose was selected on the basis of the docking score, binding energy and types of interactions.

The MD simulations study of M pro -X77 (Control) and most stable docked complexes obtained as per binding energy and docking scores i.e. M pro -22 and M pro -26 was performed by using GROMACS 5.0.5 (www.gromacs.org) program (Abraham et al., 2015) on Linux operating system to examine the structural behavior in terms of stability of M pro -22 and M pro -26 complexes in comparison to control under dynamic state. The receptor and ligand molecules from control and docked complexes were separated, and their topology files were built separately. The topology file of receptor M pro was built using optimized potentials for liquid simulations of all atoms (OPLS-AA) force field (Kaminski et al., 2001) . Whereas, topology files of ligand X77 (control), 22 and 26 were generated using an automated online PRODRG server (Van Aalten et al., 1996) . For dynamics study the simulation parameters were kept uniform and constant for all three complexes. Initially the solvation of all the complexes was carried out in the system of a cubic box with a size of 1.0 nm by using SPC (simple point-charge) water molecules to provide an aqueous environment. The systems were then neutralized by adding three Na þ ions. The periodic boundary conditions (PBC) were applied in all directions. Energy minimization of the solvated complexes was performed with the steepest descent (SD) method for 50,000 steps at 300 K to remove bad contacts until a tolerance of 1000 kJ/mol was achieved. The linear constraint solver (LINCS) algorithm was used to constrain all bond lengths (Hess et al., 1997) . A cutoff of 1.4 nm for the van der Waals interaction and 1.2 nm for electrostatic interactions were used for the simulation of all three complexes. The energy minimization of complexes was then followed by equilibration in two phases. In the first phase of equilibration, an NVT (constant number of particles, volume, and temperature) ensemble with a constant temperature of 300 K was used with a coupling constant of 0.1 picosecond for 500 picosecond, whereas in second phase an NPT (constant number of particles, pressure, and temperature) ensemble with a constant pressure of 1 bar was used with a coupling constant of 2 picosecond for 1 ns equilibration period. For both ensembles of equilibration, the coupling scheme of V-rescale was used and Particle Mesh Ewald (PME) algorithm (Essmann et al., 1995) was used to calculate long-range electrostatic interactions. After equilibration, the production MD run was performed for 100 ns to carry out structural analysis on all complexes. All resulting trajectories obtained from 100 ns MD simulations were analyzed using GROMACS utilities (Abraham et al., 2015) . The root mean square deviation (RMSD) was calculated over the 100 ns production simulation for all back-bone atoms of each complex. Similarly, root mean square fluctuations (RMSF), radius of gyration (Rg) relative to the initial structure and H bonds were also calculated for both the complexes and control to assess the stability under dynamic state. MD simulations for all complexes were performed in replicates to check the consistency and reproducibility of the results.

The binding free-energy (DG bind ) calculations for trajectories obtained from 100 ns MD simulation of complexes viz; Control, M pro -22 and M pro -26 were performed using g_mmpbsa tool which implements the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) to estimate interaction free energies (Kumari et al., 2014) . g_mmpbsa tool represents the DG bind of receptor with ligand (inhibitor) in the solvent by using following equations:

where, G complex corresponds to total free energy of the complex, G receptor and G ligand are the total energies of the only receptor and ligand respectively.

where, DG MM indicates the average molecular mechanics potential energy in vacuum and DG solv represents the free energy of solvation.

whereas, DG vdw corresponds to Van der Waals energy and DG elec represents electrostatic energy.

To calculate binding free energies, the VDW energy, electrostatic energy, molecular mechanics energy, polar solvation energy, and solvent-accessible surface area (SASA) energy were calculated for the trajectories obtained from 100 ns MD simulation of all three complexes. The contribution of each residue in binding free energy was calculated by MmPbSaDecomp.py python script (Kumari et al., 2014) . This has been helpful in determining the crucial residues and their thermodynamic state involved in interaction with the SARS-CoV-2 M proinhibitor complexes. The binding energy calculated using g_mmpbsa shows an apparent correlation of 0.80-0.85 with the experimental binding free energy (Kumari et al., 2014) .

Multiple small molecules have been reported to exhibit potent SARS-CoV M pro inhibition activity (Ghosh et al., 2020) . Moreover, Andrographolide (Enmozhi et al., 2020) and resveratrol (Wahedi et al., 2020) were found to be potential inhibitors of SARS-CoV-2 main protease. The 28 different nitric oxide donor furoxan derivatives coupled with benzhydrylpiperazine and spiro-isoquinolino-piperidine from our previous research work (Prabhuling et al., 2020; Pudukulatham et al., 2016) were selected for the current in silico analysis.

Development of small molecules with selective action on the M pro of the COVID-19 is challenging task as protease inhibition is often correlated with the peptide like structures. The development of imidazole carboxamide derivatives as potent SARS-CoV-2 M pro inhibitors open the avenue for discovery of small peptide or the peptide like structures as an attractive approach for investigation of potent inhibitors 

19 4-(4-nitrophenyl)-3-(f3-oxo-2,3-dihydro-1H-spiro [isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1,2, 5-oxadiazol-2-ium-2-olate À8.8 À43.83 GLN192 HIS41 MET49 20 4-(4-methoxyphenyl)-3-(f3-oxo-2,3-dihydro-1H-spiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1, 2,5-oxadiazol-2-ium-2-olate against SARS-CoV-2 M pro . As selected furoxan derivatives have peptide like structures and due to their structural resemblance with reported imidazole carboxamide (X77) they are virtually analyzed against SARS-CoV-2 M pro to explore their binding potential.

The results of the predicted pharmacokinetics properties of the analyzed furoxan derivatives are displayed in Figures 1-7 and Table 1 . The molecular docking revealed that the spiroisoquinolino-piperidine class of molecules have more strong binding potentials than benzhydrylpiperazine class of molecules. The docking analysis indicated that the spiro-isoquinolino-piperidine ring is responsible for the anchoring of the molecule in the binding pocket in such a way that other parts of molecule will be adopting the perpendicular orientation in the bind pocket which is responsible for the key interactions in the binding pocket of SARS-CoV-2 M pro . In case of the benzhydrylpiperazine class of molecules the 4-(2,4-dichlorophenyl)-3-(f3-oxo-2,3-dihydro-1Hspiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1, 2,5-oxadiazol-2-ium-2-olate

planner ring orientation of the piperazine pushing interactive phenyl ring in the lower plane and two phenyl ring in the co planner conformation which might be responsible for steric repulsion and subsequently misfitting in the binding pocket of SARS-CoV-2 M pro . The crystal structure of the SARS-CoV-2 M pro bound to potent broad-spectrum non-covalent inhibitor X77 (PDB ID: 6W63) was selected for the docking analysis. The resolution of the selected SARS-CoV-2 M pro was found to be 2.10 Å and as it is X-ray determined structure, which was found to be suitable for docking process. Ramachandran plot of the selected protein structure was determined to ascertain the structural geometry and to validate the selected protein structure. The Core region of the protein was found to be contain 84.16%, allowed region was found to contain 11.22% of amino acids, which indicated selected protein structure is in good shape and having excellent structural geometry (Figure 1 ).

Redocking analysis was performed to validate the docking protocol. The RMSD between docked non-covalent inhibitor X77 and co-crystalized X77 was found to be 0.11 Å which validated the docking protocol. Results of the redocking process indicated that using applied docking protocol to the selected furoxan derivatives will give similar type of binding potential to that of co-crystalized X77 (Figure 2) . Docking analysis of the studied furoxan derivatives were found to have similar binding ability to SARS-CoV-2 M pro as compared to reported inhibitors. Derivative 4-(3-nitrophenyl)-3-(f3-oxo-2,3-dihydro-1H-spiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1,2,5-oxadiazol-2-ium-2-olate (22) have binding affinity of À9.6 kcal/mol, and docking score of À90.93, with important hydrogen binding interactions with CYS145 (2.208 A 0 ) and SER144 (2.164 A 0 ) and aromatic interactions with HIS41 (5.281 A 0 ) and HIS 163 (4.628 A 0 ) as shown in Figure 3 .

Derivative 4-(3-bromophenyl)-3-(f3-oxo-2,3-dihydro-1Hspiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1,2,5-oxadiazol-2-ium-2-olate (26) shown binding affinity of À9.8 kcal/ mol and docking score of À90.91 as compared to 22.

Hydrogen binding interactions with CYS145 (2.153 A 0 ) and SER144 (1.823 A 0 ) and aromatic interactions with HIS41 (5.175 A 0 ) and HIS 163(4.760 A 0 ) as shown in Figure 4 .

Derivative 4-(4-fluorophenyl)-3-(f3-oxo-2,3-dihydro-1H-spiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1,2,5-oxadiazol-2-ium-2-olate (23) was shown binding affinity of À8.3 kcal/ mol and docking score of À89.44. Important hydrogen binding interactions was observed with GLN189 (1.926 A 0 ), GLY143 (1.565 A 0 ) and aromatic interactions with HIS41(4.521 A 0 ) and HIS 163(4.994 A 0 ) as shown in Figure 5 .

Derivative 4-(4-methoxyphenyl)-3-(f3-oxo-2,3-dihydro-1Hspiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1,2,5-oxadiazol-2-ium-2-olate (20) displays binding affinity of À8.4 kcal/mol and docking score of À88.64. Hydrogen bond interactions with GLY 143 (1.507 A 0 ), SER144 (2.533 A 0 ) and aromatic interactions with HIS41 (4.323 A 0 ) as shown in Figure 6 . 4-(2-chlorophenyl)-3-(f3-oxo-2,3-dihydro-1H-spiro[isoquinoline-4,4 0 -piperidine]-1 0 -ylgmethyl)-1,2,5-oxadiazol-2-ium-2olate (24) showed binding affinity of À8.1 kcal/mol and docking score À84.66, along with interaction with SER144

(2.364 A 0 ) via formation of hydrogen bond and with HIS163 (4.797 A 0 ) via formation of aromatic interaction as shown in Figure 7 .

The MD simulation study of most stable docked complexes viz: 2018a, 2018b). All systems are subjected to 100 ns MD simulations and their stability in dynamic state was investigated by using GROMACS utilities for calculating RMSD, RMSF, Rg and H bonds (Abraham et al., 2015) . The RMSD is a significant parameter used to analyze the stability as per output obtained from MD trajectories which is used to verify the convergence of simulations. RMSD of the protein backbone atoms are plotted as a function of time to check the stability of each complex throughout the simulation. From Figure 8A it is clear that control system reached stable equilibrium after 10 ns whereas, both M pro -22 and M pro -26 systems reached equilibrium after 15 ns. The average RMSDs for all simulated systems were calculated and shown in Table 2 . These average RMSDs vary between 0.18 nm to 0.28 nm for all systems ( Table 2 ). The RMSD values of Control and M pro -22 complex were found to be in the range of 0.1 nm to 0.2 nm and remained stable over entire simulation time after reaching equilibrium ( Figure 8A ). The RMSD values of M pro -26 complexes are in the range of 0.2 nm to 0.3 nm with slight higher deviation and remained stable throughout 100 ns simulation after attaining equilibrium ( Figure 8A ). The RMSF parameter in MD simulation is used as a mean describing flexibility differences among residues. The higher value of RMSF generally indicates more flexibility, whereas the low RMSF value indicates limited movements of residues during simulation in relation to its average position. Hence, we analyzed fluctuations of Ca atoms from both docked complexes along with control to study the flexibility of the residues from respective structures in presence of furoxan derivatives 22, 26 and X77 (control). The RMSF values of both the docked complexes along with control are shown in Figure-8B . These RMSF results clearly depicts that residues of M pro in the presence of 22 and 26 showed similar type of behavior as in presence of X77 (control) [ Figure 8B ]. Precisely, the active site residues of M pro i.e. CYS145, SER144, HIS41, HIS163 and GLU166 involved in binding with 22 and 26 showed significantly less fluctuations indicating the stable hydrogen and aromatics interactions of these residues with both the inhibitors as observed for X77 inhibitor used as control (Table 1, Figure 8B ). We also determined Rg and calculated average values as reported in Table 2 to understand the level of compactness in the structure of M pro in presence of furoxan derivatives 22, 26 and X77 (control). Interestingly, the Rg values of both the complexes and control were found to be in the range of 2.2 nm to 2.3 nm and were consistent over entire 100 ns simulation time scale ( Figure  8C ). Thus, this suggests the similarity in compactness of receptor protein M pro in the presence of 22, 26 and X77 (control) [ Figure 8C ]. Additionally, the effective binding ability of furoxan derivatives 22, 26 with SARS-CoV-2 M pro protein was illustrated by H bonds analysis (Figure 9 ). Both inhibitors 22, 26 and more prominently 26 showed strong and consistent hydrogen bonding with M pro in comparison to X77 (control) over entire 100 ns simulation run (Figure 9 ). Altogether, MD simulation results for both the complexes were in accordance with control system which significantly suggests the effective binding ability of inhibitors 22, 26 as like X77 (control) towards SARS-CoV-2 M pro receptor structure ( Figure 8A -C and Figure 9 ).

The binding free energy calculations on the basis of different energetic terms for control, M pro -22 and M pro -26 complexes confirmed the binding efficiency of X77 (control) and furoxan derivatives 22, 26 at the known inhibitor binding site region of SARS-CoV-2 M pro structure (Table 3) . As seen from Table 3 control showed involvement of higher VDW energies for stabilizing the M pro -X77 complex whereas, M pro -22 and M pro -26 complexes showed significant involvement of higher electrostatic energies for binding thus suggesting prominent role of Figure 9 ). The control system seems to be more stabilized by non-bonded contacts in comparison to hydrogen bonds which could be the reason for higher VDW energies and lower electrostatic energies observed for control (Table 3, Figure 9 ). Similarly, polar salvation energy values were also obtained higher for M pro -22 and M pro -26 complexes in comparison to control (Table 3 ). The estimated binding energies (DG bind ) for control was obtained as À130.014 ± 9.977 whereas, DG bind for M pro -22 and M pro -26 was observed to be À183.860 ± 9.106 and À171.972 ± 8.757 respectively (Table 3) . Hence, on comparative account 22, 26 showed more effective binding than X77 (control) towards SARS-CoV-2 M pro receptor on different energetic terms (Table 3) .

Furthermore, individual energy decomposition of all residues from SARS-CoV-2 M pro receptor structure was calculated to investigate the key residue of M pro involved in binding with X77 (control) and 22, 26 ( Figure 10A-C) . The results obtained from the calculation of energetic contribution of individual residues clearly confirm the effective binding ability of 22, 26 in comparison to X77 (control) towards M pro ( Figure 10A-C) . These results are totally in accordance with the results obtained for the binding free energy calculations of all three complexes on the basis of different energetic terms (Table 3, Figure 10A -C). As seen from Figure 10A -C less residue wise energy contribution is observed for control whereas, significant high values (especially electrostatics and polar solvation) are observed for complexes M pro -22 and M pro -26 thus reflecting the strong binding nature of 22, 26 towards SARS-CoV-2 M pro structure. Residue wise decomposition results confirms that the interacting residues from binding pocket of SARS-CoV-2 M pro as obtained from docking analysis (Table 1 ) also shows significantly higher energetic contribution in binding with potent furoxan derivatives 22, 26 in comparison to control ( Figure 10A -C).

Thus, overall in silico analysis performed in the current study reveal the effective binding affinity of investigated spiro-isoquinolino-piperidine-furoxan derivatives in comparison to the experimentally known inhibitor of M pro i.e. X77 (control) and also they were found to be superior to the corresponding benzhydrylpiperazine derivatives. Further modification of these spiro-isoquinolino-piperidine-furoxan derivatives can lead to discovery of potent SARS-CoV-2 M pro inhibitors having the advantage of exogenous NO release feature. Although the molecules have shown prominent results in the in silico analysis, the wet lab experiment are required to validate the results. Present results can act as starting point for development of potent selective SARS-CoV-2 M pro inhibitors.

Moroccan medicinal plants as inhibitors against SARS-CoV-2 main protease: Computational investigations

In silico study the inhibition of Angiotensin converting enzyme 2 receptor of COVID-19 by Ammoides verticillata components harvested from western Algeria

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX, 1-2, 19-25

Nitric oxide inhibits the replication cycle of severe acute respiratory syndrome coronavirus

Dual effect of nitric oxide on SARS-CoV replication: Viral RNA production and palmitoylation of the S protein are affected

Using Integrated Computational Approaches to Identify Safe and Rapid Treatment for SARS-CoV-2

Novel 2019 coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment

Inhalation of nitric oxide in the treatment of severe acute respiratory syndrome: A rescue trial in Beijing

In silico oncology drug repositioning and polypharmacology

An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective

Natural products may interfere with SARSCoV-2 attachment to the host cell

Andrographolide as a potential inhibitor of SARS-CoV-2 main protease: An in silico approach

A smooth particle mesh Ewald method

NO donors: Focus on furoxan derivatives

Drug Development and Medicinal Chemistry Efforts Toward SARS-Coronavirus and Covid-19 Therapeutics

Potential Inhibitors of Coronavirus 3-Chymotrypsin-Like Protease (3CLpro): An in silico screening of Alkaloids and Terpenoids from African medicinal plants

LINCS: A linear constraint solver for molecular simulations

A Molecular Modeling Approach to Identify Effective Antiviral Phytochemicals against the Main Protease of SARS-CoV-2

Structure of Mpro from COVID-19 virus and discovery of its inhibitors

Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease

Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides

Inhibition of SARS-coronavirus infection in vitro by Snitroso-N-acetylpenicillamine, a nitric oxide donor compound

Targeting SARS-CoV-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2'-O-ribose methyltransferase

Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach

g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations

FDA-approved thiol-reacting drugs that potentially bind into the SARS-CoV-2 main protease, essential for viral replication

Drug treatment options for the 2019-new coronavirus (2019-nCoV)

The antiviral role of nitric oxide

Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 protease against COVID-19

Peptidelike and small-molecule inhibitors against Covid-19

Molecular modeling studies to explore the binding affinity of virtually screened inhibitor toward different aminoglycoside kinases from diverse MDR strains

Insights into the antibiotic resistance and inhibition mechanism of aminoglycoside phosphotransferase from Bacillus cereus: In silico and in vitro perspective

UCSF Chimera-a visualization system for exploratory research and analysis

Synthesis and Modeling Studies of Furoxan Coupled Spiro-Isoquinolino Piperidine Derivatives as NO Releasing PDE 5 Inhibitors

Synthesis and characterization of a disubstitutedpiperazine derivative with T-type channel blocking action and analgesic properties

Does nitric oxide play a critical role in viral infections

In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain)

An antiviral mechanism of nitric oxide: Inhibition of a viral protease

Nitric oxide: State of the art in drug design

An in-silico evaluation of different Saikosaponins for their potency against SARS-CoV-2 using NSP15 and fusion spike glycoprotein as targets

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

PRODRG, a program for generating molecular topologies and unique molecular descriptors from coordinates of small molecules

Stilbene-based natural compounds as promising drug candidates against COVID-19

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved a-ketoamide inhibitors

The authors are thankful to V life Sciences, Pune for providing trial version of software for research work. The authors are thankful to the Institute of research and consulting studies at King Khalid University for funding this research through grant number 3-N-20/21 and the support of Research center for advanced materials Science is highly acknowledged.

There is no conflict of interest.