Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Aug 3;23(15):8635.
doi: 10.3390/ijms23158635.

Identification of Novel Inhibitors Targeting SGK1 via Ensemble-Based Virtual Screening Method, Biological Evaluation and Molecular Dynamics Simulation

Affiliations

Identification of Novel Inhibitors Targeting SGK1 via Ensemble-Based Virtual Screening Method, Biological Evaluation and Molecular Dynamics Simulation

Hui Zhang et al. Int J Mol Sci. .

Abstract

Serum and glucocorticoid-regulated kinase 1 (SGK1), as a serine threonine protein kinase of the AGC family, regulates different enzymes, transcription factors, ion channels, transporters, and cell proliferation and apoptosis. Inhibition of SGK1 is considered as a valuable approach for the treatment of various metabolic diseases. In this investigation, virtual screening methods, including pharmacophore models, Bayesian classifiers, and molecular docking, were combined to discover novel inhibitors of SGK1 from the database with 29,158 compounds. Then, the screened compounds were subjected to ADME/T, PAINS and drug-likeness analysis. Finally, 28 compounds with potential inhibition activity against SGK1 were selected for biological evaluation. The kinase inhibition activity test revealed that among these 28 hits, hit15 exhibited the highest inhibition activity against SGK1, which gave 44.79% inhibition rate at the concentration of 10 µM. In order to further investigate the interaction mechanism of hit15 and SGK1 at simulated physiological conditions, a molecular dynamics simulation was performed. The molecular dynamics simulation demonstrated that hit15 could bind to the active site of SGK1 and form stable interactions with key residues, such as Tyr178, ILE179, and VAL112. The binding free energy of the SGK1-hit15 was -48.90 kJ mol-1. Therefore, the identified hit15 with novel scaffold may be a promising lead compound for development of new SGK1 inhibitors for various diseases treatment.

Keywords: SGK1 inhibitor; biological evaluation; molecular dynamics; virtual screening.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

