Skip to main content

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 29 May 2023
Sec. Virus and Host
This article is part of the Research Topic Evolution, Pathogenesis, Host Interactions and Therapeutic Strategies Against Monkeypox Virus View all 10 articles

Computer-assisted drug repurposing for thymidylate kinase drug target in monkeypox virus

  • 1Department of Biochemistry, Computational Medicinal Chemistry Laboratory, Abdul Wali Khan University, Mardan, Pakistan
  • 2Center for Medical Genetics and Hunan Key Laboratory of Medical Genetics, School of Life Sciences, Central South University, Changsha, Hunan, China
  • 3Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Al-Quwayiyah, Shaqra University, Riyadh, Saudi Arabia
  • 4Department of Life Sciences, School of Science, University of Management and Technology (UMT), Lahore, Pakistan
  • 5Department of Pharmacology and Toxicology, College of Pharmacy, Umm Al-Qura University, Makkah, Saudi Arabia
  • 6Institute of Biomedical Sciences, Shanxi University, Taiyuan, China
  • 7Department of Obstetrics and Gynecology, Nanfang Hospital, Southern Medical University, Guangzhou, China
  • 8Department of Central Laboratory, SSL Central Hospital of Dongguan City, Affiliated Dongguan Shilong People’s Hospital of Southern Medical University, Dongguan, China

Introduction: Monkeypox is a zoonotic disease caused by brick-shaped enveloped monkeypox (Mpox) virus that belongs to the family of ancient viruses known as Poxviridae. Subsequently, the viruses have been reported in various countries. The virus is transmitted by respiratory droplets, skin lesions, and infected body fluids. The infected patients experience fluid-filled blisters, maculopapular rash, myalgia, and fever. Due to the lack of effective drugs or vaccines, there is a need to identify the most potent and effective drugs to reduce the spread of monkeypox. The current study aimed to use computational methods to quickly identify potentially effective drugs against the Mpox virus.

Methods: In our study, the Mpox protein thymidylate kinase (A48R) was targeted because it is a unique drug target. We screened a library of 9000 FDA-approved compounds of the DrugBank database by using various in silico approaches, such as molecular docking and molecular dynamic (MD) simulation.

Results: Based on docking score and interaction analysis, compounds DB12380, DB13276, DB13276, DB11740, DB14675, DB11978, DB08526, DB06573, DB15796, DB08223, DB11736, DB16250, and DB16335 were predicted as the most potent. To examine the dynamic behavior and stability of the docked complexes, three compounds—DB16335, DB15796, and DB16250 —along with the Apo state were simulated for 300ns. The results revealed that compound DB16335 revealed the best docking score (-9.57 kcal/mol) against the Mpox protein thymidylate kinase.

Discussion: Additionally, during the 300 ns MD simulation period, thymidylate kinase DB16335 showed great stability. Further, in vitro and in vivo study is recommended for the final predicted compounds.

1 Introduction

The Mpox virus is widespread primarily in West and Central Africa and is a member of the Orthopoxvirus genus, which also include cowpox, vaccinia, variola, and smallpox (MacNeil et al., 2009). In Africa, various strains of Mpox have fatality rates ranging from 3.6% to 10.6% (Bunge et al., 2022), raising the possibility of increasing lethality. Recently, travelers have confirmed cases of monkeypox exported to Singapore (Yong et al., 2020), the United Kingdom (MacNeil et al., 2009), Israel, and the United States (Erez et al., 2019). Mpox is spread through direct contact with an animal’s body secretions or through animal bites. However, it can also be transmitted through breathing droplets during close and prolonged face-to-face contact, through direct contact with an infected person’s body fluids, or via objects contaminated with virus particles (Farahat et al., 2022; Lam et al., 2022).

Monkeypox is considered a sexually transmitted disease (STD) in the current outbreak. The smallpox vaccine is 85% effective at preventing monkeypox (Vallée et al., 2022). Malaise, rash, headaches, and fever between 38.5 and 40.5 degrees Celsius are some of the symptoms of monkeypox. The presence of hard, deep, and umbilicated lesions along with swollen lymph nodes are also some of the symptoms (MacNeil et al., 2009). The incubation period for Mpox lasts from one week to 17 days, with the fever going away after three days following the appearance of the rash. The lesions are described as painful, stiff, and swollen. It has been proposed that lymphadenopathy, which occurs in Mpox but not smallpox, causes a higher immune response than smallpox (Lam et al., 2022). Sepsis caused by lesions has been recorded, but it is generally believed to be uncommon (Reynolds et al., 2017). The central regions of the Mpox genome, which contain crucial enzymes and proteins, have 96.3% similarity with the smallpox genome, according to a prior study (Shchelkunov et al., 2001). The genome of Mpox is linear and has double double-stranded DNA with a genome size of 197 kb (Kugelman et al., 2014). A48R is a thymidylate kinase that has previously been found to bind with thymidine diphosphate. It is a unique target because no recognized drugs currently target it. The substantial structural difference from human thymidylate kinase at the active site makes it an attractive target (Caillat et al., 2008). According to reports, smallpox immunizations are 85% effective against monkeypox (Fine et al., 1988) (Grant et al., 2020) (Reynolds and Damon, 2012),,. Since the illness was eradicated in 1980, smallpox immunization programs have ceased (Belongia and Naleway, 2003). Although the use of smallpox medications against monkeypox is advised (Rizk et al., 2022), their safety and efficacy in human beings have not yet been determined (Sherwat et al., 2022). The present surge, therefore, indicates the urgent need for developing Mpox-specific drugs and therapies.

