key: cord-0802705-0ajxde2x
authors: Garg, Saksham; Roy, Arpita
title: In silico analysis of selected alkaloids against main protease (M(pro)) of COVID-19
date: 2020-11-09
journal: Chem Biol Interact
DOI: 10.1016/j.cbi.2020.109309
sha: 8e7bf8fd03d0eb46a4a174393a91a2110ab15a69
doc_id: 802705
cord_uid: 0ajxde2x

In the present situation, COVID-19 has become the global health concern due to its high contagious nature. It initially appeared in December 2019 in Wuhan, China and now affected more than 190 countries. As of now preventive measures are the sole solution to stop this disease for further transmission from person to person transmissions as there is no effective treatment or vaccine available to date. Research and development of new molecule is a laborious process; therefore, drug repurposing can be an alternative solution that involves the identification of potential compounds from the already available data. Alkaloids are potential source of therapeutic agents which might be able to treat novel COVID-19. Therefore, In the present study, twenty potential alkaloid molecules that possess antiviral activity against different viral diseases have taken into consideration and scrutinized using Lipinski's rule. Then out of twenty compounds seventeen were further selected for docking study. Docking study was performed using Autodock software and the best four molecule which provides maximum negative binding energy was selected for further analysis. Two alkaloids namely thalimonine and sophaline D showed potential activity to inhibit the M(pro) but to confirm the claim further in-vitro studies are required.

In late December of 2019, few patients were admitted to the hospitals of Wuhan city of China. The patients claimed to have acute respiratory distress and initially were diagnosed to have pneumonia like symptoms (Houlston et al. 2010 ) (Bogoch et al. 2020) . The tests revealed presence of an unknown infectious agent which was later identified to be a member of coronavirus family (Zhao et al. 2020) . Coronaviruses (CoVs) are reported to cause major disorders in respiratory and digestive tracts. Infection is seen in both animals and humans. Some reports suggest that only certain species are infected by CoVs which include mammals, avian species and reptiles (Malik et al. 2020) . Human civilization had also seen earlier the outbreaks of CoVs such as severe acute respiratory syndrome (SARS)-CoV and M pro contains eleven cleavage sites on the large polyprotein and the recognition sequence at most sites is Leu-Gln (cleavage site) (Ser, Ala, Gly) (Mirza and Froeyen 2020) . Inhibition of this enzyme's activity would block the replication of virus because there are no known human proteases which show similar cleavage specificity (Zhang et al. 2020) . A lot of investigations are going on to tackle and curb the virus infection. The secondary metabolites of plants are considered to be highly potential and able to prevent various diseases such as cancer, diabetics, viral infections, etc. Alkaloids are secondary metabolites which possess various pharmacological activities (Roy 2017) . Over the past decade, due to shifted interest towards the natural compounds, various alkaloids are being tested for different purposes including the potential treatment against viral infections. Therefore, in this study, twenty alkaloid molecules (ligand) were selected based on their antiviral potential reported in the literature. Lycrorine is found in the members of Amaryllidaceae family and shows various medical importance (Lamoral-Theys et al. 2009 ). Various studies reported its potential antiviral activity against Human Immunodeficiency Virus (HIV) (Szlávik et al. 2004 ), Herpes Simplex Virus type 1 (HSV1) (Ieven M, van den Berghe DA 1983) , etc. Hemanthamine is another phytochemical present in Lycoris radiate and showed potential antiviral activity against influenza virus (H5N1) (He et al. 2013) . Berberine, an isoquinoline alkaloid found in a significant number of plants such as Berberis vulgaris (barberry), Berberis aristata (tree turmeric), Mahonia aquifolium (Oregon grape) and Hydrastis canadensis (goldenseal). It possesses antiviral activity against influenza virus (H1N1) (Yan et al. 2018) . Thalimonine another alkaloid presents in Thalictrum simplex showed potential to inhibit influenza virus (Serkedjieva and Velcheva 2003) . Hippeastrine is derived from the member of Amaryllidaceae family and found in Lycoris radiate. It has been found to inhibit influenza virus (H5N1) (He et al. 2013) . Hirsutine is an indole alkaloid which is extracted from Uncaria rhynchophylla. It is a potential inhibitor of influenza virus (H3N2) (Takayama et al. 1997) . 5alpha-Hydroxysophocarpine is isolated from Sophora flavescens var. angustifolia and possess anti-viral activity against hepatitis B virus (Ding et al. 2006) . Fangchinoline is found in S. tetrandrae and possess diverse biological activities. Its antiviral activity reported against HIV type 1 (Wan et al. 2012) . Tetrandrine is derived from Stephania tetrandra and act as a calcium channel blocker. Its antiviral activity reported against herpes simplex virus type-1, human coronavirus OC43 (Hu et al. 1997; Kim et al. 2019) . Cepharanthine (CEP) is extracted from Stephania cepharantha and possess potential antiviral activity against human coronavirus OC43 (Kim et al. 2019) . Skimmianine and methoxydihydronitidine both can be isolated from roots of Zanthoxylum nitidum and are being tested for potential activity against HBV (Yang and Chen 2008) . Sophaline D again has shown to inhibit the pathogenicity of HBV virus and can be isolated from Sophora alopecuriodes (Zhang et al. 2017) . A steroidal alkaloid; tomatidine derived from unripe stem and leaves green tomatoes showed reasonable inhibition of Chikungunya virus in various cell lines (Troost et al. 2020) . Emetine can be isolated from Pyschotria ipecacuanha. It has been considered a potential compound against the HIV-1 virus by interfering with the replication machinery of the virus (Valadão et al. 2015) . Dendrobine was collected from Dendrobium nobile and has been investigated for anti-influenza A virus treatment (H1N1, H3N2) . Homonojirimycin was first isolated from Omphalea diandra. Recently homonojirimycin extracted from Commelina communis was studied for antiviral activity against H1N1 virus (ZHANG et al. 2013) . 11-hydroxy Vittatine has been tested for potential activity against H5N1 (Moradi et al. 2018) . Aloperine, a lupine alkaloid exhibited potential activity for anti-influenza A virus treatment. The alkaloid was isolated from Sophora sp. plants belonging to Fabaceae family (Dang et al. 2014) . Atropine was found competent against both DNA HSV-1 and RNA parainfluenza type-3 (PI-3) virus (Özçelik et al. 2011) . The present study aims to find a potential natural compound using computational approach i.e. repurposing the alkaloids against active site of M pro of COVID-19 and which can further utilize for effectively treating COVID-19. Selected alkaloid molecules were subjected to molecular modelling techniques including ADME analysis, molecular docking and molecular dynamic simulation to find out their interactions and stability with M pro of COVID-19.

