The analysis of the physicochemical properties of the ligands offers important insights into their potential bioactivity, pharmacokinetics, and overall drug-likeness. Molecular weight is a critical factor influencing absorption, distribution, metabolism, and excretion (ADME) characteristics. Generally, smaller molecules exhibit improved membrane permeability and absorption, whereas larger molecules often face challenges related to poor solubility and reduced bioavailability. Among the ligands evaluated, Curcuminol, with a molecular weight of 566.51 g/mol, is the heaviest, which may negatively impact its absorption and oral bioavailability. Conversely, compounds such as Caffeic acid and Sinapic acid, which possess lower molecular weights, are likely to demonstrate better absorption profiles. Most ligands fall within the typical molecular weight range for drug-like compounds (200–500 g/mol), with Curcuminol as the notable exception, suggesting potential limitations for its oral administration [19]. These findings align with established knowledge that smaller molecules generally achieve superior membrane permeability and absorption, while larger compounds may be hindered by poor solubility and permeability issues [20]. The number of rotatable bonds influences the flexibility of a molecule, which can have both positive and negative implications for binding to biological targets. A higher number of rotatable bonds can increase the conformational diversity of a molecule, potentially improving its ability to bind to a range of targets. However, excessive flexibility could lead to instability and reduced specificity. Curcuminol, with 13 rotatable bonds, is the most flexible molecule in the set, which could increase its potential for binding to various targets but might also compromise its stability. On the other hand, compounds like Curdione (1 RB) and Caffeic acid (2 RB) are more rigid, which may offer better specificity for well-defined targets, although it could limit their binding diversity [21, 22]. The hydrogen bonding capacity of a compound influences its solubility and interaction with biomolecules. Compounds with a high number of HBA and donors (HBD) can exhibit strong interactions with biological targets, but excessive hydrogen bonding may reduce membrane permeability [23]. Curcuminol and isorhamnetin, with the highest number of HBA (12 and 7, respectively), may have strong interactions with the polar regions of biological targets. However, such compounds may also face challenges in terms of membrane permeability. Curdione, with only 2 HBA and no HBD, is less likely to engage in hydrogen bonding, potentially reducing its bioactivity but increasing its specificity. Ribavirin, with a balanced hydrogen bonding profile, could be seen as an optimized structure for maintaining bioactivity while potentially having favorable ADME properties [24].
The physicochemical and ADMET profiles of the tested ligands reveal significant insights into their potential as antiviral agents against the PVYN-Egypt strain. Curcuminol, despite its high molecular weight and flexibility, violates Lipinski’s rule and shows low gastrointestinal absorption, which may limit its oral bioavailability and clinical application, consistent with previous observations that large, flexible molecules often face bioavailability challenges [25]. In contrast, ligands like curcumin, bisdemethoxycurcumin, and isorhamnetin exhibit favorable drug-like properties with high GI absorption and compliance with Lipinski’s rule, highlighting their potential as orally active candidates [26]. The toxicity profiles further refine candidate selection, as compounds such as curdione and curcuminol show potential toxicity, warranting caution in further development [27].
Virtual screening results demonstrated strong binding affinities of curcuminoids and flavonoids to key viral proteins, with isorhamnetin and bisdemethoxycurcumin showing particularly promising interactions with the P1 protease and helper component proteinase, suggesting their role as potent inhibitors in viral replication and transmission processes [28]. The control drug ribavirin, despite moderate binding affinity, confirms the validity of the docking approach and underscores the improved efficacy of these natural ligands in targeting PVY proteins [29]. Detailed interaction analyses reveal that hydrogen bonding, Pi-Alkyl, and salt bridge formations are critical in ligand stabilization within the active sites, aligning with prior studies that emphasize the importance of multi-modal interactions for effective protease inhibition [30, 31].
Moreover, the inability of most ligands to permeate the blood-brain barrier aligns with the peripheral targeting profile required for plant viral infections, reducing the risk of central nervous system side effects [32]. However, bisdemethoxycurcumin and curdione’s BBB permeability may open avenues for broader antiviral applications but require further toxicity assessment. Collectively, these findings suggest that the curcuminoids and flavonoids under study merit further in vitro and in vivo validation as potential antiviral agents, especially given their multifaceted interaction profiles and favorable ADMET characteristics, supporting the growing body of evidence on natural products as valuable scaffolds in antiviral drug discovery [33, 34].
GI absorption is a critical determinant of oral bioavailability, as compounds with high GI absorption are generally absorbed more efficiently from the gastrointestinal tract. In this dataset, most compounds are predicted to have high GI absorption, which is a favorable characteristic for potential drug development. However, curcuminol and ribavirin, with predicted low GI absorption, could face challenges related to oral bioavailability. These compounds might experience poor intestinal permeability or high first-pass metabolism, which can limit their therapeutic effectiveness [35].
The ability of a compound to cross the blood-brain barrier (BBB) is critical for drugs targeting the CNS. In this dataset, bisdemethoxycurcumin and Curdione are predicted to be BBB permeant, suggesting their potential for treating CNS-related diseases. On the other hand, most other ligands, such as Caffeic acid, curcumin, and demethoxycurcumin, are not predicted to cross the BBB, suggesting they may be more effective for treating peripheral diseases [36].
Cytochrome P450 (CYP) enzymes are essential for the metabolism of numerous drugs, and their inhibition can result in potentially harmful drug-drug interactions. Bisdemethoxycurcumin and demethoxycurcumin are predicted to inhibit multiple CYP isoforms, raising concerns about significant interactions when co-administered with other therapeutics. Likewise, curcumin, isorhamnetin, and ribavirin inhibit CYP3A4, a key enzyme involved in the metabolism of a wide variety of drugs [37, 38]. Therefore, the potential for CYP enzyme inhibition must be carefully evaluated during drug development, especially when these compounds are intended for use alongside other medications that depend on CYP-mediated metabolism [39].
Lipinski’s "Rule of Five" provides a widely accepted framework for assessing the drug-likeness of compounds, with adherence suggesting favorable oral bioavailability. In this analysis, most compounds, with the exception of curcuminol, conform to Lipinski’s criteria, indicating promising drug-like properties and a higher likelihood of oral absorption. Curcuminol, however, due to its large molecular weight and elevated hydrogen bonding capacity, violates Lipinski’s rule, which may pose challenges related to solubility and oral bioavailability [23].
Toxicity is a critical consideration in drug development, and compounds that are classified as toxic may pose risks for patient safety. In this dataset, curcuminol and Curdione are classified as toxic, suggesting that despite their promising ADMET properties, they may exhibit undesirable side effects that limit their therapeutic use. Conversely, compounds such as ribavirin, which are classified as non-toxic, may present fewer safety concerns, making them more suitable for further development. Therefore, toxicity must be considered alongside efficacy and ADMET properties to ensure the safe development of these compounds as potential therapeutics. In a summary, the ligands in this virtual screening dataset demonstrate varied ADMET properties, which can be linked to their pharmacological potential [40].
Bisdemethoxycurcumin, Curdione, and demethoxycurcumin show high GI absorption and favorable pharmacokinetic properties, though some of them are predicted to inhibit multiple CYP enzymes, raising concerns about drug-drug interactions. Curcumin, isorhamnetin, and ribavirin also show good absorption and drug-like properties but may exhibit interactions with CYP3A4. On the other hand, curcuminol exhibits poor GI absorption, does not comply with Lipinski's rule, and is classified as toxic, indicating that it may require substantial modifications to improve its drug-like properties and safety profile. Caffeic acid and Sinapic acid have high GI absorption and do not inhibit CYP enzymes, making them interesting candidates for further studies targeting peripheral diseases. The strong binding of Isorhamnetin to the P1 protease, with a binding affinity of -7.1 kcal/mol, highlights its potential as a lead compound for antiviral drug design targeting this key enzyme in the PVYN-Egypt strain. Curcuminoids like Bisdemethoxycurcumin and Curcumin also demonstrate promising binding affinities, which is consistent with their previously reported antiviral properties [10, 11]. Curcuminoids have shown broad-spectrum activity against various viruses, including plant viruses, suggesting that they could be used to inhibit P1 protease activity in PVY [16, 41].
Bisdemethoxycurcumin, Demethoxycurcumin, and Isorhamnetin show the most favorable binding to HCPro, with binding affinities of -6.4 kcal/mol. This is significant because HCPro plays an essential role in viral transmission and replication. The observed high binding affinities of curcuminoids and flavonoids like Isorhamnetin align with studies that report their broad-spectrum antiviral properties, particularly against plant viruses [16, 41, 42]. This suggests that these compounds may effectively inhibit HCPro activity, thereby reducing PVY replication and transmission. Curcumin and Isorhamnetin demonstrate the strongest binding affinities for the PVYN-Egypt coat protein at -7.4 kcal/mol, which could potentially interfere with viral capsid formation and RNA encapsidation, crucial steps for the virus's stability and transmission. The observed results are consistent with previous findings that curcuminoids can disrupt viral capsid formation [11, 41, 43], and Isorhamnetin, a flavonoid, has shown promising antiviral activity in several studies [44]. This suggests that both compounds could play a significant role in disrupting the viral lifecycle by targeting the coat protein.
The RMSD (Root Mean Square Deviation) values of 0.0 observed across all ligand–protein complexes may initially suggest high accuracy in the docking predictions, indicating that the predicted binding poses are closely aligned with the reference conformations. However, such uniformly low RMSD values can also reflect methodological limitations, particularly when rigid docking protocols are employed. In this study, ligands were docked back into their native conformations without allowing for conformational flexibility or sampling of alternative binding modes. While this approach confirms pose reproducibility, it does not fully account for the dynamic nature of ligand–protein interactions in realistic biological environments [45, 46]. We acknowledge this as a limitation of the current methodology and recommend that future studies incorporate flexible docking protocols or molecular dynamics simulations to enhance the predictive accuracy [47]. Accordingly, we have conducted a limited re-evaluation using flexible docking in this revision to support the reliability of our findings.
Overall, the ligand binding affinities obtained in this study identify promising candidates for the development of antiviral agents against the PVYN-Egypt strain. Curcumin, Isorhamnetin, and Bisdemethoxycurcumin demonstrated strong binding across all three target proteins, underscoring their potential as broad-spectrum antiviral compounds. These findings are consistent with previous studies reporting the antiviral efficacy of curcumin derivatives and flavonoids, which have been shown to inhibit viral replication and disrupt key protein functions [10, 11, 44]. Additionally, Ribavirin, employed here as a standard antiviral control, exhibited comparatively strong binding affinities, reaffirming its established antiviral activity [24]. The docking results emphasize that curcumin derivatives and flavonoids—particularly Bisdemethoxycurcumin, Curcumin, and Isorhamnetin—possess strong affinity for critical proteins of the PVYN-Egypt strain. This aligns with the well-documented antiviral properties of curcuminoids, which are known to suppress viral proliferation by targeting essential viral enzymes and protein functions [45–47].
The high binding affinity of Isorhamnetin, a flavonoid, further supports the potential of flavonoids in antiviral drug discovery [48, 49]. Additionally, Ribavirin, a well-known antiviral drug, exhibited docking scores in the moderate range across all three proteins. While its binding affinity was lower than that of curcumin derivatives and Isorhamnetin, Ribavirin's established antiviral activity [50–53], makes it a valuable control in this study. Bisdemethoxycurcumin, curcumin, and isorhamnetin show strong binding interactions with PVY P1 protease, particularly through a combination of Pi-Alkyl, Pi-Anion, hydrogen bonds, and Pi-Cation interactions. These ligands seem to have a well-established binding profile with significant contributions from both hydrophobic and polar interactions, which may contribute to their high binding affinity [54–56].
The salt bridge formation observed with GLU (J:2974) in isorhamnetin could indicate a more stable binding compared to other ligands, potentially making isorhamnetin a better candidate for inhibition [57]. Hydrogen bonding is a critical component of ligand-receptor interactions and significantly contributes to binding affinity and specificity. The ligands bisdemethoxycurcumin, curcumin, and isorhamnetin form multiple hydrogen bonds with residues like GLN, SER, and ARG, which help stabilize the ligand-receptor complex [48–50]. Ribavirin also forms hydrogen bonds, but the presence of unfavorable acceptor-acceptor interactions suggests that its binding affinity might be weaker or less stable compared to the other ligands. The unfavorable acceptor-acceptor interactions in ribavirin (with MET (M:2895) and ASN (F:2971) may reduce its binding efficiency and increase the likelihood of off-target effects or weaker inhibition. This suggests that ribavirin may require structural modifications to improve its binding profile and reduce unfavorable interactions [61, 62].
Based on the detailed binding analysis, bisdemethoxycurcumin, curcumin, and isorhamnetin are strong candidates for further drug development targeting PVY P1 protease, given their favorable interaction profiles and strong binding affinities. These ligands could be optimized through structural modifications to improve potency, selectivity, and minimize toxicity [63]. Ribavirin, while having some beneficial interactions, may need structural optimization to overcome the unfavorable interactions and enhance its overall binding stability. Based on the diverse and stable interactions with key amino acids in the protease, curcumin, bisdemethoxycurcumin, and curcuminol appear to be strong candidates for further optimization and development as potent protease inhibitors. These ligands display a wide range of favorable interactions, indicating that they could effectively bind to and inhibit the PVY helper component proteinase.
The variety of interactions, including Pi-Alkyl, Pi-Sulfur, Pi-Pi T-Shaped, and hydrogen bonding, suggests that these compounds are engaging in both hydrophobic and electrostatic interactions that contribute to their high binding affinity [64, 65]. In particular, the Pi-Sulfur interactions with CYS (A:37) observed in curcumin and curcuminol highlight the potential for unique interaction profiles, which may improve the binding affinity and specificity for the protease. These interactions may provide these ligands with an advantage in terms of potency and selectivity for the target enzyme [66]. Moreover, the diversity of bonding types seen in bisdemethoxycurcumin and curcuminol-such as Pi-Alkyl, Pi-Pi T-Shaped, and hydrogen bonds-suggests that these compounds could engage the protease in a stable, multifaceted manner, likely enhancing their inhibitory potential [67]. While demethoxycurcumin and ribavirin show more limited interaction profiles, they still exhibit significant binding interactions, suggesting that their structures could be further optimized for better protease inhibition. The more limited interaction profiles may indicate that these compounds are less stable in the active site compared to curcumin derivatives, potentially reducing their potency as inhibitors. Structural modifications could improve these interactions and enhance the binding efficiency [56, 68].
Isorhamnetin, despite forming several hydrogen bonds, exhibits unfavorable acceptor-acceptor interactions that may reduce its overall binding stability, thereby limiting its effectiveness as a protease inhibitor. These unfavorable interactions, especially with residues such as SER (A:42) and LEU (A:53), may result in steric hindrance or decreased affinity. Therefore, isorhamnetin could benefit from structural modifications aimed at overcoming these unfavorable interactions to enhance binding efficiency and improve its inhibitory potency [69, 70].
Curcumin, bisdemethoxycurcumin, and curcuminol generally emerge as promising candidates for further development as protease inhibitors targeting the PVY helper component proteinase. Although isorhamnetin and ribavirin exhibit some potential, they may require structural optimization to mitigate unfavorable interactions and enhance their overall inhibitory efficacy [42, 70]. Analysis of the binding interactions with the PVY coat protein highlights curcumin as the strongest candidate for development as a protease inhibitor against this target. The compound’s robust and stable binding affinity is likely driven by a diverse array of interactions, including Pi-Cation, Pi-Alkyl, and hydrogen bonds [67, 71].
Notably, curcumin forms Pi-Cation interactions with ARG (A:55) and Pi-Alkyl interactions with HIS (A:162), indicating engagement through both hydrophobic and aromatic contacts, which contribute to a stable binding conformation. Furthermore, hydrogen bonds with residues such as GLN (A:140), GLY (A:61), and LYS (A:60) further reinforce the stability of the curcumin-protein complex, enhancing its affinity for the PVY coat protein [47, 72]. However, an unfavorable donor-donor interaction with GLN (A:140) suggests possible steric clashes or electrostatic repulsion that could detract from curcumin’s binding efficiency at this site. Addressing this issue in future ligand optimization—potentially through structural modifications to reduce steric hindrance or electrostatic conflict—may improve binding potency [73].
Isorhamnetin, while promising as a ligand for the PVY coat protein, forms fewer interaction types compared to curcumin. It mainly establishes hydrogen bonds with key residues such as ARG (A:55) and ASP (A:91), suggesting effective binding, though possibly with lower affinity due to the more limited interaction profile [74, 75].
Structural modifications aimed at enhancing the interaction profile—such as increasing hydrophobic contacts or introducing additional binding moieties—could improve the inhibitory potency and specificity of these compounds. Ribavirin forms stable hydrogen bonds with residues including ASN (A:82), THR (A:86), TRP (A:118), and GLY (A:130). Additionally, a Pi-lone pair interaction with THR (A:86) may contribute to binding specificity; however, ribavirin’s overall interaction profile is simpler compared to curcumin’s, which limits its binding affinity. Consequently, ribavirin may benefit from optimization to enhance its Pi-lone pair interactions and overall binding conformation [76]. Among the compounds studied, curcumin stands out as the most promising candidate for further development, owing to its diverse and stabilizing interactions, particularly the Pi-Cation and Pi-Alkyl interactions [77, 78]. Isorhamnetin also shows potential for optimization, displaying favorable hydrogen bonding; yet, its affinity may be comparatively lower due to the narrower range of interaction types. Ribavirin, while forming stable interactions, is likely to require structural refinement to improve both its binding affinity and specificity for the PVY coat protein.
Limitations
While the docking results indicate promising binding interactions with RMSD values close to zero, it is important to recognize the inherent limitations of molecular docking simulations. Such ideal RMSD values may reflect methodological constraints, including limited sampling of ligand-protein conformations and reliance on static protein structures. Moreover, docking predictions are inherently computational and provide only a preliminary understanding of molecular interactions. Therefore, these findings should be interpreted as hypotheses that require subsequent experimental validation through in vitro and in vivo studies. Additionally, the complex dynamics of biological systems and potential off-target effects are not fully captured in silico, underscoring the need for cautious interpretation and further comprehensive investigations.
Future Directions
Building on the promising in silico findings of this study, future research should focus on experimental validation through in vitro assays to confirm the antiviral activity of the identified turmeric-derived compounds against PVY. Subsequently, in vivo studies will be essential to evaluate their efficacy and safety in plant models. Additionally, structural modifications of these bioactive molecules could be explored to enhance their binding affinity, bioavailability, and specificity. Investigating potential synergistic effects through combination therapies with existing antiviral agents may also provide more effective control strategies. Ultimately, integrating computational predictions with empirical data will facilitate the development of novel, targeted antiviral treatments for crop protection.