key: cord-324166-6ydn2bvy
authors: Kumar, Neeraj; Awasthi, Amardeep; Kumari, Anchala; Sood, Damini; Jain, Pallavi; Singh, Taru; Sharma, Neera; Grover, Abhinav; Chandra, Ramesh
title: Antitussive noscapine and antiviral drug conjugates as arsenal against COVID-19: a comprehensive chemoinformatics analysis
date: 2020-08-20
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1808072
sha: 
doc_id: 324166
cord_uid: 6ydn2bvy

Coronavirus pandemic has caused a vast number of deaths worldwide. Thus creating an urgent need to develop effective counteragents against novel coronavirus disease (COVID-19). Many antiviral drugs have been repurposed for treatment but implicated minimal recovery, which further advanced the need for clearer insights and innovation to derive effective therapeutics. Strategically, Noscapine, an approved antitussive drug with positive effects on lung linings may show favorable outcomes synergistically with antiviral drugs in trials. Hence, we have theoretically examined the combinatorial drug therapy by culminating the existing experimental results with in silico analyses. We employed the antitussive noscapine in conjugation with antiviral drugs (Chloroquine, Umifenovir, Hydroxychloroquine, Favlplravir and Galidesivir). We found that Noscapine-Hydroxychloroquine (Nos-Hcq) conjugate has strong binding affinity for the main protease (Mpro) of SARS-CoV-2, which performs key biological function in virus infection and progression. Nos-Hcq was analyzed through molecular dynamics simulation. The MD simulation for 100 ns affirmed the stable binding of conjugation unprecedentedly through RMSD and radius of gyration plots along with critical reaction coordinate binding free energy profile. Also, dynamical residue cross-correlation map with principal component analysis depicted the stable binding of Nos-Hcq conjugate to Mpro domains with optimal secondary structure statistics of complex dynamics. Also, we reveal the drugs with stable binding to major domains of Mpro can significantly improve the work profile of reaction coordinates, drug accession and inhibitory regulation of Mpro. The designed combinatorial therapy paves way for further prioritized in vitro and in vivo investigations for drug with robust binding against Mpro of SARS-CoV-2.

