Next Article in Journal
The Accuracy of Patient-Specific Spinal Drill Guides Is Non-Inferior to Computer-Assisted Surgery: The Results of a Split-Spine Randomized Controlled Trial
Previous Article in Journal
Qualitative Interview Study of Gynecologic Oncologist Utilization of Recommended Same-Day Discharge Following Minimally Invasive Hysterectomy
Previous Article in Special Issue
Bioinformatics Analyses Identify the Therapeutic Potential of ST8SIA6 for Colon Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 12
Issue 7
10.3390/jpm12071083
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovering Tuberosin and Villosol as Potent and Selective Inhibitors of AKT1 for Therapeutic Targeting of Oral Squamous Cell Carcinoma

by
Mohd Adnan
1,
Deeba Shamim Jairajpuri
2,
Muskan Chaddha
3,
Mohd Shahnawaz Khan
4,
Dharmendra Kumar Yadav
5,*,
Taj Mohammad
3,
Abdelbaset Mohamed Elasbali
6,
Waleed Abu Al-Soud
7,8,
Salem Hussain Alharethi
9 and
Md. Imtaiyaz Hassan
3,*
1
Department of Biology, College of Science, University of Hail, Hail P.O. Box 2440, Saudi Arabia
2
Department of Medical Biochemistry, College of Medicine and Medical Sciences, Arabian Gulf University, Manama 26671, Bahrain
3
Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, Jamia Nagar, New Delhi 110025, India
4
Department of Biochemistry, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
5
College of Pharmacy, Gachon University of Medicine and Science, Hambakmoeiro, Yeonsu-gu, Incheon 21924, Korea
6
Department of Clinical Laboratory Science, College of Applied Sciences-Qurayyat, Jouf University, Sakaka 72388, Saudi Arabia
7
Department of Clinical Laboratory Sciences, Faculty of Applied Medical Sciences, Jouf University, Sakaka 72388, Saudi Arabia
8
Health Sciences Research Unit, Jouf University, Sakaka 72388, Saudi Arabia
9
Department of Biological Science, College of Arts and Science, Najran University, Najran 66252, Saudi Arabia
*
Authors to whom correspondence should be addressed.
J. Pers. Med. 2022, 12(7), 1083; https://doi.org/10.3390/jpm12071083
Submission received: 16 May 2022 / Revised: 23 June 2022 / Accepted: 27 June 2022 / Published: 30 June 2022
(This article belongs to the Special Issue Novel Targeted Cancer Therapeutics and Personalized Medicine)
Browse Figures
Versions Notes

Abstract

:
Oral squamous cell carcinoma (OSCC) is a major cause of death in developing countries because of high tobacco consumption. RAC-alpha serine-threonine kinase (AKT1) is considered as an attractive drug target because its prolonged activation and overexpression are associated with cancer progression and metastasis. In addition, several AKT1 inhibitors are being developed to control OSCC and other associated forms of cancers. We performed a screening of the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) database to discover promising AKT1 inhibitors which pass through various important filters such as ADMET (absorption, distribution, metabolism, excretion, and toxicity) properties, physicochemical properties, PAINS (pan-assay interference compounds) filters, PASS (prediction of activity spectra for substances) analysis, and specific interactions with AKT1. Molecules bearing admirable binding affinity and specificity towards AKT1 were selected for further analysis. Initially, we identified 30 natural compounds bearing appreciable affinity and specific interaction with AKT1. Finally, tuberosin and villosol were selected as potent and selective AKT1 inhibitors. To obtain deeper insights into binding mechanism and selectivity, we performed an all-atom molecular dynamics (MD) simulation and principal component analysis (PCA). We observed that both tuberosin and villosol strongly bind to AKT1, and their complexes were stable throughout the simulation trajectories. Our in-depth structure analysis suggested that tuberosin and villosol could be further exploited in the therapeutic targeting of OSCC and other cancers after further clinical validations.

1. Introduction

Cancer is a conglomeration of diseases that develop over a period of time and progress into a state of uncontrolled division of cells [1]. Such cells are functionally impaired and therefore cannot behave in the required manner. Such a state is harmful to physiological and molecular processes that sustain a living organism [2]. The occurrence of cancer is attributed to the alterations that take place in one or many processes that regulate cell division and cell death (apoptosis) [3]. Cancer requires constant communication between different factors and proteins that are involved in the regulation of tumor growth and vascular permeability that support its spread. Src family kinases (SFKs) control the vascular permeability of tumors along with tumor growth, metastasis, and endothelial-barrier regulation [4]. Similar to SFKs, protein kinase B also plays a crucial role in establishing communication between the factors promoting tumor growth and metastasis. Protein kinase B, or AKT, is a type of serine-threonine kinase [4,5].
The serine-threonine kinases have a key role in the management of processes such as programmed cell death, cell differentiation, cell proliferation, and embryonic development [6]. The protein kinase B exists in three isoforms that are AKT1, AKT2, and AKT3 [7]. Among them, AKT1 plays crucial functions in various activities related to the cell cycle, including, cell proliferation, development of blood vessels (angiogenesis), cell growth, and survival [8]. It also participates in maintaining glucose levels through insulin-induced translocation of the specific glucose transporter, SLC2A4/GLUT, to the surface of the cell [9]. The overexpression of AKT1 is linked to the development of oral squamous cell carcinoma (OSCC) [10,11] and many other complex diseases [12,13,14]. However, although it is involved in providing resistance against radiation, its inhibition is linked to the cessation of important pathways required for cell survival and overactivation of signal transduction pathways such as the receptor tyrosine kinase [15]. Studies show that reduced expression of AKT1 in knockdown mice also reduced the expression of TGFBR1 (transforming growth factor beta receptor 1) which is expressed in malignant phenotypes of OSCC [10].
The three AKT isoforms have also been studied and their role in the pathogenesis of cancer has been found through rigorous evaluation [16]. The intricate molecular network through which these isoforms work in triggering the oral squamous cell carcinoma can provide a means to effectively target and treat such challenging diseases [17]. The AKT1 protein isoforms are phosphorylated through the intervention of phosphoinositide 3-kinase, or PI3K [18]. The AKT1/PI3K signaling complex is involved in the regulation of processes such as angiogenesis [19]. Thus, this signaling complex forms a crucial part of various signaling cascades including G protein-coupled receptor (GPCR) signaling, receptor tyrosine pathway, and integrin-linked kinase [20]. AKT, therefore, has a crucial role in the induction of tumorigenesis by acting on different pathways that can lead to cancer development. Owing to its remarkable role in cancer progression, AKT1 represents itself as a potential druggable target in anticancer therapeutics [21,22].
High oxidative stress is linked to an increase in the rate of cell proliferation, metastasis, cell survival, and angiogenesis [23]. Cellular oxidative stress also hampers cell death signaling and increases resistance to biochemical drugs [24]. Therefore, reducing the reactive oxygen species and maintaining mitochondrial membrane potential can aid in decreased oxidative stress [25]. Phytoconstituents progressively emerged as safe and effective agents against anticancer drug development [26,27,28]. Tuberosin (5 hydroxy 3,4,7,3′,4′ pentamethoxy flavone) is extracted from the Indian medicinal plants Calopogonium mucunoides and Pueraria tuberosa [29]. It exhibits antioxidant properties, acts as a biochemical scavenger against hydroxyl radicals and superoxide radicals, and performs metal chelation [30]. Villosol, obtained from Tephrosia purpurea, possesses anti-inflammatory and anticancer properties and can serve as an important inhibitor of inflammation and oxidative damage that is linked to the development of cancer [31].
Although there are a few inhibitors already known for AKT1, still, new discoveries are required for the development of novel, safe, and effective drugs against AKT1-mediated malignancies [32]. Computational methods involving the use of virtual screening and molecular dynamics (MD) simulation are commendable for discovering potential inhibitors of predefined targets [33,34,35,36]. We performed a systematic structure-guided screening of natural compounds from the IMPPAT (Indian Medicinal Plants, Phytochemistry and Therapeutics) database to find out the most promising compounds that can help in the inhibition of AKT1, which is involved in the pathways leading to increased cell survival in cancer development.

