The emergence of COVID-19 has put the whole world in a standstill for more than 60 days. The disease originated from Wuhan (Hubei, China) in December 2019 and spread to the rest of the world, and was then declared a global pandemic in February 2020. COVID-19 is caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), a member of the coronavirus family. They are the largest known RNA viruses made up of positive single-stranded RNA and are subdivided into alpha, beta, gamma and delta coronaviruses (Shereen et al., 2020; Huynh et al., 2012) . This set of viruses have been documented to have a zoonotic origin from bats, mice or domestic animals and are linked with severe respiratory illness in humans Moore and June, 2020) . The SARS-CoV was first implicated as the causative agent of the severe acute respiratory syndrome outbreaks between 2002/2003 in the Guangdong area of China (Zhong et al., 2003; Zhong, 2004) , and the Middle East respiratory coronavirus (MERS-CoV) in 2012 (de Groot et al., 2013) . SARS-CoV-2, another variant of this devastating virus family has forced the human race to be on lockdown as a preventive measure to avoid disease transmission (contraction and spreading of the disease). The disease symptoms include but are not limited to fever, dry cough, sore throat and difficulty breathing. As of 29 May 2020, there are 5 701337 confirmed cases and 357 688 mortalities due to this alarming infection have been reported from 216 countries (WHO, 2020) .
Unfortunately, there is neither a specific drug nor vaccine approved for the treatment or modulation of immune response against this infection, yet the infection and mortality rates are steadily increasing. Angiotensin-converting enzyme-2 (ACE2) receptor in the host was identified as one of the receptors that the spike protein of SARS-CoV-2 interact with during infection, and was explored as a potential target for treatment of COVID-19. Chloroquine and its derivative Hydroxychloroquine were the first drugs to be repurposed for use in the treatment of COVID-19 (Omar et al., 2020) . Although these drugs were potent in inhibiting SARS-CoV, and showed no serious health effects when used for treatment of malaria and autoimmune diseases (Vincent et al., 2005) , they are now red flagged for use in COVID-19 (Mehra et al., 2020) . The efficacy of these drugs (Colson et al., 2020) together with Lopinavir/Ritonavir (Lim et al., 2020) , and Remdesivir are currently undergoing investigation and clinical trials. However, the World Health Organization (WHO) have suspended the use of Chloroquine and Hydroxychloroquine after they were reported to have serious side effects including cardiac toxicity (Mehra et al., 2020) . As the world is gradually easing their lockdown restrictions, there are concerns about the possibility of the infections increasing exponentially, hence there is an urgent need to identify and develop drug candidates and vaccines for the treatment and containment of SARS-CoV-2 infections. Notably, traditional drug discovery and development pipeline are time-consuming, very costly and often associated with high clinical failure (Aruleba, 2018) . Therefore, employing computational drug discovery tools for this raging infectious agent is of paramount importance. In silico techniques can augment the opportunities for designing, verifying activity, and developing new drug and vaccine candidates within a short space of time. In this regard, bioinformatic tools have been utilized extensively in many interesting studies; and recently for the SARS-CoV-2 (Elmezayen et al., 2020; Gupta et al., 2020; Bhattacharya et al., 2020) .
The virus main protease (M pro ) is one of the best characterized drug targets because its inhibition blocks viral replication . Holistically, the M pro is highly promising as no human protease has a similar cleavage pathway, thus making it a non-probable toxic agent . The interaction of this enzyme and several compounds as inhibitors have shown promising results (Khaerunnisa et al., 2020; Muralidharan et al., 2020; Ton et al., 2020) . Furthermore, M pro is highly conserved among the coronavirus family showing 40-44% of sequence homology (Muralidharan et al., 2020) . Inhibiting the virus replication and maturation in the host cells by blockade of the enzymatic activity of the SARS-CoV-2 M pro could potentially give the patients a fighting chance. Hence, we have employed in silico molecular docking to explore whether Lamivudine, Pradimicin A, Plerixafor, and Lopinavir could inhibit the M pro of the SARS-CoV-2 and used molecular dynamics simulation to validate the stable and strong binding interaction between the protein and the ligands with the aim of providing structural basis for rational drug discovery.
The SARS-CoV-2 M pro receptor with Protein Database (PDB) ID: 6W63 and the potential ligands: Pradimicin A (PubChem ID: 5479145), Lamivudine (PubChem ID: 60825), Plerixafor (PubChem ID: 65015), and Lopinavir (PubChem ID: 92727) were used in this study. The ligands were downloaded from NCBI PubChem database in 3D conformer (SDF) format at (https://www.ncbi.nlm.nih.gov/pccompound) while the crystallized 3D structure of the receptor in complex with co-crystallized ligand: X77, PDB ID:6W66 was retrieved from protein data bank at http://www.rcsb.org/pdb/home/home.do. The M pro selection is as a result of the presence of an inhibitor ligand (X77) at the active site. The antiviral drug, Lopinavir served as control in this study.
To determine whether the ligands can easily gain access to the target site in M pro after entering the blood stream, the QikProp module in the Schrodinger-2019-4 software package was used to evaluate the 3D molecular structures of the ligands' pharmacokinetic properties. With respect to Lipinski's rule of five, the absorption, distribution, metabolism, excretion and toxicity (ADME/T) properties were predicted (Ligprep and Macromodel, 2011).
These compounds were further subjected to ligand preparation (LigPrep) prior to docking. LigPrep is a robust collection of tools designed to prepare high quality, all-atom 3D structures for large numbers of drug-like molecules, starting with 2D or 3D structures in spatial data file (SDF) or Maestro format (Release, 2017) . In order to generate a single, lowenergy, 3D structure with correct chiralities for each successfully processed input structure, the ligands were subjected to preparation by LigPrep tool. This tool is also capable of generating multiple structures from input structure with various ionization states, tautomers, stereochemistries, and ring conformations, and eliminate molecules using specified criteria. The ligands for this study were also prepared by optimizing geometries through OPLS_2005 force field and Ionization of possible state were generated at pH 7.0 ± 2.0 by Epik (Greenwood et al., 2010) . Desalt and generate tautomers were also selected and the stereoisomer computation was checked to retain specific chiralities (vary other chiral centers) and to generate at most 32 conformations per ligand.
Excellent docking studies require both accurate software and starting structures . Experimentallyderived structures can waste time and resources or produce false result if not corrected. Schr€ odinger's Protein Preparation Wizard is designed to ensure structural correctness at the onset of a project, equipping them with a high-confidence structure ideal for use with a wide variety of modeling applications (Sastry et al., 2013) . Because the typical structure file of the receptor obtained from the PDB is not suitable for immediate use in molecular modeling calculations, the protein preparation wizard was used to prepare the receptor (PDB ID: 6W63). Following the previously reported method , the Protein Preparation Wizard was used for the following: remove alternate conformation, remove HetAtoms from the protein structure, add hydrogen atoms, correct missing or incorrectly specified residues, detect and correct valence violations, optimize protonation state and hydrogen positions.
An energy-optimized pharmacophore (e-pharmacophore) model combines the advantages of structure-based and ligand-based drug-design theories, and can be used to rapidly screen ligands based on pharmacophore properties (Salam et al., 2009; Negi et al., 2014) . Using the 'Develop Pharmacophore hypothesis' option in Maestro v12.2 tasks view mode, pharmacophore sites were generated from the receptor with redocked co-crystalized ligand (M pro -X77) complex, preserving a maximum of seven pharmacophore features as default. Pharmacophore chemical properties include: hydrogen-bond acceptor (A), represented as vectors, hydrogen-bond donor (D) as projected points, aromatic ring (R) as ring, positive ionizable (P), and negative ionizable (N) (Veeramachaneni et al., 2015) . Explicit matching was required in the e-pharmacophore approach for generation of the most energetically favorable site. The hypothesis settings were configured to treat atoms as projected points with a radii scaling factor of 0.50 and limit excluded volume shell to 5.0 Å.
Glide uses a hierarchical series of filters to search for possible locations of the ligand in the active-site region of the receptor. This is done by ligand ranking via high-throughput virtual screening (HTVS). Three modes of sampling ligand conformational and positional degrees of freedom are available to determine the optimal ligand orientation relative to a rigid protein receptor geometry (Repasky et al., 2007) . Flexible ligand docking with Glide, optionally includes ligand constraints or ligand molecular similarities. Docking grid of the receptor's active site was detected using the PDB file of the coordinates with the receptor grid generation tool in maestro v12.2. This site defines the area around the active site in term of co-ordinates x, y and z. The receptor grid box resolution was centered at coordinates À20.57, 18.10, and À26.99 corresponding to x, y and z-axis, respectively. Docking and calculations were run in the extra precision (XP) mode of Glide (Schr€ odinger, 2018).
The free binding energies of ligand docked complexes were computed using the molecular mechanic-generalized Born surface area (MM-GBSA) (Huang et al., 2006) with specific parameters. Based on the docking score and MM/GBSA binding-free energy, Jin et al. (Jin et al., 2011) developed correlation model between docking scores or calculated binding-free energies and experimental pIC50 values. The Prime module in Maestro was employed to calculate the MM-GBSA energy of Glide XP docked complex. The OPLS_2005 force field in conjunction with GBSA continuum model (Yu et al., 2006 ) was used to calculate energies of complexes of ligands following the equation reported by Lyne et al. (Lyne et al., 2006) and according to our previous study .
The docked complexes were subjected to MD simulation using Desmond module of Schr€ odinger software with OPLS 2005 force field. The protein-ligand complex was bounded with a predefined TIP3P water model (Jorgensen et al., 1983) in orthorhombic box. The volume of the box was minimized and the overall charge of the system was neutralized by adding Na þ and Clions. The temperature and pressure were kept constant at 310 Kelvin and 1.01325 bar using Nose-Hoover thermostat (Hoover, 1985) and Martyna-Tobias-Klein barostat (Zhong, 2004) methods. The simulations were performed using NPT ensemble by considering number of atoms, pressure and timescale. During simulations, the long-range electrostatic interactions were calculated using Particle-Mesh-Ewald method (Essmann et al., 1995) . RMSD plots for the backbone atoms for both the protein and ligand bound protein were generated to decipher the relative stability of the ligand in the binding pocket of the M pro . The results were analyzed and visualized by simulation interaction diagram and MS-MD trajectory analysis.
The chemical 2D structures of the four ligands ( Figure 1 ) were obtained from PubChem database. In essence, the database is a public repository composed of numerous validated chemical entities and their biological activities (Kim et al., 2016) .
Based on the Lipinski's rules of five, the pharmacokinetic properties (ADME) of Pradimicin A, Lamivudine, Plerixafor, and Lopinavir were calculated using the QikProp module in the Schrodinger-2019-4-software package. Here, Lipinski stated that a drug/compound would be orally bioavailable if it follows the following criteria such as molecular weight < 500 amu, Hydrogen bond acceptor sites < 10, Hydrogen bond donor sites < 5, and Lipophilicity value LogP 5. The result of the present study showed that only lamivudine satisfied Lipinski's rules, which indicated good oral bioavailability (Table 1) .
Extra precision docking was used to investigate favorable interactions between the selected ligands and the SARS-CoV-2 M pro . The redocked X77 replicated the binding pose and orientation of the crystal ligand, and indicated that it was successful in mimicking the native pose ( Figure 2 ). Using the same docking approach, the M pro receptor was docked with the four ligands ( Figure 1 ). Pradimicin A, Lamivudine, and Plerixafor were considered for docking studies. An FDAapproved anti-viral drug Lopinavir was also docked to correlate the docking scores and interactions (Figure 3) .
The results revealed that Pradimicin A was properly positioned into the binding cavity constructed by polar amino acids: GLN192 THR190 GLN189 ASN142 THR45 SER46 THR26 THR25 THR24; hydrophobic amino acids: ALA191 PRO168 LEU167 MET165 CYS44 MET49; and charged amino acids: GLU166 with binding energy of À2.4 kcal/mol (Figure 3(A) ). The binding parameters of the ligands together with the number of hydrogen bonds and the residues within the distance of 4 A o are presented in Table 2 . From the 2D interaction result, CYS44 exhibited bi-hedral H-bonding with hydroxyl group of oxane ring (tetrahydropyran) of Pradimicin A and other oxygen group showed H-bonding with GLN189 (Figure 3(A) ).
Lamivudine was properly positioned into the binding cavity constructed by polar amino acids: ASN 142 SER144, hydrophobic amino acids: LEU141 PHE140 CYS145 MET165 LEU167, charged amino acids: GLU166; with binding energy À5.3 kcal/mol. The hydroxyl and Oxo groups of lamivudine exhibited bi-hedral H-bonding with GLU166. The amine group and keto groups of the ring structure were also involved in hydrogen bonding with PHE140 (Figures 3(B) and Ye et al., 2020and Table 2) .
Lopinavir was properly docked into the binding cavity of the receptor as constructed by polar amino acids: GLN192 THR190 GLN189 ASN142 SER144 THR25 THR45 SER46; hydrophobic amino acids: PRO168 LEU167 MET165 PHE140 LEU141 CYS145 VAL42 CYS44 TYR54 PRO52 MET49; Charged amino acids: ASP187 GLU166; and glycine: GLY143 with binding energy of À9.2 kcal/mol (Figure 3) . Among the residues of Lopinavir and M pro , three H-bonds where involved between the residues GLU166 and GLN189 and the interacting atoms of Lopinavir (Figure 4 and Table 2 ).
Plerixafor interacted with the binding cavity constructed by polar amino acids: ASN142 THR25 THR169; hydrophobic amino acids: CYS145 LEU27 MET49 CYS44 VAL42 MET165 LEU167 PRO168; charged amino acids: GLU166 and glycine amino acids: GLY170 with a binding energy of À3.7 kcal/mol. Two amine groups showed bilateral H-bonding with GLU166 and other bonds with CYS145 and ASN142 (Figure 3(D) and Figure 4 ).
To predict the binding mode and binding free energy (DG bind ), the Prime MM-GBSA simulation was calculated for M pro -ligands and M pro -co-crystalized ligand complexes utilizing Maestro v12.2 ( Table 2 ). The DGbind (kcal/mol) of all the interactions were reported in order to understand the binding affinities between the ligands and the receptor to provide an insight into the stability of the ligands in the active site of M pro using the molecular mechanics-generalized Born surface area (MMGBSA) method. We computationally estimated the binding energy of the co-crystalized ligand to be À48.56 kcal/mol. The binding energy of lopinavir (À85.83 kcal/mol) to the Mpro appears to be the lowest among the ligands estimated followed by plerixafor (À75.78 kcal/mol). Pradimicin and lamivudine showed relatively higher (À45.89 and À34.94 kcal/mol, respectively) when compared to the co-crystalized ligand. The binding free energy determination, based on Prime MM-GBSA, established the stability of M pro -ligands complexes.
The simulation interaction diagram of Schrodinger suit was used to analyze and survey the protein-ligand interactions for the desmond molecular dynamic trajectory carried out on M pro and the studied ligands' complexes. (3) 2 GLU166, GLN189 and molecular dynamics (MD) simulation are a widely used method for the exploration of novel inhibitors against a target protein (Jani and Dalafave, 2012) . MD simulation provides information regarding the time dependent behavior of any molecular system by integrating Newton's laws of motion (Mandlik and Singh, 2016) . The stability of the docked complexes were determined by performing an MD simulation as previously described (Al-Shabib et al., 2018; AlAjmi et al., 2018) . The protein-ligand root mean square deviation (PL-RMSD) was used to measure the scaler distance between M pro and the ligands throughout the trajectory (100 nsecs) and also to evaluate the evolvement of the Ca protein backbone as well as the protein and ligand RMSD over the course of 100 nsecs period. All the protein frames were initially aligned on the reference frame backbone and then the RMSDs were calculated based the Ca, the protein and the ligands. Monitoring the RMSD of M pro provided insights into the structural conformation throughout the 100 nsecs. Furthermore, it was used to evaluate the equilibration of the complexes and its fluctuation directly to simulation within the sum of thermal average.
shows the PL-RMSD of the co-crystalized ligand (X77) and studied ligands, respectively, in the active site of M pro over the 100 nsecs simulation. This illustrates the super positioning or fitting as well as measuring the RMSD using the protein and ligand's heavy atoms. The ligand RMSD shows the stability of X77 as well as the studied ligands with respect to the protein binding pocket that is, the internal fluctuation within the ligands' structures is low. MD simulation of the docked complexes protein backbone RMSD plots indicate the stability of the ligands in the active site of the M pro . The ligands' backbone RMSD plots indicate that all the four ligands are comparable to the co-crystalized ligand and maintained their interactions with M pro .
The characterization of the local fluctuation of the protein was further investigated using the protein root mean square fluctuation (P-RMSF) as depicted in Figures 5(B) and 7. The peaks indicate the areas of the protein that fluctuate the most during the simulation period with respect to residue location (the X-axis). Additionally, the ligands' contacts were shown in bars. This corresponds to the protein residues interacting with the ligands. 31 residues of the M pro were involved in the contacts of pradimicin A and plerixafor. With respect to 14 contacts with the X77, lamivudine and lopinavir interacted with 19 and 21 Mpro residues, respectively. Figure 5(C) shows the Ligand RMSF which explains how the ligands fragment interact with M pro and the entropic role in their binding events. Most of the atoms of the X77 (Figure 5(C) ), plerixafor and pradimicin A (not shown) showed a high degree of fluctuation due to the fact that the atoms are solvent exposed and are able to rotate easily around the atom bonds and interact with water. Atoms of lamivudine and lopinavir (not shown) are buried deeply in the binding pocket of M pro with less chance to rotate freely and thus, have lower fluctuation.
The timeline of the contacts of the residues of the M pro to X77 over the 100 nsecs period as well as the time of contact (X-axis), the specific residues (Y-axis) and number of contacts (color intensity) were shown in Figure 5 (D). The M pro -ligands interaction was monitored throughout the simulation as presented in Figures 5(E) and 8 for the co-crystalized ligand as well as the studied ligands, respectively. These interactions or the amount of times the ligands made contact with M pro were categorized into four types: Hydrogen Bonds, Hydrophobic, Ionic and Water Bridges (Ojo et al., 2020) . The result of the M pro -ligands complexes were compared with Mpro-co-crystalized ligand and observed that the contacts made with M pro were three namely; Hydrogen Bonds, Hydrophobic, and Water Bridges. GLU166 was maintained almost throughout the period of simulation with hydrogen bond and water bridges. In comparison with the four ligands used in this study, several amino acids were constantly in contact with M pro but the residues ASN142, GLY143, GLN189, GLN192, and GLN166 were significantly in contact with M pro . Residues SER46 and ASN142 were significantly common to Lopinavir, Pradimicin A, and Plerixafor maintained 90%contacts M pro over the period of the simulation time.
Computational biology entails in silico investigation which is an essential arm of biotechnology focused at enhancing a deeper insight of biomolecular interactions in order to address cellular disease pathogenesis whilst having immense contribution towards design and development of possible therapeutic candidates (Abdullahi et al., 2017) . Indeed, this technique have assisted in identifying lead compounds for various diseases Shanmuga Priya et al., 2018; . On the other hand, the emergence and sporadic spread of COVID-19 has called for the urgent identification of novel drug and/or repurposing of existing available chemotherapeutic agents against this pandemic. The virus polyprotein encodes two proteases, the 3-C-like protease (main protease) and a papain-like protease, both are vital targets for drug discovery platforms against coronaviruses (Kandeel and Al-Nazawi, 2020) . Accordingly, the M pro is the most probable antiviral candidate due to its crucial role in self-maturation and ensuing development of polyproteins (Muralidharan et al., 2020) . Hence, we explored the molecular implications of blocking the activity of M pro using compounds (Pradimicin A, Lamivudine, Plerixafor and Lopinavir) that have been shown to possess pharmacological activity against viruses.
In view of the aforementioned, we redocked the 6W63 with its co-crystalized ligand using the Glide XP docking protocol.
With the receptor grid and extra precision option in the precision mode, X77 was docked into the active site of the SARS-CoV-2 M pro receptor. The M pro -X77 redocked complex was used as an input for the e-pharmacophore to generate a hypothesis. The initial number for pharmacophore site generation was set to default; however, five pharmacophore sites were predicted, and the final hypothesis consisted of three aromatic rings (R9, R10, and R11) and two H-bond acceptors (A3 and A4) as shown in Figure 2 (E). The quantitative structureactivity relationship (QSAR) of this enzyme (M pro ) also known as novel 3C-like protease (3CLpro) was developed using a multi-linear regression (MLR) based 2D-QSAR model for the identification of important structural features responsible for the inhibition of Mpro in SARS-CoV-2 infections using eight significant descriptors classified as topological, 2D atom pair, functional group count, and atom entered fragment descriptors (Kumar and Roy, 2020) . The study reported that molecules or ligands consisting of features such as O-S fragments and thiophene rings could contribute to higher M pro enzyme inhibitory activity, and the presence of pyridine rings, methylene groups connected with the highly electronegative atoms, imides (including thioimides) group, and N-O fragment may lower the inhibitory activity for 3CLpro enzyme. Based on this context, the features that negatively influence the inhibitory activity against M pro are not contained in the studied ligands' structures and thus the ligands may have high M pro enzyme inhibitory activities.
Pharmacophore-based screening was performed against the ligands in order to determine M pro inhibitors with desired chemical features and to compensate for the QSAR due to limited dataset in this study. The generated hypothesis can be used in a virtual-based screening study to generate additional inhibitors of M pro if they satisfy a minimum match of four sites on the generated e-pharmacophore hypothesis.
Subsequently, molecular docking using the Glide molecular tool revealed that Lopinavir bound to the active site with the highest affinity and lowest free energy with a Glide score of À9.2 kcal/mol; suggesting that it could be a vital inhibitor of protease activity of coronavirus . Following Lopinavir in terms of excellent binding affinity for M pro is Lamivudine; having a glide score of À5.3 kcal/mol. Earlier studies have shown this compound to be a potent anti-viral agent and it has been proposed for COVID-19 treatment due to its ability to impair DNA synthesis (Lai et al., 2020) . Jarvis et al. (Jarvis and Faulds, 1999) highlighted Lamivudine as an excellent inhibitor of hepatitis B virus (HBV) replication and suppresses viral replication in HIV patients. Interestingly, Lamivudine fulfilled the five Lipinski's rule in our study. In light of the anti-viral activity of Plerixafor, we repurposed it for the management of COVID-19. Plerixafor had a glide score of À3.7 kcal/mol against the SARS-CoV-2 M pro , a score lower than the co-crystallized ligand. Nevertheless, this compound is an antagonist to C-X-C chemokine receptor type-4 (CXCR4) (Keating, 2011) , with anti-HIV activities. It inhibits the attachment of the virus to the host cell. Amongst all the docked M pro -ligand complexes, the M pro -Pradimicin A complex had the lowest binding affinity of À2.4 kcal/mol. However, it has been shown to possess inhibitory activities against virus infection/replication. Importantly, the anti-viral mechanism of action of Pradimicin A is through binding to mannose residue (Tanabe-Tochikura et al., 1990) and recent studies have documented anti-HIV and anti-coronavirus ( Van der Meer et al., 2007) properties for this compound.
Finally, docking of the ligands with the Mpro revealed their binding modes. To account for the flexibility of Mpro and ligands and to evaluate their binding affinity with Mpro, the MD simulation of the docked complexes were carried out at three different time points (1.2 nsecs, 10 nsecs and 100 nsecs) in order to avoid false positive result of the stability result (only the 100 nsecs was reported in this study). Binding mode analysis revealed that the binding modes obtained after MD simulation were more or less similar to that obtained after docking. The presence of a large number of H bond acceptors, H bond donors as well as hydrophobic groups in the ligands account for the stability of the ligand inside the binding pocket of Mpro. Based on the RMSD of the ligand-protein complexes, it was confirmed that ligands maintained their interaction with Mpro with lower root mean square fluctuations.
Globally, various research institutions and hospitals are on a quest to identify suitable drug targets to stop/manage SARS-CoV-2. However, using the traditional methods for identification of drug targets for COVID-19 might not be viable as the disease is already crippling the health care resources. Hence, anti-viral drugs are repurposed for immediate attention. This have proven to be the most expensive exercise as the drugs did not behave as hypothesized. As such, Bioinformatics could serve as a tool to screen potential drug activity before use. Accordingly, computational drug design and repurposing has laid a foundation for researchers to investigate various chemotherapeutic agents towards identifying potent drug candidates for clinical trials.
Our study revealed that all the antiviral drugs had satisfactory binding affinity for the SARS-CoV-2 M pro . Amongst all the inhibitors, Lopinavir had the best binding score while pradimicin A had the lowest binding score for the SARS-CoV-2 M pro . These results were comparable after we used molecular dynamics simulation to validate the stability of all the ligand-receptor complexes. The ligands bound steadily to M pro active site and restricted the target movement. GLU166 was found to be an important residue in the active site of M pro necessary for binding. Furthermore, Lamivudine warrants an in-depth investigation taking into account its binding affinity for M pro and fulfilment of the Lipinski's rule of five. It is worth noting that this is the first report that indicated the inhibitory effect of Pradimicin A and Plerixafor against the M pro , and could serve as potential drug candidates for COVID-19 treatment. Overall, these data can be utilized in drug repurposing pipelines in the continued search for SARS-CoV-2 therapy with high efficacy. The QSAR of these ligands alongside other dataset is recommended for further analysis in order to provide the structure features against the inhibitory activities of the SARS-CoV-2 main protease (M pro ).