Identification of Effective and Nonpromiscuous Antidiabetic Drug Molecules from <i>Penicillium</i> Species

Saleem, Shahzad; Bibi, Shabana; Yousafi, Qudsia; Hassan, Tehzeem; Khan, Muhammad Saad; Hasan, Mohammad Mehedi; Chopra, Hitesh; Moustafa, Mahmoud; Al-Shehri, Mohammed; Khalid, Mohammad; Kabra, Atul

doi:https://doi.org/10.1155/2022/7040547

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Plant-Derived Bioactive Compounds as an Antidiabetic Agent

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 7040547 | https://doi.org/10.1155/2022/7040547

Identification of Effective and Nonpromiscuous Antidiabetic Drug Molecules from Penicillium Species

Shahzad Saleem,¹Shabana Bibi,^2,3Qudsia Yousafi ,¹Tehzeem Hassan,¹Muhammad Saad Khan,¹Mohammad Mehedi Hasan,⁴Hitesh Chopra,⁵Mahmoud Moustafa,^6,7Mohammed Al-Shehri,⁶Mohammad Khalid,⁸and Atul Kabra⁹

Academic Editor: Fernanda Tonelli

Received15 Apr 2022

Accepted08 May 2022

Published08 Jun 2022

Abstract

Diabetes mellitus (DM) is a very common metabolic disorder/disease. The deterioration of β-cells by autoimmune system is the hallmark of this disease. Thioredoxin-Interacting Protein (TXNIP) is responsible for β-cells degradation by T-cells in the pancreas. This protein had been declared a good drug target for controlling DM. Lots of side effects have been reported as a result of long-time consumption of conventional antidiabetic drugs. The development of new and effective drugs with the minimal side effects needs time. TXNIP was selected as a target for Computer-Aided Drug Design. The antidiabetic fungal metabolite compounds were selected from the literature. The compounds were screened for their drug-likeness properties by DruLiTo and DataWarior tools. Twenty-two drug-like fungal compounds were subjected to Quantitative Structure-Activity Relationship (QSAR) analysis by using CheS-Mapper 2.0. The lowest (0.01) activity cliff was found for three compounds: Pinazaphilone A, Pinazaphilone B, and Chermesinone A. The highest value for apol (81.76) was shown by Asperphenamate, while Albonoursin and Sterenin L showed highest score (40.66) for bpol. The lowest value (0.46) for fractional molecular frame (FMF) was calculated for Pinazaphilone A and Pinazaphilone B. TPSA for Pinazaphilone A and Pinazaphilone B was 130.51 Å². was observed for all the twenty-two compounds. Molecular docking of fungal compounds with TXNIP was done by AutoDock Vina. The binding energy for complexes ranged between −9.2 and −4.6 kcal/mol. Four complexes, TXNIP-Pinazaphilone A, TXNIP-Pinazaphilone B, TXNIP-Asperphenamate, and TXNIP-Sterenin L, were selected for MD simulation to find out the best lead molecule. Only one complex, TXNIP-Pinazaphilone B, showed a stable conformation throughout the 80 ns run of MD simulation. Pinazaphilone B derived from the Penicillium species fungi was selected as the lead molecule for development of antidiabetic drug having the least side effects.

1. Introduction

Diabetes mellitus (DM) is a very common disease related to the metabolic disorders. It is prevalent around the world irrespective of the economic conditions [1]. Type 1 DM is more common among the children from infancy to 14 years of age [2]. Type II DM was found to be more prevalent in adults but it has also been reported in children recently [2]. Diabetes had affected 8.5% of persons aged 18 years and above in 2014. Diabetes had been reported as the direct cause of 1.5 million deaths in 2019, with 48% of deaths occurring before the age of 70 years. The premature mortality due to diabetes was found to be more than 5% during the period from 2010 to 2016, with the highest rate being in poor countries. It has been estimated that more than 629 million people will be affected by this disease by 2049 [3]. Type I DM is induced by deterioration of pancreatic β-cell by T-cells of human immunity system [4]. Type II DM is related to insulin resistance and is more common in obese people.

The hallmark of the type 1 and type 2 diabetes is the destruction of pancreatic β-cells [5, 6]. The treatment strategies for this disease are focused on enhancing β-cells functioning by reducing their apoptosis. Thioredoxin-Interacting Protein (TXNIP) has been identified as an important glucose-induced gene [7]. The role of TXNIP as a β-cell destructor has been revealed in recent decades [8]. The overexpression of this protein had been reported in diabetic patients [9]. TXNIP is a promising target that has just emerged in this field of study [10]. A good number of antidiabetic drugs have been developed to reduce the TXNIP level in blood [11]. TXNIP inhibition has been shown to have beneficial effects beyond the β-cell including various extra pancreatic tissues [12], such as the cardiovascular system [9], kidney [13], and retina [14]. TXNIP has been identified as a crucial regulator of pancreatic β-cell malfunction and death, both of which are important factors in the etiology of type 1 and type 2 diabetes. TXNIP appears to be a promising therapeutic target for diabetes, according to growing evidence based on basic, preclinical, and retrospective epidemiological research [15].

