The Binding mechanism of Ivermectin and levosalbutamol with spike protein of SARS-CoV-2

doi:10.21203/rs.3.rs-160254/v1

Download PDF

Research Article

The Binding mechanism of Ivermectin and levosalbutamol with spike protein of SARS-CoV-2

https://doi.org/10.21203/rs.3.rs-160254/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 11 Apr, 2021

Read the published version in Structural Chemistry →

You are reading this latest preprint version

In this study, we have investigated the binding mechanism of two FDA approved drugs (ivermectin and levosalbutamol) with the spike protein of SARs-CoV-2 using three different computational modeling techniques. Molecular docking results predict that ivermectin shows a large binding affinity for spike protein (-9.0 kcal/mol) compared to levosalbutamol (-4.1 kcal/mol). Ivermectin binds with GLN492, GLN493, GLY496 and TRY505 residues in the spike protein through hydrogen bonds and levosalbutamol binds with TYR453 and TYR505 residues. Using density functional theory (DFT) studies, we have calculated the binding energies between ivermectin and levosalbutamol with residues in spike protein which favor their binding are − 17.8 kcal/mol and − 20.08 kcal/mol, respectively. The natural bond orbital (NBO) charge analysis has been performed to estimate the amount of charge transfer that occurred by two drugs during interaction with residues. Molecular dynamics (MD) study confirms the stability of spike protein bound with ivermectin through RMSD and RMSF analyses. Three different computer modeling techniques reveal that ivermectin is more stable than levosalbutamol in the active site of spike protein where hACE2 binds. Therefore, ivermectin can be a suitable inhibitor for SARS-CoV-2 to enter into the human cell through hACE2.

Medicinal Chemistry

Computational Chemistry

SARS-CoV-2

covid19

hACE2

docking

ivermectin

MD simulation

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infests the respiratory tract and causes coronavirus disease 2019 (COVID-19). It is a new member of the corona virus family which is closely related to SARS-CoV in genetic sequencing.[1, 2] SARS-CoV-19 was first recognized in December 2019 which was globally spread and infected 80,773,033 people including 1,783,619 deaths according to a report of “World Health Organization” (9.32 pm pm GMT+6, 30 December 2020). The cause of the worldwide outbreak of this virus is the ability to transmit man to man by contact with infected droplets of saliva, cough, or sneeze.[3] Its high infective potential is disrupting the economic and social life of every single people worldwide. It is reported that in 20% of infected people, COVID-19 spreads in the lungs and other organs of the body. But no drug or vaccine has been approved against COVID-19. A few FDA-approved drugs are being used to COVID-19 patients at the early stage of infection which may cause less severe infection due to reduced entry of COVID-19 into the cell. Ivermectin, an FDA-approved antiparasitic drug can reduce the replication of SARS-CoV-19 in vitro.[4] Many research groups have proved that SARS-CoV-2 utilizes the homotrimeric spike (S) glycoprotein to enter into the human cell through binding to the functional receptor human angiotensin-converting enzyme 2 (hACE2).[5, 6] The receptor-binding motif (RBM) in S protein of SARS-CoV-2 has two regions that form the interface between S protein and hACE2.[7] In this interface, the existing L455 and Q493 of the SARS-CoV-2 have favorable interaction with hACE2 and K31 and E35 in hACE2 form hydrogen bonds with Q493 to enhance binding S protein.[5] If any drug molecule can inhibit the binding between SARS-CoV-2 and hACE2 by blocking S protein's interface, especially Q493 and L455, it will prohibit the SARS-CoV-2 entry into the human cell.

In the present study, we investigated the interaction of ivermectin with the interface region of the spike protein through which it binds with hACE2. We also tried to understand which amino acid residues in the spike protein are favorable to bind with ivermectin and how they interact? The molecular docking was conducted to predict the active site at the S protein interface for ivermectin and their potential binding affinity and interactions. Density functional theory (DFT) was carried out to calculate the binding energy between interacting residues and ivermectin and the amount of charge transfer during an interaction. The Hirshfeld surface analysis was performed to understand the intermolecular interaction and identification of close contact atoms in the ligand-protein system.[8] We considered one more FDA-approved drug levosalbutamol, a β2 adrenergic receptor with a poor antiviral effect.[9] This drug is used for asthma and chronic obstructive pulmonary diseases (COPD).[10]

