Evidence-Based Complementary and Alternative Medicine
Volume 2018 (2018), Article ID 1943203, 12 pages
https://doi.org/10.1155/2018/1943203

Research Article

Exploring the Therapeutic Ability of Fenugreek against Type 2 Diabetes and Breast Cancer Employing Molecular Docking and Molecular Dynamics Simulations

Shailima Rampogu,¹ Saravanan Parameswaran,¹ Mary Rampogu Lemuel,² and Keun Woo Lee¹

¹Division of Life Science, Division of Applied Life Science (BK21 Plus), Plant Molecular Biology and Biotechnology Research Center (PMBBRC), Research Institute of Natural Science (RINS), Gyeongsang National University (GNU), 501 Jinju-daero, Jinju 52828, Republic of Korea
²West Thames College, London, UK

Correspondence should be addressed to Keun Woo Lee; kwlee@gnu.ac.kr

Received 9 February 2018; Revised 12 June 2018; Accepted 24 June 2018; Published 11 July 2018

Academic Editor: Almir Gonçalves Wanderley

Copyright © 2018 Shailima Rampogu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Fenugreek (Trigonella foenum-graecum) is used as a spice throughout the world. It is known for its medicinal properties such as antidiabetic, anticarcinogenic, and immunological activities. The present study shows the properties and the nutritional quality of fenugreek seed extract and focuses on screening of active compounds in drug designing for type 2 diabetes and breast cancer. Quantitative analysis was used to calculate the percentages of protein, carbohydrates moisture, fatty acid, galactomannan, oil, and amino acid. Phytochemical analysis revealed the presence of flavonoids, terpenoids, phenols, proteins, saponins, and tannins in fenugreek seed extracts. Molecular docking and molecular dynamics simulation-based computational drug discovery methods were employed to address the role of fenugreek seed constituents against type 2 diabetes and breast cancer. The computational results reveal that the compound galactomannan can be ascribed as potential drug candidate against breast cancer and type 2 diabetes rendered by higher molecular dock scores, stable molecular dynamics (MD) simulations results, and lower binding energy calculations.

1. Introduction

The legume fenugreek (Trigonella foenum-graecum) is a short annual plant from the Fabaceae family [1, 2]. The name Trigonella foenum-graecum is a Latin-Greek name as it bears a typical triangular shaped flowers and is employed as a common fodder for animals in Greece [1]. It is found in various parts of the globe and is often used as spice, condiment, and medication [3, 4]. Largely, fenugreek leaves and seeds have been used as spices in different parts of the globe. In Africa, fenugreek is used as supplement during bread preparation and the seed components of fenugreek are known to enhance the nutritional quality of the bread. In India, the leaves and seeds are utilized as favouring and seasoning agents [1]. In China, it is used as cure edema, while the ancient Egyptians employed fenugreek to incense the mummies [1, 5, 6]. Additionally, fenugreek is used as a medicine to treat several diseases besides being used as antioxidant [7], against inflammation [8, 9], as anticancer [10], as hepatoprotective agent [11, 12], as antibacterial [13–15], and as antifungal [16]. Additionally, fenugreek is also used as off-season fodder and animal food supplement [17].

Fenugreek seeds are widely studied part of the plant. The powdered fenugreek is used as condiment and the seed endosperm serves to secure fenugreek gum [1]. The seeds have a strong aroma with bitter taste [18]. The major chemical constituents found in fenugreek seed are galactomannan (fibre), diosgenin (sapogenin), trigonelline (alkaloid), and 4-hydroxyisoleucine that have the antidiabetic properties and are also employed to treat breast cancer [19].

Diabetes mellitus is a common and chronic disease concern globally associated with a ten-year shorter life expectancy [20]. According to WHO, type 2 diabetes occurs because either body does not produce enough insulin or body resists the effects of insulin [21, 22]. Type 2 diabetes is dominant in developing countries and accounts to around 85%–90% worldwide [20, 21]. Fenugreek is another promising antidiabetic drug [23]. It was also confirmed that consuming fenugreek as a dietary supplement in the prediabetic patients could efficiently reduce the outbreak of type 2 diabetes [24]. Additionally, it was further reported that the socked fenugreek seeds can act as adjuvant in mitigating the type 2 diabetes and also in noninsulin dependent diabetes [25, 26] and serum lipids in type I diabetes [27]. Additionally, it is well evidenced that the fenugreek seeds are antidiabetic in nature [24, 28, 29].

Fenugreek also possesses anticancer properties and chemical constituents of fenugreek are known to induce apoptosis [30, 31]. Furthermore, it induces dose-dependent effect on human breast cancer cell line [32]. Breast cancer is the most common cause of death in female worldwide [33, 34]. The discovery of BRCA1 and BRCA2 genes helped to understand that hereditary factor is the main cause of most cancers [35]. Chloroform seed extract studies by Khoja et al. proved the effective killing of MCF-7 human immortalized breast cells [30]. Amin et al. (2005) studies suggest that fenugreek seed chemical constituents have preventive effect against breast cancer which inhibit MDA 231-induced mammary hyperplasia [36]. However, it is not yet delineated on the most effective compound that can act on both the morbidities. Therefore, in the current investigation, we employed the computational technique such as molecular docking and molecular dynamics simulations to identify candidate compound as compared with the reference compounds.

Molecular docking is one of the widely adapted methods to predict the binding affinities between the ligand and the target protein and further the lead optimization [37]. Additionally, the molecular docking imparts knowledge on the interactions at the atomic level [37] and predicts the ideal binding mode [33, 38]. Molecular docking mechanism generally evaluates the binding conformations, its orientation, and the accommodation of the small molecule at the active site of the proteins binding site and are read as scores [39]. The molecular dynamics simulation imparts knowledge on the nature of the small molecules at the proteins binding pocket thereby affirming the appropriate binding modes [38]. The identified Hits that have demonstrated a higher dock score than the reference compounds or the known drugs, exhibiting the interactions with the key residues complemented by stable molecular dynamics simulation results, are considered the most promising candidate compounds.

In the current investigation, the quantitative analysis of fenugreek seeds was conducted to gain information on the components and further the computational analysis was performed to discover the potential compound against breast cancer and type 2 diabetes. The in silico results have illuminated galactomannan as the prospective compound against both diseases.

2. Materials and Methods

Fenugreek seeds were used as a sample to test the medicinal properties. Fenugreek seeds were sourced from a local market (Hyderabad, India) and were of high quality grade. They were shade dried, cleaned, and finely powdered and used for chemical analysis.

2.1. Biochemical Analysis

The biochemical studies were carried out to identify the protein content, total soluble carbohydrates, oil content and fatty acid values, free amino acids, and soluble fibres from the collected seed samples.

2.1.1. Estimation of Total Protein

Percentage of proteinaceous nitrogen and proteins was estimated by the Micro-kjeldahl method [40]. Proteinaceous nitrogen was calculated by the following formula.

 T is titration reading of the sample, B is blank reading of the sample, S is the amount of sample taken in grams, N is normality of hydrochloric acid (N/28).

To calculate the percentage of protein, the nitrogen value was multiplied by the factor 6.25.

2.1.2. Estimation of Total Carbohydrate

Total carbohydrate content of the seed samples was estimated by the procedure suggested by Loewis (1952) [41]. Anthrone reagent was used and the developed colour was read at 620nm in a colorimeter against blank.

2.1.3. Estimation of Oil Content

Total oil content of the said spices was estimated as suggested by Meara (1955) [42].

Percentage of oil was calculated by following formula:

 Wo is the weight of oil extracted, Ws is the weight of seed taken.

2.1.4. Estimation of Fatty Acid Value

Method used to estimate the fatty acid value is suggested by Meara (1955) [42].

Fatty acid value was calculated using the formula

 U is the volume of titration of 0.1 n KOH, W is the grams of oil taken.

2.1.5. Isolation of Amino Acids

Column chromatography was used to isolate free amino acids from fenugreek seeds [43].

To find the concentration of 4-hydroxyisoleucine, first the total amino acid content was determined by using spectrophotometric method. Then the relative concentration of 4-hydroxyisoleucine in the mixture of amino acid was determined by high performance thin layer chromatography (HPTLC).

2.1.6. Isolation of Galactomannans

Extraction and isolation of the water-soluble polysaccharides (galactomannans) from endosperm of fenugreek seeds were done using the procedure of Kooiman (1971) [44].

2.1.7. Estimation of Moisture Percentage

Moisture content of seeds was estimated by “Dry air oven” method association of official analytical chemists (AOAC) (anonymous, 1947)[45] and the percentage was calculated from the following formula:

2.2. Molecular Docking, Simulations, and Free Energy Calculations

To further assess the suitability of the compounds as antidiabetic and potential breast cancer agents, the investigation proceeds employing the computational methods such as molecular docking recruiting CDOCKER available on Discovery Studio (DS) v4.5, molecular dynamics (MD) simulations conducted using GROningen MAchine for Chemical Simulations (Gromacs) v5.0, which was followed by MM/PBSA calculations.

2.2.1. Molecular Docking

For the execution of the docking protocol, the proteins for both the diseases were imported from protein data bank (PDB) of high resolution. The protein with the PDB id 3EQM (2.9Å) was chosen for breast cancer and 1GFY (2.1Å) was elected for type 2 diabetes, respectively. These proteins were prepared on DS by initiating the clean protein module embedded with the DS and subsequently heteroatoms together with the water molecules were dislodged and the addition of hydrogens was performed adapting the CHARMm force field accessible on the DS. The active sites were selected in accordance with the co-crystal geometry, thereby, considering the residues around 10 Å radius [46, 47].

Phytochemicals along with the type 2 diabetic and breast cancer drugs, canagliflozin [48] and anastrozole [49], were used to comparatively evaluate the effect of the prospective drug molecules on the diseases labelling the latter as reference drug. These compounds were imported onto the DS to obtain their 3D structures and were subsequently minimized. The prepared proteins and the ligands were subjected to molecular docking studies employing the CDOCKER protocol.

CDOCKER available on the DS happens to be the most reliable method as it employs the CHARMm-based dynamics methods [50]. Subsequently, 30 conformations were allowed to be generated for each ligand, while the other parameters were set at default. The results were evaluated based upon the higher –CDOCKER interaction energy and higher –CDOCKER energy that significantly correspond to the favourable binding. The most appropriate binding mode was judged by the maximum clusters formed and was therefore subjected to MD simulations to understand its dynamic behaviour.

2.2.2. MD Simulations

Molecular dynamics (MD) simulations were performed for the favourable systems obtained after docking using GROMACS 5.0 with CHARMm27 force field. Ligand topologies were generated adapting the SwissParam [51]. All the parameters were attributed as described earlier [52–56]. Dodecahedron water box was generated and the systems were solvated comprising three-site transferrable intermolecular potential (TIP3P) water model, to which the counter ions were added. The system was energy minimized with steepest descent algorithm with 10000 steps which was then subjected to equilibration using constant number N, volume V, and temperature T (NVT) [57] and constant number N, pressure P, and temperature T (NPT) [58]. During this process, the protein backbone was restrained and the periodic boundary conditions were fostered to avoid bad effects. Thereafter, the MD run was conducted for 10 ns, saving the data for every one picosecond (ps). Visual molecular dynamics (VMD)[59] and DS were utilized to analyse the MD results.

2.2.3. Binding Free Energy Calculations

Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) was recruited to compute the binding free energy calculations [60, 61]. 10 snapshots were evenly extracted from the MD trajectories of the protein ligand complex. A variety of energetic values were calculated using

3. Results

3.1. Biochemical Analysis

The total seed percentage revealed that galactomannan and 4-hydroxyisoleucine were present in 26.4 and 13 percentages, respectively, as in Table 1.

Table 1: Percentage of the seed contents.

Contents of fenugreek seed extract
Average percentage of the seed extracts (%)

protein
28.5

carbohydrate
16.2

oils
5.3

fatty acid
3.8

galactomannan
26.4

moisture
6.8

4-hydroxyisoleucine
13

Further phytochemical screening of acetone seed extract of fenugreek was carried out to test the presence of tannins, phenols, terpenoids, flavonoids, saponins, and alkaloids [62] and are tabulated in Table 2.

Table 2: Summary of Phytochemicals in Acetone Extract of Fenugreek Seed.

Tests
Results

Flavanoid
+ve

Tannin
+ve

Terpenoids
+ve

Alkaloids
+ve

Saponins
+ve

Test for flavonoids: 1 ml of extract in a test tube and 5ml of diluted ammonium solution were added followed by few drops of concentrated sulphuric acid. Formation of yellow colour indicated the presence of flavonoids [62].

Test for tannins: Formation of reddish-brown colour indicated the presence of tannins (ferric chloride test) when 1% ferric chloride solution was added to 1 ml of extract of fenugreek seeds [62].

Test for terpenoids: To find out the presence of terpenoids, Salkowski test was conducted. 1 ml of extract was taken and dissolved in chloroform and then a few drops of concentrated sulphuric acid were added to it. On the inner face, a reddish-brown colour was formed that indicated the presence of terpenoids [62].

Test for alkaloids: Dragendorffa’s test results indicated the presence of alkaloids by giving orange-red precipitate, when 1 ml of Dragendroffa’s reagent was added (potassium bismuth iodide solution) to 1 ml of extract [62].

Test for saponins: Frothing test was conducted to test for saponins in the seed extract. 1ml of extract was vigorously shaken with distilled water and was allowed to stand for 10 min. Stable froth indicated the presence of saponins [62].

3.2. Molecular Docking, Simulations, and Free Energy Calculations

3.2.1. Molecular Docking Studies

Molecular docking was executed independently for diabetes and breast cancer. The ligands along with their respective proteins were docked to assess their binding affinities. It was interesting to note that 4-hydroxyisoleucine has generated a relatively lower dock score while galactomannan produced higher dock score as compared to their respective reference compounds, as in Table 3. Therefore, 4-hydroxyisoleucine was refrained from further calculations and the other systems were proceeded forward.

Table 3: Molecular dock scores between the drug targets and the compounds.

S. no.
Name of the compound
-CDOCKER interaction energy

Dock scores of diabetes mellitus

1
canagliflozin
36.55

2
galactomannan
43.19

3
4-hydroxyisoleucine
28.27

Dock scores of breast cancer

1
anastrozole
34.05

2
galactomannan
58.15

3
4-hydroxyisoleucine
23.88

3.2.2. Molecular Dynamics Simulations

To secure the results obtained from the docking, the MD simulations were performed to establish the most reliable ligand-receptor complex and additionally to understand their behaviour at proteins active site [52, 53]. The MD for 10 ns was initiated and the behaviour of each system was monitored. Accordingly, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and potential energies were calculated for each system. The RMSD for the breast cancer systems were observed to be stable after 4000 ps with no significant variation, thereafter, implying that the system is well converged, as in Figure 1. Moreover, the RMSD values were demonstrated to be less than 0.25 nm. Similar results were noted with RMSF values as well, as in Figure 2. The potential energy further states that there were no abnormal behaviours of the systems which were stable throughout the simulations, as in Figure 3. The last 5ns trajectories were retrieved to study the binding mode analysis. Upon superimposition, it was conceived that the binding mode pattern of the reference and the galactomannan were similar, as in Figure 4. The interactions of the ligand with the protein were evaluated with the key residues located at the active site. The reference compound anastrozole was seen to form a hydrogen bond with the NH atom of Met374 residue, joined by N5 atom with a bond length of 2.9 Å. Phe134 was found to form the π – π with the ligand molecule. Galactomannan was found to interact with the protein by forming 7 hydrogen bonds. The O13 atom of the ligand has interacted with the HH22 atom of Arg115 with a bond length of 2.8 Å. The HH21 atom of Arg115 has interacted with O15 atom of the ligand with a bond length of 2.5 Å. The O atom of Ile132 has joined with H62 of the ligand displaying a bond length of 2.6 Å. Another hydrogen bond was observed between the HH11 atom of Arg145 and the O14 atom of the ligand with a length of 2.0 Å. The OD2 atom of the residue Asp309 has interacted with the H57 of the ligand with a bond distance of 2.8 Å. The O atom of the key residue Met374 has interacted with the H53 atom of the ligand with a bond length of 2.5 Å. The SG atom of the Cys437 residue has interacted with the H63 atom of the ligand with a distance of 2.5 Å. The details of the interaction are represented in Figure 5 and Table 4. Furthermore, the intermolecular hydrogen bond interactions were recorded during the simulations to elucidate deposition of the ligand within the active site. It was observed that the reference molecule displayed an average of 0.3 hydrogen bonds, while those within 0.35 nm were observed to be 0.7, as in Figure 6, while the candidate molecule demonstrated an average of 1.3 hydrogen bonds and the bonds within 0.35 nm were 4.4, as in Figure 7.

Table 4: The molecular interactions between the compounds and the protein.

S. no.
Compound
Ligand Atom
Amino acid
Amino acid atom
Bond length (Å)
Hydrophobic interactions

1
anastrozole
N5
Met374
HN
2.9
Ile133,Asp309,Val370, Leu372, Val373,Pro429,Phe430,Cys437,Leu 477

2
galactomannan
O13
Arg115
HH22
2.8
Ala306, Asp309, Phe430.

O15
Arg115
HH21
2.5

H62
Ile132
O
2.6

O14
Arg145
HH11
2.0

H57
Asp309
OD2
2.8

H53
Met374
O
2.5

H63
Cys437
SG
2.5

Figure 1: RMSD plots for backbone atoms.

Figure 2: RMSF profiles for backbone atoms.

Figure 3: Potential energy graphs of the systems.

Figure 4: Binding mode assessment of the reference (cyan) and galactomannan (pink). Superimposition of the representative structures (left) and zoomed (right).

Figure 5: Depiction of hydrogen bond interactions and binding conformations. Only polar atoms are displayed for clarity.

Figure 6: Graphical depiction of number of intermolecular hydrogen bond interactions between the protein and the reference compound.

Figure 7: Graphical depiction of number of intermolecular hydrogen bond interactions between the protein and the candidate compound.

Similar types of calculations were determined for the type 2 diabetes disease target and its respective ligands. The RMSD was recorded to be stable after 7000 ps for both the reference and galactomannan. Further, it was noted that the RMSD of the reference was established to be within 0.2 nm while the drug-like molecule demonstrated a RMSD within 0.15 nm, as in Figure 8. However, no major fluctuations were noticed during the simulations referring to the stability of the systems. The same results were depicted through the RMSF, as in Figure 9, and the potential energy calculations, as in Figure 10. Therefore, to examine the binding mode of the ligand molecules, the last 5 ns trajectories were extracted and were superimposed. The results represented a similar binding mode between the reference and the galactomannan, as in Figure 11. Furthermore, intermolecular interactions were inspected with the key residues residing at the active site. It revealed that the reference molecule has formed three hydrogen bonds with the active site residues. The F2 of the ligand has interacted with the HG atom of Cys215 with bond length of 2.6 Å. The other two hydrogen bonds are formed with HN and HE atoms of Arg221 and 2.1 Å each. Tyr46 and Phe182 have been involved with the π – π interactions. On the contrary, Galactomannan on the other hand generated eight hydrogen bonds, two hydrogen bonds with Lys120 and Asp181 amino acid residues and one hydrogen bond with Arg221, Ser216, Gln262, and Gln266, respectively. The details of the interactions are tabulated in Figure 12 and Table 5. Furthermore, the intermolecular hydrogen bonds were evaluated throughout the simulations. The average hydrogen bonds were computed to be 0.09 and those within 0.35 nm were found to be 0.7, as in Figure 13. The prospective drug molecule however has represented average hydrogen bonds of 3.9 while the bonds within 0.35 nm were enumerated to be 4.4 projecting the superiority of galactomannan, as in Figure 14.

Table 5: The molecular interactions between the compounds and the protein.

S.no
Compound
Ligand Atom
Amino acid
Amino acid atom
Bond length Å
Hydrophobic interactions

1
canagliflozin
F2
Cys215
HG
2.6
Lys120,Lys116, Ser216,Gly218, Ile219,Gly220, Ala217,Gln262

F2
Arg221
HN
2.1

F2
Arg221
HE
2.1

2
galactomannan
O9
Lys120
HZ2
1.7
Tyr46,Lys116, Phe182,Gly183, Cys215,Ser216, Gly218,Ile219, Gly220

O3
Lys120
HZ1
2.0

H66
Asp181
OD1
2.3

H64
Asp181
OD1
1.9

O14
Arg221
HN
2.4

O16
Ser216
HN
2.4

H62
Gln262
OE1
2.1

O13
Gln266
HE22
2.4

Figure 8: RMSD plots for backbone atoms.

Figure 9: RMSF profiles for backbone atoms.

Figure 10: Potential energy graphs of the systems.

Figure 11: Binding mode assessment of the reference (purple) and galactomannan (orange). Superimposition of the representative structures (left) and zoomed (right).

Figure 12: Depiction of hydrogen bond interactions and binding conformations. Only polar atoms are displayed for clarity.

Figure 13: Graphical depiction of number of intermolecular hydrogen bond interactions between the protein and the reference compound.

Figure 14: Graphical depiction of number of intermolecular hydrogen bond interactions between the protein and the candidate compound.

3.3. Binding Free Energy Analysis

Binding free energies are computed after the MD simulations that inspect protein fluctuations and ligand conformations thereby ensuring a suitable positioning of the ligand within the binding site. The MM/PBSA calculations have produced a favourable ΔG that ranged between -10 to 100 kJ/mol for breast cancer target, as in Figure 15. Furthermore, the average binding energy produced by reference was -42.45 kJ/mol while that generated by galactomannan was -47.95 kJ/mol, respectively, as in Table 6.

Table 6: Comparative assessment between dock scores and the binding energies of breast cancer systems.

S. no.
Name of the compound
-CDOCKER interaction energy
Average binding energy (kJ/mol)

1
anastrozole
34.05
-42.45

2
galactomannan
58.15
-47.95

Figure 15: MM/PBSA binding energy representation of the reference and the candidate compound.

The binding free energies were subsequently calculated for canagliflozin-protein and galactomannan-protein systems for type 2 diabetes. 10 snapshots were evenly extracted and the binding energies were computed accordingly. The binding energies ranged between -15 kJ/mol and -100 kJ/mol, as in Figure 16. Additionally, it was observed that the average binding energy was calculated as -51.75 kJ/mol for the reference and -68.11 kJ/mol for galactomannan, as in Table 7.

Table 7: Comparative assessment between dock scores and the binding energies type 2 diabetes systems.

S. no.
Name of the compound
-CDOCKER interaction energy
Average binding energy (kJ/mol)

1
canagliflozin
36.55
-51.75

2
galactomannan
43.19
-68.11

Figure 16: MM/PBSA binding energy representation of the reference and the candidate compound.

From the results, it is evident that galactomannan has represented higher –CDOCKER interaction energy values and lower binding free energies than their respective reference compounds. These results demonstrate that galactomannan has stronger binding affinities than the reference inhibitors.

4. Discussion

In the present study, the seed extract showed the presence of proteins, carbohydrates, fatty acids, oils, saponin, flavonoids, tannins, terpenoids, alkaloids, soluble fibre galactomannan, and amino acid 4 hydroxyisoleusine (Tables 1 and 2). Some chemicals screened are similar to the work done by Yadav R. et al. 2014 [63].

Out of these chemicals, the special interest in this investigation is on the percentages of soluble fibre galactomannan 26.4 % and free amino acids 4 hydroxyleucine 13% and the presence saponins, as these are linked to human health benefits mainly in the reduction of plasma glucose levels and anticancer activities [64].

In order to further evaluate molecular inhibitory effect of the selected phytochemicals, the investigation proceeds in silico. Computational results have revealed that the phytochemical 4 hydroxyisoleucine could not induce the inhibitory activity against both the diseases. Although reports exist to explain its antidiabetic and antibreast cancer activity, the present finding foretells its inability as an inhibitor [31, 65]. Therefore, this amino acid was not forwarded for further studies. The other compound galactomannan has proved to be potential against both the diseases. This was represented by the RMSD, RMSF, and the potential energy values. The results were found to be unaltered as compared with the reference throughout the simulations. Moreover, the binding energies of the prospective drug molecules are found to be less, while rendering the highest –CDOCKER interaction energies. It is documented from the previous reports regarding the role of breast cancer inhibitors on diabetes mellitus as there exists a linkage between them [66, 67]. All the above results conclude that galactomannan could be considered as a potential drug for both the diseases.

Chemically, galactomannan is a polysaccharide molecule comprising a mannose backbone and the galactose side groups, hence, the name. More precisely they exist with 1-6 alpha-D-galactopyranose linkage. However, in fenugreek, mannose and galactose are linked by 1:1 linkage. Upon observing the docking conformations, it can be elucidated that the galactose side groups have involved in forming the hydrogen bond interaction with the active side residues, with the ring structures of the mannose involved in the formation of the π bond interactions.

In conclusion, the present study has examined the active components of fenugreek seeds against two common but different diseases, viz-a-viz: type-2 diabetes and breast cancer, using a well-established computational drug discovery method. The chemical composition of fenugreek seeds was assessed, and galactomannan and 4-hydroxyisoleucine were identified as major components and are similar to previous studies [68]. The therapeutic potential of these two identified active components was further assessed using molecular docking and molecular dynamics simulations. Our results identify galactomannan as a potential active component of fenugreek seeds, with a docking score compared to established drugs such as canagliflozin and anastrozole in binding simulations of therapeutics against type-2 diabetes and breast cancer, respectively. These results establish galactomannan, derived from fenugreek seeds, as a potential candidate for further drug discovery experiments in establishing their value as therapeutics against type-2 diabetes and breast cancer.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Shailima Rampogu and Saravanan Parameswaran contributed equally to this work.

Acknowledgments

This research was supported by the Pioneer Research Center Program (NRF-2015M3C1A3023028) through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning.

References

K. C. Nagulapalli Venkata, A. Swaroop, D. Bagchi, and A. Bishayee, “A small plant with big benefits: Fenugreek (Trigonella foenum-graecum Linn.) for disease prevention and health promotion,” Molecular Nutrition & Food Research, vol. 61, no. 6, Article ID 1600950, 2017.

C. Poole, B. Bushey, C. Foster et al., “The effects of a commercially available botanical supplement on strength, body composition, power output, and hormonal profiles in resistance-trained males,” Journal of the International Society of Sports Nutrition, vol. 7, 2010.

S. Shabbeer, M. Sobolewski, R. K. Anchoori et al., “Fenugreek: a naturally occurring edible spice as an anticancer agent,” Cancer Biology & Therapy, vol. 8, no. 3, pp. 272–278, 2009.

S. A. Wani and P. Kumar, “Fenugreek: A review on its nutraceutical properties and utilization in various food products,” Journal of the Saudi Society of Agricultural Sciences, vol. 17, no. 2, pp. 97–106, 2018.

E. Basch, C. Ulbricht, G. Kuo, P. Szapary, and M. Smith, “Therapeutic applications of fenugreek,” Alternative Medicine Review, vol. 8, no. 1, pp. 20–27, 2003.

D. Tiran, “The use of fenugreek for breast feeding women,” Complementary Therapies in Nursing and Midwifery, vol. 9, no. 3, pp. 155-156, 2003.

K. Szabó, R. Gesztelyi, N. Lampé et al., “Fenugreek (Trigonella Foenum-Graecum) Seed Flour and Diosgenin Preserve Endothelium-Dependent Arterial Relaxation in a Rat Model of Early-Stage Metabolic Syndrome,” International Journal of Molecular Sciences, vol. 19, no. 3, p. 798, 2018.

N. Sharma, S. Suresh, A. Debnath, and S. Jha, “Trigonella seed extract ameliorates inflammation via regulation of the inflammasome adaptor protein, ASC,” Frontiers in Bioscience - Elite, vol. 9, no. 2, pp. 246–257, 2017.

K. Pundarikakshudu, D. H. Shah, A. H. Panchal, and G. C. Bhavsar, “Anti-inflammatory activity of fenugreek (Trigonella foenum-graecum Linn) seed petroleum ether extract,” Indian Journal of Pharmacology, vol. 48, no. 4, pp. 441–444, 2016.

G. Sethi, M. Shanmugam, S. Warrier et al., “Pro-Apoptotic and Anti-Cancer Properties of Diosgenin: A Comprehensive and Critical Review,” Nutrients, vol. 10, no. 5, p. 645, 2018.

A. R. Shivashankara, A. Azmidah, R. Haniadka, M. P. Rai, R. Arora, and M. S. Baliga, “Dietary agents in the prevention of alcohol-induced hepatotoxicty: Preclinical observations,” Food & Function, vol. 3, no. 2, pp. 101–109, 2012.

S. Kaviarasan and C. V. Anuradha, “Fenugreek (Trigonella foenum graecum) seed polyphenols protect liver from alcohol toxicity: a role on hepatic detoxification system and apoptosis,” Die Pharmazie, vol. 62, no. 4, pp. 299–304, 2007.

D. Bano, H. Tabassum, A. Ahmad, A. Mabood, and I. Z. Ahmad, “The medicinal significance of the bioactive compounds of trigonella foenum-graecum: a review,” International Journal of Research in Ayurveda & Pharmacy, vol. 7, no. 4, pp. 84–91, 2016.

S. Goyal, N. Gupta, and S. Chatterjee, “Investigating therapeutic potential of trigonella foenum-graecum L. As our defense mechanism against several human diseases,” Journal of Toxicology, vol. 2016, 2016.

R. Premanath, J. Sudisha, N. L. Devi, and S. M. Aradhya, “Antibacterial and anti-oxidant activities of fenugreek (Trigonella foenum graecum L.) leaves,” Research Journal of Medicinal Plant, vol. 5, no. 6, pp. 695–705, 2011.

R. Haouala, S. Hawala, A. El-Ayeb, R. Khanfir, and N. Boughanmi, “Aqueous and organic extracts of Trigonella foenum-graecum L. inhibit the mycelia growth of fungi,” Journal of Environmental Sciences, vol. 20, no. 12, pp. 1453–1457, 2008.

A. Ahmad, S. S. Alghamdi, K. Mahmood, and M. Afzal, “Fenugreek a multipurpose crop: Potentialities and improvements,” Saudi Journal of Biological Sciences, vol. 23, no. 2, pp. 300–310, 2016.

E. Altuntaş, E. Özgöz, and Ö. F. Taşer, “Some physical properties of fenugreek (Trigonella foenum-graceum L.) seeds,” Journal of Food Engineering, vol. 71, no. 1, pp. 37–43, 2005.

S. Rizvi and N. Mishra, “Traditional Indian Medicines Used for the Management of Diabetes Mellitus,” Journal of Diabetes Research, vol. 2013, Article ID 712092, pp. 1–11, 2013.

M. S. Kirkman, V. J. Briscoe, N. Clark et al., “Diabetes in older adults,” Diabetes Care, vol. 35, no. 12, pp. 2650–2664, 2012.

H. Schneider, J. Shaw, and P. Zimmet, “Guidelines for the Detection of Diabetes Mellitus - Diagnostic Criteria and Rationale for Screening,” The Clinical Biochemist Reviews, vol. 24, no. 3, pp. 77–80, 2003.

L. Bellamy, J. P. Casas, A. D. Hingorani, and D. Williams, “Type 2 diabetes mellitus after gestational diabetes: a systematic review and meta-analysis,” The Lancet, vol. 373, no. 9677, pp. 1773–1779, 2009.

F. H. Moghadam, B. Vakili-Zarch, M. Shafiee, and A. Mirjalili, “Fenugreek seed extract treats peripheral neuropathy in pyridoxine induced neuropathic mice,” EXCLI Journal, vol. 12, pp. 282–290, 2013.

A. Gaddam, C. Galla, S. Thummisetti, R. K. Marikanty, U. D. Palanisamy, and P. V. Rao, “Role of Fenugreek in the prevention of type 2 diabetes mellitus in prediabetes,” Journal of Diabetes and Metabolic Disorders, vol. 14, no. 1, 2015.

Z. Madar, R. Abel, S. Samish, and J. Arad, “Glucose-lowering effect of fenugreek in non-insulin dependent diabetics,” European Journal of Clinical Nutrition, vol. 42, no. 1, pp. 51–54, 1988.

M. Attokaran, Effectiveness of phytotherapy in supportive treatment of type 2 diabetes mellitus II. Fenugreek (Trigonella foenum-graecum), Česká a slovenská farmacie, vol. 64, pp. 67–71, 2015.

R. D. Sharma, T. C. Raghuram, and N. S. Rao, “Effect of fenugreek seeds on blood glucose and serum lipids in type I diabetes,” European Journal of Clinical Nutrition, vol. 44, 1990.

G. S. Kumar, A. K. Shetty, K. Sambaiah, and P. V. Salimath, “Antidiabetic property of fenugreek seed mucilage and spent turmeric in streptozotocin-induced diabetic rats,” Nutrition Research, vol. 25, no. 11, pp. 1021–1028, 2005.

N. Neelakantan, M. Narayanan, R. J. De Souza, and R. M. Van Dam, “Effect of fenugreek (Trigonella foenum-graecum L.) intake on glycemia: A meta-analysis of clinical trials,” Nutrition Journal , vol. 13, no. 1, article no. 7, 2014.

K. K. Khoja, G. Shafi, T. N. Hasan et al., “Fenugreek, a naturally occurring edible spice, kills MCF-7 human breast cancer cells via an apoptotic pathway,” Asian Pacific Journal of Cancer Prevention, vol. 12, no. 12, pp. 3299–3304, 2011.

M. I. M. Khalil, M. M. Ibrahim, G. A. El-Gaaly, and A. S. Sultan, “Trigonella foenum (Fenugreek) Induced Apoptosis in Hepatocellular Carcinoma Cell Line, HepG2, Mediated by Upregulation of p53 and Proliferating Cell Nuclear Antigen,” BioMed Research International, vol. 2015, Article ID 914645, pp. 1–11, 2015.

S. Vígh, Z. Zsvér-Vadas, C. Pribac et al., “Fenugreek (Trigonella foenum-graecum l.) extracts are inducing dose-dependent hormetic response and cytotoxic effects in case of human breast cancer cell lines,” Studia Universitatis Vasile Goldis Arad, Seria Stiintele Vietii, vol. 26, no. 4, pp. 435–448, 2016.

S. Rampogu, M. Son, A. Baek et al., “Targeting natural compounds against {HER}2 kinase domain as potential anticancer drugs applying pharmacophore based molecular modelling approaches,” Computational Biology and Chemistry, vol. 74, pp. 327–338, 2018.

J. Ferlay, H. R. Shin, F. Bray, D. Forman, C. Mathers, and D. M. Parkin, “Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008,” International Journal of Cancer, vol. 127, no. 12, pp. 2893–2917, 2010.

N. Petrucelli, M. B. Daly, and G. L. Feldman, “Hereditary breast and ovarian cancer due to mutations in BRCA1 and BRCA2,” Genetics in Medicine, vol. 12, no. 5, pp. 245–259, 2010.

A. Amin, A. Alkaabi, S. Al-Falasi, and S. A. Daoud, “Chemopreventive activities of Trigonella foenum graecum (Fenugreek) against breast cancer,” Cell Biology International, vol. 29, no. 8, pp. 687–694, 2005.

G. Wang and W. Zhu, “Molecular docking for drug discovery and development: a widely used approach but far from perfect,” Future Medicinal Chemistry, vol. 8, no. 14, pp. 1707–1710, 2016.

S. Rampogu, A. Baek, A. Zeb, and K. W. Lee, “Exploration for novel inhibitors showing back-to-front approach against VEGFR-2 kinase domain (4AG8) employing molecular docking mechanism and molecular dynamics simulations,” BMC Cancer, vol. 18, no. 1, 2018.

X.-Y. Meng, H.-X. Zhang, M. Mezei, and M. Cui, “Molecular docking: a powerful approach for structure-based drug discovery,” Current Computer-Aided Drug Design, vol. 7, no. 2, pp. 146–157, 2011.

H. Wang, N. Pampati, W. M. McCormick, and L. Bhattacharyya, “Protein nitrogen determination by kjeldahl digestion and ion chromatography,” Journal of Pharmaceutical Sciences, vol. 105, no. 6, pp. 1851–1857, 2016.

F. A. Loewus, “Improvement in Anthrone Method for Determination of Carbohydrates,” Analytical Chemistry, vol. 24, no. 1, p. 219, 1952.

K. Paech and M. V. Tracey, Modern Methods of Plant Analysis / Moderne Methoden der Pflanzenanalyse, Springer Berlin Heidelberg, Berlin, Heidelberg, 1955.

V. Rolland-Fulcrand, M. Rolland, M.-L. Roumestant, and J. Martinez, “Chemoenzymatic synthesis of enantiomerically pure (2S,3R,4S)-4- hydroxyisoleucine, an insulinotropic amino acid isolated from fenugreek seeds,” European Journal of Organic Chemistry, no. 4, pp. 873–877, 2004.

P. Kooiman, “Structures of the galactomannans from seeds of Annona muricata, Arenga saccharifera, Cocos nucifera, Convolvulus tricolor, and Sophora japonica,” Carbohydrate Research, vol. 20, no. 2, pp. 329–337, 1971.

S. S. Nielsen, “Determination of Moisture Content,” in Food Analysis Laboratory Manual, Food Science Texts Series, pp. 17–27, Springer US, Boston, MA, 2010.

D. Ghosh, J. Griswold, M. Erman, and W. Pangborn, “Structural basis for androgen specificity and oestrogen synthesis in human aromatase,” Nature, vol. 457, no. 7226, pp. 219–223, 2009.

G. H. Peters, L. F. Iversen, S. Branner et al., “Residue 259 Is a Key Determinant of Substrate Specificity of Protein-tyrosine Phosphatases 1B and α,” The Journal of Biological Chemistry, vol. 275, no. 24, pp. 18201–18209, 2000.

Y. Toderika and N. Ferguson, “Canagliflozin: A new class of antidiabetic agent targeting the sodium-glucose cotransporter,” Cardiology in Review, vol. 22, no. 2, pp. 97–104, 2014.

M. Sanford and G. L. Plosker, “Anastrozole: A review of its use in postmenopausal women with early-stage breast cancer,” Drugs, vol. 68, no. 9, pp. 1319–1340, 2008.

S. Rampogu and M. Rampogu Lemuel, “Network Based Approach in the Establishment of the Relationship between Type 2 Diabetes Mellitus and Its Complications at the Molecular Level Coupled with Molecular Docking Mechanism,” BioMed Research International, vol. 2016, Article ID 6068437, pp. 1–6, 2016.

V. Zoete, M. A. Cuendet, A. Grosdidier, and O. Michielin, “SwissParam: a fast force field generation tool for small organic molecules,” Journal of Computational Chemistry, vol. 32, no. 11, pp. 2359–2368, 2011.

S. Rampogu, M. Son, C. Park, H. Kim, J. Suh, and K. Lee, “Sulfonanilide Derivatives in Identifying Novel Aromatase Inhibitors by Applying Docking, Virtual Screening, and MD Simulations Studies,” BioMed Research International, vol. 2017, pp. 1–17, 2017.

S. Rampogu, A. Baek, M. Son et al., “Computational Exploration for Lead Compounds That Can Reverse the Nuclear Morphology in Progeria,” BioMed Research International, vol. 2017, pp. 1–15, 2017.

D. van der Spoel, E. Lindahl, B. Hess, G. Groenhof, A. E. Mark, and H. J. C. Berendsen, “GROMACS: fast, flexible, and free,” Journal of Computational Chemistry, vol. 26, no. 16, pp. 1701–1718, 2005.

B. Hess, H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije, “LINCS: a linear Constraint Solver for molecular simulations,” Journal of Computational Chemistry, vol. 18, no. 12, pp. 1463–1472, 1997.

T. Darden, D. York, and L. Pedersen, “Particle mesh Ewald: an N·log(N) method for Ewald sums in large systems,” The Journal of Chemical Physics, vol. 98, no. 12, pp. 10089–10092, 1993.

H. J. C. Berendsen, J. P. M. Postma, W. F. Van Gunsteren, A. Dinola, and J. R. Haak, “Molecular dynamics with coupling to an external bath,” The Journal of Chemical Physics, vol. 81, no. 8, pp. 3684–3690, 1984.

M. Parrinello and A. Rahman, “Polymorphic transitions in single crystals: a new molecular dynamics method,” Journal of Applied Physics, vol. 52, no. 12, pp. 7182–7190, 1981.

W. Humphrey, A. Dalke, and K. Schulten, “VMD: visual molecular dynamics,” Journal of Molecular Graphics, vol. 14, no. 1, pp. 33–38, 1996.

N. A. Baker, D. Sept, S. Joseph, M. J. Holst, and J. A. McCammon, “Electrostatics of nanosystems: application to microtubules and the ribosome,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 98, no. 18, pp. 10037–10041, 2001.

R. Kumari, R. Kumar, and A. Lynn, “g_mmpbsa—A GROMACS tool for high-throughput MM-PBSA calculations,” Journal of Chemical Information and Modeling, vol. 54, no. 7, pp. 1951–1962, 2014.

M. A. Hossain, K. A. S. AL-Raqmi, Z. H. AL-Mijizy, A. M. Weli, and Q. Al-Riyami, “Study of total phenol, flavonoids contents and phytochemical screening of various leaves crude extracts of locally grown Thymus vulgaris,” Asian Pacific Journal of Tropical Biomedicine, vol. 3, no. 9, pp. 705–710, 2013.

R. Yadav, R. Tiwari, P. Chowdhary, and C. K. Pradhan, “A pharmacognostical monogroaph of Trigonella foenum-graecum seeds,” International Journal of Pharmacy and Pharmaceutical Sciences, vol. 3, pp. 442–445, 2011.

Z. Madar and I. Shomer, “Polysaccharide Composition of a Gel Fraction Derived from Fenugreek and Its Effect on Starch Digestion and Bile Acid Absorption in Rats,” Journal of Agricultural and Food Chemistry, vol. 38, no. 7, pp. 1535–1539, 1990.

M. I. Zafar and F. Gao, “4-Hydroxyisoleucine: A Potential New Treatment for Type 2 Diabetes Mellitus,” BioDrugs, vol. 30, no. 4, pp. 255–262, 2016.

S. D. V. Rampogu, “Role of breast cancer inhibitors on diabetes mellitus- an in silico approach,” Journal of Diabetes and Metabolic Disorders, vol. 14, no. 1, 11 pages, 2015.

H. Ahmadieh and S. T. Azar, “Type 2 Diabetes Mellitus, Oral Diabetic Medications, Insulin Therapy, and Overall Breast Cancer Risk,” ISRN Endocrinology, vol. 2013, pp. 1–8, 2013.

J. E. Thomas, M. Bandara, E. L. Lee, D. Driedger, and S. Acharya, “Biochemical monitoring in fenugreek to develop functional food and medicinal plant variants,” New Biotechnology, vol. 28, no. 2, pp. 110–117, 2011.