Remember me
ORIGINAL ARTICLE
Year : 2021 | Volume
: 20
| Issue : 4 | Page : 371-392
Repurposing of available antiviral drugs against SARS-CoV-2 by targeting crucial replication machinery proteins: a molecular docking study
Ahmed F El-Sayed1, Amira T Mohammed2, Wafy Hamed3, Shaymaa Abdelmalek4
1 Department of Microbial Genetic, Genetic Engineering and Biotechnology Division, National Research Centre, Giza, Egypt
2 Department of Biochemistry, Genetic Engineering and Biotechnology Division, National Research Centre, Giza, Egypt
3 Department of Microbiology and Immunology, Faculty of Veterinary Medicine, Sadat University, Menoufia, Egypt
4 Department of Microbiology and Immunology, Faculty of Veterinary Medicine, Cairo University, Giza, Egypt
Correspondence Address:
PhD Ahmed F El-Sayed
Researcher at Department of Microbial Genetics, Genetic Engineering and Biotechnology Division, National Research Centre, Giza
Egypt
Source of Support: None, Conflict of Interest: None
CheckDOI: 10.4103/epj.epj_70_21
Background and objectives The newly emerged severe acute respiratory syndrome coronavirus is spreading worldwide rapidly with increasing incidence rates. Due to lack of effective treatments and vaccines, various drug repurposing studies are being developed. Searching for available antiviral drug libraries is the best and fast option to advance to clinical trials and spread their application among infected patients.
Materials and methods Molecular docking study was performed utilizing AutoDock 4.2 system and Discovery Studio 4.5, which were utilized to predict the activity pocket of the target proteins.
Results and conclusion The results found that the interacting affinities resulted from the molecular simulation of 3CL protease with ligands Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir were −7.2, −7.4, −7.2, −6.3, −6.1 and −6.6 kcal/mol, respectively. Similarly, the interacting energies obtained from the docking of RNA helicase with ligands were −7.9, −7.4, −6.4, −7.9, −6.2, and −6.9 kcal/mol. Also, the binding energies obtained from the docking of 3′-5′ exoribonuclease with ligands were −10.6, −10.1, −6.5, −7.1, −6.1, and −9.3 kcal/mol. Finally, the binding energies score from the docking of the RNA-dependent RNA polymerase with ligands was −9.6, −6.9, −6.2, −6.6, −6.7, and −6.4 kcal/mol. Based on the binding energy score and docking result, Ledipasvir and Sofosbuvir have a higher affinity of the drug molecule such as against RNA-dependent RNA polymerase, exonuclease, and 3CL protease. Besides, Ledipasvir and Galidesivir show prominent binding interaction with severe acute respiratory syndrome coronavirus RNA helicase. The results are promising for evaluated drugs especially Ledipasvir and Sofosbuvir and could be useful in emergency treatment of coronavirus disease 2019 patients.
Keywords: 3CL protease, coronavirus disease 2019, exonuclease, helicase, molecular docking, RNA-dependent RNA polymerase
A novel human strain of coronavirus severe acute respiratory syndrome coronavirus (SARS-CoV-2) has emerged in Wuhan, China [1]. The WHO has proclaimed the outbreak of coronavirus illness [coronavirus disease 2019 (COVID-19)] a worldwide crisis and a pandemic [2]. COVID-19 is caused by SARS-CoV-2 and can cause symptoms 2–14 days after exposure, with fever, cough, difficulty breathing, general fatigue, loss of taste and/or smell, sore throat, nausea, or vomiting and diarrhea [3],[4].
Human coronaviruses, positive-sense single-stranded RNA viruses, have a genome size of about 30 000 bp, which encode the replicate complex (ORF1ab), translated in the form of polyproteins (pp), which contain four structural proteins and 20 nonstructured proteins (NSP). Some of these are NSPs, RNA-dependent RNA polymerase (RdRp) (NSP12), 3CL protease (3CLpro), 3′-5′ exoribonuclease, and RNA helicase. The last two enzymes play an important role in replication machinery, replication fidelity, and RNA capping and replication of the longest RNA genomic material of corona viruses [5],[6].
RdRp is a vital protein in the replication of RNA viruses, including coronaviruses. It focuses on several types of RNA viruses, including hepatitis C virus (HCV), Zika virus, and coronaviruses [7]. Among the coronaviruses NSPs, there are two proteases that are vital to virus replication [8]. As part of translational NSPs processing and release, the major protease is known as the 3-C-like protease (3CLpro) and the papain-like protease (PLpro), both of which are responsible for the recent coronavirus epidemics of SARS and Middle East respiratory syndrome (MERS) [9]. The main protease (3 CLpro) is a key enzyme in the processing of polyproteins pp1a and pp1ab. ORF1a and ORF1ab are cleaved by papain-like protease (PLpro, nsp3) and 3C-like protease (3CLpro, nsp5) [10]. The SARS-CoV 3CLpro has an important function and considered an active target for antiviral drugs [11]. 3CLpro is the major protease responsible for the control of many major viral functions and has a highly conserved SARS-CoV catalytic domain [12]. Some of its functions include viral replication processes, making it an ideal drug development goal [13].
The SARS epidemic and recent emergence of MERS feature the possible lethality of zoonotic coronavirus diseases in people. No particular antiviral treatment alternatives are accessible [5].
In this investigation, the virtual screening concentrate against the first known SARS-CoV-2 was performed. The obtained results will help in the repurposing of already available marketed Food and Drug Administration (FDA)-approved antiviral drugs, Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir to combat the recent dangerous SARS-CoV-2 nonstructured proteins (RdRp, 3CL protease, exonuclease, and RNA helicase). The antiviral drugs that the FDA has approved, act on nonstructured proteins (nucleotide analogs and polymerase inhibitors). This study helps the scientists to do clinical trials in the field, which gives the hope to the world to fight the newly emerged threatening outbreak.
Materials and methodsPreparation of protein receptors
The two three-dimensional (3D) structures of the coronavirus proteins with their PDB:ID 3CLpro (5r7y) and RdRp (6M71) used in our study were downloaded from the Protein Databank. The other two proteins (helicase and exoribonuclease with accession numbers YP_009725308.1 and YP_009725309, respectively) and their amino acid sequences were downloaded from NCBI and modeled to form a 3D structure by Modeler 20.v. All the 3D structures of the objective proteins were performed by erasing the local native drug and water molecules. In addition, polar hydrogen were added according to the amino acid protonation state at pH 7.0 utilizing AutoDock, version 4.2 program. Thereafter, nonpolar hydrogens were merged while polar hydrogens were added to each protein. Repeat this process for each protein, then save it in the pdbqt format and insert it for molecular binding. We then added hydrogen atoms to the molecule using the MG tool in AutoDock Vina4 the Molecular Graphics Lab. at The Scripps Research Institute.
Preparation of ligands
All FDA-approved antiviral inhibitors (Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir) were extracted and collected from ZINC database. The preparation steps of compounds and inhibitors were performed by converting from Structure Data Format to MOL2 chemical format using an Open babel. In addition, nonpolar hydrogen was deleted. Also, polar hydrogen charges were added. Furthermore, the torsions degrees were set to 0. Also using open babel, it was converted to the pdb format. Finally, inhibitor ligands were converted and saved into the dockable pdbqt format through AutoDock tools. The structures of prepared ligands were energy minimized and used as input files for AutoDock 4.2.5.
Docking studies
Docking studies were endeavored to investigate the binding method of the recommended inhibitors onto the 3D model of the protein of SARS-CoV-2 utilizing AutoDock tools 4.2. Before docking, polar-H atoms were added to the SARS-CoV-2 model followed by Gasteiger charge computation utilizing AutoDock devices. The macromolecule files were energy minimized and then saved in the pdbqt format and fit to be utilized for docking. Polar-H charges of the Gasteiger-type were allocated, and nonpolar-H were converged with the carbons. Default settings were utilized for all other parameters. PyMOL package was utilized to visualize the binding cooperation between these ligands with a 3D model of SARS-CoV-2 receptors. All the docking tests were performed utilizing AutoDock Vina as it offers (a) more precision in anticipating ligand–protein cooperation contrasted with its past AutoDock 4.2;(b) more limited running time due to its numerous center processors; and (c) more exactness for ligand handling more than 20 rotatable bonds.
Preparation of a grid box
The receptors file was saved in the pdbqt format. Ligand-centered maps were generated by the AutoGrid program with a spacing of 0.388 A° and grid dimensions of 90 A°×90 A°×90 A°. Grid box center was set to the coordinates 0.076, 0.093, −0.012 in the x, y, z format, respectively. All the docking tests for helicase and exoribonuclease were finished utilizing the strategy for blind docking (utilizing the grid box sufficiently huge to cover the entire protein structure).
Visualization
The illustration and interpretation of the docking analysis was achieved at the 3D and 2D level of interactions for the COVID protein–drug complex structure by the Discovery Studio 4.5 program. Molecular interactions of amino acids located in the activity pocket and catalytic groves with drug include hydrophobic bonds, hydrophilic bonds, hydrogen bonds, and their bond lengths in the form of graphical representation.
ResultsDocking and molecular interaction studies of severe acute respiratory syndrome coronavirus 3CL protease
Our results indicated that inhibitor drugs prepared for docking were screened for Lipinski’s rule of five as presented in [Table 1]. All the six potential inhibitors viz. Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir got docked onto the 5R7Y_A 3C-like protease as shown in [Figure 1]. The good, recorded binding energy value was obtained for Sofosbuvir (−7.4 kcal/mol). The best pose with the least energy is considered promising because it is very stable for the compound.
Table 1 Lipinski’s rule screening for physiochemical properties of molecules
Figure 1 The six drug ligands docked in SARS-CoV-2 main protease (PDB ID 5r7y) with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Ledipasivir binding with 3CLpro, (b) Sofosbuvir binding with 3CLpro. SARS-CoV-2, severe acute respiratory syndrome coronavirus.On the basis of simulation analysis, data in [Table 2] and [Figure 1]a show that the docking of Ledipasvir with the SARS-CoV-2 main 3CL protease shows high-affinity interaction −7.2 kcal/mol. Ledipasvir formed one hydrogen bond with Gln189. Also, Ledipasvir formed nine nonhydrogen bond interactions within the activity pocket; five alkyl bonds with residues Leu50, Pro168, Cys44, and Met49, two halogen (fluorine) bonds formed with Pro168 and Glu166, and only one π-anion contact with Glu166.
Table 2 Molecular contact of drugs with residues of 3C-like proteinase (residues having similar contact are red colored and in bold)Data shown in [Figure 1]b illustrated that the interaction of Sofosbuvir with the SARS-CoV-2 main 3CL protease in complex shows the highest affinity interaction (7.4 kcal/mol). Also, we noticed that four H bonds were formed with His41, Cys44, Glu166, and Gln189. Besides, six non-H-bonds were formed with Asn142 (Carbon-H), His41 (Van der Waals bond), Pro168 (Alkyl), Pro168 (Alkyl), Met165 (Pi-sulfur), and His41 (Pi-pi T shaped). Results were obtained by docking of Ribavirin with the SARS-CoV-2 main 3CL protease [Figure 2]a show high-affinity interaction −7.2 kcal/mol. Ribavirin formed four H-bonds with His41, Cys44, Glu166, and Gln189 and six hydrophobic reactions with Asn142 (carbon-H), His41 (Van der Waals), Pro168 (Alkyl), Pro168 (Alkyl), Met165 (Pi-sulfur), and His41 (Pi-pi T shaped). Also, [Figure 2]b shows that the docking of Galidesivir with the SARS-CoV-2 main 3CL protease shows moderate-affinity interaction with the chain A of protein with an affinity of −6.3 kcal/mol. Galidesivir formed four H bonds with Gln189, Gln189, Cys44, and Arg188 and two hydrophobic reactions with His164 (Pi-sigma) and Met49 (carbon-H).
Figure 2 The six drug ligands docked in SARS-CoV-2 main protease (PDB ID 5r7y) with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional interaction amino acid residues involved in the interaction (with ligand as color sticks) and 2D interaction binding of ligands with an amino acid with a hydrogen bond (green dash line) (a) Ribavirin binding with 3CLpro, (b) Galidesivir binding with 3CLpro. SARS-CoV-2, severe acute respiratory syndrome coronavirus.The docking of Tenofovir with 3CL protease in the complex has the least affinity interaction −6.1 kcal/mol. Tenofovir created five H bonds with Glu166, Arg188, His164, and Leu14, one hydrophobic reaction with Asn142 (carbon-H) ([Figure 3]a). Docking of Remdesivir with the SARS-CoV-2 main 3CL protease shows moderate affinity interaction −6.6 kcal/mol. As shown in [Figure 3]b, it has three H bonds with His41, Glu166, and Gln189. Also, four hydrophobic reactions with His41 (Van der Waals), Asn142 (Alkyl bond), Met165 (π-sulfur), and Pro168 (Alkyl). Finally, our data indicated that comparison of the docking analysis of six inhibitors find that the most common residues Gln189, Glu166, and His164 in the activity and binding pocket increase the interactions of all six compounds through the formation of H bonds.
Figure 3 The six drug ligands docked in SARS-CoV-2 main protease (PDB ID 5r7y) with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Tenofovir binding with 3CLpro, (b) Remdesivir binding with 3CLpro. SARS-CoV-2, severe acute respiratory syndrome coronavirus.Docking and molecular interaction studies of severe acute respiratory syndrome coronavirus RNA helicase
Our data obtained from this simulation supported and enhanced the suitable binding of the promising drug inhibitors Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir against the predicted 3D model of SARS-CoV-2 RNA helicase enzyme with docking score (binding energy). Our results show the binding of amino acids that take part in the drug–protein interactions. The result of the binding energy from the docking analysis of Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir to the 3D model of helicase is represented in [Table 3].
Table 3 Molecular contact of drugs with residues of three-dimensional model of severe acute respiratory syndrome coronavirus helicase (residues having similar contact are red colored and in bold)The docking of Ledipasvir with the SARS-CoV-2 RNA helicase enzyme has the highest affinity interaction with a binding of −7.9 kcal/mol. No H bonds were formed by Ledipasvir. But Ledipasvir formed 18 nonhydrogen bond interactions within the activity pocket, Ala553 (Pi-Alkyl), Ala553 (halogen), Ala553 (unfavorable bump), Val558 (Pi-Alkyl), Val558 (Alkyl), Val558 (Alkyl), Val456 (Alkyl), Thr416 (Pi-sigma), Leu455 (Alkyl), Leu581 (Alkyl), Arg560 (unfavorable bump), His554 (pi-anion), Cys556 (carbon-H), Leu412 (Alkyl), Leu412 (Alkyl) Tyr457 (Alkyl), Lys584 (Alkyl), and Lys584 (Alkyl) as shown in [Table 3] and [Figure 4]a.
Figure 4 The six drug ligands docked in three-dimensional (3D) model of SARS-CoV-2 helicase with the best binding mode in the pocket of a protein (with ligand as color sticks), 3D interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Ledipasivir binding with SARS-CoV-2 helicase, (b) Sofosbuvir binding with SARS-CoV-2 helicase. SARS-CoV-2, severe acute respiratory syndrome coronavirus.Also, the docking of Sofosbuvir with the 3D model of SARS-CoV-2 helicase shows high affinity interaction −7.4 kcal/mol. Also, we noticed that four H bonds were formed with Asn179, Arg560, Arg560, and Asp534. Besides, 11 non-H-bonds were formed with Leu412 (Alkyl), Leu417 (Alkyl), Ala407 (Amide pi-stacked), Pro406 (Alkyl), Pro406 (Pi-Alkyl), Arg409 (Alkyl), Pro408 (Alkyl), Pro408 (Pi-Alkyl), Arg560 (pi-anion), Asp534 (pi-anion), and Phe422 (Alkyl) as presented in [Figure 4]b.
Results in [Figure 5]a show the docking of Ribavirin with the 3D model of SARS-CoV-2 helicase showing moderate-affinity interaction −6.4 kcal/mol. Ribavirin have three H bonds with Asn557, Ala407, and Leu417but no hydrophobic interaction. The docking of Galidesivir with the 3D model of helicase for SARS-CoV-2 shows the same affinity as that of Ledipasvir (−7.9 kcal/mol). Galidesivir formed four H bonds with Gln189, Gln189, Cys44, and Arg188 and two hydrophobic reactions with His164 (Pi-sigma) and Met49 (carbon-H) as presented in [Figure 5]b. Also, the docking of Tenofovir with the model of helicase shows the least-affinity interaction −6.2 kcal/mol. Tenofovir formed eight H bonds with Asn557, Asn559, Leu412, Leu417, and Arg409 and four hydrophobic reactions with Leu417 (Alkyl) Pro406 (Alkyl), Phe422 (Alkyl), and Thr410 (Van der Waals) ([Figure 6]a). Data in [Figure 6]b illustrates the docking of Remdesivir with the model of helicase has a slightly better affinity interaction −6.9 kcal/mol. It formed three H bonds with His41, Glu166, and Gln189 and four hydrophobic reactions with His41 (Van der Waals), Asn142 (Alkyl bond), Met165 (π-sulfur), and Pro168 (alkyl).
Figure 5 The six drug ligands docked in three-dimensional (3D) model of SARS-CoV-2 helicase with the best binding mode in the pocket of a protein (with ligand as color sticks), 3D interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (b) Ribavirin binding with SARS-CoV-2 helicase, (b) Galidesivir binding with SARS-CoV-2 helicase. SARS-CoV-2, severe acute respiratory syndrome coronavirus.
Figure 6 The six drug ligands docked in three-dimensional (3D) model of SARS-CoV-2 helicase with the best binding mode in the pocket of a protein (with ligand as color sticks), 3D interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding interaction of ligands with an amino acid with a hydrogen bond (green dash line). (a) Tenofovir binding with SARS-CoV-2 helicase, (b) Remdesivir binding with SARS-CoV-2 Helicase. SARS-CoV-2, severe acute respiratory syndrome coronavirus.Various residues included in the ligand–helicase interactions are shown as a stick with red color. From our simulation analysis, it is concluded that the residues Arg560, Leu412, and Leu417 in the pocket cavity increase the interaction of all the six drugs by the creation of H bonds.
Docking and molecular interaction studies of severe acute respiratory syndrome coronavirus 3′-5′ exoribonuclease
Our data of the binding affinity from the simulation study of Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir to exoribonuclease for SARS-CoV-2 is represented in [Table 4].
Table 4 Molecular contact of drugs with residues of severe acute respiratory syndrome coronavirus 3′-5′ exoribonuclease (residues having similar contact are red colored and in bold)Results obtained by [Table 4], and [Figure 7]a show the docking of Ledipasvir with the SARS-CoV-2 3′-5′ exoribonuclease in the complex has the greatest affinity interaction −10.6 kcal/mol. Ledipasvir created four H bonds with Asn386, Asn386, Gln313, and Arg289 and Ledipasvir formed 22 nonhydrogen bond interactions within the activity pocket, seven alkyl bonds with residues Pro429, Pro342, Pro342, Val290, Val287, Ile338, and Arg310, three halogen (fluorine) bonds formed with Ala553, Asn386, and Trp385, five pi-alkyl bonds with residues Val290, Ile338, Arg289, Arg310, and Cys387, and only one Van der Waals contact with Lys339. Data in [Figure 7]b shows the docking of Sofosbuvir with the SARS-CoV-2 3′-5′ exoribonuclease in the complex showing very high-affinity interaction −10.1 kcal/mol. Also, we noticed that four H bonds were formed with Asn388, Asn388, Tyr368, and Ala553. Besides, 11 non-H-bonds were formed with Val290 (Alkyl), Arg289 (Alkyl), Phe426 (Pi-Pi T shape), Asp352 (Halogen), Phe367 (Pi-Pi stacked), Tyr368, (unfavorable Bond), Pro429 (Alkyl), His427 (Alkyl), Val287 (Alkyl), Ala553 (Pi-Alkyl), and Asn386 (Van der Waals).
Figure 7 The six drug ligands docked in SARS-CoV-2 3′-5′ exoribonuclease with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional (3D) interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding interaction of ligands with an amino acid with a hydrogen bond (green dash line). (a) Ledipasivir binding with SARS-CoV-2 3′-5′ exoribonuclease, (b) Sofosbuvir binding with SARS-CoV-2 3′-5′ exoribonuclease. SARS-CoV-2, severe acute respiratory syndrome coronavirus.Results obtained by [Figure 8]a show that the docking of Ribavirin with the SARS-CoV-2 3′-5′ exoribonuclease in the complex shows low-affinity interaction −6.5 kcal/mol. Ribavirin formed seven H bonds with Gln354, Asp352, Asn388, Asn386, and Tyr368 and two hydrophobic reactions with Ala553 (unfavorable bond) and Gln354 (unfavorable bond). Also [Figure 8]b shows that the docking of Galidesivir with the SARS-CoV-2 3′-5′ exoribonuclease in the complex shows moderate-affinity interaction −7.1 kcal/mol. Galidesivir made six H bonds with Ala553, Asn386, and Leu366 and three hydrophobic reactions with Val290 (pi-alkyl), Asn388 (carbon–H-bond), and His427 (pi-anion). The docking of Tenofovir with the SARS-CoV-2 3′-5′ exoribonuclease in the complex shows the least affinity interaction −6.1 kcal/mol. Tenofovir formed five H bonds with Asp331, Gln313, Trp385, Gly333, and Ile332 and two hydrophobic reactions with Trp385 (unfavorable bond) and Lys336 (alkyl), ([Figure 9]a). Data in [Figure 9]b indicate that docking of Remdesivir with the SARS-CoV-2 3′-5′ exoribonuclease in the complex shows high-affinity interaction −9.3 kcal/mol. It formed six H bonds with Asn388, Leu366, Asp352, and Asn386 and10 hydrophobic reactions with Ala353 (pi-alkyl), Val389 (pi-alkyl), Cys387 (Pi-alkyl), Phe367 (Pi-Alkyl), Pro335 (pi-alkyl), Pro335 (alkyl), Val290 (pi-alkyl), Trp292 (alkyl), Trp385 (alkyl), and Asn388 (unfavorable bond).
Figure 8 The six drug ligands docked in SARS-CoV-2 3′-5′ exoribonuclease with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional (3D) interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Ribavirin binding with SARS-CoV-2 3′-5′ exoribonuclease, (b) Galidesivir binding with SARS-CoV-2 3′-5′ exoribonuclease. SARS-CoV-2, severe acute respiratory syndrome coronavirus.
Figure 9 The six drug ligands docked in SARS-CoV-2 3′-5′ exoribonuclease with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional (3D) interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Tenofovir binding with SARS-CoV-2 3′-5′ exoribonuclease, (b) Remdesivir binding with SARS-CoV-2 3′-5′ exoribonuclease. SARS-CoV-2, severe acute respiratory syndrome coronavirus.All residues included in the drug–exonuclease complex are illustrated as bold with red color. Our studies supported the fact that residues Asn386, Tyr368, and Gln313 enhanced the binding interaction in the catalytic site by the formation of bonds. The best results were obtained for exonuclease inhibitors Ledipasvir and Sofosbuvir, with free-binding energies of −10.6 and −10.1 kcal/mol, respectively, followed by the Remdesivir and Galidesivir with −9.3 and −7.1 kcal/mol.
Docking and molecular interaction studies of RNA-dependent RNA polymerase (NSP12)
Data obtained in this molecular analysis of Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir to RdRp of the novel SARS-CoV-2 is represented in [Table 5]. Out of the six drugs, Ledipasvir and Sofosbuvir were used as promising inhibitors and represented the highest interactions with catalytic pocket residues of RdRp.
Table 5 Molecular contact of drugs with amino acids of RNA-dependent RNA polymerase (residues having similar contact are red colored and in bold)The result in [Figure 10]a showed the docking of Ledipasvir with RdRp, which have the greatest affinity interaction −9.6 kcal/mol. Ledipasvir formed two H bonds with Asp623 and Arg553. Furthermore, Ledipasvir has nine nonhydrogen bond interactions within the activity pocket, with Ala685 (Pi-Alkyl), Ala688 (Pi-Alkyl), Val560 (Alkyl), Ser682 (halogen), Ser501 (carbon H-Bond), Asp760 (carbon H-Bond), Arg555 (Alkyl), Arg624 (Alkyl), Lys621 (Alkyl), and Asp623 (carbon H-bond).
Figure 10 The six drug ligands docked in SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional (3D) interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Ledipasivir binding with SARS-CoV-2 RdRp, (b) Sofosbuvir binding with SARS-CoV-2 RdRp. SARS-CoV-2, severe acute respiratory syndrome coronavirus.Also [Figure 10]b illustrated that the docking of Sofosbuvir with RdRp has high-affinity interaction −6.9 kcal/mol. Also, we noticed that eight H bonds were formed with Asn691, Asp623, Arg553, Lys621, Ser682, Thr556, and Cys622. Besides, three non-H-bonds were formed with Asp623 (pi-anion), Lys621 (pi-cation), and Thr556 (halogen). Results obtained by [Figure 11]a showed that the docking of Ribavirin with COVID-19 (RdRp) have the least affinity interaction −6.2 kcal/mol. Ribavirin formed nine H bonds with Arg624, Thr556, Ala554, Asp452, Thr556, and Tyr456, and also two hydrophobic reactions with Asp623, (Van der Waals) and Ser682 (unfavorable bond). Besides, [Figure 11]b shows that the docking of Galidesivir with the COVID-19 (RdRp) shows moderate-affinity interaction with an affinity of −6.6 kcal/mol. Galidesivir formed four H bonds with Asp760, Asp760, Asp623, and Thr556 and three hydrophobic reactions with Thr556 (carbon H-bond), Asp760 (unfavorable bond), and Asp760 (pi-anion).
Figure 11 The six drug ligands docked in SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional (3D) interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding interaction of ligands with an amino acid with a hydrogen bond (green dash line). (a) Ribavirin binding with SARS-CoV-2 RdRp, (b) Galidesivir binding with SARS-CoV-2 RdRp. RdRp, RNA-dependent RNA polymerase; SARS-CoV-2, severe acute respiratory syndrome coronavirus.The docking of Tenofovir with COVID-19 (RdRp) shows slightly good affinity interaction −6.7 kcal/mol. Tenofovir formed nine H bonds with Thr556, Thr556, Asp452, Ala554, Arg553, Arg624, Asp623, Thr680, and Ser681 and two hydrophobic reactions with Arg555 (Alkyl) and Asp623 (carbon H-Bond) as presented in [Figure 12]a). [Figure 12]b shows that the docking of Remdesivir with (RdRp) shows moderate-affinity interaction −6.4 kcal/mol. It formed eight H bonds with Arg624, Ser682, Lys621, Asp761, Asp760, and Tyr619 and two hydrophobic reactions with Asp760 (pi-anion) and Tyr619 (carbon H-bond). Different residues included in the compound–target interactions are in red bold illustration. Our studies supported that residues Asp623, Arg624, Thr556, and Asp760 enhanced the binding interaction in the catalytic site by the formation of bonds. The best results were obtained for RdRp inhibitors Ledipasvir and Sofosbuvir −9.6 and −6.9 kcal/mol, respectively, followed by Tenofovir and Galidesivir with −6.7 and −6.6 kcal/mol, respectively.
Figure 12 The six drug ligands docked in SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) with the best binding mode in the pocket of a protein (with ligand as color sticks), three-dimensional (3D) interaction amino acid residues involved in the interaction (with ligand as color sticks) and two-dimensional interaction binding of ligands with an amino acid with a hydrogen bond (green dash line). (a) Tenofovir binding with SARS-CoV-2 RdRp, (b) Remdesivir binding with SARS-CoV-2 RdRp. SARS-CoV-2, severe acute respiratory syndrome coronavirus. Discussions Alarming outbreak of COVID-19 was caused by SARS-CoV-2. First discovered in Wuhan, Hubei, China the mode of transmission is human to human through air droplets. The WHO reported this illness as a global public health emergency on February 15, 2020 [14].
FDA has not approved any antiviral medications for the treatment of COVID-19. All clinical treatments are guided toward symptomatic treatment. Several studies have developed promising drug candidates that can reduce COVID-19 symptoms by inhibiting certain components of SARS-CoV-2. For example, Remdesivir showed adverse effects against SARS-CoV-2 [15]. There is currently no safe drug candidate against the SARS and MERS viruses because the epidemic has ended and no clinical trials have been conducted [8].
Our study is a molecular simulation study against the first known SARS-CoV-2. The obtained results will help in the repurposing of the already available antiviral drugs Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir to combat the recent dangerous SARS-CoV-2 NSPs (3CLpro, RNA helicase, and 3′-5′ exoribonuclease, and RdRp). The antiviral drugs such as Ledipasvir (NS5A inhibitor of HCV), Sofosbuvir (NS5B RNA polymerase inhibitor of HCV), nucleotide analogs as Ribavirin, Galidesivir, Remdesivir, and Tenofovir (HIV-1 and HBV reverse transcriptase inhibitor). RdRp is an essential enzyme in the replication of RNA viruses; hence, it is a main target in treatment of viral infections, including HCV, ZIKV, and CoVs [16]. The active site of RdRp is highly similar, with two successive and surface-accessible aspartates in a beta-turn structure [17].
The protein 3-chimotrypsin like SARS-CoV (3CLpro) has been shown to be important for replication and has therefore been recognized as an important drug target for SARS infection. 3CLpro inhibitors have been discovered during the last 12 years from many sources as structure-based molecular docking studies. This perspective is a good guide for researchers to use SARS-CoV 3CLpro inhibitors as anti-SARS chemotherapy [18]. The 3CLpro enzymes are the preferred routes for the immediate treatment of coronaviruses ;and other coronavirus proteins could be targeted for virtual screening [8].
RdRp is the essential enzyme for RNA virus regeneration [19]. In mammalian positive-strand RNA (+RNA) viruses, the enzyme is produced by transcription of the recorded viral genome, which provides a polyprotein precursor that is released by the RdRp-containing subunit by proteolytic cleavage. RdRp is then embedded in a membrane containing viral enzymes that drive the production of new RNA (−RNA) and genome molecules [20]. RNA viral genome is translated with low fidelity because the reproductive product is thought to be untested and unchanged [21]. This property is an important part of the evolution, modification, and proliferation of RNA viruses. 3′-5′ exoribonuclease (nsp14) is one of the replication mechanisms, which interacts with nsp7/nsp8/nsp12.
All these nonstructured proteins help in the replication of long coronavirus RNA. The nsp14 is responsible for the viral RNA capping and the viral genome proofreading to safeguard coronavirus replication fidelity, so all replication machinery proteins that are associated with the RdRp are targets for antiviral drug developments. This complex is important in RNA polymerization, proofreading, and cap-modifying into a multifunctional protein assembly. This macromolecular assembly with RNA helicase (nsp13) shares to form the core of the coronavirus replication machinery, which can engage in coordinated RNA synthesis and processing activities. The RNA helicase (nsp13) is a multifunctional protein, which is able to use the dsDNA and dsRNA as substrates with 5′-3′ polarity, plus it can also work with nsp12 in the viral genome replication machinery and included in viral RNA assembly and associates with nucleoproteins in the membrane complex [5].
We concentrated on the action, receptors, and side effects of six available drugs, Ledipasvir, Sofosbuvir, Ribavirin, Galidesivir, Tenofovir, and Remdesivir. Ledipasvir is approved in combination with Sofosbuvir for the treatment of chronic HCV genotype 1 (Harvoni). Sofosbuvir nowadays is considered as an inhibitor for SARS-CoV-2 based on the similarity between RdRp of HCV and coronaviruses. Both are positive-sense RNA viruses and similar in replication machinery [22]. The dual-component HCV drugs, Harvoni, may be attractive candidates to repurpose because they may inhibit two coronavirus enzymes. Targeting two viral proteins reduces viral resistance. Harvoni has very few side effects and is easy to administer [8]. No toxicity measurements are required because they were tested previously by FDA and no dose recommendation for patients with renal failure [22].
Ribavirin is a nucleotide analog FDA-approved drug for the treatment of chronic HCV in combination with interferon alfa-2b, but this drug is teratogenic, embryonic, and not suitable for cardiac patients. Galidesivir is an approved drug for the treatment of HCV and Ebola virus, which shows broad-spectrum antiviral action against many RNA viruses. Ribavirin is screened virtually against SARS and MERS suggesting strong effectiveness [23]. Tenofovir is an approved drug for HIV-1 and HBV treatment, which is used in adults and the young (2 years of age). But Tenofovir side effects are decreased bone mineral density and immune reconstitution syndrome. COVID-19 treatment with Remdesivir showing improved patients has just been reported, and clinical trial is now underway [24]. The presence of FDA-approved anti-RdRp drugs helps treat patients and reduces the risk of new SARS-CoV-2 infection. Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir are potent drugs against SARS-CoV-2 because they strongly bind to RdRp [25].
The current study found that the molecular docking and interaction between the six available drug candidates with 3CLpro of SARS-CoV-2 show a good interaction and better fit to the 3CLpro enzyme. It is an excellent drug target as it has more than one binding sites. Also, the results obtained by this study showed that Sofosbuvir is the most suitable ligand for 3CLpro with the highest interactions (−7.4 kcal/mol). The interaction affinity of the used six drug ligands in the descending manner is Sofosbuvir, Ledipasvir, and Ribavirin with high interaction (−7.2 kcal/mol), then Galidesivir and Remdesivir are with moderate interaction (−6.3 kcal/mol) and (−6.6 kcal/mol), respectively. Tenofovir shows the least interaction affinity (−6.1 kcal/mol); 3CLpro targeting reported that Harvoni (Ledipasvir /Sofosbuvir) could be promising medication against COVID-19 [8].Another important finding is that all six available drug candidates interact with specific residues in activity groves of SARS-CoV-2 RNA helicase enzyme. The results show that Ledipasvir and Galidesivir have the highest binding affinity (−7.9 kcal/mol). Sofosbuvir shows high affinity (−7.4 kcal/mol) and Remdesivir has slightly good binding affinity of −6.9 kcal/mol. Ribavirin has moderate interaction (−6.4 kcal/mol), but Tenofovir has the least interaction (−6.2 kcal/mol). A surprising finding is that the best results were obtained for SARS-CoV-2 3′-5′ exoribonuclease inhibitors, Ledipasvir and Sofosbuvir, with free binding energies of −10.6 and −10.1 kcal/mol, respectively, followed by Remdesivir and Galidesivir with −9.3 and −7.1 kcal/mol. Ribavirin shows moderate interaction affinity of −6.5 kcal/mol and Tenofovir has the least interaction affinity of −6.1 kcal/mol.
Our study proved the effectiveness of Ledipasvir and Sofosbuvir to inhibit RdRp protein. The promising data were obtained for RdRp inhibitors, Ledipasvir and Sofosbuvir, with free interacting energies of −9.6 and −6.9 kcal/mol, respectively, followed by Tenofovir of slightly good binding affinity of −6.7 kcal/mol; however, Galidesivir and Remdesivir were with moderate binding of −6.6 and −6.4 kcal/mol, respectively. Ribavirin has the least binding power of −6.2 kcal/mol. Modeling and docking for SARS-CoV-2 RdRp and targeting by anti-polymerase drugs such as Sofosbuvir, Ribavirin, and Remdesivir were performed, and it has been reported that these drugs are effective against the newly emerged viral infections such as COVID-19. Many of the FDA-approved anti-RdRp drugs could inhibit COVID-19 with no toxicity measurements required [25]. According to the data explained above, our study found that all six drug candidates bind to the used four proteins effectively and the Harvoni (Ledipasvir and Sofosbuvir) gave the best binding affinity with all proteins in this study. That is means that Ledipasvir and Sofosbuvir have broad-spectrum effect on replication machinery proteins of new emerged SARS-CoV-2.
The strength points of our work: (a) using in-silico prediction of the effect of our six selected available safe FDA-approved antiviral candidates to inhibit four nonstructured proteins of the SARS-CoV-2 replication machinery. (b) These proteins are so important in coronaviruses’ replication, so inhibition of these proteins could stop viral multiplication and infectivity. (c) The selected drugs are marketed, available, and safe drugs especially Ledipasvir and Sofosbuvir with no side effects and no toxicity according to FDA approval sheets. (d) Most of these selected drugs are multiple viral protein inhibitors that prevent the coronavirus antiviral resistance. In-silico study results need to be applicable through experimental and clinical trials. Infected patients should be recruited for clinical trials. Our future study is using the same selected drug candidates and other approved antivirals to be trailed in silico by molecular docking against more nonstructured SARS-CoV-2 proteins such as nsp7 and nsp8.
ConclusionThe disease caused by the novel coronavirus SARS-CoV-2 originated at Wuhan city of China and has been spreading at a fast pace all over the world. Because of the absence of successful treatment choices for SARS-CoV-2, different procedures are being considered. Our study proved in silico that Harvoni (Ledipasvir and Sofosbuvir) is a broad-spectrum anti-SARS-CoV-2 drug with minimal side effects and no toxicity and safe for chronic disease patients with hepatic and renal failure. It is suitable for old and young alike. We recommend using Harvoni in emergency treatment of COVID-19 patients.
Acknowledgements
Authors’ contributions: A.F.E., the corresponding author, molecular docking and in silico study specialist, manuscript drafting. A.T.M., coauthor, manuscript drafting, manuscript revising. S.A., coauthor design of the work, interpretation of the results, manuscript drafting, and revising. W.H., coauthor, design of the study, manuscript revising. All authors have read and approved the manuscript.
Financial support and sponsorship
Nil.
Conflicts of interest
There are no conflicts of interest.
References
Comments (0)