Finding new drugs requires several steps, including target identification, lead identification, animal studies, and clinical trials, which can typically take at least ten to twelve years to complete. Consequently, drug repurposing is a desirable substitute in the epidemic’s situation (Yong et al., 2020). The strategy of repurposing drugs has numerous advantages, including a significant decrease in testing time (Sherwat et al., 2022). In particular, drugs that have already received approval for treating different illnesses have undergone comprehensive toxicity testing and can thus be given safely to the general public (Yong et al., 2020). The drug discovery process now includes computational methods that allow the screening of small molecule libraries to identify lead candidates that can be optimized to find potential medications for clinical testing (Ashburn and Thor, 2004) (Sliwoski et al., 2014),. The crystal structure of the TMPK enzyme of Mpox was not available in the PDB database so we developed a homology model of the TMPK drug target. Furthermore, a total of 9000 FDA-approved drugs retrieved from the drug bank database were screened against the drug target. In order to computationally evaluate the stability of the ligand-protein complexes, we then carried out molecular dynamics (MD) simulations for the top three inhibitors predicted by the docking studies. In summary, our in-silico research explores whether drug that have already been approved can be good candidates against Mpox.

2 Materials and method

2.1 Structure prediction

In computational biology, predicting protein structure has been a significant scientific challenge for a long time. Regardless of the size of a protein, homology modeling has proved a relatively quick method for predicting a protein’s structure based on experimental structures available in the Protein Data Bank. AlphaFold is a trained neural network-based deep learning method for homology modeling. The most recent version of this software, AlphaFold 2, revealed extremely high accuracy in predicting protein structures (Cramer, 2021). In the present study, Google Colab was used for the prediction of the 3D structure of thymidylate kinase (TK). For the purpose of evaluating Ramachandran plots to confirm the stereochemical accuracy of the predicted protein structure, PROCHECK (Suganya et al., 2014) was used. ERRAT (Zheng et al., 2022) was also used for model validation.

2.2 Structure-based virtual screening and docking

2.2.1 Structure preparation

It is crucial to employ precise protein structures in structure-based molecular modeling. In this study, the model developed by AlphaFold 2 was used for virtual screening. The structure preparation wizard of MOE (Molecular Operating Environment 2016) software was used to prepare the structure (Ahmed et al., 2022). The protein structure was 3D protonated to add hydrogen atoms, and then the MOE software default parameters were utilized to minimize the energy of the structure.

2.2.2 Ligands preparation

Over 9000 compounds documented in the drug-bank database were retrieved from the drug-bank database. All of the compounds in the drug-bank database were subjected to three-dimensional protonation using the MMFF94x force field and energy minimization using MOE software with an RMS gradient of 0.05 (Rauf et al., 2022).

2.2.3 Molecular docking

Structure-based virtual screening (SBVS) techniques require the 3D structures of the drug target (receptor) and ligands in the database. We used molecular docking to evaluate the binding patterns of drug-target proteins and ligands in order to discover novel potential inhibitors. The interactions between ligands and receptors can be predicted by molecular docking approaches. The molecular docking investigations were carried out using the MOE docking program (Taj et al., 2022). The retrieved compounds were docked with the active site of thymidylate kinase to predict the active compounds against thymidylate kinase (A48R). For each compound, a total of five conformers were generated, and the top-ranked conformation for each compound was analyzed for interaction analysis. Using the rigid receptor docking protocol and GBVI/WSA scoring function, molecular docking (Taha et al., 2023) of all FDA-approved drugs was performed. AutoDock Vina software was used to carry out the docking operation (Ali et al., 2023) of the lowest docking score hits obtained from the MOE software to validate the docking protocol of the MOE software. Adding polar hydrogen atoms and assigning partial charges to each atom were all done using the AutoDock software, and the structure of ligands and receptors were saved in PDBQT format. The grid spacing was set to 1 Å and the box size was 16×20 ×24. Twenty conformations were set in the docking output. The binding poses of the molecule with the best docking energies were visualized using the PyMol software (Du et al., 2023).

2.3 MD simulation

