In-silico study for African plants with possible beta-cell regeneration effect through inhibition of DYRK1A

1. Introduction

Diabetes is a life-threatening global health issue as a result of its high incidence [1], associated disability, and mortality [2]. The pancreatic beta-cell deficit is a significant part of the pathophysiological mechanism [3]. Beta-cells considerable damage leads to long-lasting endocrine insufficiency with the possibility of a permanent diabetic state. On the bright side, pancreatic beta-cell regeneration is a promising pharmacological strategy for recovering Beta-cells. In adults, it is known that the endocrine pancreas has a regulated ability for self-regeneration [4]. Consequently, approaches for stimulating beta-cell restoration have insightful inferences for the treatment and management of diabetes, particularly for type 1 diabetes and late-type 2 diabetes with considerable beta-cell loss.

Two possible approaches exist through which pancreatic beta-cells can be regenerated. The first approach is by preventing beta-cell loss precisely through the inhibition of beta-cell apoptosis and dedifferentiation. The second approach is to stimulate new endogenous regeneration and exogenous supplementation. For about a century, researchers have attempted pancreatic beta-cells regeneration. Under specific physiological environments, such as pregnancy, obesity, and conditions of insulin resistance, the adaption of islet and improved beta-cell mass take place in animal models and humans [5-8]. Contemporary advances in new technologies have offered additional substantiation on the generation of beta-cells. Single-cell RNA sequencing available data have revealed that human islets comprise four discrete subtypes of beta-cells [9] and probably transitional phases [10]. These suggest that beta-cells can acclimatize and undergo transdifferentiation or neogenesis. Physiological restoration research can make available data on the development of medication targeted toward beta-cell regeneration. Several approaches have been reported to be utilized in the promotion of beta-cells regeneration. The strategies include pancreatectomy, partial duct ligation, and chemical-induced massive beta-cell loss [11-15]. Molecular routes that cause multiplications in the mass of beta-cells have been comprehensively explored. Only a small portion of the materials studied have been found to have clinical, pre-clinical, or clinical potential as medicines. Thousands of materials have been studied, and hundreds are effective in the course of beta-cell restoration.

The CMGC (CDK, MAPK, CDC-like kinases, GSK3 kinase) family of eukaryotic protein kinases has been demonstrated to play crucial roles in neurodegenerative illnesses [16, 17], cell death, and tumorigenesis.  Dual-specificity tyrosine phosphorylation-regulated kinase A (DYRK1A) is a member of this family [18, 19]. A regulator of regeneration pathways essential to human insulin-producing pancreatic -cells has recently been discovered as DYRK1A [20-23]. Numerous studies have explored the development of DYRK1A inhibitor scaffolds, given the involvement of DYRK1A in these diseases [17-20, 22-24]. Several DYRK1A inhibitors from natural sources like harmine and small molecules have been identified and characterized [22, 25-48]. Harmine and its analogs (-carbolines) are the most often researched DYRK1A inhibitors, and they continue to be among the most effective and readily available inhibitor families that can be taken orally [17, 49]. The presence of harmine in the hallucinogenic infusion of ayahuasca and its affinity for serotonin, tryptamine, and other receptors in the central nervous system, in addition to its kinase inhibitory activity, have led to the hypothesis that harmine is a hallucinogen [50, 51]. Harmine and its analogs had been reported to block the DYRK1A-mediated phosphorylation of tau proteins in the CNS [52]. They also showed anti-proliferative action, including inhibition of topoisomerase I [53, 54], inhibition of CDKs [55], activation of cell apoptosis [56], and DNA intercalation [57].

This research aims to determine druggable enzyme/target/receptor that is vital in the pathogenesis of beta-cell apoptosis, identify phytocompounds with high binding affinity against the identified target using molecular docking simulation, and determine the drug-likeness of promising phytocompounds based on Lipinski’s rule, determine the toxicity of the phytocompounds in-silico, undertake bioactivity prediction of the phytocompounds on Molinspiration platform and identify the plant sources of the frontrunner compounds.

2.    Experimental Section

2.1. Materials

Personal computer, African Natural Compounds Database, PubChem (http://Pubchem.ncbi.nlm.nih.gov) [58], Linux operating system (Ubuntu desktop 18.04), Protein data bank (https://www.rcsb.org/) [59], DataWarrior software [60], PyMol software [61], AutoDockTools-1.5.6 software [62], Autodockvina 1.1.2 software [63], on Ubuntu operating system, Molinspiration Chemoinformatics web tool (https://www.molinspiration.com/cgi-bin/properties) [64].

2.2. Literature mining

Literature was mined to identify the target/receptor for possible induction of beta-cell regeneration. This was done to check the importance of the target/receptors in the onset and pathophysiology of beta-cell destruction. This gives more information about the receptor, functions, properties, and its druggability.

2.3. Selection and Preparation of the Receptor

After the identification of several targets/receptors, literature mining, and analysis of the target/receptor, Dual-specificity tyrosine phosphorylation-regulated kinase A in 3D format was obtained from Protein Data Bank (PDB) with the respective PDB code; 6UWY. The initial preparation of the PDB file to select the required chains, and delete multiple ligands was done with PyMol software. The PyMol software was employed to gain insight into the ligands binding to the receptors. The receptor was prepared for molecular docking simulations with the AutoDockTool. In the preparation, polar hydrogens and Kollman’s charges were added to the receptors and they were saved in the pdbqt file format. The pdbqt file format is the structural format required for the protein and ligand for carrying out molecular docking simulation. The electrostatic grid boxes and the 3-dimensional affinity of different sizes and centers, as indicated in Table 1 below, were created around the active site of the protein.

Table 1 Grid box parameters used for the molecular docking simulations

 

2.4. Selection, drug-likeness, and toxicity assessment of ligands (Phytocompounds)

A total number of 6511 isolated phytocompounds were obtained from the African Natural Products Database (African-compounds.org) [65, 66] in the 3D-structure data file format. The phytocompounds were loaded onto the DataWarrior software. Molecular properties such as molecular weight, hydrogen bond donor, hydrogen bond acceptor, partition coefficient (log P), and topological polar surface area (TPSA) were determined. Lipinski's rule of five violations was noted. The phytocompounds were also screened for toxicity (mutagenicity, carcinogenicity, tumorigenicity, and reproductive effect).

2.5. Selection and preparation of ligands

Phytocompounds following Lipinski’s rule of 5 with no toxicity in-silico were prepared for the molecular docking simulation. Reference ligands were identified from the literature including the compound co-crystallized with the receptor/protein on the PDB database. In preparation of the ligands for molecular docking simulation, all rotatable bonds, Torsions, and Geistegers charges were assigned and saved in the pdbqt file format.

2.6. Validation of docking protocol

To validate the molecular docking simulations protocol for the 6UWY (DYRK1A) protein, the PDB structure of this protein in complex with a reference inhibitor was reproduced in-silico. The deletion of the reference compound from the protein was done with the PyMol software. Polar hydrogen, Kollman charges, grid box sizes, and centers at a grid space of 1.0 Å were determined with the AutoDockTools-1.5.6 [62, 63]. The protein was saved in the pdbqt file format. The reference compound was prepared for molecular docking simulation with the AutoDockTools-1.5.6. All rotatable bonds, including torsions, were permitted to remain rotatable. Then, output was produced with the pdbqt file extension. Molecular docking simulation of the protein and reference compound was implemented locally with the AutoDockVina® [63] on a Linux platform using the centers and sizes with a virtual screening shell script. Docked conformations were visualized in the PyMol-1.4.1 software and the binding poses of the co-crystal inhibitors were compared with the re-docked co-crystal structures of the reference compound.

2.7. Molecular docking of the phytocompounds on Dual-specificity tyrosine phosphorylation-regulated kinase A

The phytocompounds were prepared in batches for molecular docking simulations using virtual screening scripts against the dual-specificity tyrosine phosphorylation-regulated kinase A. Following the validation of docking methods, four replicates of Molecular Docking Simulations were performed on a Linux platform using AutoDockVina® and related tools. To determine the leading phytocompounds, binding free energy values (kcal/mol SD) were ranked.

2.8. Bioactivity prediction of phytocompounds

The online Molinspiration web tool version 2011.06 (www.molinspiration.com) was supplied SMILES notations of the leading phytocompounds to forecast the bioactivity scores for kinase inhibition.

2.9. Calculation of the predicted percentage of absorption

The predicted percentage of absorption (% ab) of the frontrunner phytocompounds was calculated with the method reported by Zhao et al. (2002) [67]. The following formula was used: %ab = 109 – (0.345 x TPSA).

3. Results

3.1. Drug-likeness, and toxicity assessment of ligands (Phytocompounds)

The drug-likeness assessment of the 6511 phytocompounds based on Lipinski’s rule of five was done to screen out phytocompounds with violations of the rules. After the screening, a total number of 3814 phytocompounds had no violation of Lipinski’s rule, while 2697 phytocompounds violated the rules. Toxicity assessment on the 3814 phytocompounds that did not violate Lipinski's rule was carried out with DataWarrior in to identify phytocompounds that might be mutagenic, tumorigenic, irritating, or have reproductive consequences. A total number of 1897 phytocompounds were found to have none of the listed toxicities in-silico. Total polar surface area (TPSA) was also analyzed for all the phytocompounds.

3.2. Validation of docking protocol

The docking protocol validation was done to ensure in-silico reproducibility of the experimental protein-ligand interactions obtained from the protein data bank. The results obtained from the docking validations are presented below in Figures 1 and 2. Figure 1 represents the structural conformation and superimposition of the docked ligand (blue) and co-crystallized ligand (green) in the Dual-specificity tyrosine phosphorylation-regulated kinase A binding site. Figure 2A shows the 2D representation of the co-crystallized ligand-protein interaction, while Figure 2B shows the 2D representation of the docked ligand-protein interaction. Comparative analysis of the docked ligand and co-crystallized ligand-protein interaction reveals a 90.9% match.

Figure 1 Superimposed view of DYKR1A reference compound in blue and docked reference compound in green

Figure 2 2D representation of the (A) co-crystallized ligand-protein interaction and (B) the docked ligand-protein interaction

3.3. Molecular docking of the phytocompounds DYKR1A protein 

The molecular docking of the phytocompounds was performed on DYKR1A to identify phytocompounds with better in-silico inhibitory activity against DYKR1A than the reference compounds. The reference compounds are listed in the last three rows of Table 2. The docking was also performed to study the phytocompounds-proteins interaction pattern at the binding sites of these proteins. Phytocompounds with better binding affinities/energies than the reference compounds as can be observed from the mean binding affinity, are presented in Table 2.

Table 2 Phytocompounds with better binding energy values on DYRK1A than reference compounds

3.4. Bioactivity prediction of phytocompounds

Results of the bioactivity prediction of the 85 phytocompounds with better binding affinities than the reference compounds are presented in Table 3. The phytocompounds were screened for kinase inhibitory activity because the protein of interest (DYRK1A) is a kinase. Twelve phytocompounds were found to possess kinase inhibitory activity based on the scores. Some of the phytocompounds have better inhibitory scores than the reference compounds as can be observed in Table 3.

Tables 3 Bioactivity scores of DYRK1A active phytochemicals with their plant sources

3.5. Calculation of the predicted percentage of absorption

The results of the predicted percentage absorption of the frontrunner phytocompounds with that of the reference compounds are presented in Table 4. The prediction is based on the TPSA values.

Table 4 Predicted percentage of absorption

Figure 3 Structures of the frontrunner phytocompounds and reference compounds

4. Discussion

The study was set out to determine the binding affinities of phytocompounds from the African natural product database to DYRK1A compared to the reference compounds INDY, 4-(7-methoxy-1-methyl-9H-beta-carbolin-9-yl) butanamide and GNF4877 with in-silico molecular docking simulation. The process of developing new drugs has proven to be difficult due to the enormous expense of drug discovery and development, and the time needed. Modern research now relies heavily on computer-aided drug design, also known as the "in-silico" approach to drug discovery and design. Drug discovery and development are sped up by computer-aided drug design elements including molecular docking, molecular dynamics, QSAR, ADMET tool, and their accurate predictions.

On the other hand, medicines and medicinal substances have historically been derived from nature, primarily plants. The majority of medicines on the market today are either isolated or created from isolates derived from natural sources. Based on their use in conventional medical procedures, the majority of these currently utilized medications are made from natural sources [68]. To date, more novel compounds are being isolated from plants [69] and deposited in chemical databases. There are also general biological databases and specialized databases on which thousands of proteins are deposited to aid scientific research [70].

In this work, we retrieved 6511 phytocompounds from the African natural database that were purported to be isolated from African plants. Lipinski's rule of five was utilized to determine the drug-likeness of these phytocompounds. Pharmaceutical chemists frequently utilize Lipinski's rule of five to forecast the oral bioavailability of possible lead or therapeutic compounds during drug design and development. Lipinski's "rule of five" states that a candidate molecule is more likely to be orally active if it meets the following criteria: a) has a molecular weight below 500; b) has an estimated octanol/water partition coefficient (Log P is less than 5); c) has fewer than five hydrogen bond donors; and d) has fewer than ten hydrogen bond acceptors [71-73].

The DataWarrior software uses a precomputed collection of structural pieces that, when found in the structures under study, trigger toxicity alerts. All compounds from the Registry of Toxic Effects of Chemical Substances (RTECS) database [74] that are acknowledged to be active in a certain toxicity class were thoroughly destroyed to compile these fragment lists. During the shredding, compounds were first severed, with each rotating link leading to a set of core fragments. These, in turn, were utilized to reconstitute all possible significant substructures of the original molecule. The frequency of any fragment (core and created fragments) within all chemicals in that toxicity class was then determined using a substructure search process. Additionally, it identified these fragment frequencies in more than 3000 traded medications' structural data. Any fragment was viewed as a risk factor if it frequently occurs as the substructure of dangerous chemicals but never or very occasionally in the traded pharmaceuticals. This assumption was made based on the notion that sold drugs are primarily free of toxic effects. A total of 1897 phytocompounds did not exhibit any in-silico mutagenicity, tumorigenic, irritant, or reproductive impacts based on this explained fragments search. No fragments or fragments known to have any of the toxicities listed in the Registry of Toxic Effects of Chemical Substances were present in these phytocompounds.

From the molecular docking result, 85 phytocompounds were obtained with better binding affinity than the reference compounds, as shown in Table 2. Lower binding affinity suggests better ligand binding. The importance of binding affinity values is determined by the most significant magnitude negative value, representing the most favorable conformation of the complex formed when the ligand involved efficiently binds with the protein's active site. As observed, the mean binding affinity scores are in negative values. This is because protein-ligand binding only occurs spontaneously when the free energy change is negative, and the difference in ∆G levels of complexed and unbound free states is proportional to the stability of the protein-ligand interaction. Both protein folding and protein-ligand binding occur when ∆G is low in the system [75, 76]. Hence, negative ∆G scores indicate the stability of the resulting complexes with the receptor molecules, which is an essential characteristic of efficacious drugs [77].  

From the molinspiration bioactivity prediction, twelve compounds were found to be very active kinase inhibitors. Based on the prediction, two of the three reference compounds used were also very active kinase inhibitors. One of the reference compounds was predicted to be a moderately active kinase inhibitor. In molinspiration, biological activity is measured by a bioactivity score that is categorized as active (0.00 to 0.5), moderately active (0.00 to -0.5), and inactive (less than -0.5) [64].

The calculated percentage absorption (%ABS) of the frontrunner phytocompounds ranged between 66.12% and 101.85%, indicating that these phytocompounds have good permeability in the cellular membrane. The percentage absorption was calculated from the topological polar surface area (TPSA). The frontrunner phytocompounds exhibited computational TPSA values between 20.72 and 124.29 Ų and have good intestinal absorption. As a guide, orally active drugs transported by the transcellular route should not exceed a PSA of about 120 Ų [78,79]. Similarly, for good brain penetration of CNS drugs, this number should even be tailored to PSA<100Ų [79] or even smaller, <60–70 Ų [78].

Finally, observation of the frontrunner phytocompounds' structures compared with reference compounds, as presented in Figure 3, reveals some structural activity relationships that might be necessary for the inhibition of DYRK1A. The frontrunner compounds are composed of phenolics and alkaloids. From the 2D structure of the PDB reference compound presented in Figure 2, it can be observed that nitrogen, oxygen, and hydrogen atoms (which are all components of the frontrunner phytocompounds) are necessary for the protein-ligand interactions. Previous in-vitro research has shown that some natural products such as alkaloids and polyphenolic compounds act as inhibitors of DYRK1A. Epigallocatechin gallate, a major catechin component of green tea, when tested in a panel of 28 kinases structurally and functionally related to DYRK1A, showed selective inhibitory activity against DYRK1A (IC₅₀ 330 nM [ATP] = 100 μM) [80]. Acaninol B, which was isolated from the Leguminosae plant Acacia nilotica [81], exhibited only moderate action against DYRK1A (IC50 19 M [ATP] = 15 M) [82]. The already known CDK inhibitor flavopiridol was discovered to be a powerful DYRK1A inhibitor (IC50 0.3 M) [85] after being screened against a panel of five kinases using a set of natural and synthetic flavonoids and flavonoidal alkaloids. A strong DYRK1A inhibitor with an IC50 of 19 nM, staurosporine is an indolecarbazole derived from Streptomyces staurosporeus [86]. However, it is highly nonselective toward other kinases [87,88]. The L-rhamnulose-modified staurosporine analog was similarly significantly effective against DYRK1A (IC50 4 nM) [88]. Alkaloid acrifoline has a very strong DYRK1A inhibitory effect (IC50 = 0.075 M). Atalaphyllidine and chlorospermine B are both moderately effective DYRK1A inhibitors [89]. Recent studies have demonstrated the potency of two granulatimide analogs as DYRK1A inhibitors, with IC50 values of 0.26 and 0.09 M, respectively [90,91].

The results of in-silico studies translate well during in-vitro or in-vivo studies depending on the computational study design. In synthetic chemistry, the bioavailability, toxicity, and potential bioactivity of compound libraries are often studied with in-silico techniques before synthesis [92]. Similarly, we have prepared the study design of the present in-silico study to maximize the chance of obtaining fruitful results in bioassays. Molecular docking is one of the fastest and most reliable in-silico techniques to study the binding affinity, binding pose, and molecular interactions between a ligand and a protein [93]. Molecular docking studies are sometimes combined with in-vivo studies to study the molecular interaction and binding affinity of compounds with studied biomarkers [94, 95]. To tackle the severity of the ongoing coronavirus disease 2019 pandemic, cancer, and other infectious diseases, several researchers have designed their in-silico studies to maximize the chance of obtaining fruitful results in bioassays [96-104]. This discussion highlights the robustness of in-silico studies in drug discovery and development, corroborates and validates the study design used in the present in-silico study, and justifies that the African natural product database harbors promising phytocompounds as DYRK1A inhibitors.

5. Conclusions

Small-scale suppression of the DYRK1A molecules can offer a remedy for the pharmaceutical intervention of beta-cell regeneration in diabetes since options for treating beta-cell regeneration are a significant unmet therapeutic need. However, due to the conventional function of DYRK1A in controlling multiple signaling pathways vital to neuronal development and functions, caution should be used while trying to modulate it so that its activity is reduced to that which is typically seen in healthy individuals. These current in-silico tests' findings imply that 3-[6- (3-methyl-but-2-enyl) -1H-indolyl] -6-(3-methyl-but-2-enyl) -1H-indole, Lanceolatin B, Lysicamine, Pratorinine, Pratorimine, Lanceolatin A, Lanuginosine, Hippacine, (-)-Semiglabrin, Aegyptinone B, 3'-Prenylnaringenin and 8-C-p-hydroxybenzylluteolin are candidate ligands for activating beta-cells regeneration. The phytocompounds exhibit good intestine absorption, according to computational analyses of drug-likeness, TPSA, and % absorption. Finally, these phytocompounds have been recognized by the in-silico analysis as prospective novel medication candidates. Validating this in-silico work requires further thorough research using different models, such as in-vivo assays using the phytocompounds or extracts containing the phytocompounds.

Funding

Not applicable

Acknowledgment

We express our gratitude to the Principal Investigator, Drug Design and Informatics Group (DDIG) laboratory, for granting us free access to the facility. We appreciate the contribution of the CURIES research team of the Faculty of Pharmaceutical Sciences, Nnamdi Azikiwe University, Awka to this ongoing research.

Conflict of Interest

The authors declare no conflicting interests.

Author contributions

IMC: Study design and experiments; ICI: Figures and critical review; NMC: Figures and critical review; EII: Study design, experiments, and supervision

References