2.1 Molecular Docking: The crystal structure of the spike protein of COVID-19 was recovered from the protein data bank (PDB:6 m0j) at 2.45 Å resolution. The water molecules and heteroatoms were removed and polar hydrogen atoms were added to the amino acid residues. The initial structure of ivermectin was obtained from the Pubchem database and levosalbutamol was constructed through Gaussview. All prepared PDB files of protein and drug molecules were processed to PDBQT format for reckoning energy grid map.

The binding affinities of these two drugs with S protein were elucidated by AutoDock Vina. The docking grid was set around the region1 and region2 of spike protein (Fig. 1) with center X=-37.33, Y = 31.93 and Z = 2.68 and the dimension were X = 29.50, Y = 44.97 and Z = 19.33. The docking process was repeated 10 times and the mean affinity score was considered as binding affinity. To visualize the docking results PyMol and Discovery Studio Visualizer were utilized.

2.2 DFT calculations: Using docking results all amino acid residues of the spike protein were selected out which favored two drugs to bind. The coordinates of selected amino acid residues of S protein and drug molecule were considered as an initial configuration to calculate binding energies. Equilibrium structures with local minima in potential energy profiles were computed using density functional theory (DFT) with Becke three-parameters exchange functional[11] and the Lee-Yang-Parr functional[12] (B3LYP) with the all electron 6–31 g(d) basis set. The vibrational frequencies were calculated and no imaginary frequencies were found for the optimized structures. To get binding energies, all selected amino acid residues and drug molecules were optimized separately. The binding energy was calculated using Eq. 1;[13]

$${E}_{BE}={E}_{res-lig}-{E}_{res}-{E}_{lig}$$

where ${E}_{BE}$ is the binding energy between selected amino acid residues and drug molecules. ${E}_{res-lig}$ is the total energy of residues-drug complex after binding. ${E}_{res}$ and ${E}_{lig}$ are total energies of residues and drug molecules, respectively. All calculations are carried out at the gas phase with Gaussin16 program[14].

2.3 Molecular dynamics: Molecular dynamics (MD) simulations were performed to the conformations of S protein bound with ivermectin and levosalbutamol with the highest docking score, separately. Prior to MD simulation, the complexes were immersed in a cubic periodic box and solvated in 19,533 simple point charge (SPC) water molecules during the topology process. Two chlorine ions were added to neutralize the system. This system maintained 1.0 nm distance from the box edge and was subjected to energy minimization using 5000-step steepest descent minimization. 100 ps NVT followed by 100 ps NPT ensemble at 300 K and 1.0 bar were performed using Berendsen thermostat and barostat[15] to equilibrate water molecules around the complex. The particle mesh Ewald (PME)[16] algorithm was used to consider the long-range Coulomb interactions. For van der Waals interactions, Lennard-Jones potentials were used. A cutoff 1.2 nm was applied for considering both long and short range interactions. The LINCS algorithm was applied to constrain all bonds.[17] A 100 ns long MD simulation was performed in the NVT ensemble at 300 K with a time step of 2.0 fs and 10.0 ps interval final coordinates were saved. MD simulation was performed using the Gromcas 5.1.4[18] and force field Charmm36[19].

3.1 Molecular Docking: The highest binding affinity of ivermectin and levosalbutamoland at the interface of S protein were calculated − 9.0 and − 4.1 kcal/mol, respectively. Figure 2 shows the binding pocket of ivermectin with S protein. Here, residues GLN493 and LEU492 are involved in hydrogen bond (HB) formation with ivermectin. The hydrogen atom of NH₂ group of GLN493 forms two HBs with the oxygen atom of -CO and oxygen atom of –OH groups in ivermectin and the oxygen atom of –CO group in LEU492 forms one HB with the hydrogen atom of –OH of ivermectin. The HB plays a vital role in stabilizing the protein-ligand complex. Besides these, oxygen atom of –CO groups in TYR505 and GLY496 form HBs with the hydrogen atom of alkyl group in ivermectin. GLU484, PHE497, ASN501, GLN498, TYR453, SER494, AGR403 and TYR495 interact with ivermectin through van der Waals interaction. Another stacking interaction between aromatic rings of residues (TYR449, TYR489, PHE456 and PHE490) and alkyl group of ivermectin plays a significant role to strengthen the binding of ivermectin with S protein, including hydrogen bonds and van der Waals interaction. Among all interacting residues, GLN493 and LEU455 residues which have favorable interactions with ACE2, can bind with ivermectin. Therefore, before the interaction between S protein and hACE2, it can block the S protein's binding pocket, which may inhibit SARS-CoV-2 binding with hACE2.