Molecular dynamics simulations were performed using the AMBER 20 software suite (Ajmal et al., 2022). The best docking score complexes were used as the starting structures for the MD simulations. The AMBER force field FF14SB was used for the protein, and the generic AMBER force field (GAFF) was used for the ligands (Wadood et al., 2022). The topology files and atomic charges of ligands were generated using the Antechamber suite in the AMBER 20 package (Mahmood et al., 2023). The topology and coordinate data for the entire system were generated using the tleap module of the AMBER 20 software. The entire system was soaked in a TIP3P water box at a margin distance of 8 Å. An appropriate quantity of chloride ions was introduced in order to neutralize the system’s charge. During the simulation, the particle mesh Ewald (PME) was used to handle the long-range electrostatic interactions, and the cut-off distance for nonbonded interactions was adjusted to 10 Å. The SHAKE algorithm (Poli et al., 2022) was employed to constrain the hydrogen-containing bonds. Each system underwent two stages of energy minimization: the first stage involved the use of the algorithms (10,000 steps of the steepest descent and 10,000 steps of the conjugate gradient) with restraint, and the second stage involved the use of the same algorithms without restraint. Then, each system was gradually heated from 0 to 300 K. The system was then brought up to equilibrium at 300 K and constant pressure. Finally, a production process lasting 200 ns was carried out under conditions of constant temperature and pressure (NTP). Finally, the CPPTRAJ was used for trajectory analysis (Junaid et al., 2018).

2.3.1 Principal component analysis

Proteins’ high-amplitude movements can be captured by using principal component analysis (PCA). The PCA analysis in this study was performed using the cpptraj package (Ajmal et al., 2022). The covariance matrix was generated based on the Cartesian coordinates of the carbon alpha of each system to evaluate the dynamic behavior of each system. The covariance matrix’s eigenvalues and eigenvectors were obtained by diagonalization. The eigenvectors and eigenvalues, respectively, showed the direction of high-amplitude motion and their mean square fluctuation. The PC1 and PC2, for each system, were calculated and plotted to track their movements.

2.4 Binding energy calculation

The free energy of binding between ligands and protein complexes was computed by the MMGBSA and MMPBSA techniques. For the calculation of binding energy, the last 100 frames were used. Two effective methods for analyzing the free energy of binding are MM/GBSA and MM/PBSA. The values of MM/PBSA strongly correlate with experimental methods (33). Here, we calculated the binding free energy using both the MMPBSA and MMGBSA method. The following equation was employed to calculate the free energy:

Gbind=ΔGcomplex[ΔGreceptor+ΔGligand]

In the above equation, Gbind stands for the total binding free energies, ΔGcomplex for complex free energies, and the remaining terms stand for the corresponding free energies of the receptor protein and ligand.

3 Results

3.1 Protein sequence retrieval

The protein sequence of the Mpox virus thymidylate kinase was retrieved from the NCBI database.

3.2 Homology modeling

In addition, the 3D structure was predicated by the AlphaFold2 server as presented in Figure 1. A per-residue confidence rating between 0 and 100 was produced by AlphaFold2. There may be some disorder regions in the predicted structure which can be represented by a low pLDDT value. In the generated model, most of the residues have very high confidence scores that represent the model accuracy (pLDDT > 90, Figure 1). The dark green region in the predicted aligned error plot is ideal, indicating the accuracy of the model, but light green is bad (indicating a high error).

FIGURE 1
www.frontiersin.org

Figure 1 3D structure of thymidylate kinase, a drug target in Mpox, generated by AlphaFold 2 server.

3.3 Structure validation

The structure was validated by PROCHECK and ERRAT servers. PROCHECK is used for three-dimensional model validation. The Ramachandran plot is shown in Figure S1A. In the plot, residues were determined in three various regions such as favored region (94%), allowed region (4.3%), and outlier region (0.5%). Furthermore, the ERRAT server was used to check the quality of the model. ERRAT validated the model by using statistical relation of non-bonded interactions between various atom types on basis of atomic interaction. The values around 95% and higher than 95% indicate standard high-resolution structures; however, values around 91% indicate low-resolution structures. Figure S2B illustrates a 97% ERRAT score, indicating the best model according to the ERRAT plot.

3.4 Molecular docking analysis

