Next Article in Journal
FLASH Radiotherapy: Mechanisms of Biological Effects and the Therapeutic Potential in Cancer
Previous Article in Journal
The Role of Hyperuricemia in Cardiac Diseases: Evidence, Controversies, and Therapeutic Strategies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 14
Issue 7
10.3390/biom14070755
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Comparative Analysis of Ivermectin and Nirmatrelvir Inhibitors Interacting with the SARS-CoV-2 Main Protease

by
Yuri Alves de Oliveira Só
1,
Katyanna Sales Bezerra
2,
Ricardo Gargano
1,*,
Fabio L. L. Mendonça
3,
Janeusa Trindade Souto
4,
Umberto L. Fulco
2,
Marcelo Lopes Pereira Junior
3 and
Luiz Antônio Ribeiro Junior
1,5
1
Institute of Physics, University of Brasília, Brasília 70910-900, Brazil
2
Department of Biophysics and Pharmacology, Federal University of Rio Grande do Norte, Natal 59078-570, Brazil
3
Department of Electrical Engineering, College of Technology, University of Brasília, Brasília 70910-900, Brazil
4
Department of Microbiology and Parasitology, Biosciences Center, Federal University of Rio Grande do Norte, Natal 59064-741, Brazil
5
Computational Materials Laboratory, University of Brasília, LCCMat, Brasília 70919-970, Brazil
*
Author to whom correspondence should be addressed.
Biomolecules 2024, 14(7), 755; https://doi.org/10.3390/biom14070755
Submission received: 25 April 2024 / Revised: 4 June 2024 / Accepted: 6 June 2024 / Published: 25 June 2024
(This article belongs to the Section Chemical Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Exploring therapeutic options is crucial in the ongoing COVID-19 pandemic caused by SARS-CoV-2. Nirmatrelvir, which is a potent inhibitor that targets the SARS-CoV-2 Mpro, shows promise as an antiviral treatment. Additionally, Ivermectin, which is a broad-spectrum antiparasitic drug, has demonstrated effectiveness against the virus in laboratory settings. However, its clinical implications are still debated. Using computational methods, such as molecular docking and 100 ns molecular dynamics simulations, we investigated how Nirmatrelvir and Ivermectin interacted with SARS-CoV-2 Mpro(A). Calculations using density functional theory were instrumental in elucidating the behavior of isolated molecules, primarily by analyzing the frontier molecular orbitals. Our analysis revealed distinct binding patterns: Nirmatrelvir formed strong interactions with amino acids, like MET49, MET165, HIS41, HIS163, HIS164, PHE140, CYS145, GLU166, and ASN142, showing stable binding, with a root-mean-square deviation (RMSD) of around 2.0 Å. On the other hand, Ivermectin interacted with THR237, THR239, LEU271, LEU272, and LEU287, displaying an RMSD of 1.87 Å, indicating enduring interactions. Both ligands stabilized Mpro(A), with Ivermectin showing stability and persistent interactions despite forming fewer hydrogen bonds. These findings offer detailed insights into how Nirmatrelvir and Ivermectin bind to the SARS-CoV-2 main protease, providing valuable information for potential therapeutic strategies against COVID-19.

Graphical Abstract

1. Introduction

Emerging in Wuhan, China, toward the end of 2019, a new acute respiratory illness caused by the novel coronavirus (SARS-CoV-2) was officially named COVID-19. The World Health Organization (WHO) declared it a global pandemic on 11 March 2020 [1]. This virus is known for its high transmissibility. It primarily spreads through human-to-human contact via contaminated respiratory secretions released during coughing, sneezing, or talking [2,3]. When in close proximity, aerosol particles can also transmit the infection, initiating the viral replication process in a new host [4].
Previously, smaller-scale epidemics, such as those in 2002, 2003, and 2012, were linked to the coronavirus, impacting various hosts, including bats, rodents, mammals, and birds [5,6], and by 2024, just over four years after the initial case in Wuhan, SARS-CoV-2 infections have affected nearly 704 million people, with approximately 7 million resulting in fatalities [7]. It is worth noting the well-known challenges in official reporting, which may lead to underestimating these values [8].
SARS-CoV-2 comprises three structural proteins and 16 non-structural proteins, among which the Mpro, or 3CLpro, holds a central role in the virus’s replication mechanism, acting as a cysteine protease [9,10]—positioned strategically at the active site of Mpro, a cysteine acts like molecular scissors, facilitating the cleavage of viral polyproteins. These lengthy protein chains, which are synthesized from the viral genome, are broken down into smaller, functionally independent units crucial for viral replication, assembly, and release. Given the vital importance of Mpro, various research groups are currently actively engaged in the development of targeted inhibitors, aiming for potential therapeutic applications for COVID-19 [11]. Several promising candidates are in the advanced stages of clinical trials, raising optimism about their eventual incorporation into the antiviral arsenal against this persistent pathogen [12].
In this context, Alugubelli and colleagues explored the potential of Boceprevir, which is a known inhibitor of HCV’s NSP3, as a repurposed treatment for COVID-19 [13]. This compound exhibited promising inhibitory effects on the Mpro of SARS-CoV-2, featuring an α -ketoamide group as a warhead, along with specific structural elements, including a P1 β -cyclobutylalanine segment, a P2 dimethyl cyclopropyl proline motif, a P3 test-butyl glycine unit, and a P4 N-terminal start-butyl carbamide. Through modifications at all four positions, twenty Mpro inhibitors based on Boceprevir were synthesized [13]. One such derivative, Nirmatrelvir (PF-07321332), underwent characterization for its potency against Mpro, both in vitro and in cellulo, using test tubes and 293T cells, respectively. The study involved analyzing the crystal structures of Mpro bound to ten of these inhibitors, alongside assessing the cytotoxicity and antiviral efficacy of four selected inhibitors.
By substituting the P1 site with a β -(S-2-oxopyrrolidin-3-yl)-alanil (Opal) residue and altering the warhead to an aldehyde, a high in vitro potency was achieved [13]. Additionally, the original structural components in P2, P3, and the P4 N-terminal positions of Boceprevir demonstrated superior performance compared with alternative chemical groups tested for their in vitro potency. In crystal structures, all inhibitors formed a covalent adduct, with the cysteine at the active site of Mpro [13]. Specifically, the Opal P1 residue, P2 dimethylcyclopropylproline, and P4 N-terminal tert-butylcarbamide established robust hydrophobic interactions with Mpro, elucidating the high in vitro potency associated with compounds containing these structural features.
Andi and colleagues conducted a study that focused on developing novel antivirals and vaccines that target Mpro [14]. Their findings, which are supported by crystallographic evidence and binding assay data, reveal that three drugs approved for hepatitis C virus treatment, along with two drug-like compounds, establish covalent bonds with the catalytic residue Cys145 within the active site of Mpro. Furthermore, molecular docking studies offered additional insights, aiding in the formulation of new antiviral inhibitors for SARS-CoV-2 by leveraging these drugs as lead compounds. It is worth mentioning that numerous other investigations into potential Mpro inhibitors have been documented in the literature [15,16,17,18].
The substrate binding site of Mpro, which is a region known for its remarkable conservation across various coronavirus genera, was investigated by Li and colleagues [19]. They specifically examined Nirmatrelvir, which is an oral inhibitor developed by Pfizer that targets the Mpro of SARS-CoV-2, with demonstrated effectiveness against other coronaviruses as well. Through their analysis of these structures, the researchers unveiled a conserved binding site shared between coronaviruses, offering valuable insights into the mechanism underlying the inhibition of viral replication [19].
Viewing it from this angle, Nirmatrelvir specifically targets Mpro [20]. By inhibiting this protease at a critical stage of proteolysis before viral RNA replication, Nirmatrelvir effectively curbs the proliferation of SARS-CoV-2. It has exhibited potent antiviral activity against various human coronaviruses and notable off-target selectivity [20]. Additionally, animal experiments confirmed its favorable safety profile and demonstrated its capacity to reduce the lung viral load in a mouse model engineered to simulate SARS-CoV-2 infection [20]. Preclinical investigations found no evidence of mutagenic interactions with DNA [21]. Moreover, a phase I clinical trial involving healthy volunteers showed plasma concentrations surpassing the threshold of in vitro cellular antiviral potency [21].
Conversely, Ivermectin is a broad-spectrum antiparasitic medication that initially garnered attention due to its in vitro ability to inhibit SARS-CoV-2 replication, sparking interest in its potential as a COVID-19 treatment [22,23,24]. However, its off-label use for this purpose is cautiously advised. Available in both human and veterinary formulations, Ivermectin is widely used clinically to combat parasites, like Onchocerca volvulus (which causes river blindness) and Strongyloides stercoralis. Nonetheless, at elevated doses, it can lead to potentially severe adverse effects, such as nausea, vomiting, diarrhea, dizziness, and even fatalities. Large-scale clinical trials that investigated its efficacy against COVID-19 yielded mixed results. While some studies failed to demonstrate significant benefits in reducing the symptoms or duration, others suggest potential advantages [25,26]. Consequently, a consensus based on consistent, high-quality evidence regarding Ivermectin’s role in treating COVID-19 has yet to emerge in the literature.
As presented, although initially developed as an antiparasitic, Ivermectin gained notoriety for its supposed antiviral effects against COVID-19. However, robust scientific evidence is crucial to support its clinical use. Nirmatrelvir, which is supported by solid scientific evidence, shows antiviral promise with proven efficacy and a favorable safety profile. Several works are found in the literature involving Ivermectin and SARS-CoV-2, both experimental [27,28,29] and theoretical [30,31,32], although there are no theoretical investigative works that compared the inhibitory potential of Ivermectin against Mpro. The choice of Mpro as a promising therapeutic target against SARS-CoV-2 is based on comprehensive and consistent experimental evidence due to its importance in viral replication and the potential of Mpro inhibitors to block disease progression [33,34,35]. In view of this, we conducted an in silico study to investigate the potential of interaction of Ivermectin with Mpro, using Nirmatrelvir as an inhibitor reference. Employing molecular docking, we first examined how Ivermectin bound to Mpro, specifically focusing on its interaction with protomer A, denoted as Mpro(A). Subsequently, we conducted molecular dynamics simulations for both Ivermectin and Nirmatrelvir bound to Mpro(A), facilitating a comprehensive analysis of their interaction profiles at the molecular level.

2. Methodology

2.1. Electronic Properties

We investigated the electronic properties of Ivermectin and Nirmatrelvir using density functional theory (DFT) and frontier molecular orbital (FMO) analyses within the Gaussian 16 software suite, Gaussian, Inc., Wallingford, CT, USA [36]. Molecular geometries were optimized through the semi-empirical Parameterization Method 6 (PM6) coupled with a steepest descent algorithm [37]. Following optimization, single-point energy calculations utilizing restricted DFT [38] were conducted to determine the energies of the highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), and their energy gap [39]. The analysis of HOMO/LUMO orbitals is presented here to provide contextual information about the electronic and structural differences between the investigated molecules. For a comprehensive DFT analysis, Becke’s three-parameter Lee–Yang–Parr (B3LYP) exchange–correlation functional [40] and the def2-TZVP (valence triple-zeta polarization) basis set [41] were employed.