Figure 1
Figure 1
The ensemble-based virtual screening strategy was developed for discovering potent inhibitors targeting for SGK1.
Figure 2
Figure 2
Representative protein conformers based on 3HDM (A) and 3HDN (B) were extracted from the MD trajectory by using structural clustering.
Figure 3
Figure 3
The best pharmacophore models of SGK1 inhibitors. (A) The 3D spatial relationship and geometric parameters of Hypo1. (B) The Hypo1 aligned with cpd1 (IC50 = 0.04 μM). (C) The Hypo1 aligned with lower activity of cpd6 (IC50 = 6.78 μM). (D) The 3D spatial relationship and geometric parameters of Hypo2. (E) The Hypo2 aligned with GMG. (F) The Hypo2 mapped onto the compound with lower activity (IC50 = 2.70 μM).
Figure 4
Figure 4
Distributions of the docking scores of the inhibitors/noninhibitors for each SGK1 protein.
Figure 5
Figure 5
The ROC curve of the naive Bayesian classifier. (A) Training set; (B) test set.
Figure 6
Figure 6
Dynamics of hit15, GMG, and STSP binding to the SGK1 active site. (A) Comparison of orientations of protein backbone for three SGK1-ligand systems. (B) Comparison of orientations of ligands for the three SGK1-ligand systems. (C) Time evolution of radius of gyration (Rg) values during 100 ns MD simulation. (D) Cα atomic fluctuations root-mean-square fluctuation (RMSF) plots. Blue, purple, and red represent the SGK1-hit15, SGK1-GMG, and SGK1-STSP, respectively. Root-mean-square deviation (RMSD) plots as a function of time. (EG) represented the structural superposition before (blue) and after (green) MD for SGK1-hit15 (E), SGK1-GMG (F), and SGK1-STSP (G). A porcupine plot was generated and the movements of flexible parts of the SGK1 indicated with red arrows, and cone represents the direction.
Figure 7
Figure 7
Molecular dynamics results of SS and SASA during 100 ns MD simulation. (AC) represented the SS changes during MD simulation for SGK1-hit15 (A), SGK1-GMG (B), and SGK1-STSP (C). (DF) exhibited the total, hydrophobic, and hydrophilic area changes during MD simulation for SGK1-hit15 (D), SGK1-GMG (E), and SGK1-STSP (F).
Figure 8
Figure 8
Number of hydrogen bond interactions formed during MD simulation in the case of SGK1-hit15 (A), SGK1-GMG (B), and SGK1-STSP (C).
Figure 9
Figure 9
Molecular dynamics results of tertiary contact map for SGK1-STSP (A), SGK1-GMG (B), and SGK1-hit15 (C).
Figure 10
Figure 10
Molecular dynamics results of PCA analysis from the last 5 ns trajectories. Projection of the motion of SGK1-hit15, SGK1-GMG, and SGK1-STSP in phase space along the PC1 and PC2.
Figure 11
Figure 11
Molecular dynamics results of clustering and binding mode analysis. (A) Conformation clustering from the last 3 ns of SGK1-hit15. The Cluster_1 was colored by blue. (B) Conformation clustering from the last 3 ns of SGK1-GMG. The Cluster_1 was colored by purple. (C) Conformation clustering from the last 3 ns of SGK1-STSP. The Cluster_1 was colored by red. (DF) represents the binding modes of SGK1-hit15 (D), SGK1-GMG (E), and SGK1-STSP (F). Hydrogen bonds are depicted as green dashed lines, Alkyl and Pi-Alkyl are represented as pink dashed lines, Pi-Pi Stacked are depicted as purple-pink dashed lines, Pi-Cation are exhibited as orange dashed lines.
Figure 12
Figure 12
Binding free energy calculated and decomposed by MM/PBSA method. (A) The binding free energies of SGK1-hit15, SGK1-GMG, and SGK1-STSP. (B) The component energies of binding free energy for SGK1-hit15, SGK1-GMG, and SGK1-STSP. (CE) indicate per-residue binding free-energy decomposition for SGK1-hit15 (C), SGK1-GMG (D), and SGK1-STSP (E). (F) The radar map of 6 important residues binding free-energy decomposition for SGK1-hit15, SGK1-GMG, and SGK1-STSP.
Figure 13
Figure 13
The chemical structures of SGK1 inhibitors together with their biological activity data (IC50) for HipHop running.

Similar articles

Cited by

References

    1. Basnet R., Gong G.Q., Li C., Wang M.W. Serum and glucocorticoid inducible protein kinases (SGKs): A potential target for cancer intervention. Acta Pharm. Sin. B. 2018;8:767–771. doi: 10.1016/j.apsb.2018.07.001. - DOI - PMC - PubMed
    1. Vanhaelen Q., Mamoshina P., Aliper A.M., Artemov A., Lezhnina K., Ozerov I., Labat I., Zhavoronkov A. Design of efficient computational workflows for in silico drug repurposing. Drug Discov. Today. 2017;22:210–222. doi: 10.1016/j.drudis.2016.09.019. - DOI - PubMed
    1. Whitlock G., Lewington S., Sherliker P. Body-mass index and cause-specific mortality in 900,000 adults: Collaborative analyses of 57 prospective studies. Lancet. 2009;373:1083–1096. - PMC - PubMed
    1. Lang F., Böhmer C., Palmada M., Seebohm G., Strutz-Seebohm N., Vallon V. (Patho)physiological significance of the serum-and glucocorticoid-inducible kinase isoforms. Physiol. Rev. 2006;86:1151–1178. doi: 10.1152/physrev.00050.2005. - DOI - PubMed
    1. Kobayashi T., Cohen P. Activation of serum- and glucocorticoid-regulated protein kinase by agonists that activate phosphatidylinositide 3-kinase is mediated by 3-phosphoinositide-dependent protein kinase-1 (PDK1) and PDK2. Biochem. J. 1999;339:319–328. doi: 10.1042/bj3390319. - DOI - PMC - PubMed

Grants and funding

This work was supported by the National Natural Science Foundation of China (Grant No. 81903543 and 81660589) and Science and Technology Program Project of Gansu Province (20JR5RA534).

LinkOut - more resources