key: cord-1029229-2e4bj7zf authors: Ullah, Anwar; Ullah, Kifayat title: Inhibition of SARS-CoV-2 3CL M(pro) by Natural and Synthetic Inhibitors: Potential Implication for Vaccine Production Against COVID-19 date: 2021-04-12 journal: Front Mol Biosci DOI: 10.3389/fmolb.2021.640819 sha: 1e03e0f91d3f66195996106d66824e3261df751e doc_id: 1029229 cord_uid: 2e4bj7zf COVID-19 has created a pandemic situation all over the world. It has spread in nearly every continent. Researchers all over the world are trying to produce an effective vaccine against this virus, however; no specific treatment for COVID-19 has been discovered -so far. The current work describes the inhibition study of the SARS-CoV-2 main proteinase or 3CL M(pro) by natural and synthetic inhibitors, which include 2S albumin and flocculating protein from Moringa oleifera (M. oleifera) and Suramin. Molecular Docking study was carried out using the programs like AutoDock 4.0, HADDOCK2.4, patchdock, pardock, and firedock. The global binding energy of Suramin, 2S albumin, and flocculating proteins were −41.96, −9.12, and −14.78 kJ/mol, respectively. The docking analysis indicates that all three inhibitors bind at the junction of domains II and III. The catalytic function of 3CL M(pro) is dependent on its dimeric form, and the flexibility of domain III is considered important for this dimerization. Our study showed that all three inhibitors reduce this flexibility and restrict their motion. The decrease in flexibility of domain III was further confirmed by analysis coming from Molecular dynamic simulation. The analysis results indicate that the temperature B-factor of the enzyme decreases tremendously when the inhibitors bind to it. This study will further explore the possibility of producing an effective treatment against COVID-19. A new virus named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified in patients in China in December 2019 (Kotta et al., 2020) . It spread throughout the country and world quickly and infected millions of people all over the world (Kneller et al., 2020) . Till now (November 2020), 55.6 million people have been detected with this virus of which 35.86 million have been recovered and 1.34 million have died (Johns Hopkins University). The disease produced by SARS-CoV-2 is termed COVID-19 (Hussin et al., 2020; , which is a short name given to this disease by the World Health Organization (WHO, 2020) . The coronavirus spread through the air and physical contact among people (Graham Carlos et al., 2020; Helmy et al., 2020; Zhao et al., 2020) . The usual symptoms of COVID-19 include mild fever, cough, lethargy, dyspnea (difficulties in breathing), and anosmia (loss of smell) and taste (ageusia) (Kotta et al., 2020) . These symptoms usually appear after 5 days of infection by the virus . Interestingly some of these symptoms like mild fever, cough, lethargy, and dyspnea are common among both betacoronavirus and COVID-19 ; however, COVID-19 displays some distinctive clinical symptoms like sore throat, a runny nose (rhinorrhea), and sneezing (sternutation) (Lee et al., 2003; Assiri et al., 2013) . One way to prevent the COVID-19 from spreading among people is to keep a suitable distance of 1.5-2 meter as recommended by WHO Kotta et al., 2020) , although a recent study has suggested that the virus can travel more than 2 m in the air (Setti et al., 2020; van Doremalen et al., 2020) . The lockdown option is used in all countries of the world to achieve this social distancing and it has worked tremendously like in China Wu et al., 2020) . Currently, there is no specific treatment for COVID-19 (Barati et al., 2020) , though some antiviral drugs like redeliver, oseltamivir, lopinavir/ritonavir, ganciclovir chloroquine, and hydroxychloroquine are used that can produce some relief to the patients (Costanzo et al., 2020; Kumar et al., 2020; Ledford, 2020) . The genomics and proteomics of SARS-COV2 have been described in the literature (Vandelli et al., 2020) . The structure of this new virus is composed of single-stranded ribonucleic acid (RNA) and displays high sequence identity to other betacoronavirus such as SARS-CoV and MERS-CoV (Middle East respiratory syndrome coronavirus) (Cascella et al., 2020) . These viruses use a specific protein named spike (S) protein to adhere specifically to the angiotensin-converting enzyme (ACE2) on the host cell (Park et al., 2019; Turoòová et al., 2020) . Besides spike glycol protein, the SARS-COV2 contains proteins like 3CL M pro [also called the main proteinase (M pro )] and RNA-dependent RNA polymerase (RdRp) (Jeong et al., 2020) . The life cycle of SARS-COV2 begins when the virus infects the host cell through the interaction of S protein with the angiotensin I-converting enzyme-2 (ACE2) (V'kovski et al., 2020) . The S protein has two subunits called S1 and S2 , S1 it uses to attach to the N-terminal of ACE2, and the S2 subunit assists in the binding of the protein to the host membrane. This results in the binding of the virus to the membrane of the host cell. Consequently, the disruption of the membrane of the host cell occurs and endocytosis takes place (V'kovski et al., 2020) . The furin proteinase and transmembrane serine proteinase 2 of the host cells cause the cleavage of S protein at the S1/S2 boundary position (V'kovski et al., 2020) , which allow the entry of transmembrane serine proteinase 2-dependent entrance to the host cells (Belouzard et al., 2009; Hoffmann et al., 2020; Walls et al., 2020) . The polycistronic RNA of the virus is released into the cytoplasm. The ribosomal-1 frameshifts then translates the replicase gene either into replicase polyprotein pp1a or into pp1ab (∼750 kDa, nsp1-16). This process occurs near the 3 -end of ORF 1a. This autoproteolytic cleavage results into 16 non-structural proteins (NS) by two ORF1a encoded proteinase domains (Brierley et al., 1989; Herold et al., 1993; Thiel et al., 2001 Thiel et al., , 2003 Harcourt et al., 2004; Prentice et al., 2004; Ziebuhr, 2004) . The two other proteinases assisting in these proteolytic cleavages include the main proteinase M pro (3CL M pro ) and papain-like proteinase (PL pro ) (Hegyi and Ziebuhr, 2002) . The polyprotein pp1ab is cleaved by M pro (Ziebuhr et al., 2000; Hegyi and Ziebuhr, 2002) . The replication (production of the entire genome) or transcription (synthesis of intermittent mRNAs) is intervened by cytoplasmic enzyme complex termed replicase-transcriptase complex (Gorbalenya et al., 2006; Pasternak et al., 2006; Sawicki et al., 2007) . The key proteins (structural and accessory) are translated from these transcripts; consequently the viruses are released into the cell (V'kovski et al., 2020) . The two important proteins in the life cycle of SARS-CoV-2, are the S protein and 3CL M pro (Kneller et al., 2020; V'kovski et al., 2020) . As discussed earlier, the S protein help the virus to binds to the host cell and to facilitate its entry to the host cell (Duan et al., 2020) , while 3CL M pro or the main proteinase assists in the processing of the polyproteins (Kneller et al., 2020) . Owing to the main roles of these two proteins, researchers from all over the world are targeting these proteins to find out a new treatment for COVID-19 (Kotta et al., 2020) . Taking this into consideration, the current work has been designed to test the efficacy of natural and synthetic inhibitors (2S albumin and flocculating proteins of Moringa oleifera and Suramin), against 3CL M pro and discover a new treatment for this pandemic disease. The atomic coordinates of SARS-CoV-2 3CL M pro , Suramin, and 2S albumin were retrieved from the Protein Data Bank (PDB), with PDB IDs: 6WQF (Kneller et al., 2020) , 6CE2 (SVR) (Salvador et al., 2018) , and 5DOM (Ullah et al., 2015) . The structure of flocculating protein was obtained as a model using the Swiss Model (Waterhouse et al., 2018) . The three-dimensional atomic structure of 2S albumin from M. oleifera was used as a template (74% sequence identity). The ligands and crystallographic water molecules were removed from the protein and the H-atoms were added. The ionization states of the atoms were kept in the ligand as mentioned in the database. The optimization of the ligand geometry was done using the AM1 method (Dewar et al., 1985) . The partial charges of the ligands were calculated by AM1-BCC method (Jakalian et al., 2002) . The atoms type, bond angle, dihedral, and van der Waals parameters for the ligands were assigned using the general AMBER force field (GAFF) method (Wang et al., 2004) . The programs used for molecular docking include AutoDock 4.0 (Morris et al., 2009) , pardock (Gupta et al., 2007) , patchdock (Schneidman-Duhovny et al., 2005) , HADDOCK2.4 (van Zundert et al., 2016) , and the refinement of the docked ligands with protein was carried out using firedock (Mashiach et al., 2008) . The binding affinity of the docked ligands were find out using Kdeep web server (Jiménez et al., 2018) . The interactions (hydrogen bonds and hydrophobic contacts) between 3CL M pro was determined using LigPlot (Wallace et al., 1995) from PDBsum web server (Laskowski et al., 2018) . The MDMoby and MDweb programs (Hospital et al., 2012) , GROMACS (Berendsen et al., 1995) , AMBER16 (Case et al., 2005; Maier et al., 2015) were used for Molecular Dynamic Simulation as described previously (Ullah et al., 2019; Ullah and Masood, 2020) . The all-atom-protein interaction was found out using FF14SB force field (Darden et al., 1993) . The online server H + + (Anandakrishnan et al., 2012) was used for the determination of the protonation states of the amino acid side chain at pH 7.0. The neutralization of the system was carried out using Cl-. The minimization of the simulation system was carried out in order to remove the clashes in the atomic position, structural errors (bond length and bond angle). This minimization was done by executing a 500-step descent (SD) minimization, accompanied by a 2 ns position restricted MD simulation with NVT and NPT ensemble separately (Zhang et al., 2013) . Subsequently, it was put in a rectangular box of TIP3P water, and extended to a minimum of 20 Å from any protein atom. The system was heated gradually from 0 to 350 K for 250 ps with a constant atom number and volume. The protein was kept with a constant force of 10 kcal/mol.Å 2 . A constant atom number, pressure, and temperature (NPT) ensemble was conducted for 500 ps to attain the equilibration step. The simulation was executed for 100 ns with a 4-fs time step. The pressure was kept at 1 atm using Nose ì-Hoover Langevin Piston algorithm (Tu et al., 1995) and the temperature was kept at 300 K, using Langevin coupling (Washio et al., 2018) . The long-range electrostatic interactions were calculated using the particle-mesh Ewald (PME) method (Darden et al., 1993) , by retaining the cutoff distance of Van der Waals interactions at 10 Å. The protein and ligands were prepared for surface charge distribution using PDB2PQR (Dolinsky et al., 2007) and the charges were visualized using ABS Tools from PyMOL (DeLano, 2000) . The Overall Structure of SARS-CoV-2 3CL M pro The three-dimensional structure of SARS-CoV-2 3CL M pro has been described by Kneller et al. (2020) with PDB ID: 6M03. The structure is composed of 306 amino acid residues and these amino acid residues fold into distinct three domains, named domains I, II, and III ( Figure 1A) . Domain I is composed of amino acid residues, from, Phy8-Tyr101, and has four α-helices and seven beta-strands. Domain II (amino acid residues, Lys102-Pro184) comprises seven beta strands only, whereas domain III (amino acid residues Thr201-Val303) contains five alpha-helices only. The enzyme active site is situated at the junction of domains I and II and comprises the amino acid residues His41 and Cys145 (Figure 1B) , which make a dyad (Cys145-His41) instead of the triad (His47-Asp102-Ser195) as in the case of classical serine proteinases (Ullah et al., 2018) . A catalytic water molecule is also bound to His41 and helps in the catalytic process of this enzyme ( Figure 1B) . The enzyme is active in the dimeric state and the flexibility of domain III is required for its dimerization (Kneller et al., 2020) . The binding energy calculated for interaction between SARS-CoV-2 3CL M pro and Suramin was ∼ −42 kcal/mol ( Table 1) . All the other form of bond energies are listed in Tables 1, 2. Global binding energy, attractive VdW (van der Walls forces), Repulsive VdW, ACE (atomic contact energy) and HB (contribution of the hydrogen bonds to global energy) for the interaction between SARS-CoV-2 3CL M pro , Suramin, 2S albumin and flocculating protein). Suramin binding site is between the two domains (Domains II and III) of SARS-CoV-2 3CL M pro (Figures 2A-C) . The amino acid residues of SARS-CoV-2 3CL M pro that interact with Suramin include Lys102, Pro108, Gln110, Asp155, Glu240, and His246 (Figures 1D,E) . The Kdeep results indicate that both equilibrium dissociation constant (pKd) and Gibb's free energy ( G) are large (Table 3) , which further confirmed the binding between SARS-CoV-2 3CL M pro and the three ligands (Suramin, Figure 3 and Supplementary Table 1) . Suramin is a drug that is used to treat African sleeping sickness and river blindness (Lima et al., 2009) . Suramin has been shown to inhibit Human α-thrombin (Lima et al., 2009) , snake venom phospholipases A2 (Salvador et al., 2018) , snake venom serine proteinases (Ullah et al., 2018) , severe Fever with thrombocytopenia syndrome virus nucleocapsid protein (Jiao et al., 2013) , murine Norovirus RNA-dependent RNA polymerase (Mastrangelo et al., 2012) , and Leishmania mexicana pyruvate kinase (Morgan et al., 2011) . In most of these cases, the Suramin binds toward the C-terminal of the proteins and restrict the motion of the C-terminal (Lima et al., 2009; Ullah et al., 2018 ). In the current study, Suramin binds toward the N-terminal of SARS-CoV-2 3CL M pro (Figures 2A-E) . The binding energies for SARS-CoV-2 3CL M pro , 2S albumin and flocculating protein were ∼ −9.12 and ∼ −15 kJ/mol, Frontiers in Molecular Biosciences | www.frontiersin.org respectively ( Table 1) . The other form of bond energies come from docking as indicated in Tables 1, 2. The amino acid residues involved in these interactions, include S139, T139, G302, Q299 (SARS-CoV-2 3CL M pro ), R143, Q97 (2S albumin) and Q15, and Q38 (Flocculating protein). The interactions between SARS-CoV-2 3CL M pro , 2S albumin, and flocculating protein are Frontiers in Molecular Biosciences | www.frontiersin.org largely electrostatic (Figures 3A-E and 4A-E) . In both cases, the ligands binding site is between the two domains (Domains II and III) of SARS-CoV-2 3CL M pro . The LigPlot analysis shows a total of three hydrogen bonds between SARS-CoV-2 3CL M pro and both 2S albumin and flocculation protein, while the number of hydrophobic interactions were 130 and 152 for 2S albumin and flocculating protein, respectively (Supplementary Figures 4, 5) . The MD simulation analysis indicates that the flexibility of SARS-CoV-2 3CL M pro decreases tremendously when the ligands bind to it (Supplementary Figure 1) . For Suramin as an inhibitor, the fluctuation increases a little bit (temperature B-factor increases from 14 to 16) ( Supplementary Figure 1 and Figures 2A-D) , while in the case of 2S albumin and flocculating proteins the fluctuation decreases (temperature B-factor decreases from 12 to 10, respectively (Supplementary Figure 1, Figures 3A-D and 4A-D) . The RMSD vs. time graph indicates that the interaction between SARS-CoV-2 3CL M pro and the three ligands was stable throughout the simulation process (Supplementary Figure 6) . Suramin can make 1-5 hydrogen bonds, while both 2S albumin and flocculating protein can make 2-5 hydrogen bonds according to 100 ns MD simulation analysis (Supplementary Figure 7) . The flexibility analysis from PyMOL also indicates that all the ligands decrease the flexibility of SARS-CoV-2 3CL M pro upon binding ( Supplementary Figures 2A-D) . Inhibition Study of SARS-COV-2 3CL M pro by Other Researchers Teli et al. (2021) , have screened ten compounds namely, Procyanidin A3, Rutin, Solanine, Procyanidin A4, Procyanidin B4, Hypericin, Quercetagetin, Procyanidin, and Astragalin for inhibition of SARS-COV-2 3CL M pro . In that study they have shown that most of these compounds binds in the active site cavity of SARS-COV-2 3CL M pro (Teli et al., 2021) . Chourasia et al. (2020) have used a potential peptide (with amino acid sequence, KFVPKQPNMIL) from soy cheese for effective Inhibition of SARS-CoV-2 Main Protease and S1 Glycoprotein (Chourasia et al., 2020) . This peptide binds specifically to the amino acid residues that are important for the host cell entry and multiplication (3CL pro ) of SARS-CoV-2. Abdusalam and Murugaiyah (2020) have used ZINC database to identify zinc containing compounds as inhibitors of SARS-COV-2 3CL M pro , and have encountered four active zinc compounds (ZINC32960814, ZINC12006217, ZINC03231196, and ZINC33173588) which shows high binding affinity for 3CL pro pocket (Abdusalam and Murugaiyah, 2020) . Vincent et al. (2020) have used Compounds From Kabasura Kudineer on SARS-CoV-2 3CL pro and have shown that Acetoside, Luteolin 7, rutinoside, rutin, Chebulagic acid, Syrigaresinol, Acanthoside, Violanthin, Andrographidine C, myricetin, Gingerenone -A, Tinosporinone, Geraniol, Nootkatone, Asarianin, and Gamma sitosterol are the natural compounds in Kabasura Kudineer extracts, which can used as effective inhibitors against SARS-CoV-2 3CL pro (Vincent et al., 2020) . • The inhibition of 3CL M pro by natural (2S albumin and flocculating protein from M. oleifera) and synthetic inhibitor (Suramin) was demonstrated in this study. • The interaction between 3CL M pro and the inhibitors are largely through electrostatic force of attraction and with the interactions of amino acid residues from both sides. • All the three inhibitors bind between domain II and III (3CL M pro amino acid residues, Lys102, Pro108, Gln110, Asp155, Glu240, and His246, with Suramin and S139, T139, G302, Q299 with 2S albumin and flocculating protein. These interactions restrict the moment in domain III, which is important for dimerization and further for the function of SARS-COV2 3CL M pro . • Here we proposed that these inhibitors will inhibit 3CL M pro by preventing this enzyme from dimerization. • The current study will lead to the production of a new vaccine against COVID-19. The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author. AU designed the project and reviewed the manuscript. KU drafted and proofread the manuscript, and did English language corrections in the manuscript revision stage. Both authors contributed to the article and approved the submitted version. The Identification of Potential Inhibitors of 3CL Protease of SARS-CoV-2 From ZINC Database by Molecular Docking-Based Virtual Screening H++ 3.0: automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulation Epidemiological, demographic, and clinical characteristics of 47 cases of Middle East respiratory syndrome coronavirus disease from Saudi Arabia: a descriptive study Potential Drugs and Remedies for the Treatment of COVID-19: a Critical Review Activation of the SARS coronavirus spike protein via sequential proteolytic cleavage at two distinct sites GROMACS: A message-passing parallel molecular dynamics implementation Characterizatio,.n of an efficient coronavirus ribosomal frameshifting signal: requirement for an RNA pseudoknot Novel wuhan (2019-nCoV) coronavirus Features, evaluation and treatment coronavirus(COVID-19) The Amber biomolecular simulation programs A Potential Peptide From Soy Cheese Produced Using Lactobacillus delbrueckii WS4 for Effective Inhibition of SARS-CoV-2 Main Protease and S1 Glycoprotein SARS-CoV-2: Recent Reports on Antiviral Therapies Based on Lopinavir/Ritonavir, Darunavir/Umifenovir, Hydroxychloroquine, Remdesivir, Favipiravir and other Drugs for the Treatment of the New Coronavirus Particle mesh Ewald: an N log (N) method for Ewald sums in large systems The PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC AM1: a new general purpose quantum mechanical molecular model PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations The SARS-CoV-2 Spike Glycoprotein Biosynthesis, Structure, Function, and Antigenicity: Implications for the Design of Spike-Based Vaccine Immunogens Nidovirales: evolving the largest RNA virus genome Novel Wuhan (2019-NCoV) coronavirus ParDOCK: An All Atom Energy Based Monte Carlo Docking Protocol for Protein-Ligand Complexes Identification of severe acute respiratory syndrome coronavirus replicase products and characterization of papain-like protease activity Conservation of substrate specificities among coronavirus main proteases The COVID-19 Pandemic: A Comprehensive Review of Taxonomy, Genetics, Epidemiology, Diagnosis, Treatment, and Control Nucleotide sequence of the human coronavirus 229E RNA polymerase locus SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor MDWeb and MDMoby: an integrated web-based platform for molecular dynamics simulations Clinical features of patients infected with 2019 novel coronavirus in Wuhan Structural and functional properties of SARS-CoV-2 spike protein: potential antivirus drug development for COVID-19 The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak Fast, efficient generation of highquality atomic charges. AM1-BCC model: II. Parameterization and validation Therapeutic Strategies Against COVID-19 and Structural Characterization of SARS-CoV-2: A Structure of severe fever with thrombocytopenia syndrome virus nucleocapsid protein in complex with suramin reveals therapeutic potential K DEEP : Protein-Ligand Absolute Binding Affinity Prediction via 3D-Convolutional Neural Networks Structural plasticity of SARS-CoV-2 3CL M pro active site cavity revealed by room temperature X-ray crystallography Combating the Pandemic COVID-19 Repurposing Antiviral Protease Inhibitors Using Extracellular Vesicles for Potential Therapy of COVID-19 PDBsum: Structural summaries of PDB entries The race to make COVID antibody therapies cheaper and more potent A major outbreak of severe acute respiratory syndrome in Hong Kong Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB FireDock: a web server for fast interaction refinement in molecular docking Structure-based inhibition of Norovirus RNA-dependent RNA polymerases The trypanocidal drug suramin and other trypan blue mimetics are inhibitors of pyruvate kinases and bind to the adenosine site Autodock4 and AutoDockTools4: automated docking with selective receptor flexiblity Structures of MERS-CoV spike glycoprotein in complex withsialoside attachment receptors Nidovirus transcription: how to make sense Identification and characterization of severe acute respiratory syndrome coronavirus replicase proteins The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak Molecular aspects of COVID-19 differential pathogenesis Structural and functional characterization of suramin-bound MjTX-I from Bothrops moojeni suggests a particular myotoxic mechanism A contemporary view of coronavirus transcription PatchDock and SymmDock: servers for rigid and symmetric docking Airborne Transmission Route of COVID-19: Why 2 Meters/6 Feet of Inter-Personal Distance Could Not Be Enough In silico Screening of Natural Compounds as Potential Inhibitors of SARS-CoV-2 Main Protease and Spike RBD: Targets for COVID-19 Infectious RNA transcribed in vitro from a cDNA copy of the human coronavirus genome cloned in vaccinia virus Mechanisms and enzymes involved in SARS coronavirus genome expression Constant-pressure andtemperature molecular-dynamics simulations of crystals of the lecithinfragmentsglycerolphophorylcholine and dilauroylglycerol In situ structural analysis of SARS-CoV-2 spike reveals flexibility mediated by three hinges The Sequence and Three-Dimensional Structure Characterization of Snake Venom Phospholipases B. Front Crystal structure of mature 2S albumin from Moringa oleifera seeds Thrombinlike enzymes from snake venom: Structural characterization and mechanism of action The Sequence and a Three-Dimensional Structural Analysis Reveal Substrate Specificity Among Snake Venom Phosphodiesterases The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes Structural analysis of SARS-CoV-2 genome and predictions of the human interactome Molecular Docking Studies on the Anti-viral Effects of Compounds From Kabasura Kudineer on SARS-CoV-2 3CL pro Coronavirus biology and replication: implications for SARS-CoV-2 LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein Development and testing of a general AMBER force field Coupling Langevin Dynamics With Continuum Mechanics: Exposing the Role of Sarcomere Stretch Activation Mechanisms to Cardiac Function SWISS-MODEL: homology modelling of protein structures and complexes WHO/Europe | Coronavirus disease (COVID-19) outbreak -WHO announces COVID-19 outbreak a pandemic Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan novel antibody humanization method based on epitopes scanning and molecular dynamics simulation 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 Molecular biology of severe acute respiratory syndrome coronavirus Virus-encoded proteinases and proteolytic processing in the Nidovirales The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.Copyright © 2021 Ullah and Ullah. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.