key: cord-0959008-5ozf18lg
authors: Al-Khafaji, Khattab; AL-DuhaidahawiL, Dunya; Taskin Tok, Tugba
title: Using integrated computational approaches to identify safe and rapid treatment for SARS-CoV-2
date: 2020-05-15
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1764392
sha: a8bd381318f31069b5dff2ec23c112b299a3399f
doc_id: 959008
cord_uid: 5ozf18lg

SARS-CoV-2 is a new generation of coronavirus, which was first determined in Wuhan, China, in December 2019. So far, however, there no effective treatment has been found to stop this new generation of coronavirus but discovering of the crystal structure of SARS-CoV-2 main protease (SARS-CoV-2 Mpro) may facilitate searching for new therapies for SARS-COV-2. The aim was to assess the effectiveness of available FDA approved drugs which can construct a covalent bond with Cys145 inside binding site SARS-CoV-2 main protease by using covalent docking screening. We conducted the covdock module MMGBSA module in the Schrodinger suite 2020-1, to examine the covalent bonding utilizing. Besides, we submitted the top three drugs to molecular dynamics simulations via Gromacs 2018.1. The covalent docking showed that saquinavir, ritonavir, remdesivir, delavirdine, cefuroxime axetil, oseltamivir and prevacid have the highest binding energies MMGBSA of –72.17, −72.02, −65.19, −57.65, −54.25, −51.8, and −51.14 kcal/mol, respectively. The 50 ns molecular dynamics simulation was conducted for saquinavir, ritonavir and remdesivir to evaluate the stability of these drugs inside the binding pocket of SARS-CoV-2 main protease. The current study provides a powerful in silico results, means for rapid screening of drugs as anti-protease medications and recommend that the above-mentioned drugs can be used in the treatment of SARS-CoV-2 in combined or sole therapy. Communicated by Ramaswamy H. Sarma