2.2. Molecular Docking

Ligands and proteins were prepared using AutoDock Tools 1.5.6 (ADT), The Scripps Research Institute, San Diego, CA, USA [42]. Molecular docking simulations were employed to examine the non-covalent binding interactions between the Mpro(A) [20] (PDB: 7VH8) complex and Ivermectin. The crystallographic structure from the PDB contained the complex of Mpro and the inhibitor Nirmatrelvir. Autodock Vina, which is integrated into PyRx, The Scripps Research Institute, San Diego, CA, USA [43], was integrated with the genetic algorithm (GA) set as the default configuration for the docking simulations. This method employs a random walk simulation to explore the conformational space of the ligand while the protein remains in a fixed structure. The characterization of interactions between Ivermectin and Mpro was facilitated using ChimeraX, Resource for Biocomputing, Visualization, and Informatics, San Francisco, CA, USA [44] and Discovery Studio 2016, Dassault Systèmes, San Diego, CA, USA [45] programs. The simulation box encompassing the entire protein had dimensions of 60 × 40 × 65 Å3, with a center point at (−1.74, 18.92, −26.05) Å. Ligand screening was performed with an accuracy of ±2 Å for ligand positions and ±0.01 kcal/mol for binding affinities.

2.3. Molecular Dynamics

Nanoscale Molecular Dynamics (NAMD) 2.9 package, Beckman Institute for Advanced Science and Technology at the University of Illinois at Urbana-Champaign, Urbana, IL, USA [46] with the CHARMM36 force field [47] was employed for simulations. A time step of 1.0 fs was used with particle mesh Ewald (PME) electrostatics under NPT conditions (constant pressure and temperature). Langevin dynamics maintained a temperature of 300 K and pressure of 1 atm. The system was solvated with TIP3P water and neutralized with Na+ and Cl− ions to reach 0.15 M NaCl, mimicking a physiological environment. Initially, the protein–ligand complex was restrained, followed by energy minimization and equilibration for 10 ns. A subsequent 100 ns production run allowed the system to evolve freely. The radius of gyration (Rg) and root-mean-square deviation (RMSD) and fluctuation (RMSF) analyses using Visual Molecular Dynamics (VMD) 1.9.1, Beckman Institute for Advanced Science and Technology at the University of Illinois at Urbana-Champaign, Urbana, IL, USA [48] were performed to evaluate the system behavior over time. Finally, re-docking was performed based on interactions derived from the initial docking poses.