Lifelong diabetes medications have some side effects. The drugs might have side effects due to their off-target activities. Angioedema is a disorder that results in larynx, uvula, and tongue swelling, which leads to breathing problems [16]. This is due to off-target activity of antidiabetic drugs against Angiotensin Converting Enzyme (ACE) [17]. Most of the drugs are associated with skin problems [18]. Scars formation is a very common side effect of antidiabetic drugs [19]. Some of the side effects are allergy-related, especially to insulin intake [20]. Sulfonylurea drugs were found to be involved in hyperexcitation of β-cells resulting in their apoptosis [21]. This side effect leads to fatigue and exhaustion in patients. About one-third of cases of allergy are due to insulin itself [22] and the rest are related to the components of insulin formulation, that is, zinc protamine and meta-cresol [23]. In addition to side effects, treatment failure in some was reported in some patients using very promising drugs [24]. Some important therapies with their side effects are as follows. First is biguanides: Metformin is usually the first prescribed medicine by the doctors to treat type II diabetes. This improves body’s way to use insulin and lowers blood sugar. It may cause serious lactic acidosis after continues use over long time. Second is sulfonylureas: Glipizide is another important drug which controls pancreas to produce more insulin. The patients feel shaky, dizzy, sweaty, and confused due to drastically lowering blood sugar level. Third is meglitinides: nateglinide also controls pancreas and more insulin is produced. This drug is not retained in the body and is excreted rapidly, but it still has side effects like lowering blood sugar and weight gain. Fourth is thiazolidinediones (TZDs) : Pioglitazone works by boosting the insulin use by body. The water retention in body increases, which causes swelling. Weight gain and increase of bad cholesterol are its side effects. Some life-threating effects are bone fracture heart failure and bladder cancer. Fifth is alpha-glucosidase inhibitors: Acarbose is a drug taken with first loaf of meal which slows down the breakdown of carbohydrates. This affects the metabolism of carbohydrates and digestion; some common problems faced by using this medicine are gases, stomach pain, and diarrhea. Sixth is DPP-4 inhibitors: alogliptin helps the pancreas to make more insulin. But associated problems are stuffy nose, sore throat, and stomach pain. Seventh is insulin therapy, which is the last resort for diabetes treatment. The insulin is injected in the body with the help of syringe or pen injector. This treatment is not free of side effects. Some common insulin therapy issues are cough, mouth dryness, dizziness, headache, and rashes.

The fungi have a very unique metabolism system and produce metabolites with diversified nature [25]. The discovery of penicillin by Alexander Fleming in 1929 has proven a to be breakthrough in the fungal metabolites research [25]. Besides, penicillin (antibacterial), echinocandin B (antifungal), cyclosporine A (immunosuppressive), and lovastatin (cholesterol-lowering) are all fungal originated and marketed pharmaceutical drugs. Furthermore, fungal metabolites have been reported to be used in human diet and health improvement [26].

New drug development is a time-consuming and expensive approach, which might be the reason for scarce research new drug development [27]. The approval of new drug molecules from FDA takes more than fourteen years [28] with an estimated cost of production of more than 800 million dollars [29]. Most of the drugs fail in clinical trials, which waste entire efforts and expenses of initial screening [30]. It is important to curtail the cost of production of new drugs and expedite this process at successful and smooth pace. More than half of these problems can be solved by using Computer Aided Drug Design (CADD) methods [31]. A large number of molecules from different databases can be screened out in terms of their toxicity and drug-likeness properties by using CADD techniques [32]. The side effects cause a serious problem in the lifetime of using life-saving drugs. This problem can also be addressed by computational Quantitative Structure-Activity Relationship (QSAR) analysis of selected drug molecules [33]. These CADD approaches can be used in combination with wet lab methods for rapid development of new and effective drugs with lesser side effects [34].

Miglitol and Acarbose are two commonly used commercial antidiabetic drugs [35]. These drugs were found to be associated with the serious liver injuries and gastrointestinal damage in patients using them regularly [36]. Identification and development of new and safer antidiabetic medications with low toxicity and few side effects are of crucial importance. Secondary metabolites from plants have been extensively studied for their possible antidiabetic activities. Different chemicals produced by fungi might be used as active ingredients in new antidiabetic drug development. In this study, we used CADD approaches for designing antidiabetic drugs, targeting TXNIP, from fungal metabolites. The QSAR analysis was done for predicting the least promiscuous drug molecules to avoid side effects.

2. Materials and Methods

2.1. Retrieval and Preparation of Ligands

Different fungal metabolites with antidiabetic properties were selected from literature survey. Three-dimensional (3D) structures of selected compounds were downloaded from PubChem and ChEBI in SDF format. The ligand molecules were prepared, for docking, by MGL tools in AutoDock Vina. The roots were detected and torsion count was done. The prepared molecules were saved in PDBQT format.

2.2. Retrieval and Preparation of Protein

The Thioredoxin-Interacting Protein (TXNIP) was selected as target protein. The three-dimensional structure of the protein’s N-terminal domain was retrieved in PDB format from Protein Data Bank (PDB ID : 4GEI) (Figure 1). The protein was purified by removing water molecules and het-atoms through UCSF Chimera. Polar hydrogen and charges, Gasteiger charge (−9.9935) and Kollman’s charge (−211.153), were added by MGL tools of AutoDock Vina. The 3D grid box, as ligand molecule binding site, was adjusted as 39.708 × 33.059 × 4.623 Å. The prepared protein structure was saved in PDBQT format.

2.3. Drug-Likeness and Toxicity Analysis of Ligands

The compounds were screened for their drug-likeness by DruLiTo and nontoxicity by DataWarior tools. Mutagenic, tumorigenic, nasty function, irritant, carcinogenic, and reproductive effects were predicted by DataWarrior and the compounds with none of these properties were selected. The Drug-likness parameters like number of hydrogen bond donors/acceptors, , and molecular weight were predicted by DruLiTo. The compounds following Lipinski rule of five were selected.

2.4. Binding Affinity and Molecular Docking

Molecular docking of N-terminal domain of TXNIP was done with selected (nontoxic and drug-like) fungal compounds (ligands). Protein-ligand blind docking was done using AutoDock Vina tool. The PDBQT files of protein and ligands were used in this process. The results were generated in the form of out.pdbqt file. The out.pdbqt and protein complexes were generated by using PyMOL software. The complex files were saved in PDB format and visualized by LIGPLOT (2D) and BIOVIA Discovery Studio (3D). The protein-ligand complexes showing lowest binding affinities and good hydrogen bond interactions were selected for further studies.

3. Quantitative Structure-Activity Relationship (QSAR) Analysis of Selected Compounds

The compounds that qualified the drug-likness and nontoxicity parameters were further analyzed for their Quantitative Structure-Activity Relationship (QSAR) by using CheS-Mapper 2.0 [36]. Statistical approaches were used to compare the predicted and experimental results of compounds. Different QSAR properties were selected for analyses and groups of compounds were created according to their structural features and physical properties.

4. Molecular Dynamics (MD) Simulation

Molecular Dynamics (MD) simulation study is a useful strategy towards the understanding of protein-ligand interactions in detail at the atomic level [37]. The top four docked complexes, in terms of lower binding energies, that is, TXNIP-Pinazaphilone A (CID_122182011; Penicillium species), TXNIP-Pinazaphilone B (CID_122182012; Penicillium species), TXNIP-Asperphenamate (CID_173952; Penicillium spathulatum), and TXNIP-Sterenin L (CID_77461068; Stereum hirsutum), were selected for MD simulation to find out the best lead molecule. Molecular Dynamics (MD) simulations were performed, for 80 ns time interval, using the Amber v18 software program. For the addition of hydrogen atoms to the protein-ligand combination, the LeaP module was used. To maintain the system’s neutrality, counterions (Na⁺ and Cl⁻) were also introduced. For system solvation, a truncated octahedral box of the TIP3P water model with 15 buffers was utilized. For the treatment of long-range electrostatic interactions, the PME (Particle Mesh Ewald) method was used. For simulation, the ff14SB force field was used [38]. Bonds involving hydrogen atoms were constrained using the SHAKE method [39]. To speed up simulations, the PMEMD of the CUDA version of the PMEMD was used [40]. In the NVT ensemble, the steepest descent approach was used to minimize the solvated system for roughly 20000 steps, followed by heating for up to 400 ps and then equilibration for 200 ps. The manufacturing runs at 298 kelvin with a 2 ps time step in Berendsen's thermostat and barostat under the NPT ensemble. The trajectories were analyzed using the CPPTRAJ program.

5. Results

5.1. Binding Affinity and Interactions in Molecular Docking

The overall binding affinity among docked complexes ranged from −9.2 to −4.6 kcal/mol (Table 1). The lowest binding affinity (−9.2) was observed for one compound (CID_10068406), but no hydrogen bonds were found, so no interaction of ligands was observed with the protein. Only fifteen ligand molecules showed hydrogen bonding with the Thioredoxin-Interacting Protein (TXNIP). Four ligand-receptor complexes formed one hydrogen bond, six ligand-receptor complexes formed two hydrogen bonds, and three complexes formed three hydrogen bond interactions, whereas only two ligand-receptor complexes formed four hydrogen bonds.

5.2. Quantitative Structure-Activity Relationship (QSAR) Analysis of Selected Compounds

Only twenty-two nontoxic and drug-like fungal compounds were selected for QSAR analysis (Table 2). The compounds were clustered in two groups (Table 3). In this study, we selected seven features (Table 4) for QSAR analysis of the selected compounds. Activity cliffs (ACs) are pairs or groups of structurally identical compounds which are functional against the same target but have significant potency variances. The two groups of compounds showed more subgroups for activity cliff. Two compounds in cluster I, CID_54436672 (2-methoxy-4,5-dihydroxybenzaldehyde) and CID_528089 (abscisic acid), showed similar value of 0.02. Three compounds from cluster II, CID_122182012 (Pinazaphilone B), CID_122182011 (Pinazaphilone A), and CID_53355009 (Chermesinone A), showed similar value of 0.01. The sum of the atomic polarizabilities (including implicit hydrogens) is calculated as apol [41] and the sum of the absolute values of the difference between atomic polarizabilities of all bonded atoms in the molecule (including implicit hydrogens) is calculated as bpol. The highest value of 81.73 for apol was recorded for CID_173952 (Asperphenamate). For cluster II, two structurally close compounds, CID_77461068 (Sterenin L) and CID_77461069 (Sterenin M), showed the highest value of 75.71, while their other relative CID_77461067 (Sterenin K) showed slightly lower value, 69.52. The highest bpol in cluster I (38.91) was in the case of CID_173952 (Asperphenamate) followed by 29.73 for ChEBI: 171609 (Albonoursin). In cluster II, an overlapping cluster of 3 structurally related compounds showed the highest values of bpol. Two compounds, CID_77461069 (Albonoursin) and CID_77461068 (Sterenin L), shared the same value, that is, 40.66, and CID_77461067 (Sterenin K) showed a value of 36.29. A cluster of two other compounds, CID_122182011 (Pinazaphilone A) and CID_122182012 (Pinazaphilone B), share similar values (28.57) for bpol. Fractional molecular frame (FMF), the key descriptor utilized in the analysis, is based on Bemis and Murcko's hierarchical molecular classification method [42, 43]. Fractional molecular frame (FMF) was calculated for the compounds in both groups. A similar trend to other descriptors was observed. In cluster I, four compounds, CID_173952 (Asperphenamate), CID_10068406 (Benzomalvin A), CID_64435512 (Benzomalvin B), and CID_7090 (Phenidone), shared the same value, that is, 0.92. Similarly, in cluster II, two structurally related compounds shared the same value of 0.47, while two other compounds, CID_12218211 (Pinazaphilone A) and CID_12218212 (Pinazaphilone B), showed 0.46 FMF value. Total polar surface area (TPSA) showed a trend similar to that in other QSAR descriptors. The three structurally related compounds CID_77461067, CID_77461068, and CID_77461069 showed a similar value of 128.16 Å². Two other similar compounds, CID_122182011 and CID_122182012, showed the same value of 130.51 Å². The value of was observed for all the compounds. The number of hydrogen bond acceptors <10 and the number of hydrogen bond donors <5 were recorded. These conditions fulfill the requirements for Lipinski rule of 5.

5.3. Molecular Dynamics Simulation

The top four selected docked complexes, that is, TXNIP-Pinazaphilone A (CID_122182011; Penicillium species), TXNIP-Pinazaphilone B (CID_122182012; Penicillium species), TXNIP-Asperphenamate (CID_173952; Penicillium spathulatum), and TXNIP-Sterenin L (CID_77461068; Stereum hirsutum), were simulated in an explicit water environment for a total 80 ns time period. The complex TXNIP-Pinazaphilone B showed stable behavior throughout the run with four ligand-receptor hydrogen bonds at each time slot (Figures 2(a)–2(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

5.4. Root Mean Square Deviation (RMSD)

The simulation process stability was tested measuring the backbone atoms deviation using Root Mean Square Deviation (RMSD). In general, the RMSD result revealed that the curve showed a progressive increase starting at 1 and oscillating around 2.2 to 3.2 (Figure 3). A peak of 4.3 Å was observed at 49 ns and after that there was a decline, which favors the stability and reliability of the complex.

5.5. Radius of Gyration (Rg)

We investigated the complex’s overall compactness further using MD simulation and radius of gyration (Rg). The result revealed a very consistent behavior in terms of Rg value below 19.25 Å throughout the MD simulation time (Figure 4).

5.6. Root Mean Square Fluctuations (RMSF)

We looked into Root Mean Square Fluctuations (RMSF) to better understand the fluctuation of individual residues in the docked complex, which gives us more information about each residue’s flexibility. The results of this analysis revealed fluctuation peaks for glutamine 53, glutamic acid 74, tyrosine 113, and valine 140 at 4.5 Å, 3.8 Å, 2.7 Å, and 4.1 A, respectively. The high fluctuation in the docked complex residues might be due to the free movement of the residues, while most of the residues in docked protein showed a steady behavior which might indicate the interaction with the ligand and that is why that region showed a stable behavior with less RMSF curve (Figure 5). The RMSF analysis indicates that the binding of ligand to protein showed a very consistent association.

6. Discussion

Diabetes mellitus is a metabolic syndrome which affects millions of people around the world. It is a metabolic disorder of glucose level dysregulation. Different types of drugs are available in the market, among which Miglitol and Acarbose are more commonly prescribed. These drugs hinder the carbohydrate absorption from the gut and can be taken alone or in combination with insulin [44]. As these medicines have to be used for the rest of patients’ life, they pose serious side effects in the patients [45]. These associated problems tend to restrict the prescription and uses of certain drugs [36]. This threatening situation of use of conventional antidiabetic drugs leads to an urge for development of some new safer and least promiscuous drugs.

This study was conducted to predict some nontoxic and effective lead molecules, from fungal metabolites, for inhibition of Thioredoxin-Interacting Protein (TXNIP) with minimal side effects. The overexpression of TXNIP causes death of β-cells, which is a major cause of diabetes mellitus [46]. Inhibition of this protein had been recommended for targeted antidiabetic drug development. Fungal metabolites have been reported to have many medicinal properties, especially antidiabetic and anticancer nature [47]. These fungal metabolites can be used for synthesizing antidiabetic drugs [48]. The structures of fungal metabolite compounds downloaded from PubChem and ChEBI were screened for their drug-likeness and nontoxic nature. The pharmacokinetic properties, that is, hydrogen bod donors, hydrogen bond acceptors, atomic molar refractivity, number of rotatable bonds, and [49], were determined by using Drug Likeness Tool (DruLiTo) 2.1. Lipinski postulated five rules for a molecule to become a drug. He declared that a drug has poor absorption if it has more than 5 hydrogen bond donors, 10 hydrogen bond acceptors, greater than 5, and molecular weight above 500 Da. The Lipinski rule of five is popularly followed to define the druggability of the chemical compounds [50]. Twenty-six out of thirty-three compounds passed Lipinski rule of five filter. These compounds were further screened out for their nontoxic nature by submitting to DataWarior tool.

The binding affinities of selected fungal metabolite molecules against TXNIP were investigated using Molecular Docking simulation in this work. Molecular docking is an integral component of Computer Aided Drug Design (CADD) pipeline [51]. This approach can be used to generate protein-ligand interaction complexes with high precision [52] and also for ligand target site identification [41]. The AutoDock Vina tool was used to generate best docked models of protein (TXNIP) and ligands (fungal metabolite molecules) complex. This tool uses scoring functions, in terms of binding affinity, to predict protein-ligand complex interaction strength. Ten different poses, with respective binding energies and RMSD, were generated, which are ranked by scoring functions [53]. The lower the binding energy is, the higher the protein-ligand bonding interaction is. Twenty-two compounds showed good binding affinity between −9.2 and −5.2 kcal/mol.

The selected 22 antidiabetic fungal compounds were also submitted to CheS-Mapper for Quantitative Structure-Activity Relationship (QSAR) Analysis. The pharmacokinetic study of chemical compounds done through QSAR is a very important study in cheminformatics [54]. This approach is used in medicinal chemistry and CADD for lead optimization. In this study, we have studied six descriptors, activity cliffs (ACs), apol, bpol, total polar surface area (TPSA), , and fractional molecular frame (FMF), of the selected compounds. Activity cliff of the set of compounds shows their structural closeness but potency might be different [55]. The polarizability is defined as a linear response of electronic charge distribution with respect to an external applied electric field [56]. The study’s electronic interaction of chemical compounds is very crucial in the field of medicinal chemistry for proper chemical-protein interaction [57]. Two descriptors of polarizability, apol and bpol, have been calculated for selected compounds in this study. As the antidiabetic drugs are to be taken for the rest of patients’ life, the side effects of these medicines are a very serious concern. Different skin disease, heart diseases, and kidney dysfunctions have been reported as antidiabetic-drug-related issues. Promiscuity of a drug can be defined as its ability to interact with multiple proteins and affect their functioning, which results in unrequired pharmacological effects [58]. The promiscuous drugs, if used for a long time, can result in undesired drug metabolism pharmacokinetics [59]. This is an emerging problem in the field of medicinal chemistry and special emphasis has been raised for development of specific and least promiscuous drugs [60]. Three QSAR descriptors, that is, molecular weight, , and FMF, are used to describe the promiscuity of a drug molecule [61]. Molecular weight of a drug molecule is very critical for its absorption through the membrane. It is also reported that the promiscuity of the drug molecule is directly proportional to its promiscuity [62]. The molecular weight of the selected lead molecule was 416.40 g/mol, which shows that it has better absorption and has low promiscuity tendency. The quantification of the molecular lipophilicity is generally done by value. It is a very critical characteristic of the molecule which should be considered at the time of designing a drug, as higher results in higher promiscuity and least transport across biological membrane [57]. The compounds having are more promiscuous, while the least promiscuity was reported in the case of [61]. The selected compound for this study showed , while the lead molecule showed , which ensures its least promiscuous nature. The intestinal absorption and blood brain barrier (BBB) penetration of the drug molecule is measured by total polar surface area (TPSA) descriptor [63]. Total polar surface area more than 140 Å² renders least permeability across biological membrane [64]. A molecule with TPSA < 140 Å² has good biological membrane permeability and bioavailability. Most of the compounds in our dataset follow this rule and selected lead molecule (Pinazaphilone B) has TPSA value of 130.51 Å², which reflects its good permeability and bioavailability.

Top four compounds, Pinazaphilone A (CID_122182011), Pinazaphilone B (CID_122182012), Asperphenamate (CID_173952), and Sterenin L (CID_77461068), showing minimum binding energy and good hydrogen bonding with target protein were submitted for Molecular Dynamics (MD) simulation, for lead selection. Molecular simulation protocols are popularly used in CADD pipeline to study detailed structural dynamics at atomic level. The most important use of MD simulation is to test the accuracy of the predicted models in comparison with the experimental results [65]. The standard way to measure the stability of the docked complex is the RMSD of all heavy atoms with respect to the X-ray structure. The RMSD is an important tool that is used to characterize the conformational changes of ligand-protein complexes [66]. The RMSD for the TXNIP-Pinazaphilone B complex showed very little fluctuation in the start and then started to decline; this trend confirms the stability of the complex [67]. The radius of gyration (Rg) is a basic measurement of the overall size of a chain molecule [68]. Each residue’s RMSF is simple to calculate across the trajectory, but they can also be translated to temperature factors, which are present for each atom in a PDB file. The Pinazaphilone B-TXNIP complex showed a stable behavior in terms of different MD simulation descriptors.

Pinazaphilone B (CID_122182012) has been identified as an effective inhibitor of TXNIP. This can be qualified as lead molecule for development of effective and least promiscuous antidiabetic drugs. Pinazaphilone B secondary metabolite was extracted from fungus Penicillium sp. (Phylum: Ascomycota). It is a well-known fungus in medical history with the discovery of its antibiotic potential [69]. Its adaptable nature increases its presence in diversified environment, that is, air, soil, water, food, and vegetation. This is an economically important genus with its use in medicinal industry and food production and preservation. Some species are environmentally important because of their role as decomposers [70]. Some species are pathogenic in nature and produce a wide range of mycotoxins [71].

7. Conclusion

Fungal metabolites can be used in the development of drugs for many diseases. The side effects and toxicity of the medicines can be reduced by using natural compounds like fungal metabolites as active ingredients of drug. The patients of chronic and lifelong diseases are more prone to the harmful side effects of the medicines. Diabetes is one of the most common diseases in which drugs’ side effects have been observed in patients, which sometimes lead to life-threating conditions. In this study, fungal metabolites are used to test their efficacy against TXNIP protein. The Computer Aided Drug Designing techniques have been used to develop targeted drug. Further QSAR analysis confirmed the nonpromiscuous nature of the lead molecule. Only one fungal compound, Pinazaphilone B, showed a stable interaction with TXNIP. Pinazaphilone B derived from Penicillium fungi was selected as lead molecule in the development of antidiabetic drug having the least side effects. This fungal metabolite can be tested in vitro and in vivo to confirm efficacy for inhibition of TXNIP. This compound might be used solely or in combination with some commercially available antidiabetic compounds for development of safe, effective, and nonpromiscuous drugs.

Data Availability

Data will be available upon request to the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

All authors have contributed to the paper and read and approved the final manuscript.

Acknowledgments

This research was supported by the Deanship of Scientific Research at King Khalid University, Saudi Arabia, under Grant no. RGP. 2/90/43.

References

N. G. Forouhi, A. Koulman, S. J. Sharp et al., “Differences in the prospective association between individual plasma phospholipid saturated fatty acids and incident type 2 diabetes: the EPIC-InterAct case-cohort study,” Lancet Diabetes and Endocrinology, vol. 2, no. 10, pp. 810–818, 2014.
View at: Publisher Site | Google Scholar
K. Sahoo, B. Sahoo, A. Choudhury, N. Sofi, R. Kumar, and A. Bhadoria, “Childhood obesity: causes and consequences,” Journal of Family Medicine and Primary Care, vol. 4, no. 2, p. 187, 2015.
View at: Publisher Site | Google Scholar
E. English and E. Lenters-Westra, “HbA1c method performance: the great success story of global standardization,” Critical Reviews in Clinical Laboratory Sciences, vol. 55, no. 6, pp. 408–419, 2018.
View at: Publisher Site | Google Scholar
A. L. Burrack, T. Martinov, and B. T. Fife, “T cell-mediated beta cell destruction: autoimmunity and alloimmunity in the context of type 1 diabetes,” Frontiers in Endocrinology, vol. 8, 2017.
View at: Publisher Site | Google Scholar
J. Montane, L. Cadavez, L. Cadavez-Trigo, and A. Novials, “Stress and the inflammatory process: a major cause of pancreatic cell death in type 2 diabetes,” Diabetes, Metabolic Syndrome and Obesity: Targets and Therapy, vol. 7, p. 25, 2014.
View at: Publisher Site | Google Scholar
V. Poitout and R. P. Robertson, “Minireview: secondary β-cell failure in type 2 diabetes—a convergence of glucotoxicity and lipotoxicity,” Endocrinology, vol. 143, pp. 339–342, 2002.
View at: Publisher Site | Google Scholar
A. Shalev, C. A. Pise-Masison, M. Radonovich et al., “Oligonucleotide microarray analysis of intact human pancreatic islets: identification of glucose-responsive genes and a highly regulated TGFβ signaling pathway,” Endocrinology, vol. 143, no. 9, pp. 3695–3698, 2002.
View at: Publisher Site | Google Scholar
A. H. Minn, C. Hafele, and A. Shalev, “Thioredoxin-interacting protein is stimulated by glucose through a carbohydrate response element and induces β-cell apoptosis,” Endocrinology, vol. 146, no. 5, pp. 2397–2405, 2005.
View at: Publisher Site | Google Scholar
J. Chen, M. E. Young, J. C. Chatham, D. K. Crossman, L. J. Dell’Italia, and A. Shalev, “TXNIP regulates myocardial fatty acid oxidation via miR-33a signaling,” American Journal of Physiology—Heart and Circulatory Physiology, vol. 311, no. 1, pp. H64–H75, 2016.
View at: Publisher Site | Google Scholar
A. Shalev, “Minireview: thioredoxin-interacting protein: regulation and function in the pancreatic β-cell,” Molecular Endocrinology, vol. 28, no. 8, pp. 1211–1220, 2014.
View at: Publisher Site | Google Scholar
M. Shaked, M. Ketzinel-Gilad, Y. Ariav, E. Cerasi, N. Kaiser, and G. Leibowitz, “Insulin counteracts glucotoxic effects by suppressing thioredoxin- interacting protein production in INS-1E beta cells and in Psammomys obesus pancreatic islets,” Diabetologia, vol. 52, no. 4, pp. 636–644, 2009.
View at: Publisher Site | Google Scholar
W. A. Chutkow, P. Patwari, J. Yoshioka, and R. T. Lee, “Thioredoxin-interacting protein (Txnip) is a critical regulator of hepatic glucose production,” Journal of Biological Chemistry, vol. 283, no. 4, pp. 2397–2406, 2008.
View at: Publisher Site | Google Scholar
Y. Shi, Y. Ren, L. Zhao et al., “Knockdown of thioredoxin interacting protein attenuates high glucose-induced apoptosis and activation of ASK1 in mouse mesangial cells,” FEBS Letters, vol. 585, no. 12, pp. 1789–1795, 2011.
View at: Publisher Site | Google Scholar
M. F. El-Azab, B. R. B. Baldowski, B. A. Mysona et al., “Deletion of thioredoxin-interacting protein preserves retinal neuronal function by preventing inflammation and vascular injury,” British Journal of Pharmacology, vol. 171, no. 5, pp. 1299–1313, 2014.
View at: Publisher Site | Google Scholar
L. Thielen and A. Shalev, “Diabetes pathogenic mechanisms and potential new therapies based upon a novel target called TXNIP,” Current Opinion in Endocrinology Diabetes and Obesity, vol. 25, no. 2, pp. 75–80, 2018.
View at: Publisher Site | Google Scholar
M. A. Wu, F. Perego, A. Zanichelli, and M. Cicardi, “Angioedema phenotypes: disease expression and classification,” Clinical Reviews in Allergy and Immunology, vol. 51, no. 2, pp. 162–169, 2016.
View at: Publisher Site | Google Scholar
S. I. Scott, M. F. Andersen, L. Aagaard, C. V. Buchwald, and E. R. Rasmussen, “Dipeptidyl peptidase-4 inhibitor induced angioedema—an overlooked adverse drug reaction?” Current Diabetes Reviews, vol. 14, no. 4, pp. 327–333, 2018.
View at: Publisher Site | Google Scholar
P. M. Jedlowski, C. H. Te, R. J. Segal, and M. T. Fazel, “Cutaneous adverse effects of diabetes mellitus medications and medical devices: a review,” American Journal of Clinical Dermatology, vol. 20, no. 1, pp. 97–114, 2019.
View at: Publisher Site | Google Scholar
P. Ross, A. Gray, J. Milburn et al., “Insulin pump-associated adverse events are common, but not associated with glycemic control, socio-economic status, or pump/infusion set type,” Acta Diabetologica, vol. 53, no. 6, pp. 991–998, 2016.
View at: Publisher Site | Google Scholar
J. Jacquier, C. L. Chik, and P. A. Senior, “A practical, clinical approach to the assessment and management of suspected insulin allergy,” Diabetic Medicine, vol. 30, no. 8, pp. 977–985, 2013.
View at: Publisher Site | Google Scholar
A. Rosengren, X. Jing, L. Eliasson, and E. Renström, “Why treatment fails in type 2 diabetes,” PLoS Medicine, vol. 5, no. 10, p. e215, 2008.
View at: Publisher Site | Google Scholar
D. Kim and J. Baraniuk, “Delayed-type hypersensitivity reaction to the meta-cresol component of insulin,” Annals of Allergy, Asthma, and Immunology, vol. 99, no. 2, pp. 194-195, 2007.
View at: Publisher Site | Google Scholar
H. Gin and J. Aubertin, “Generalized allergy due to zinc and protamine in insulin preparation treated with insulin pump,” Diabetes Care, vol. 10, no. 6, pp. 789-790, 1987.
View at: Publisher Site | Google Scholar
B. Abubakar, R. T. Giaze, and A. Haliru, “Clinical investigation of treatment failure in type 2 diabetic patients treated with metformin and glibenclamide at a hospital in northwestern Nigeria,” Tropical Journal of Pharmaceutical Research, vol. 13, no. 9, p. 1521, 2014.
View at: Publisher Site | Google Scholar
A. K. Srivastava, “The role of fungus in bioactive compound production and nanotechnology,” Role of Plant Growth Promoting Microorganisms in Sustainable Agriculture and Nanotechnology, 2019.
View at: Google Scholar
L. Dufossé, M. Fouillaud, and Y. Caro, “Fungi and fungal metabolites for the improvement of human and animal nutrition and health,” Journal of Fungi, vol. 7, no. 4, p. 274, 2021.
View at: Publisher Site | Google Scholar
J. Drews, “Drug discovery today—and tomorrow,” Drug Discovery Today, vol. 5, no. 1, pp. 2–4, 2000.
View at: Publisher Site | Google Scholar
C. M. Song, S. J. Lim, and J. C. Tong, “Recent advances in computer-aided drug design,” Briefings in Bioinformatics, vol. 10, pp. 579–591, 2009.
View at: Publisher Site | Google Scholar
A. Lavecchia and C. Giovanni, “Virtual screening strategies in drug discovery: a critical review,” Current Medicinal Chemistry, vol. 20, no. 23, pp. 2839–2860, 2013.
View at: Publisher Site | Google Scholar
A. Lavecchia and C. Cerchia, “In silico methods to address polypharmacology: current status, applications and future perspectives,” Drug Discovery Today, vol. 21, no. 2, pp. 288–298, 2016.
View at: Publisher Site | Google Scholar
T. Cheng, Q. Li, Z. Zhou, Y. Wang, and S. H. Bryant, “Structure-based virtual screening for drug discovery: a problem-centric review,” The AAPS Journal, vol. 14, no. 1, pp. 133–141, 2012.
View at: Publisher Site | Google Scholar
S. Bibi and K. Sakata, “Current status of computer-aided drug design for type 2 diabetes,” Current Computer-Aided Drug Design, vol. 12, pp. 167–177, 2016.
View at: Publisher Site | Google Scholar
A. Gissi, D. Gadaleta, M. Floris et al., “An alternative QSAR-based approach for predicting the bioconcentration factor for regulatory purposes,” ALTEX, vol. 31, pp. 23–36, 2014.
View at: Publisher Site | Google Scholar
W. Yu and A. D. Mackerell, “Computer-aided drug design methods,” Methods in Molecular Biology, vol. 1520, 2017.
View at: Google Scholar
A. Bhatia, B. Singh, R. Arora, and S. Arora, “In vitro evaluation of the α-glucosidase inhibitory potential of methanolic extracts of traditionally used antidiabetic plants,” BMC Complementary and Alternative Medicine, vol. 19, no. 1, p. 74, 2019.
View at: Publisher Site | Google Scholar
B. Usman, N. Sharma, S. Satija et al., “Recent developments in alpha-glucosidase inhibitors for management of type-2 diabetes: an update,” Current Pharmaceutical Design, vol. 25, no. 23, pp. 2510–2525, 2019.
View at: Publisher Site | Google Scholar
M. S. Khan, B. Mehmood, Q. Yousafi et al., “Molecular docking studies reveal rhein from rhubarb (rheum rhabarbarum) as a putative inhibitor of ATP-binding cassette super-family G member 2,” Medicinal Chemistry, vol. 17, no. 3, pp. 273–288, 2021.
View at: Publisher Site | Google Scholar
I. Andricioaei and M. Karplus, “Particle mesh Ewald: an N log (N) method for Ewald sums in large systems,” Statistical Mechanics of Fluid Mixtures The Journal of Chemical Physics, vol. 115, 2001.
View at: Google Scholar
J. P. Ryckaert, G. Ciccotti, and H. J. C. Berendsen, “Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes,” Journal of Computational Physics, vol. 23, no. 3, pp. 327–341, 1977.
View at: Publisher Site | Google Scholar
R. Salomon-Ferrer, A. W. Götz, D. Poole, S. Le Grand, and R. C. Walker, “Routine microsecond molecular dynamics simulations with AMBER on GPUs 2 Explicit solvent particle mesh ewald,” Journal of Chemical Theory and Computation, vol. 9, pp. 3878–3888, 2013.
View at: Publisher Site | Google Scholar
M. F. Adasme, D. Parisi, K. Van Belle et al., “Structure-based drug repositioning explains ibrutinib as VEGFR2 inhibitor,” PLoS One, vol. 15, no. 5, Article ID e0233089, 2020.
View at: Publisher Site | Google Scholar
G. W. Bemis and M. A. Murcko, “Properties of known drugs 2 Side chains,” Journal of Medicinal Chemistry, vol. 42, no. 25, pp. 5095–5099, 1999.
View at: Publisher Site | Google Scholar
G. W. Bemis and M. A. Murcko, “The properties of known drugs 1 Molecular frameworks,” Journal of Medicinal Chemistry, vol. 39, no. 15, pp. 2887–2893, 1996.
View at: Publisher Site | Google Scholar
H. Y. Hung, K. Qian, S. L. Morris-Natschke, C. S. Hsu, and K. H. Lee, “Recent discovery of plant-derived anti-diabetic natural products,” Natural Product Reports, vol. 29, no. 5, p. 580, 2012.
View at: Publisher Site | Google Scholar
J. J. Salazar, W. J. Ennis, and T. J. Koh, “Diabetes medications: impact on inflammation and wound healing,” Journal of Diabetes and Its Complications, vol. 30, no. 4, pp. 746–752, 2016.
View at: Publisher Site | Google Scholar
K. Hong, G. Xu, T. B. Grayson, and A. Shalev, “Cytokines regulate β-Cell Thioredoxin-interacting Protein (TXNIP) via distinct mechanisms and pathways,” Journal of Biological Chemistry, vol. 291, no. 16, pp. 8428–8439, 2016.
View at: Publisher Site | Google Scholar
J. R. Abima Shazhni, A. Renu, and P. Vijayaraghavan, “Insights of antidiabetic, anti-inflammatory and hepatoprotective properties of antimicrobial secondary metabolites of corm extract from Caladium x hortulanum,” Saudi Journal of Biological Sciences, vol. 25, no. 8, pp. 1755–1761, 2018.
View at: Publisher Site | Google Scholar
S. Agrawal, S. Samanta, and S. K. Deshmukh, “The antidiabetic potential of endophytic fungi: future prospects as therapeutic agents,” Biotechnology and Applied Biochemistry, 2021.
View at: Publisher Site | Google Scholar
M. Jablonsky, A. Haz, Z. Burcova, F. Kreps, and J. Jablonsky, “Pharmacokinetic properties of biomass-extracted substances isolated by green solvents,” Bioresources, vol. 14, 2019.
View at: Publisher Site | Google Scholar
H. Chopra, S. Bibi, S. Kumar, M. S. Khan, P. Kumar, and I. Singh, “Preparation and evaluation of chitosan/PVA based hydrogel films loaded with honey for wound healing application,” Gels, vol. 8, no. 2, 2022.
View at: Publisher Site | Google Scholar
H. Gohlke, M. Hendlich, and G. Klebe, “Knowledge-based scoring function to predict protein-ligand interactions,” Journal of Molecular Biology, vol. 295, no. 2, pp. 337–356, 2000.
View at: Publisher Site | Google Scholar
M. Batool, B. Ahmad, and S. Choi, “A structure-based drug discovery paradigm,” International Journal of Molecular Sciences, vol. 20, no. 11, p. 2783, 2019.
View at: Publisher Site | Google Scholar
Z. Liu, Y. Liu, G. Zeng et al., “Application of molecular docking for the degradation of organic pollutants in the environmental remediation: a review,” Chemosphere, vol. 203, pp. 139–150, 2018.
View at: Publisher Site | Google Scholar
P. Ertl and B. Rohde, “The Molecule Cloud - compact visualization of large collections of molecules,” Journal of Cheminformatics, vol. 4, no. 1, p. 12, 2012.
View at: Publisher Site | Google Scholar
D. Stumpfe, H. Hu, and J. Bajorath, “Evolving concept of activity cliffs,” ACS Omega, vol. 4, no. 11, pp. 14360–14368, 2019.
View at: Publisher Site | Google Scholar
K. D. Bonin and V. V. Kresin, Electric-Dipole Polarizabilities of Atoms, Molecules, And Clusters, World Scientific, Singapore, 2010.
C. Hansch, W. E. Steinmetz, A. J. Leo, S. B. Mekapati, A. Kurup, and D. Hoekman, “On the role of polarizability in chemical - biological interactions,” Journal of Chemical Information and Computer Sciences, vol. 43, no. 1, pp. 120–125, 2003.
View at: Publisher Site | Google Scholar
P. Schneider and G. Schneider, “A computational method for unveiling the target promiscuity of pharmacologically active compounds,” Angewandte Chemie International Edition, vol. 56, no. 38, pp. 11520–11524, 2017.
View at: Publisher Site | Google Scholar
Z. R. Guo, “Drug promiscuity,” Yaoxue Xuebao, vol. 46, pp. 361–369, 2011.
View at: Google Scholar
S. K. Mencher and L. G. Wang, “Promiscuous drugs compared to selective drugs (promiscuity can be a virtue),” BMC Clinical Pharmacology, vol. 5, no. 1, p. 3, 2005.
View at: Publisher Site | Google Scholar
Y. Yang, H. Chen, I. Nilsson, S. Muresan, and O. Engkvist, “Investigation of the relationship between topology and selectivity for druglike molecules,” Journal of Medicinal Chemistry, vol. 53, no. 21, pp. 7709–7714, 2010.
View at: Publisher Site | Google Scholar
K. Azzaoui, J. Hamon, B. Faller et al., “Modeling promiscuity based on in vitro safety pharmacology profiling data,” ChemMedChem, vol. 2, no. 6, pp. 874–880, 2007.
View at: Publisher Site | Google Scholar
N. Strazielle and J. F. Ghersi-Egea, “Factors affecting delivery of antiviral drugs to the brain,” Reviews in Medical Virology, vol. 15, no. 2, pp. 105–133, 2005.
View at: Publisher Site | Google Scholar
H. Pajouhesh and G. R. Lenz, “Medicinal chemical properties of successful central nervous system drugs,” NeuroRx, vol. 2, no. 4, pp. 541–553, 2005.
View at: Publisher Site | Google Scholar
M. Karplus and J. A. McCammon, “Molecular dynamics simulations of biomolecules,” Nature Structural Biology, vol. 9, pp. 646–652, 2002.
View at: Publisher Site | Google Scholar
I. Aier, P. K. Varadwaj, and U. Raj, “Structural insights into conformational stability of both wild-type and mutant EZH2 receptor,” Scientific Reports, vol. 6, no. 1, Article ID 34984, 2016.
View at: Publisher Site | Google Scholar
U. Kalathiya, M. Padariya, and M. Baginski, “Structural, functional, and stability change predictions in human telomerase upon specific point mutations,” Scientific Reports, vol. 9, no. 1, p. 8707, 2019.
View at: Publisher Site | Google Scholar
Z. Jiang, L. You, W. Dou, T. Sun, and P. Xu, “Effects of an electric field on the conformational transition of the protein: a molecular dynamics simulation study,” Polymers, vol. 11, no. 2, p. 282, 2019.
View at: Publisher Site | Google Scholar
B. Hoff, S. Pöggeler, and U. Kück, “Eighty years after its discovery, Fleming’s Penicillium strain discloses the secret of its sex,” Eukaryotic Cell, vol. 7, no. 3, pp. 465–470, 2008.
View at: Publisher Site | Google Scholar
J. C. Frisvad and R. A. Samson, “Polyphasic taxonomy of Penicillium subgenus Penicillium: a guide to identification of food and air-borne terverticillate Penicillia and their mycotoxins,” Studies in Mycology, vol. 2004, 2004.
View at: Google Scholar
J. C. Frisvad, J. Smedsgaard, T. O. Larsen, and R. A. Samson, “Mycotoxins, drugs and other extrolites produced by species in Penicillium subgenus Penicillium,” Studies in Mycology, vol. 2004, 2004.
View at: Google Scholar

Copyright

Copyright © 2022 Shahzad Saleem 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.

PDF Download Citation

Download other formats

Order printed copies

Views

610

Downloads

531

Citations