1. Introduction
The Lumpy skin disease (LSD) is a contagious viral disease that affects ruminant animals such as cattle, water buffaloes, and giraffes, posing a significant threat across borders [
1,
2,
3,
4]. This disease is attributed to the Lumpy skin disease virus (LSDV), a member of the Capripoxvirus genus within the Poxviridae family that exhibits close kinship with the sheep pox virus (SPPV) and goat pox viruses (GTPV) [
5]. It has a genome spanning 151 kilobase pairs (kb) and 156 putative genes [
6]. Similar to other poxviruses, LSDV possesses several conserved genes involved in its fundamental replicative mechanisms. Seven homologues of genes found in chordopoxviruses are associated with DNA replication. These genes, namely, LSDV039, LSDV077, LSDV082, LSDV083, LSDV112, LSDV133, and LSDV139, are either known or potentially involved in DNA replication processes [
6]. However, the virus is most likely spread mechanically by blood-sucking arthropods such as flies, mosquitoes, and ticks. Direct contact between animals can also transmit the virus to a lesser extent [
7,
8]. Additionally, infected animals can spread the virus through milk, blood, nasal secretions, and saliva, which serve as alternative transmission routes through feeding or drinking [
9]. Affected animals primarily exhibit symptoms such as fever, nodular skin lesions, dramatic decrease in milk production, and weight loss [
10]. The introduction of LSDV into a herd can result in a high incidence rate, ranging from 5 to 45%, and a case death rate from 0.5% to 7.0% [
11,
12]. Consequently, LSD poses a significant economic threat to the livestock industry worldwide because it causes vast economic losses, including abortions in females and sterility in males [
13].
It was first reported in Africa (Zambia) in 1929 and was identified as a communicable disease in the 1940s [
14]. It has also spread to countries in Central and Eastern Africa, the Middle East, Asia, and Eastern Europe [
15]. LSD outbreak data were accessed from the World Organization for Animal Health (WOAH) World Animal Health Information System (WAHIS) database in Southeast Asian countries during the study period between October 2020 and October 2021. During the epidemic period, 866 LSD outbreaks were reported in six Southeast Asian countries, including 1,758,923 susceptible cattle, 93,465 cases, 5,936 deaths, and 1,117 culled cattle [
16]. In Bangladesh, the outbreak first emerged in the Chattogram region in July 2019 and then quickly spread throughout the entire country [
17].
To date, the scientific community has preferred some homologous and heterologous vaccines to develop immunity owing to their cross-protection against lumpy skin diseases, such as Kenyan sheep and goat pox (KSGP) O-180 strain vaccines and Gorgan goat pox (GTP) vaccines. While the KSGPO-180 vaccine failed to protect cattle against LSDV, the Gorgan GTP vaccine successfully prevented the clinical symptoms of LSD in all vaccinated calves [
18]. Certain SPPV vaccines provide only partial or incomplete protection against LSDV. Despite the use of goat pox vaccine virus KSGPO-240 and Romanian SPPV strains, their attenuation levels were insufficient to ensure the safety of cattle. Consequently, vaccinated animals remain ill [
19]. The Neethling vaccine was developed in 2016 and tested in six Balkan countries. Its effectiveness was found to be average at 79.8% (95% CI: 73.2–84.7) in the Balkan region, with varying rates in different countries. Albania had a 62.5% effectiveness rate, whereas Bulgaria and Serbia had effectiveness rates of 97% [
20].
Over the past three decades, many diseases have been treated using repurposed drugs. Zidovudine (AZT), the first successful repurposed drug used to treat Human Immunodeficiency Virus (HIV) in 1987, was originally developed to treat cancer [
21,
22]. Recently, a drug named Remdezivir (RDV), a repurposed drug approved by Food and Drug Administration (FDA) in 2020, was used to treat Covid-19 patients [
23]. Additionally, a recent in vitro study showed that ivermectin, an antiparasitic drug, is effective against capripoxviruses [
24]. This study found that ivermectin demonstrated significant inhibitory effects on viral replication and the attachment and penetration stages of the LSDV virus. Specifically, the results showed that ivermectin reduced viral replication by 99.82 and 99.87% at the replication stage, respectively. It also exhibited inhibitory effects of 68.38 and 25.01% at the attachment stage and 57.83 and 0.0% at the penetration stage [
24]. Nevertheless, plant-based phytocompounds with antiviral activity are suitable natural solutions and show high efficacy in inhibiting several viral diseases, such as HIV and Covid-19 [
25]. For the alternative treatment approach, propolis-alginate nanoparticles (Propolis-ALg NPs) have a potential therapeutic approach through different routes, including eye drops, oral routes, and topical spray. Transmission electron microscopy was used to characterize the propolis-ALg NPs. Propolis-ALg NPs effectively treat infected animals by reducing fever and boosting overall health [
13]. Furthermore, colchicine, a potent natural alkaloid, has the potential to effectively treat a range of skin diseases, either as a standalone therapy or in combination with other drugs [
26]. In this study, ivermectin was used as a control drug to serve as a benchmark for comparison with repurposed drugs and plant-based compounds investigated for their antiviral properties against LSDV. Although these have demonstrated some level of effectiveness, the efficacy and safety of these treatments for LSDV remain unclear. Considering these constraints, it is crucial to identify a specific drug that can effectively inhibit or reduce LSDV infection [
27,
28].
Currently, there is no cure for LSD, and available treatments mainly focus on relieving the symptoms of the disease [
13]. Researchers are exploring new drug targets for LSD to develop more effective treatments. Computer-aided drug design is a valuable tool in this process because it allows for the design of new drugs based on the structure of protein targets within the virus [
29]. Computational methods have been employed to analyze the protein targets of the LSD virus and design novel lead molecules or repurposed drugs for drug development. Target selection is a critical step in drug discovery, particularly for viruses with evasion mechanisms [
29]. The identification of highly conserved and essential targets for viral survival and replication is crucial for successful drug development. Therefore, we targeted the DNA replication associated with DNA replication, including the DNA polymerase LSDV039, which is potentially involved in DNA replication.
Considering the morbidity and mortality, a truly effective antiviral drug against LSDV has emerged. Repurposed drugs [
27,
28,
30], which can be rapidly used, are readily available, and show significance as antivirals, are nothing but a great choice to eradicate LSDV. Hence, this study aimed to identify LSDV inhibitors that could effectively inhibit viral replication. Virtual screening and molecular modeling techniques enable the identification of compounds (plant-based compounds and repurposed drugs) with high binding affinities for viral proteins. Further experimental validation will provide valuable insights into the effectiveness and safety of potential inhibitors and pave the way for the development of new therapeutic strategies against LSDV.
4. Discussion
LSD is a viral disease that affects livestock and is considered a transboundary animal disease because of its ability to cross national borders, affecting multiple countries and regions. LSDV has significantly expansion of LSDV beyond its original endemic regions in Africa. LSDV have been reported in various geographical locations worldwide, including Asian, European, and Middle Eastern countries [
70]. This disease was first reported in Bangladesh in July 2019 [
71], China in August 2019 [
72], India in November 2019 [
73], Nepal in June 2020 [
74], and other Asian countries [
75]. In recent years, multiple cluster outbreaks of LSD have occurred on the Asian subcontinent, leading to substantial cattle mortality and posing a serious concern for the agricultural and livestock sectors [
76]. For instance, in India alone, over 155,000 cattle deaths will be reported by 2022, underscoring the urgent need for disease control and management [
77].
According to a report by the Department of Livestock Services (DLS) in 2019, 553,528 cattle were affected by LSD in all eight divisions of Bangladesh. The highest incidences were observed in Chattogram (8.26%) and Khulna (6.52%), whereas the lowest was observed in Sylhet (0.01%) (DLS, 2019) [
78]. In one study, the LSDV collected from two Bangladeshi strains, BD-V392.1 and BD-V395.1, were distinct from contemporary field strains found in other Asian countries such as Hong Kong, China, Taiwan, and Vietnam. This suggests a unique genetic lineage or origin for Bangladeshi LSDV strains, which differs from strains circulating in neighboring Asian countries [
79]. The Livestock Research Institute (LRI) in Bangladesh provides vaccines that are used as facultative measures, meaning that it is not mandatory but optional for farmers to vaccinate their livestock against LSDV. However, vaccination coverage is limited, resulting in many small holdings remaining unvaccinated [
79]. Furthermore, the genetic variation of LSDV is one of the reasons why vaccines are not effective in a specific geographical region. Currently, no subunit vaccinations or chemotherapeutic medications are available for the treatment of LSD. It is currently feasible to prevent or treat this disease using other therapies such as phytochemical-based medicines or repurposed medications.
In this study, we used computer-aided methods to identify potential DNA polymerase inhibitors. The top four compounds, two drug compounds (canagliflozin and tepotinib) and two phytocompounds (rhein and taxifolin), were identified as highly promising LSDV DNA polymerase protein inhibitors. Molecular docking results reveal that canagliflozin poses the highest Glide Gscore (−9.858 kcal/mol) and the binding affinity (−45.68 kcal/mol) among the derived compounds. tepotinib, rhein, and taxifolin exhibited higher Glide Gscores and binding affinities than those of the control, Ivermectin B1a (−5.630 kcal/mol) (
Table 3). These compounds form hydrogen bonds, pi-pi interactions, salt bridges, and cation-pi bonds with the LSDV DNA polymerase protein.
The formation of hydrogen bonds between the ligand and protein is highly selective and specific [
80]. This depends on the spatial arrangement of the atoms in both the binding site of the protein and the ligand. The complementarity of these structures allows precise interactions to occur. The presence of hydrogen bonds stabilizes the ligand-protein complex, enhancing binding specificity [
81]. Our results showed that each of the selected compounds formed hydrogen bonds with specific active-site residues of the LSDV DNA polymerase protein (
Table 4), maintaining a specific distance. Taxifolin formed the highest number of hydrogen bonds. The positive control, ivermectin B1a, formed three hydrogen bonds. The hydrogen bond distance between Asn 659 and canagliflozin is 1.80 Å and 1.73 Å, respectively. The difference in the bond distance was due to the differences in the molecules within the same amino acid residue to which they bonded. The minimum bond distance between the ligand and specific amino acid residues indicates the maximum hydrogen bond energy.
Pi-Pi interactions are frequent in protein crystal structures, which help in the interaction of proteins and small molecules. The geometry of aromatic compounds and electrostatic interactions play a role in Pi-Pi interactions [
82]. Besides, Pi-Pi stacking is also involved in the binding energy of receptor ligands [
83]. From our findings, only Tepotinib formed three Pi-Pi interactions with the Tyr 477, Phe 499, and Tyr 663 amino acid residues, maintaining bond distances of 5.49, 4.89, and 5.38 Å, respectively.
Salt bridges are one of the strongest non-covalent interactions in nature. It plays crucial roles in protein folding, mediates protein-protein interactions, and facilitates molecular recognition [
84]. Our results suggest that canagliflozin, tepotinib, and rhein, but not taxifolin, form salt bridges. Canagliflozin forms a salt bridge with Lys 483 with a bond distance of 2.61 Å. Tepotinib forms a salt bridge with Glu 339 with a bond distance of 1.93 Å. Moreover, rhein formed a salt bridge with Lys 337 by maintaining a bond distance of 2.71 Å. The Positive control ivermectin B1a did not exhibit such bonding interactions.
Cation-pi bond interactions play a crucial role in protein stability and structure. It involves attractive forces between the aromatic ring of a small molecule and a positively charged cation in a protein [
85,
86,
87,
88]. In terms of energy, this interaction is comparable to or sometimes stronger than a hydrogen bond [
89]. Research suggests that cation-pi interactions are essential for protein-ligand recognition and have valuable implications for predicting drug-receptor interactions [
90]. Our study reveals that tepotinib and taxifolin from our selected compounds only formed cation-pi bonded interactions with Lys 483 having bond distances of 4.01 and 3.90 Å, respectively. In contrast, positive control did not form these bonds.
A similar binding energy pattern was observed for the non-bonded interactions (
Table 5). Coulomb energy is a form of electrostatic energy and a significant non-bonded energy found in protein-ligand complexes. Rhein has the highest coulomb energy (−77.35 kcal/mol), and ivermectin B1a possesses the lowest coulomb energy without any protein-ligand strains. Lipophilic energy is important for the evaluation of drug uptake and metabolism. Canagliflozin possesses the highest lipophilic energy (−25.21 kcal/mol) without strain. The Van der Waals interaction energies showed the opposite results. The Van der Waals interaction energy without receptor and ligand strains for the positive control was the highest (−63.90 kcal/mol) compared to the selected compounds.
The drug-like properties of the compounds satisfied Lipinski's Rule of Five, indicating their potential as drug candidates. Pharmacokinetics and toxicology offer important insights into the interactions between drug molecules in the human body [
91]. The properties of these compounds were calculated using ADMET analysis tools such as QikProp. These tools predict the results of various chemical and physical properties of drug candidates, including molecular weight (MW), surface area, hydrogen bond (HB) interactions, lipophilicity, and oral absorption rates in humans [
92]. Similarly, Qplog Po/w, QPlogS, Qplog HERG, and Human oral absorption were the major parameters used for the assessment of chemical compounds and for identifying the pharmacokinetic features of drugs (
Table 4).
In addition, most drugs on the market have a molecular mass between 200 and 600 Da, with the majority being less than 500 Daltons [
93]. To determine whether Canagliflozin, Tepotinib, rhein, and taxifolin are viable options, Lipinski's Rule of Five was considered. Lipinski's Rule of Five is a set of criteria used to evaluate a compound's drug-like behavior, including assessments such as the AMES test, Veber rule, and bioavailability [
94]. Canagliflozin, tepotinib, rhein, and taxifolin satisfied most of the requirements, specifically Lipinski's Rule of Five, ensuring their drug-likeliness behavior (i.e., AMES test, Veber rule with no violations, better bioavailability value, and other parameters), which could be used as drug candidates (
Table 1 and
Supplementary Table S1).
Molecular simulations are a powerful approach for understanding the stability and dynamics of protein-ligand complexes [
95,
96]. Where the Mean Square Deviation indicates the average atomic displacement, with higher RMSD values indicating structural deviation or flexibility and a lower value indicating structural stability. The Square Fluctuation quantifies the average atomic positional variance from their mean positions within a molecule, with higher RMSF values suggesting greater atomic fluctuation or flexibility, and lower RMSF values indicating less atomic fluctuation and greater stability. The radius of gyration in molecular dynamics is a measure of the compactness or spread of a molecule's structure, where a higher Rg indicates a more extended or less compact structure, and a lower Rg signifies a more compact or folded structure. The solvent-accessible surface area quantifies the surface area of a molecule that is accessible to solvent molecules, where a higher SASA value suggests a greater exposed surface area, and a lower SASA value indicates a reduced exposed surface area. These parameters collectively provide essential insights into the molecular behavior in simulations [
63,
97,
98].
The stability and dynamics of the compounds in protein-ligand complexes were studied. From the average mean value of the molecular dynamic simulation trajectory and analysis of the whole simulation data within 100 ns, it can be determined that the phytocompound rhein showed the highest stability in RMSD and RMSF. This fluctuated slightly in comparison with that of the positive control. The Rg and SASA values indicated less compactness and less exposure to the surface area compared to the positive control. It exhibited structural stability from 30 to 63 ns, whereas the rest of the simulation fluctuated. Taxifolin maintained its compactness in the RMSD. It showed stability at 10–89 ns; after that, it fluctuated from 1 to 10 ns and 91 to 100 ns, which is comparatively less than that of the positive control and gives a slightly higher RMSF value. The Rg and SASA values provided more compactness and a greater exposed surface area, respectively, than the other three compounds, but were marginally higher than those of the positive control. Tepotinib had the highest RMSD value among all compounds, indicating a lower level of compactness and even much higher compactness than the positive control, although this drug gave a minimum value in RMSF compared to the positive control and the rest of the compounds. It exhibited compactness in the range of 1–77 ns, and instability in the range of 78–100 ns. Canagliflozin showed the same RMSD score as the positive control. The RMSF of this compound was the highest among all the other compounds, including the positive control. Rg and SASA values followed the same trends. Canagliflozin exhibited compactness in 15 to 63 ns.
Canagliflozin and tepotinib are repurposed drugs, which means they have already been approved for use in treating specific diseases, and currently, there is interest in exploring their potential use for different diseases, including LSD. Canagliflozin (Invokana) is primarily used to the treatment of type 2 diabetes. It belongs to a class of drugs known as sodium-glucose co-transporter 2 (SGLT2) inhibitors. Canagliflozin is typically prescribed as a third-line treatment option after metformin (first-line medication) and other second-line alternatives in cases where these options are insufficient to effectively manage blood sugar levels [
99,
100,
101]. Tepotinib belongs to a class of medications known as kinase inhibitors. It works by blocking the action of an abnormal protein that signals cancer cells to multiply [
102,
103,
104].
Similarly, in this study, two specific phytochemicals, rhein and taxifolin, were identified as top-screened compounds with drug-like properties in inhibiting the macro domain of the chikungunya virus
[105], with a focus on finding phytochemical inhibitors (drug-like properties of the top-screened phytochemicals rhein and taxifolin) against a unique macro domain found in the conserved N-terminal region of the non-structural protein nsP3 of the chikungunya virus. Medicinal plants have long been used as a source of potential therapeutic compounds and may offer a promising approach to combat LSD. This study investigated the binding affinities of two phytocompounds (rhein and taxifolin) and two repurposed drugs (canagliflozin and tepotinib) to the active sites of the DNA polymerase protein LSDV039, which are important targets for antiviral drug development. Post-MM/GBSA analysis, which calculated the binding free energy, indicated that the complexes exhibited a better binding free energy at 100 ns than the pre-MM/GBSA binding free energy. This indicates that the binding of phytocompounds to target proteins became more favorable over time. This suggests that they have the potential to interact effectively with these proteins, possibly by inhibiting their activity and interfering with viral replication.
However, to develop these selected compounds as LSD-specific antivirals, it is crucial to understand their behavior in their natural hosts, confirm their ability to inhibit replication, and provide antiviral prophylaxis. Notably, the identified compounds (canagliflozin, tepotinib, rhein, and taxifolin) are publicly available, which facilitates their rapid utilization and further research. Further experimental studies are necessary to validate their efficacy as novel compounds against LSD.
Overall, LSD continues to pose significant challenges to the livestock industry, and more comprehensive control measures, including improved vaccination coverage and alternative treatment options, are required to mitigate its impact. This study emphasized the importance of pharmacokinetic properties, efficacy, and safety levels in determining therapeutic efficacy and releasing new drugs. Based on in-silico screening of natural compounds, that is, rhein and taxifolin, and repurposed drugs, that is, canagliflozin and tepotinib, combinatorial docking, molecular dynamic simulation, drug-likeliness analysis, and all other experimental data were used. Taxifolin and canagliflozin have been suggested as drug candidates for the treatment of LSD.