3. Results

Initially, DFT calculations were performed to understand the behavior of the molecules when isolated, particularly concerning the calculated frontier molecular orbitals (HOMO and LUMO). Figure 1 displays, on the left side, the results for the HOMO and LUMO of the Ivermectin molecule, while, on the right side, it presents the same characteristics for Nirmatrelvir.
The figure analysis reveals favorable features for the interaction of Ivermectin with its target protein, indicating a potential enhancement of the bioactivity. The high HOMO energy of −6.08 eV, coupled with a low LUMO energy of −0.96 eV, resulted in a significant energy gap of 5.12 eV for Ivermectin. In the case of Nirmatrelvir, the HOMO energy was −6.75 eV, and the LUMO energy was −1.65 eV. Consequently, the energy gap was 5.10 eV. These wide gaps show efficient charge transfers of both Ivermectin and Nirmatrelvir, potentially facilitating stable interactions with the target protein. Moreover, the high energy gap implies a strong inhibition efficiency due to the minimum energy required for electron removal from the HOMO.
It was observed that ideal inhibitory molecules readily accept and donate electrons, exhibiting an electron-rich character associated with superior inhibition potential. Consistent with these principles, the previously discussed calculated LUMO energy for Ivermectin showed a good electron-accepting capacity, further supporting its potential as an effective inhibitor. Regarding the overall electronic characteristics of both Ivermectin and Nirmatrelvir, we found they were similar in both molecules. This trend implies that Ivermectin, like Nirmatrelvir, can potentially be used as an inhibitor ligand for Mpro(A), hindering the replication process of SARS-CoV-2.
In molecular docking, we prioritized ligand poses with more favorable binding affinity, as the lower ΔG free energy scores indicate, reflecting the most probable and stable binding modes [49]. The docking simulations revealed favorable binding affinities for Ivermectin, with the top-scoring pose presenting a ΔG of −9.0 kcal/mol. Table 1 represents the values of all binding affinities and the RMSD upper bounds and lower bounds.
From this perspective, Figure 2 details the ligand–Mpro(A) interactions for the best score pose of the molecular docking. At the top, we present Nirmatrelvir. This complex was obtained from Zhao and colleagues [20]. We observed that the molecule established interactions near the N-terminus of Mpro(A), encompassing methionines 49 and 165 (MET), histidine 41 (HIS), phenylalanine 140 (PHE), glycine 143 (GLY), threonine 190 (THR), cysteine 145 (CYS), glutamate 166 (GLU), glutamines 189 and 192 (GLN), and proline 168 (PRO). On the other hand, at the bottom of this figure, Ivermectin, in the reported best energy configuration, occupied a distinct pocket closer to the C-terminus of Mpro(A), interacting with alanine 154 (ALA), aspartic acid 197 (ASP), leucines 286 and 287 (LEU), methionine 176 (MET), and proline 168 (PRO). Remarkably, while Nirmatrelvir and Ivermectin were bound to separate pockets near the N- and C-termini, respectively, they converged toward the central region. The dominant interactions of Nirmatrelvir involved conventional H-bonds. In contrast, those of Ivermectin were dominated by alkyl and pi-alkyl interactions, showing that Nirmatrelvir interacted more strongly with Mpro(A).
Building upon the molecular docking results for Ivermectin and the Nirmatrelvir/Mpro(A) complex discussed earlier, MD simulations of these systems were conducted to explore the interaction dynamics of both molecules with Mpro(A) under ambient temperature and pressure. This aimed to verify their stability in the presence of external factors and their persistence as inhibitors throughout the simulation.
Figure 3 illustrates the temporal evolution of the RMSD (a), RMSF (b), and radius of gyration (Rg) (c) simultaneously for Mpro(A) in the Apo form (green curves), Nirmatrelvir (blue curves), and Ivermectin (red curves) over the 100 ns simulation time. The middle panel also presents the binding free energy ( Δ G ) (d), and the lower panel indicates the number of H-bonds for Nirmatrelvir (e) and Ivermectin (f). The average RMSD values for Mpro(A) in the Apo form was 2.15 Å; for Mpro(A) with Nirmatrelvir, it was 2.00 Å; and with Ivermectin, it was 1.87 Å. These values, hovering around 2.00 Å, suggest that the reference structures were highly similar to the initial structures, signifying that both ligands remained securely bound to the cavity throughout the simulation time. The RMSD curves exhibited notable similarity, characterized by low fluctuations, showing that Ivermectin has the potential to act as an effective ligand against Mpro(A). Moreover, in the interval between 25 ns and 45 ns of the simulation, our results indicate that Nirmatrelvir experienced more pronounced RMSD fluctuations than Ivermectin, followed by increased stability beyond this point. These results imply a more pronounced stability of Mpro(A) in complex with Ivermectin. On the other hand, RMSF analysis was employed to assess the dynamic behavior of the Mpro(A) protein when bound to Nirmatrelvir and Ivermectin. Both complexes exhibited minimal amino acid fluctuations throughout the 100 ns simulations when compared with Mpro(A) in the Apo form, as illustrated by their RMSF plots. Thus, both ligands significantly stabilized the Mpro(A) protein. The Rg measured for the two systems (protein with Ivermectin and protein with Nirmatrelvir) varied between 22.01 Å and 22.66 Å, respectively. This result indicates that the protein atoms were distributed in a similar way along its axis during the simulation in the presence of both ligands. The results indicate that the ΔG for the binding of Ivermectin and Nirmatrelvir with the target protein had negative average values of −8.25 kcal/mol and −8.50 kcal/mol, respectively. These values indicate that the interaction between the ligands and the protein was spontaneous and favorable, suggesting a high affinity between them. It is worth noting that the average affinity value for Ivermectin (−8.25 kcal/mol) was close to the molecular docking value of 9.0 kcal/mol, which corroborated the strength of the interaction. Nirmatrelvir, on the other hand, had a slightly higher average affinity (−8.50 kcal/mol), indicating an even more stable bond with the protein.
In terms of the H-bond numbers, Figure 3e,f shows that the Mpro(A) complex with Ivermectin exhibited superior stabilization over 100 ns compared with the complex with Nirmatrelvir. The temporal progression of hydrogen bond formation between the protein and ligands was also scrutinized to glean insights into the molecular interactions. The analysis uncovered that Nirmatrelvir initially formed four hydrogen bonds with the Mpro(A) protein, with two and three persisting and alternating throughout the simulation. In contrast, Ivermectin formed two to three initial hydrogen bonds, with only one remaining stable throughout the simulation.
From the above-discussed MD simulations, we emphasize the snapshots at 0 ns, 50 ns, and 100 ns (see Figure 4). These visualizations unveil a notable difference in conformational flexibility, wherein both ligands maintained binding to Mpro(A) throughout the simulations. However, Nirmatrelvir exhibited a high conformational change compared with Ivermectin. Notably, the flexibility of the loop regions in the Mpro(A) protein played a critical role in retaining the ligands, like Ivermectin and Nirmatrelvir. These flexible loops can adapt their shapes, forming tight interactions with the ligands and creating a more favorable binding pocket. This flexibility allows for the formation of additional hydrogen bonds and other stabilizing forces, which is crucial for strong ligand retention.
Finally, Figure 5 delineates the interaction profiles of Nirmatrelvir and Ivermectin with Mpro(A) at 0 ns, 50 ns, and 100 ns, which represent the molecular dynamics snapshots discussed earlier. At 0 ns, Nirmatrelvir (a) interacted with the following amino acids: methionines 49 and 165 (MET); histidines 41, 163, and 164 (HIS); glutamate 166 (GLU); serine 144 (SER); threonine 189 (THR); glutamine 192 (GLN); and leucine 141 (LEU). At the 50 ns mark, Nirmatrelvir (c) engaged with amino acids, including methionines 49 and 165 (MET); histidines 41, 163, and 164 (HIS); phenylalanine 140 (PHE); cysteine 145 (CYS); glutamate 166 (GLU); and asparagine 142 (ASN). In the subsequent 100 ns interval (e), the Nirmatrelvir/Mpro(A) complex maintained interactions with methionines 49 and 165 (MET), histidine 41, glutamate 166 (GLU), asparagine 142 (ASN), glutamine 192 (GLN), and leucine 167 (LEU). Conversely, Ivermectin’s interactions persisted at 0 ns (b), 50 ns (d), and 100 ns (f) with the amino acids tyrosine 237 (THR) and leucine 287 (LEU). Notably, Ivermectin exhibited a more enduring interaction pattern with the same amino acids over time, while Nirmatrelvir interacted with approximately half of the amino acids at 0 ns and 50 ns compared with 100 ns. Consistent with the docking results, Nirmatrelvir predominantly relied on conventional hydrogen bonds, whereas Ivermectin primarily engaged in alkyl and pi-alkyl interactions. Consequently, the more robust nature of hydrogen bonds contributed to Nirmatrelvir forming a higher number of interactions with Mpro(A).