In the present study, molecular docking was performed to evaluate the interaction of FDA-approved drugs against the drug target thymidylate kinase. By employing the site finder option of MOE software, active site residues were identified. The docking scores of all the compounds were arranged in ascending order. Twenty percent of the data was further used for interaction analysis. Among these compounds, five compounds exhibited the most active compounds against thymidylate kinase. As noted in Figure 2A, compound DB15796 mediates a docking score of -7.49 kcal/mol. Compound DB15796 made six hydrogen bonds with the active site residues of thymidylate kinase (Asp13, Asp92, Glu145, Lys17, Arg41, Arg93), and one π-H (Thr18) interaction (Table 1). Moreover, in Figure 2B, the compound DB08223 exhibited six hydrogen bonds with Lys17, Gly16, Thr18, and Arg93 with a docking score of -7.12 kcal/mol. The docking complex of DB11736 and A48R exhibits hydrogen bonds with Thr18, Lys14, and Gly16, one ionic (Glu142), and one Pi-cation (Arg93) interaction with a docking score of -7.29 kcal/mol as depicted in Figure 2C. In addition, the complex of DB16250 and A48R protein, as shown in Figure 2D, mediates one hydrogen bond donor (Asp13), three hydrogen bonds acceptor (Gly16, Lys17, Arg93), one ionic bond (Glu142), and one π-cation (Thr18) interactions. The docking score of the compound DB16250 was predicted as -8.65 kcal/mol. Finally, the compound DB16335 exhibits the best docking score of -9.57 kcal/mol and made seven hydrogen bonds with Lys 14, lys17, Thr 18, Asn37, Arg93, Phe38, and Glu142 residues of thymidylate kinase (Figure 2E) (Table 1). A number of studies reported that the success rate of virtual screening can be increased by combining the results of different docking software (36). To increase the success rate of virtual screening, the top best docking score hits obtained from the MOE software were further docked by the AutoDock software. The binding energy values obtained from AutoDock Vina are presented in Table 2. The binding energy obtained from AutoDock and the S score obtained from the MOE software is almost similar and ranges from -7.5 to -9.4 kcal/mol. The docking results also revealed that by performing the docking with MOE and AutoDock software, the residues including Asp13, Lys17, Thr18, Arg93, Arg41, and Gly16 were involved in the hydrogen bond formations with the ligands atoms.

FIGURE 2
www.frontiersin.org

Figure 2 3D interactions of top five final compounds: (A) DB15796; (B) DB08223; (C) DB11736; (D) DB16250; and (E) DB16335 in complex with thymidylate kinase. The green dotted line represents the bond.

TABLE 1
www.frontiersin.org

Table 1 Docking score and protein-ligands interaction of best-scored compounds.

TABLE 2
www.frontiersin.org

Table 2 The binding energy and interacted residues of the docked compounds obtained from the AutoDock vina software.

3.5 Molecular dynamic simulation analysis

Molecular dynamic simulation helps to study conformational changes of compounds (ligands) into the binding pocket of a protein. On the basis of docking score and interactions, the top three compounds—DB16335, DB15796, DB16250—and Apo state of thymidylate kinase were subjected to MD simulation in which dynamic behavior and stability of complexes were examined by utilizing Amber 20 software. Further, RMSD values were analyzed to evaluate the stability of complexes.

3.5.1 Stability analysis

In the Apo state, initially, RMSD was drastically high at 2.0 Å up to 70 ns. After that, the RMSD value decreases to 1.4 Å from 80 to 230 ns. The RMSD value then gradually increases to 2.2 Å up to 300 ns. The average RMSD of the Apo-state was found to be 2.2 Å. The thymidine kinase DB16335 compound has a minimum RMSD value of 1.2-1.5 Å up to 160 ns, representing maximum stability in trajectories, however during 160 to 170 ns, small fluctuations were observed, and then after 180 to 300 ns, the complex revealed great stability and remained stable throughout the 300 ns (Figure 3A). The thymidine kinase DB15796 initially showed stable behavior with an RMSD value of 1.2 Å up to 70 ns, then the RMSD value gradually increases to 2.2 Å up to 80-250ns, and finally, becomes stable at 300 ns (Figure 3B). The thymidylate kinase DB16250 shows a steady state behavior with an RMSD value of ∼1.5 Å at 25 ns up to 230 ns. After that, the RMSD increases to 2.4 Å up to 270 ns and then reaches a stable condition at 300 ns (Figure 3C). Among all the complexes, the RMSD of the thymidine kinase DB16335 complex was more stable.

FIGURE 3
www.frontiersin.org

Figure 3 Black-colored lines in all graphs represent the Apo state. (A) Root-mean-square deviation (RMSD) graph of thymidylate kinase DB16335, (B) thymidylate kinaseDB15796, (C) graph of thymidylate kinase DB16250.

3.5.2 Residues flexibility index analysis

RMSF analysis was performed in order to evaluate the contributions of every amino acid to the stability of complexes. The residues from 50-53 and -150-153 show major fluctuations in all complexes during simulation (Figure 4). Smaller fluctuations represent a stable complex. The overall result of RMSF indicates that thymidine kinase DB16335 and thymidine kinase DB15796 show low fluctuations compared with Apo; however, thymidylate kinase DB16250 shows more similar fluctuations to that of the Apo state. Residues Asp 50, Asp51, Tyr 52, Leu 53, Gln 151, Lys 152, and Val 153 revealed high fluctuations during the 300 ns MD simulation. The RMSF plot is shown in Figure 4.

FIGURE 4
www.frontiersin.org

Figure 4 Black-colored lines in all graphs represent the Apo state. (A) Root-mean-square deviation (RMSD) graph of thymidylate kinase DB16335; (B) thymidylate kinase DB15796; (C) thymidylate kinase DB16250.

3.5.3 Compactness analysis