Figure 3 shows the interaction between S protein and levosalbutamol. TYR453 and TYR505 are involved to form two HBs with levosalbutamol. The oxygen atoms of –OH and –CO group in TYR453 and TYR505 form HBs with hydrogen atoms of –OH groups in levosalbutamol, respectively. Another stacking interaction is found between aromatic rings of residues (TYR453 and TYR505) and alkyl carbon of levosalbutamol. GLY496 residue interacts with levosalbutamol through the pi-donor hydrogen bond.

3.2 DFT Calculations: The initial geometries of ivermectin and levosalbutamol were optimized separately at the above mentioned level of theory. We considered only those amino acid residues of S protein that interact with either ivermectin or levosalbutamol through hydrogen bonds to reduce simulation time. Figure 4 shows the optimized configuration of four residues (GLN492, GLN493, GLY496 and TRY505), which interact with ivermectin through hydrogen bonds. The calculated binding energy of these four residues with ivermectin is -17.8 kcal/mol. Oxygen and hydrogen atoms of –CO and –OH groups of ivermectin form hydrogen bonds with –NH group of GLN493 and –OH group of LEU492, respectively. Alkyl hydrogens of ivermectin form two hydrogen bonds with GLY496 and TYR505. Hydrogen bonds (–CO- -NH-) and (–OH- -OC-) are shorter in distance compared to those forms with alkyl hydrogens of ivermectin. The hydrogen bond distances are shown in Fig. 4.

Similarly, two residues TYR453 and TYR505 with levosalbutamol were considered for DFT calculation because they formed hydrogen bonds with levosalbutamol. Their calculated binding energy is -20.08 kcal/mol. The hydrogen atom of –OH and the oxygen atom of –CO group of TYR453 and TYR505 form hydrogen bonds with –OH groups in levosalbutamol, respectively. Their hydrogen bond distances are 1.94 Å and 1.88 Å, respectively (Fig. 5).

The amount of charge transfers between two interacting species that occurs for electron transfer establishes the strength of the interaction between them. We have analyzed NBO charge transfer to compute the number of charges transferred by ivermectin and levosalbutamol from or to amino acid residues. Ivermectin donates 0.00845e⁻ to four residues (GLN492, GLN493, GLY496 and TRY505) of S protein. On the other hand, levosalbutamol donates 0.01586e⁻ to two residues (TYR453 and TYR505).

3.3 Stability of spike protein bound with ivermectin and levosalbutamol during MD simulation: To study the stability of the initial conformation of docked complex in an aqueous solution, the dynamics of the complex such as RMSD (root mean square deviation) and RMSFs (root mean square fluctuations) were examined. For this, a 100 ns long simulation was performed on the docked structure of the spike protein-ivermectin complex.

The RMSD computed for the backbone of protein (Fig. 6) shows the relative stability of the spike protein bound with ivermectin and levosalbutamol. The RMSD values for all structures were calculated using the mass-weighted least-squares fitting method. In Fig. 6, the black color represents the RMSD of the backbone of the protein bound with ivermectin, which reaches equilibrium within 2.0 ns and fluctuates between 0.9 nm to 0.12 nm up to 70 ns. During this period, the average RMSD value is ~ 0.11 nm. After 70 ns, the protein depicts a sudden drift towards a little higher deviation at ~ 0.18 nm. During this period, the RMSD fluctuates from 0.12 nm to 0.17 nm and the average RMSD is 0.15 nm. This lower deviation indicates the stability of protein bound with ivermectin in the active site throughout the simulation. On the other hand, the RMSD (red color) shows a large fluctuation for the protein bound with levosalbutamol. This is beacuse levosalbutamol moves away from the active site of the spike protein. Therefore, it confirms that ivermectin can bind strongly in the active site of spike protein than levosalbutamol.

