key: cord-0759181-hb8psokr
authors: Mahmud, Shafi; Biswas, Suvro; Kumar Paul, Gobindo; Akter Mita, Mohasana; Afrose, Shamima; Robiul Hasan, Md.; Sharmin Sultana Shimu, Mst.; Abu Raihan Uddin, Mohammad; Salah Uddin, Md.; Zaman, Shahriar; Kaderia Kibria, K. M.; Arif Khan, Md.; Bin Emran, Talha; Abu Saleh, Md.
title: Antiviral peptides against the main protease of SARS-CoV-2: A molecular docking and dynamics study
date: 2021-07-14
journal: Arabian journal of chemistry
DOI: 10.1016/j.arabjc.2021.103315
sha: 8e4fd61838ddea5e3ca1ac127201764425505dd1
doc_id: 759181
cord_uid: hb8psokr

The recent coronavirus outbreak has changed the world’s economy and health sectors due to the high mortality and transmission rates. Because the development of new effective vaccines or treatments against the virus can take time, an urgent need exists for the rapid development and design of new drug candidates to combat this pathogen. Here, we obtained antiviral peptides obtained from the data repository of antimicrobial peptides (DRAMP) and screened their predicted tertiary structures for the ability to inhibit the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using multiple combinatorial docking programs, including PatchDock, FireDock, and ClusPro. The four best peptides, DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804, had binding energies of −1125.3, −1084.5, −1005.2, and −924.2 Kcal/mol, respectively, as determined using ClusPro, and binding energies of −55.37, −50.96, −49.25, −54.81 Kcal/mol, respectively, as determined using FireDock, which were better binding energy values than observed for other peptide molecules. These peptides were found to bind with the active cavity of the SARS-CoV-2 main protease; at Glu166, Cys145, Asn142, Phe140, and Met165, in addition to the substrate-binding sites, Domain 2 and Domain 3, whereas fewer interactions were observed with Domain 1. The docking studies were further confirmed by a molecular dynamics simulation study, in which several descriptors, including the root-mean-square difference (RMSD), root-mean-square fluctuation (RMSF), solvent-accessible surface area (SASA), radius of gyration (Rg), and hydrogen bond formation, confirmed the stable nature of the peptide–main protease complexes. Toxicity and allergenicity studies confirmed the non-allergenic nature of the peptides. This present study suggests that these identified antiviral peptide molecules might inhibit the main protease of SARS-CoV-2, although further wet-lab experiments remain necessary to verify these findings.