Folded and unfolded states of protein complexes were predicated by RoG. Our result revealed that the final three complexes along with the Apo state represent change behaviors throughout molecular dynamic simulations. In addition, we also analyze the radius of gyration (RoG) that helps to measure the compactness of the complexes. The average ROG values of thymidine kinase DB16335, thymidine kinase DB15796; thymidylate kinase DB16250; and Apo were 16.5 ± 16.7; 16.5 ± 16.9; 16.6 ± 17.0; and 16.7± 17.1 A, respectively (Figure 5). Among all the complexes as well as the Apo state, the complex thymidine kinase DB16335 revealed a highly compact complex.

FIGURE 5
www.frontiersin.org

Figure 5 Radius of gyration (RoG) plot of Apo (black) in all graphs. (A) ROG graph of thymidylate kinase DB16335; (B) thymidylate kinase DB15796; (C) thymidylate kinase DB16250.

3.5.4 Principal component analysis

A dimensionality reduction technique such as principal component analysis was performed to examine the mobility of the proteins and the clusters of the relevant structural frames.

A statistical technique known as principal component analysis (PCA) combines several correlated variables with a smaller set of uncorrelated variables called principal components. The PCA was calculated and is displayed in Figure 6 in order to properly evaluate the effect of the drug binding on protein mobility. The first two principal components, known as PC1 and PC2, were plotted against one another. Colors ranging from red to blue reflect how one conformation changes into another. Each dot in Figure 6 corresponds to a single frame of the trajectory. The coordinate covariance matrix, which was derived from the time series of 3D positional coordinates of the complexes over the course of the 200 ns MD simulation, served as the input for principal component analysis. The Apo state displayed a slightly dispersed form of motion, while the drug complexes displayed a cluster type of motion. Among all the complexes TK DB16335 revealed more cluster types of motion.

FIGURE 6
www.frontiersin.org

Figure 6 PCA plot for (A) thymidylate kinase-Apo state; (B) thymidylate kinase DB16335; (C) thymidylate kinase DB15796; (D) thymidylate kinase DB16250.

3.6 Binding energy calculations

Finally, we calculated ΔG binding in order to refine ligands predicted through virtual screening by using MM/PBSA and MM/GBSA techniques. The MMGBSA analysis composed of electrostatic energy, van der Waals energy, surface area energy, and electrostatic contribution to the solvation-free energy was calculated as presented in Table 3, while Table 4 represents the MMPBSA analysis. Among all the complexes, DB16335 in complex with thymidylate kinase revealed lower binding energy and high binding affinity for thymidylate kinase, a drug target in the Mpox virus predicted by both MMPBSA and MMGBSA methods.

TABLE 3
www.frontiersin.org

Table 3 MMGBSA analysis reported in kcal/mol.

TABLE 4
www.frontiersin.org

Table 4 MMPBSA analysis reported in kcal/mol.

4 Discussion

Mpox virus, a member of Orthopoxvirus genus that includes smallpox, has become endemic to Africa and throughout the world. No vaccines against Orthopoxviruses have been developed in the last four decades, after the eradication of smallpox. The reemergence of monkeypox in unprepared and unvaccinated populations is an emergency demanding attention from the scientific community. Instead of infeasible mass vaccination campaigns in a short time period, specific therapeutics against Mpox can provide sustainable solutions for affected patients (Potter et al., 2007). Approximately 49 genes are common among all the members of the provirus family from the total 150 genes encoded by poxviruses (Llanos et al., 2021). Currently, no approved drug is available for the treatment of monkeypox virus (Anwar et al., 2023).

Drug repurposing is best way to determine effective therapeutics against Mpox. Therefore, the present study aims to identify the best drugs against Mpox by using different computational approaches such as homology modeling, virtual screening, molecular docking, and MD simulation. For this purpose, we collected FDA-approved drugs, and 3D structure of protein was modeled. The model structure was subjected to molecular docking against a curated library of FDA-approved drugs. Afterwards, by identifying binding energy and interaction, various compounds were selected in which drugs DB12380, DB13276, DB13276, DB11740, DB14675, DB11978, DB08526, DB06573, DB15796, DB08223, DB11736, DB16250, and DB16335 mediate significant binding energy and binding interaction with thymidylate kinase (A48R). Among these compounds, the compounds having the lowest binding energy—DB15796, DB08223, DB11736, DB16250, and DB16335, with -7.49, -7.12, -7.29, 8.65, and -9.57 kcal/mol, respectively—were recognized as the best compounds (Table 1).

Finally, compounds DB16335, DB15796, and DB16250 along with Apo were subjected to 300 ns MD simulation to evaluate the stability of the compounds. Among the Apo and selected compounds, DB16335 was found to be active and showed strong binding against Mpox. DB16335 is an antibiotic drug that is under a clinical trial NCT03354598) and is used to treat complicated and uncomplicated urinary tract infections, intra-abdominal infections, pneumonia, and acute pyelonephritis. In addition, it is also active against Gram-positive and Gram-negative bacteria (Dunne et al., 2023). The selected compounds have low RMSD and low amino acid residue fluctuations compared with Apo, which indicates that the binding of selected compounds to the active site of Mpox protein makes the protein more stabilized. In addition, the ROG also indicates that thymidylate kinase DB16335 has a low value of RoG (Figure 5). The binding of these compounds to Mpox can destroy its activity and provide an unfavorable situation for Mpox.