To understand the dynamics in more details, we compute the time-averaged RMSF (Fig. 6) for the spike protein bound with ivermectin (black color) and levosalbutamol (red color) at 300 K. The RMSF is used to evaluate the flexibility of the residues of protein which reveals the extent to which fluctuations vary along the protein chain before and after binding with ivermectin. The RMSF values of GLN492, GLN493, GLY496 and TRY505 are found close to the baseline (~ 0.06 nm) except GLN493 and TRY505 (< 0.15 nm). This small fluctuation of these four important amino acids will be beneficial for binding ivermectin in S protein's active site. The large fluctuation (~ 0.45 nm) is found for amino acids 369–371 which are relatively far from the active site of ivermectin. In the case of spike protein bound with levosalbutamol, the RMSF values for TYR453 and TYR505 are found 0.06 and 0.1, respectively.

Using three different computational modeling techniques, we have investigated the binding mechanism of two FDA approved drugs (ivermectin and levosalbutamol) in the active site of the spike protein of SARs-CoV-2 where hACE2 binds. Molecular docking results predict that ivermectin shows a large binding affinity for spike protein (-9.0 kcal/mol) compared to levosalbutamol (-4.1 kcal/mol). Ivermectin binds with GLN492, GLN493, GLY496 and TRY505 residues in spike protein which favor hACE2 to bind in the active site. Ivermectin binds with these residues through hydrogen bonds. Levosalbutamol binds with TYR453 and TYR505 residues. Using density functional theory (DFT) studies, we have calculated the binding energy between ivermectin and four residues in the active site -17.8 kcal/mol. In the case of levosalbutamol with two residues, it is -20.08 kcal/mol. The natural bond orbital (NBO) charge analysis has been performed to estimate the amount of charge transfer that occurred by two drugs during interaction with residues. Ivermectin donates 0.00845e^- to four residues (GLN492, GLN493, GLY496 and TRY505) of the spike protein. On the other hand, levosalbutamol donates 0.01586e^- to two residues (TYR453 and TYR505). Molecular dynamics (MD) study confirms the spike protein's stability with ivermectin through RMSD and RMSF analyses. During 100 ns long MD simulation, ivermectin was bound in the spike protein's active site. On the other hand, levosalbutamol moves away from that site. Three different computer modeling techniques confirm that ivermectin is more stable than levosalbutamol in the active site of spike protein where hACE2 binds. Therefore, ivermectin can be a suitable inhibitor for SARS-CoV-2 to enter into the human cell through hACE2.

Funding: This study was supported by Research Grant (No. 16-008 RG/CHE/AS_I – FR3240293317) funded by The World Academy of Sciences (TWAS). The authors thank TWAS for this generous support.

Conflict of interest: The authors declare that they have no conflict of interest.

Availability of data and material: Not applicable

Code availability: Not applicable

Authors’ contribution: Md. Jahir Raihan performed data analysis and figure preparation. Study conception, design and manuscript writing were performed by Joyanta K. Saha.