The current coronavirus disease 2019 pandemic has spread rapidly across the world, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which was first identified in December 2019 in Wuhan, Hubei Province, China Y. C. Wu et al., 2020; Zu et al., 2020) . COVID-19 has spread globally across 200 countries (Naqvi et al., 2020) , and the most affected countries are the USA, Brazil, India, Russia, Colombia, Spain, the UK, Peru, Argentina, Mexico, France, and Italy. According to the most recent update, released on 2 July 2021 by the World Health Organization (WHO), 182,319,261 confirmed cases and 3,954,324 deaths have been recorded worldwide (https://covid19.who.int/). The consequences of COVID-19 appear to be exceedingly severe among people older than 55 years (Emergency and Team, 2020; Guan et al., 2020; Report, 2020; Uddin et al., 2020) . By contrast, the fatality rate has been the lowest among individuals younger than 19 years, at 0%-0.1%, and a low fatality rate has been observed among those aged 20-54 years, at 0.1%-0.8%. The fatality rate increases gradually, ranging from 1.4% to 4.9% among those aged 55-74 years and reaching 4.3%-10.5% among those aged 75-84 years, whereas those older than 85 years are associated with the highest fatality rate of 10.4%-27.3% (Emergency and Team, 2020; Guan et al., 2020; Li et al., 2020; Report, 2020; Uddin et al., 2020) . Individuals with delicate immune systems or chronic diseases, such as diabetes, kidney disease, liver disease, malignant tumors, and cardiovascular disease, are at the highest risk of COVID-19 and experience a higher mortality rate (Emergency and Team, 2020; Gao et al., 2020; Guan et al., 2020; Li et al., 2020; Report, 2020; Uddin et al., 2020) .

Coronaviruses are single-stranded, positive-sense, enveloped ribonucleic acid viruses that are 60-140 nm in diameter and feature 9-12-nm long spike-like surface projections that manifest a crown or solar corona-like appearance (Singhal, 2020; Zu et al., 2020) . Coronaviruses have been categorized into four genera, including α-, β-, γ-, and δ-coronaviruses, and SARS-CoV-2 belongs to the β-coronavirus genus, of the order Nidovirales, the family Coronaviridae, and the subfamily Coronavirinae Naqvi et al., 2020; . Mammals can typically be infected by α-and β-coronaviruses, whereas avian species are typically infected by γ-coronaviruses, and both mammals and avians can be infected by δcoronaviruses (Naqvi et al., 2020) . SARS-CoV-2 represents the seventh coronavirus known to infect humans, including HCoV-229E, HKU-NL63, HCoV-OC43, HCoV-HKU1, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV) (Chan et al., 2015 (Chan et al., , 2013 Chen et al., 2020; Liu et al., 2020; .

The SARS-CoV-2 genome contains 13-15 open reading frames (ORFs), which are straddled by the 3′-untranslated region (UTR) and the 5′-UTR Mohammad et al., 2020) . The ORF 1ab encodes the pp1a and pp1ab polyproteins, which are subsequently cleaved by the main protease (Mpro) and the papain-like protease (PLpro) to form 16 nonstructural proteins (nsp1 to nsp16). Mpro, which is also known as 3CLpro, and PLpro play significant roles in the post-translational modification of the two replicase polyproteins, pp1a and pp1ab (Anand et al., 2003; Mohammad et al., 2020; Naqvi et al., 2020; Ton et al., 2020; Zhang et al., 2020) , and are required to generate functional subunits. These abundant functional subunits are required for the flexing and packaging of newly formed virions and broadly assist in the viral replication process, which is necessary to propagate the infection (Du et al., 2004; Hegyi and Ziebuhr, 2002; Naqvi et al., 2020; Ullrich and Nitsche, 2020; Ziebuhr et al., 2000) .

Mpro plays a central role in proteolysis, viral replication processes, transcriptional processes and is essential for the life cycle of the virus (Dai et al., 2020; Naqvi et al., 2020; . Because Mpro is the key enzyme necessary for SARS-CoV-2 replication and propagation, its inhibition can block these processes, preventing viral spread. Therefore, Mpro is considered to be a more promising drug target than PLpro (Estrada, 2020; Zhang et al., 2020) . Some drugs have been examined explored for their ability to be repurposed for use in the treatment of COVID-19. Among these, remdesivir and favipiravir have been identified as RNA-dependent RNA polymerase (RdRp)-targeted drugs that exhibit antiviral activity against SARS-CoV-2. In addition, antimalarial drugs, including hydroxychloroquine and chloroquine, have been demonstrated to suppress SARS-CoV-2 replication (Kiplin Guy et al., 2020) . Other antimalarial drugs, including halofantrine, doxycycline, mefloquine, lumefantrine, atovaquone, primaquine, and sulfonamides, have also been used to treat COVID-19 (Schlagenhauf et al., 2020) . Typically, vaccines are required to successfully demonstrate efficacy and safety in three clinical trial phases (I, II, and III). As of 2 July 2021, according to the most recent update reported by the WHO, approximately 50, 37, and 32 vaccine candidates are currently undergoing Phase I, II, and III clinical trials, respectively (https://www.nytimes.com/interactive/2020/science/coronavirus-vaccine-tracker.html).

To date, several vaccines have been approved for use in different countries. On 31 December 2020, a COVID-19 vaccine named 'BNT162b2/COMIRNATY Tozinameran (INN)', manufactured by Pfizer, was granted an Emergency Use Listing (EUL) by the WHO. Subsequently, the WHO permitted the EUL use of two versions of the vaccine named 'AZD1222', manufactured by AstraZeneca and SKBio. In addition, the COVID-19 vaccine named 'Covishield' was manufactured by the Serum Institute of India, and the 'mRNA-1273' vaccine, manufactured by Moderna, became available near the end of February 2021. Two additional vaccine manufacturers, Sinopharm and Sinovac, launched two identically named vaccines, called 'SARS-CoV-2 Vaccine,' in early March. Additionally, on 12 March 2021, Janssen (Johnson & Johnson) developed the 'Ad26.COV2.S' vaccine. Besides, 'Sputnik V,' 'Ad5-nCoV', and 'EpiVacCorona' vaccines were developed by various manufacturers and have been issued EULs by the WHO (https://www.who.int/news-room/q-adetail/coronavirus-disease-(covid-19)-;(Organización Mundial de la Salud, 2021)). Moreover, on March 15, China endorsed a vaccine for exigency use named "ZF2001" which has been made by two companies entitled-Anhui Zhifei Longcom and the Institute of Medical Biology at the Chinese Academy of Medical Sciences as copartners. As research continues, "NVX-CoV2373" is a vaccine developed by the company Novavax which is Maryland-based. Besides, the Beijing Institute of Biological Products manufactured the "BBIBP-CorV" vaccine which is authorized as exigence use by WHO on May 7. Nevertheless, the "CoronaVac" vaccine is developed by Sinovac Biotech, an unofficial Chinese company and on June 1 it gets permission for emergency use by WHO (https://www.nytimes.com/interactive/2020/science/coronavirus-vaccine-tracker.html).

Peptides allocate themselves as fastidious, spiffing, benign signaling molecules that can able to affix to distinct cell surface receptors. Having these potential and underlying characteristics, peptides personate as novel therapeutics (Fosgerau and Hoffmann, 2015) . As a drug, peptides propound sundry precedence comprehend with exuberant biological alacrity, bluffy specificity, and mangy toxicity (Bruckdorfer et al., 2005; Lien and Lowman, 2003) . The comprehensive obstacle of being peptides therapeutics are petty active period due to impetuous renal clearance, scantiness of selective receptor subtype, minuscule half-life span. Peptides accommodate exorbitant affinity or specificity to a particular target and beneath toxic silhouette which makes assorted superiority over small molecules (Nestor Jr., 2009; Sato et al., 2006) . US Food and Drug Administration (FDA) endorse peptide therapeutics against the Hepatitis C virus in 2011 and also ratify as the first peptide drug of Human immunodeficiency virus (HIV) for clinical attestation. In addition, several peptides are undergoing phase trials against the Influenza virus, Hepatitis B virus, and Hepatitis D virus. Contrariwise, Chikungunya virus (CHIKV), Zika virus, and Dengue virus (DENV) have peptide therapeutics exhibiting antiviral activity (Agarwal and Gabrani, 2021) .

With the aim of identifying plausible antiviral therapeutic agents to treat SARS-CoV-2, researchers all over the world have conducted virtual screening strategies that include molecular docking (Abdusalam and Murugaiyah, 2020; Abo-zeid et al., 2020; Cheke, 2020; Hagar et al., 2020; Joshi et al., 2020; Khaerunnisa et al., 2020; Nag et al., 2021; Rutwick Surya and Praveen, 2021) or molecular dynamics simulations combined with molecular docking (Al-Karmalawy et al., 2021; Chikhale et al., 2021; Das et al., 2020; Keretsu et al., 2020; S. Kumar et al., 2020; Majumder and Mandal, 2020; Mhatre et al., 2020; Sepay et al., 2021; Teli et al., 2021) . SARS-CoV-2 has exhibited almost 1,516 nucleotide-level variations across various genome-wide annotations, and due to frequent adaptive mutations through which the viral genome evolves viral resistance, no effective or plausible antiviral treatment has yet been developed (M. R. Xu et al., 2020) . However, peptide-oriented molecular docking and molecular dynamics techniques can provide a rapid platform for reliable drug development in the race to combat SARS-CoV-2.

Potential antiviral peptides were retrieved from the data repository of antimicrobial peptides (DRAMP), (Kang et al., 2019) which contains 22,209 entries, including 5,841 antimicrobial peptides. In our study, we only examined those peptide molecules with known antiviral activity. Based on antiviral activity, we screened 215 antiviral peptides from 5820 peptides. Additional screening procedures were executed according to the following parameters. First, we eliminated 62 antiviral peptides longer than 50 amino acids in length to facilitate the application of a new approach called PEP-FOLD 3.5 (Maupetit et al., 2010) , which can accurately predict the structures of peptides containing from 5 to 50 amino acids. Two additional antiviral peptides with unknown amino acids in their sequences, denoted as 'X,' were also eliminated. Approximately 31 antiviral peptides were removed due to a lack of authentic information regarding their activity or function within the database. Another 16 antiviral peptides were removed due to a lack of any plausible antiviral activity or that presented with inefficient antiviral activities and were therefore not suitable for combating any viral disease. Consequently, among the initially identified 5,841 antimicrobial peptides, we analyzed 104 antiviral peptides, each consisting of fewer than 50 amino acids and with evidence of antiviral activity, which were identified by an in-depth analysis of the database. The peptide sequence information was used to estimate the peptide structures using the PEP-FOLD 3.5 webserver. This webserver predicts the peptide structure using a Hidden Markov Model suboptimal sampling algorithm. The energy minimizations of the peptides were implemented in Avogadro software with the aid of MMFF94 force field, using the steepest descent algorithm and 500 steps.

The three-dimensional structure of the SARS-CoV-2 Mpro (PDB ID: 6LU7) was retrieved from the Protein Data Bank (Protein Data Bank, 2019) with a resolution of 2.16 Å. The protein structure was initially prepared in Discovery Studio (San Diego: Accelrys Software Inc., 2012), where the water molecules and heteroatoms were removed. The cleaned protein structure was further minimized in YASARA tools by employing the AMBER14 force field (Land and Humble, 2018) . The minimized protein structure was used for further molecular docking and dynamics studies.

Peptide and protein docking was implemented as described in our previous study (Mahmud, S., Paul, G. K., Biswas, S., Afrose, S., Mita, M. A., Hasan, M., ... & Saleh, 2021). The protein structure of Mpro was input as the receptor molecule, and the peptides were used as the ligands in the PatchDock (Schneidman-Duhovny et al., 2005) webserver. The final clustering was selected on the basis of the root-mean-square deviation (RMSD) value. The top ten solutions were taken from PatchDock and further docked using FireDock (Mashiach et al., 2008) tools. The FireDock program optimizes lateral string conformations and rigid body formations to provide larger refinement. These web tools were employed to obtain rigid protein-peptide docking. The best protein-peptide complex was selected from among the top ten conformers based on energy scoring. All the 104 peptides were then further docked using the ClusPro program (Comeau et al., 2004) to obtain more accurate binding energy profiles from the corresponding peptide and protein complexes. Based on the combination of multiple docking programs, four peptide-protein complexes were selected for non-bonded structure analysis. The interaction analysis was conducted using PyMOL (DeLano, 2002) and Discovery Studio (San Diego: Accelrys Software Inc., 2012) software package.

The docked peptide-protein molecules and their conformational variability were assessed through molecular dynamics simulations in YASARA dynamics software (Krieger et al., 2004) . The AMBER14 force field (Maier et al., 2015) was used in this study, and the docked complexes were initially cleaned, along with hydrogen bond orientation and optimization for simulation. The cubic simulation box with a cell size of 110×110×110Å 3 was filled with water molecules, with a 0.9899 g/cm 3 density, and the system was neutralized with 0.9% NaCl (Krieger, Elmar and Spronk, 2013) . The TIP3P water model was used to solvate the complex. The acid dissociation constant value (pKa) was calculated for the amino acids in the complex. The SCWRL algorithm combined with hydrogen bond network optimization was applied to maintain the correct protonation state of each amino acid residue. The particle mesh Ewald method (Darden et al., 1993) was applied to calculate long-range electrostatic interactions, using a cutoff radius of 8 Å. The simulation cell box was 20 Å larger than the peptide-protein complex to allow for maximum free motion, and the periodic boundary condition was maintained (Krieger and Vriend, 2015) . The energy of each simulated system was minimized by using the steepest descent gradient approaches (5,000 cycles) via simulated annealing methods. An 80-ps position-restrained molecular dynamics simulation was performed for each system at constant pressure (1 atm) and temperature (310 K). Temperature coupling was performed using the Berendsen thermostat, with a temperature coupling constant of 0.1 ps, and a manometer method was used for pressure coupling, with a reference pressure of 1 atm. For the simulation study, a time step of 1.25 fs was used. The chemical bond lengths involving hydrogen bond atoms were fixed using the SHAKE algorithms (Brooks et al., 1983) . After reaching an equilibrium state at 1 ns, the simulation was allowed to run for 250 ns, and the simulation trajectories were saved after every 100 ps. The simulation was conducted three times, and the average values from the trajectories were used to analyze the RMSD, root-mean-square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), and hydrogen bonds (Bappy et al., 2020; Khan et al., 2020; Swargiary et al., 2020) .

Furthermore, the simulation snapshots were subjected to binding free energy calculations from MM-PBSA approaches by following equations (Mitra and Dash, 2018) .

The YASARA macro was used to calculate the MM-PBSA binding energy where positive energy indicates the better bindings (Srinivasan and Rajasekaran, 2016) .

The physiochemical properties of the four best peptides were evaluated using ProtParam (Gasteiger et al., 2005) tools. This tool provides relevant properties, including molecular weight, net charge at pH 7, volume, peptide properties, stability, and charge.

The allergenicity and toxicity of the peptides were calculated using AllergenFP (Dimitrov et al., 2014b) and AllerTOP (Dimitrov et al., 2014a) webservers. The AllergenFP webserver predicts peptide allergenicity by employing a five E-descriptor-based fingerprinting method, whereas the AllerTOP webserver uses amino acid E-descriptors and k-nearest neighbor (kNN) machine learning approaches.

Molecular docking approaches can be used to identify the binding affinities of various ligands for the target protein structure to assist in the performance of rational drug design (Jones and Willett, 1995) . This approach can contribute toward understanding the interaction dynamics and potential binding mechanisms, which can be applied for more rigorous inhibition (Liu et al., 2018) . By combining several docking algorithms and programs, our approach provided enhanced structural accuracy and the precise calculation of binding interactions and affinities (Bartuzi et al., 2017) . The peptide model structure derived from PEP-FOLD is presented in Figure S1 . The structures of 104 total peptides were modeled, and all 104 peptide structures were used to dock against the target Mpro enzyme. Howbeit, we selected the top ten elevated peptide molecules in conformity with the combined docking scores procured from PatchDock, ClusPro, and FireDock (Table 1, Table S1 ). Furthermore, we opted top four plausible peptides amid the top ten peptide molecules predicated on the calculated binding energy for non-bonded structure analysis (Figure 1) The non-bonded interactions between the peptide molecules and the SARS-CoV-2 Mpro are shown in figure 2. The DRAMP00877 peptide and Mpro complex was stabilized by five hydrogen bonds at GLN107, ASP245, GLU240, ILE106, and PRO108, five hydrophobic bonds at VAL140, VAL297, HIS246, PHE294, and ILE249, two electrostatic bonds at ASP153 and GLU240 in addition to two unfavorable bonds at ASP248 and GLN110 ( Figure  2 and Table 2 ). The DRAMP02333 peptide and Mpro enzyme formed four hydrogen bonds at LYS137, LYS236, SER123, and TYR237, one electrostatic interaction at ASP197, and five hydrophobic bonds at ALA7, VAL125, CYS128, LEU286, and LYS5 together with four unfavorable bonds at TYR126, GLN127, ASP289, and ARG131 ( Figure 2 and Table 2 ). The DRAMP02669 peptide and Mpro complex formed five hydrophobic bonds at PHE294, VAL297, VAL104, ILE106 and PRO252, four hydrogen bonds at LYS100, LYS102, PRO108 and ASP245 and lastly four electrostatic bonds at LYS100, ASP153, ASP248 and ASP245. (Figure 2 and Table 2 ). The DRAMP03804 peptide and Mpro complex created nine hydrogen bonds at ASN142, CYS145, GLN189, THR24, THR25, SER144, SER46, HIS172 and GLU166 and five hydrophobic bonds at MET165, Pro168, LEU141, Met49, and Leu50 and two electrostatic bonds at HIS172 and GLU166 (Figure 2 and Table 2 ). The DRAMP03804 peptide binds in the active groove of Mpro, at Glu166, Cys145, Asn142, and Met165 . Binding to the active sites may lead to the possible inhibition of the target molecule (Daddam et al., 2020) .

The binding interaction and the conformational variations that occur during binding can be understood using a docking study to provide insight into the binding sites to a limited extent. Further changes in the ligand or the bound complex can be studied by varying the temperature and pressure in a molecular dynamics simulation, which can minimize the cost of performing such experiments in the laboratory (Arfin et al., 2018; Durrant and McCammon, 2011; Nair and Miners, 2014) . A combined docking and molecular dynamics approach can be used to validate the docking-derived results (Gioia et al., 2017) .

Molecular dynamics simulations were implemented for the peptide and SARS-CoV-2 Mpro complexes to confirm the structural rigidity and validate the docking outcomes for the complexes. The RMSD values of the C-alpha atoms were explored to understand the structural rigidity. Figure 3 (a) demonstrated that the DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804 complexes showed initial RMSD increases due to instability. The increase in RMSD was higher for the DRAMP00877 peptide complex than for the other complexes. The upward trend in RMSD was maintained until 40 ns, after which all four complexes achieved stability, which was maintained throughout the remainder of the simulation period. Although the DRAMP00877 complex had a higher RMSD trend during the initial phase, the RMSD profile decreased from 40-100 ns. The DRAMP02333 complex had a lower RMSD than the other complexes, indicating increased structural integrity. The average RMSD profile of the complexes formed with DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804 were 1.88 Å, 1.58 Å, 1.82 Å, and 1.81 Å, respectively. The overall RMSD trend did not exceed 2.5 Å, which indicates overall structural stability as higher RMSD profiles for C-alpha atoms are indicators of low stability .

The SASA values of the complexes were analyzed to evaluate changes in the surface of the SARS-CoV-2 Mpro in response to binding with the peptide molecules. The DRAMP00877, DRAMP02333, and DRAMP03804 complex had similar SASA trends throughout the entire simulation trajectory, indicating no change in protein volume. However, the DRAMP02669 complex showed an initial upward trend in SASA until 50 ns, indicating an expansion of the surface area following initial complex formation [ Figure 3 (b) ]. The complex then stabilized and showed a trend similar to that of the DRAMP02333 peptide complex. The average SASA profiles of the four complexes were 15663.42Å 2 , 16302.79Å 2 ,16066.2Å 2 and 15590.45Å 2 for DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804, respectively. The SASA profiles of the peptide-protein complexes had similarly stable trends, with little fluctuation, indicating a lack of expansion or contraction for the protein complexes during the simulation time (Rakib et al., 2021) .

The compactness of the protein complexes was evaluated using Rg, for which higher values indicate a more labile nature and a lower value indicates the firmness of the protein. The DRAMP02669 complex displayed fluctuations in the Rg value until 20 ns, indicating an initially loose packaging system for the protein [ Figure 3 (c) ]. However, this complex maintained stability for the remainder of the simulation time. The DRAMP03804 complex had a larger Rg trend than all other complexes but was associated with a lower degree of deviation, demonstrating the achievement of a steady state. The comparative higher Rg of this peptide may be responsible for relatively more flexible state and mobile nature of the complexes. The average Rg values were 22.504, 22.527, 22.501, and 22 .867 Å for DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804, respectively. The higher Rg profile of the DRAMP03804 complex than other complexes may defines the more mobile nature and folding and unfolding of the protein complexes. In addition, hydrogen bond formation between the protein and peptide in complex plays a vital role in maintaining complex stability. Figure 3 (d) demonstrated that the DRAMP02333 complex contained the highest number of hydrogen bonds throughout the entire simulation time. All of the evaluated peptide and Mpro complexes were stable, and high levels of fluctuation were generally absent for these complexes.

The RMSF of the peptide and Mpro complexes were explored to understand the flexibility across the amino acid region of Mpro. Figure 3 (e) demonstrated that almost all residues had low RMSF, except Ser1 (helix-strand), Glu47 (beta-turn), Asp48 (beta-turn), Arg60 (helixstrand), Lys97 (beta-turn), Lys137 (beta-turn), Arg222 (beta-turn), Phe294 (helix-strand), Ser301 (beta-turn), Gly302 (beta-turn), Val303 (beta-turn), Thr304 (beta-turn), Phe305 (betaturn), and Gln306 (beta-turn).

The pre-and post-molecular dynamics docked structures were superimposed to explore variations among these structures. The RMSD of the DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804 peptides were found to be 2.076, 1.658. 1.793, and 1.778 Å, respectively, indicating limited flexibility among the structures (Figure 4) . Moreover, the binding free energy of the peptide-protein complexes were calculated from the MM-PBSA method ( Figure 5 ) where more positive score indicates the better binding of the complexes. The average free energy of the DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804 peptide complexes were -245.612, -243.441, 19.51237, -51 .5827KJ/mol respectively. The DRAMP02669 peptide complex exhibited more favorable binding than the other complexes.

The flexibility of the amino acid residues within the peptides was also assessed through RMSF analysis, and all four peptide molecules were found to display RMSF values below 2.5 Å, except Lys23 in the DRAMP00877 peptide; Tyr29, Asp33, and Val34 in the DRAMP02333 peptide; Glu13, Gly17, Thr18, Tyr21, and Gly22 in the DRAMP02669 peptide; and Val1, Val2, and Arg16 in the DRAMP03804 peptide. Lower RMSF values among the amino acid sequence indicate a less flexible and more stable nature of the peptide molecule ( Figure S6 ). The secondary structures of the peptide molecules were assessed to explore changes in the secondary structures, included the formation of α helices, turns, random coils, 310 helices, and pi helices ( Figure S7 ). The DRAMP00877 peptide displayed stable proportions of α-helices, β-sheets, turns, random coils, and pi helices, although some deviations were observed for 310 helices; however, 310 helices represented a lower percentage of secondary structure contents for the DRAMP00877 peptide. Random coils and helices were primarily noted for the DRAMP02333 peptide, which displayed a stable trend throughout the entire simulation times, with minor deviations in the proportions of 310 helices and turns. Random coils were primarily noted for the DRAMP02669 peptide, and stable profiles were observed for all secondary structures observed for this peptide throughout the simulation. Turns and random coils were primarily observed for the DRAMP03804 peptide, which featured stable trends, with some aberrations for α helix and 310 helix formation.

The SARS-CoV-2 Mpro consists of three domains; Domain 1, Domain 2, and Domain 3. Domain 1 contains residues 8-101; Domain 2 is formed by residues 102-184, which features an antiparallel β-barrel structure; and Domain 3 consists of residues 201-303 . Domain III (residues 201-303) comprises of five α-helices and is liable for the dimerization of enzymes exhibiting involvement in dimer formation Kneller et al., 2020) . In the Mpro dimerization, the troublesome the foreword of domain III has been adjusted by fragment elimination demonstrating that a decollated enzyme inexistent domain III stays as a monomer which is catalytically passive (Suárez and Díaz, 2020) . Thus, Mpro constructs a functioning dimer via intermolecular interactivity, predominantly mesne the helical domains eliciting indispensable accosts of Domain III in the protease functionality (Kneller et al., 2020) . The substrate-binding position is residing in the cleft interim domains I and II and the protomers, that cohere with every other via N-terminus residues 1-7, which are lying middle domains II and III with playing roles in the formation of the substrate-binding site. Four subsites namely; S1', S1, S2, and S4 comprised the substrate-binding cleft (Mengist et al., 2021) . The S1 subsite is formed with Phe140, Gly143, His163, His172, Cys145 and Glu166 while S2 comprised with Thr25, His41, and Cys145 amino acid residues; chiefly inlaid in electrostatic and hydrophobic interactions. The perfunctory subsites S3-S5 consist of Met49, His41, Met165, Glu166, and Gln189 amino acid residues. These shallow subsites can endure varied performance Yang et al., 2005) .

The binding interactions of the docked complexes were analyzed after taking snapshots at the 100, 200, and 250ns time points during the simulation to understand interaction stability and deviations that occur across the simulation time. The DRAMP00877 peptide had several stable interactions at Domain 3 (Val 297, Phe294, His246, Glu240, Ile249, and Asp248), and equal stable interactions were found with Domain 2 (Gln107, Ile106, Pro108, Val104, Asp153, and Gln110) after 100 ns of simulation (Table S2 , Figure S2 ). Meanwhile, the DRAMP00877 peptide had the most stable interactions at Domain 3 (Phe294, Glu240, Ile249, His246, Asp248, and Asp245), and two stable interactions were observed with Domain 2 (Gln110 and Asp153) after 200 ns of simulation. Howbeit, the DRAMP00877 peptide had two stable interactions at Domain 2 (Gln110 and Asp153) and diverse stable interactions at Domain 3(Phe294, Ile249, Glu240, His246, and Asp245).

In case of the DRAMP02333 peptide molecule, there were diverse interactions in the simulation environment at Domain 3 (Lys236, Asn238, Leu287, Tyr237, Asp289 and Leu286) and four stable interactions at Domain 2 (Lys137, Cys128, Tyr126 and Ser123 in conjunction with an interaction at Domain 1 (Ala7). Domain2 and Domain 3 of Mpro are connected by a long loop region, and DRAMP02333 formed stable contacts with Asp197 within the loop region. Furthermore, the DRAMP02333 peptide molecule exhibited several interactions at Domain 2 (Val125, Gly124, Ser139, Lys137, Ser123, and Tyr126) and three stable interactions at Domain 3 (Asn238, Lys236, and Leu286) in tandem with an interaction at Domain 1 (Glu14 and Ala7) in the simulation environment after 200 ns of simulation (Table S3 , Figure S3 ), The DRAMP02333 constituted stable contacts with Asp197 within the loop region that connects Domain 2 and Domain 3 of Mpro. Howbeit, after 250 ns of simulation, the DRAMP02333 peptide molecule demonstrated various interactions in the simulation environment, including nine static interactions were noticed at Domain 2 (Val125, Lys137, Ser139, Gly124, Gly138, Tyr126, Pro122, Ser123 and Glu166) and six stable interactions at Domain 3 (Asn238, Lys236, Tyr237, Leu286, Leu287, and Met235) in conjunction with two rigid interactions at Domain 1 (Glu14, and Ala7). Here, the DRAMP02333 formed stable contacts with Thr199 within the loop region connecting Domain2 and Domain 3 of Mpro simulation.

For the DRAMP02669 peptide, sole one interaction audited at Domain 1 (Lys100), in addition with six interactions, which observed at Domain 2 (Lys102, Ile106, Val104, Gln110, Pro108, Asp153) and five interactions observed at Domain 3 (Phe294, Val297, Pro252, Asp248, Asp245) after the 100ns simulation (Table S4, Figure S4 ). Additionally, four interactions noticed at Domain 2 (Lys102, Val104, Gln107, Asp153), with five interactions observed at Domain 3 (Asp245, His246, Asp248, Phe294, Val297) after the 200ns simulation. Besides, after the 250ns simulation only three interactions were seen at Domain 2 (Lys102, Val104, Asp153), with notably seven interactions observed at Domain 3 (Pro252, His246, Phe294, Val297, Asp248, Ile249, Asp245).

There are six interactions observed for the DRAMP03804 peptide were formed with Domain 1 (His41, Ser46, Leu50, Thr25, Thr24, Met49), with six interactions observed at Domain 2 (Asn142, Met165, Pro168, Glu166, Phe140, Leu141) after the 100 ns simulation (Table S5 , Figure S5 ). Moreover, afterwards the 200ns simulation, five interactions were observed at Domain 1 (Ser46, Met49, Leu50, Thr25, Thr24), including seven interactions observed at Domain 2 (Asn142, Met165, Pro168, Leu141, Phe140, Glu166, His172) . Furthermore, only three interactions were observed at Domain 1 (Cys44, Ser46, Leu50), with six interactions observed at Domain 2 (Asn142, Glu166, Met165, Pro168, Leu141, His163) after the 250ns simulation and this peptide constructed static contacts with GLN 189, ALA 191 among the loop region that appends Domain 2 and Domain 3 of Mpro. The interactions formed by the DRAMP03804 peptide primarily involve the active residues of the SARS-CoV-2 Mpro, which may interfere with the Mpro function. The co-crystalized ligand obtained with the Mpro structure binds in the substrate-binding pocket in extended conformations, and the inhibitor backbone atoms from an antiparallel sheet that interacts with residues 164-168 . Interestingly, the DRAMP03804 peptide molecules also featured three conserved interactions within these regions, at Met165, Pro168 and Glu166. It is cognizant that, DRAMP03804 peptide has Asn142, Gln189, Met165, Pro168, Leu141 and Glu166 residues that elicit resemblance with sundry of the abuzz residues of SARS-CoV-2 Mpro .

The peptide molecule DRAMP00877 interacted with the SARS-CoV-2 Mpro through conserved residues Phe294, His246, Glu240, Ile249, Asp153, and Gln110, whereas the DRAMP02333 peptide interaction involved Lys137, Tyr126, Ser123, Lys236, Ala7 and Leu286 throughout the simulation at 0ns, 100ns, 200ns and 250ns. Moreover, in all respects of the entire simulation trajectory, the amino acid residues Lys102, Phe294, Asp153, Val297, Val104, Asp248 and Asp245 of DRAMP02669 were involved in the interaction with SARS-CoV-2 Mpro and the Asn142, Gln189, Met165, Pro168, Leu141, Ser46, Glu166 and Leu50 residues contributed to the interaction between Mpro and DRAMP03804.

An allergenic antigen can trigger the activity of the Th2 response, causing B cells to produce IgE, which binds with the receptor molecule FcεRI and activates eosinophils. Activated eosinophils can lead to inflammation and tissue damage (Brock, 1995; Dimitrov et al., 2013; Huby et al., 2000) . The AllerTOP tools predict allergens based on descriptors associated with amino acid properties and has been demonstrated to predict allergenicity with a high level of accuracy and robustness (Dimitrov et al., 2014a) . The toxicity and allergenic profiles of the peptide molecules were evaluated using multiple web tools (Table 3 ). The DRAMP00877, DRAMP03804, and DRAMP02333 peptide molecules were non-allergic and were characterized as non-toxic, whereas the DRAMP02669 peptide was identified as a probable allergen in the AllerTOP tool, although it was reported as a non-allergen by the AllergenFP webserver.

We assessed the physicochemical properties of the best four peptide molecules ( Table 4 ). The largest peptide contained 44 amino acids, whereas the smallest peptide only contained 30 amino acids. The molecular weights and theoretical isoelectric point (pI) were higher for the DRAMP02333 and DRAMP03804 peptides than for the other peptides. The peptides DRAMP00877, DRAMP02333, DRAMP02669, and DRAMP03804 contained 1,5,1, and 1 negatively charged residues, respectively. The total number of atoms included in the entire peptide molecule ranged from 400-600, except for DRAMP02333.

Antiviral peptides can serve as effective leads for the development of new therapeutic options against SARS-CoV-2. In this study, peptide molecules were tested against the SARS-CoV-2 Mpro, and the four best peptides were identified based on the binding affinities with SARS-CoV-2 Mpro. This study revealed that the best peptide molecules bind to the active amino acid residues in Mpro. The molecular dynamics simulation study confirmed the docked conformations and structural properties that were identified in the rigid state. In addition, allergenicity profiling confirmed the non-allergic properties of the peptide molecules selected in this current study. The present study may aid the development of effective drugs against this deadly virus. Author Agreement Statement I the undersigned on behalf of the all authors, declare that this manuscript is original, has not been published before and is not currently being considered for publication elsewhere.

I confirm that the manuscript has been read and approved by all named authors and that there are no other persons who satisfied the criteria for authorship but are not listed. I further confirm that the order of authors listed in the manuscript has been approved by all of us.

Identification of Potential Inhibitors of 3CL Protease of SARS-CoV-2 From ZINC Database by Molecular Docking-Based Virtual Screening

A molecular docking study repurposes FDA approved iron oxide nanoparticles to treat and control COVID-19 infection

Antiviral Peptides: Identification and Validation

Molecular docking and dynamics simulations reveal the potential of anti-HCV drugs to inhibit COVID-19 main protease

Coronavirus main proteinase (3CLpro) Structure: Basis for design of anti-SARS drugs. Science (80-. )

Flexible docking-based molecular dynamics simulation of natural product compounds and Ebola virus Nucleocapsid ( EBOV NP ): a computational approach to discover new drug for combating Ebola

Extensive immunoinformatics study for the prediction of novel peptide-based epitope vaccine with docking confirmation against envelope protein of Chikungunya virus: a computational biology approach

Recent advances and applications of molecular docking to g protein-coupled receptors

Immunobiology: the immune system in health and disease

CHARMM: A program for macromolecular energy, minimization, and dynamics calculations

From Production of Peptides in Milligram Amounts for Research to Multi-Tons Quantities for Drugs of the Future

Middle East Respiratory syndrome coronavirus: Another zoonotic betacoronavirus causing SARS-like disease

Interspecies transmission and emergence of novel viruses: Lessons from bats and birds

The Molecular Docking Study of Potential Drug Candidates Showing Anti-COVID-19 Activity by Exploring of Therapeutic Targets of SARS-CoV-2

Emerging coronaviruses: Genome structure, replication, and pathogenesis

Computational assessment of saikosaponins as adjuvant treatment for COVID-19: molecular docking, dynamics, and network pharmacology analysis

Quaternary structure of the severe acute respiratory syndrome (SARS) coronavirus main protease

ClusPro: An automated docking and discrimination method for the prediction of protein complexes

Designing, docking and molecular dynamics simulation studies of novel cloperastine analogues as anti-allergic agents: Homology modeling and active site prediction for the human histamine H1 receptor

Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease

Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems

In-Silico approach for identification of effective and stable inhibitors for COVID-19 main protease (Mpro) from flavonoid based phytochemical constituents of Calendula officinalis

Pymol: An open-source molecular graphics tool. CCP4 Newsl

AllerTOP v.2 -A server for in silico prediction of allergens

AllerTOP -a server for in silico prediction of allergens

AllergenFP: Allergenicity prediction by descriptor fingerprints

Polyprotein cleavage mechanism of SARS CoV M pro and chemical modification of the octapeptide

Molecular dynamics simulations and drug discovery

Osong Public Health and Research Perspectives Early Epidemiological and Clinical Characteristics of 28 Cases of

Topological analysis of SARS CoV-2 main protease

Peptide therapeutics: Current status and future directions

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

The Proteomics Protocols Handbook

Dynamic docking: A paradigm shift in computational drug discovery

Clinical characteristics of coronavirus disease 2019 in China

Investigation of some antiviral N-heterocycles as COVID 19 drug: Molecular docking and DFT calculations

Conservation of substrate specificities among coronavirus main proteases

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Why are some proteins allergens?

Genome-wide analysis of SARS-CoV-2 virus strains circulating worldwide implicates heterogeneity

Identification and in silico molecular modelling study of newly isolated Bacillus subtilis SI-18 strain against S9 protein of Rhizoctonia solani

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

Docking small-molecule ligands into active sites

In silico screening of natural compounds against COVID-19 by targeting Mpro and ACE2 using molecular docking

DRAMP 2.0, an updated data repository of antimicrobial peptides

Rational approach toward COVID-19 main protease inhibitors via molecular docking, molecular dynamics simulation and free energy calculation

Potential Inhibitor of COVID-19 Main Protease ( M pro ) from Several Medicinal Plant Compounds by Molecular Docking Study

Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: a molecular docking study

Rapid repurposing of drugs for COVID-19. Science (80-. )

Structural plasticity of SARS-CoV-2 3CL Mpro active site cavity revealed by room temperature X-ray crystallography

YASARA -Yet Another Scientific Artificial Reality Application

Making optimal use of empirical energy functions: Force-field parameterization in crystal space

New ways to boost molecular dynamics simulations

Promising inhibitors of main protease of novel corona virus to prevent the spread of COVID-19 using docking and molecular dynamics simulation

Discovery of New Hydroxyethylamine Analogs against 3CLproProtein Target of SARS-CoV-2: Molecular Docking, Molecular Dynamics Simulation, and Structure-Activity Relationship Studies

YASARA: A tool to obtain structural guidance in biocatalytic investigations

Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science (80-. )

Therapeutic peptides

COVID-19: The first documented coronavirus pandemic in history

Application of molecular docking for the degradation of organic pollutants in the environmental remediation: A review

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

Prospective role of peptide-based antiviral therapy against the main protease of SARS-CoV-2. Front

Virtual screening and molecular dynamics simulation study of plant-derived compounds to identify potential inhibitors of main protease from SARS-CoV-2

ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB

Screening of plant-based natural compounds as a potential COVID-19 main protease inhibitor: an in silico docking and molecular dynamics simulation approach

FireDock: a web server for fast interaction refinement in molecular docking

A fast method for large-scale de novo peptide and miniprotein structure prediction

Structural Basis of Potential Inhibitors Targeting SARS-CoV-2 Main Protease

Elsevier has created a COVID-19 resource centre with free information in English and Mandarin on the novel coronavirus COVID-19 . The COVID-19 resource centre is hosted on Elsevier Connect , the company ' s public news and information

Structural dynamics and quantum mechanical aspects of shikonin derivatives as CREBBP bromodomain inhibitors

Identification of high-affinity inhibitors of SARS-CoV-2 main protease: Towards the development of effective COVID-19 therapy

Phytochemicals as potential drug candidates for targeting SARS CoV 2 proteins, an in silico study

Molecular dynamics simulations: from structure function relationships to drug discovery

Insights into SARS-CoV-2 genome, structure, evolution, pathogenesis and therapies: Structural genomics approach

The Medicinal Chemistry of Peptides

Status of COVID-19 Vaccines within WHO EUL/PQ evaluation process

RCSB PDB: Homepage. Rcsb Pdb

A molecular modelling approach for identifying antiviral selenium-containing heterocyclic compounds that inhibit the main protease of SARS-CoV-2: An in silico investigation

Bialek -2020 Severe-outcomes-among-patients-withcoronavirus-disease-2019-COVID19--United-States-February-12march-16-20202020Morbidity-and-Mortality-Weekly-ReportOpen-Access

A molecular docking study of SARS-CoV-2 main protease against phytochemicals of Boerhavia diffusa Linn. for novel COVID-19 drug discovery

Discovery Studio Modeling Environment

Therapeutic peptides: technological advances driving peptides into development

Repurposing antimalarials and other drugs for COVID-19

PatchDock and SymmDock: Servers for rigid and symmetric docking

Anti-COVID-19 terpenoid from marine sources: A docking, admet and molecular dynamics study

Glecaprevir and Maraviroc are high-affinity inhibitors of SARS-CoV-2 main protease: possible implication in COVID-19 therapy

A Review of Coronavirus Disease-2019 (COVID-19)

Computational investigation of curcumin, a natural polyphenol that inhibits the destabilization and the aggregation of human SOD1 mutant (Ala4Val)

SARS-CoV-2 Main Protease: A Molecular Dynamics Study

Screening of phytochemicals as potent inhibitor of 3-chymotrypsin and papain-like proteases of SARS-CoV2: an in silico approach to combat COVID-19

In silico Screening of Natural Compounds as Potential Inhibitors of SARS-CoV-2 Main Protease and Spike RBD: Targets for COVID-19

Rapid Identification of Potential Inhibitors of SARS-CoV-2 Main Protease by Deep Docking

SARS-CoV-2/COVID-19: Viral genomics, epidemiology, vaccines, and therapeutic interventions

The SARS-CoV-2 main protease as drug target

Genome Composition and Divergence of the Novel Coronavirus (2019-nCoV) Originating in China

The outbreak of COVID-19: An overview

Systematic Comparison of Two Animal-to-Human Transmitted Human Coronaviruses: SARS-CoV-2 and SARS-CoV

Design of wide-spectrum inhibitors targeting coronavirus main proteases

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

Virus-encoded proteinases and proteolytic processing in the Nidovirales

Pr e ss In Pr e

The authors declare that they have no conflicts of interest.

The authors did not receive any external funding.