5 Concluding remarks

Monkeypox has rapidly spread through the world, but still, no suitable vaccine or drug is available. The present study is designed to determine a potent drug against Mpox by using computational methods. Our result shows that all three compounds were notable because of good interaction and binding energy. Among all the drugs, DB16335 was identified as the most potent compound in terms of docking score, interactions, and binding energy calculation. Overall, compound DB16335 could bind and inhibit the drug target thymidylate kinase, which can be helpful to reduce infections associated with the Mpox virus. However, further in vitro and in vivo experiments are required to validate the findings of our study.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

AA, AM and CH performed experiments, drafted the manuscript. MAH performed some of the experiments and revised the manuscript and discussion. BSA, PL, MU, ANA and PH drafted the manuscript and critically revised the manuscript. AW and JH conceptualized, designed, and supervised the study and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: 23UQU4331128DSR001. This study is also funded by the Incubation Fund of SSL Central Hospital of Dongguan (YJ002).

Acknowledgments

The authors thank Dr. Ajmal, University of Nizwa, Sultanate Oman for providing technical support for the computational study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1159389/full#supplementary-material

References

Ahmed, A., Saeed, A., Ejaz, SA., Aziz, M., Hashmi, M. Z., Channar, P. A., et al. (2022). Novel adamantyl clubbed iminothiazolidinones as promising elastase inhibitors: design, synthesis, molecular docking, ADMET and DFT studies. RSC Adv. 12 (19), 11974–11991. doi: 10.1039/D1RA09318E

PubMed Abstract | CrossRef Full Text | Google Scholar

Ajmal, A., Ali, Y., Khan, A., Wadood, A., Rehman, A. U. (2022). Identification of novel peptide inhibitors for the KRas-G12C variant to prevent oncogenic signaling. J. Biomolecular Structure Dynamics p, 1–10. doi: 10.1080/07391102.2022.2138550

CrossRef Full Text | Google Scholar

Ali, Y., Imtiaz, H., Tahir, M. M., Gul, F., Saddozai, U. A., ur Rehman, A., et al. (2023). Fragment-based approaches identified tecovirimat-competitive novel drug candidate for targeting the F13 protein of the monkeypox virus. Viruses 15 (2), 570. doi: 10.3390/v15020570

PubMed Abstract | CrossRef Full Text | Google Scholar

Anwar, F., Haider, F., Khan, S., Ahmad, I., Ahmed, N., Imran, M., et al. (2023). Clinical manifestation, transmission, pathogenesis, and diagnosis of monkeypox virus: a comprehensive review. Life 13 (2), 522. doi: 10.3390/life13020522

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburn, T. T., Thor, K. B. (2004). Drug repositioning: identifying and developing new uses for existing drugs. Nat. Rev. Drug Discov. 3 (8), 673–683. doi: 10.1038/nrd1468

PubMed Abstract | CrossRef Full Text | Google Scholar

Belongia, E. A., Naleway, A. L. (2003). Smallpox vaccine: the good, the bad, and the ugly. Clin. Med. Res. 1 (2), 87–92. doi: 10.3121/cmr.1.2.87

PubMed Abstract | CrossRef Full Text | Google Scholar

Bunge, E. M., Hoet, B., Chen, L., Lienert, F., Weidenthaler, H., Baer, L. R., et al. (2022). The changing epidemiology of human monkeypox–a potential threat? a systematic review. PloS Negl. Trop. Dis. 16 (2), e0010141. doi: 10.1371/journal.pntd.0010141

PubMed Abstract | CrossRef Full Text | Google Scholar

Caillat, C., Topalis, D., Agrofoglio, L. A., Pochet, S., Balzarini, J., Deville-Bonne, D., et al. (2008). Crystal structure of poxvirus thymidylate kinase: an unexpected dimerization has implications for antiviral therapy. Proc. Natl. Acad. Sci. 105 (44), 16900–16905. doi: 10.1073/pnas.0804525105

CrossRef Full Text | Google Scholar

Cramer, P. (2021). AlphaFold2 and the future of structural biology. Nat. Struct. Mol. Biol. 28 (9), 704–705. doi: 10.1038/s41594-021-00650-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, J., Wang, S., Zhang, X., Liu, C., Zhang, Y., Zhang, H., et al. (2023). The inhibitory mechanism of 7H-pyrrolo [2, 3-d] pyrimidine derivatives as inhibitors of P21-activated kinase 4 through molecular dynamics simulation. Molecules 28 (1), 413. doi: 10.3390/molecules28010413

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunne, M. W., Aronin, S. I., Das, A. F., Akinapelli, K., Breen, J., Zelasky, M. T., et al. (2023). Sulopenem for the treatment of complicated urinary tract infections including pyelonephritis: a phase 3, randomized trial. Clin. Infect. Dis. 76 (1), 78–88. doi: 10.1093/cid/ciac704