2. Materials and Methods

2.1. Virtual Screening

The crystal structure of a protein is a prerequisite for molecular docking and, thus, the PDB file for the crystal structure of AKT1 was retrieved from the RCSB Protein Data Bank (PDB ID: 4GV1). The protein was assembled for virtual screening with the help of MGLTools [37], Swiss-PdbViewer [38], and InstaDock [39]. This step was crucial for filling gaps in the protein structure, adding hydrogen atoms to polar groups, and incorporating appropriate atoms. The molecular docking was conducted through InstaDock software under a blind search space. The visualization and interaction analysis was performed through DS Visualizer [40].
The phytoconstituents library was extracted from the IMPPAT database [41] after applying a filter following Lipinski’s rule of five [42]. The IMPPAT database contains a plethora of naturally derived plant compounds that possess pharmacological activity. The process of molecular docking was carried out to screen compounds that exhibit higher binding affinity towards the target protein AKT1. The compounds were filtered after using various segregating criteria, including binding affinity and ligand efficiency, and the most promising 30 candidates were taken ahead for further investigation. The compounds that exhibited acceptable affinity scores were processed and finalized to be investigated through InstaDock.

2.2. PAINS Patterns and ADMET Properties of Compounds

Safety and efficacy are two very important pillars in drug discovery and development [43]. The compounds that are to be used for treating and curing any ailment should possess specificity and selectivity so that they only target the process or the protein they are designed to affect [44,45,46,47]. Multiple attacking sites can pose a threat to the various vital processes in the body and can do more damage than inferring the desired effect. The compounds selected from the molecular-docking-based screening of the IMPPAT database were further investigated to carefully analyze their physicochemical properties. Top hits were collected from InstaDock and were taken for PAINS analysis [48]. The PAINS filter helps in identifying pan-assay interference compounds which might have a higher tendency to bind with multiple targets.
The compounds from the PAINS filter were further passed through an ADMET analysis which evaluates the absorption, digestion, metabolism, excretion, and toxicity of a chemical molecule. To carry out the ADMET analysis, the pkCSM [49] and SwissADME [50] servers were utilized. These software predict the pharmacokinetic and toxic properties of small molecules through graph-based signatures [51,52,53]. The software provides information on the efficacy of the compound through parameters such as total clearance, AMES toxicity, water solubility, skin permeability, hepatic toxicity, and so on. AMES toxicity helps to assess the probable carcinogenic effect of a compound through the use of the bacteria S. typhi. It is an effective and reliable method to estimate the genotoxicity of the compound and its propensity to reverse the mutations at the selected target in the bacterial strain [54,55,56]. We selected GI absorption, BBB permeation, CYP2D inhibitor, OCT2 substrate, AMES, and hepatotoxicity to finalize the compounds that showed the most promising results for effective molecules. The file format for pkCSM was selected as the SMILES format that was obtained from the Discovery Studio Visualizer software.

2.3. Biological Activity Prediction through PASS

The PASS (prediction of activity of spectra for substance) analysis is crucial for identifying the biological properties of a chemical compound [57]. The PASS analysis carefully assesses the biological and pharmacological properties of compounds that are already selected after PAINS and ADMET screening. The results that were obtained from the PASS analysis were represented in the form of two parameters, “Pa” or probability of being active and “Pi” or probability of being inactive. The ratio between the two values was also displayed to indicate the biological potential of the molecule. The molecules with a higher Pa value indicated a greater and better biological property.

2.4. Interaction Analysis

PyMOL was utilized for protein visualization and for exploring the possible interactions between the ligands and the shortlisted compounds. Apart from PyMOL, DS Visualizer served as an equally important tool for identifying probable interactions between the protein and ligands and in the identification of the important amino acid residues that participated in the formation of the binding pocket of the target protein.

2.5. MD Simulations

The key to a successful understanding of the protein and ligand relationship is found through the all-atom simulation studies [34,36,58]. The GROMACS simulation suite [59] enabled performance force-field-aided simulation in the AKT1 kinase and its complexes with the selected ligands villosol and tuberosin. Out of all possible conformations from the docked compounds, the top conformation of each compound was selected for simulation as they exhibited the most promising interactive results with the active sites of the target protein, AKT1. The ligand complex topologies were generated through the PRODRG server [60]. For solvation, we placed each system in a cubic box at a distance of 10Å from the center to the edges while utilizing the simple point-charge (SPC216) water model. The steepest descent algorithm provided energy minimization and position refinement, followed by NVT and NPT ensembles. The MD run was performed for 100 ns and resultant trajectories were analyzed through the QtGrace software [61].

3. Results and Discussion

3.1. Virtual Screening

The IMPPAT database served as a source of ~6000 compounds that conform to Lipinski’s rule of five (RO5) [41]. These compounds were selected after careful screening through advanced search filters such as Lipinski’s rule of five, the number of heavy atoms and heteroatoms, and stereochemical complexity. A molecular-docking-based virtual screening of filtered compounds resulted in the identification of 30 compounds selected based on their binding affinity scores (Table 1). High-affinity compounds might have a high tendency to bind in the active site pocket of AKT1 for its subsequent inhibition [55,56,62]. The docking study showed that all the selected hits exhibited appreciable binding affinity and ligand efficiency towards AKT1.

3.2. ADMET Properties and PAINS of the Compounds

The high-affinity compounds were subjected to the PAINS filter, where 16 out of 30 compounds were passed and further subjected to ADMET analysis. ADMET analysis aided in selecting the molecules based on their drug-like properties without any toxic moieties. The selected molecules that have an acceptable range of ADMET properties also possess a greater probability of being successful in clinical trials. Different ADMET parameters were obtained from the pkCSM software where only six compounds were selected based on their AMES and hepatotoxicity. The ADMET values of these six molecules are tabulated in Table 2. The molecules that show no toxicity qualify to be investigated further as AKT1 kinase inhibitors.

3.3. Screening for Biological Activity Prediction: PASS Analysis

The six selected compounds were subjected to PASS analysis for their biological properties. The PASS server (Way2drug) accepts the SMILES format for ligand files and displays results based on Pa and Pi values [57]. The activity predicted also served as an important criterion for selecting the compounds that possessed anticarcinogenic or antimutagenic properties. The six compounds showed acceptable results; however, the two compounds tuberosin and villosol exhibited the most promising outcomes from PASS analysis (Table 3). Both compounds possessed anticarcinogenic properties and, hence, were finalized for detailed interaction analysis.

