key: cord-323824-74xvvwrw
authors: de Oliveira, Osmair Vital; Rocha, Gerd B.; Paluch, Andrew S.; Costa, Luciano T.
title: Repurposing approved drugs as inhibitors of SARS-CoV-2 S-protein from molecular modeling and virtual screening
date: 2020-06-02
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1772885
sha: 
doc_id: 323824
cord_uid: 74xvvwrw

Herein, molecular modeling techniques were used with the main goal to obtain candidates from a drug database as potential targets to be used against SARS-CoV-2. This novel coronavirus, responsible by the COVID-19 outbreak since the end of 2019, became a challenge since there is not vaccine for this disease. The first step in this investigation was to solvate the isolated S-protein in water for molecular dynamics (MD) simulation, being observed a transition from “up” to “down” conformation of receptor-binding domain (RBD) of the S-protein with angle of 54.3 and 43.0 degrees, respectively. The RBD region was more exposed to the solvent and to the possible drugs due to its enhanced surface area. From the equilibrated MD structure, virtual screening by docking calculations were performed using a library contained 9091 FDA approved drugs. Among them, 24 best-scored ligands (14 traditional herbal isolate and 10 approved drugs) with the binding energy below –8.1 kcal/mol were selected as potential candidates to inhibit the SARS-CoV-2 S-protein, preventing the human cell infection and their replication. For instance, the ivermectin drug (present in our list of promise candidates) was recently used successful to control viral replication in vitro. MD simulations were performed for the three best ligands@S-protein complexes and the binding energies were calculated using the MM/PBSA approach. Overall, it is highlighted an important strategy, some key residues, and chemical groups which may be considered on clinical trials for COVID-19 outbreak.

At the end of December 2019, it was reported an outbreak of many pneumonia cases in Wuhan, China . Hereafter, according to the Chinese Center for Disease Control, this anomalous amount of cases was caused by a novel coronavirus. Initially, this coronavirus was called as 2019-nCoV (2019-novel coronavirus) by World Health Organization (WHO) in January 2020. Now, it has been officially renamed to SARS-CoV-2 (severe acute respiratory syndrome 2) (Guarner, 2020) and the disease caused by this coronavirus is known as COVID-19. Due to the COVID-19 outbreak rapidly spread worldwide, it was declared on 11 March 2020 as a pandemic disease by the WHO. To date at 14 March 2020, there are 216 countries affected by the COVID-19, which 4,307,287 cases and 295,101 deaths were confirmed (https://www.who.int/emergencies/diseases/novel-coronavirus-2019).

The symptoms of COVID-19 may appear between 2 and 14 days after exposure to the SARS-CoV-2 with main signs: cough, runny nose, fever, sore throat, headache and difficulty to breath . The SARS-CoV-2 is belong the SARS-CoV family as verified by Xu et al. (Xu et al., 2020) from comparison of its sequence with others coronavirus. Therefore, SARS-CoV-2 is composed by glycosylated spike protein (S-protein), membrane protein, envelope protein and nucleocapsid protein (Chen et al., 2006; Boopathi et al., 2020) . Likewise, the SARS family, the SARS-CoV-2 infect the cell human using the angiotensin-converting enzyme 2 (ACE2) as receptor, which was confirmed by Ou et al. (Ou et al., 2020; Shang et al., 2020) .Recent molecular modeling study developed by Xu et al. (2020) revealed that the S-protein from SARS-CoV-2 presents strong binding affinity with the ACE2 receptor.

The S-protein has the functional S1 and S2 subunits and they are responsible by the human cell attachment and membrane fusion, respectively (Boopathi et al., 2020) . In the S1 subunit is located the receptor-binding domain (RBD) which binds to human ACE2 and mediates the viral fusion and cellular membranes (Song et al., 2018) . However, the RBD from SARS-CoV-2 has highest affinity with the ACE2 than SARS-CoV as confirmed and evidenced by experimental (Wrapp et al., 2020; Wan et al., 2020) , and theoretical data . In these studies, the authors suggest that the large spread from human to human can be attributed the higher affinity between the SARS-CoV-2 and human ACE2. Moreover, the optimized mutation in the RBD residues has evidenced high-affinity binding with ACE2 and suggest that the SARS-CoV-2 surged from a natural selection process instead of a human manipulation (Andersen et al., 2020) .

Thus, S-protein emerges an interesting target to drug design for the SARS-CoV-2 inhibition, which blocking the RBD by a drug can prevent coupling between the S-protein and ACE2. Consequently, the initial step of fusion of the SARS-CoV-2 and human cell may be avoided inducing the virus death or at least no replication.

Although the S-protein is very important in the initial stage of infection, there is not experimental work reporting their use on drug design. On the other hand, there are various theoretical studies in the literature using the S-protein as target to provide its inhibition. For instance, Zhang et al. carried out docking calculations using Chinese medical herbs to S-protein inhibition, and they found 13 compounds with potential SARS-CoV-2 activity. In another study, the Smith et al. (Smith & Smith, 2020) used molecular dynamics simulation to generate an ensemble docking of the S-protein@ACE2 complex. In their docking calculations, it was found 19 potent candidates to bind at the S-protein:ACE2 interface and 30 top-scoring compounds were predicted to bind in the isolated S-protein. In same way, the Sandeep & McGregor (Sandeep & McGregor, 2020) used the S-protein:ACE2 interface to predict the binding of the hydroxychloroquine and azithromycin using docking calculations. Then Utomo et al. (Utomo et al., 2020) used the Sprotein RBD domain as target in docking calculations using compounds presented in Curcuma sp., Citrus sp., Alpinia galanga, and Caesalpinia sappan. Likewise, herein we used virtual screening and molecular modeling to study the SARS-CoV-2 S-protein with the main insight to obtain potent candidates to inhibit the SARS-CoV-2.

Our approach adopted here differs from those studies in the following way: the isolated S-protein receptor for docking calculations will be obtained from molecular dynamics (MD) simulation, and not directly from crystal structure or S-protein@ACE2 complex. Where, in our understanding, the coupling of full structural flexibility accounted from MD technique and flexible docking calculations provide an efficient strategy to be used in virtual screening. In addition, the investigation of approved drugs is a straightforward way to repurpose quickly a good candidate against the virus. It is interesting to point that, up to date, there is not drug or vaccine against the COVID-19 and the social isolation is one way to control this outbreak as recommended by WHO.

The initial structure of the SARS-CoV-2 viral S-protein was obtained by deleting the ACE2 enzyme from the S-protein@ACE2 complex elaborated by Shah et al. (Shah et al., 2020) . The isolated S-protein was merged into a box with edge of 16.1 Â 16.7 Â 18.5 nm with 147,803 water molecules included. Six chloride ions were added to keep the system neutralized. The GROMOS54a7 force field (Schmid et al., 2011) was used to describe the structural and energetic parameters for the S-protein, and for water the SPC model (Berendsen et al., 1981) was considered. Periodic boundary conditions were considered within 1.0 nm cut-off for nonbonded interactions. The NpT ensemble was used to keep the constant pressure at 1 bar using the Berendsen barostat (Berendsen et al., 1984) (with coupling time of 5 ps) and the temperature at 310 K using the stochastic velocity rescaling method (Bussi et al., 2007) (with coupling time of 0.1 ps). Long-range electrostatic forces were taken into account using the smooth Particle-Mesh Ewald (PME) method (Berendsen, 2007 , Deserno & Holm, 1998 ) with a real space interactions truncated at 1.0 nm cut-off. The system was minimized by the steepest descent method (Morse & Feshbach, 1953) to avoid unfavorable contacts between the atoms, and the convergence was archived at potential energy below 500 kJ/molÁnm. The water molecules were relaxed for 300 ps of MD simulation keeping S-protein rigid by means of a position restraint potential with a force constant of 1000 kJ/molÁnm. Afterwards, the system was equilibrated for 18 ns. All MD simulation steps were performed with the GROMACS package version 2019.2 (Berendsen et al., 1995; Hess et al., 2008) . The total energy, temperature, pressure and the root mean square deviation (RMSD) of the Sprotein were used to monitor the system equilibration. The equations of motions were integrated using the leap-frog algorithm (Berendsen, 2007) with an integration step of 2 fs. The covalent bonds to the hydrogen atoms were fixed using the P-LINCS method (Hess et al., 1997; Hess et al., 2008) .

MD equilibrated structure of the SARS-CoV-2 viral S-protein was used as input for docking calculations. The AutoDockTools software (ADT) (Morris et al., 2009 ) was employed to build the protein-ligand complexes. The AutoDock Vina method (Trott & Olson, 2010) was used in virtual screening calculations for docking the ligands into the S-protein RDB domain. The searching for ligands by docking calculations was done in a grid size of 1.8 Â 4.8 Â 2.2 nm, which was centered at 9.0226 Â 9.3000 Â 2.9897 nm in the RDB region. The S-protein and ligands structures were considered rigid and flexible, respectively, along the docking calculations. The ligands were retrieved from the SWEETLEAD library (Novick et al., 2013) and this was chosen because it was elaborated using approved drugs in USA, India, China, Brazil, WHO Essential Medicines List and others.

Moreover, the use of approved drug in virtual screening calculations is an efficient way to decrease the drug discovery costs and the time spend in research. Recently, Smith et al. (Smith & Smith, 2020) has been used this library to identify small molecules to inhibit the SARS-CoV-2 using a different approach than adopted here. Smith et al. (Smith & Smith, 2020) converted the ligands from SWEETLEAD library into the format (PDBQT) which is accepted in the Vina software. In this conversion, the authors have considered the same isomers of the respective ligands. Thus, here 9091 of these ligands and some their respectively isomers with PDBQT format were downloaded from the supplementary material available by Smith et al. (Smith & Smith, 2020) and used in docking calculations. For each ligand, 20 binding modes were generated. The visualization software PyMol (DeLano, 2009) and Discovery Studio Visualizer (Discovery Studio Visualizer Software, 2020) were used to prepare the figures and to analyze the results.

MD simulations were carried out for the three best-scored ligands complexed to S-protein previously obtained from docking calculations. The ATB version 3.0 (http://compbio. chemistry.uq.edu.au/atb/) (Stroet et al., 2018) was used to generate the GROMOS force field for all ligands. Herein, we used the same protocol adopted in the Molecular Dynamics simulation section. Afterwards, the binding free energy (DG bind. ) was calculated using the MM/PBSA method (Kollman et al., 2000) for each complex along 18 ns of equilibrated MD simulation. For each system, 60 snapshots were extracted at interval of 300 ps along the trajectory.

The DG bind is calculated with the following equation,

where G complex is the Gibbs free-energy of the ligand@S-protein complex, while G S-protein and G ligand are the total G of Sprotein and the ligands in the presence of the solvent, respectively. The individual G for complex, S-protein, and ligands was calculated with equation:

where hE MM i is the average of the potential energy obtained from the MD simulation in vacuum. The hG solv. i is the Gibbs free energy of solvation and its was calculated with the equation,

where G polar is the electrostatic contributions and it is calculated by solving the Poisson-Boltzmann equation (Honig & Nicholls, 1995; Srinivasan et al., 1998; Baker et al., 2001) . The G nonpolar is calculated from the solvent-acessible surface area (SASA). The APBS program (Baker et al., 2001 ) was used to calculate the energy components (G polar and G nonpolar ). The entropy contributions were not taken in account in the equation 3 due its high computational cost. The g_mmpbsa subroutine implemented in the GROMACS package by Kumari (Kumari et al., 2014) was used in the calculations.

To our best knowledge, this work is shedding light on the potential candidates to inhibit the action of the SARS-CoV-2 S-protein on binding to ACE2 and then starting the replication in the human cell.

In this section, relevant results obtained from the MD trajectory are described. As shown in Figure S1 , the potential energy, temperature and pressure were equilibrated. Figure  1 presents the RMSD of the C-alpha atoms of the S-Protein along the MD trajectory, showing that 10 ns was enough to reach a stable and well behaved structure. RMSD values were calculated by the superposition of the S-protein structure along the trajectory with their initial structure (at 0 ps), increasing rapidly until 6 ns with a fluctuation around 0.65 nm and standard deviation of 0.015 nm at end of the simulation. Considering the large size of the protein, the RMSD value of 0.65 nm implies that the S-protein does not suffer significant structural change along the simulation. The inset in Figure 1 shows the RMSF (root mean square fluctuation) values, which indicate lowest fluctuation of the amino acids of the S-protein. The highest RMSF values are attributed the N-and C-terminals, and loop regions. For best comparison, Figure 2 shows the molecular representation of the initial and final snapshots obtained in the MD simulation. Figure 2 reveals that the initial overall S-protein structure remains quite similar after 18 ns of MD simulation, in agreement with the low RMSD values (Figure 1 ). In general, the main difference observed is located in the loop regions, which is expected due to its high flexibility compared to helix motif.

From MD trajectory was observed an important change in RBD region (Figure 2 , highlighted), which is dislocated along the simulation from the "up" to "down" conformation into direction of the S-protein center. Therefore, it indicates that when the RBD is in "down" configuration and it approximates to ACE2, a change is needed to the "up" conformation for an efficient interaction with the ACE2 receptor. In this transition, the angle formed by the C-alpha atoms of the D405, V622 and V991 residues decreases from 54.3 ("up" at 0 ns) to 43.0 degrees ("down" at 18 ns). These data are in agreement with the Peng et al. , which showed that the "up" angle between 84.8 and 52.2 degrees is accessible by ACE receptor and the "down" angle that ranges from 31.6 to 52.2 degrees is not accessible by the ACE2 human receptor. The RMSD value of the RBD region (C336-C525) is 0.33 nm (standard deviation of 0.02 nm) after 6 ns of simulation. In this region, the RMSF values are lowest than 0.32 nm, except for G482, V483, E484, G485 and F486, which remain the RMSF of 0.42, 0.64, 0.60 and 0.50 nm, respectively.

It is important to emphasize again that the initial S-protein model was obtained from the S-protein@ACE2 complex (Shah et al., 2020) . Therefore, the results show that there is a significant structural difference between the RBD in isolated S-protein and that one when complexed with the human ACE2 receptor. It implies that RBD domain is very sensitive to the chemical environment, temperature and solvent effects. Thus, in drug design these differences must be considered in the search of compounds to inhibit the S-protein.

Considering this, in the next section we used the final structure generated in the MD simulation to find good candidates for the SARS-CoV-2 S-protein inhibition using docking calculations.

AutoDock Vina was used for docking calculations of 9091 approved drugs to obtain potential candidates for inhibition of the SARS-CoV-2 S-protein. Differently from other works Smith & Smith, 2020; Sandeep & McGregor, 2020; Utomo et al., 2020) , the isolated S-protein receptor was previously equilibrated by MD simulation, providing a realistic structure when the virus infects humans. Then, the solvent and temperature effects are taking into account in the S-protein conformation. Figures 3 presents the binding affinity or energy for all ligands obtained in Docking calculations. Overall, all ligands docked in the RDB region has favorable energy. It is interesting to notice that the Smith et al. (Smith & Smith, 2020) has been used the present ligands library in docking calculation for S-protein. Therein the isolated protein receptor was retrieved from the S-protein@ACE2 complex previously obtained by their MD simulation. However, our MD results showed that the RBD region suffer significantly structural change from the S-protein@ACE2 complex to the isolated S-protein. For instance, the Smith et al. (Smith & Smith, 2020) identify only 30 compounds with binding affinity between -6.3 and -7.2 kcal/mol. However, as shown in Figure 3 it is possible to identify thousand compounds in this interval of energy, which indicates that our approach is more efficient to explore the S-protein RBD region than the strategy adopted by Smith et al. (Smith & Smith, 2020) . There are 536 ranking docking poses of ligands (including isomers) with binding energy between -7.2 and -8.0 kcal/ mol. However, the ligands with the lowest binding affinity is an important criterion in drug design. In this way, we selected 33 ligands with binding affinity bellow -8.1 kcal/mol, and among them, the isomers with highest energy were removed. In this process, 24 ligands were selected and analyzed here. In the Table 1, we present the results obtained  from the docking calculations. Recently, the azithromycin, hydroxychloroquine and chloroquine have been tested as possible treatment of the COVID-19 outbreak with some potential results (Juurlink, 2020; Gautret et al., 2020; Touret & de Lamballerie, 2020) , but with undesirable side effects. However, our results shows high binding energy of the -6.8, À5.1 and -4.7 kcal/mol, respectively. Therefore, these values indicate that these drugs may act in the S-protein RDB region, but with lowest efficiency than our best-scored ligands in the same concentration. In Table 1 , we can see that 14 compounds are derived from the traditional herbal isolate. The advantages of this type of therapeutics is due the cheap cost, little side effects, good availability and local cultural aspects. On the other hands, 10 approved drugs were selected from the docking calculations. Thus, the laboratories specialized can buy these drugs and test those to verify their efficiency against COVID-9 outbreak. The most scored-ligands presented in Table 1 are rich in hydroxyl group and eight (Lig6982, Lig6982, Lig7089, Lig8121, Lig8283, Lig8402, Lig8541 and Lig8757) has the presence at least one monosaccharides group attached in their structures.

This finding is very interesting because the S-protein is a glycoprotein, which has high affinity with oligosaccharides and various sugar chains are covalently linked. Therefore, we suggest that a drug against SARS-CoV-9 S-protein may have the presence of the oligosaccharide groups, which is important in recognition of S-protein. Analyzing the interacting residues with the S-protein (Table 1) , we found that the R403, R405, Y449, L455, G496 and Y505 residues are common in the most ligands. Therefore, theses residues play an important role in the ligand-S-protein interactions. Thereby, we suggest that in drug design these residues can be used do draw a candidate to S-protein inhibition.

Although the lopinavir and emetine drugs does not appear in the Table 1 , the binding energy obtained in our docking calculations were -7.0 and -6.5 kcal/mol, respectively, see supporting information. These results are interesting because recently the Choy et al. (Choy et al., 2020) showed that these drugs inhibit the SARS-CoV-2 replication with EC 50 under 100 lM. On the other hands, Muralidharan et al. (Muralidharan et al., 2020) obtained from docking calculations a binding energy of -4.1 kcal/mol using the SARS-CoV-2 protease and lopinavir drug, respectively, as receptor and ligand. Thus, the lopinavir may act preferentially in the S-protein than protease, due the lowest binding energy (-7.0 kcal/mol), predicted using the S-protein as receptor.

For best visualization of the interacting residues of the Sprotein and the best-scored ligand, in the Figure 4 is presented the 2 D ligand interaction diagram and molecular electrostatic map (MEP) for the four ligand with lowest binding affinity. In the supporting information ( Figure S2 ), we present the 2 D diagram for all compounds in the Table 1 .

Overall, in the 2 D diagram, the p-p stacking, p-cation and hydrogen bond are the most interactions observed in the ligands@S-protein RBD complexes. For best visualization, in Figure 2 , the MEP was showed only for the RDB region and its position is in front compared with RDB highlighted in the Figure 1 . In the MEP, it can see that all ligands interact with a negative (intense red color) of the RDB region. Contrary, the ligands bind with RBD in a region with partial positive (light blue color) charge.

Along the manuscript written, the ivermectin drug has been tested as anti-SARS-CoV-2 (Caly et al., 2020) and their results showed a reduction of $500 times the cell culture of the virus in 48 h. The authors also highlighted the anti-viral activity of the ivermectin against a broad range of viruses in vitro (Tay et al., 2013; Wagstaff et al., 2012) . Thus, our docking calculation using this drug gives a binding energy of -8.1 kcal/mol. This finding is very promising since the second best-scored ligand also an anti-parasitic agent like the Lig8301 (selamectin) and Lig8568 (doramectin) drugs, Table  1 . Moreover, the main interacting residues with the S-protein are: R403, I418, Y489 and F490. In this way, our calculations indicate that the ivermectin drug may bind in the RBD region, inhibiting the coupling of the SARS-CoV-2 S-protein with the human ACE2 receptor. Recently, Fan et al. (Fan et al., 2020) showed that the selamectin (lig8301) inhibits completely the spike protein of coronavirus GX_P2V, which shares 92.2% amino acid identity with the 2019-nCoV.

Although docking calculations are useful to predict the binding of a ligand into a protein, they are limited due to the lack of temperature and solvent effects, which have an important role in the stability and conformational flexibility of a system. Therefore, we carried out 18 ns of MD simulation for the three best-scored ligands complexed with the S-protein. These simulations were carried out using the Ohio Supercomputer Center (OSC, 1987) . Figure 5 depicts some relevant properties obtained along the MD trajectories.

Overall, Figure 5 highlights that the S-protein is equilibrated at the last 5 ns of MD simulation for all system, as shown by RMSD ( Figure 5(A) ) of the S-protein, in which the values were $0.28 nm with a standard deviation of 0.02 nm for the three complexes. This imply that the presence of the ligands does not change significantly the S-protein structure along the MD simulation. Likewise, the RMSF per residues using their C-alpha atoms ( Figure 5(B) ) was calculated to determine and to identify the flexibility of the S-protein regions. Comparing the RMSF of the ligands@S-protein complexes with the apo S-protein (Figure 1) , it is possible to see that the residues have lowest flexibility in the presence of the ligands.

However, it is important to empathize that we used an equilibrated S-protein as an initial configuration in the MD simulation of the complexes. In addition, all ligands are positioned in the S-protein RDB domain, which is composed by the atoms CA26100 (C336)-CA28081 (C525 residue). In this region is observed low residue fluctuations, keeping the RMSF values below $0.3 nm (Figure 5(B) ). The compactness of the S-protein was achieved by the radius of gyration ( Figure 5(C) ). The S-Protein complexed with Lig8522 and Lig8970 ligands present a similar compactness along the entire MD simulation. For Lig6843, the Rg of the protein remains low fluctuation at the last 7 ns of the trajectory. Overall, the S-protein in all systems (apo e holo forms) is compact with Rg value of $4.35 nm. In general, the RMSD, RMSF and Rg analyzes show that the S-protein has low conformation change at the last 6 ns of the MD trajectories. Thus, as expected, the S-protein present a similar solvent accessible surface area ( Figure 5(C) ), retaining its accessibility in the last 6 ns of the MD simulations.

In Figure 6 , we present the representation of the initial and final structures obtained in simulations.

As one can see in Figure 6 , the ligands are attached on the S-protein RBD along the MD trajectory, in which only the Lig6843 differs from its initial configuration, but it remains in the RBD region. Table 2 shows the binding energetic components obtained in the MM/PBSA calculations.

Although the binding energy obtained from MM/PBSA method differs from the docking calculations, they follow the same energetic order. Lig8522 and Lig8970 present the same binding energy with similar van der Waals and apolar solvation (SASA) energies. However, the electrostatic and polar solvation energies are very different due to the large negatively charge (-6e) of the Lig8970 salt. The positive value of the binding energy for the Lig6843 was an unexpected result. However, this result is a consequence of the Lig6843 reorientation along the MD simulation, see Figure 6 . In this process, the initial position of the Lig6843 (charge þ2e) is favorable for interactions with two negative regions of the RBD (Figure 4) . However, in the MD trajectory it is dislocated for an unfavorable positive region. Moreover, the large error in the binding energy of this ligand indicate that is necessary more time to be simulated in the MD simulation. Therefore, the Lig6841 remain binding in the RBD region by the van der Waals interactions along the 18 ns of the MD simulation. Nowadays, we are performing long MD simulations for the best-scored ligands presented in the Table 1 . Forthcomingly, when in vitro activity were disposable in the literature, we will do receiver operating characteristic (ROC) curve analysis from machine learning to identify possible true positive ligands.

Currently, we are surviving at a global health crisis caused by the SARS-CoV-2 virus, which is responsible by the COVID-19 outbreak since the end of the 2019. To date, there is not a vaccine against COVID-19 and the social isolation is the main strategy adopted to avoid the spread of this novel coronavirus. Therefore, herein we used molecular dynamics (MD) simulation and docking calculations to study the SARS-CoV-2 S-protein with the main goal to obtain possible drugs candidates for repurposing them against to COVID-19. Our approach was study the SARS-CoV-2 S-protein before it bounded with the human angiotensin-converting enzyme 2 receptor preventing the human cell infection. Overall, the Sprotein does not undergo significantly structural change along the 18 ns MD simulation equilibration in water solvent. The main change occurs in the receptor-binding domain (RBD) which was converted from "up" (angle of 54.3 degrees) to "down" (angle of 43.0 degrees) conformation. In this transition, the RBD increased surface area enhanced the possibility of interactions with the solvent and also with the possible drugs. In this way, we carried out virtual screening by docking calculations using the S-protein at 18 ns of the MD simulation as receptor and 9091 approved drugs as ligands. In these calculations, 14 traditional herbal isolate and 10 approved drugs were obtained as a potent candidates to inhibit the SARS-CoV-2 S-protein with binding affinity bellow -8.1 kcal/mol. In addition, our calculations were successful to predict the binding of the ivermectin in ACE2, which this drug was used recently with high efficiency to control viral replication in vitro. Overall, the R403, R405, Y449, L455, G496 and Y505 residues are the most interacting with ligands. MD simulations of the three ligands@complexes show that the ligands remain bounded in the RBD region along the simulation. Where the binding energies obtained from MM/MBSA calculations were -38.51 (± 1.59), À40.43 (± 1.92) and 2.67 (± 1.18) kcal/mol for Lig8522, Lig8970 and Lig6843, respectively. Finally, herein we highlight some important issues to be taken into account for repurposing known drugs for the COVID-19 outbreak. 

The proximal origin of SARS-CoV-2

Electrostatics of nanosystems: Application to microtubules and the ribosome

Intermolecular Forces

Molecular dynamics with coupling to an external bath

GROMACS: A message-passing parallel molecular dynamics implementation

Simulating the physical world: hierarchial modeling from quantum mechanics to fluid dynamics

Novel 2019 coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment

Canonical sampling through velocity rescaling

The FDA-approved Drug Ivermectin inhibits the replication of SARS-CoV-2 in vitro

An overall picture of SARS coronavirus (SARS-CoV) genome-encoded major proteins: Structures, functions and drug development

Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro

How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines

Repurposing of clinically approved drugs for treatment of coronavirus disease 2019 in a 2019-novel coronavirusrelated coronavirus model

Hydroxychloroquine and azithromycin as a treatment of COVID-19: Results of an open label non-randomized clinical trial

Three emerging coronaviruses in two decades

LINCS: A linear constraint solver for molecular simulations

GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation

P-LINCS: A parallel linear constraint solver for molecular simulation

Classical electrostatics in biology and chemistry

Safety considerations with chloroquine, hydroxychloroquine and azithromycin in the management of SARS-CoV-2 infection

Calculating structures and free energies of complex molecules: Combining molecular mechanics and continuum models

g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations

New coronavirusinfected pneumonia Engulfs Wuhan

Autodock4 and AutoDockTools4: Automated docking with selective receptor flexibility

Asymptotic series, method of steepest descent

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

SWEETLEAD: An in silico database of approved drugs, regulated chemicals, and herbal isolates for computer-aided drug discovery

Ohio supercomputer center

Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV

Exploring the binding mechanism and accessible angle of SARS-CoV-2 spike and ACE2 by molecular dynamics simulation and free energy calculation

Energetics based modeling of hydroxychloroquine and azithromycin binding to the SARS-CoV-2 spike (S)Protein-ACE2 Complex

Definition and testing of the GROMOS force-field versions 54A7 and 54B7

Structural basis of receptor recognition by SARS-CoV-2

Sequence variation of SARS-CoV-2 spike protein may facilitate stronger interaction with ACE2 promoting high infectivity

Repurposing therapeutics for COVID-19: Supercomputer-based docking to the SARS-CoV-2 viral spike protein and viral spike protein-human ACE2 interface

Cryo-EM structure of the SARS coronavirus spike glycoprotein in complex with its host cell receptor ACE2

Continuum solvent studies of the stability of DNA, RNA, and phosphoramidate -DNA helices

Automated topology builder version 3.0: Prediction of solvation free enthalpies in water and hexane

Nuclear localization of dengue virus (DENV) 1-4 non-structural protein 5; protection against all 4 DENV serotypes by the inhibitor Ivermectin

Of chloroquine and COVID-19

The PyMOL Molecular Graphics System

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading

Revealing the potency of citrus and galangal constituents to halt SARS-CoV-2 infection

Ivermectin is a specific inhibitor of importin a/b-mediated nuclear import able to inhibit replication of HIV-1 and dengue virus

Receptor recognition by the novel coronavirus from Wuhan: An analysis based on decade-long structural studies of SARS coronavirus

Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation

The SARS-CoV-2 outbreak: What we know

Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission

In silico screening of Chinese herbal medicines with the potential to directly inhibit 2019 novel coronavirus

The authors also acknowledge the computing support from the Ohio Supercomputing Center for the binding free energy analysis calculations. We are grateful for the researcher Charly Empereur-mot for his valuable discussions.

No potential conflict of interest was reported by the author(s).