Research Article |
Corresponding author: Ahmed L. Alaofi ( ahmedofi@ksu.edu.sa ) Academic editor: Georgi Momekov
© 2025 Ahmed L. Alaofi, Mudassar Shahid, Mohd Abul Kalam.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Alaofi AL, Shahid M, Kalam MA (2025) Inhibition of MERS-CoV papain-like protease by sunitinib: In vitro and in silico investigations. Pharmacia 72: 1-12. https://doi.org/10.3897/pharmacia.72.e141634
|
Middle East Respiratory Syndrome Coronavirus (MERS-CoV) remains a significant public health threat, with high mortality rates and no approved antiviral therapies. The papain-like protease (PLpro) of MERS-CoV plays a critical role in viral replication and immune evasion, making it a key target for drug discovery. This study evaluated the inhibitory effects of four anticancer drugs (sunitinib, olaparib, mitoxantrone, and bicalutamide) on recombinant MERS-CoV PLpro using a combination of in vitro and in silico techniques. Protease inhibition assays revealed that sunitinib displayed potent, dose-dependent inhibition of PLpro activity, with an IC50 of 1.75 µM, while olaparib, mitoxantrone, and bicalutamide exhibited negligible inhibition. Thermal shift assays confirmed the strong interaction of sunitinib with PLpro, showing a ΔTm of 26.64 °C, indicative of increased protein stability. Furthermore, molecular dynamics (MD) simulations and docking studies provided structural insights into the mechanism of inhibition. Sunitinib bound within the thumb domain of PLpro, forming stable interactions with residues such as D76, R82, and F79. Binding induced significant stabilization of PLpro’s structure, reducing flexibility in critical regions, including the thumb and catalytic domains, as indicated by a decreased radius of gyration and alterations in the free energy landscape. Importantly, the stabilization of PLpro by sunitinib was consistent between in vitro and in silico analyses, highlighting its robust inhibitory potential. These findings position sunitinib as a promising inhibitor of MERS-CoV PLpro, with strong binding affinity and the ability to disrupt enzymatic function. Further preclinical studies are warranted to explore its therapeutic potential against MERS-CoV. This study underscores the utility of repurposing existing drugs for emerging viral threats and contributes to the development of targeted antiviral strategies.
Middle East respiratory syndrome coronavirus, molecular dynamics simulations, viral inhibition
Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is an epidemic and zoonotic virus that was first identified in Saudi Arabia in 2012 (
The virus is primarily transmitted from dromedary camels to humans, although human-to-human transmission, particularly in healthcare settings, is also possible and poses a significant risk (
Despite its significant public health impact, there are currently no approved vaccines or specific antiviral treatments for MERS-CoV. Treatment options are primarily supportive, focusing on alleviating symptoms and providing respiratory support in severe cases (
In targeting viral infections, protease inhibitors have proven to be a vital therapeutic class, particularly for viruses like HIV, Hepatitis C (HCV), and SARS-CoV-2. For instance, HIV protease inhibitors such as Ritonavir and Lopinavir disrupt the viral replication cycle by inhibiting the HIV protease enzyme (
In literature, repurposing of small molecules as inhibitors for coronavirus enzymes is extensively studied. Sunitinib is a multi-targeted receptor tyrosine kinase (RTK) inhibitor that is primarily used in the treatment of renal cell carcinoma and gastrointestinal stromal tumors. It works by inhibiting multiple signaling pathways that promote tumor growth, angiogenesis, and metastasis (
These drugs highlight the diversity in therapeutic strategies employed in oncology, ranging from targeted kinase inhibition and DNA repair disruption to androgen receptor antagonism and DNA intercalation. Each of these drugs plays a critical role in the tailored treatment of various cancers, reflecting the complexity and specificity required in modern cancer therapies. These drugs were shown to have promising antiviral activity against SARS-CoV-2 and have varying effects on inhibiting SARS-CoV-2 replication, but no study has been conducted on MERS-CoV. Therefore, we conducted experimental and in silico studies to evaluate the inhibitory activities of sunitinib, olaparib, mitoxantrone, and bicalutamide against MERS-CoV PLpro.
Sunitinib, olaparib, mitoxantrone, and bicalutamide were purchased from Biosynth, Carbosynth Ltd. (United Kingdom). Z-Arg-Leu-Arg-Gly-Gly-AMC Acetate was purchased from Bachem (CA, United States). SYPROTM orange protein gel stain was purchased from Invitrogen by Thermo Fisher Scientific (Rockford, Illinois, United States). All other chemicals were analytical grade and were obtained from standard commercial suppliers.
E. coli BL21(DE3)pLysS harboring MERS-CoV-Plpro (Genscript) was used to express Plpro as described in our previous work (
MERS-CoV Plpro inhibition and IC50 values for anti-cancer drugs (sunitinib, olaparib, mitoxantrone, and bicalutamide) were determined using a fluorescence-based assay in a 384-well plate format. The hydrolysis of the fluorogenic peptide Z-RLRGG-AMC leads to increased fluorescence of the AMC moiety, allowing precise measurement of enzyme activity. The reactions were carried out in 50 µL total volume, consisting of 20 mM Tris buffer (pH 8.0), 30 µM Z-RLRGG-AMC, and various concentrations of inhibitors (0.78–100 µM). After the addition of 60 nM Plpro, the plates were shaken for 30 seconds, and fluorescence was monitored using a BiotekHT Microplate Reader at excitation 360 nm and emission 460 nm. IC50 values were calculated using the dose-response inhibition function in GraphPad Prism. The experiments were performed in triplicate for accuracy according to our previous methods (
To assess the interaction between anti-cancer drugs (sunitinib, olaparib, mitoxantrone, and bicalutamide) and the MERS-CoV PLpro protein, Differential Scanning Fluorimetry (DSF) was performed using the Applied Biosystems 7500 Real-Time PCR System. The assay involved incubating varying concentrations of tested drugs, ranging from 0.78 to 100 µM, with a reaction buffer containing a final concentration of 2 µM MERS-CoV PLpro. The buffer without the test compound acted as the control. After the mixtures were incubated for 30 minutes at 25 °C, the fluorescence signal was recorded using SYPRO Orange dye, which binds to hydrophobic regions exposed during protein unfolding. The fluorescence was monitored as the temperature was increased incrementally at a rate of 1 °C per minute, from 25 °C to 95 °C. This method allowed for the determination of the protein’s melting temperature (Tm), which is an indicator of protein stability (
The structure of PLpro of MERS-CoV was obtained from the Protein Data Bank (PDB ID 4RNA), while sunitinib was abstracted from a ligand-protein complex (PDB ID 3MIY). The small molecule, sunitinib, was prepared as follows: The Avogadro website was used to assign all hydrogen atoms, and the output file (.mol2) of the sunitinib ligand was obtained (
Molecular dynamics (MD) simulations were used to assess the binding interactions of sunitinib toward MERS-CoV PLpro. Sunitinib was placed with PLpro without predefined binding sites, i.e., before docking to thoroughly sample the PLpro’s conformational space, increasing the likelihood of identifying productive binding poses (
Utilizing simulated structures takes into account the dynamics and flexible trajectories of both ligand and protein. Therefore, the trajectories at 100 ns produced from MD simulations of MERS-CoV PLpro and sunitinib (i.e., PDB files) were used to perform blind docking using the HADDOCK.2.4 server (
To assess the in vitro inhibition potency of sunitinib, olaparib, mitoxantrone, and bicalutamide (Fig.
Chemical structures of sunitinib, olaparib, mitoxantrone, and bicalutamide (a). MERS protease inhibition by sunitinib, olaparib, mitoxantrone, and bicalutamide (b). Dose response curve B. Percent inhibition, dose response curve of MERS-CoV Plpro inhibition identified from the IC50 protease inhibition assay with dissociation of RLRGG-AMC substrate. Values represent the average ±standard deviation of three replicates.
Non-linear regression analysis employing a simple inhibition model allowed estimating apparent inhibition at different concentrations. Sunitinib can be assumed to bind to an active site on MERS-CoV Plpro, functioning as a competitive inhibitor. This apparent inhibition can be observed in the IC50 value of 1.75 μM, while there is no significant inhibition for olaparib, mitoxantrone, and bicalutamide (Table
While TSA is a widely accepted method for evaluating protein-ligand interactions by measuring shifts in protein melting temperature, it has limitations, including the potential for false negatives and not always correlating directly with inhibitory activity. However, methods measuring actual enzymatic activity offer more reliable evidence of inhibition. In our assay, olaparib (ΔTm 1.4), mitoxantrone (ΔTm 5.45), and bicalutamide (ΔTm 0.94) exhibited no binding affinity with MERS-CoV PLpro, and their presence did not significantly affect the melting temperature of the protein, indicating no major interaction (Table
Melting points (in Celsius; °C) and melting point differences ΔTm (ΔTm = Tmmolecule -TmPLpro) calculated for small molecules incubated with MERS-CoV PLpro.
Sunitinib | Olaparib | Mitoxantrone | Bicalutamide | MERS-CoV PLpro | |
---|---|---|---|---|---|
Tm (°C) | 69.6 | 41.56 | 48.41 | 42.02 | 42.96 |
ΔTm (°C) | 26.64 | 1.4 | 5.45 | 0.94 | 0.00 |
Sunitinib, however, demonstrated consistent inhibition across a range of concentrations (0.78–100 µM), with a corresponding thermal shift in the TSA (ΔTm ranging from 16.26 °C to 23.4 °C) (Fig.
The melting temperature (Tm) of MERS-CoV-2 PLpro by thermal shift assay (TSA). Representative image of Tm profile of sunitinib (a), olaparib (b), mitoxantrone (c), bicalutamide (d), and MERS-CoV-2 PLpro (e). The raw and fitted data figure was plotted by TSA-CRAFT (Suppl. material
Cα-root mean square deviation (Cα-RMSD) of free and bound PLpro reached a plateau, indicating stable equilibrium systems (Fig.
The Cα-RMSD as a function of time (a) was plotted for free PLpro (black color) and bound PLpro (red color). Both systems reached a plateau indicated by the stabilities of the simulations. In (b), the Cα-RMSF as a function of residues was plotted for free PLpro (black color) and bound PLpro (red color). The thumb domain (A61-L178) showed a significant rigidity for bound PLpro in comparison to the free one. In (c), the left panel showed a surface representation of free PLpro (smudge color) and bound PLpro (violet color). In the right panel, the boxes showed cartoon representations of the alignment of free and bound PLpro. The upper box was rotated 180 °C in the lower box. Dark blue arrows showed the significant differences between bound and free PLpro.
The use of radius of gyration (Rg) is an indicator of the conformational compactness of a protein. We assessed the conformational compactness using the Rg of the backbone over time for the free and sunitinib-bound PLpro. Sunitinib binding to MERS-CoV PLpro decreased the conformational compactness of the PLpro structure (Fig.
The docking techniques are commonly used to identify binding sites for ligands in receptor pockets. Our docking results suggested sunitinib binding to the thumb domain of MERS-CoV PLpro (Fig.
The left panel shows the docking results of sunitinib (stick representation; indicated by a small arrow) to the MERS-CoV PLpro (cartoon representation). MERS-CoV PLpro domains were labeled next to the corresponding domain. The right panel shows a surface representation of the binding pocket of sunitinib (indicated by a small arrow) in the upper box. The lower box shows the interaction residues of PLpro with the sunitinib molecule.
The COVID-19 pandemic has underscored the urgency of developing effective pharmacological interventions against other coronaviruses, particularly Middle East Respiratory Syndrome Coronavirus (MERS-CoV), which poses a significant threat due to its high mortality rate. The primary challenge in developing treatments for MERS-CoV is the virus’s potential for rapid mutation and drug resistance. To address this, there is a need for a repository of active compounds with diverse mechanisms of action that maintain efficacy while minimizing susceptibility to resistance (Zumla et al. 2015;
The screening process involved assessing the efficacy of four drugs: sunitinib, olaparib, mitoxantrone, and bicalutamide against MERS-CoV PLpro. Using a peptidic FRET substrate (Z-RLRGG-AMC), we monitored the enzyme’s activity in a continuous enzymatic assay. Following initial evaluations of each drug’s ability to inhibit PLpro in a dose-dependent manner, we further examined their effects on the protease’s thermal stability using a fluorescence-based Thermal Shift Assay (TSA). Olaparib, mitoxantrone, and bicalutamide showed insignificant inhibition of MERS-CoV PLpro at any concentration. However, sunitinib demonstrated consistent inhibition across a range of concentrations (0.78–100 µM), with IC50 1.75 µM, indicating a significant inhibition of MERS-CoV PLpro. This shift in the protein,s unfolding temperature due to the drugs, precise and effective binding with the protease suggests their potential as candidates for further development and preclinical trials.
It is interesting that in silico studies were consistent with in vitro studies since sunitinib increased the MERS-CoV PLpro stability observed by elevated melting temperature and stable structural conformation of bound PLpro (i.e., sunitinib-PLpro complex). The interaction between active compounds or ligands and a protein can lead to changes in the protein’s stability, particularly in response to thermal or chemical denaturation. When a ligand preferentially binds to a specific conformational state of the protein, this interaction either stabilizes or destabilizes the protein depending on whether the native or non-native state is favored (
There are limited studies on the role of olaparib’s potential application against viral infections, including viruses like MERS-CoV or SARS-CoV-2, which may also counteract SARS‐CoV‐2‐induced and inflammation‐induced cell death and support cell survival (
The study demonstrates a strong correlation between in silico and in vitro findings regarding sunitinib’s role as a PLpro inhibitor for MERS-CoV. In vitro assays, including protease inhibition and thermal shift analysis, revealed a dose-dependent inhibitory effect of sunitinib, with a substantial thermal stability increase (ΔTm of 26.64 °C) indicating a strong interaction with PLpro. Complementing this, in silico molecular docking and dynamic simulations identified sunitinib binding at the thumb domain of PLpro, stabilizing the protein structure by reducing flexibility in key regions such as the M137-D142 segment. This is consistent with previous findings where sunitinib has shown inhibitory interactions with key viral targets, including HCV NS5B polymerase virus envelope protein and spike protein, by stabilizing regions and disrupting their functional activity. These findings highlight sunitinib’s versatility in targeting conserved viral domains, supporting its repurposing for antiviral therapy (
Computational tools are extensively used in the literature for studying proteins and macromolecules. Our MD simulations and docking results were consistent with experimental data as sunitinib was able to inhibit the catalytic activity of the MERS-CoV PLpro enzyme. The MD simulation results indicate a significant stabilization of the PLpro structure in the complex form compared to free PLpro. Based on our observation, the overall flexibility of the PLpro was impacted by the presence of sunitinib, indicating a specific and potent binding of sunitinib to MERS-CoV PLpro. There was a clear stable structure of bound Plpro, suggesting favorable binding to the sunitinib molecule. This was consistent with the binding site of sunitinib, as the thumb domain showed higher rigidity than the palm and finger domains of PLpro. Computational analysis predicts a strong binding affinity between sunitinib and the thumb domain of MERS-CoV PLpro, mediated by both hydrophobic and hydrogen bond interactions. Experimental validation of this binding affinity would strengthen these findings.
MERS-CoV PLpro is a significant drug target with various inhibitors discovered through different approaches. However, quality science requires thorough validation, including detailed biochemical and binding studies. Drug candidates should show consistent IC50 values in FRET assays and align with the results of thermal shift assay results. This research aims to enhance understanding of PLpro inhibitors while ensuring scientific rigor. It supports further drug development without compromising quality in the rush for effective treatments. The identification of sunitinib as a potential inhibitor of MERS-CoV PLpro through enzymatic, thermal stability assays, and MD simulations underscores the importance of repurposing existing drugs and offers a promising avenue for therapeutic development. This approach not only addresses current treatment gaps but also enhances global preparedness for future coronavirus outbreaks.
Conflict of interest
The authors have declared that no competing interests exist.
Ethical statements
The authors declared that no clinical trials were used in the present study.
The authors declared that no experiments on humans or human tissues were performed for the present study.
The authors declared that no informed consent was obtained from the humans, donors or donors’ representatives participating in the study.
The authors declared that no experiments on animals were performed for the present study.
The authors declared that no commercially available immortalised human and animal cell lines were used in the present study.
Funding
This study was funded by the Research Supporting Project number (RSPD2025R560), King Saud University, Riyadh, Saudi Arabia.
Author contributions
Ahmed L. Alaofi: Conceptualization, Funding acquisition, Data curation, Writing – original draft, Writing – review & editing, Investigation, Validation, Formal analysis, Methodology, Supervision, Resources, Project administration, Software. Mudassar Shahid: Conceptualization, Data curation, Writing – original draft, Writing – review & editing, Visualization, Investigation, Validation, Formal analysis, Methodology, Software. Mohd Abul Kalam: Writing – review & editing, Methodology.
Author ORCIDs
Ahmed L. Alaof https://orcid.org/0000-0001-8967-173X
Mudassar Shahid https://orcid.org/0000-0003-3714-4772
Mohd Abul Kalam https://orcid.org/0000-0002-5713-8858
Data availability
All of the data that support the findings of this study are available in the main text or Supplementary Information.
Supplementary information
Data type: docx