4. Summary and Conclusions

In conclusion, this study aimed to investigate the distinct binding modes exhibited by Nirmatrelvir and Ivermectin with the SARS-CoV-2 Mpro(A) from the perspective of understanding the nature of these molecules as inhibitors in the viral replication process. In this regard, Nirmatrelvir demonstrated engagement with the Mpro(A), forming numerous interactions involving amino acids, such as methionines 49 and 165 (MET); histidines 41, 163, and 164 (HIS); phenylalanine 140 (PHE); cysteine 145 (CYS); glutamate 166 (GLU); and asparagine 142 (ASN). Through molecular dynamics simulations, the RMSD analysis revealed a stable binding of Nirmatrelvir, with an RMSD value of around 2.0 Å, indicating its secure attachment to the binding pocket over the 100 ns simulation period.
On the other hand, Ivermectin exhibited good stability, maintaining prolonged interactions with specific amino acids, such as tyrosines 237 and 239 (THR) and leucines 271, 272, and 287 (LEU). The RMSD for the Ivermectin-Mpro(A) complex was 1.87 Å, reflecting its consistent and secure binding throughout the simulation. It is worth noting that Ivermectin interacted with approximately the same set of amino acids at both 50 ns and 100 ns, demonstrating enduring interaction patterns over time.
This diversified exploration of ligand–protein interactions emphasized the intricate nature of these binding profiles. Nirmatrelvir, which relied on conventional hydrogen bonds, engaged in a significant number of stronger interactions with Mpro(A), contributing to its strong binding. However, Ivermectin primarily utilized alkyl and pi-alkyl interactions, indicating a different binding mode.

Author Contributions

Y.A.d.O.S., K.S.B., M.L.P.J. and F.L.L.M.: data curation, formal analysis, methodology, prepared figures, and writing—original draft preparation. J.T.S., R.G., U.L.F. and L.A.R.J.: conceptualization, funding acquisition, and writing—reviewing and editing. All authors read and agreed to the published version of this manuscript.

Funding

The financial support from the Brazilian Research Council FAP-DF grants 00193.00001808/2022-71, 00193-00001247/2021-20, 00193-00000857/2021-14, 00193-00001247/2021-20, and 00193-00000811/2021-97, and the CNPq grants 302236/2018-0 and 350176/2022-1. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 88887.691997/2022-00. The financial support from FAP-DF grant 00193-00001807/2023-16.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors gratefully acknowledge the financial support from the Brazilian research agencies CNPq, CAPES, and FAP-DF. L.A.R.J. acknowledges the financial support from the Brazilian Research Council FAP-DF and the CNPq. M.L.P.J. acknowledges the financial support from FAP-DF. R.G., M.L.P.J., and L.A.R.J acknowledge CENAPAD-SP for providing computational facilities. F.L.L.M. and L.A.R.J. acknowledge the General Attorney of the Union—AGU grant 697.935/2019. L.A.R.J. acknowledges Núcleo de Computação de Alto Desempenho (NACAD) for providing the computational facilities through the Lobo Carneiro supercomputer. L.A.R.J. acknowledges the National Laboratory for Scientific Computing (LNCC/MCTI, Brazil) for providing HPC resources of the SDumont supercomputer (URL: http://sdumont.lncc.br (accessed on 7 November 2023)), which contributed to the research results reported in this paper. The authors acknowledge the technical and computational support from the LATITUDE Laboratory, University of Brasília; the TED 01/2019 grant from Advocacia Geral da União (grant AGU 697.935/2019); the Administrative Council for Economic Defense—CADE (TED 08700.000047/2019-14); the TED 01/2021 grant from the National Social Assistance Secretariat—SNAS/DGSUAS/CGRS; the Dean of Research and Innovation—DPI/UnB; the SISTER City Project—Safe and Real-Time Effective Intelligent Systems for Smart Cities (grant 625/2022); and the Foundation for Research Support of the Federal District—FAP-DF.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Atzrodt, C.L.; Maknojia, I.; McCarthy, R.D.; Oldfield, T.M.; Po, J.; Ta, K.T.; Stepp, H.E.; Clements, T.P. A Guide to COVID-19: A global pandemic caused by the novel coronavirus SARS-CoV-2. FEBS J. 2020, 287, 3633–3650. [Google Scholar] [CrossRef] [PubMed]
  2. Jayaweera, M.; Perera, H.; Gunawardana, B.; Manatunge, J. Transmission of COVID-19 virus by droplets and aerosols: A critical review on the unresolved dichotomy. Environ. Res. 2020, 188, 109819. [Google Scholar] [CrossRef] [PubMed]
  3. Louten, J. Virus transmission and epidemiology. In Essential Human Virology; Academic Press: Cambridge, MA, USA, 2016; pp. 71–92. [Google Scholar] [CrossRef]
  4. Asadi, S. Airborne Infectious Disease Transmission via Expiratory Aerosol Particles and Aerosolized Fomites; University of California: Davis, CA, USA, 2020. [Google Scholar]
  5. Weiss, R.A.; Sankaran, N. Emergence of epidemic diseases: Zoonoses and other origins. Fac. Rev. 2022, 11. [Google Scholar] [CrossRef] [PubMed]
  6. Rajeev, R.; Prathiviraj, R.; Kiran, G.S.; Selvin, J. Zoonotic evolution and implications of microbiome in viral transmission and infection. Virus Res. 2020, 290, 198175. [Google Scholar] [CrossRef] [PubMed]
  7. Worldometer. COVID-19 Coronavirus Pandemic. 2024. Available online: https://www.worldometers.info/coronavirus (accessed on 5 March 2024).
  8. Pellis, L.; Scarabel, F.; Stage, H.B.; Overton, C.E.; Chappell, L.H.; Fearon, E.; Bennett, E.; Lythgoe, K.A.; House, T.A.; Hall, I.; et al. Challenges in control of COVID-19: Short doubling time and long delay to effect of interventions. Philos. Trans. R. Soc. B 2021, 376, 20200264. [Google Scholar] [CrossRef] [PubMed]
  9. Cannalire, R.; Cerchia, C.; Beccari, A.R.; Di Leva, F.S.; Summa, V. Targeting SARS-CoV-2 proteases and polymerase for COVID-19 treatment: State of the art and future opportunities. J. Med. Chem. 2020, 65, 2716–2746. [Google Scholar] [CrossRef] [PubMed]
  10. Banerjee, R.; Perera, L.; Tillekeratne, L.V. Potential SARS-CoV-2 main protease inhibitors. Drug Discov. Today 2021, 26, 804–816. [Google Scholar] [CrossRef] [PubMed]
  11. Jin, Z.; Du, X.; Xu, Y.; Deng, Y.; Liu, M.; Zhao, Y.; Zhang, B.; Li, X.; Zhang, L.; Peng, C.; et al. Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors. Nature 2020, 582, 289–293. [Google Scholar] [CrossRef] [PubMed]
  12. Agost-Beltrán, L.; de la Hoz-Rodríguez, S.; Bou-Iserte, L.; Rodríguez, S.; Fernández-de-la Pradilla, A.; González, F.V. Advances in the development of SARS-CoV-2 Mpro inhibitors. Molecules 2022, 27, 2523. [Google Scholar] [CrossRef] [PubMed]
  13. Alugubelli, Y.R.; Geng, Z.Z.; Yang, K.S.; Shaabani, N.; Khatua, K.; Ma, X.R.; Vatansever, E.C.; Cho, C.C.; Ma, Y.; Xiao, J.; et al. A systematic exploration of boceprevir-based main protease inhibitors as SARS-CoV-2 antivirals. Eur. J. Med. Chem. 2022, 240, 114596. [Google Scholar] [CrossRef] [PubMed]
  14. Andi, B.; Kumaran, D.; Kreitler, D.F.; Soares, A.S.; Keereetaweep, J.; Jakoncic, J.; Lazo, E.O.; Shi, W.; Fuchs, M.R.; Sweet, R.M.; et al. Hepatitis C virus NS3/4A inhibitors and other drug-like compounds as covalent binders of SARS-CoV-2 main protease. Sci. Rep. 2022, 12, 12197. [Google Scholar] [CrossRef]
  15. Yamamoto, K.Z.; Yasuo, N.; Sekijima, M. Screening for inhibitors of main protease in SARS-Cov-2: In silico and in vitro approach avoiding peptidyl secondary amides. J. Chem. Inf. Model. 2022, 62, 350–358. [Google Scholar] [CrossRef] [PubMed]
  16. Lin, M.; Zeng, X.; Duan, Y.; Yang, Z.; Ma, Y.; Yang, H.; Yang, X.; Liu, X. Molecular mechanism of ensitrelvir inhibiting SARS-CoV-2 main protease and its variants. Commun. Biol. 2023, 6, 694. [Google Scholar] [CrossRef] [PubMed]
  17. Di Chio, C.; Previti, S.; Amendola, G.; Ravichandran, R.; Wagner, A.; Cosconati, S.; Hellmich, U.A.; Schirmeister, T.; Zappalà, M.; Ettari, R. Development of novel dipeptide nitriles as inhibitors of rhodesain of Trypanosoma brucei rhodesiense. Eur. J. Med. Chem. 2022, 236, 114328. [Google Scholar] [CrossRef] [PubMed]
  18. Miura, T.; Malla, T.R.; Owen, C.D.; Tumber, A.; Brewitz, L.; McDonough, M.A.; Salah, E.; Terasaka, N.; Katoh, T.; Lukacik, P.; et al. In vitro selection of macrocyclic peptide inhibitors containing cyclic γ2, 4-amino acids targeting the SARS-CoV-2 main protease. Nat. Chem. 2023, 15, 998–1005. [Google Scholar] [CrossRef] [PubMed]
  19. Tyndall, J.D. S-217622, a 3CL protease inhibitor and clinical candidate for SARS-CoV-2. J. Med. Chem. 2022, 65, 6496–6498. [Google Scholar] [CrossRef] [PubMed]
  20. Zhao, Y.; Fang, C.; Zhang, Q.; Zhang, R.; Zhao, X.; Duan, Y.; Wang, H.; Zhu, Y.; Feng, L.; Zhao, J.; et al. Crystal structure of SARS-CoV-2 main protease in complex with protease inhibitor PF-07321332. Protein Cell 2022, 13, 689–693. [Google Scholar] [CrossRef] [PubMed]
  21. Ahmad, B.; Batool, M.; Ain, Q.u.; Kim, M.S.; Choi, S. Exploring the binding mechanism of PF-07321332 SARS-CoV-2 protease inhibitor through molecular dynamics and binding free energy simulations. Int. J. Mol. Sci. 2021, 22, 9124. [Google Scholar] [CrossRef] [PubMed]
  22. Zaidi, A.K.; Dehgani-Mobaraki, P. The mechanisms of action of ivermectin against SARS-CoV-2—An extensive review. J. Antibiot. 2022, 75, 60–71. [Google Scholar] [CrossRef] [PubMed]
  23. Hill, A.; Garratt, A.; Levi, J.; Falconer, J.; Ellis, L.; McCann, K.; Pilkington, V.; Qavi, A.; Wang, J.; Wentzel, H. Retracted: Meta-analysis of randomized trials of ivermectin to treat SARS-CoV-2 infection. In Open Forum Infectious Diseases; Oxford University Press: New York, NY, USA, 2021; Volume 8, p. ofab358. [Google Scholar]
  24. Sharun, K.; Dhama, K.; Patel, S.K.; Pathak, M.; Tiwari, R.; Singh, B.R.; Sah, R.; Bonilla-Aldana, D.K.; Rodriguez-Morales, A.J.; Leblebicioglu, H. Ivermectin, a new candidate therapeutic against SARS-CoV-2/COVID-19. Ann. Clin. Microbiol. Antimicrob. 2020, 19, 23. [Google Scholar] [CrossRef]
  25. Cobos-Campos, R.; Apiñaniz, A.; Parraza, N.; Cordero, J.; García, S.; Orruño, E. Potential use of ivermectin for the treatment and prophylaxis of SARS-CoV-2 infection. Curr. Res. Transl. Med. 2021, 69, 103309. [Google Scholar] [CrossRef] [PubMed]
  26. Buonfrate, D.; Chesini, F.; Martini, D.; Roncaglioni, M.C.; Fernandez, M.L.O.; Alvisi, M.F.; De Simone, I.; Rulli, E.; Nobili, A.; Casalini, G.; et al. High-dose ivermectin for early treatment of COVID-19 (COVER study): A randomised, double-blind, multicentre, phase II, dose-finding, proof-of-concept clinical trial. Int. J. Antimicrob. Agents 2022, 59, 106516. [Google Scholar] [CrossRef] [PubMed]
  27. Arévalo, A.; Pagotto, R.; Pórfido, J.L.; Daghero, H.; Segovia, M.; Yamasaki, K.; Varela, B.; Hill, M.; Verdes, J.; Duhalde Vega, M.; et al. Ivermectin reduces in vivo coronavirus infection in a mouse experimental model. Sci. Rep. 2021, 11, 7132. [Google Scholar] [CrossRef] [PubMed]
  28. Peña-Silva, R.; Duffull, S.B.; Steer, A.C.; Jaramillo-Rincon, S.X.; Gwee, A.; Zhu, X. Pharmacokinetic considerations on the repurposing of ivermectin for treatment of COVID-19. Br. J. Clin. Pharmacol. 2021, 87, 1589. [Google Scholar] [CrossRef] [PubMed]
  29. Dixit, A.; Yadav, R.; Singh, A.V. Ivermectin: Potential role as repurposed drug for COVID-19. Malays. J. Med. Sci. MJMS 2020, 27, 154. [Google Scholar] [CrossRef] [PubMed]
  30. Choudhury, A.; Das, N.C.; Patra, R.; Bhattacharya, M.; Ghosh, P.; Patra, B.C.; Mukherjee, S. Exploring the binding efficacy of ivermectin against the key proteins of SARS-CoV-2 pathogenesis: An in silico approach. Future Virol. 2021, 16, 277–291. [Google Scholar] [CrossRef]
  31. Aminpour, M.; Cannariato, M.; Preto, J.; Safaeeardebili, M.E.; Moracchiato, A.; Doria, D.; Donato, F.; Zizzi, E.A.; Deriu, M.A.; Scheim, D.E.; et al. In silico analysis of the multi-targeted mode of action of ivermectin and related compounds. Computation 2022, 10, 51. [Google Scholar] [CrossRef]
  32. Kaur, H.; Shekhar, N.; Sharma, S.; Sarma, P.; Prakash, A.; Medhi, B. Ivermectin as a potential drug for treatment of COVID-19: An in-sync review with clinical and computational attributes. Pharmacol. Rep. 2021, 73, 736–749. [Google Scholar] [CrossRef] [PubMed]
  33. Hu, Q.; Xiong, Y.; Zhu, G.H.; Zhang, Y.N.; Zhang, Y.W.; Huang, P.; Ge, G.B. The SARS-CoV-2 main protease (Mpro): Structure, function, and emerging therapies for COVID-19. MedComm 2022, 3, e151. [Google Scholar] [CrossRef] [PubMed]
  34. Liang, S.; Liu, X.; Zhang, S.; Li, M.; Zhang, Q.; Chen, J. Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations. Phys. Chem. Chem. Phys. 2022, 24, 1743–1759. [Google Scholar] [CrossRef]
  35. Katre, S.G.; Asnani, A.J.; Pratyush, K.; Sakharkar, N.G.; Bhope, A.G.; Sawarkar, K.T.; Nimbekar, V.S. Review on development of potential inhibitors of SARS-CoV-2 main protease (MPro). Future J. Pharm. Sci. 2022, 8, 36. [Google Scholar] [CrossRef] [PubMed]
  36. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16 Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  37. Stewart, J.J. Application of the PM6 method to modeling proteins. J. Mol. Model. 2009, 15, 765–805. [Google Scholar] [CrossRef] [PubMed]
  38. Parr, R.G. Density functional theory. Annu. Rev. Phys. Chem. 1983, 34, 631–656. [Google Scholar] [CrossRef]
  39. Sham, L.J.; Schlüter, M. Density-functional theory of the energy gap. Phys. Rev. Lett. 1983, 51, 1888. [Google Scholar] [CrossRef]
  40. Becke, A.D. A new mixing of Hartree–Fock and local density-functional theories. J. Chem. Phys. 1993, 98, 1372–1377. [Google Scholar] [CrossRef]
  41. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  42. Huey, R.; Morris, G.M.; Forli, S. Using AutoDock 4 and AutoDock vina with AutoDockTools: A tutorial. Scripps Res. Inst. Mol. Graph. Lab. 2012, 10550, 1000. [Google Scholar]
  43. Dallakyan, S.; Olson, A.J. Small-molecule library screening by docking with PyRx. In Chemical Biology: Methods and Protocols; Humana Press: New York, NY, USA, 2015; pp. 243–250. [Google Scholar]
  44. Meng, E.C.; Goddard, T.D.; Pettersen, E.F.; Couch, G.S.; Pearson, Z.J.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Tools for structure building and analysis. Protein Sci. 2023, 32, 4792. [Google Scholar] [CrossRef] [PubMed]
  45. Studio, Discovery. Discovery studio. Accelrys [2.1], 2016; Volume 420. Available online: https://www.3ds.com/products/biovia/discovery-studio (accessed on 24 April 2024).
  46. Phillips, J.C.; Zheng, G.; Kumar, S.; Kalé, L.V. NAMD: Biomolecular simulation on thousands of processors. In Proceedings of the SC’02: Proceedings of the 2002 ACM/IEEE Conference on Supercomputing, Baltimore, MD, USA, 16–22 November 2002; pp. 36–36. [Google Scholar]
  47. Klauda, J.B.; Venable, R.M.; Freites, J.A.; O’Connor, J.W.; Tobias, D.J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A.D., Jr.; Pastor, R.W. Update of the CHARMM all-atom additive force field for lipids: Validation on six lipid types. J. Phys. Chem. B 2010, 114, 7830–7843. [Google Scholar] [CrossRef]
  48. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef] [PubMed]
  49. Fukunishi, Y.; Yamashita, Y.; Mashimo, T.; Nakamura, H. Prediction of protein- compound binding energies from known activity data: Docking-score-based method and its applications. Mol. Inform. 2018, 37, 1700120. [Google Scholar] [CrossRef]