3.4. Interaction Analysis

The protein AKT1 and the selected compounds tuberosin and villosol were subjected to interaction analysis. The interaction data for the reference compound AZD5363 [63] were utilized for pose selection and the different conformations of docked tuberosin and villosol were superposed on the structure of AKT1. The results showed that tuberosin and villosol interacted with the crucial residues of the AKT1 that formed the binding pocket for ATP and known inhibitors (Figure 1A). The ATP binding site Lys179 and the active site residue Asp274 interacted with both compounds. The magnified view of the protein–ligand interactions is illustrated in Figure 1B. Following the known sites, tuberosin and villosol make various close interactions with the critical sites on AKT1. Both compounds block the binding site for ATP and the active site with virtuous complementarity (Figure 1C). These interaction outcomes indicate that the stability of both the compounds obstructs ATP accessibility to the protein and, therefore, they function as potent inhibitors of the protein.
The data from the interaction analysis is taken further for a detailed study against the target protein. The elaborate studies in the interaction unveil different types of noncovalent interactions that are formed between AKT1 and tuberosin or villosol. The two-dimensional plots of the probable interactions were obtained for both the compounds using Discovery Studio Visualizer. Both compounds interact with the Lys179 and Asp 274 residues on the ATP binding site and active site, respectively (Figure 2). Both compounds share common features of binding and, thus, can act as ATP-competitive inhibitors of AKT1. The selected conformations of both compounds were taken for all-atom MD simulations for 100 ns.

3.5. MD Simulations

The MD simulations provide insights into the flexibility of proteins with and without the presence of ligands [34,36,58,64]. Along with evaluating the interaction and flexibility of proteins and ligands, MD simulation fills in the structural details and dynamic behavior of the proteins which were not accomplished through experimental interventions [65]. The docked complexes with AKT1 were subjected to simulation in an explicit solvent environment to analyze their structural dynamic properties and stability index as a function of time for a time period of 100 ns.

3.5.1. AKT1 Structural Changes through Ligand Binding

In order to assess the structural changes in the protein, the root-mean-square deviations (RMSDs) are crucial as they help to estimate the deviations in protein over time [35]. The alterations present in the backbone of the AKT1 were characterized by the RMSD analysis. The outcomes of RMSD analysis of the protein complexes AKT1-tuberosin and AKT1-villosol are represented in the form of time evolution plots (Figure 3). These plots show deviations occurring in the protein backbone and the results of these simulations are shown in Figure 3A. The data on simulation trajectory were utilized for analysis of the stability of the complexes. The data derived from the RMSD plot showed a fair consistency with small fluctuations. The deviations were reduced in the ligand–protein systems, especially in the AKT1-yuberosin complex. The three systems of apo-AKT1, AKT1-yuberosin, and AKT1-villosol reached equilibrium and endured stable conformation through the simulation trajectory of 100 ns (Figure 3A). The presence of random alterations in the AKT1-villosol complex’s RMSD of approximately 0.1 nanometers is distributed at 40–50 ns, which indicates a small adjustment in the system. The RMSD plots as a probability density function (PDF) depict lower RMSD and stabilization of AKT1 dynamics when the protein is present in the form of the complex. Consequently, the inference from the RMSD analysis is that the ligand-bound protein complex is folded and stable during the simulation.
The root-mean-square fluctuations (RMSFs) assist in the estimation of the residual vibrations present in the protein molecule at the time of the MD simulation [66,67]. The protein backbone was subjected to RMSF studies to evaluate the flexibility of the residues independently (Figure 3B). The figure depicts that the three systems show a similar pattern of RMSF and the residual fluctuations are minimized and stabilized when AKT1 is bound to tuberosin and villosol. Through RMSF intervention it was observed that residues in the protein binding pocket are almost stable and only show slight fluctuation at the time of the simulation. Only a few residues showed a deviation higher than the other systems, implying initial adjustment of the binding as seen in the RMSD analysis. Any major fluctuations in the AKT1 structure are only local to the loop regions and, comprehensively, tuberosin and villosol show stable interactions with AKT1.
The radius of gyration (Rg) is the parameter associated with the tertiary structure of a protein [68]. It indicates the root-mean-square distance of the collection of atoms from the collective center of mass. The parameter is extensively used as the means to evaluate the compactness seen in the tertiary structure of the protein. The determination of Rg is conducted on the basis of a time scale for apo and protein–ligand complex systems. After a rigorous assessment of the three systems, the plot indicates that the protein AKT1 in the presence of the compounds tuberosin and villosol shifted from the stability range of 2.05 nm to 1.90 nm during the simulation trajectory. The shift in the values is depicted in Figure 4.
Following interaction analysis and estimation of Rg, the studies on the surface area of the protein and its accessibility to the adjacent solvent form an important part of the protein structure analysis [58]. The SASA, or solvent-accessible surface area, helps to study the stability of the protein structure and its folding patterns. Upon SASA analysis, we can understand the number of native contacts that are formed by the protein. The time-based evaluation of the SASA values against the protein AKT1 in its free and complex form were plotted (Figure 4B). The plot depicts no significant fluctuations in the SASA values during the entire simulation, even when bound to the two ligand molecules (Figure 4B). The overall compactness and the protein folding remain undisturbed as suggested by SASA analysis and, while correlating the values with Rg, the SASA distribution illustrates a similar equilibration pattern. The PDF distribution pattern of the AKT1 and its complex with tuberosin and villosol are shown in the lower panel of Figure 4B.

3.5.2. Dynamics of Hydrogen Bonds

Protein structures are dynamic in nature and are stabilized by the formation of crucial bonded interactions, i.e., intramolecular hydrogen bonds (H-bonds) and hydrophobic interactions [69,70]. The hydrogen bonds between the two participate in the folding of the protein and determining the stable protein conformation. The MD trajectories helped to estimate the formation of H-bonds in time evolution form. This analysis was essential in understanding the consistency of the intramolecular H-bonds in AKT1 when it is present in a bounded state. The changes in the number of H-bonds before and after binding with both compounds show that these bonds are more persistent and are significant in the geometry of the AKT1 structure. A slight increment in the number of intramolecular hydrogen bonds from the apostate to the bound state of protein AKT1 with tuberosin and villosol indicates that there is an increase in the compactness of the protein structure’s inbound state, thereby making it dynamically active. The PDF for the comparative change in the hydrogen bonds is depicted in Figure 5.
The intramolecular H-bonds were studied through time evolution and the persistence of these bonds between the selected compounds and AKT1 was further investigated. The average number of H-bonds formed between the protein–ligand complexes is one (Figure 6). The presence of acceptable consistency in any system is shown by a higher PDF value and an average of one H-bond; such values are observed in AKT1-tuberosin and AKT1-villosol. There was no major change observed in the initial docking position of tuberosin and villosol, which proves the stability of both complexes.

3.6. Principal Components Analysis