PubMed Abstract | CrossRef Full Text | Google Scholar

Erez, N., Achdout, H., Milrot, E., Schwartz, Y., Wiener-Well, Y., Paran, N., et al. (2019). Diagnosis of imported monkeypox, Israel, 2018. Emerging Infect. Dis. 25 (5), 980. doi: 10.3201/eid2505.190076

CrossRef Full Text | Google Scholar

Farahat, R. A., Sah, R., El-Sakka, A. A., Benmelouka, A. Y., Kundu, M., Labieb, F., et al. (2022). Human monkeypox disease (MPX). Le infezioni in Medicina 30 (3), 372–391. doi: 10.53854/liim-3003-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fine, P., Jezek, Z., Grab, B., Dixon, H. (1988). The transmission potential of monkeypox virus in human populations. Int. J. Epidemiol. 17 (3), 643–650. doi: 10.1093/ije/17.3.643

PubMed Abstract | CrossRef Full Text | Google Scholar

Grant, R., Nguyen, L.-B. L., Breban, R. (2020). Modelling human-to-human transmission of monkeypox. Bull. World Health Organ. 98 (9), 638. doi: 10.2471/BLT.19.242347

PubMed Abstract | CrossRef Full Text | Google Scholar

Junaid, M., Shah, M., Khan, A., Li, C. D., Khan, M. T., Kaushik, A. C., et al. (2018). Structural-dynamic insights into the h. pylori cytotoxin-associated gene a (CagA) and its abrogation to interact with the tumor suppressor protein ASPP2 using decoy peptides. J. Biomolecular Structure Dynamics.

Google Scholar

Kugelman, J. R., Johnston, S. C., Mulembakani, P. M., Kisalu, N., Lee, M. S., Koroleva, G., et al. (2014). Genomic variability of monkeypox virus among humans, democratic republic of the Congo. Emerging Infect. Dis. 20 (2), 232. doi: 10.3201/eid2002.130118

CrossRef Full Text | Google Scholar

Lam, H. Y. I., Guan, J. S., Mu, Y. (2022). In silico repurposed drugs against monkeypox virus. Molecules 27 (16), 5277. doi: 10.3390/molecules27165277

PubMed Abstract | CrossRef Full Text | Google Scholar

Llanos, M. A., Gantner, M. E., Rodriguez, S., Alberca, L. N., Bellera, C. L., Talevi, A., et al. (2021). Strengths and weaknesses of docking simulations in the SARS-CoV-2 era: the main protease (Mpro) case study. J. Chem. Inf. Modeling 61 (8), 3758–3770. doi: 10.1021/acs.jcim.1c00404

CrossRef Full Text | Google Scholar

MacNeil, A., Reynolds, M. G., Braden, Z., Carroll, D. S., Bostik, V., Karem, K., et al. (2009). Transmission of atypical varicella-zoster virus infections involving palm and sole manifestations in an area with monkeypox endemicity. Clin. Infect. Dis. 48 (1), e6–e8. doi: 10.1086/595552

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahmood, A., Samad, A., Shah, A., Wadood, A., Alkathiri, A., Alshehri, M. A., et al. (2023). A novel biallelic variant in the popeye domain-containing protein 1 (POPDC1) underlies limb girdle muscle dystrophy type 25. Clin. Genet. 103 (2), 219–225. doi: 10.1111/cge.14238

PubMed Abstract | CrossRef Full Text | Google Scholar

Poli, G., Barravecchia, I., Demontis, G. C., Sodi, A., Saba, A., Rizzo, S., et al. (2022). Predicting potentially pathogenic effects of h RPE65 missense mutations: a computational strategy based on molecular dynamics simulations. J. Enzyme Inhibition Medicinal Chem. 37 (1), 1765–1772. doi: 10.1080/14756366.2022.2090547

CrossRef Full Text | Google Scholar

Potter, M. A., Sweeney, P., Iuliano, A.D., Allswede, M. P. (2007). Performance indicators for response to selected infectious disease outbreaks: a review of the published record. J. Public Health Manage. Pract. 13 (5), 510–518. doi: 10.1097/01.PHH.0000285205.40964.28

CrossRef Full Text | Google Scholar

Rauf, A., Ibrahim, M., Muhammad, N., Naz, S., Wadood, A., Khan, B., et al. (2022). Enzyme inhibitory activities of extracts and carpachromene from the stem of ficus benghalensis. BioMed. Res. Int. 2022.

Google Scholar