SARS-CoV-2 also called 'severe acute respiratory syndrome coronavirus 2' abbreviated SARS-CoV-2 was recognized to be the causative of atypical pneumonia (Joshi, 2020; Pant et al., 2020) outbreak in Wuhan, China (Hasan, 2020; S. A. Khan et al., 2020) . The virus belongs to the family known as 'coronaviruses' because of the crown-like appearance of spikes glycoproteins on the envelope under an electron microscope . World health organization (WHO) recently announced that the virus transforms from epidemic to pandemic, which requires urgent intervention to prevent the growing spread of the virus across the globe (Chan et al., 2020) . The total confirmed cases 2,347,884 with 738,923 cases in the United State of America (USA) alone and the total death of 161,138 (as of April 19), with mortality estimated within 2% and about 3.4%, according to estimates of approved cases and death worldwide (N. . The most familiar is a virus that arose from the Rhinolophus bat which is > 96% homologous with the modern SARS-CoV-2 virus and it is just 79% homologous with the initial SARS-CoV (Fisher & Heymann, 2020) . The fast-growing number of infected cases globally urged the World Health Organization to announce a state of global health emergency to correlate scientific and medical disciplines to develop rapidly an effective treatment for patients, (Morse et al., 2020; Sarma et al., 2020) elderly patients and people with severe underlying health diseases like heart diseases, lung illness, and diabetic patients, for instance, appear to be at greater risk of revealing severe SARS-CoV-2 requires immediate intervention rather than waiting virus vaccine which may require 1 year to be available (Enayatkhani, 2020) . While drug repurposing could be a short-term and fast resolution to handle SARS-CoV-2 patients (Elfiky, 2020 ; R. J. Kumar et al., 2019) , repurposing existing drugs can offer a good choice to overcome the virus and offer better risk-versus trade-off as compared with discovering new drug and can help overcome time waiting for new therapy rather than use the available resources (Elmezayen, 2020; Muralidharan, 2020) One successful repurposing drug story includes duloxetine which originally developed for depression and FDA approved as the first-inclass choice for stress urinary incontinence (Sweeney & Chancellor, 2005) , duloxetine initially created as antidepressant also is now passed to Phase III clinical trials as a first-in-class treatment for premature ejaculation (McMahon, 2012) and thalidomide, which had a tragic start as an over-the-counter sedative for morning sickness in pregnancy is now being applied to manage leprosy and multiple myeloma (Hideshima & Anderson, 2002) .

As a result of, a newly issued X-ray crystal of SARS-CoV-2 Main protein (Mpro), we planning to use computational approaches (Cameron et al., 2013) to contribute to find an effective treatment for SARS-Cov-2.Thereby, computational analyses speed up these approaches since they allow to handle millions of data simultaneously (Gupta et al., 2020) . Molecular docking includes a set of computational methods and algorithms that aimed to identify novel relationships between chemical ligands and targets through using the modelling of their direct physical interaction (Aanouz, 2020; Ekins et al., 2007) . In present study, we attend to evaluate some of approved drugs to be as covalent binders, irreversible interactions, which can provide a powerful strategy to fight against epidemic viruses. And molecular dynamics simulations can give a more detail for the image which got from molecular covalent docking.

The crystal structure of FDA-available covalent drugs which available in Table 1 were selected based on the review of Kumalo et al. (2015) and some of the antiviral drugs that can form a covalent bond to the target protein. And we aim also here to redirect them for other indications specially to see their possibility to fight against SARS-CoV-2. Thence, we searched about the chosen drugs in PubChem (https://pubchem.ncbi.nlm.nih.gov/) to identify the possibility of the selected drugs to be as covalent binders toward SARS-CoV-2 Mpro. Where, PubChem provides detailed information about the selected drugs rather than other repositories, especially, the property of drug to form covalent bond. Before starting covalent docking, we downloaded the selected drugs one by one and optimize them by using the Ligprep (Lim et al., 2020) based on the OPLS_2005 force field and generated possible state employing Epik in the Schrodinger 2020 (Elfiky 2020; S. A. Khan et al.,2020) . In the next step, the structure of SARS-CoV-2 Mpro (6LU7) was downloaded from Protein Data Bank (http://www.rcsb.org/pdb) (Jin et al. , 2020) . The protease structure was optimized by adding hydrogens, removing water molecules and optimizing charge using the Protein Preparation Wizard module (Kumalo et al., 2015) in Schrodinger suite 2020-1. The covalent docking protocol was preferred since the cysteine 145 residue which is available in the binding site of SARS-CoV-2 Mpro considered as a vital residue that can form a covalent bond with the drug if it can interact covalently. The different mechanisms for the Cys145-FDA drug were mentioned in Table 1 based on the nature of the drug so that the reactive functional group on the ligand and receptor residue are identified and the bond is formed between the correct atoms. These covalently docked complexes were created using Covdock in Schrodinger suite 2020-1. Finally, we selected the lowest the MMGBSA value for each drug as a propriate conformation of the drug inside the binding pocket.

The chosen drugs binding energies were calculated using Prime MM-GBSA modules (Vijayakumar et al., 2014) in the Schr€ odinger (2020). The best poses of selected drugs-SARS-CoV-2 Mpro-were chosen to obtain the binding free energy calculation. Prime MMGBSA is a method that combines Optimized Potential for Liquid Simulations-All Atoms (OPLSAA) force field, molecular mechanics energies (EMM), an SGB solvation model for polar solvation (GSGB), and a non-polar solvation term (GNP) composed of the non-polar solvent accessible surface area and van der Waals interactions. The total binding free energy:

Molecular dynamics simulations are a decision-making process for inspections of protein-drug complexes' stabilities (Al-Khafaji & Taskin Tok, 2020b) . It is used to clarify the dynamic behavior at an atomic level of biological systems, which is hard to handle in labs (Shukla et al., 2019) . In the current study, we conducted molecular dynamics simulations for the top three drugs based on MMGBSA values. The got protein-drug complex structures from covalent docking were submitted to MD simulations (saquinavir, ritonavir, and remdesivir with SARS-CoV-2 Mpro). We hired Gromacs 2018.1 to run 50 ns MD simulations (Abraham et al., 2015) . Charmm 27 force field for all atoms were chosen to run MD simulation (Bjelkmar et al., 2010) . We used Swiss PARAM to produce the topologies of drugs (Zoete et al., 2011) . All protein-drug systems were solvated with three-point transferable intermolecular potential (TIP3P) and their charges were neutralized via adding Na or Cl ions. In the following step, the protein-drugs systems were energetically minimized through the steepest descent algorithm at a tolerance value of 1000 kJ/mol.nm. Then the equilibration with position restraint on the protein molecules for 0.1 ns using NVT and NPT ensembles were done. Electrostatic interactions were evaluated by Particle Mesh Ewald summation (Darden et al., 1993) . We performed the molecular dynamics simulation with no restraint on the protein molecules or ligand to determine the stability in the final step (time step of 0.015 ns). RMSD, RMSF, Rg, and number of hydrogen bonds were chosen to analyze MD trajectories by using GROMACS utilities.

The principal component analysis (PCA) approach was employed to calculate eigenvectors and eigenvalues and their projection along with the first two principal components (Al-Khafaji & Taskin Tok, 2020a). This approach is based on the protocol of GROMACS 2018.1 (Abraham et al., 2015) . It is used to simplify the effect of drugs on the dynamic motion of the targeted protein where it extracts the dynamic motions in simulations that are required for their biological function (Amadei et al., 1993) . We got the principal component analysis from the MD trajectories. A series of eigenvectors and eigenvalues were generated by diagonalizing the matrix. We chose trajectories of the protein backbone of the complexes to get 2 D-projection of motion of trajectory.

To assess the possibility of selected FDA drugs to work as treatment of SARS-CoV-2, the covalent docking was utilized to screen the selected library and rank them according to their binding affinities. The calculated binding free energies of some available drugs using docking score, Glide Gscore, and ensemble-average MM/GBSA are shown in Table 1 .

To validate our covalent docking results, the correlation between MMGBSA and docking score was constructed ( Figure  1 ). The binding energy MMGBSA-docking score correlation shows a good correlation (R 2 ¼ 0.6299). Based on this correlation, we chose the top three ranked-MMGBSA and docking score values of FDA drugs for dissection their binding modes inside the binding site. Covalent docking showed that saquinavir, ritonavir, remdesivir, delavirdine, cefuroxime axetil, oseltamivir, and prevacid have À72.17, À72.02, À65.19, À57.65, À54.25, À51.8, and À51.14, respectively. Where all the topeight FDA drugs show a higher affinity to form a covalent, irreversible bond with Cys145 of SARS-CoV-2 Mpro. Here we investigated the role of molecular weight upon the affinity of selected drugs to bind covalently to Cys145, whereas 57% of selected drugs which have molecular weight over than 600 g/ mol ( Figure 2 ). and have higher free binding energy MMGBSA than À50 kcal/mol. While the ratio of selected drugs decreased to be 16% of drugs which can form covalent bonding with over than À50 kcal/mol. Surprisingly, none of the selected drugs which have a molecular weight below than 300 g/mol can form a good affinity of binding energy. This indicates that higher molecular weight covalent warheads can form stable and efficient binding energy.

What stands out in the Table 1 is saquinavir has the highest binding affinity (lowest binding energy MMGBSA of 72.17 kcal/mol). Therefore, the deep examination of saquinavir is needed. Figure 3 shows that saquinavir not only formed covalent bond of 1.81 Å with Cys145 but also formed five hydrogen bonds inside the pocket of SARS-CoV-2 Mpro. In the second rank ritonavir as presented in Table 1 has À72.02 kcal/mol. This high affinity resulted from covalent bond between ritonavir and Cys145 of 1.82 Å (Figure 4 ) through nucleophilic addition to double bond reaction besides it interacted within binding site by forming three hydrogen bonds. Despite remdesivir formed covalent bond with Cys145 of 1.82 Å and three hydrogen bonds ( Figure 5 ) and this is similar to Ritonavir, but the MMGBSA value is lower than that of ritonavir this may due to nature of reaction which in remdesivir nucleophilic addition to triple bond.

The effect of drug-protein interactions upon dynamics of biological system is a fundamental in drug discovery thereby we used RMSD to investigate the influence of saquinavir, ritonavir, and remdesivir upon the stability of SARS-CoV-2 Mpro. We utilized Gromacs to execute the MD simulations of 50 ns for three drug-protein systems besides of apo protein.

RMSD fluctuations for both apo and hollo forms were measured and presented. The RMSD was calculated to assess the overall dynamics, stability, and convergence of the various systems and the results are presented in Figure 6 (a). Figure 6 (a) shows there is a significant decrease in the RMSD value when SARS-CoV-2 main protease whether it bound to saquinavir, ritonavir, or remedisivir. Further analysis revealed that the RMSD average of apo SARS-CoV-2 main protease was 0.294 nm but when it bound to saquinavir, ritonavir, and remedisivir the RMSD averages were.01865, 0.2130, and 0.2053 nm, respectively. Another significant aspect of MD simulation is the flexibility of protein's backbone which can be assessed through measuring RMSF value. The results of the comparative analysis between these drugs and their effects upon SARS-CoV-2 main protease are illustrated in Figure 6 (b). Closer scrutiny of Figure 6 (b) exhibits the binding of saquinavir diminished the fluctuations of the protein's backbone. And behavior is also can be seen from Figure 6 (b) where the binding of both ritonavir and remedisivir led to reducing the flexibility of the protein. The radius of gyration (Rg) is a definition of system's density, and substantially influences the folding rate and stability of proteins. Rg was employed to assess the compactness of all complexes. In this work, Rg values are in agreement with RMSF values where there are no significant differences between apo form and hollo forms as presented in Figure 6 (c). This reveals that protein remained stable and compact all through the 50 ns time. The number of intermolecular hydrogen bonds is a vital aspect to give an impression about the stability of drug and protein. Also, the number of hydrogen bonds is relevant to the binding scores of molecular docking process. We calculated the hydrogen bonds over time to validate the stable interactions between top three drugs and their corresponding target. As seen in the Figure 6 (d) the Ritonavir has highest order of hydrogen bonding (the average ¼3.35), while saquinavir has an average of 1.68.

The essential dynamic method is a tool to explore the dynamical behavior in the space of SARS-CoV-2 main protease combined with saquinavir, ritonavir, and remedisivir. Basically, the comparison of the drug-bound SARS-CoV-2 Mpro and drug-unbound was made as reference. In order to further understand the configurational space, we selected the first two principal components (PC1, PC2) to analyze their projection of trajectories during the simulations of ligand free and ligand bound SARS-CoV-2 Mpro of the phase space (shown in Figure 7) . During the four system simulations, the results clearly show that the unbound ligand protein covered a wider region of phase space, while all three drug-protein system occupied a smaller region of phase space. Especially, saquinavir reduced the essential dynamics to lowest degree of functional motions as compared with another drugs. Moreover, the PCA results suggest that the drug-bound SARS-CoV-2 Mpro is more stable than ligandunbound SARS-CoV-2 Mpro form of SARS-CoV-2 Mpro. In short, the PCA results are also in agreement with the RMSD and RMSF results, which enhance the validity of the performed analysis.

SARS-CoV-2 causes major pandemic health issue since its spread across the world and can infect humans mainly respiratory system causing severe pneumonia with no vaccine and drug treatment available. Prior study that have referred to the significance of molecular docking to determine effective treatment in short time (Wu, et al., 2020) . In reviewing the literature, we took the advantage of possibility of FDA available drugs that can be as a covalent warhead to inhibit the SARS-CoV-2 Mpro with Cys145. An initial objective of the project was to identify effective and applicable treatment. The present study focused on the main protease (Mpro), especially PDB ID (6LU7) as a potential target for several marketed drugs as possible therapeutic option to combat the virus to see the capability of these drugs to bind with the cysteine 145 residue which is available in the binding site of SARS-CoV-2 main protease. The Mpro in SARS-CoV-2 is necessary for the proteolytic maturation of the virus targeting this protein to limit the expansion of infection by hindering the cleavage of the viral polyprotein. The most interesting finding was that the both saquinavir and ritonavir have the same affinity (binding energy MMGBSA % À72 kcal/mol) to block binding site of SARS-CoV-2 with irreversible interactions. This study supports evidence from previous observations that lopinavir/ritonavir can inhibit SARS-CoV-2 (Lim et al., 2020) . Another important outcome was that the remdesivir comes in second rank with binding affinity (MMGBSA ¼ À65 kcal/mol ). Whereas previous research has established that remdesivir can inhibit SARS-CoV-2 M proenzyme through docking results. Another significant finding, delavirdine computationally showed the ability to form irreversible covalent bond of À57.65 kcal/mol MMGBSA binding energy. Moreover, indinavir and cefuroxime axetil exhibited the possibility to halt the pocket of SARS-CoV-2 Mpro by forming stable interaction through covalent bonding and hydrogen bonding (MMGBSA binding energies % À54 Kcal/ mol). Furthermore, oseltamivir and prevacid exhibited same affinity to bind with SARS-CoV-2 protease through one covalent bond and one hydrogen bond (% À51 Kcal/mol), which is a good explanation for the activity of oseltamivir (Peeri et al., 2020) . In this study, results obtained show the binding affinity of drugs to an active site depends on several factors mainly the ability of a compound to form a covalent bond with amino acid residues of the Mpro (Cys145) and length of a covalent bond, the number of H-bonds that can form with a pocket of the active site and type of nucleophilic addition of unsaturated bonds. In this concept, the structural features in saquinavir and ritonavir like free amine group (-NH2), hydroxyl groups (-OH), carbonyl groups (C¼O) in addition to ether group play key structural feature to form H-bond. Results show the promising activities for antiretroviral drugs saquinavir that used for HIV/AIDS more than ritonavir and remdesivir followed by lopinavir as the best drugs to bind covalently toward SARS-CoV-2 Mpro with lowest energy of binding. As discussed above, saquinavir, ritonavir, remdesivir, delavirdine, cefuroxime axetil, oseltamivir, and prevacid showed the ability to block the binding site of SARS-CoV-2 Mpro by forming covalent bond and stabilized with hydrogen bonding.

These outcomes are contrary to that of Kandeel and Al-Nazawi (2020), they found that ribavirin, telbivudine, vitamin B12, and nicotinamide can be can form non-bonding interactions. Whilst our results showed that, saquinavir, ritonavir, and remdesivir can form irreversible interactions, which are considered an effective way for viral infections.

In the present work, comparative molecular dynamics simulations were conducted to evaluate the effects of saquinavir, ritonavir, and remdesivir on the conformational and dynamical demeanor in either apo or hollo forms to understand the inhibitory possibility at atomic level. An initial mission of the project was to identify rapid, effective, safe, and available for large proportions of people, so we run covalent docking and MMGBSA to sift the possibility of these drug to form irreversible interactions. Therefore, we had chosen these drugs because they carry covalent warheads that can bind covalently to the target. But the docking results are not adequate, so we run MD simulations examine how much these drugs are able to form stable interactions with targeted protein. One interesting finding is the RMSD of protein backbone in apo status has higher average value, whereas the binding with top three drugs diminished these RMSD average. Another important finding was that the RMSF of protein's backbone be less flexible when it compared to apo form of SARS-CoV-2 main protease. It is somewhat surprising that the binding of three drugs have not noted in Rg values. It is not yet clear whether the top three drugs can show an evidence to inhibit the targeted protein, thereby we lean on the principal component analysis to analyze the MD trajectories to judge without dispute. The most obvious finding to emerge from 2D PCA analysis is that the binding of under investigated drugs caused stately impact on essential dynamics of protein by reducing its essential dynamics to its least possible motions. These results seem to be consistent with RMSD results. A possible explanation for this might be that binding of the drugs make the binding site much narrower due to hydrogen bindings this makes the 3D structure of targeted protein more rigid so it will lose its biological functions. It is possible, therefore, that using these available drugs in two ways: either alone or in combined way. These findings suggest the possible use of nominated drugs against SARS-CoV-2 in short time as approved by covalent docking screening. Also, the present results are significant in at least directing the clinicians to use these safe drugs to stop development of Corona virus and second giving hope that available drugs that can be efficient treatment.

Coronavirus today emphasizes as a potential threat to all people worldwide. Although extensive researches have been directed to stop SARS-CoV-2, but till now there is no medication. Meanwhile, the spreading with complicated crisis requires immediate therapy to overcome the spread and minimize mortality of SARS-CoV-2. The aim of the present study was to discover effective treatment through repurposing some of available FDA-approved drugs against SARS-CoV-2 Mpro. Where, they can provide covalent warheads in virtual screening. The most prominent finding to emerge from this study is that the affinity of covalent binder toward SARS-CoV-2 Mpro is ranked: saquinavir > ritonavir > remdesivir > delavirdine > cefuroxime axetil > oseltamivir ¼ prevacid. One of the more noteworthy findings in this study is that MD simulation analysis that saquinavir, ritonavir, and remdesivir can form stable interaction inside the binding site of SARS-CoV-2 Mpro. Also, they restricted the essential motions of protein. Overall, the results of screening toward Mpret encourage for further clinical evaluations. They can be easy to reach and exploit as persuasive treatments for SARS-CoV-2.

Moroccan Medicinal plants as inhibitors of COVID-19: Computational investigations

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

Amygdalin as multi-target anticancer drug against targets of cell division cycle: Double docking and molecular dynamics simulation

Understanding the mechanism of Amygdalin's multifunctional anti-cancer action using computational approach

Essential dynamics of proteins

Implementation of the CHARMM force field in GROMACS: Analysis of protein stability effects from correction maps, virtual interaction sites, and water models

A graph-based recovery and decomposition of Swanson's hypothesis using semantic predications

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Emerging coronaviruses: Genome structure, replication, and pathogenesis

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

In silico pharmacology for drug discovery: Methods for virtual ligand screening and profiling

Anti-HCV, nucleotide inhibitors, repurposing against COVID-19

Drug repurposing for coronavirus (COVID-19): In silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes

Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: An in silico study

Two-dimensional projection of motion of trajectory of SARS-CoV-2 Mpro bound with drugs over the PC1 and PC2

Q&A: The novel coronavirus outbreak causing COVID-19

In-silico approaches to detect inhibitors of the human severe acute respiratory syndrome coronavirus envelope protein ion channel

A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin

Molecular mechanisms of novel therapeutic approaches for multiple myeloma

Structure of Mpro from COVID-19 virus and discovery of its inhibitors

Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease

Virtual screening and repurposing of FDA approved drugs against COVID-19 main protease

Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach

Targeting SARS-CoV-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2 0 -O-ribose methyltransferase

Theory and applications of covalent docking in drug discovery: Merits and pitfalls

Extrapolation of phenolic compounds as multi-target agents against cancer and inflammation

Case of the index patient who caused tertiary transmission of coronavirus disease 2019 in Korea: The application of lopinavir/ritonavir for the treatment of COVID-19 pneumonia monitored by quantitative RT-PCR

Dapoxetine: A new option in the medical management of premature ejaculation

Learning from the past: Possible urgent prevention and treatment options for severe acute respiratory infections caused by 2019-nCoV

Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 protease against COVID-19

Peptide-like and small-molecule inhibitors against Covid-19

The SARS, MERS and novel coronavirus (COVID-19) epidemics, the newest and biggest global health threats: What lessons have we learned?

In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain)

Identification of novel small molecules against GSK3b for Alzheimer's disease using chemoinformatics approach

Treatment of stress urinary incontinence with duloxetine hydrochloride

Identification of natural inhibitors against angiotensin I converting enzyme for cardiac safety using induced fit docking and MM-GBSA studies

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

SwissParam: A fast force field generation tool for small organic molecules

Authors thankful for staff membered of Gaziantep University Institute of Health Sciences, Department of Bioinformatics and Computational Biology.

The authors declare no conflict of interest.