Biomolecules 14 00755 g001
Figure 1. Analysis of the frontier molecular orbitals (FMOs) of Ivermectin and Nirmatrelvir, depicting their LUMOs and HOMOs, and the energy GAPs (eV) between them.
Figure 1. Analysis of the frontier molecular orbitals (FMOs) of Ivermectin and Nirmatrelvir, depicting their LUMOs and HOMOs, and the energy GAPs (eV) between them.
Biomolecules 14 00755 g001
Biomolecules 14 00755 g002
Figure 2. Crystal structure of Nirmatrelvir and Mpro(A) from RCSB PDB Protein Data Bank (up) and structure of Ivermectin and Mpro(A) from molecular docking (down). The blue-to-red color band represents the orientation from the N-terminus to the C-terminus of Mpro(A).
Figure 2. Crystal structure of Nirmatrelvir and Mpro(A) from RCSB PDB Protein Data Bank (up) and structure of Ivermectin and Mpro(A) from molecular docking (down). The blue-to-red color band represents the orientation from the N-terminus to the C-terminus of Mpro(A).
Biomolecules 14 00755 g002
Biomolecules 14 00755 g003
Figure 3. Root-mean-square deviation (a), root-mean-square fluctuation (b), radius of gyration (c), binding free energy (d), and number of hydrogen bonds (e,f) of Mpro(A) in the Apo form and complex with Nirmatrelvir and Ivermectin throughout the 100 ns simulation.
Figure 3. Root-mean-square deviation (a), root-mean-square fluctuation (b), radius of gyration (c), binding free energy (d), and number of hydrogen bonds (e,f) of Mpro(A) in the Apo form and complex with Nirmatrelvir and Ivermectin throughout the 100 ns simulation.
Biomolecules 14 00755 g003
Biomolecules 14 00755 g004
Figure 4. By employing the optimal docking pose (see Figure 2), we examined snapshots at 0 ns (a,b), 50 ns (c,d), and 100 ns (e,f). These depictions reveal a significant disparity in the conformational flexibility: while both ligands consistently maintained their binding to Mpro(A) throughout the simulations, Nirmatrelvir exhibited a markedly more pronounced conformational variation compared with Ivermectin. The red, black, cyan, blue, and grey spheres represent the oxygen, carbon, fluorine, nitrogen, and hydrogen atoms, respectively.
Figure 4. By employing the optimal docking pose (see Figure 2), we examined snapshots at 0 ns (a,b), 50 ns (c,d), and 100 ns (e,f). These depictions reveal a significant disparity in the conformational flexibility: while both ligands consistently maintained their binding to Mpro(A) throughout the simulations, Nirmatrelvir exhibited a markedly more pronounced conformational variation compared with Ivermectin. The red, black, cyan, blue, and grey spheres represent the oxygen, carbon, fluorine, nitrogen, and hydrogen atoms, respectively.
Biomolecules 14 00755 g004
Biomolecules 14 00755 g005
Figure 5. Interactions of Nirmatrelvir and Ivermectin with Mpro(A) at 0 ns (a,b), 50 ns (c,d), and 100 ns (e,f). The red, black, cyan, and blue spheres represent the oxygen, carbon, fluorine, and nitrogen atoms, respectively.
Figure 5. Interactions of Nirmatrelvir and Ivermectin with Mpro(A) at 0 ns (a,b), 50 ns (c,d), and 100 ns (e,f). The red, black, cyan, and blue spheres represent the oxygen, carbon, fluorine, and nitrogen atoms, respectively.
Biomolecules 14 00755 g005
Table 1. AutoDock Vina binding affinity (ΔG) and RMSD upper bound/lower bound for Ivermectin/Mpro(A) complex.
Table 1. AutoDock Vina binding affinity (ΔG) and RMSD upper bound/lower bound for Ivermectin/Mpro(A) complex.
Ranking Score Δ G (kcal/mol)Distance from Best Mode (Å)
RMSD Upper BoundRMSD Lower Bound
1−9.00.00.0
2−8.710.8715.250
3−8.711.2025.632
4−8.715.5748.398
5−8.713.6332.868
6−8.418.25312.631
7−8.323.13815.835
8−8.323.12515.631
9−8.19.7875.067
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

de Oliveira Só, Y.A.; Bezerra, K.S.; Gargano, R.; Mendonça, F.L.L.; Souto, J.T.; Fulco, U.L.; Pereira Junior, M.L.; Junior, L.A.R. In Silico Comparative Analysis of Ivermectin and Nirmatrelvir Inhibitors Interacting with the SARS-CoV-2 Main Protease. Biomolecules 2024, 14, 755. https://doi.org/10.3390/biom14070755

AMA Style

de Oliveira Só YA, Bezerra KS, Gargano R, Mendonça FLL, Souto JT, Fulco UL, Pereira Junior ML, Junior LAR. In Silico Comparative Analysis of Ivermectin and Nirmatrelvir Inhibitors Interacting with the SARS-CoV-2 Main Protease. Biomolecules. 2024; 14(7):755. https://doi.org/10.3390/biom14070755

Chicago/Turabian Style

de Oliveira Só, Yuri Alves, Katyanna Sales Bezerra, Ricardo Gargano, Fabio L. L. Mendonça, Janeusa Trindade Souto, Umberto L. Fulco, Marcelo Lopes Pereira Junior, and Luiz Antônio Ribeiro Junior. 2024. "In Silico Comparative Analysis of Ivermectin and Nirmatrelvir Inhibitors Interacting with the SARS-CoV-2 Main Protease" Biomolecules 14, no. 7: 755. https://doi.org/10.3390/biom14070755

APA Style

de Oliveira Só, Y. A., Bezerra, K. S., Gargano, R., Mendonça, F. L. L., Souto, J. T., Fulco, U. L., Pereira Junior, M. L., & Junior, L. A. R. (2024). In Silico Comparative Analysis of Ivermectin and Nirmatrelvir Inhibitors Interacting with the SARS-CoV-2 Main Protease. Biomolecules, 14(7), 755. https://doi.org/10.3390/biom14070755

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop