Computational screening and molecular dynamics reveal curcumin III and taxifolin as potential thyroid receptor modulators for hypothyroidism therapy

Hypothyroidism is a condition marked by inadequate thyroid hormone production. It is typically treated with Levothyroxine, which, despite its effectiveness, can cause adverse effects on metabolism and the cardiovascular system. In this study, an in silico approach was used to screen phytochemicals from Curcuma longa, Moringa oleifera, and Nigella sativa for their potential to activate the thyroid receptor. A total of 439 compounds were docked against Thyroid Receptor Beta 1 (TRβ1) and Thyrotropin-Releasing Hormone Receptor (TRHR) using AutoDock Vina. Among them, valoneic acid dilactone, curcumin III, taxifolin, luteolin, lophenol, and stigmastanol were identified from the three plants as performing better than others. These compounds were further evaluated with ADMET predictions and molecular dynamics simulations to assess their drug-like properties and stability. Lophenol from Nigella sativa showed the highest binding affinity to TRβ1 (-13.62 kcal/mol) with an inhibition constant of 0.1037 nM, surpassing levothyroxine (-11.60 kcal/mol and 6.63 nM). Also, stigmastanol from Nigella sativa showed the highest binding affinity to TRHR (-8.71 kcal/mol) with an inhibition constant of 415 nM, surpassing levothyroxine (-7.64 kcal/mol and 2530 nM). However, the two compounds exhibited some unfavourable ADMET properties that could be improved through formulation science. Curcumin III and taxifolin demonstrated favourable ADMET properties, including optimal solubility, lipophilicity, and polar surface area, suggesting good oral bioavailability. Molecular dynamics simulations were carried out using Desmond. The results revealed that curcumin III, lophenol, and taxifolin maintained strong stability, with RMSD values of 1.70 ± 0.005 Å, 2.17 ± 0.012 Å, and 2.50 ± 0.014 Å, respectively, for TRβ1, and RMSD values of 3.82 ± 0.014 Å, 4.29 ± 0.019 Å, and 3.31 ± 0.017 Å, respectively, for TRHR. This study identified Curcumin III from Curcuma longa and taxifolin from Moringa oleifera as promising natural compounds for the treatment of hypothyroidism. Further validation through in vitro, in vivo, and ex vivo studies is recommended.

Introduction

Hormones are chemical substances produced by endocrine glands and enter the bloodstream to control several bodily functions. These endocrine glands include the pituitary, adrenal, gonad, thyroid, parathyroid, and pancreatic glands (1). Scientific research has shown that the primary endocrine gland to develop during embryonic development is the thyroid gland (2, 3). Thus, it is regarded as a basic endocrine gland (4). The inability of the thyroid gland to produce adequate thyroid hormones - thyroxine (T4) and triiodothyronine (T3) (which are essential for the metabolic needs of the body) results in the disease condition known as hypothyroidism or under-active thyroid (5). This is brought about by several factors, namely, (a) insufficient stimulation of the thyroid gland by the hypothalamus/pituitary gland, (b) inability of the thyroid gland to receive proper messages from the brain to produce more hormones, (c) thyroid surgery or exposure to radiation (6). Hypothyroidism is a widespread condition with the capacity to culminate in destructive health consequences that affect humans globally (7). This condition alters nearly the entire body system, and it can present as asymptomatic or symptomatic (5). More so, it could be transient or permanent (8). The manifestations and indicators indicative of thyroid impairment are characterized by their nonspecific and non-diagnostic nature, particularly in the initial stages of disease manifestation. (9). If hypothyroidism is detected in its early stages, it is easy to treat, but if proper treatment is neglected, it can trigger life-threatening complications across many organ systems, where the circulatory system is the most strongly investigated target (10). These complications include infertility, hypertension, neuromuscular dysfunction, dyslipidemia, and cognitive impairment (6). The occurrence of hypothyroidism increases with age, particularly in people above 50 years, with a higher occurrence in women than in men. (11) It has been assessed that about 5% of the general population of humans are affected by hypothyroidism, while a further estimation of 5% of the human population is not identified through diagnosis (12).

There are different types of hypothyroidism, including primary, secondary, tertiary, and peripheral types. Primary hypothyroidism is a product of a deficiency of thyroxine, while the deficiency of thyroid-stimulating hormone (TSH) results in secondary hypothyroidism. Peripheral hypothyroidism occurs due to a shortage of thyrotropin-releasing hormone (5). Both the secondary and tertiary hypothyroidism make up central hypothyroidism. Both the central hypothyroidism and peripheral hypothyroidism (extrathyroidal) are uncommon, with only about less than 1% cases (5). In days gone by, a deficiency of iodine was the common cause of hypothyroidism because iodine mainly comes from diet, and the body does not make iodine (13). Nowadays, the major origin of primary hypothyroidism is Hashimoto’s disease, an autoimmune disorder in which the body produces antibodies that harm thyroid tissues, causing the thyroid gland to deteriorate gradually and suppress the production of the thyroid hormone (14). More than 99% of affected patients are diagnosed with primary hypothyroidism (12). Other causes of underactive thyroid include stress, some medications and hormonal fluctuations owing to pregnancy and menopause Signs (15) and symptoms of underactive thyroid include fatigue, depression and apathy, weight gain, constipation, insomnia, hair loss, infertility, anxiety, panic attack, tingling hands and feet, hoarse voice, lowered immunity, cold intolerance, lethargy, low libido, headache, slow healing and so on (8, 9).

Majority of individuals with hypothyroidism will require protracted thyroid hormone therapy (12). Levothyroxine, the recommended drug for the treatment of hypothyroidism, unfortunately has unfriendly backsides, though it normalises the quantity of the thyroid-stimulating hormone (TSH) (6). Excessive or under-treatment with levothyroxine usually has negative side effects on glucose and lipid metabolism and cardiovascular function. It decreases physical capability owing to impaired cardiac output, causing muscle weakness and reducing the proper functioning of the lungs, hampering good growth and development, endangering reproductive function, cognitive function, bone metabolism, and gastrointestinal function (16). And for that reason, there is an accelerating interest in alternative treatments such as the use of bioactive compounds in medicinal plants. These would have the potential to support thyroid function, normalize thyroid hormones and inhibit side effects. Certain bioactive compounds may support patients on thyroid hormone replacement therapy by improving their overall health and reducing the need for higher doses of medication.

This study aims to utilize computational tools to investigate the potential therapeutic effects of bioactive compounds from these three medicinal plants - Moringa oleifera, Curcuma longa, and Nigella sativa that can be used as alternatives to the conventional drug (levothyroxine) with minimal or zero side effects, having the potential to target the Thyroid receptor beta 1 and Thyrotropin-Releasing Hormone receptor.

Moringa oleifera Lam. belonging to the Moringaceae family, sometimes called the ‘miracle tree’ as every part of the plant has nutritional and pharmacological benefits (17, 18). Previous reports showed its medicinal effects in hypertension, diabetes, inflammatory conditions, cancer, osteoporosis and thyroid dysfunctions (19–22). M. oleifera, especially its leaf extract has been shown to regulate thyroid hormones in different animal studies (23; –20, 24). Curcuma longa (F. Zingiberaceae), commonly called turmeric, is a rhizomatous plant with origin in India and cultivated widely in China, West Africa, East Africa and some other tropical areas (25). Tumeric is used in Ayurvedic medicine (India) and Chinese Traditional Medicine (China) for the treatment of a wide range of illness including cough, cancer, diabetes, ulcers inflammation. The dry rhizome is also used as flavouring agent, having a characteristic yellow colour due to the phenolic compounds ‘curcuminoids’ present (25, 26). Nigella sativa (F. Ranunculaceae), an annual herb plant found to grow in different parts of the world though native to Asia, has been referred to as an important medicinal plant. It is commonly referred to as Black seed, and long been used to relieve different ailments as used in Unani, African, Eastern, Arabic, Chinese and Ayurvedic medicine (27). N. sativa parts seeds, oils, and extracts have been reported in studies for treatment of asthma, inflammatory disorders, cancer, arthritis, diabetes, conditions resultant from oxidative stress e.t.c (27–29). Also, a number of animal studies have been conducted to evaluate its potential as therapeutics in hormonal dysfunction including hypothyroidism (30–36).

Furthermore, and more specifically, this study is aimed at identifying bioactive compounds that have therapeutic effects on dual-targets (TRβ1 and TRHR), because there are no medicinal plants previously identified in scientific literatures that specifically target both TRβ1 and TRHR simultaneously for hypothyroidism treatment (37). In addition, no medicinal plants have been reported to specifically target TRHR for hypothyroidism (38), thus alluding to the novelty of this study. Moreso, research into taxifolin's role in hypothyroidism is limited (39) and reported findings on curcumin are primarily based on animal studies, rather than computational. Based on available data, very few natural products with a dual-target modulatory effect have been reported for hypothyroidism therapy (40, 41). This is the insight that this study provides computational empirical evidence for, thus providing a compass for further in vivo or in vitro studies.

Materials and methods

Phytochemical library Generation and Ligand preparation

A library of phytochemicals of Curcuma longa (158), Moringa oleifera (106), and Nigella sativa (175), previously reported in published articles, was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/) and IMPPAT online databases (42, 43) in two-dimensional structure (2D) SDF format. The phytochemicals were converted to pdbqt formats using the Open Babel LigPrep command lines in AutoDock Vina (44). The ligands were protonated at physiological pH (~7.4) using Open Babel. Gasteiger–Marsili partial charges, which were well-matched with the AutoDock Vina scoring function were assigned to the ligands.

Receptor selection, retrieval and preparation

Two proteins involved in hypothyroidism were selected for this study. Three-dimensional (3D) Experimental structures of Thyroid receptor beta 1 (PDB ID: 1NAX) (45) and Thyrotropin-Releasing Hormone receptor (PDB ID: 7X1T) (46) were retrieved from the Protein Data Bank (https://www.rcsb.org) and visualised using PyMOL molecular graphics system (version 2.0, Schrödinger, LLC., New York, NY, USA). Polar hydrogen was added while water molecules, other atoms and the co-crystallised ligands were removed. The proteins were protonated at physiological pH (~7.4) using Open Babel and gasteiger–Marsili partial charges, compatible with the AutoDock Vina scoring function were assigned to the receptors before the molecular docking (44).

Binding site definition and receptor grid generation

The coordinates of active sites on the receptors (receptor grid), where ligand interaction occurs, were defined by using the co-crystallised ligands to specify the dimensions on Autodock (47). The user-specified grid map point was set to 60 X 60 X 60. The XYZ grid coordinates are 4.01 Å, 21.1 Å, 31.86 Å for 1NAX and 98.97Å, 117.55 Å, and 76.52 Å for 7X1T, respectively.

Molecular docking

The study employed the PyMOL/Autodock Plugin tool for docking of the phytochemicals against the proteins (47, 48). The co-crystallised ligands were extracted from each receptor on PyMOL and redocked on Audock Vina, while the other ligands were docked using script where the receptor (protein target) was treated as static (rigid) while the ligands’ rotatable bonds were set to be free (flexible). The exhaustiveness was set at 8, and one pose with the lowest RMSD was generated as output. The 3D protein-ligand interaction was generated using the Discovery Studio Visualizer version 2024. The online tool, PLIP, was also used to determine the interacting amino acid residues (49). The ligands from each plant were ordered based on binding affinity and inhibition constant values. The selected ligands were subjected to molecular dynamics simulation.

Drug likeness and in silico pharmacokinetic profile

Selected pharmacokinetics (ADMET) and drug-likeness properties were conducted using the canonical smile notations of selected lead compounds in ADMET Lab 3.0 (50). These in silico pharmacokinetic and drug-likeness properties included hydrogen bond donor/acceptor, molecular weight, Lipinski’s rule of five, Blood-brain barrier penetration, partition coefficient, etc.

Molecular dynamics simulation and trajectory analysis

Molecular dynamics (MD) simulation and Trajectory analysis were carried out to validate the docking results obtained earlier and evaluate the stability of lead complexes using the previously described method reported in (51). The docked complexes of lead compounds, from each plant, were subjected to MD simulation alongside the standard drug, levothyroxine, and co-crystallized ligands on each receptor. The system was prepared using the TIP3P solvent model of the System builder module of the Maestro Schrödinger suite (52) while each complex was bound in an orthorhombic box with 10 X 10 X 10 Å dimensions, minimized, overall charges neutralized by adding sodium and chloride ions, and the salt concentration was set to 0.15M, mimicking a physiologic state. Thereafter, each system was allowed to undergo a 200ns simulation in an NPT ensemble at 300K and 1 atm using the OPLS3e force field, following the standard protocols of Desmond simulation package of Schrödinger LLC (53). The data outputs of the simulations, including root mean square deviation (RMSD) and radius of gyration (rGyr) were plotted using Origin version 6.0 (54) while statistical analyses were done using GraphPad Prism version 9.0.0 for Windows, GraphPad Software, San Diego, CA, USA (www.graphpad.com, accessed on 1st March 2025).

Binding free energy calculations

The binding free energy was calculated using the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) approach for the molecular dynamic simulation trajectory analysis. (55, 56) This was done using the Prime module in the Maestro environment. In the MM/GBSA methods, the free binding free energy is calculated thus:

$\begin{array}{ll}G=E_{bind}+ E_{el}+E_{vdw}+G_{pol}+G_{np}-TS& (1)\end{array}$

$\begin{array}{ll}\Delta G_{bind}=\Delta G_{complex}–\Delta G_{protein}–\Delta G_{ligand }& (2)\end{array}$

$\begin{array}{ll}\Delta G_{bind}=\Delta H–T\Delta S& (3)\end{array}$

$\begin{array}{ll}\Delta G_{bind}={\Delta }_{gas}+\Delta G_{solv}–T\Delta S& (4)\end{array}$

$\begin{array}{ll}\Delta E_{gas}=\Delta E_{int}+\Delta E_{vdw}+\Delta E_{ele}& (5)\end{array}$

$\begin{array}{ll}\Delta G_{sol}=\Delta G_{GB}+\Delta G_{SA}& (6)\end{array}$

$\begin{array}{ll}{G}_{SA}{=}_{\gamma }SASA& (7)\end{array}$

Where the first three terms in Equation 1 represent energy from bonded, electrostatic and van der Waals interactions respectively while G_pol,and G_np represents polar and non-polar contributions to solvation free energies. ΔE_gas, ΔG_solv, and – TΔS denotes changes in gas phase energy, free solvation energy and total entropy upon ligand binding to the receptor as T is temperature of the systems. ΔE_gasrepresents change in gas phase internal energy which includes changes in internal energies ΔE_int(bond, angle and dihedral energies), van der Waals energy ΔE_vdw and electrostatic (Coulomb) energies ΔE_eleMeanwhile ΔG_solis calculated the summation of the polar energy contribution: electrostatic solvation energy ΔG_GBand nonpolar solvation energy contribution ΔG_SA usually estimated using the solvent-accessible surface area (SASA) between the solute and the solvent (57–59).

The plant materials used in this study were Curcuma longa, Moringa oleifera, and Nigella sativa, while the methods engaged in the study involved the generation and preparation of proteins and ligands, generation of the grid site, molecular docking of the prepared ligands into the active sites, pharmacokinetics analysis of the hit compounds, molecular dynamics simulation of the protein-ligand complex of the hit compounds and the binding energy calculations from the MD trajectories of the hit compounds.

Result

Binding affinities and inhibition constant

In this study, the therapeutic of phytochemicals from Curcuma longa, Moringa oleifera, and Nigella sativa as agonists at the binding sites of Thyroid receptor beta 1 (PDB ID: 1NAX) and Thyrotropin-Releasing Hormone receptor (PDB ID: 7X1T) were assessed. On 1NAX, lophenol (T5) showed the highest binding affinity of -13.62 kcal/mol, while taxifolin is the least performing hit on the receptor with a docking score of -8.28 kcal/mol. The binding affinity of valoneic acid dilactone (T1) is -10 kcal/mol, whereas the co-crystallized ligand, IH5 {3, 5-dichloro-4-[4-hydroxy-3-(propan-2-yl)phenoxy]phenyl}acetic acid) and standard drug, levothyroxine, are -11.44 and -11.60 kcal/mol.

On 7X1T, the co-crystallized ligand has a binding affinity of -8.67kcal/mol, but lophenol, valoneic acid dilactone, and taxifolin have decreasing binding affinities of -8.07, -7.91, and -7.29 kcal/mol, respectively. Figure 1 depicts the heatmaps of the screened phytochemicals with an arrow (red) on the potential leads from each plant. The top two (2) lead compounds with lower binding affinities on both receptors for each plant were identified as shown in Table 1.

Figure 1

Figure 1. Heatmap representation of the molecular docking scores of phytochemicals identified in selected medicinal plants against the studied target proteins. (A) Moringa oleifera, (B) Nigella sativa, and (C) Curcuma longa.

Table 1

Table 1. The binding affinities (docking scores) and the inhibition constants of the selected compounds.

Docking validation

The validation to check the accuracy of the docking procedure was done by re-docking the co-crystallised ligands into the binding sites of the receptor. The low root-mean-square deviation (RMSD) values of 0.000 Å (1NAX) and 0.000Å (7X1T) established the consistency of the grid parameters. Superimposed ligand structures are shown in Figure 2.

Figure 2

Figure 2. Superimposed structures of the co-crystallised ligands in their co-crystallised (magenta) and re-docked poses (yellow) at the active site of the two receptors 1NAX and 7X1T.

Ligand interaction network

Common key network interactions of the hit compounds, valoneic acid, taxifolin, and lophenol with the amino residues of 1NAX include Arg320, Arg282, Arg316, Leu330, and His435. While hydrogen bonds and hydrophobic interactions are the major interaction networks, some formed water bridges, while others formed salt bridges (Table 2; Figures 3a, b).

Table 2

Table 2. Ligand interaction types and amino acid residues at 1NAX AND 7X1T.

Figure 3

Figure 3. (A) The interaction networks between phytochemicals (in magenta) and receptor 1NAX. (B) The interaction networks between phytochemicals (in magenta) and receptor 7X1T.

In silico ADMET properties and drug likeness

The ADMET properties of hit compounds and standards were evaluated using ADMET Lab 3.0. Valoneic acid dilactone and Taltirelin violated the Lipinski rule of 5 because of the number of hydrogen bond acceptors, which exceeded 10. Details are given in Table 3 and Figure 4.

Table 3

Table 3. The in silico ADMET properties of top leads.

Figure 4

Figure 4. The radar chart of selected physicochemical properties of hits from M. oleifera, N. sativa, and C. longa phytochemicals.

Molecular dynamics simulation analysis

The best ligand from each plant was selected based on the low binding affinity, key molecular interactions, and ADMET properties. The two major targets that were of concern in hypothyroidism conditions are Thyroid receptor beta 1 and Thyrotropin-Releasing Hormone receptor.

For the simulation period of 200ns, the protein and ligand properties, including protein-ligand root mean square deviation (PL_RMSD), protein root mean square fluctuation (P_RMSF), radius of gyration (rGyr), molecular surface area (MolSA), polar surface area (PSA), solvent-accessibility surface area (SASA) were analysed with average values calculated and recorded as shown in Table 4 and the plots were represented in Figures 5a, b.

Table 4

Table 4. The simulation results of the complexes of the hit compounds on the receptors.

Figure 5

Figure 5. (a) The MD simulation results for hit compounds and co-crystallised ligands on the receptor 1NAX. Note: Plots (A) Protein-Ligand Root Mean Square Deviation (PL_RMSD); (B) Protein Root Mean Square Fluctuation (P_RMSF); (C) Radius of Gyration (rGyr); (D) Molecular Surface Area (MolSA); (E) Polar Surface Area (PSA); (F) Solvent-Accessibility Surface Area (SASA). (b) Shows MD simulation findings for hit compounds on the receptor 7X1T Note: Plots (A) Protein-Ligand Root Mean Square Deviation (PL_RMSD); (B) Protein Root Mean Square Fluctuation (P_RMSF); (C) Radius of Gyration (rGyr); (D) Molecular Surface Area (MolSA); (E) Polar Surface Area (PSA); (F) Solvent-Accessibility Surface Area (SASA).

Protein and ligand properties of hit compounds across the 1NAX receptors are shown in Table 5 and Figure 4a. An RMSD value measures the system stability throughout the simulation time. The co-crystallized ligand (IH5) exhibited the lowest RMSD value, closely similar to curcumin III RMSD value of 1.70 ± 0.006Å and 1.70 ± 0.005 Å, respectively. The lead compounds, taxifolin and lophenol, have close values with the standard drug, levothyroxine (2.50 ± 0.014 Å, 2.17 ± 0.012 Å, and 2.10 ± 0.012 Å, respectively). These values, though higher than that of the co-crystallized ligand, still fall within the acceptable RMSD fluctuation range of 1 and 3Å. While the leads, taxifolin and lophenol, displayed fluctuations within the first 40ns, curcumin III was quite stable during the simulation time.

Table 5

Table 5. The MM/GBSA binding free energy estimation of the protein–ligand complex using the MD trajectories of the hit compounds.

The stability of individual amino acid residues at the active sites of a protein is measured by RMSF analysis (60). The average RMSF value from the lowest, 0.92 ± 0.043 Å of IH5, similar to that of curcumin III, 0.924 ± 0.035 Å to the highest value exhibited by taxifolin (1.09 ± 0.049 Å), all indicated minimal fluctu ations in the active binding site region, hence revealing stable ligand-binding. In Figure 4a (B), unstable fluctuations of the amino acid residues were observed in the different ligand systems within the regions 40-55, where lophenol showed the highest fluctuation, while curcumin III has the lowest.

The plots of radius of gyration (rGyr), which analyse protein stability and compactness, revealed that the compounds were stable. The average value of lophenol is the highest (5.13 ± 0.004 Å) and taxifolin exhibited the lowest rGyr (3.74 ± 0.001 Å). For Curcumin III and Lophenol, their Gyr values are quite higher than the values of the standard drug and the co-crystallized ligand (5.08 ± 0.008 Å, 5.13 ± 0.004 Å, 4.40 ± 0.001 Å, 4.21 ± 0.002 Å, respectively).

On the 7X1T, as shown in Figure 4b, PL-RMSD values of co-crystallized ligand (taltirelin) and taxifolin are quite similar, and the highest value recorded (4.29 Å), while lophenol, levothyroxine, and curcumin III were 3.31 ± 0.017 Å, 3.67 ± 0.012 Å, and 3.82 ± 0.014 Å, respectively. These values exceeded 3.0 Å, showing that the ligands are not very stable on this receptor. The RMSF values of curcumin III, taxifolin, and lophenol are lower (1.58 ± 0.067 Å, 1.52 ± 0.051 Å, 1.59 ± 0.056 Å) than that of the co-crystallized ligand (1.87 ± 0.079 Å). These values are much lower than levothyroxine (2.09 ± 0.070 Å). The plots of rGyr showed minimal fluctuations; hence, the compounds were stable in the binding pocket.

The values of PSA and MolSA are relatively similar for each of the ligands, but the values of SASA on the receptors. For instance, curcumin III’s PSA value on 1NAX and 7X1T are 180.7 ± 0.079 Å and 181 ± 0.125 Å, respectively, its MolSA is 311 ± 0.0702 Å and 313 ± 0.066 Å, respectively, but its SASA values are 10.3 ± 0.232 Å and 143 ± 1.84 Å, respectively.

The MM/GBSA binding free energy calculations of the complex, as shown in Table 5, provides a broad overview of the energy breakdown including van der Waals interaction energy (ΔE_vdW) as well as binding free energy (ΔG_bind) for Curcumin III, Taxifolin, and Lophenol compared to the co-crystallized compounds for the two different protein targets 1NAX and 7X1T. The positive ΔE_vdW values (72 ± 103 kcal/mol) and exceptionally high positive ΔG_bind values (574 ± 612 kcal/mol) for Lophenol showed the unfavourable and unstable binding to 7X1T. Curcumin III showed optimal and stable binding with the two targets, having ΔE_vdW -52 ± 0.29 kcal/mol, ΔG_bind -80 ± 0.69 kcal/mol, on 1NAX and ΔE_vdW -40 ± 0.30, kcal/mol, ΔG_bind -55 ± 0.56 kcal/mol on 7X1T. Taxifolin has its ΔE_vdW -44 ± 0.23 kcal/mol, ΔG_bind -67± 0.40 kcal/mol and 1NAX, and ΔE_vdW -27 ± 0.34 kcal/mol, ΔG_bind -35 ± 0.64 kcal/mol on 7X1T.

The phytochemicals of Curcuma longa, Moringa oleifera, and Nigella sativa were screened as potential agonists of the TRβ1 (1NAX) and TRHR (7X1T) receptors where curcumin III was found to be the most promising as it formed stable interactions on both receptors, while Taxifolin showed moderate binding properties in the computational study discussed extensively.

Discussion

In silico analysis has provided great advancement in the discovery of drugs and the study of diseases at the molecular level using computational skills. It performs a key function in the comprehension of protein–ligand interaction and in the prediction of drug-likeness properties of any compound, which is crucial for bench-to-bedside drug development (61). Hypothyroidism treatment with conventional drugs, namely, levothyroxine and liothyronine, is accompanied by side effects. Therefore, the discovery of novel anti-hypothyroidism drugs has become crucial, birthing the invaluable need for evaluation of small molecules as a substitute for hypothyroidism treatment. One of the novel options in the last decade for the enhancement and management of several disease conditions has been the use of bioactive plant compounds, so that there can be minimal side effects of drugs on humans (62). Studies on the treatment of hypothyroidism using medicinal plants are limited, (41). This present study focuses on investigating the activation potential of the hit bioactive compounds across board from Moringa oleifera, Curcuma longa, and Nigella sativa against the active site of the clinical targets. Thyroid receptor beta 1 and Thyrotropin-Releasing Hormone receptor. Molecular docking, molecular dynamics simulation, ADMET properties, among others, are a few of the methods employed.

Systematically, valoneic acid dilactone, curcumin III, taxifolin, luteolin, lophenol and stigmastanol interacted with basic residues such as Arg282, Asn331, Arg320, His435, and Arg316, Gly332, Thr329, Ala279, similar to known activators of Thyroid receptor beta 1 (TRβ). The co-crystallised ligand of the TRβ 1 N-{3-[(4s, 5r)-2-Amino-5-Fluoro-4-(Fluoromethyl)-5, 6-Dihydro-4h-1, 3-Oxazin-4-Yl]-4-Fluorophenyl}-5-Cyanopyridine-2-Carboxamide formed hydrogen bonds with Asn331 and His435 amino acid at the pocket of the active site (45). This finding agrees with prior scholarly work on the structural analysis of selective agonists targeting thyroid hormone receptor β using 3D-QSAR and molecular docking methodologies. It has been established that the specific binding pocket in which the agonists bind contains both polar (Arg282, Ser314, Arg316, Asn331) and nonpolar (Phe269, Phe272, Ile275, Ile276, Ala279, Met310, Ile312, Met313, Ala317, Leu330, Leu346, His435, Phe459) amino acids (63). The presence of a hydrophobic interaction observed between the co-crystallised ligand and the amino residues at the pocket site of the TRβ 1 (N-{3-[(4s, 5r)-2-Amino-5-Fluoro-4-(Fluoromethyl)-5, 6-Dihydro-4h-1, 3-Oxazin-4-Yl]-4-Fluorophenyl}-5-Cyanopyridine-2-Carboxamide) was absent in the interaction occurring between the standard drug (levothyroxine) and the binding site residues. The drug's pharmacokinetics are profoundly influenced by hydrophobic interactions, which are indispensable for the drug's engagement with its biological targets. The affinity and specificity of the drug-receptor interaction depend on hydrophobic interactions, which are the primary features stabilising the drug within this binding site (64). Hence, this elucidates the underlying reason for the adverse effects of Levothyroxine.

For taxifolin, the three hydrogen bonds to Arg282, Arg316 and Gly332, the water bridge formation with Asn331, and the hydrophobic interaction with Phe272 collectively improved binding affinity by making best use of corresponding electrostatic and hydrophobic contacts in the TRβ pocket. Asn331, through a structured hydrogen bonding network, facilitates isoform-specific binding, which is a critical switch for β-selective thyromimetics. The promising stabilisation of the active conformation, facilitated by hydrogen bond and hydrophobic interaction, etc., enables the agonist to activate the receptor and initiate drug action in TRβ, thereby aiding in TRβ-targeted drug development for hypothyroidism (64). This binding mode implies that taxifolin has the potential to be developed to ameliorate hypothyroidism symptoms largely through TRβ activation, contingent upon in vitro and in vivo data supporting that this docking pattern relates to functional β-selective agonism.

Curcumin III also has two well-positioned oxygen atoms in an aromatic ring that can form a bidentate hydrogen bonding network with amino acids,  Arg316 and Arg320. It participates in a π–π interaction with Phe272. Together, these interactions increase ligand affinity and drive a receptor conformation that is agonistic only in TRβ, as long as the pocket geometry in TRβ allows this arrangement to be accommodated. This structural correspondence provides the most elementary mechanism for selective activation of TRβ and provides clear directions toward its rational optimisation and experimental verification.

Also, the basic residues akin to activating thyrotropin-releasing hormone receptors (TRHR) include Gln105, Tyr106, Tyr282, Asn298, Arg306, Arg185, Tyr181, and Tyr310 (46), whereas the co-crystallized ligand- taltirelin forms a hydrogen bond with Gln105, Asn289, Arg306, Tyr181, and Arg185 and hydrophobic interaction with Leu302. A comparable interaction was also noted for valoneic acid dilactone, curcumin III, taxifolin, luteolin, lophenol, and stigmastanol within the active site. The likelihood of Curcumin III to exhibit selective TRHR stimulation is due to its β-diketone and phenolic oxygens forming a firm, multi-point hydrogen bond network with Arg316 and Tyr106, while its aromatic rings take part in π-π stacking with Phe272. These amino acids arrange the ligand optimally to effectively stabilise the active conformation of TRHR, resulting in selective activation of the receptor. This categorises Curcumin III as a promising compound for developing therapeutic agents aimed at TRHR within the context of hypothyroidism because activation of TRHR promotes physiological TSH release and reinstates normal thyroid hormone synthesis. This implies that the potency of the interaction increases due to the numerous hydrogen bond interactions, which facilitate unique hydrogen-bonded interactions when compared to levothyroxine.

Taxifolin possesses the capability to selectively activate TRHR owing to the presence of phenolic hydroxyl and carbonyl groups, allowing a multi-point H-bond network with Gln105, Tyr106, Tyr282, and Arg306. The above residues offer a unique polar-aromatic microenvironment within the receptor’s binding pocket. The contact and interaction with Gln105 help to stabilise the molecule; a dual-pronged contact by the Tyr106-Tyr282 module secures, in an effective way, the flavanonol ring system through H- and π-interactions, and a solid directional H-bond by Arg306 helps to anchor the receptor firmly into its active conformation. The described topology together culminates in TRHR selectivity, thus indicating that Taxifolin or its derivatives may represent a new class of upstream therapeutics for managing hypothyroidism by triggering normal TSH secretion through targeted activation of TRHR.

From the ADMET results, two of our lead compounds, Curcumin III and taxifolin from Curcuma longa and Moringa oleifera, respectively, gave good results in relation to how effectively they would be absorbed, distributed, metabolized, and excreted. From the literature, the LogS result above (-5) indicates an overall efficacy of the compounds. From the results, Valoneic acid dilactone (-3.14), Curcumin III (-3.34), taxifolin (-3.35), and luteolin (-4.02) are within the accepted range, while lophenol (-7.78) and stigmastanol (-8.22) could not suffice. Their logP results also show that valoneic acid dilactone, curcumin III, taxifolin, and luteolin have their values less than 3.0. This shows a better result than the standard drug and the co-crystallized ligand (Taltirelin), whose liposolubility is a little higher than 3.0 and less than -1, respectively. This result suggests that they can penetrate the lipid membrane, and since any compound that will be effective for treating hypothyroidism must be able to penetrate the lipid membrane (65, 66). The selected compounds in this study give a better result. The high liposolubility of lophenol and stigmastanol makes it difficult to penetrate the lipid membrane. TPSA range for an effective hypothyroidism drug is between 60-140 Å (67). Fortunately, Curcumin III has 74.6 Ų, taxifolin 127.5 Ų, luteolin 111.1Ų, levothyroxine 92.78 Å^2, and co-crystallized compound (1H5) 113.4Ų compared to the co-crystallized compound (taltirelin) with 170.6 Å. All the selected compounds are more likely to have good oral bioavailability, having followed the Lipinski rule of five, excluding valoneic acid dilactone and taltirelin (68). The ADMET result, however, confirms that curcumin and taxifolin will give better results compared to levothyroxine.

Additionally, the stability of protein-activator complexes was estimated by performing Molecular Dynamics (MD) simulations. The principal tool in computational chemistry to acquire a thorough report on the fluctuations and conformational shifts of proteins and ligands in a time-bound approach is MD simulation (69). The results of the molecular dynamic simulation, particularly, the RMSD, RMSF, protein-ligand contact, and ligand properties show that the choice compounds are quite stable on both receptors as their RMSD values are within the accepted range of 2-3Å, Curcumin III (1.70 ± 0.005), taxifolin (2.50 ± 0.014 Å), and lophenol (2.17 ± 0.012Å) (70) The P_RMSF results (< 2 Å) are indicative of a stabilized protein. The stabilised binding mode implies that curcumin (0.924 ± 0.035 Å), taxifolin (1.09 ± 0.049 Å), and lophenol (1.06 ± 0.058 Å) provide strong interaction with TRβ 1 and TRHR, thereby improving its activator potential. The low SASA values, high PSA values, and rGyr results of the curcumin and taxifolin on the two receptors suggest the need for optimisation of the compound to enhance its bioavailability.

The MMGBSA results reveal that on the 1NAX protein, all three ligands produce positive van der Waals properties, with negative ΔE_vdW values. It is appropriate to note that van der Waals activity in Lophenol is the highest (-69 ± 0.34 kcal/mol) compared with Curcumin III (-52 ± 0.29 kcal/mol), Taxifolin (-44 ± 0.23 kcal/mol), and the co-crystallized ligand (-52 ± 0.22 kcal/mol). This trend is evident in the analysis of binding free energy, too. Lophenol has the best ΔG_bind (−96 ± 0.62 kcal/mol), which indicates that Lophenol has the best binding ability for the 1NAX target molecule. Curcumin III has a ΔG_bind of (−80 ± 0.69 kcal/mol), which is better when compared to the co-crystallized molecule with a ΔG_bind of (−75 ± 0.45 kcal/mol). However, compared to the co-crystallized molecule, Taxifolin has a relatively weaker binding free energy. The binding ability of Curcumin III shows a more stable binding when compared to the binding ability of the co-crystallized ligand (−52 ± 0.32 kcal/mol).

It is significant to note that Lophenol has profoundly unfavourable and unstable binding to 7X1T, as made clear by its positive ΔE_vdW values (72 ± 103 kcal/mol) and its exceptionally high positive ΔGbind values (574 ± 612 kcal/mol). Additionally, the large error bars reflect a lack of convergence of the interactions that occurred during the simulation due to steric issues or the inability of Lophenol to adopt a stable bindable conformation in the 7X1T bindable site.

Several compounds, including Taxifolin and Curcumin III, exhibited interactions similar to known receptor agonists. 2-D Ligand interactions revealed that hydrogen bonding and hydrophobic interactions show selective activation of TRβ1 and TRHR. ADMET results ranked Curcumin III and Taxifolin higher because of their solubility, membrane permeability, and oral bioavailability over the other phytocompounds and levothyroxine. Molecular dynamics and MMGBSA analyses confirmed the stability of the binding of Curcumin III and Taxifolin, making them leading dual-target compounds.

Conclusion

The present study highlights the potential of curcumin III from Curcuma longa and taxifolin from Moringa oleifera as the most effective activators of the selected dual target of hypothyroidism - thyroid receptor beta 1 and thyrotropin-releasing hormone receptors. The dynamics and steadiness of curcumin III and taxifolin complexes have been proven by carrying out molecular dynamics simulations. Besides, the physicochemical and ADMET properties of curcumin III and taxifolin revealed that they have almost all the properties of a drug-like molecule. The MMGBSA results showed that lophenol has exceptional binding affinity for 1NAX, but it does not bind favorably with 7X1T, which clearly signifies selectiveness. This disqualifies it from being a potential dual target, which this study aims at. Curcumin III shows optimal and stable binding with all the targets and performs better than the ligand used in the co-crystallization process, clearly signifying its potential to be a universal binder. On the other hand, taxifolin shows relatively weaker binding despite its stable nature when compared with the other ligands. The binding stability in taxifolin is an important strength, which clearly implies strong compatibility with the receptor pocket. The relatively low affinity could be explained based on an inadequately addressed hydrophobic/disperse interaction, which could be improved by hydrophobic site modification, rigidity improvement, and site-directed hybridization. By optimizing taxifolin, it has the potential to be developed into a highly effective and selective lead compound for the corresponding molecular targets for hypothyroidism. Nevertheless, the outcomes of this study necessitate in vitro, in vivo and ex vivo endorsement for the verified natural compounds beforehand clinical procedure.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Author contributions

BO: Conceptualization, Methodology, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing, Data curation, Software, Validation. AA: Conceptualization, Methodology, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing, Data curation, Software, Validation. SA: Conceptualization, Methodology, Project administration, Supervision, Visualization, Writing – original draft, Writing – review & editing, Data curation, Software, Validation. OF: Investigation, Methodology, Writing – original draft. OSO: Conceptualization, Investigation, Project administration, Writing – review & editing. OMO: Investigation, Writing – original draft. AMA: Methodology, Writing – review & editing. ABA: Investigation, Methodology, Writing – original draft. OI: Investigation, Supervision, Writing – original draft. AI: Conceptualization, Investigation, Writing – original draft. FO: Investigation, Methodology, Supervision, Writing – original draft. BA: Conceptualization, Supervision, Writing – review & editing.

Funding

The author(s) declared that financial support was not received for this work and/or its publication.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

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

References

1. Hosseinzadeh H. Endocrine system. In: Fundamentals of medicine for biomedical engineering. Springer Nature Switzerland, Cham (2025). p. 409–33. doi: 10.1007/978-3-031-74085-5_12

Crossref Full Text | Google Scholar

2. Benvenga S, Tuccari G, Ieni A, and Vita R. Thyroid gland: anatomy and physiology. In: Encyclopedia of endocrine diseases. New York: Elsevier (2018). p. 382–90. doi: 10.1016/B978-0-12-801238-3.96022-7

Crossref Full Text | Google Scholar

3. Rosen RD and Sapra A. Embryology, Thyroid. StatPearls. (2023). Available online at: https://www.ncbi.nlm.nih.gov/books/NBK551611/ (Accessed August 12, 2025).

Google Scholar

4. Bhalerao SD and Varadharajulu G. Awareness of role of physiotherapy in women with hypothyroidism. J Ecophysiology Occup Health. (2024) 24:117–21. doi: 10.18311/jeoh/2024/36335

Crossref Full Text | Google Scholar

5. Zamwar UM and Muneshwar KN. Epidemiology, types, causes, clinical presentation, diagnosis, and treatment of hypothyroidism. Cureus. (2023) 15. doi: 10.7759/cureus.46241

PubMed Abstract | Crossref Full Text | Google Scholar

6. Gaitonde DY, Rowley KD, and Sweeney LB. Hypothyroidism: an update. Am Fam Physician. (2012) 86:244–51. doi: 10.1080/20786204.2012.10874256

PubMed Abstract | Crossref Full Text | Google Scholar

7. Olanrewaju OA, Asghar R, Makwana S, Yahya M, Kumar N, Khawar MH, et al. Thyroid and its ripple effect: impact on cardiac structure, function, and outcomes. Cureus. (2024) 16:1–13. doi: 10.7759/cureus.51574

PubMed Abstract | Crossref Full Text | Google Scholar

8. Almandoz JP and Gharib H. Hypothyroidism: etiology, diagnosis, and management. Med Clinics North America. (2012) 96:203–21. doi: 10.1016/j.mcna.2012.01.005

PubMed Abstract | Crossref Full Text | Google Scholar

9. Wilson SA. Hypothyroidea. 1290–1290: Encyclopedia of Clinical Neuropsychology (2011). doi: 10.1007/978-0-387-79948-3_3711

Crossref Full Text | Google Scholar

10. Chaker L, Razvi S, Bensenor IM, Azizi F, Pearce EN, and Peeters RP. Hypothyroidism. Nat Rev Dis Primers. (2022) 8:30. doi: 10.1038/s41572-022-00357-7

PubMed Abstract | Crossref Full Text | Google Scholar

11. McLeod DSA and Cooper DS. The incidence and prevalence of thyroid autoimmunity. Endocrine. (2012) 42:252–65. doi: 10.1007/s12020-012-9703-2

PubMed Abstract | Crossref Full Text | Google Scholar

12. Chiovato L, Magri F, and Carlé A. Hypothyroidism in context: where we’ve been and where we’re going. Adv Ther. (2019) 36:47–58. doi: 10.1007/s12325-019-01080-8

PubMed Abstract | Crossref Full Text | Google Scholar

13. Zimmermann MB and Boelaert K. Iodine deficiency and thyroid disorders. Lancet Diabetes Endocrinol. (2015) 3:286–95. doi: 10.1016/S2213-8587(14)70225-6

PubMed Abstract | Crossref Full Text | Google Scholar

14. Carkit S. Hashimoto’s thyroiditis. In: Hypothyroidism - causes, screening and therapeutic approaches. London, United Kingdom: IntechOpen (2024). doi: 10.5772/intechopen.1005431

Crossref Full Text | Google Scholar

15. Roop JK. Hormone imbalance—A cause for concern in women. Res J Life Sci Bioinform Pharm Chem Sci. (2018) 4:237–51. doi: 10.26479/2018.0402.18

Crossref Full Text | Google Scholar

16. Feldt-Rasmussen U, Effraimidis G, Bliddal S, and Klose M. Consequences of undertreatment of hypothyroidism. . Endocrine. (2024) 84:301–8. doi: 10.1007/s12020-023-03460-1

PubMed Abstract | Crossref Full Text | Google Scholar

17. Sagona WCJ, Chirwa PW, and Sajidu SM. The miracle mix of Moringa: Status of Moringa research and development in Malawi. South Afr J Bot. (2020) 129:138–45. doi: 10.1016/J.SAJB.2019.03.021

Crossref Full Text | Google Scholar

18. Pareek A, Pant M, Gupta MM, Kashania P, Ratan Y, Jain V, et al. Moringa oleifera: an updated comprehensive review of its pharmacological activities, ethnomedicinal, phytopharmaceutical formulation, clinical, phytochemical, and toxicological aspects. Int J Mol Sci. (2023) 24:2098. doi: 10.3390/IJMS24032098

PubMed Abstract | Crossref Full Text | Google Scholar

19. Kaur A, Verma SK, Kalsi S, and Neha N. Hypothyroidism: management based on ayurvedic and modern therapeutic perspective. Int J Pharmaceutics Drug Anal. (2016) 4:281–8.

Google Scholar

20. Alam MA, Quamri MA, and Haider N. Efficacy and safety of Barg-e-Sahajna (Moringa oleifera Lam.) in primary hypothyroidism. Drug Metab Pers Ther. (2022) 37:21–6. doi: 10.1515/DMPT-2021-0136

PubMed Abstract | Crossref Full Text | Google Scholar

21. Abdelazim AM, Afifi M, Abu-Alghayth MH, and Alkadri DH. Moringa oleifera: recent insights for its biochemical and medicinal applications. J Food Biochem. (2024) 2024:1270903. doi: 10.1155/2024/1270903

Crossref Full Text | Google Scholar

22. Camilleri E and Blundell R. A comprehensive review of the phytochemicals, health benefits, pharmacological safety and medicinal prospects of Moringa oleifera. Heliyon. (2024) 10:e27807. doi: 10.1016/j.heliyon.2024.e27807

PubMed Abstract | Crossref Full Text | Google Scholar

23. Tahiliani P and Kar A. Role of moringa oleifera leaf extract in the regulation of thyroid hormone status in adult male and female rats. . Pharmacol Res. (2000) 41:319–23. doi: 10.1006/PHRS.1999.0587

PubMed Abstract | Crossref Full Text | Google Scholar

24. Hadidy AAA and Dawood ST. The effect of moringa oleifera leaves powder on some hormones to prevent the development of experimental hypothyroidism in rabbits. Indian J Forensic Med Toxicol. (2021) 15:1040–7. doi: 10.37506/IJFMT.V15I3.15454

Crossref Full Text | Google Scholar

25. Iweala EJ, Uche ME, Dike ED, Etumnu LR, Dokunmu TM, Oluwapelumi AE, et al. Curcuma longa (Turmeric): Ethnomedicinal uses, phytochemistry, pharmacological activities and toxicity profiles—A review. Pharmacol Res - Modern Chin Med. (2023) 6:100222. doi: 10.1016/J.PRMCM.2023.100222

Crossref Full Text | Google Scholar

26. Fuloria S, Mehta J, Chandel A, Sekar M, Rani NNIM, Begum MY, et al. A comprehensive review on the therapeutic potential of curcuma longa linn. In: Relation to its major active constituent curcumin. Front Pharmacol. (2022) 13:820806. doi: 10.3389/fphar.2022.820806

PubMed Abstract | Crossref Full Text | Google Scholar

27. Majeed A, Muhammad Z, Ahmad H, Rehmanullah Hayat SSS, Inayat N, and Siyar S. Nigella sativa L.: Uses in traditional and contemporary medicines – An overview. Acta Ecologica Sin. (2021) 41:253–8. doi: 10.1016/J.CHNAES.2020.02.001

Crossref Full Text | Google Scholar

28. Kooti W, Hasanzadeh-Noohi Z, Sharafi-Ahvazi N, Asadi-Samani M, and Ashtary-Larky D. Phytochemistry, pharmacology, and therapeutic uses of black seed (Nigella sativa). Chin J Nat Med. (2016) 14:732–45. doi: 10.1016/S1875-5364(16)30088-7

PubMed Abstract | Crossref Full Text | Google Scholar

29. Dabeer S, Rather MA, Rasool S, Rehman MU, Alshahrani S, Jahan S, et al. History and traditional uses of black seeds (Nigella sativa). In: Black seeds (Nigella sativa): pharmacological and therapeutic applications. Elsevier (2022). p. 1–28. doi: 10.1016/B978-0-12-824462-3.00016-0

Crossref Full Text | Google Scholar

30. Tajmiri S, Farhangi MA, and Dehghan P. Nigella Sativa treatment and serum concentrations of thyroid hormones, transforming growth factor β (TGF-β) and interleukin 23 (IL-23) in patients with Hashimoto’s Thyroiditis. Eur J Integr Med. (2016) 8:576–80. doi: 10.1016/J.EUJIM.2016.03.003

Crossref Full Text | Google Scholar

31. Beheshti F, Hosseini M, Shafei MN, Soukhtanloo M, Ghasemi S, Vafaee F, et al. The effects of Nigella sativa extract on hypothyroidism-associated learning and memory impairment during neonatal and juvenile growth in rats. Nutr Neurosci. (2017) 20:49–59. doi: 10.1179/1476830514Y.0000000144

PubMed Abstract | Crossref Full Text | Google Scholar

32. Pakdel R, Hadjzadeh H, Sadegh MM, Hosseini M, Emami B, and Hadjzadeh M-AR. The Effects of Hydroalcoholic Extract of Nigella sativa Seeds on Serum Estradiol and Prolactin Levels and obstetric Criteria due to Hypothyroidism in Rat. Adv BioMed Res. (2017) 6:166. doi: 10.4103/2277-9175.221860

PubMed Abstract | Crossref Full Text | Google Scholar

33. Hosseini M, Ghasemi S, Hadjzadeh MAR, Ghorbani A, Aghili S, Aghaei A, et al. Administration of Nigella sativa during neonatal and juvenile growth period improved liver function of propylthiouracil-induced hypothyroid rats. J Maternal-Fetal Neonatal Med. (2020) 33:718–25. doi: 10.1080/14767058.2018.1500540

PubMed Abstract | Crossref Full Text | Google Scholar

34. Avci G, Ulutas E, Ozdemir V, Kivrak I, and Bulbul A. The positive effect of black seed (Nigella sativa L.) essential oil on thyroid hormones in rats with hypothyroidism and hyperthyroidism. J Food Biochem. (2022) 46:e13801. doi: 10.1111/JFBC.13801

PubMed Abstract | Crossref Full Text | Google Scholar

35. Wani JA, Tsagkaris C, Majid S, Ganie MA, Akhter R, Ahmad SB, et al. Therapeutic effects of Nigella sativa on hormonal dysfunctions. In: Black seeds (Nigella sativa): pharmacological and therapeutic applications. Elsevier (2022). p. 217–38. doi: 10.1016/B978-0-12-824462-3.00001-9

Crossref Full Text | Google Scholar

36. Hassan HY and Fadel HK. Comparison between the efficacy of nigella sativa aqueous extract and its oil on methimazole-induced hypothyroidism in albino mice. Jordan J Pharm Sci. (2023) 16:103. doi: 10.35516/JJPS.V16I1.1072

Crossref Full Text | Google Scholar

37. Al Zarzour RH, Kamarulzaman EE, Saqallah FG, Zakaria F, Asif M, and Abdul Razak KN. Medicinal plants’ proposed nanocomposites for the management of endocrine disorders. Heliyon. (2022) 8:e10665. doi: 10.1016/j.heliyon.2022.e10665

PubMed Abstract | Crossref Full Text | Google Scholar

38. Taïbi K, Ait Abderrahim L, Helal F, and Hadji K. Ethnopharmacological study of herbal remedies used for the management of thyroid disorders in Algeria. Saudi Pharm J. (2021) 29:43–52. doi: 10.1016/j.jsps.2020.12.004

PubMed Abstract | Crossref Full Text | Google Scholar

39. Liu Y, Shi X, Tian Y, Zhai S, Liu Y, Xiong Z, et al. An insight into novel therapeutic potentials of taxifolin. Front Pharmacol. (2023) 14:1173855. doi: 10.3389/fphar.2023.1173855

PubMed Abstract | Crossref Full Text | Google Scholar

40. Proschak E, Stark H, and Merk D. Polypharmacology by design: A medicinal chemist’s perspective on multitargeting compounds. J Med Chem. (2019) 62:420–44. doi: 10.1021/acs.jmedchem.8b00760

PubMed Abstract | Crossref Full Text | Google Scholar

41. El-Akad RH, Elosaily AH, Gamil NM, Merghany RM, El-Shiekh RA, Almagharbeh WT, et al. Natural product-based interventions for thyroid disorders: mechanisms and applications. Naunyn Schmiedebergs Arch Pharmacol. (2025). doi: 10.1007/s00210-025-04604-y

PubMed Abstract | Crossref Full Text | Google Scholar

42. Mohanraj K, Karthikeyan BS, Vivek-Ananth RP, Chand RPB, Aparna SR, Mangalapandi P, et al. IMPPAT: A curated database of Indian Medicinal Plants, Phytochemistry and Therapeutics. Sci Rep. (2018) 8:1–17. doi: 10.1038/s41598-018-22631-z

PubMed Abstract | Crossref Full Text | Google Scholar

43. Vivek-Ananth RP, Mohanraj K, Sahoo AK, and Samal A. IMPPAT 2.0: an enhanced and expanded phytochemical atlas of Indian medicinal plants. ACS Omega. (2023) 8:8827–45. doi: 10.1021/acsomega.3c00156

PubMed Abstract | Crossref Full Text | Google Scholar

44. Agyapong O, Asiedu SO, Kwofie SK, Miller WA, Parry CS, Sowah RA, et al. Molecular modelling and de novo fragment-based design of potential inhibitors of beta-tubulin gene of Necator americanus from natural products. Inform Med Unlocked. (2021) 26:100734. doi: 10.1016/j.imu.2021.100734

PubMed Abstract | Crossref Full Text | Google Scholar

45. Ye L, Li Y-L, Mellström K, Mellin C, Bladh L-G, Koehler K, et al. Thyroid receptor ligands. 1. Agonist Ligands Selective for the Thyroid Receptor β 1. J Med Chem. (2003) 46:1580–8. doi: 10.1021/jm021080f

PubMed Abstract | Crossref Full Text | Google Scholar

46. Yang F, Zhang H, Meng X, Li Y, Zhou Y, Ling S, et al. Structural insights into thyrotropin-releasing hormone receptor activation by an endogenous peptide agonist or its orally administered analogue. Cell Res. (2022) 32:858–61. doi: 10.1038/s41422-022-00646-6

PubMed Abstract | Crossref Full Text | Google Scholar

47. Trott O and Olson AJ. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. (2010) 31:455–61. doi: 10.1002/jcc.21334

PubMed Abstract | Crossref Full Text | Google Scholar

48. Seeliger D and de Groot BL. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J Comput Aided Mol Des. (2010) 24:417–22. doi: 10.1007/s10822-010-9352-6

PubMed Abstract | Crossref Full Text | Google Scholar

49. Salentin S, Schreiber S, Haupt VJ, Adasme MF, and Schroeder M. PLIP: Fully automated protein-ligand interaction profiler. Nucleic Acids Res. (2015) 43:W443–7. doi: 10.1093/NAR/GKV315

PubMed Abstract | Crossref Full Text | Google Scholar

50. Fu L, Shi S, Yi J, Wang N, He Y, Wu Z, et al. ADMETlab 3.0: an updated comprehensive online ADMET prediction platform enhanced with broader coverage, improved performance, API functionality and decision support. Nucleic Acids Res. (2024) 52:W422–31. doi: 10.1093/nar/gkae236

PubMed Abstract | Crossref Full Text | Google Scholar

51. Agbebi EA, Adeyemi SO, Adewale AI, Ajofoyinbo OS, Olugbogi EA, Oyinloye OM, et al. Uvarinol and Dichamanetin Derived from Uvaria chamae as Potential Dual-Site Inhibitors Against PBP2a in Methicillin Resistant Staphylococcus aureus : An In Silico Study. Pharmaceuticals. (2025) 18(4):529.

PubMed Abstract | Google Scholar

52. Bowers KJ, Chow DE, Xu H, Dror RO, Eastwood MP, Gregersen BA, et al. (2006). Scalable algorithms for molecular dynamics simulations on commodity clusters, in: ACM/IEEE SC 2006 Conference (SC’06), pp. 43–3. Tampa, Florida, USA: IEEE. doi: 10.1109/SC.2006.54

Crossref Full Text | Google Scholar

53. Roos K, Wu C, Damm W, Reboul M, Stevenson JM, Lu C, et al. OPLS3e: extending force field coverage for drug-like small molecules. J Chem Theory Comput. (2019) 15:1863–74. doi: 10.1021/acs.jctc.8b01026

PubMed Abstract | Crossref Full Text | Google Scholar

54. Mitchell KE. Origin update. Sci (1979). (2000) 288:1982–2. doi: 10.1126/science.288.5473.1982a

Crossref Full Text | Google Scholar

55. Kollman PA, Massova I, Reyes C, Kuhn B, Huo S, Chong L, et al. Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc Chem Res. (2000) 33:889–97. doi: 10.1021/AR000033J

PubMed Abstract | Crossref Full Text | Google Scholar

56. Hou T, Wang J, Li Y, and Wang W. Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J Chem Inf Model. (2010) 51:69–82. doi: 10.1021/CI100275A

PubMed Abstract | Crossref Full Text | Google Scholar

57. Genheden S and Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. (2015) 10:449. doi: 10.1517/17460441.2015.1032936

PubMed Abstract | Crossref Full Text | Google Scholar

58. Wang E, Sun H, Wang J, Wang Z, Liu H, Zhang JZH, et al. End-point binding free energy calculation with MM/PBSA and MM/GBSA: strategies and applications in drug design. Chem Rev. (2019) 119:9478–508. doi: 10.1021/ACS.CHEMREV.9B00055

PubMed Abstract | Crossref Full Text | Google Scholar

59. dos Santos Nascimento IJ, Santana Gomes JN, de Oliveira Viana J, de Medeiros e Silva YMS, Barbosa EG, and de Moura RO. The power of molecular dynamics simulations and their applications to discover cysteine protease inhibitors. Mini-Reviews Medicinal Chem. (2023) 24:1125–46. doi: 10.2174/1389557523666230901152257/CITE/REFWORKS

PubMed Abstract | Crossref Full Text | Google Scholar

60. Shoaib TH, Almogaddam MA, Andijani YS, Saib SA, Almaghrabi NM, Elyas AF, et al. Marine-derived compounds for CDK5 inhibition in cancer: integrating multi-stage virtual screening, MM/GBSA analysis and molecular dynamics investigations. Metabolites. (2023) 13:1090. doi: 10.3390/METABO13101090

PubMed Abstract | Crossref Full Text | Google Scholar

61. Onikanni SA, Lawal B, Munyembaraga V, Bakare OS, Taher M, Khotib J, et al. Profiling the Antidiabetic Potential of Compounds Identified from Fractionated Extracts of Entada africana toward Glucokinase Stimulation : Computational Insight. Molecules (Basel, Switzerland) (2023) 28(15):5752. doi: 10.3390/molecules28155752

PubMed Abstract | Crossref Full Text | Google Scholar

62. McAninch EA and Bianco AC. The history and future of treatment of hypothyroidism. Ann Intern Med. (2016) 164:50–6. doi: 10.7326/M15-1799

PubMed Abstract | Crossref Full Text | Google Scholar

63. Wang F, Yang W, Shi Y, and Le G. Structural analysis of selective agonists of thyroid hormone receptor β using 3D-QSAR and molecular docking. J Taiwan Inst Chem Eng. (2015) 49:1–18. doi: 10.1016/j.jtice.2014.11.009

Crossref Full Text | Google Scholar

64. Sahu VK, Khan AKR, Singh RK, and Singh PP. Hydrophobic, polar and hydrogen bonding based drug-receptor interaction of tetrahydroimidazobenzodiazepinones. American Journal of Immunology. (2008) 4(3):33–42. doi: 10.3844/ajisp.2008.33.42

Crossref Full Text | Google Scholar

65. Han SY, Qiao JQ, Zhang YY, Lian HZ, and Ge X. Determination of n-octanol/water partition coefficients of weak ionizable solutes by RP-HPLC with neutral model compounds. Talanta. (2012) 97:355–61. doi: 10.1016/J.TALANTA.2012.04.045

PubMed Abstract | Crossref Full Text | Google Scholar

66. Korinth G, Wellner T, Schaller KH, and Drexler H. Potential of the octanol–water partition coefficient (logP) to predict the dermal penetration behaviour of amphiphilic compounds in aqueous solutions. Toxicol Lett. (2012) 215:49–53. doi: 10.1016/j.toxlet.2012.09.013

PubMed Abstract | Crossref Full Text | Google Scholar

67. Kralj S, Jukič M, and Bren U. Molecular filters in medicinal chemistry. Encyclopedia. (2023) 3:501–11. doi: 10.3390/encyclopedia3020035

Crossref Full Text | Google Scholar

68. Ivanović V, Rančić M, Arsić B, and Pavlović A. Lipinski’s rule of five, famous extensions and famous exceptions. Chemia Naissensis. (2020) 3:171–81. doi: 10.46793/CHEMN3.1.171I

Crossref Full Text | Google Scholar

69. Grewal S, Deswal G, Grewal AS, and Guarve K. Molecular dynamics simulations: Insights into protein and protein ligand interactions. Advances in pharmacology (San Diego, Calif.) (2025) 103:139–62. doi: 10.1016/bs.apha.2025.01.007.

PubMed Abstract | Crossref Full Text | Google Scholar

70. Rudrapal M, Issahaku AR, Agoni C, Bendale AR, Nagar A, Soliman MES, et al. In silico screening of phytopolyphenolics for the identification of bioactive compounds as novel protease inhibitors effective against SARS-CoV-2. J Biomol Struct Dyn. (2022) 40:10437–53. doi: 10.1080/07391102.2021.1944909

PubMed Abstract | Crossref Full Text | Google Scholar