Gorbalenya AE, Baker SC, Baric RS, de Groot RJ, Drosten C, Gulyaeva AA, Haagmans BL, Lauber C, Leontovich AM, Neuman BW, Penzar D, Perlman S, Poon LLM, Samborskiy DV, Sidorov IA, Sola I, Ziebuhr J (2020) V. Coronaviridae Study Group of the International Committee on Taxonomy of, The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nature Microbiology 5:536–544
Yang P, Wang X (2020) COVID-19: a new challenge for human beings. Cell Mol Immunol 17:555–557
Hagar M, Ahmed HA, Aljohani G, Alhaddad OA, Investigation of Some Antiviral N-Heterocycles as COVID 19 Drug: Molecular Docking and DFT Calculations, International Journal of Molecular Sciences, 21 (2020)
Caly L, Druce JD, Catton MG, Jans DA, Wagstaff KM (2020) The FDA-approved drug ivermectin inhibits the replication of SARS-CoV-2 in vitro. Antiviral Res 178:104787
Yi C, Sun X, Ye J, Ding L, Liu M, Yang Z, Lu X, Zhang Y, Ma L, Gu W, Qu A, Xu J, Shi Z, Ling Z, Sun B (2020) Key residues of the receptor binding motif in the spike protein of SARS-CoV-2 that interact with ACE2 and neutralizing antibodies. Cell Mol Immunol 17:621–630
Ragia G, Manolopoulos VG, Inhibition of SARS-CoV-2 entry through the ACE2/TMPRSS2 pathway: a promising approach for uncovering early COVID-19 drug therapies, European Journal of Clinical Pharmacology, (2020)
Morse JS, Lalonde T, Xu S, Liu WR, Learning from the Past: Possible Urgent Prevention and Treatment Options for Severe Acute Respiratory Infections Caused by 2019-nCoV, ChemBioChem, 21 (2020) 730–738
Spackman PR, Yu L-J, Morton CJ, Parker MW, Bond CS, Spackman MA, Jayatilaka D, Thomas SP (2019) Bridging Crystal Engineering and Drug Discovery by Utilizing Intermolecular Interactions and Molecular Shapes in Crystals. Angew Chem Int Ed 58:16780–16784
Touret F, Gilles M, Barral K, Nougairède A, van Helden J, Decroly E, de Lamballerie X, Coutard B (2020) In vitro screening of a FDA approved chemical library reveals potential inhibitors of SARS-CoV-2 replication. Sci Rep 10:13093
Jat KR, Khairwa A (2013) Levalbuterol versus albuterol for acute asthma: a systematic review and meta-analysis. Pulm Pharmacol Ther 26:239–248
Becke AD (1993) Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys 98:5648–5652
Lee C, Yang W, Parr RG (1988) Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B 37:785–789
Saha JK, Hossain MS, Ghosh MK (2019) DFT study of response mechanism and selectivity of poly(3,4-ethylenedioxythiophene) towards CO2 and SO2 as gas sensor. Struct Chem 30:1427–1436
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Petersson GA, Nakatsuji H, Li X, Caricato M, Marenich AV, Bloino J, Janesko BG, Gomperts R, Mennucci B, Hratchian HP, Ortiz JV, Izmaylov AF, Sonnenberg JL, Williams F, Ding F, Lipparini F, Egidi J, Goings B, Peng A, Petrone T, Henderson D, Ranasinghe VG, Zakrzewski J, Gao N, Rega G, Zheng W, Liang M, Hada M, Ehara K, Toyota R, Fukuda J, Hasegawa M, Ishida T, Nakajima Y, Honda O, Kitao H, Nakai T, Vreven K, Throssell JA, Montgomery JE Jr, Peralta F, Ogliaro MJ, Bearpark JJ, Heyd EN, Brothers KN, Kudin VN, Staroverov TA, Keith R, Kobayashi J, Normand K, Raghavachari AP, Rendell JC, Burant SS, Iyengar J, Tomasi M, Cossi JM, Millam M, Klene C, Adamo R, Cammi JW, Ochterski RL, Martin, K. Morokuma, O. Farkas, J.B. Foresman, D.J. Fox, Gaussian 16, in, Wallingford, CT, 2016
Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR (1984) Molecular dynamics with coupling to an external bath. J Chem Phys 81:3684–3690
Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG (1995) A smooth particle mesh Ewald method. J Chem Phys 103:8577–8593
Yuan X-H, Wang Y-C, Jin W-J, Zhao B-B, Chen C-F, Yang J, Wang J-F, Guo Y-Y, Liu J-J, Zhang D, Gong L-L, He Y-W (2012) Structure-Based High-Throughput Epitope Analysis of Hexon Proteins in B and C Species Human Adenoviruses (HAdVs). PLOS ONE 7:e32938
Pronk S, Páll S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, Shirts MR, Smith JC, Kasson PM, van der Spoel D, Hess B, Lindahl E (2013) GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics 29:845–854
Best RB, Zhu X, Shim J, Lopes PEM, Mittal J, Feig M, MacKerell AD (2012) Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone ϕ, ψ and Side-Chain χ1 and χ2 Dihedral Angles. J Chem Theory Comput 8:3257–3273

Download PDF

Journal Publication

published 11 Apr, 2021

Read the published version in Structural Chemistry →

Editorial decision: Minor revisions
21 Mar, 2021
Reviewers invited by journal
26 Feb, 2021
Reviews received at journal
26 Feb, 2021
First submitted to journal
14 Jan, 2021

You are reading this latest preprint version

The Binding mechanism of Ivermectin and levosalbutamol with spike protein of SARS-CoV-2

Status:

Journal Publication

Version 1

Abstract

Figures

1. Introduction

2. Computational Methodology

3. Results And Discussions

4. Conclusion

Declarations

References

Status:

Journal Publication

Version 1