Reynolds, M. G., Damon, I. K. (2012). Outbreaks of human monkeypox after cessation of smallpox vaccination. Trends Microbiol. 20 (2), 80–87. doi: 10.1016/j.tim.2011.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Reynolds, M. G., McCollum, A. M., Nguete, B., Shongo Lushima, R., Petersen, B. W. (2017). Improving the care and treatment of monkeypox patients in low-resource settings: applying evidence from contemporary biomedical and smallpox biodefense research. Viruses 9 (12), 380. doi: 10.3390/v9120380

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizk, J. G., Lippi, G., Henry, B. M., Forthal, D. N., Rizk, Y. (2022). Prevention and treatment of monkeypox. Drugs 82 (9), 957–963. doi: 10.1007/s40265-022-01742-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Shchelkunov, S. N., Totmenin, A. V., Babkin, I. V., Safronov, P. F., Ryazankina, O. I., Petrov, N. A., et al. (2001). Human monkeypox and smallpox viruses: genomic comparison. FEBS Lett. 509 (1), 66–70. doi: 10.1016/S0014-5793(01)03144-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Sherwat, A., et al. (2022). Tecovirimat and the treatment of monkeypox–past, present, and future considerations. New Engl. J. Med. 387 (7), 579–581. doi: 10.1056/NEJMp2210125

CrossRef Full Text | Google Scholar

Sliwoski, G., Kothiwale, S., Meiler, J., Lowe, E. W. (2014). Computational methods in drug discovery. Pharmacol. Rev. 66 (1), 334–395. doi: 10.1124/pr.112.007336

PubMed Abstract | CrossRef Full Text | Google Scholar

Suganya, P. R., Sudevan, K., Kalva, S., Saleena, L. M. (2014). Homology modeling for human adam12 using prime, i-tasser and easymodeller. Int. J. Pharm. Pharm. Sci. 6 (2), 782–786.

Google Scholar

Taha, M., Rahim, F., Khan, I. U., Uddin, N., Farooq, R. K., Wadood, A., et al. (2023). Synthesis of thiazole-based-thiourea analogs: as anticancer, antiglycation andantioxidant agents, structure activity relationship analysis and docking study. J. Biomolecular Structure Dynamics p, 1–16. doi: 10.1080/07391102.2023.2171134

CrossRef Full Text | Google Scholar

Taj, S., Ahmad, M., Alshammari, A., Alghamdi, A., Ashfaq, U. A. (2022). Exploring the therapeutic potential of benzothiazine-pyrazole hybrid molecules against alpha-glucosidase: pharmacological and molecular modelling based approach. Saudi J. Biol. Sci. 29 (3), 1416–1421. doi: 10.1016/j.sjbs.2021.11.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Vallée, A., Farfour, E., Zucman, D. (2022). Monkeypox virus: a novel sexually transmitted disease? a case report from France. Travel Med. Infect. Dis. 49, 102394. doi: 10.1016/j.tmaid.2022.102394

PubMed Abstract | CrossRef Full Text | Google Scholar

Wadood, A., Ajmal, A., Junaid, M., Rehman, A. U., Uddin, R., Azam, S. S., et al. (2022). Machine learning-based virtual screening for STAT3 anticancer drug target. Curr. Pharm. Design 28 (36), 3023–3032. doi: 10.2174/1381612828666220728120523

CrossRef Full Text | Google Scholar

Yong, S. E. F., Ng, OT., Ho, Z. J., Mak, T. M., Marimuthu, K., Vasoo, S., et al. (2020). Imported monkeypox, Singapore. Emerging Infect. Dis. 26 (8), 1826. doi: 10.3201/eid2608.191387

CrossRef Full Text | Google Scholar

Zheng, W., Su, R., Lin, X., Liu, J. (2022). Nanochannel array modified three-dimensional graphene electrode for sensitive electrochemical detection of 2, 4, 6-trichlorophenol and prochloraz. Front. Chem. 10. doi: 10.3389/fchem.2022.954802

CrossRef Full Text | Google Scholar

Keywords: monkeypox, homology modeling, molecular docking, MD simulation, drugs development

Citation: Ajmal A, Mahmood A, Hayat C, Hakami MA, Alotaibi BS, Umair M, Abdalla AN, Li P, He P, Wadood A and Hu J (2023) Computer-assisted drug repurposing for thymidylate kinase drug target in monkeypox virus. Front. Cell. Infect. Microbiol. 13:1159389. doi: 10.3389/fcimb.2023.1159389

Received: 05 February 2023; Accepted: 10 May 2023;
Published: 29 May 2023.

Edited by:

Curtis Brandt, University of Wisconsin-Madison, United States

Reviewed by:

Yasir Ali, Quaid-i-Azam University, Pakistan
Sobia Ahsan Halim, University of Nizwa, Oman

Copyright © 2023 Ajmal, Mahmood, Hayat, Hakami, Alotaibi, Umair, Abdalla, Li, He, Wadood and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Abdul Wadood, awadood@awkum.edu.pk; Junjian Hu, hujunjian79@163.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.