Protein/Macromolecule 3-dimensional (3D) structure of COVID-19's main protease (M pro ) (pdb id 6LU7) was obtained from PDB website ("Protein Data Bank"). The structure of 6LU7 is formed of two A chains hence, making the molecule a homodimer. Along with the primary chains, a hetatom named N3 is found as an in-activator of the protein.

PubChem repository ("PubChem") was used to obtain the structure of alkaloids required for the analysis in .sdf format. Later, Biovia Discovery Studio Visualizer v19.1.0.18287 was used to convert the structures into .pdb format. 3D structures of following compounds were used in this study: Lycorine ( 

Lipinski's Rule of Five was applied to 10 compounds that were selected for this investigation. Filters like Molecular weight of the ligand (<500Da), high lipophilicity (LogP<5), Number of hydrogen bonds donors (<5), Number of hydrogen bond acceptors (<10) and Molar refractivity (40-130) (Ghose Rule) were used to carry out the further selection of ligand molecules. Violation of more than 2 of the above stated parameters debarred further analysis of particular molecule (Lipinski 2004) . Parameter details were taken from PubChem portal ("PubChem") except molar refractivity. Molar refractivity was calculated with the use of an online tool SwissADME ("SwissADME").

Autodock v4.2.6 was used to carry out the molecular docking of each ligand with 6LU7. Both protein and ligand optimization were carried out by removing water, adding polar hydrogens and charges including Kollman and Gasteiger. Finally, the removal of N3 i.e. hetatom from macromolecule concluded the preparation of ligand and protein. A grid box of 60×60×60 with a spacing of 0.375Å was used for active site residues. Lamarckian GA output was obtained from docking studies. For each ligand, docking procedure was repeated thrice and out of 10 conformations obtained for each docking run, the final best conformation with minimum binding energy was taken into consideration and converted into a 2-dimensional diagram showing interaction of ligand with active site residues using Biovia Discovery Studio Visualizer v19.1.0.18287.

A more comprehensible analysis of physiochemical properties was used to continue the filtration of potent ligand molecule. Bioavailability radars of ligands having better results than control compound were obtained using SwissADME web based tool. A total of 6 parameters were used to scrutinize the compounds: solubility, size, polarity, lipophilicity, flexibility and saturation. Ligands deviating from the standardized values suggested non oral bioavailability and hence were debarred from further testing (Daina et al. 2017) .

Simulation of receptor-ligand docked complex was done using Desmond-Maestro 2020 module. By default, Desmond uses algorithms which are entitled to perform high speed and precision-based MD simulations. For conformational stability of protein-ligand complex, a system setup based on default algorithms was used. The protein-ligand complex was solvated in orthorhombic TIP3P water model. System was neutralized with three counter ions of Na + at 0.15M salt concentration. Optimized potentials for liquid simulations-AA (OPLS-AA) 2005 force field was assigned for the system. Employed SHAKE/RATTLE algorithm limited the movement of all the covalent bonds with default tolerance value of 10 -7 Å. NVT was utilized as an ensemble class to relax the system before the simulation begins. The simulation was set at 300K and 1 bar pressure for 50ns. Thermostat method was set on Nose-Hoover chain with relaxation time of 1ps while Martyna−Tuckerman−Tobias−Klein (MTTK) barostat method was employed with 2ps relaxation time to maintain temperature and pressure. The approaches taken to set up the system was combined by RESPA integrator. A trajectory of 50ns was setup with 1000 frames to show the interaction between protein and ligand based on the literature reports ("User Manual, Impact 6.7" 2015; Kumar et al. 2020 ).

Biological activity prediction Prediction of biological activity was done using the PASS webserver ("PASS webserver"). PASS webserver allows us to predict the potential biological effect of the compound using MNA (multilevel neighbors of atoms) descriptors. This particular method prediction is entirely based on the chemical structure of the compound (Goel et al. 2011; Gulzar et al. 2019) .

The selected 20 compounds were subjected to a preliminary comparison between few properties of ligand molecule (Table 1 ). This comparison helped to rule out a few compounds which were less likely to be a drug candidate. Along with all the properties as mentioned in Lipinski's Rule of 5, one of the Ghose Rule was also taken into consideration as filter for comparison among the different alkaloids. Qualifying majority of the parameters doesn't guarantee a particular compound to be a drug; it only tells about drug-likeliness and helps in eliminating the weak compounds in the preclinical phase only. The result of analysis showed that among the selected 20 alkaloids, tetrandrine, cepharanthine and fangchinoline showed violation of more than two parameters, hence they are eliminated from further studies. The remaining 17 compounds were considered to be as a potential drug candidate and molecular docking studies were performed to predict the interaction of protein-ligand complex.

Molecular docking is a computational method used in the field of drug designing. Docking results predict the residues of protein that are interacting with ligand's atoms and generates a lowest energy conformation complex. All the seventeen selected alkaloids were docked at the active site of M pro (PDB ID: 6LU7) protein of COVID-19. THR24, THR26, PHE140, ASN142, GLY143, CYS145, HIS163, HIS164, GLU166, HIS172 ( Figure 1A ) were found at the active site of the M pro protein. The binding energy was taken as the comparison parameter and N3 (Native inhibitor) was considered as control molecule. Binding energy of N3 with 6LU7 was -8.15 kcal/mol (Table  2 and Figure 2E ). Thalimonine, Emetine, Sophaline D and Tomatidine were showing better binding energy than N3.

Thalimonine-M pro complex showed the least binding energy of -8.39 kcal/mol. SER144 and CYS145 were forming conventional H bonds with CYS145 forming 2 H bonds along with Pi-donor and Alkyl interactions. Electrostatic bond formation was seen in case of HIS41 with both benzene rings present in the compound. Other interactions were also observed with HIS41 residue like Pi-Pi cation ( Figure 2A) and Pi-alkyl. Two sulfur bond formation was taking place with sulfur containing residues i.e. MET165 and MET49. One more alkyl bond formation takes place with MET49. HIS164 and ASP187 were forming a Pi-donor H-bond with the ligand. Sophaline D-M pro docked complex gave -8.79 kcal/mol as the minimum binding energy. Interaction can be classified into 6 categories namely: Conventional H-bond, Alkyl, pi-alkyl, carbon hydrogen and van der Waal bonds. Among the various residues, HIS163 was seen bonding via conventional hydrogen bond. HIS41 formed pi-alkyl bond and CYS145 and MET165 can were bound to ligand by alkyl bond formation. HIS164 residue was found interacting via carbon hydrogen bond. 12 residues stabilizing the complex by van der Waals attractive force can be observed too ( Figure 2B ).

Tomatidine-M pro complex was found to have minimum binding energy of -9.58 kcal/mol. Only 1 conventional hydrogen bond was observed with HIS164 residue. GLN189 bonded using carbon hydrogen bond. LEU141, MET49 and CYS145 interacted with the ligand by forming alkyl bonds while due to presence of ring structure in the residues HIS41, HIS163 and HIS172 formed pi-alkyl type of bond with tomatidine. Along with these major interaction 8 more van der Waal attractive interaction were seen in the bound protein-ligand complex ( Figure 2C ).

Emetine-M pro complex showed minimum binding energy of -10.17 kcal/mol. A total of 6 types of bond formation was observed in the docked complex. HIS163 and GLU166 were forming conventional H-bond with emetine. GLU166 was also interacting using a pi-anion interaction. MET165 alone was forming alkyl as well as pi alkyl bonds. LEU141 and HIS164 formed carbon hydrogen bond with the ligand. In the surrounding 11 other residues can be seen interacting through weak van der Waals interaction thereby, adding to the stability of the overall conformation ( Figure 2D ).

Comparison with control (N3) left us with 4 ligands viz., thalimonine, emetine, tomatidine and sophaline D.

Observing bioavailability radars of four compounds revealed that all the ligands except tomatidine were found orally bioavailable. Tomatidine failed to follow the standard values of lipophilicity and insolubility. Other three ligands can be seen to lying in the shaded pink region of standard values.

Three compounds were shortlisted, thalimonine, sophaline D and emetine. 6LU7 in complex with mentioned three ligands was subjected to a 50ns MD simulation individually. The average potential energy of the thalimonine-M pro , sophaline D-M pro and emetine-M pro complexes was found to be -115,805 kcal/mol, -117,514 kcal/mol and -117,360 kcal/mol respectively. Along with potential energy total energy parameter was also calculated during simulation time, -93,656 kcal/mol, -95,415 kcal/mol and -95,270 kcal/mol respectively were the total energy values obtained.

Upon binding of a small molecule, the stability of the conformation of the protein is a fundamental property to explore and thus, gain insight about changes and dynamics of the protein structure. Root mean square deviation (RMSD), Root mean square fluctuation (RMSF), Radius of Gyration (rGyr) and Solvent accessible surface area (SASA) plots can be utilized to estimate the properties like compactness of the protein.

RMSD plot ( Figure 3A ) of thalimonine-M pro complex can be observed getting stabilized at the end of the simulation suggesting protein stability after binding of the ligand. Though, a lot of steep peaks can be seen during the simulation which are a result of internal vibration of molecule. A maximum deviation of 1Å was seen which is considered acceptable for the protein stability. RMSD graph ( Figure 4A ) of 6LU7 in complex with sophaline D ligand can be observed ending with a lower RMSD indicating stabilized protein. A sudden peak with deviation of nearly 1Å can be seen at about 45ns which might be the result of sudden internal vibration in the protein but henceforth from 45-50ns RMSD plot can be seen going towards lower numbers in the graph. RMSD plot ( Figure  5A ) of M pro with emetine was obtained and during the most time of the simulation the values can be seen fluctuating around a fixed number. It was only after 40ns of the simulation of RMSD plot went a little higher than before. A structural deviation was seen in the trajectory and the deviation is nearly 1Å only hence, the elevated graph can be attributed to internal vibrations. RMSF plot helped in understanding the residual vibrations in three different complexes. Since, alpha helixes and beta strands are rigid hence, less fluctuation should be seen with the residues involved in such secondary structures. RMSF plot (( Figure 3B ) of 6LU7 with thalimonine can be observed to follow the same principle with high peaks at the unstructured part of the protein. The residues of the binding pocket can be seen obtaining lower values of the graph indicating less conformational change at the binding pocket. Similarly, RMSF plot ( Figure 4B ) of M pro with sophaline D can be observed to follow the same lines, residues making alpha helixes and beta sheets with the lower values. Relatively, the RMSF graph obtains less values with less peaks. RMSF graph of 6LU7-emetine ( Figure 5B ) complex can be seen obtaining relatively higher values and peaks. Though the residues at binding site and secondary structure obtained smaller values. Also, in all three graphs the B-factor graph can be running mostly parallel to the RMSF plot indicating towards protein stability.

Radius of gyration (rGyr) is yet another parameter which is linked with the tertiary structure and general conformational state defining our understanding of compactness and folding of the protein. All three complexes were allowed to yield the rGyr plot. A fluctuation of 0.1Å, 0.1Å and 0.4Å around the average values were observed for M pro -thalimonine ( Figure 3C ), M pro -sophaline D ( Figure 4C ) and M pro -emetine ( Figure 5C ) respectively. Though, in case of M pro -emetine complex a sudden drop in rGyr was observed between 15ns to 25ns but immediately after 25ns the graph was restored to its average value. This sudden deviation might be due to protein's packing. All three plots attained the equilibrium around their average value thus, suggesting complex stoutness through the simulation course.

SASA is the surface area accessible to the water molecule and study of SASA can be used to investigate the conformational dynamics. Average value of SASA observed were: 200Å 2 , 120Å 2 and 450Å 2 respectively for 6LU7thalimonine ( Figure 3D ), 6LU7-sophaline D ( Figure 4D ) and 6LU7-emetine ( Figure 5D ) complex. In each case result, the equilibrium is reached around the average values. Emetine showed highest SASA value suggesting larger surface are in contact with water molecules i.e., more number of inner residues coming in contact with the solvent while thalimonine and sophaline D complexes are more stable with less inner residues interacting with surroundings. Though the attained equilibrium in each case suggests that the complexes are structurally stable.

Computing the 6LU7-thalimonine complex, a new residue, THY25 interacted with ligand almost 100% of the simulation time via water bridge formation. Out of all residues, four residues i.e. CYS145, ASN142, HIS41 and TYR54 formed H-bond with ligand. HIS41 showed pi-pi stacking majority of the simulation time and forming water bridges for the rest of the remaining simulation time. Majority of the amino acids including THR24, CYS44, TYR54, ASN119, ASN142, GLY143, SER144, CYS 145, VAL 186, ASP187, ARG188 and GLN189 were involved in water bridge bond formation which were slightly less weak than conventional H-bonds. Other residues like LEU27, MET49, PRO52, MET165 were involved in hydrophobic interaction which include pi-pi stacking, pi-cation and other non-specific interactions ( Figure 3E ). Figure 3F represents a timeline representation of number of contacts made by residues with ligand. The upper panel showed total number of contacts at any instant throughout the trajectory course. On average at any given time 3-4 contacts were made by the 6LU7 residues. The lower panel showed which residue makes contact with ligand molecule at any instantaneous moment of the trajectory. During the middle and later stages of the simulation THR26 was interacting with ligand in a significant amount, often making more than 4 contacts. HIS41 was also stabilizing the ligand over 50ns time. At the later stages of the trajectory CYS44 also plays role in stabilizing the overall interaction with TYR54, CYS145 and ASP187.

Considering the 6LU7-sophaline D complex, in total of 5 H-bonds can be observed. THR190 and GLN192 both residues interacted with the residue for 100% of the simulation time. Along with these residues, GLY143, CYS145 and SER144 were also seen to form H-bond with the ligand for a brief period of time. Few residues were also seen to change the nature of bond formation. GLY143 for majority of interaction time formed H-bond while for a brief time period also interacted through water bridge formation while SER144 interacting in a vice-versa manner. CYS145 was seen to deviate along all three types of interaction, for almost equal amount of fraction, H-bond and hydrophobic interactions can be observed while water bridge relatively for a lower fraction. On the other and MET145, MET49, PRO168 interacted with hydrophobic interactions. Other residues like HIS41, LEU141, ASN142, HIS163, GLU166, GLN189 all contributed towards the stability of the complex via water bridge bond formation ( Figure 4E ). Figure 4F shows the timeline representation of the residues interacting with the ligand at any specific time. An average of 4-6 bonds were observed with the ligand at any given point of time during the 50ns simulation. GLN192 and THR190 can be seen to hold the ligand in the binding pocket during the full 50ns time. While other residues like HIS41, CYS145, GLY143 and MET165 binding with the ligand in a very scattered manner. Throughout the trajectory at no given period of time the ligand is devoid of any interaction.

In the case of 6LU7-emetine complex, the highest interaction fraction can be observed to be only 40%. Many residues formed 2 types of bonds. GLN189 and THR190 formed H-bond and water bridge interaction for almost equal time. ASN142, GLU166 and SER46, all interacted via water bridge formation for majority of their respective times and very briefly by forming H-bond with the ligand. PRO168 and ALA199 both formed hydrophobic interactions along with water bridges. 12 other residues individually formed water bridges and hydrophobic interactions for very short fraction of time contributing to the complex stability ( Figure 5E ). The timeline of the residues in contact with ligand can be seen in Figure 5F . Observations from the upper panel suggested 2 contacts to an average number of contacts. We can also see significant number of intervals were no contact with the ligand was made. The lower panel confirms the same. THR190 and GLN189 were seen interacting most with emetine though, contacts are scattered. Other residues also make contact with ligand but overall, the contacts are distributed unevenly with intermediate complete absence of any interactions.

Secondary structural dynamics can be used to understand the conformational changes in the protein and folding behavior of the protein. It is a result of changes in the residue upon binding of different ligands. Table 4 represents the percentage composition of secondary structure of the whole protein. Analysis showed that with thalimonine as the ligand, 41.87% of all the residues were involved in secondary structure formation. This number was observed slightly decreased with sophaline D as the ligand (41.56%) and with the emetine as the ligand molecule further drop was observed (40.19%). Both alpha-helixes and beta-strands were reduced.

Biological activity prediction PASS server was being used to make the effective biological activity predictions. The analysis resulted in similar biological activity for the selected three compounds. Thalimonine, sophaline D and emetine have shown antitussive, antiviral (Rhinovirus) and antitussive predictions respectively. All the biological activities predicted were related to respiratory ailments and Pa ranged from 0.325 to 0.827 when Pi < Pa. Table 3 shows the predicted activities of the three compounds.

M pro of COVID19 is predicted to contain 306 amino acid residues and structural analysis showed that it has total of 3 domains. The binding site is located in-between domain I and domain II. Both domains contribute a single residue for the formation of a catalytic site which is essential for substrate binding. Domain I contain His41 residue and domain II contains Cys145 residue which are conserved in SARS-CoV family (Mirza and Froeyen 2020) . From literature survey it was found that His41 acts as a base and -SH group of Cys145 mediates the electrostatic trigger to start the chemical reaction and in collaboration it is known as Cys-His dyad. Functional M pro is an essential requirement for maturation of the virus and the catalytic dyad is essential to carry out substrate conversion by M pro . Even a slight mutation in either of the residues leads to enzyme inactivation (Chang 2010) . Hence, considering the genomic and structural similarities the dyad becomes a potential target for M pro of Covid 19. Comparing all the filters, parameters and mechanisms, thalimonine and sophaline D can be considered for the further analysis to inhibit the M pro of COVID-19. Both were predicted to have respiratory track related biological activity and can be seen interacting with the catalytic dyad of the 6LU7. Emetine showed the maximum negative binding energy of all docked complexes but the lack of interaction with the catalytic dyad residues and frequent absence of any contact with ligand during simulation makes emetine a hard choice.

For combating the pandemic caused by the COVID19 many studies are aiming to find antiviral therapies. Most of the studies as of now are concentrated towards finding a cure in chemically synthesized or earlier used drugs (drug repurposing) for other diseases. In this study, the aim was to find a natural solution which can effectively solve the problem. Considering that as a fact, a five-step filter study was conducted consisting of ADME analysis, followed by molecular docking, bioavailability radar analysis, molecular dynamic simulation and prediction of biological activity on different alkaloids which are known to show some antiviral activity. In this study, all the twenty alkaloids which were chosen filtered through ADME parameters. Further, the compounds with binding energy lower than N3, 

After the MD simulation biological activity was also predicted. Overall, the study concludes two alkaloids namely thalimonine and sophaline D having potential activity to inhibit the M pro but to confirm the claim further in-vitro studies

Pneumonia of unknown aetiology in Wuhan, China: potential for international spread via commercial air travel

Quaternary Structure of the SARS Coronavirus Main Protease

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

Identification and synthesis of quinolizidines with anti-influenza A virus activity

Quinolizidine Alkaloids with Anti-HBV Activity from Sophora tonkinensis

PASS-assisted exploration of new therapeutic potential of natural products

Elucidation of interaction mechanism of ellagic acid to the integrin linked kinase

Amaryllidaceae alkaloids inhibit nuclearto-cytoplasmic export of ribonucleoprotein (RNP) complex of highly pathogenic avian influenza virus H5N1

Meta-analysis of three genome-wide association studies identifies susceptibility loci for colorectal cancer at 1q41, 3q26.2, 12q13.13 and 20q13

Tetrandrine potently inhibits herpes simplex virus type-1-induced keratitis in BALB/c mice

Plant antiviral agents. IV. Influence of lycorine on growth pattern of three animal viruses

Natural Bis-Benzylisoquinoline Alkaloids-Tetrandrine, Fangchinoline, and Cepharanthine, Inhibit Human Coronavirus OC43 Infection of MRC-5 Human Lung Cells

Withanone and caffeic acid phenethyl ester are predicted to interact with main protease (M pro ) of SARS-CoV-2 and inhibit its activity

Lycorine, the Main Phenanthridine Amaryllidaceae Alkaloid, Exhibits Significant Antitumor Activity in Cancer Cells That Display Resistance to Proapoptotic Stimuli: An Investigation of Structure−Activity Relationship and Mechanistic Insight

Anti-influenza A Virus Activity of Dendrobine and Its Mechanism of Action

Lead-and drug-like compounds: the rule-of-five revolution

Emerging novel coronavirus (2019-nCoV)-current scenario, evolutionary perspective based on genome analysis and recent developments

Structural elucidation of SARS-CoV-2 vital proteins: Computational methods reveal potential drug candidates against main protease, Nsp12 polymerase and Nsp13 helicase

Alkaloids as the natural anti-influenza virus agents: a systematic review

Cytotoxicity, antiviral and antimicrobial activities of alkaloids, flavonoids, and phenolic acids

PASS webserver

A Review on the Alkaloids an Important Therapeutic Compound from Plants

Vitro Anti-Influenza Virus Activity of the Pavine Alkaloid

Alkaloids from Leucojum vernum and Antiretroviral Activity of Amaryllidaceae Alkaloids

Discovery of anti-influenza A virus activity of a corynanthe-type indole alkaloid, hirsutine, in vitro and the structure-activity relationship of natural and synthetic analogs

Tomatidine, a natural steroidal alkaloid shows antiviral activity towards chikungunya virus in vitro

Natural plant alkaloid (Emetine) inhibits HIV-1 replication by interfering with reverse transcriptase activity

Fangchinoline Inhibits Human Immunodeficiency Virus Type 1 Replication by Interfering with gp160 Proteolytic Processing

Antiinfluenza activity of berberine improves prognosis by reducing viral replication in mice

Alkaloids from the roots of Zanthoxylum nitidum and their antiviral and antifungal effects

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

Four Matrine-Based Alkaloids with Antiviral Activities against HBV from the Seeds of Sophora alopecuroides

Homonojirimycin, an alkaloid from dayflower inhibits the growth of influenza A virus in vitro

Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak

Figure 2. (A) Interaction of Thalimonine with M pro (B) Interaction of Sophaline D with M pro (C) Interaction of Tomatidine with M pro (D) Interaction of Emetine with M pro (E) Interaction of N3 (control) with M pro

√☐ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.√☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:Saksham Garg Arpita Roy