The recent outbreak of the novel Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is unprecedented. The World Health Organization (WHO) has given the guidelines that the current pandemic might go long, and the virus may become endemic in our society. As per the statistics of WHO, seven million people (as of 8 th June 2020) have already been infected with coronavirus disease-2019 and expected to increase exponentially in some of the regions (WHO, 2020) . The coronavirus is a single-stranded RNA virus that typically affects the vertebrates and causes lethal effects (Fung & Liu, 2019) . The transmission of the SARS-CoV-2 virus between different host species is generally through respiratory droplets, in aerosol form, and fomites . Due to the genetic proximity of this novel Betacoronavirus, it is similar to the previous outbreaks of severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) (Zhu et al., 2020) . The diverse spike protein of SARS-CoV-2 makes it more lethal than previous SARS and MERS family viruses. Among the various routes for its treatment/ prevention by the genetic code information of Coronavirus, reported by the Chinese research group in March 2020 (Ren et al., 2020) provided an early opportunity for scientists across the globe to develop the vaccine for the COVID-19 disease at earliest. However, it is too early to assume success for the COVID-19 vaccine because the development of vaccination is tricky; it may take time to manufacture an effective and safe vaccine on a large scale. Alternatively, to the current situation, drug repurposing seems to be an effective drug discovery strategy by utilizing the existing drugs and natural plant products for the treatment of the SARS-CoV-2 outbreak. Drug repurposing strategy can potentially shorten the time and as well as reduce the cost compared to develop novel drug discovery and their clinical trials (Cheng, 2019; Cheng et al., 2016 Cheng et al., , 2017 , . Overall, the world research community is repurposing known FDA approved drugs and improving them by conjugation or synthesizing potential analogs. Remdesivir (Beigel et al., 2020) , Hydroxychloroquine (Colson et al., 2020) , Chloroquine , Favipiravir (T-705) (Coomes & Haghbayan, 2020) , and Galidesivir (BCX4430) (Li & De Clercq, 2020) are some of the most successful repurposed small molecule drugs which has shown broad-spectrum activity against Coronavirus. These Repurposed drugs had initially shown good results; however, there were some major limitations as well, likewise potency, effectiveness, specificity to viral targets, and in some cases, side effects due to non-specificity, etc. Interestingly, combination therapy of approved drugs suggested for being an effective strategy to treat corona infection. Recently, Hydroxychloroquine (HCQ), a well-known drug for malaria, has been repurposed, which brought attention to potentially helping patients to recover from COVID-19. Also, it has been used in the combination therapy and reaches to the trials such as; Hydroxychloroquine with Mivermectin (Patr ı & Fabbrocini, 2020), Hydroxychloroquine with Azithromycin (Magagnoli et al., 2020) , Hydroxychloroquine with Nitazoxanide (Amawi et al., 2020) , etc. However, these have not shown expected outcomes and associated with issues like high mortality rates and increased frequency of irregular heartbeat. As per WHO reports, these HCQ based combination therapies took major setbacks after not performing significantly in clinical trials. Hence in present work, we propose a drug conjugate of Noscapine along with major antiviral drugs in repurposing trials (Chloroquine, Hydroxychloroquine, Umifenovir, Favlplravir, and Galidesivir) . Noscapine is an approved antitussive drug for the treatment of cough through simulating the lung linings. It is known to inhibit bradykinin enhanced cough response (Ebrahimi et al., 2003) . We are aware of the fact that the primary target for SARS-CoV is the lungs, and noscapine has been widely used to inhibit bradykinin, eventually inhibiting inflammation and lung damage (Ebrahimi, 2020) . Hence it makes it sense to conjugate Noscapine with antiviral drugs, which can decrease the severe symptoms due to the lung damage; this can overcome the setbacks of previous combinations therapies, including HCQ conjugations, by attenuating the high mortality rates due to COVID-19. Noscapine is known to obtain from the naturally occurring Opium Poppy (Papaver Sominiferum) medicinal plant, and it is characteristically a non-toxic phthalide-isoquinoline alkaloid (Aneja et al., 2006) .

We aimed to design the combinatorial therapy against the main protease (Mpro or 3CLpro) enzyme, a key regulator of the life cycle of the coronavirus. Mpro is essential for the replication of the virus as it processes the polyproteins that are translated from the viral RNA to yield functional viral proteins (Chen et al., 2020) . This makes the Mpro enzyme a potential drug target to combat coronaviruses. This fact has led to some of the research groups to develop a potential inhibitor to the protease (Mengist et al., 2020; Zhang et al., 2020) , . These inhibitors prevent the replication of the virus after entering into the host cells. Mpro protease inhibitors are reported to modify the polypeptides essential for SARS-COV-2 infection and progression in the host. Also, our group has recently reported the potent interaction of Noscapine analogs and subsequent inhibition of the Mpro protease of the SARS-CoV-2 (Kumar, Kumari, et al., 2020) .

With detailed literature, we designed and analyzed the noscapine conjugates with antiviral drugs (chloroquine, hydroxychloroquine, favipiravir, and galidesivir that are currently being clinically tested (Table 1 ). The molecular binding affinity of conjugates to inhibit the Mpro has been evaluated using high throughput computational assays. The molecular dynamics studies have been widely used in modern drug discovery and development. Computer-aided drug design has been employed to annotate the molecular binding mechanism of drugs to their specific targets (Singh et al., 2019; , . In many studies, advanced molecular docking approach has been harnessed to screen an extensive library, define the binding affinity of drugs and to elucidate the involved molecular interactions in the complexation of the ligand with receptor protein Chaudhary et al., 2020) , . Molecular dynamics (MD) analysis helps to insight into the binding mechanism of the drugs through the Root-mean-square deviation (RMSD), radius of gyration binding energy calculations, and relevant calculations (Kumar, Sood, Sharma, et al., 2020; Kumar, Sood, Tomar et al., 2020; Singh et al., 2016; ., . With employing these advanced computational approaches, deigned noscapine conjugations with antiviral drugs were examined to derive the highly efficient combination against the Mpro of SARS-CoV-2 for further validation with experimental lab assays.

2.1. Drug target main protease (mpro) of SARS-CoV-2 retrieval and structural assessments

The main protease regulatory enzyme of SARS-CoV-2 has been reported to play an essential role in coronavirus infection and progression in the host (Xue et al., 2008) . To design the drug therapy against Mpro, a three-dimensional crystal structure of Mpro was retrieved from the protein data bank (PDB ID 6LU7 

We have designed the Noscapine based antiviral cognates with antiviral drugs in trials (Chloroquine, Umifenovir, Hydroxychloroquine, Favlplravir, and Galidesivir) . Noscapine is reported possess the appeased effects on lung and respiratory infection with its antitussive properties (Kim et al., 2019) . Noscapine conjugates were drawn using the Chemdraw software and assessed for their compatibility and deformation.

The molecular interaction analyses of all conjugates were performed with the Mpro to illustrate their binding affinity. To perform the docking calculations, the target Mpro was prepared and optimized by the Whatif server (Vriend, 1990) . Mpro was confined to add missing hydrogen atom, amino acid chain, any other unwanted chain and removing the water molecules. Similarly, ligand files of noscapine-antiviral conjugates were prepared and optimized for any missing atom and chain. The molecular docking was performed using the Hex 8.0 fast Fourier algorithm-based module (Macindoe et al., 2010) . By incorporating the binding sites, grid box in all axes to Mpro was formed and ligand files were inserted to the server. The receptor grid was set to 0.6 Å for X, Y, and Z coordinates of Mpro structure. The translational step, twist range of 360 and protein flip ranges were set as per the molecular docking manual.

Moreover, to strengthen our results, we redocked the high-rank noscapine conjugate with another molecular docking tool with the Swissdock server. Swissdock works on the basis of the EADock DSS algorithm with curated scripts of receptors and ligands (Grosdidier et al., 2011) . From the obtained various docked conformation of Mpro-conjugates, the model with the lowest energy was analyzed for the binding mechanism by protein-ligand profiles (Laskowski & Swindells, 2011; Salentin et al., 2015) , . The model with the lowest docked energy score in all drug conjugates was selected for molecular dynamics simulation analysis.

To mechanistic insights into the binding of the conjugate (Nos-Hcq), the molecular dynamics simulation analysis was performed. MD simulations were implemented through GROMACS v5.0 under the force field GROMOS96 54a7 having water model SPC216 along with the time step of 1 fs for 100 ns . We investigated the atomic motions, conformational changes, structural stability, and binding energy contributing to the binding of the ligand to the target molecule. In an explicit water solvent, internal atomic motion simulation was initiated. The binding complex of Nos-Hcq conjugate was solvated in the octahedron system and stretched to 10 Å in all directions. The complex operation was minimized through four steps, and all atoms were fixed. In the minimization process 2000 steepest descent steps were executed, and all atoms were relaxed at 300 K with one atmospheric pressure in boundary states. NPT ensembles, along with periodic boundary conditions, were utilized to carry out MD simulations. A cut-off of about 12 Å was used in order to manage the Vander Waals forces. The Particle Mesh Ewald model having a cut-off of 14 Å was further utilized to calculate the electrostatic interactions. The obtained trajectory for 100 ns MD simulation run was analyzed for root mean square deviation (RMSD), root mean square fluctuations (RMSF), hydrogen bond analyses, solvent accessible surface area analysis (SASA) and radius of gyration (Rg) plot analysis (Grant et al., 2006) . Also, the native contacts, along with the principal component analysis (PCA), cross-correlation map analyses were performed throughout the simulation run (McGibbon et al., 2015) . Moreover, the binding energy calculation lead conjugate with Mpro was investigated through the reaction coordinates analysis. It determined the ligand compatibility with threshold energy and the progression of ligand for stabilized binding with the target receptor. A total of 1000 frames of interactions complex trajectory was analyzed for binding energy computation. The equation employed for the calculation is

DGcomplex, DGreceptor, DGligand represents the binding energy of complex, binding energy receptor, binding energy of ligand.

3.1. Sars CoV-2 mpro enzyme as a potential target SARS coronavirus is an enveloped positive-stranded RNA virus, which majorly affects the respiratory system and enteric system (Graham et al., 2013) . SARS-CoV-2 infection and progression in the host are regulated by multiple structural proteins. Among various structural proteins, the main protease enzyme of coronavirus is reported to play significant roles in viral replication through proteolytic machinery and involved in transcription, translation, and amplification of viral proteins (Paules et al., 2020) . Mpro enzyme is of size 306 amino acids and possesses high similarity with protease enzyme from different human and animal by sequence analysis. The protein annotation results depicted Mpro has various roles in coronavirus viz; a nucleic acid-binding domain of SARS coronavirus, polyprotein cleavage domain, coronavirus endopeptidase, corona NSP16, NSP7 and NSP9 replication implication regulatory roles, reported through Pfam annotation database.

Moreover, 3 D structural analysis of Mpro by the CATH server showed it consisted of beta-barrel and alpha domain (Figure 1 ). The secondary structural analysis by the DSSP server showed it comprised of 24% alpha-helices (10 helices, 75 residues) and 27% beta-barrel sheets (19 strands). These significant analyses along with the literature survey, suggested the Mpro enzyme as a potential target to design the COVID-19 therapeutic drugs.

Moreover, prior to molecular binding analyses the 3 D structure of the Mpro enzyme was assessed for structural conformations and physicochemical by Ramachandran plot and Saves server. The Ramachandran plot analyses depicted 90.6% residues are in the favorable region, and 9.1% residues and 0.4% residues are in the outlier region of the plot. The local quality assessment showed good structural quality factors with high homology (0.1 to 0.7 Å) with the native structures in database of 160,954 known protein structures.

With Noscapine significant roles in lung linings and dilation and Coronavirus movements into the lung cells evokes the cytokine storm and leads to neutrophilic infiltration and other implication (Martin & Ernst, 2003) . Hence, it was essential to investigate the effectiveness of antitussive noscapine in conjugation with antiviral drugs against Mpro. We chose and designed the noscapine conjugation with potential antiviral drugs (Chloroquine, Umifenovir, Hydroxychloroquine, Favlplravir and Galidesivir) using the Chemdraw and simultaneously assessed their applicability and conformations. For molecular docking, both receptor Mpro and ligand files of noscapine conjugates were prepared for molecular docking using the Whatif prepdock module and ligprep server. For performing the molecular docking, binding site of Mpro enzyme-N3 inhibitor crystal structure, recently submitted to the protein data bank, was analyzed for binding sites and insights into the molecular mechanism of inhibition by drugs . Similar binding sites were employed for performing the molecular docking of the noscapine conjugations with the target Mpro of coronavirus using the Hex 8.0 and SwissDock servers. The molecular docking was performed using Hex 8.0, which works based on a fast Fourier transform by SPF shape-density correlations, and results were recalculated with electrostatic association additionally. First, we have performed the Molecular docking analysis of drug alones which depicted the high binding score for Hydroxychloroquine with À301.08 kJ/mol, among all drugs. Further we assessed their drug likeliness properties of all drugs, none of them violated Lipinski rule of five except Galidesivir which has six hydrogen donor groups (Table 2) . Taking these results altogether, we performed the docking for all conjugates. Interestingly, obtained results showed the strong binding of noscapine conjugates with Mpro of coronavirus. The Nos-Hcq conjugate depicted the most substantial binding with an energy score of À410.09 kJ/mol, among all conjugates. Besides, to assess the effectivity of noscapine (singlet state) against the Mpro, it was also docked using the same grid, which showed the binding score of À292.42 kJ/mol, larger than the known molecule N3 molecule bound to Mpro with binding score of À260.78 kJ/mol. It can be seen from the molecular binding analyses that hydroxychloroquine depicted the high binding score in drug alone analyses, whose binding efficiency was enhanced by conjugation of noscapine by about double with a docking score of À410.09 kJ/mol to Mpro (Table 3) . Moreover, lead conjugation of Nos-Hcq does not found to follow the drug likeliness Lipinski rule of five for high molecular weight and hydrogen bond acceptor groups with combinations of two drugs. In addition, to strengthen our results, we redocked the all conjugates with Mpro using the SwissDock module. The outcomes affirmed the highest efficacy of Nos-Hcq conjugate with binding DG energy score of À10.01 kcal/mol and combined full fitness score of À1025.72 kcal/mol, among all noscapine conjugates (Figure 2 ). With these significant results, it can be attributed that Nos-Hcq conjugate has a high potential to bind the target Mpro enzyme and further can be used as effective therapeutics for SARS-CoV-2.

The detailed characterization of contacts of noscapine conjugates is vital to recognize the interaction pattern of ligands for stringent binding to target Mpro. The protein-ligand profiler was employed to determine the non-covalent contacts of complex systems. The Mpro structural domain analysis by the CATH server showed it forms majorly three ligand-binding domains and one additional domain. The globular structure analysis by Chimera modeling application suggested the plausibly first domain from Ser1-Pro96, second binding remained confined to domain-2 and domain-3 ( Figure 4 ). From these contact analyses of all conjugates, it can be attributed that Nos-Hcq conjugate has high affinity and wide molecular contacts for major domains of Mpro, among all conjugates. Moreover, it was further studied through molecular dynamics to get insight into inhibitory patterns and efficacy to combat the Mpro enzyme of coronavirus.

The binding affinity of the Nos-Hcq conjugate was investigated by the molecular dynamics simulation analyses for 100 ns. The stable binding of Nos-Hcq along with protein dynamics in the target Mpro was assessed through the GROMOS96 54a7 force field with water model SPC216 consisting time step of 1 fs. The obtained trajectory was solvated in the octahedron box in an explicit water state to provide sustainable conformations. The system minimization was done for 2000 steps by the steepest descent program and further 300 steps by conjugated gradient method through inserted force fields. The root means square deviation (RMSD), root means square fluctuation (RMSF), and Radius of gyration (Rg) were computed for the resulting trajectory to assess the stable binding of lead combination (Nos-Hcq) with Mpro of coronavirus. Moreover, to compare and determine the boosting power and synergetic contribution, MD simulation was also performed for the singlet state of Mpro. The steep RMSD variation can be considered an implication of a malleable and instinctive protein. The RMSD (average backbone range of complex system) for the frames of native contacts with target Mpro depicted the stable binding of lead Nos-Hcq conjugate to Mpro with fluctuation in the range of 0.2-0.55 nm and for Mpro singlet state fluctuation in the range of 0.15-0.81 nm ( Figure 5) . Interestingly, RMSD of the complex system was more stabilized than the native form of target receptor Mpro. Moreover, the secondary structure analysis showed the compelling flexibility in terminal residues and other loops of beta linker segments of the Mpro (residues 189-2998) to about 0.21 nm and Mpro-conjugate (residues 130-140 and 189-299) to about 0.1-0.2 nm. Also, we found the residues 130 to 299 have comparatively lesser discrepancies in RMSF statistics, which are corresponding to the interacting segment of the Mpro-conjugate complex after MD simulation. Interestingly, the binding region for lead conjugation before MD simulation was also within the residues 137 to 298, which remains the same after simulation, suggesting a strong binding of Nos-Hcq towards the Mpro and provides stability to the Mpro-conjugate. After that, the obtained trajectory was computed to evaluate the Radius of gyration of the Mpro and compared with the noscapine singly. The lead combination with the Mpro protein showed the Rg fluctuation in the range of 2.12-2.26 nm and for Mpro in range of 2.14-2.33 nm and demonstrated the stable and robust binding of Nos-Hcq with Mpro ( Figure 5) .

Furthermore, similar outcomes were obtained by SASA computation, which represented the solvent-accessible protein surface, and various orientations through the folding and proceeding to alterations in buried and exposed regions of the surface area of the protein complex. Similar observations were determined through SASA analysis representing the solvent defined protein surface and its orientation through folding, making the alterations in the exposed and remains same as before till 100 ns ( Figure 5 ). Here, Mpro-NosHcq solvation profile shows a convincing SASA value suggesting a stable structure and sturdy binding. Also, the secondary structure analysis of complex trajectory and native Mpro trajectory was performed by employing the DSSP tool of Gromacs. The cluster analyses of both depicted the conformational changes before and after the simulation run. Mpro-NosHcq and Mpro were consisted of mainly conserved b-sheets along with loops of b linker segments and other random coil regions infused with the various bend, a-helices, turn and b bridges. The observations of both plot analyses showed the Mpro rationale unfolded and implemented surface with b-sheets to bind stably with lead conjugation Nos-Hcq. Binding energy landscape MM/PB (GB) surface area (SA) calculations were performed to estimate the lead conjugation binding affinities to Mpro protein. The system was found to be the most favorable in terms of the sum of various binding energy computations. MM/GBSA free binding energy calculations were executed with two fifty frames dynamically equilibrated for the stabilized trajectory of the Mpro-NosHcq complex system between 75 ns to 80 ns. The electrostatic energy distribution in molecular mechanics was obtained through the high Vander Waal energy (VDW) À203.87 ± 15.81 kJ/mol, electrostatic energy À1077.91 ± 171.34 kJ/mol and total binding energy of À1281.78 ± 179.25 kJ/mol, which are considered as most significant contributor in the binding of ligand (Nos-Hcq conjugation) to the target protein. The polar solvation energy (PBSOL) was calculated to be À8915.14 ± 210.92 kJ/ mol, which added a vital contribution to the stable binding of the ligand to target protein.

Further, the stability and flexibility of the binding of Nos-Hcq with Mpro were strengthening with the stabilized run of trajectory with set NVT parameters. The average temperature for both native Mpro and in conjugation was close to 300 Kelvin; average pressure was found to be close to 1.4 bars and an average density of 993.5 kg/m3 ( Figure 6 ). These outcomes through 100 ns MD simulation suggested the stable binding of lead combination with significant contribution of binding energy to the target Mpro of the coronavirus. Also, the energy work profile of ligand binding and Mpro protein dynamics was investigated by SMD trajectory based on Jarzynski relationships in molecular dynamics simulation. The reaction coordinate work profile demonstrated the contribution of binding energies in potential interaction. The calculated work profile for ligand binding through the reaction coordinates stated the stabilized coupling with energy contribution reached a free energy profile of 100 kcal/mol at an RC value of 33 Å and sustained throughout binding with the less flexible atomic motions. In the early stages, it reached to the maximum value and stabilized at 100 kcal/mol with hydrogen bonds and reached to work profile of 130 kcal/mol at RC 70 Å in later stages. The atomic energy profiles for Mpro-Nos-Hcq reaction coordinates showed an average radius of 1.88 Å, with overall Priority of conjugated reaction progression about 91.84%. Along with the binding energy calculations, hydrogen bond consistency and dynamic contribution of bound Nos-Hcq revealed the involvement of a high number of hydrogen bonds throughout the simulation run for both native Mpro and complex system ( Figure 6 ). Mpro formed an average of 550 hydrogen bonds with the lead combination during the simulation and played an essential role in the selectivity of intermolecular interactions. These outcomes through the MD simulation signified for the stable binding of Nos-Hcq conjugation in correlation with molecular docking results.

The principal component analysis (PCA) was performed to insight into the binding clusters variances analyses of the complex system. The binding cluster frames were categorized into two coordinate clusters in black and red color.

PCA calculations of the atomic backbone of the complex system were computed by three conformations PC1, PC2, and PC3 by normal mode molecular dynamics. Obtained results showed that the PCA cluster possesses the highest variability of 30.2% by binding of lead combination to Mpro in terms of internal motions of trajectory. After that, PC2 statistics depicted the minimal variability with 14.43%, and consequently, PC3 calculations 9.64% and suggested the stabilized binding with Eigen scores and minimal atomic motions variability of Nos-Hcq in complexation with Mpro (Figure 7) . Both the conformations in two groups in different PCA analyses were superimposed and aligned, as shown in Figure 7d . The PCA trajectory analyses showed a conformational change in clusters from the first blue cluster to white cluster progressively to red color cluster, similarly recovered from first cluster (black color) to second cluster (red color) through Simple Clustering in PC Subspace (Figure 8) . Furthermore, the dynamical residue cross-correlation map (DCCM) analysis was performed, which depicted the pairwise correlation of atomic coordinates of Nos-Hcq to the Mpro through the magnitude of pairwise cross-correlation coefficient. The map showed the correlation residues (>0.8) of the complex system in blue color, and anti-correlated residues (<-0.4) depicted in red color. The highly pairwise correlated residues (light blue color) of complex trajectory in the map showed the stable binding of ligand file Nos-Hcq. The outcomes were also found in correlation with optimal secondary structures of complex protein dynamics during the MD simulation run. The binding ability of Nos-Hcq towards the Mpro and its essence for coronavirus for replication and regulation suggests the potency of Nos-Hcq and possible role in COVID-19 clinical drug development.

With the computational rationalization of lead conjugate Nos-Hcq against the Mpro, we propose the plausible scheme of chemical synthesis. The drug conjugates of Noscapine 1 and selected Re-purposed small molecules drugs depict the feasible tentative synthetic methodology (Scheme 1). Earlier in 2006, our research group developed a series of microtubule-interfering halogenated noscapine analogs that perturb mitosis in cancer cells followed by cell death. Iodo-Nos 2 was one of the crucial analogs of halogenated Noscapine. 9-Iodo-Nos was synthesized by the iodination of noscapine 1 using pyridine-iodine chloride in acetonitrile solvent with a 75% yield (Aneja et al., 2006) . The target compound 5 can be achieved by the application of a well-known Stille metal cross-coupling reaction. Compound 3, organostannanes, which can act as a precursor for the Stille cross-coupling reaction, can be achieved via radical pathways in the presence of SnMe 3 , DIPEA (base) and MeCN (solvent) (Chen et al., 2016) . The presence of the chloro group at Hydroxychloroquine 4 facile stille coupling reaction to take place at ease. The Pd(PPh 3 ) 4 catalytic system with copper(I) iodide and cesium fluoride in DMF is most effective for coupling iodides. Copper(I) iodide and cesium fluoride combination facile the synthesis of sterically hindered systems such as 5 (Figure 9 ).

We report the efficient combinatorial therapy by conjugating the noscapine (antitussive drug) with potential hydroxychloroquine (Nos-Hcq) against the SARS-CoV-2, through the computational assays with insights into the experimental results. The conjugate Nos-Hcq showed the most definite binding affinity towards the Mpro, among the various noscapine based conjugates with (Chloroquine, Umifenovir, Favlplravir, and Galidesivir). Noscapine (natural source opium poppy) is deemed to form the conjugation with antiviral drugs in trials against SARS-CoV-2. The molecular docking assays depicted that by conjugation binding affinity score almost doubled with Nos-Hcq than with noscapine in the singlet state. The binding mechanism analyses of all conjugations illustrated that the lead conjugation (Nos-Hcq) has strong binding to three major domains of Mpro. Moreover, MD simulation analyses of the trajectory (Mpro-Nos-Hcq) for 100 ns in explicit solvent depicted the stable binding of Nos-Hcq conjugate to Mpro with minimal conformational variations in RMSD and Rg plots. We reveal the drug binding with three major domains of Mpro has high reaction coordinate completion efficacy (energy work profiles) with a significant contribution of immense binding energies and meticulous atomic level stabilization. The DCCP and PCA analyses confirmed the pivotal interactions and stabilize the binding of lead conjugation with Mpro. Besides rationalizing the combinatorial drug therapy, our study provides the fundamental plausible route of chemical synthesis of Nos-Hcq conjugate. This work paves the way for boosting the current medications into the effective combinational strategy to bind target Mpro strongly, which can enhance the inhibitory action of drugs and contributes to fighting against the pandemic COVID-19 with saving time, cost with minimal time screening in the current emergency. Further, the warranted lead drug conjugation needs to be validated experimentally in laboratories and clinical trials.

The authors declare no competing financial interest. 

COVID-19 pandemic: An overview of epidemiology, parthenogenesis, diagnostics and potential vaccines and therapeutics

Synthesis of microtubule-interfering halogenated noscapine analogs that perturb mitosis in cancer cells followed by cell death

Remdesivir for the treatment of Covid-19-preliminary report

4-Bromo-4'-chloro pyrazoline analog of curcumin augmented anticancer activity against human cervical cancer, HeLa cells: In silicoguided analysis, synthesis, and in vitro cytotoxicity

Synthesis of aryl trimethylstannanes from aryl halides: An efficient photochemical method

Emerging coronaviruses: Genome structure, replication, and pathogenesis

In silico oncology drug repositioning and polypharmacology

Individualized networkbased drug repositioning infrastructure for precision oncology in the panomics era

Drug repurposing: New treatments for zika virus infection?

Chloroquine and hydroxychloroquine as available weapons to fight COVID-19

Favipiravir, an antiviral for COVID-19

Noscapine, a possible drug candidate for attenuation of cytokine release associated with SARS-CoV-2. Drug Development Research

Interaction of noscapine with the bradykinin mediation of the cough response

Human coronavirus: Host-pathogen interaction

Breakthrough: Chloroquine phosphate has shown apparent efficacy in treatment of COVID-19 associated pneumonia in clinical studies

A decade after SARS: Strategies for controlling emerging coronaviruses

Bio3d: An R package for the comparative analysis of protein structures

SwissDock, a proteinsmall molecule docking web service based on EADock DSS

Structure of Mpro from COVID19 virus and discovery of its inhibitors

Natural Bis-benzylisoquinoline alkaloids-tetrandrine, fangchinoline, and cepharanthine, inhibit human coronavirus OC43 infection of MRC-5 human lung cells

Understanding the binding affinity of Noscapine with protease of COVID-19 using MD simulation at different temperature

Multiepitope subunit vaccine to evoke immune response against acute encephalitis

Antimicrobial peptide designing and optimization employing large-scale flexibility analysis of protein-peptide fragments

Preclinical evaluation and molecular docking of 1, 3-benzodioxole propargyl ether derivatives as novel inhibitor for combating the histone deacetylase enzyme in cancer

LigPlot þ : Multiple ligand-protein interaction diagrams for drug discovery

PROCHECK: A program to check the stereochemical quality of protein structures

Therapeutic options for the 2019 novel coronavirus (2019-nCoV)

Research and development on therapeutic agents and vaccines for COVID-19 and related human coronavirus diseases

HexServer: An FFT-based protein docking server powered by graphics processors

Outcomes of hydroxychloroquine usage in United States veterans hospitalized with Covid-19

Antiviral agents from plants and herbs: A systematic review

MDTraj: A modern open library for the analysis of molecular dynamics trajectories

Designing of improved drugs for COVID-19: Crystal structure of SARS-CoV-2 main protease M pro

Stereochemical quality of protein structure coordinates

Hydroxychloroquine and ivermectin: A synergistic combination for COVID-19 chemoprophylaxis and treatment?

Coronavirus infections-More than just the common cold

Identification of a novel coronavirus causing severe pneumonia in human: A descriptive study

PLIP: Fully automated protein-ligand interaction profiler

Designing of a novel indoline scaffold based antibacterial compound and pharmacological evaluation using chemoinformatics approach

A novel peptide thrombopoietin mimetic designing and optimization using computational approach

Mechanistic interaction study of bromo-noscapine with bovine serum albumin employing spectroscopic and chemoinformatics approaches

Antibacterial and pharmacological evaluation of fluoroquinolones: A chemoinformatics Approach

WHAT IF: A molecular modeling and drug design program

Structures of two coronavirus main proteases: implications for substrate binding and antiviral drug design

LARMD: Integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor

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

Network-based drug repurposing for novel coronavirus 2019-nCoV/ SARS-CoV-2

A novel coronavirus from patients with pneumonia in China

£-NK and AA contributed equally. N.K, A.A., D.S., R.C., and N.S. designed the studies. N.K., D.S., and A.A, A.K., carried out the in silico experiments. N.K, A.A., A.G., and D.S. wrote the manuscript.

No potential conflict of interest was reported by the authors.