The principal component analysis (PCA) provides information on the conformational sampling and the constitutive motions of the protein as seen in the simulated trajectories [71]. The PCA studies were performed on AKT1 and its complexes through the essential dynamics approach. The utilization of Cα atoms during the sampling resulted in generating the contour maps, as illustrated in Figure 7. The plot conveys that the AKT1-tuberosin and AKT1-villosol complexes are present in similar subspaces, as shown in both the EVs. The restricted movement in the essential subspace of AKT1-tuberosin with reduced flexibility can be seen on the EV2, which further supports the observations that the docked complexes are stable during the simulation studies (Figure 7).

3.7. Free Energy Landscapes

The free energy landscapes, or the FELs, are a vital means to explain the mechanism and patterns of protein folding [72]. The plots help to examine the levels of protein stability in the presence and absence of ligands, as well as under different solvent conditions during the MD simulation studies. The plots for the FELs for the three systems were generated and show the energy minima and conformational landscape in two principal components, i.e., PC1 and PC2. The plots are shown in the form of contoured maps represented in different colours, indicating the levels of the protein folding state. Figure 8 depicts that the binding of AKT1 to tuberosin and villosol causes fluctuations in the size and the positions of the stable global minima. The areas shown in dark blue denote the stable conformations with lower energy as the conformation is close to the native state of the protein (Figure 8). The unbound protein AKT1 is restricted to a large single global minima which expands to two or three basins. The AKT1-tuberosin and AKT1-villosol (Figure 8B–C) also show similar areas of global minima but with varying populations. Comprehensively, results obtained through the essential dynamics of AKT1 in free and bound states with tuberosin and villosol portray stability with minimal conformational switching during the simulations.

4. Conclusions

Cancer has become one of the most challenging diseases, as it requires constant advancement in medical sciences and a better understanding of the underlying mechanisms of disease development and progression. Anticancer therapies are constantly evolving towards more efficient and targeted specific treatments. Various proteins and metabolic processes play essential roles in the development of cancer. AKT1 acts as a positive regulator of different types of cancers, especially OSCC. The pharmacological inhibition of AKT1 using tuberosin and villosol can be a novel strategy to develop more effective anticancer molecules. The computational analysis using state-of-the-art approaches has helped in identifying these two phytochemicals that have potential therapeutic features and show a significantly high affinity for AKT1. Tuberosin and villosol possess anticancer properties and further investigation of these compounds can contribute to the management of cancer. Thus, we present a proposition that tuberosin and villosol should be taken forward for further investigation in in vitro and in vivo settings for the development of better anticancer therapeutic regimes. This study also aimed to draw attention to the anticancer effects of traditional medicinal plants and their active components, promoting their use in developing potent anticancer therapies.

Author Contributions

