Skip to main content

Computational investigation of the anticancer potential of Sorghum bicolor and Setaria italica phytochemicals against dihydrofolate reductase (DHFR) enzyme

Abstract

Breast and prostate cancer holds the position of foremost contributors to mortality. Dietary therapies for accompanied by medication are widely recognized as a potential method to successfully tackle cancer. Millet grains are the most ancient food, a perfect combination of proteins, carbohydrates, fiber, macronutrients, micronutrients, and vitamins. This study aims to examine the anticancer potential of Sorghum bicolor (Sorghum) and Setaria italica (Foxtail) phytochemicals. The 50 phytochemicals of sorghum and foxtail millets were retrieved through a literature survey and docked to the Dihydrofolate reductase (DHFR), an enzyme essential for cell growth and proliferation. The top-scoring phytochemicals were filtered and further investigated with active-site residue interaction, drug-likeness, and pharmacokinetics analysis. The ligand stability with the DHFR was evaluated through density functional theory (DFT) based HOMO and LUMO calculations. The results show that caffeic acid, ferulic acid, hesperetic acid, stigmasterol, Cis-p-Coumaric acid, and luteolinidin attained greater stability within the active site of DHFR. These phytochemicals showed a docking score of − 6.4 kcal/mol, − 6.4 kcal/mol, − 6.1 kcal/mol, − 6.4 kcal/mol, − 5.4 kcal/mol, and − 6.7 kcal/mol with DHFR (PDB ID:1BOZ) and flutamide and capecitabine have docking score of − 7.5 and − 8.1 for 1BOZ and − 7.4 and − 7.1 with DHFR (PDB ID:1OHK) respectively. The dynamic interaction at the molecular level validated the stability of these phytochemicals against both DHFR target proteins along with excellent drug-likeness and pharmacokinetic properties. However, the current findings were proven and validated through in-silico experiments to validate above identified phytochemicals as DHFR inhibitors, so millets are used as therapeutics for breast and prostate cancer.

Introduction

The food security is among the looming challenges for the twenty-first century which is flattering precariously through population growth and climate change. The United Nations forecasts a nine billion population growth by 2050 globally. The continents are witnessing a day by day increment in hunger, drought, and climate change conditions jeopardizing food security. The Norman Borlaug-led miracle seed green revolution nowadays needs another agricultural revolution. The Millets contribute as a reliable food source option for combating global hunger in the context of rapidly shifting climate conditions. The Food and Agricultural Organization of the United Nations defined Millet is as a cumulative term denoting to a number of small-seeded annual grasses cultivated as grain crops predominantly on marginal dry lands of temperate, tropical, and subtropical regions. Millets come into view as climate resilient and superior photosynthetic efficiencies crops over staple ones (Mehta et al. 2024).

During its 75th session in March 2021, the United Nations General Assembly officially designated the year 2023 as the International Year of Millets (IYM 2023) to shed light on millets nutritional and ecological importance. The extended purpose of IYM2023 is to draw policymakers, researcher and layman attention to millets potential for food security. Millets represent a viable source for nations seeking to enhance their self-sufficiency and diminish their dependence on imported cereal grains (Poshadri et al. 2023). Millets, staple crops in semiarid tropics, are thriving due to their high nutrient content and tolerance to extreme weather conditions. With concerns about lifestyle diseases arising due to refined diets, modern consumers are increasingly choosing millet as an alternative to wheat and rice, especially during COVID-19.

Millets, poaceae family members are ancient staple sustainable cereal crops foremost domesticated and cultivated by prehistoric humans Millets and other cereal grains are major food source for economically challenged parts of Asia and Africa, accounting for 75% of total calorie consumption. These regions produce 14.2 and 12.4 million tonnes of millet on average per year (Belton and Taylor 2004). India is the leading worldwide millet producer, accounting for almost 80% of overall millet output (Farsund et al. 2015).

Based on economic significance, cultivation, and consumption share, millets are grouped into two categories minor millets and major millets. The Pearl Millet (Bajra) (Pennisetum glaucum (L.), Little Millet (Sorghum) (Panicum sumatrense L), Amarnath (Proso Millet) (Amaranthus cruentus), Kodo Millet (Paspalums crobiculatum L.), Finger Millet (Ragi) (Eleusine coracana L), Barnyard Millet (Sawan) (Echinochola frumentacea L), Foxtail Millet (Kakun) (Setaria italic L), Teff Millet (Abyssinian lovegrass) (Eragrostis tef), and Browntop Millet (Panicum ramose L) are imperative millets cultivated globally. Millets are disseminated as Nutri cereals or superfoods due to their richness in proteins, dietary fibers, antioxidants, vitamins bioactive phytochemicals, polyphenols and minerals. The genomic diversity of millets helps to thrive in diverse agricultural domains (Maharajan et al. 2024).

Sorghum (Sorghum bicolor) is a staple energy source for resource-limited inhabitants of Globe. It is a genotypic and phenotypic diverse, climate-resilient, drought-tolerant and multipurpose (food, fiber, fuel, and feed) crop (Bakari et al. 2023) (Okoye et al. 2014).

Foxtail millet (Seteria italica) is the oldest cultivated crop whose bio-functional phytochemicals show antihyperglycemic, anticancerous, antioxidant, anti-inflammatory, antihyperlipidemic medical prospects properties along with substantial proportions of carbohydrate, protein, fibers, and minerals (Kalsi and Bhasin, 2023).

Cancer is a cumulative disease term in which cells start uncontrolled cell division followed by invasion and spread to other cells and tissues. Sometimes disruptions in gene expression can cause erratic cellular proliferation, leading to the abnormal increase of cells. Cancer has been an ongoing public health concern throughout history, appearing as a top cause of death in advanced and developing countries worldwide. Current data analytics indicate a projected increase, with an estimated 23.6 million additional cancer cases per year by 2030 (Bray et al. 2018).

The phytochemicals are biologically active nutritive of plants. The alkaloids, phenolic, nitrogen-containing, organosulfur, and carotenoids are major groups of bioactive phytochemicals. The phytochemicals are potential agents to combat diabetes and cardiovascular diseases The FMBP, class III peroxidase protein extracted from foxtail millet was shown anti-colon cancer properties in-vitro and in-vivo experiments on mice. The FMBP arrests cancerous colon epithelial cells in the G1 phase whereas normal epithelial cells remain unaffected. The FMBP has shown effects on epithelial mesenchymal transition by down regulation of JAK1 signalling pathway followed by accumulation of more ROS in cancerous colon epithelial cells (Shan et al. 2014).

Another similar research proved that flavonoids, carotenoids and phenolics extract of foxtail millet arrest breast cancer cell lines in G2/M phase along with DNA fragmentation. Thus extracts show antiproliferative potential. The Kodo and pearl millet bran extracts showed anti-cancerous property in HT-29 cancer cell lines. The research study concluded that phytochemicals extracted from Proso and barnyard millets induce apoptosis and cell proliferation in the HT-29 human colon cancer cell line and arrest cancerous cells in Go/G1 phase (Ramadoss and Sivalingam, 2020) (Gupta et al. 2023).

Phytate phytochemical extracted from millets shows antiproliferative activity and downregulates NF-kB, P13K/MAPK signaling pathways in breast cancer cell lines (Shukla et al. 2023) Several systematic and review articles concluded that millets, Wonder grains encompass nutritional, health-promoting, anti-inflammatory and anticancerous properties. Millets are accessible functional, low-cost foods that gaining the attention of industry to common people (Samtiya et al. 2022) (Nani et al. 2015).

Growth factor receptors (GFRs) are proteins that play important roles in many aspects of cell behavior, including tumor development, metastasis, angiogenesis, cell survival, apoptosis, cell motility, cellular differentiation, organ creation, neovascularization, and chemotherapy resistance(Tiash and Chowdhury 2015). These receptors are activated when they bind to certain ligands, commonly known as growth factors (Ornitz and Itoh 2015).

DHFR is a ubiquitous enzyme found in all organisms that catalyzes the conversion of 7,8-dihydrofolate (DHF) to 5,6,7,8-tetrahydrofolate (THF). This process requires a precise hydride transfer from the NADPH cofactor to the pterin ring’s C6 atom, followed by protonation at N5. DHFR plays a vital role in controlling cellular reservoirs of THF and its derivatives, which are required for the synthesis of purines and thymidylate, which are required for cell proliferation and growth(Schnell et al. 2004). Notably, this enzyme has received a lot of attention as a potential target for anticancer medicines and antibiotics(Chen et al. 2009). Anticancer medication development includes strategically disrupting certain proteins or pathways that are intricately engaged in cancer growth. As stated in (Siegel et al. 2023), breast cancer is one of the main causes of cancer-related death in women, while prostate cancer is similarly prevalent in men. Research has shown that DHFR is overexpressed in certain kinds of breast cancer. Furthermore, studies have revealed that inhibiting DHFR in prostate cancer cells may cause cell cycle arrest and death, thus slowing tumour development. Visual presentation of our work is shown below in Fig. 1

Fig. 1
figure 1

Graphical abstract of research

Materials and methods

Identification of the Plant’s compounds, and their preparation

A thorough review of the available literature resulted in the identification of around fifty unique plant millet chemicals within Sorghum bicolor (Sorghum) (Kauther et al. 2020) (de Morais Cardoso et al. 2017) (Shih et al. 2007) and Setaria italica (Foxtail) (Sharma et al. 2021) (Zhang et al. 2021) (Nivedita et al. 2022). The 2D/3D structures of these identified phytocompounds were obtained in SDF format from the PubChem (https://pubchem.ncbi.nlm.nih.gov/) chemical database, resulting in a complete library of chemical compounds as depicted in.

Fig. 2
figure 2

Figure 3D Structures of 50 phytochemicals used in this study. (1) 3,4-Dihydroxybenzoic acid, (2) Vanillic Acid, (3) Truxillic acid, (4) Naringenin, (5) Luteolinidin, (6) Ferulic acid, (7) Squalene,(8) cis-p-Coumaric acid, (9)Stigmasterol, (10) Kaempferol,(11) Ergostanol,(12) trans-p-Coumaric acid (13) 4-glucoside,Phytic acid,(14) 4-Hydroxybenzaldehyde,(15) Protocatechualdehyde,(16) Luteoforol,(17) Taxifolin,(18) Peonidin,(19) Apigeninidin,(20) CaffeicAcid,(21) Quercetin,(22) Sinapaldehyde,(23) Zeaxanthin,(24) Glucosyringic acid,(25) 2-Methylisocitric acid,(26) Fisetinidin,(27) 4-Hydroxybenzoic acid, (28) Syringic acid, (29) Campesterol, (30) Pelargonidin, (31) Apiforol,(32) p-Coumaric acid, (33) Isoferulic acid, (34) Vitexin, (35) Isoquercetin, (36) Lutein,(37) (2S)-eriodictoyl-7-O-beta-D-glucopyranosiduronic acid, (38) Di-p-coumaroylspermidine, (39) 2-Methylcitric acid, (40) Homocitric acid, (41) Beta-Sitosterol,(42) 2-(4-Hydroxyphenyl)chromenylium-5,7-diol, (43) Cinnamic acid,(44) Sinapinic acid,(45) cis-Sinapic acid, (46) Apigenin,(47) Rutin,(48) Ellagic Acid, (49) 3,3′,4′,7-Tetrahydroxyflavylium, (50) N-caffeoyl-N'-feruloyl spermidine

Fig. 3
figure 3

2D interaction of the selected phytochemicals with 1BOZ and 1OHK

Figure 2Furthermore, two established medications, flutamide linked to prostate cancer and capecitabine to breast cancer, were obtained from the PubChem database with unique identities PubChem ID CID_3397 and CID_60953, respectively, for comparative comparison with the phytochemicals. These medications were also obtained in 3D SDF format. The ligands were converted into the PDB format using the PyRx Software (https://pyrx.sourceforge.io/downloads), making them suitable for further docking studies. This was followed by the construction of torsion requirements required for correct binding, which was accomplished through the use of AutoDock 4.2.6 parameters.

Source of protein crystal structures and their preparation

The protein X-Ray structures of DHFR a ubiquitous enzyme, were obtained from the protein data bank repository (http://www.rcsg.org). The unique identification numbers are PDB ID:1BOZ (Gangjee et al. 1998) and PDB ID: 1OHK(Cody et al. 1997).An array of preliminary processes were performed on the protein structures using AutoDock 4.2.6. This included removing water molecules and adding polar hydrogen atoms, as well as Kollman’s charges. Following that, the changed protein structures were stored in the PDBQT format to ensure compatibility with further docking investigations.

Molecular docking studies and virtual screening

Following the careful preparation of proteins and ligands, a molecular docking method was carried out using Autodock 4.2.6.. Specific dimensions were defined along the grid box of 1OHK was used with center_x = 1.079, centre_y = − 12.26, centre_z = 37.468 and 1BOZ with center_x = 24.397, centre_y = 12.818, centre_z = − 1.563, and the size_x = 20, size_y = 22, size_z = 20 for both receptor proteins respectively, while keeping a resolution of 0.500, to promote successful binding interactions.. The grid layout was precisely centered to maximize the creation of favourable docking conformations (Mendie and Hemalatha 2022). This grid setup was then saved as a (.gpf) file and executed by Autogrid.

Following that, docking calculations were started using the Lamarckian genetic method, with a standard of ten separate runs. The docking computations that resulted were saved in a (.dpf) file. Following the conclusion of the Autodock computations, the definitive docked results were obtained and stored in a (.dlg) file format. Tables 1 and 2 provides useful information, such as binding residues, binding energy, and inhibition constants. The structures representing the interactions between ligands and proteins were visualized using Discovery Studio (https://discover.3ds.com/discovery-studio-visualizer-download) to visually examine these interactions.

Table 1 Binding parameters between ligands and target protein 1BOZ
Table 2 Binding parameters between ligands and target protein 1OHK

In silico ADME analysis

The chemical compound’s drug-likeness, physicochemical qualities, and pharmacokinetic properties were comprehensively evaluated. This assessment included key pharmacokinetic characteristics such as absorption, distribution, metabolism, and excretion, which were examined in the ligands using the SwissADME platform (http://www.swissadme.ch/index.php). This detailed study provides significant insights that may be used to guide and improve the medication research and development process.

Single-point calculation to determine the reactivity and inertness of compounds

The fundamental ideas of density functional theory (DFT) are derived from the Hohenberg–Kohn theorem, which states that the energy of a system may be described as a function that is entirely dependent on the electron density. To investigate this, chosen chemical components derived from Sorghum bicolor (Sorghum) and Setaria italica (Foxtail) were meticulously geometrically optimized in Maestro using the Jaguar fast engine. The functional set used in this optimization procedure was Becke’s three-parameter exchange–correlation functional hybrid (B3LYP), which was supplemented with a 6-31G basis set, resulting in DFT as the preferred level of theory.

Following this optimization, the investigation focused on the geometrically improved compounds surface characteristics. This involved computing molecular orbitals in the gaseous phase, which provided important insights into the electrical structure and behavior of the compounds.

Results

Docking and virtual screening analysis

That growth factor receptors play in the start and progression of cancer, it is critical to develop therapeutic candidates that can effectively inhibit growth factors at receptor sites while causing minimum or insignificant effects elsewhere. The intensity of the interaction between ligands and receptors is governed by binding affinity, which is measured by binding energy. Notably, lower binding energy indicates a stronger binding affinity (McConkey et al. 2002).

To achieve this goal, 50 bioactive chemicals were subjected to docking interactions with the DHFR protein in a thorough virtual screening approach. The Flutamide (for Prostate cancer) and Capecitabine (for breast cancer) is docked as positive control to the DHFR receptor protein. The phytochemicals were sorted out on the basis of binding energy. The post-docking investigation were performed for binding energy (in kcal/mol), the number of hydrogen bonds, the inhibitory constant (expressed in M/nM), and the amino acids intricately involved in hydrogen bonding interactions between phytochemical ligands and target protein.

At the end of the molecular docking analysis this study identifies the following compounds with high binding affinity with DHFR; caffeic Acid, ferulic Acid, hesperetic Acid, stigmasterol, cis-p-coumaric acid, and luteolinidin. These compounds showed a docking score of − 6.4 kcal/mol, − 6.4 kcal/mol, − 6.1 kcal/mol, − 6.4 kcal/mol, − 5.4 kcal/mol, and − 6.7 kcal/mol with Dihydrofolate reductase and flutamide and capecitabine have docking score of − 7.5 and − 8.1 for 1BOZ and − 7.4 and − 7.1 for 1OHK respectively.

Pharmacokinetic and drug likeness screening of phytochemical

The path to drug discovery is inherently difficult. The financial investment necessary for drug research is significant, and the potential of developing medicinal molecules with low toxicity is particularly difficult. It is worth noting that only a small percentage of prospective lead compounds get approved by the FDA (Wishart 2007). A drug lead can be terminated at any point of the development process owing to a variety of issues such as inefficacy, unwanted effects, toxicity concerns, poor absorption, or inefficient clearance. Unfortunately, it has been discovered that the more promising a medication lead seems, the more expensive it is to cease its research. Prediction or modelling of ADME (absorption, distribution, metabolism, and excretion) is one potential route in this complex process. ADME data provides useful information on how medicine will interact with the body. While a pharmacologically promising lead may demonstrate great performance in vitro, unfavourable ADME data typically result in the lead's termination. Robust ADME prediction, in essence, serves as a vital filter, impacting the development route of possible drug candidates.

The results reported in Table 3 offer a complete summary of the chemical compounds drug-likeness and pharmacokinetic properties. The examination of physicochemical parameters, which include important chemical and physical factors such as molecular weight, hydrogen bond donors, hydrogen bond acceptors, and octanol/water partition coefficient (QPlogPo/w), helps to determine the drug-likeness of the compounds. The calculated molecular weights of the plant compounds vary from 164.20 to 412, with hydrogen bond donors and acceptors ranging from 1 to 4 and 3 to 8, respectively. Furthermore, the octanol/water partition coefficients range between 0.79 and 5.08, which are all within acceptable limits.

Table 3 Lipinski’s rule of five and radar chart of the selected phytochemicals

Significantly, all of these compounds fulfil more than two of Lipinski's rule of 5 (RO5) qualifying requirements (Lipinski et al. 2001), as they have a maximum of 5 hydrogen bond donors, a maximum of 5 hydrogen bond acceptors, a molecular weight less than 500 kDa, and log P values less than 5. As a result, all the crystallized ligand that binds to the proteins exhibit drug-like properties.

This assumption is supported further by Human Oral Availability (HOA) testing, which demonstrates that these drugs have medium to high availability. With the exception of Stigmasterol, all ligands had high GIA values. It is worth noting that crucial therapeutic drugs usually operate as p-glycoprotein substrates, which are frequently linked with decreased drug absorption, permeability, oral bioavailability, and total drug retention. P-glycoproteins are highly expressed in cancer cells, posing a substantial challenge to cancer therapy. This overexpression causes drug efflux, which significantly reduces the efficiency of treatment (Bansal et al. 2009). As a result, in the context of cancer treatment, ligands that do not operate as p-glycoprotein substrates stand out as the most attractive possibilities.

The majority of phytochemicals, as shown in Table 4, have non-inhibitory effects on CYP3A4 and CYP1A2, which are key members of the cytochrome P450 family of drug-metabolizing enzymes. Cytochrome P450 enzymes are essential in the metabolism of several medications. Their interactions with pharmaceuticals can result in two distinct outcomes: rapid metabolism if the drugs function as substrates for certain CYP enzymes, resulting in drug induction; or drug accumulation if the drugs act as inhibitors of CYP enzymes, resulting in inhibition. Both of these possibilities are undesired (Ji et al. 2020). As a result, using in silico research to identify possible interactions between chemicals or medications and certain CYP isoenzymes is critical during the drug development process. The substances values for these parameters are impressively close to established standards, indicating their outstanding pharmacokinetic properties. This validates their ability to efficiently transfer metabolites to specific locations, indicating a distinct pharmacological significance.

Table 4 In silico pharmacokinetics of ligands using Swiss ADME

The Bioavailability Radar determines a compound's drug-likeness in real-time. The pink region in figure as depicted in Table 3 represents the best parameter range. A phytocompound that falls inside this pink region on the radar map is considered drug-like. This estimates the potential for oral bioavailability of the ligands. Flexibility (FLEX) and polarity (polar) are important factors influencing chemical bioavailability. Rotatable bonds impact flexibility, with substances with more than ten such bonds having limited oral bioavailability. Polarity, as assessed by topological polar surface area (TPSA), reveals that drugs with TPSA between 20Å2 and 130 Å2 have a high oral bioavailability (Khan et al. 2017).Therefore majority of phytochemicals meet radar plot requirements, showing their potential for oral bioavailability.

Homo lumo calculation

The frontier molecular orbital (FMO), which derives from the molecular orbitals with the greatest occupied (HOMO) and lowest unoccupied (LUMO) energies, plays an important role in showing compounds active sites during molecular interactions. This smart visualization assists in understanding the dynamic interaction of molecules. Figure 4 show the HOMO and LUMO areas for the compounds from the nine plants. LUMO refers to orbitals with the energy capacity to receive electrons, whereas HOMO refers to the area inside compounds with the greatest energy and most reactive electrons. As a result, it has been proposed that these bonding orbitals form strong interactions between ligands and receptors at these precise sites (Iwaloye et al. 2022).

Fig. 4
figure 4

The HOMO and LUMO diagrams of (1) Caffeic acid, (2) Ferulic Acid, (3) Hesperetic acid, (4) Stigamasterol, (5) Cis-p-coumaric acid, (6) Leuteolindin

Table 5 shows the HOMO/LUMO values predicted by DFT for the substances. The HOMO orbitals of the compounds have energy levels ranging from − 0.26751 to − 0.33045 eV, whereas the LUMO orbitals have energy levels ranging from 0.06266 to 0.18089 eV.

Table 5 Calculation of frontier molecular orbital

As a consequence, Caffeic acid was found to have the lowest LUMO energy level at − 0.06266, while Stigamasterol had the highest LUMO energy level at 0.18089. To acquire a better understanding, the Frontier Molecular Orbital (FMO) band gap (measured in eV) was calculated by subtracting the LUMO value from the HOMO value. This calculation is particularly important since it identifies both highly chemically reactive and moderately inert substances. It is generally understood that decreasing the energy gap increases a molecule’s chemical reactivity, whereas increasing the energy gap reduces reactivity (Olawale et al. 2022). This guiding concept was applied to the energy gap values calculated from the chemicals, resulting in Hesperetic acid being identified as the most chemically reactive molecule in the gaseous phase Table 5 Stigamasterol and Ferulic acid, on the other hand, appeared as the least reactive molecules among the tested candidates. The energy gap displayed by these plant-derived chemicals has the potential to reveal new insights into their interactions with biomolecules such as proteins. This study gives information on the processes driving these interactions and their possible metabolic consequences.

Discussion

Millet are staple diet source for developing nations since decades. Based on the nutritional profile, anti-inflammatory, anti-oxidant, hypocholesterolemic, hypoglycemic and anti-carcinogenic attributes of millet seed and derivative bioactive are considered miracle functional foods. Numerous research findings proved that millet seed derived phytochemicals slows down the metabolic disorder progression and some chronic diseases. In this article, Sorghum bicolor and Seteria italica phytochemicals anticancerous potential are discussed. The inhibitors of dihydrofolate reductase (DHFR), thymidylate synthase and serine hydroxyl methyltransferase of folate metabolic cycle are crucial targets for cancer treatment. Moreover, DHFR is involved in both DNA synthesis and amino acid metabolism, highlighting its importance as a therapeutic target in a variety of illnesses (Luo et al. 2006).

The folate metabolism inhibitors have received widespread attention for their effectiveness in treating cancer, rheumatoid arthritis, and microbial infections. This is primarily owing to their capacity to affect important enzymes such as thymidylate synthase, dihydrofolate reductase (DHFR), and serine hydroxymethyltransferase, all of which are linked to folate metabolism disorders. Notably, DHFR is involved in both DNA synthesis and amino acid metabolism, highlighting its importance as a therapeutic target in a variety of illnesses (Luo et al. 2006).

Rao et al. promising DHFR inhibitors were found utilizing in silico analysis. It was discovered that in human dihydrofolate reductase, Glu30 plays an important function in the active site by creating hydrogen bonds with inhibitors (Rao and Tapale, 2013). In our study we found similar hydrogen bond interaction in Hesperetic acid, Cis-p Coumaric acid obtained from Setaria italica which is also reported in (Zhang et al. 2023).The study found similar hydrogen bond interactions between Glu30 and both reference and newly discovered hit phytocompounds. It was discovered that the carboxylate oxygen of Glu30 establishes hydrogen bond interactions with the pyridopyrimidine ring found in inhibitor compounds (Rana et al. 2019).

The Sorghum bicolor and Setaria italica possess numerous bioactive phytochemical with diverse role from excellent nutrition to medicated properties. Keeping in view, the anticancerous activity of millet phytochemicals explored to cure cancer. To the best of our research investigation, this is first study reporting anticancerous activities via DHFR inhibition of hespertic acid, ferulic acid, stigmasterol, and caffeic acid.

Conclusion

In this study, molecular docking methods were used to identify six compounds isolated from Sorghum bicolor (Sorghum) and Setaria italica (Foxtail) for their inhibitory actions on Dihydrofolate reductase (DHFR) enzyme. Among these are Caffeic acid, Ferulic acid, Hesperetic acid, Stigmasterol, Cis-p-courmaric acid and Leuteolindin. The docking scores underwent a thorough post-docking examination to confirm the validity of these results. This study also used ADME predictions to show that these compounds had admirable drug-like qualities and pharmacokinetic characteristics. Hesperetic acid was found to be the chemical among the six that was found to be the most reactive when evaluated for reactivity in the gaseous phase. This finding emphasises the ability of bioactive substances to interact with growth factor receptors (GFRs) and block growth factors. In turn, this inhibition serves as a barrier to the growth of cancer cells against dihydrofolate reductase (DHFR) enzyme. This study implies that the lead compounds it has found have the potential for outstanding protein inhibitory activities, calling for more in vivo and in vitro research. These results also corroborate earlier in vivo studies on these phytochemicals, highlighting their promise as possible anti-cancerous options, as observed by other researchers in the area.

Data availability

Data will be available on request.

Code availability

Not applicable.

References

  • Bakari H, Djomdi Ruben ZF, Roger DD, Cedric D, Guillaume P, Pascal D, Philippe M, Gwendoline C. Sorghum (Sorghum bicolor L. Moench) and its main parts (By-Products) as promising sustainable sources of value-added ingredients. Waste Biomass Valorization. 2023;14(4):1023–44. https://doi.org/10.1007/s12649-022-01992-7.

  • Bansal T, Jaggi M, Khar R, Talegaonkar S. Emerging significance of flavonoids as P-glycoprotein inhibitors in cancer chemotherapy. J Pharm Pharm Sci. 2009;12(1):46–78.

    Article  CAS  PubMed  Google Scholar 

  • Belton PS, Taylor JRN. Sorghum and millets: protein sources for Africa. Trends Food Sci Technol. 2004;15(2):94–8.

    Article  CAS  Google Scholar 

  • Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–442.

    Article  PubMed  Google Scholar 

  • Chen Z, Li H, Zhang Q, Bao X, Yu K, Luo X, Zhu W. Pharmacophore-based virtual screening versus do- cking-based virtual screening : a benchmark compa- rison against eight targets. Acta Pharmacol Sin. 2009. https://doi.org/10.1038/aps.2009.159.

    Article  PubMed  PubMed Central  Google Scholar 

  • Cody V, Galitsky N, Luft JR, Pangborn W, Rosowsky A, Blakley RL. Comparison of two independent crystal structures of human dihydrofolate reductase ternary complexes reduced with nicotinamide adenine dinucleotide phosphate and the very tight-binding inhibitor PT523. Biochemistry. 1997;36(45):13897–903.

    Article  CAS  PubMed  Google Scholar 

  • de Morais Cardoso L, Pinheiro SS, Martino HSD, Pinheiro-Sant’Ana HM. Sorghum (Sorghum bicolor L.): Nutrients, bioactive compounds, and potential impact on human health. Crit Rev Food Sci Nutr. 2017;57(2):372–90. https://doi.org/10.1080/10408398.2014.887057.

  • Farsund AA, Daugbjerg C, Langhelle O. Food security and trade: reconciling discourses in the food and agriculture organization and the world trade organization. Food Secur. 2015;7:383–91.

    Article  Google Scholar 

  • Gangjee A, Vidwans AP, Vasudevan A, Queener SF, Kisliuk RL, Cody V, Li R, Galitsky N, Luft JR, Pangborn W. Structure-based design and synthesis of lipophilic 2, 4-diamino-6-substituted quinazolines and their evaluation as inhibitors of dihydrofolate reductases and potential antitumor agents. J Med Chem. 1998;41(18):3426–34.

    Article  Google Scholar 

  • Gupta M, Asfaha DM, Ponnaiah G. Millets: A Nutritional Powerhouse With Anti-cancer Potential. Cureus, 2023;15(10), e47769. https://doi.org/10.7759/cureus.47769.

  • Iwaloye O, Elekofehinti OO, Olawale F, Chukwuemeka PO, Kikiowo B, Folorunso IM. Fragment-based drug design, 2D-QSAR and DFT calculation: scaffolds of 1, 2, 4, triazolo [1, 5-a] pyrimidin-7-amines as Potential Inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase. Lett Drug Des Discov. 2022;20(3):317–34. https://doi.org/10.2174/1570180819666220422120707.

    Article  CAS  Google Scholar 

  • Ji D, Xu M, Udenigwe CC, Agyei D. Physicochemical characterisation, molecular docking, and drug-likeness evaluation of hypotensive peptides encrypted in flaxseed proteome. Curr Res Food Sci. 2020;3:41–50. https://doi.org/10.1016/j.crfs.2020.03.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Kalsi R, Bhasin JK. Nutritional exploration of foxtail millet (Setaria italica) in addressing food security and its utilization trends in food system. EFood. 2023;4(5), e111. https://doi.org/10.1002/efd2.111.

  • Kauther A, Elsheikh AM, Ahmed HM, Hassan OHA, Hamza NB, Osman MG, Daffalla HM. Phytochemical screening and antibacterial activities of sorghum bicolor leaves derived from in vitro culture. GSC Biol Pharm Sci. 2020;10(01):65–72.

    Article  Google Scholar 

  • Khan T, Dixit S, Ahmad R, Raza S, Azad I, Joshi S, Khan AR. Molecular docking, PASS analysis, bioactivity score prediction, synthesis, characterization and biological activity evaluation of a functionalized 2-butanone thiosemicarbazone ligand and its complexes. J Chem Biol. 2017;10(3):91–104. https://doi.org/10.1007/s12154-017-0167-y.

    Article  PubMed  PubMed Central  Google Scholar 

  • Lipinski CA, Lombardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings1PII of original article: S0169–409X(96), 00423–1. The article was originally published in advanced drug delivery reviews 23 (1997) 3. Adv Drug Deliv Rev. 2001;46(1):3–26. https://doi.org/10.1016/S0169-409X(00)00129-0.

    Article  CAS  PubMed  Google Scholar 

  • Luo W, Yu Q, Kulkarni SS, Parrish DA, Holloway HW, Tweedie D, Shafferman A, Lahiri DK, Brossi A, Greig NH. Inhibition of human acetyl-and butyrylcholinesterase by novel carbamates of (−)-and (+)-tetrahydrofurobenzofuran and methanobenzodioxepine. J Medic Chem. 2006;49(7):2174–85.

  • Maharajan T. Krishna TPA, Krishnakumar NM, Vetriventhan M, Kudapa H, Ceasar SA. Role of genome sequences of major and minor millets in strengthening food and nutritional security for future generations. Agriculture 2024;14(5). https://doi.org/10.3390/agriculture14050670.

  • Mehta D, Vyas S, Dudhagara D, Patel A, Parmar V. Significance of Indian millets in enhancing global food security: A comprehensive review. Trends Food Sci Technol. 2024;149, 104527. https://doi.org/10.1016/j.tifs.2024.104527.

  • McConkey BJ, Sobolev V, Edelman M. The performance of current methods in ligand–protein docking. Curr Sci. 2002;83:845–56.

    CAS  Google Scholar 

  • Mendie LE, Hemalatha S. Molecular docking of phytochemicals targeting GFRs as therapeutic sites for cancer: an in silico study. Appl Biochem Biotechnol. 2022;194(1):215–31. https://doi.org/10.1007/s12010-021-03791-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Nani A, Belarbi M, Ksouri-Megdiche W, Abdoul-Azize S, Benammar C, Ghiringhelli F, Hichami A, Khan NA. Effects of polyphenols and lipids from Pennisetum glaucum grains on T-cell activation: modulation of Ca2+ and ERK1/ERK2 signaling. BMC Complement Altern Med. 2015;15(1):426. https://doi.org/10.1186/s12906-015-0946-3.

  • Nivedita PS, Hoskeri JH, Minormillets K, Hoskeri JH. Minor Millet Phytochemicals and their Pharmacological Potentials. Pharmacogn Rev. 2022;16(32):100–6.

    Article  CAS  Google Scholar 

  • Okoye TC, Uzor PF, Onyeto CA, Okereke EK. 18 - Safe African Medicinal Plants for Clinical Studies. In: VBT-TS of AMP. Kuete, editors. 2014. p. 535–555. Elsevier. https://doi.org/10.1016/B978-0-12-800018-2.00018-2.

  • Olawale F, Iwaloye O, Olofinsan K, Ogunyemi OM, Gyebi GA, Ibrahim IM. Homology modelling, vHTS, pharmacophore, molecular docking and molecular dynamics studies for the identification of natural compound-derived inhibitor of MRP3 in acute leukaemia treatment. Chem Pap. 2022;76(6):3729–57.

    Article  CAS  Google Scholar 

  • Ornitz DM, Itoh N. The fibroblast growth factor signaling pathway. Wiley Interdiscip Rev Dev Biol. 2015;4(3):215–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Poshadri A, Deshpande HW, Kshirsagar RB. The international year of millets-2023, millets as nutricereals of 21st centenary for health and wellness. AATCC Review, 7. 2023.

  • Ramadoss DP, Sivalingam N. Vanillin extracted from Proso and Barnyard millets induce apoptotic cell death in HT-29 human colon cancer cell line. Nutr Cancer. 2020;72(8):1422–37. https://doi.org/10.1080/01635581.2019.1672763.

  • Rana RM, Rampogu S, Zeb A, Son M, Park C, Lee G, Yoon S, Baek A, Parameswaran S, Park SJ, Lee KW. In silico study probes potential inhibitors of human dihydrofolate reductase for cancer therapeutics. J Clinic Medic. 2019;8(2):1–21. https://doi.org/10.3390/jcm8020233.

  • Rao AS, Tapale SR. A study on dihydrofolate reductase and its inhibitors: A review. IJPSR, 2013;4(7): 2535.

  • Samtiya M, Aluko RE, Dhaka N, Dhewa T, Puniya AK. Nutritional and health-promoting attributes of millet: current and future perspectives. Nutr Rev. 2022;81(6):684–704. https://doi.org/10.1093/nutrit/nuac081.

  • Schnell JR, Dyson HJ, Wright PE. Structure, dynamics, and catalytic function of dihydrofolate reductase. Annu Rev Biophys Biomol Struct. 2004;33:119–40.

    Article  CAS  PubMed  Google Scholar 

  • Shan S, Li Z, Newton IP, Zhao C, Li Z, Guo M. A novel protein extracted from foxtail millet bran displays anti-carcinogenic effects in human colon cancer cells. Toxicol Lett. 2014;227(2), 129–38. https://doi.org/10.1016/j.toxlet.2014.03.008.

  • Sharma R, Sharma S, Dar BN, Singh B. Review Millets as potential nutri-cereals : a review of nutrient composition, phytochemical profile and techno-functionality. Int J Food Sci Technol. 2021. https://doi.org/10.1111/ijfs.15044.

    Article  Google Scholar 

  • Shih C-H, Siu S-O, Ng R, Wong E, Chiu LCM, Chu IK, Lo C. Quantitative analysis of anticancer 3-deoxyanthocyanidins in infected sorghum seedlings. J Agric Food Chem. 2007;55(2):254–259. https://doi.org/10.1021/jf062516t.

  • Shukla V, Srivastava S, Singh S, Mursal M, Hussain S. Unveiling the intricacies of phytate antinutrients in millets and their therapeutic implications in breast cancer. Intelligent Pharmacy. 2023. https://doi.org/10.1016/j.ipha.2023.12.005.

  • Siegel RL, Miller KD, Wagle NS, Jemal A. Cancer statistics, 2023. CA Cancer J Clin. 2023;73(1):17–48. https://doi.org/10.3322/caac.21763.

    Article  PubMed  Google Scholar 

  • Tiash S, Chowdhury EH. Growth factor receptors: promising drug targets in cancer. J Cancer Metastasis Treat. 2015;1:190–200.

    Article  CAS  Google Scholar 

  • Wishart DS. Improving early drug discovery through ADME modelling: an overview. Drugs R D. 2007;8:349–62.

    Article  CAS  PubMed  Google Scholar 

  • Zhang L, La X, Tian J, Li H, Li A, Liu Y, Wu C, Li Z. The phytochemical vitexin and syringic acid derived from foxtail fillet bran inhibit breast cancer cells proliferation via GRP78 / SREBP-1 / SCD1 signaling axis. J Funct Foods. 2021;85(June):104620. https://doi.org/10.1016/j.jff.2021.104620.

    Article  CAS  Google Scholar 

  • Zhang L, Liu Y, La X, Li S, Wen L, Liu T, Li H, Li A, Wu H, Wu C, Li Z. The bound polyphenols of foxtail millet (Setaria italica) inner shell inhibit breast cancer by promoting lipid accumulation-induced autophagic death. Food Chem Toxicol. 2023;177:113855. https://doi.org/10.1016/j.fct.2023.113855.

Download references

Funding

The authors declare that no funds, grants or other support were received during the research work and preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Kalpana Katiyar conceived and designed research. Anshika Gupta, Kalpana Katiyar and Akriti Verma performed all experimental analysis. All authors equally contributed in manuscript writing, editing and proof reading.

Corresponding author

Correspondence to Kalpana Katiyar.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors read and approved for publication.

Competing interests

The Authors declare no competing interest for this research work.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Verma, A., Gupta, A. & Katiyar, K. Computational investigation of the anticancer potential of Sorghum bicolor and Setaria italica phytochemicals against dihydrofolate reductase (DHFR) enzyme. CABI Agric Biosci 5, 66 (2024). https://doi.org/10.1186/s43170-024-00265-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s43170-024-00265-8

Keywords