Conceptualization, M.A., D.S.J. and T.M.; methodology, A.M.E., M.C., D.K.Y. and T.M.; software, W.A.A.-S., S.H.A. and A.M.E.; validation, M.A., D.K.Y. and D.S.J.; formal analysis, M.C. and M.I.H.; investigation, D.S.J., M.S.K., T.M. and S.H.A.; resources, M.A. and M.I.H.; data curation, T.M. and M.S.K.; writing—original draft preparation, M.A., T.M. and D.S.J.; writing—review and editing, M.S.K., D.K.Y. and M.I.H.; visualization, M.I.H.; supervision, M.I.H.; project administration, M.A. and D.K.Y.; funding acquisition, D.S.J., M.S.K. and D.K.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Indian Council of Medical Research (grant number ISRM/12(22)/2020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

M.S.K. acknowledges the generous support from the Research Supporting Project (RSP-2021/352)) by King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Giat, E.; Ehrenfeld, M.; Shoenfeld, Y. Cancer and autoimmune diseases. Autoimmun. Rev. 2017, 16, 1049–1057. [Google Scholar] [CrossRef] [PubMed]
  2. Carmeliet, P.; Jain, R.K. Angiogenesis in cancer and other diseases. Nature 2000, 407, 249–257. [Google Scholar] [CrossRef] [PubMed]
  3. Lawen, A. Apoptosis—An introduction. Bioessays 2003, 25, 888–896. [Google Scholar] [CrossRef] [PubMed]
  4. Alwhaibi, A.; Verma, A.; Adil, M.S.; Somanath, P.R. The unconventional role of Akt1 in the advanced cancers and in diabetes-promoted carcinogenesis. Pharmacol. Res. 2019, 145, 104270. [Google Scholar] [CrossRef]
  5. Alwhaibi, A.; Verma, A.; Artham, S.; Adil, M.S.; Somanath, P.R. Nodal pathway activation due to Akt1 suppression is a molecular switch for prostate cancer cell epithelial-to-mesenchymal transition and metastasis. Biochem. Pharmacol. 2019, 168, 1–13. [Google Scholar] [CrossRef]
  6. Mishra, A.P.; Salehi, B.; Sharifi-Rad, M.; Pezzani, R.; Kobarfard, F.; Sharifi-Rad, J.; Nigam, M. Programmed cell death, from a cancer perspective: An overview. Mol. Diagn. Ther. 2018, 22, 281–295. [Google Scholar] [CrossRef]
  7. Cozzone, D.; Fröjdö, S.; Disse, E.; Debard, C.; Laville, M.; Pirola, L.; Vidal, H. Isoform-specific defects of insulin stimulation of Akt/protein kinase B (PKB) in skeletal muscle cells from type 2 diabetic patients. Diabetologia 2008, 51, 512–521. [Google Scholar] [CrossRef] [Green Version]
  8. Somanath, P.R.; Razorenova, O.V.; Chen, J.; Byzova, T.V. Akt1 in endothelial cell and angiogenesis. Cell Cycle 2006, 5, 512–518. [Google Scholar] [CrossRef] [Green Version]
  9. Tremblay, F.; Lavigne, C.; Jacques, H.; Marette, A. Defective insulin-induced GLUT4 translocation in skeletal muscle of high fat–fed rats is associated with alterations in both Akt/protein kinase B and atypical protein kinase C (ζ/λ) activities. Diabetes 2001, 50, 1901–1910. [Google Scholar] [CrossRef] [Green Version]
  10. Nakashiro, K.-I.; Tanaka, H.; Goda, H.; Iwamoto, K.; Tokuzen, N.; Hara, S.; Onodera, J.; Fujimoto, I.; Hino, S.; Hamakawa, H. Identification of Akt1 as a potent therapeutic target for oral squamous cell carcinoma. Int. J. Oncol. 2015, 47, 1273–1281. [Google Scholar] [CrossRef] [Green Version]
  11. Sharif Siam, M.K.; Sarker, A.; Sayeem, M.M.S. In silico drug design and molecular docking studies targeting Akt1 (RAC-alpha serine/threonine-protein kinase) and Akt2 (RAC-beta serine/threonine-protein kinase) proteins and investigation of CYP (cytochrome P450) inhibitors against MAOB (monoamine oxidase B) for OSCC (oral squamous cell carcinoma) treatment. J. Biomol. Struct. Dyn. 2021, 39, 6467–6479. [Google Scholar] [PubMed]
  12. Wang, Y.W.; Lin, K.T.; Chen, S.C.; Gu, D.L.; Chen, C.F.; Tu, P.H.; Jou, Y.S. Overexpressed-eIF3I interacted and activated oncogenic Akt1 is a theranostic target in human hepatocellular carcinoma. Hepatology 2013, 58, 239–250. [Google Scholar] [CrossRef] [PubMed]
  13. Priolo, C.; Pyne, S.; Rose, J.; Regan, E.R.; Zadra, G.; Photopoulos, C.; Cacciatore, S.; Schultz, D.; Scaglia, N.; McDunn, J. AKT1 and MYC induce distinctive metabolic fingerprints in human prostate cancer. Cancer Res. 2014, 74, 7198–7204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Sun, M.; Wang, G.; Paciga, J.E.; Feldman, R.I.; Yuan, Z.-Q.; Ma, X.-L.; Shelley, S.A.; Jove, R.; Tsichlis, P.N.; Nicosia, S.V. AKT1/PKBα kinase is frequently elevated in human cancers and its constitutive activation is required for oncogenic transformation in NIH3T3 cells. Am. J. Pathol. 2001, 159, 431–437. [Google Scholar] [CrossRef]
  15. Paoli, P.; Giannoni, E.; Chiarugi, P. Anoikis molecular pathways and its role in cancer progression. Biochim. Biophys. Acta-Mol. Cell Res. 2013, 1833, 3481–3498. [Google Scholar] [CrossRef] [Green Version]
  16. Alzahrani, A.S. PI3K/Akt/mTOR inhibitors in cancer: At the bench and bedside. Semin. Cancer Biol. 2019, 59, 125–132. [Google Scholar] [CrossRef]
  17. Oudit, G.Y.; Sun, H.; Kerfant, B.-G.; Crackower, M.A.; Penninger, J.M.; Backx, P.H. The role of phosphoinositide-3 kinase and PTEN in cardiovascular physiology and disease. J. Mol. Cell. Cardiol. 2004, 37, 449–471. [Google Scholar] [CrossRef]
  18. Liu, P.; Cheng, H.; Roberts, T.M.; Zhao, J.J. Targeting the phosphoinositide 3-kinase pathway in cancer. Nat. Rev. Drug Discov. 2009, 8, 627–644. [Google Scholar] [CrossRef] [Green Version]
  19. Karar, J.; Maity, A. PI3K/AKT/mTOR pathway in angiogenesis. Front. Mol. Neurosci. 2011, 4, 51. [Google Scholar] [CrossRef] [Green Version]
  20. Jiang, B.H.; Liu, L.Z. PI3K/PTEN signaling in angiogenesis and tumorigenesis. Adv. Cancer Res. 2009, 102, 19–65. [Google Scholar] [CrossRef]
  21. Mitsiades, C.S.; Mitsiades, N.; Koutsilieris, M. The Akt pathway: Molecular targets for anti-cancer drug development. Curr. Cancer Drug Targets 2004, 4, 235–256. [Google Scholar] [CrossRef] [PubMed]
  22. Cheng, J.Q.; Lindsley, C.W.; Cheng, G.Z.; Yang, H.; Nicosia, S.V. The Akt/PKB pathway: Molecular target for cancer drug discovery. Oncogene 2005, 24, 7482–7492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Reuter, S.; Gupta, S.C.; Chaturvedi, M.M.; Aggarwal, B.B. Oxidative stress, inflammation, and cancer: How are they linked? Free. Radic. Biol. Med. 2010, 49, 1603–1616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yun, H.R.; Jo, Y.H.; Kim, J.; Shin, Y.; Kim, S.S.; Choi, T.G. Roles of autophagy in oxidative stress. Int. J. Mol. Sci. 2020, 21, 3289. [Google Scholar] [CrossRef]
  25. Trachootham, D.; Alexandre, J.; Huang, P. Targeting cancer cells by ROS-mediated mechanisms: A radical therapeutic approach? Nat. Rev. Drug Discov. 2009, 8, 579–591. [Google Scholar] [CrossRef]
  26. Dutt, R.; Garg, V.; Khatri, N.; Madan, A.K. Phytochemicals in anticancer drug development. Anti-Cancer Agents Med. Chem. 2019, 19, 172–183. [Google Scholar] [CrossRef]
  27. Anjum, F.; Mohammad, T.; Almalki, A.A.; Akhtar, O.; Abdullaev, B.; Hassan, M.I. Phytoconstituents and Medicinal plants for anticancer drug discovery: Computational identification of potent inhibitors of PIM1 kinase. OMICS J. Integr. Biol. 2021, 25, 580–590. [Google Scholar] [CrossRef]
  28. Ashraf, M.A. Phytochemicals as potential anticancer drugs: Time to ponder nature’s bounty. BioMed Res. Int. 2020, 2020, 8602879. [Google Scholar] [CrossRef]
  29. Joshi, B.S.; Kamat, V.N. Tuberosin, a new pterocarpan from Pueraria tuberosa DC. J. Chem. Soc. Perkin Trans. 1973, 907–911. [Google Scholar] [CrossRef]
  30. Pandey, N.; Tripathi, Y.B. Antioxidant activity of tuberosin isolated from Pueraria tuberose Linn. J. Inflamm. 2010, 7, 47. [Google Scholar] [CrossRef] [Green Version]
  31. Li, Q.; Dong, D.-D.; Huang, Q.-P.; Li, J.; Du, Y.-Y.; Li, B.; Li, H.-Q.; Huyan, T. The anti-inflammatory effect of Sonchus oleraceus aqueous extract on lipopolysaccharide stimulated RAW 264.7 cells and mice. Pharm. Biol. 2017, 55, 799–809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Zhao, Z.; Leister, W.H.; Robinson, R.G.; Barnett, S.F.; Defeo-Jones, D.; Jones, R.E.; Hartman, G.D.; Huff, J.R.; Huber, H.E.; Duggan, M.E. Discovery of 2, 3, 5-trisubstituted pyridine derivatives as potent Akt1 and Akt2 dual inhibitors. Bioorg. Med. Chem. Lett. 2005, 15, 905–909. [Google Scholar] [CrossRef] [PubMed]
  33. Mohammad, T.; Khan, F.I.; Lobb, K.A.; Islam, A.; Ahmad, F.; Hassan, M.I. Identification and evaluation of bioactive natural products as potential inhibitors of human microtubule affinity-regulating kinase 4 (MARK4). J. Biomol. Struct. Dyn. 2019, 37, 1813–1829. [Google Scholar] [CrossRef] [PubMed]
  34. Mohammad, T.; Siddiqui, S.; Shamsi, A.; Alajmi, M.F.; Hussain, A.; Islam, A.; Ahmad, F.; Hassan, M. Virtual screening approach to identify high-affinity inhibitors of serum and glucocorticoid-regulated kinase 1 among bioactive natural products: Combined molecular docking and simulation studies. Molecules 2020, 25, 823. [Google Scholar] [CrossRef] [Green Version]
  35. Ali, S.; Khan, F.I.; Mohammad, T.; Lan, D.; Hassan, M.; Wang, Y. Identification and evaluation of inhibitors of lipase from Malassezia restricta using virtual high-throughput screening and molecular dynamics studies. Int. J. Mol. Sci. 2019, 20, 884. [Google Scholar] [CrossRef] [Green Version]
  36. Naqvi, A.A.; Mohammad, T.; Hasan, G.M.; Hassan, M. Advancements in docking and molecular dynamics simulations towards ligand-receptor interactions and structure-function relationships. Curr. Top. Med. Chem. 2018, 18, 1755–1768. [Google Scholar] [CrossRef]
  37. Morris, G.M.; Huey, R.; Olson, A.J. Using autodock for ligand-receptor docking. Curr. Protoc. Bioinform. 2008, 24, 8–14. [Google Scholar] [CrossRef]
  38. Guex, N.; Peitsch, M.C. SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis 1997, 18, 2714–2723. [Google Scholar] [CrossRef]
  39. Mohammad, T.; Mathur, Y.; Hassan, M.I. InstaDock: A single-click graphical user interface for molecular docking-based virtual high-throughput screening. Brief. Bioinform. 2021, 22, 279. [Google Scholar] [CrossRef]
  40. Discovery Studio; Version 2.1; Accelrys: San Diego, CA, USA, 2008.
  41. Mohanraj, K.; Karthikeyan, B.S.; Vivek-Ananth, R.; Chand, R.B.; Aparna, S.; Mangalapandi, P.; Samal, A. IMPPAT: A curated database of I ndian M edicinal P lants, P hytochemistry A nd T herapeutics. Sci. Rep. 2018, 8, 4329. [Google Scholar] [CrossRef] [Green Version]
  42. Lipinski, C.A. Lead-and drug-like compounds: The rule-of-five revolution. Drug Discov. Today Technol. 2004, 1, 337–341. [Google Scholar] [CrossRef] [PubMed]
  43. Reichel, A.; Lienau, P. Pharmacokinetics in drug discovery: An exposure-centred approach to optimising and predicting drug efficacy and safety. In New Approaches to Drug Discover; Springer: Cham, Switzerland, 2015; pp. 235–260. [Google Scholar]
  44. Yadav, D.K.; Kumar, S.; Choi, E.-H.; Chaudhary, S.; Kim, M.-H. Computational modeling on aquaporin-3 as skin cancer target: A virtual screening study. Front. Chem. 2020, 8, 250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Waseem, R.; Anwar, S.; Khan, S.; Shamsi, A.; Hassan, M.; Anjum, F.; Shafie, A.; Islam, A.; Yadav, D.K. MAP/Microtubule affinity regulating kinase 4 inhibitory potential of irisin: A new therapeutic strategy to combat cancer and Alzheimer’s disease. Int. J. Mol. Sci. 2021, 22, 10986. [Google Scholar] [CrossRef]
  46. Cheon, J.H.; Kim, K.S.; Yadav, D.K.; Kim, M.; Kim, H.S.; Yoon, S. The JAK2 inhibitors CEP-33779 and NVP-BSK805 have high P-gp inhibitory activity and sensitize drug-resistant cancer cells to vincristine. Biochem. Biophys. Res. Commun. 2017, 490, 1176–1182. [Google Scholar] [CrossRef] [PubMed]
  47. Yadav, D.K. Recent advances on small molecule medicinal chemistry to treat human diseases-Part III. Curr. Top. Med. Chem. 2021, 21, 1517–1518. [Google Scholar] [CrossRef]
  48. Baell, J.B. Feeling nature’s PAINS: Natural products, natural product drugs, and pan assay interference compounds (PAINS). J. Nat. Prod. 2016, 79, 616–628. [Google Scholar] [CrossRef]
  49. Pires, D.E.; Blundell, T.L.; Ascher, D.B. pkCSM: Predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures. J. Med. Chem. 2015, 58, 4066–4072. [Google Scholar] [CrossRef]
  50. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  51. Tiwari, M.K.; Coghi, P.; Agrawal, P.; Shyamlal, B.R.K.; Yang, L.J.; Yadav, L.; Peng, Y.; Sharma, R.; Yadav, D.K.; Sahal, D. Design, synthesis, structure-activity relationship and docking studies of novel functionalized arylvinyl-1, 2, 4-trioxanes as potent antiplasmodial as well as anticancer agents. ChemMedChem 2020, 15, 1216–1228. [Google Scholar] [CrossRef]
  52. Kalani, K.; Yadav, D.K.; Alam, S.; Khan, F.; Kashyap, M.P.; Srivastava, S.K.; Pant, A.B. In-silico studies and wet-lab validation of camptothecin derivatives for anti-cancer activity against liver (HepG2) and lung (A549) cancer cell lines. Curr. Top. Med. Chem. 2021, 21, 908–919. [Google Scholar] [CrossRef]
  53. Teli, M.K.; Kumar, S.; Yadav, D.K.; Kim, M.-h. In silico identification of hydantoin derivatives: A novel natural prolyl hydroxylase inhibitor. J. Biomol. Struct. Dyn. 2021, 39, 703–717. [Google Scholar] [CrossRef] [PubMed]
  54. Shyamlal, B.R.K.; Yadav, L.; Tiwari, M.K.; Mathur, M.; Prikhodko, J.I.; Mashevskaya, I.V.; Yadav, D.K.; Chaudhary, S. Synthesis, bioevaluation, structure-activity relationship and docking studies of natural product inspired (Z)-3-benzylideneisobenzofuran-1 (3H)-ones as highly potent antioxidants and antiplatelet agents. Sci. Rep. 2020, 10, 2307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Teli, M.K.; Kumar, S.; Yadav, D.K.; Kim, M.h. In silico identification of prolyl hydroxylase inhibitor by per-residue energy decomposition-based pharmacophore approach. J. Cell. Biochem. 2021, 122, 1098–1112. [Google Scholar] [CrossRef] [PubMed]
  56. Yadav, D.K.; Kumar, S.; Teli, M.K.; Kim, M.H. Ligand-based pharmacophore modeling and docking studies on vitamin D receptor inhibitors. J. Cell. Biochem. 2020, 121, 3570–3583. [Google Scholar] [CrossRef] [PubMed]
  57. Lagunin, A.; Stepanchikova, A.; Filimonov, D.; Poroikov, V. PASS: Prediction of activity spectra for biologically active substances. Bioinformatics 2000, 16, 747–748. [Google Scholar] [CrossRef]
  58. Anjum, F.; Ali, F.; Mohammad, T.; Shafie, A.; Akhtar, O.; Abdullaev, B.; Hassan, M.I. Discovery of Natural Compounds as Potential Inhibitors of Human Carbonic Anhydrase II: An Integrated Virtual Screening, Docking, and Molecular Dynamics Simulation Study. OMICS J. Integr. Biol. 2021, 25, 513–524. [Google Scholar] [CrossRef]
  59. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  60. Schüttelkopf, A.W.; Van Aalten, D.M. PRODRG: A tool for high-throughput crystallography of protein–ligand complexes. Acta Crystallogr. Sect. D Biol. Crystallogr. 2004, 60, 1355–1363. [Google Scholar] [CrossRef] [Green Version]
  61. Turner, P. XMGRACE: Center for Coastal and Land-Margin Research; Version 5.1. 19; Oregon Graduate Institute of Science and Technology: Beaverton, OR, USA, 2005. [Google Scholar]
  62. Yadav, D.K.; Kumar, S.; Misra, S.; Yadav, L.; Teli, M.; Sharma, P.; Chaudhary, S.; Kumar, N.; Choi, E.H.; Kim, H.S. Molecular insights into the interaction of RONS and Thieno [3,2-c] pyran analogs with SIRT6/COX-2: A molecular dynamics study. Sci. Rep. 2018, 8, 4777. [Google Scholar] [CrossRef]
  63. Addie, M.; Ballard, P.; Buttar, D.; Crafter, C.; Currie, G.; Davies, B.R.; Debreczeni, J.; Dry, H.; Dudley, P.; Greenwood, R. Discovery of 4-Amino-N-[(1S)-1-(4-chlorophenyl)-3-hydroxypropyl]-1-(7 H-pyrrolo [2,3-d] pyrimidin-4-yl) piperidine-4-carboxamide (AZD5363), an Orally Bioavailable, Potent Inhibitor of Akt Kinases. J. Med. Chem. 2013, 56, 2059–2073. [Google Scholar] [CrossRef]
  64. Khan, A.; Mohammad, T.; Shamsi, A.; Hussain, A.; Alajmi, M.F.; Husain, S.A.; Iqbal, M.A.; Hassan, M.I. Identification of plant-based hexokinase 2 inhibitors: Combined molecular docking and dynamics simulation studies. J. Biomol. Struct. Dyn. 2021, 2021, 1942217. [Google Scholar] [CrossRef] [PubMed]
  65. Amir, M.; Mohammad, T.; Prasad, K.; Hasan, G.M.; Kumar, V.; Dohare, R.; Islam, A.; Ahmad, F.; Imtaiyaz Hassan, M. Virtual high-throughput screening of natural compounds in-search of potential inhibitors for protection of telomeres 1 (POT1). J. Biomol. Struct. Dyn. 2020, 38, 4625–4634. [Google Scholar] [CrossRef]
  66. Naz, F.; Khan, F.I.; Mohammad, T.; Khan, P.; Manzoor, S.; Hasan, G.M.; Lobb, K.A.; Luqman, S.; Islam, A.; Ahmad, F. Investigation of molecular mechanism of recognition between citral and MARK4: A newer therapeutic approach to attenuate cancer cell progression. Int. J. Biol. Macromol. 2018, 107, 2580–2589. [Google Scholar] [CrossRef]
  67. Fatima, S.; Mohammad, T.; Jairajpuri, D.S.; Rehman, M.T.; Hussain, A.; Samim, M.; Ahmad, F.J.; Alajmi, M.F.; Hassan, M.I. Identification and evaluation of glutathione conjugate gamma-l-glutamyl-l-cysteine for improved drug delivery to the brain. J. Biomol. Struct. Dyn. 2020, 38, 3610–3620. [Google Scholar] [CrossRef] [PubMed]
  68. Lobanov, M.Y.; Bogatyreva, N.; Galzitskaya, O. Radius of gyration as an indicator of protein structure compactness. Mol. Biol. 2008, 42, 623–628. [Google Scholar] [CrossRef]
  69. Hubbard, R.E.; Haider, M.K. Hydrogen Bonds in Proteins: Role and Strength; eLS: Cologne, Germany, 2010. [Google Scholar]
  70. Dahiya, R.; Mohammad, T.; Roy, S.; Anwar, S.; Gupta, P.; Haque, A.; Khan, P.; Kazim, S.N.; Islam, A.; Ahmad, F. Investigation of inhibitory potential of quercetin to the pyruvate dehydrogenase kinase 3: Towards implications in anticancer therapy. Int. J. Biol. Macromol. 2019, 136, 1076–1085. [Google Scholar] [CrossRef]
  71. Maisuradze, G.G.; Liwo, A.; Scheraga, H.A. Principal component analysis for protein folding dynamics. J. Mol. Biol. 2009, 385, 312–329. [Google Scholar] [CrossRef] [Green Version]
  72. Altis, A.; Otten, M.; Nguyen, P.H.; Hegger, R.; Stock, G. Construction of the free energy landscape of biomolecules via dihedral angle principal component analysis. J. Chem. Phys. 2008, 128, 6B620. [Google Scholar] [CrossRef]
Jpm 12 01083 g001 550
Figure 1. The schematic representation of the structure of AKT1 in a bound state with tuberosin and villosol. (A) Cartoon illustration of AKT1 with the elucidated molecules. (B) The representation of close interaction between AKT1 with tuberosin (red salmon) and villosol (purple). (C) The magnified illustration shows the surface potential view of the binding pocket of AKT1 with elucidated compounds.
Figure 1. The schematic representation of the structure of AKT1 in a bound state with tuberosin and villosol. (A) Cartoon illustration of AKT1 with the elucidated molecules. (B) The representation of close interaction between AKT1 with tuberosin (red salmon) and villosol (purple). (C) The magnified illustration shows the surface potential view of the binding pocket of AKT1 with elucidated compounds.
Jpm 12 01083 g001
Jpm 12 01083 g002 550
Figure 2. The two-dimensional representation of binding pocket residues of AKT1 and its respective interactions with compounds tuberosin (A) and villosol (B). The tuberosin shows the noncovalent interactions formed between the active site and ATP-binding residues of protein AKT1. The hydrogen bonds are formed between sites Lys179, Glu191, and Asp 292. The villosol shows noncovalent interactions between crucial residues of protein Glu234 and Asp292, with target protein AKT1.
Figure 2. The two-dimensional representation of binding pocket residues of AKT1 and its respective interactions with compounds tuberosin (A) and villosol (B). The tuberosin shows the noncovalent interactions formed between the active site and ATP-binding residues of protein AKT1. The hydrogen bonds are formed between sites Lys179, Glu191, and Asp 292. The villosol shows noncovalent interactions between crucial residues of protein Glu234 and Asp292, with target protein AKT1.
Jpm 12 01083 g002
Jpm 12 01083 g003 550
Figure 3. The structural dynamics of AKT1 after binding with tuberosin and villosol based on the function of time. (A) RMSD graph showing AKT1 when bound with tuberosin and villosol. (B) The average residual fluctuations (RMSF) graph of AKT1 is associated with tuberosin and villosol. The plots shown in the lower panels depict the probability distribution as PDF.
Figure 3. The structural dynamics of AKT1 after binding with tuberosin and villosol based on the function of time. (A) RMSD graph showing AKT1 when bound with tuberosin and villosol. (B) The average residual fluctuations (RMSF) graph of AKT1 is associated with tuberosin and villosol. The plots shown in the lower panels depict the probability distribution as PDF.
Jpm 12 01083 g003
Jpm 12 01083 g004 550
Figure 4. The structural compactness of AKT1 when present with tuberosin and villosol binding as a function of time. (A) The radius of gyration (Rg) through time evaluation. (B) SASA plot of AKT1 as a function of time before and after binding with selected compounds. Lower graph plots depict the probability density function (PDF).
Figure 4. The structural compactness of AKT1 when present with tuberosin and villosol binding as a function of time. (A) The radius of gyration (Rg) through time evaluation. (B) SASA plot of AKT1 as a function of time before and after binding with selected compounds. Lower graph plots depict the probability density function (PDF).
Jpm 12 01083 g004
Jpm 12 01083 g005 550
Figure 5. Formation of hydrogen bonds in time evolution plot formed within 0.35 nm of intra-AKT1 (A). The probability density function (PDF) of the intramolecular hydrogen bonds present within the structure of AKT1 (B).
Figure 5. Formation of hydrogen bonds in time evolution plot formed within 0.35 nm of intra-AKT1 (A). The probability density function (PDF) of the intramolecular hydrogen bonds present within the structure of AKT1 (B).
Jpm 12 01083 g005
Jpm 12 01083 g006 550
Figure 6. The time evolution of intermolecular hydrogen bonds formed within 0.35 nm between AKT1 and tuberosin (A) and villosol (B). The panels on the lower side depict the probability density function distribution of H-bonds.
Figure 6. The time evolution of intermolecular hydrogen bonds formed within 0.35 nm between AKT1 and tuberosin (A) and villosol (B). The panels on the lower side depict the probability density function distribution of H-bonds.
Jpm 12 01083 g006
Jpm 12 01083 g007 550
Figure 7. The conformational sampling in principal compound analysis (PCA). The two-dimensional projections illustrate the conformational sampling of AKT1-tuberosin and AKT1-villosol complexes.
Figure 7. The conformational sampling in principal compound analysis (PCA). The two-dimensional projections illustrate the conformational sampling of AKT1-tuberosin and AKT1-villosol complexes.
Jpm 12 01083 g007
Jpm 12 01083 g008 550
Figure 8. The free energy landscapes (FELs) of (A) AKT1, (B) AKT1-tuberosin complex, and (C) AKT1-villosol complex.
Figure 8. The free energy landscapes (FELs) of (A) AKT1, (B) AKT1-tuberosin complex, and (C) AKT1-villosol complex.
Jpm 12 01083 g008
Table
Table 1. Top 30 hits against AKT1 and their docking scores.
Table 1. Top 30 hits against AKT1 and their docking scores.
S. No. Compound IDPhytochemical NameAffinity (kcal/mol)pKiLigand Efficiency
1.174880Lactupicrin−10.17.410.34
2.120698Capnoidine−10.07.330.37
3.308140Diospyrin−10.07.330.36
4.101277405Arboreol−10.07.330.36
5.104940Withanolide S−9.97.260.39
6.3822Cubebol−9.97.260.34
7.5315739N-Acetyldehydroanonaine−9.87.190.43
8.20055371-−9.77.110.40
9.14630495Tuberosin−9.77.110.39
10.146680(2,2′-Binaphthalene)-5,5′,8,8′-tetrone, 1,1′-dihydroxy-6,6′-dimethyl-−9.77.110.35
11.167455134-Hydroxysesamin−9.77.110.36
12.53172975-(6,7-Dihydroxy-2-oxochromen-5-yl)-6,7-dihydroxychromen-2-one−9.77.110.37
13.101746Sesamolin−9.67.040.35
14.12309624Papaverrubine B−9.67.040.34
15.196979Frangulin A−9.67.040.32
16.2913Acora-4,10-diene−9.56.970.43
17.73393Lagerine−9.56.970.36
18.21723446Pseudostrychnine−9.56.970.36
19.5316097Corylin−9.56.970.39
20.3295843,3′-Bisjuglone−9.56.970.36
21.5281806Psoralidin−9.56.970.38
22.119204Roemerine−9.46.890.45
23.5320486Gomezine−9.46.890.47
24.5459059Ribasine−9.46.890.36
25.91572Phaseolin−9.46.890.39
26.94577Cepharadione A−9.36.820.40
27.101595Vomicine−9.36.820.33
28.5320772Psoralenol−9.36.820.37
29.5490819Villosol−9.36.820.31
30.11438278Cryptodorine−9.26.750.40
Table
Table 2. ADMET properties of the selected compounds.
Table 2. ADMET properties of the selected compounds.
S. No.Compound IDPhytochemicalAbsorptionDistributionMetabolismExcretionToxicity
GI Absorption (%)BBB PermeationCYP2D6 InhibitorOCT2 SubstrateAMES and Hepatotoxicity
1.3295843,3′-BisjugloneHighNoNoNoNo
2.5320772PsoralenolHighNoYesNoNo
3.5459059RibasineHighYesYesNoNo
4.5490819VillosolHighNoYesNoNo
5.12309624Papaverrubine BHighYesYesNoNo
6.14630495TuberosinHighYesYesNoNo
Table
Table 3. PASS analysis and molecular structures of the selected compounds.
Table 3. PASS analysis and molecular structures of the selected compounds.
S. No.PhytochemicalMolecular StructurePaPiActivity
1.3,3′-Bisjuglone Jpm 12 01083 i00108660003Antimutagenic
08790016Membrane integrity agonist
08280009Antineoplastic
07990006Kinase inhibitor
07730007Caspase-3 stimulant
2.Psoralenol Jpm 12 01083 i00207040008Anticarcinogenic
06310023Apoptosis agonist
05350013Antimutagenic
05640053Antineoplastic
05060033Kinase inhibitor
3.Ribasine Jpm 12 01083 i00308780002Neurotransmitter uptake inhibitor
05120005Antineoplastic alkaloid
04740023Caspase-8 stimulant
03450079Caspase-3 stimulant
02760029Antineoplastic (lymphoma)
4.Villosol Jpm 12 01083 i00409930001TP53 expression enhancer
09020004Apoptosis agonist
08540007Antineoplastic
06620007Antineoplastic (breast cancer)
05560006Prostate cancer treatment
5.Papaverrubine B Jpm 12 01083 i00506830029Antineoplastic
05260028Antineoplastic (non-Hodgkin lymphoma)
05030006Antineoplastic alkaloid
04740009Antimitotic
03760027Antineoplastic (lung cancer)
6.Tuberosin Jpm 12 01083 i00607920013Antineoplastic
07700005Antineoplastic (breast cancer)
07270021TP53 expression enhancer
06270023Apoptosis agonist
05710005Antineoplastic (ovarian cancer)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Adnan, M.; Jairajpuri, D.S.; Chaddha, M.; Khan, M.S.; Yadav, D.K.; Mohammad, T.; Elasbali, A.M.; Abu Al-Soud, W.; Hussain Alharethi, S.; Hassan, M.I. Discovering Tuberosin and Villosol as Potent and Selective Inhibitors of AKT1 for Therapeutic Targeting of Oral Squamous Cell Carcinoma. J. Pers. Med. 2022, 12, 1083. https://doi.org/10.3390/jpm12071083

AMA Style

Adnan M, Jairajpuri DS, Chaddha M, Khan MS, Yadav DK, Mohammad T, Elasbali AM, Abu Al-Soud W, Hussain Alharethi S, Hassan MI. Discovering Tuberosin and Villosol as Potent and Selective Inhibitors of AKT1 for Therapeutic Targeting of Oral Squamous Cell Carcinoma. Journal of Personalized Medicine. 2022; 12(7):1083. https://doi.org/10.3390/jpm12071083

Chicago/Turabian Style

Adnan, Mohd, Deeba Shamim Jairajpuri, Muskan Chaddha, Mohd Shahnawaz Khan, Dharmendra Kumar Yadav, Taj Mohammad, Abdelbaset Mohamed Elasbali, Waleed Abu Al-Soud, Salem Hussain Alharethi, and Md. Imtaiyaz Hassan. 2022. "Discovering Tuberosin and Villosol as Potent and Selective Inhibitors of AKT1 for Therapeutic Targeting of Oral Squamous Cell Carcinoma" Journal of Personalized Medicine 12, no. 7: 1083. https://doi.org/10.3390/jpm12071083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop