1
|
Piña MDLN, Sahu AK, Frontera A, Biswal HS, Bauzá A. Tetrel bonds involving a CF 3 group participate in protein-drug recognition: a combined crystallographic and computational study. Phys Chem Chem Phys 2023; 25:12409-12419. [PMID: 37093130 DOI: 10.1039/d3cp00839h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
In this study, the ability of CF3 groups to bind to the electron-rich side chains and backbone groups of proteins has been investigated by combining a Protein Data Bank (PDB) survey and ab initio quantum mechanics calculations. More precisely, an inspection of the PDB involving organic ligands containing a CF3 group and electron-rich atoms (A = N, O and S) in the vicinity revealed 419 X-ray structures exhibiting CF3⋯A tetrel bonds (TtBs). In a posterior stage, those hits that exhibited the most relevant features in terms of directionality and intermolecular distance were selected for theoretical calculations at the RI-MP2/def2-TZVPD level of theory. Also, Hammett's regression plots of several TtB complexes involving meta- and para-substituted benzene derivatives were computed to shed light on the substituent effects. Moreover, the TtBs were characterized through several state-of-the-art computational techniques, such as the Quantum Theory of Atoms in Molecules (QTAIM) and Noncovalent Interactions plot (NCIplot) methodologies. We believe that the results gathered from our study will be useful for rational drug design and biological communities as well as for further expanding the role of this interaction to biomedical applications.
Collapse
Affiliation(s)
- María de Las Nieves Piña
- Department of Chemistry, Universitat de les Illes Balears, Crts de Valldemossa km 7.6, 07122 Palma de Mallorca, Baleares, Spain.
| | - Akshay Kumar Sahu
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), PO-Bhimpur-Padanpur, Via-Jatni, Khurda, 752050, Bhubaneswar, India.
- Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Antonio Frontera
- Department of Chemistry, Universitat de les Illes Balears, Crts de Valldemossa km 7.6, 07122 Palma de Mallorca, Baleares, Spain.
| | - Himansu S Biswal
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), PO-Bhimpur-Padanpur, Via-Jatni, Khurda, 752050, Bhubaneswar, India.
- Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Antonio Bauzá
- Department of Chemistry, Universitat de les Illes Balears, Crts de Valldemossa km 7.6, 07122 Palma de Mallorca, Baleares, Spain.
| |
Collapse
|
2
|
Li Y, Liu S, Wang J, Rui X, Tian H, Li C, Guo C. In Silico Studies of Piperidine Derivatives as Protein Kinase B Inhibitors
through 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation. LETT DRUG DES DISCOV 2022. [DOI: 10.2174/1570180818666211207105516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Background:
Protein kinase B (Akt) is a serine/threonine-protein kinase that drives the diverse
physiological process. Akt is a promising therapeutic target, which involves cancer cell growth, survival,
proliferation and metabolism.
Objective:
The study aims to design highly active Akt inhibitors, and to elucidate the structural requirements
for their biological activity, we analyzed the key binding features and summarized the structural
determinants for their bioactivities.
Methods:
A series of piperidine derivatives have been investigated employing three-dimensional quantitative
structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics simulation.
Results:
The statistics of the comparative molecular field analysis (CoMFA) model (Q2=0.631, R2=0.951)
and the comparative molecular similarity index analysis (CoMSIA) model (Q2=0.663, R2=0.966) indicated
that our 3D-QSAR model was accurate and reliable. Besides, the stability of receptor-ligand interactions
under physiological conditions was then evaluated by molecular dynamics simulation, in agreement
with the molecular docking results.
Conclusion:
Our study provided valuable insights for the discovery of potent Akt inhibitors.
Collapse
Affiliation(s)
- Yaxin Li
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
- Hebei Key Laboratory of Neuropharmacology,
Hebei North University, Zhangjiakou 075000, P.R. China
| | - Sisi Liu
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
- Hebei Key Laboratory of Neuropharmacology,
Hebei North University, Zhangjiakou 075000, P.R. China
| | - Jin Wang
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
- Hebei Key Laboratory of Neuropharmacology,
Hebei North University, Zhangjiakou 075000, P.R. China
| | - Xue Rui
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
| | - Haobo Tian
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
| | - Chenshuo Li
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
| | - Chunyan Guo
- Department of Pharmacy, Hebei North University, Zhangjiakou 075000, P.R. China
- Hebei Key Laboratory of Neuropharmacology,
Hebei North University, Zhangjiakou 075000, P.R. China
| |
Collapse
|
3
|
Keretsu S, Ghosh S, Cho SJ. Computer aided designing of novel pyrrolopyridine derivatives as JAK1 inhibitors. Sci Rep 2021; 11:23051. [PMID: 34845259 PMCID: PMC8630053 DOI: 10.1038/s41598-021-02364-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2021] [Accepted: 11/10/2021] [Indexed: 02/07/2023] Open
Abstract
Janus kinases (JAKs) are a family of non-receptor kinases that play a key role in cytokine signaling and their aberrant activities are associated with the pathogenesis of various immune diseases. The JAK1 isoform plays an essential role in the types 1 and II interferon signaling and elicits signals from the interleukin-2, interleukin-4, gp130, and class 2 receptor families. It is ubiquitously expressed in humans and its overexpression has been linked with autoimmune diseases such as myeloproliferative neoplasm. Although JAK1 inhibitors such as Tofacitinib have been approved for medical use, the low potency and off-target effects of these inhibitors have limited their use and calls for the development of novel JAK1 inhibitors. In this study, we used computational methods on a series of pyrrolopyridine derivatives to design new JAK1 inhibitors. Molecular docking and molecular dynamics simulation methods were used to study the protein-inhibitor interactions. 3D-quantitative structure–activity relationship models were developed and were used to predict the activity of newly designed compounds. Free energy calculation methods were used to study the binding affinity of the inhibitors with JAK1. Of the designed compounds, seventeen of the compounds showed a higher binding energy value than the most active compound in the dataset and at least six of the compounds showed higher binding energy value than the pan JAK inhibitor Tofacitinib. The findings made in this study could be utilized for the further development of JAK1 inhibitors.
Collapse
Affiliation(s)
- Seketoulie Keretsu
- Department of Biomedical Sciences, College of Medicine, Chosun University, 375 Seosuk-dong, Dong-gu, Gwangju, 501-759, Republic of Korea
| | - Suparna Ghosh
- Department of Biomedical Sciences, College of Medicine, Chosun University, 375 Seosuk-dong, Dong-gu, Gwangju, 501-759, Republic of Korea
| | - Seung Joo Cho
- Department of Biomedical Sciences, College of Medicine, Chosun University, 375 Seosuk-dong, Dong-gu, Gwangju, 501-759, Republic of Korea. .,Department of Cellular Molecular Medicine, College of Medicine, Chosun University, Gwangju, 501-759, Republic of Korea.
| |
Collapse
|
4
|
Ghosh S, Keretsu S, Cho SJ. Designing of the N-ethyl-4-(pyridin-4-yl)benzamide based potent ROCK1 inhibitors using docking, molecular dynamics, and 3D-QSAR. PeerJ 2021; 9:e11951. [PMID: 34434664 PMCID: PMC8359802 DOI: 10.7717/peerj.11951] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Accepted: 07/20/2021] [Indexed: 01/20/2023] Open
Abstract
Rho-associated kinase-1 (ROCK1) has been recognized for its pivotal role in heart diseases, different types of malignancy, and many neurological disorders. Hyperactivity of ROCK phosphorylates the protein kinase-C (PKC), which ultimately induces smooth muscle cell contraction in the vascular system. Inhibition of ROCK1 has been shown to be a promising therapy for patients with cardiovascular disease. In this study, we have conducted molecular modeling techniques such as docking, molecular dynamics (MD), and 3-Dimensional structure-activity relationship (3D-QSAR) on a series of N-ethyl-4-(pyridin-4-yl)benzamide-based compounds. Docking and MD showed critical interactions and binding affinities between ROCK1 and its inhibitors. To establish the structure-activity relationship (SAR) of the compounds, 3D-QSAR techniques such as Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) were used. The CoMFA (q 2 = 0.774, r 2 = 0.965, ONC = 6, and r p r e d 2 = 0.703) and CoMSIA (q 2 = 0.676, r 2 = 0.949, ONC = 6, and r p r e d 2 = 0.548) both models have shown reasonable external predictive activity, and contour maps revealed favorable and unfavorable substitutions for chemical group modifications. Based on the contour maps, we have designed forty new compounds, among which, seven compounds exhibited higher predictive activity (pIC50). Further, we conducted the MD study, ADME/Tox, and SA score prediction using the seven newly designed compounds. The combination of docking, MD, and 3D-QSAR studies helps to understand the coherence modification of existing molecules. Our study may provide valuable insight into the development of more potent ROCK1 inhibitors.
Collapse
Affiliation(s)
- Suparna Ghosh
- Department of Biomedical Sciences, College of Medicine, Chosun University, Gwangju, South Korea
| | - Seketoulie Keretsu
- Department of Biomedical Sciences, College of Medicine, Chosun University, Gwangju, South Korea
| | - Seung Joo Cho
- Department of Biomedical Sciences, College of Medicine, Chosun University, Gwangju, South Korea.,Department of Cellular and Molecular Medicine, College of Medicine, Chosun University, Gwangju, South Korea
| |
Collapse
|
5
|
Ghosh S, Keretsu S, Cho SJ. Computational Modeling of Novel Phosphoinositol‐3‐kinase γ Inhibitors Using Molecular Docking, Molecular Dynamics, and
3D‐QSAR. B KOREAN CHEM SOC 2021. [DOI: 10.1002/bkcs.12305] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Suparna Ghosh
- Department of Biomedical Sciences, College of Medicine Chosun University Gwangju 501‐759 Republic of Korea
| | - Seketoulie Keretsu
- Department of Biomedical Sciences, College of Medicine Chosun University Gwangju 501‐759 Republic of Korea
| | - Seung Joo Cho
- Department of Biomedical Sciences, College of Medicine Chosun University Gwangju 501‐759 Republic of Korea
- Department of Cellular and Molecular Medicine, College of Medicine Chosun University Gwangju 501‐759 Republic of Korea
| |
Collapse
|
6
|
Keretsu S, Ghosh S, Cho SJ. Molecular Modeling Study of c-KIT/PDGFRα Dual Inhibitors for the Treatment of Gastrointestinal Stromal Tumors. Int J Mol Sci 2020; 21:ijms21218232. [PMID: 33153146 PMCID: PMC7662224 DOI: 10.3390/ijms21218232] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Revised: 10/26/2020] [Accepted: 10/30/2020] [Indexed: 02/06/2023] Open
Abstract
Gastrointestinal stromal tumors (GISTs) are the most common Mesenchymal Neoplasm of the gastrointestinal tract. The tumorigenesis of GISTs has been associated with the gain-of-function mutation and abnormal activation of the stem cell factor receptor (c-KIT) and platelet-derived growth factor receptor alpha (PDGFRα) kinases. Hence, inhibitors that target c-KIT and PDGFRα could be a therapeutic option for the treatment of GISTs. The available approved c-KIT/PDGFRα inhibitors possessed low efficacy with off-target effects, which necessitated the development of potent inhibitors. We performed computational studies of 48 pyrazolopyridine derivatives that showed inhibitory activity against c-KIT and PDGFRα to study the structural properties important for inhibition of both the kinases. The derivative of phenylurea, which has high activities for both c-KIT (pIC50 = 8.6) and PDGFRα (pIC50 = 8.1), was used as the representative compound for the dataset. Molecular docking and molecular dynamics simulation (100 ns) of compound 14 was performed. Compound 14 showed the formation of hydrogen bonding with Cys673, Glu640, and Asp810 in c-KIT, and Cys677, Glu644, and Asp836 in PDGFRα. The results also suggested that Thr670/T674 substitution in c-KIT/PDGFRα induced conformational changes at the binding site of the receptors. Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed based on the inhibitors. Contour map analysis showed that electropositive and bulky substituents at the para-position and the meta-position of the benzyl ring of compound 14 was favorable and may increase the inhibitory activity against both c-KIT and PDGFRα. Analysis of the results suggested that having bulky and hydrophobic substituents that extend into the hydrophobic pocket of the binding site increases the activity for both c-KIT and PDGFRα. Based on the contour map analysis, 50 compounds were designed, and the activities were predicted. An evaluation of binding free energy showed that eight of the designed compounds have potential binding affinity with c-KIT/PDGFRα. Absorption, distribution, metabolism, excretion and toxicity (ADMET) and synthetic feasibility tests showed that the designed compounds have reasonable pharmaceutical properties and synthetic feasibility. Further experimental study of the designed compounds is recommended. The structural information from this study could provide useful insight into the future development of c-KIT and PDGFRα inhibitors.
Collapse
Affiliation(s)
- Seketoulie Keretsu
- Department of Biomedical Sciences, College of Medicine, Chosun University, Gwangju 501-759, Korea; (S.K.); (S.G.)
| | - Suparna Ghosh
- Department of Biomedical Sciences, College of Medicine, Chosun University, Gwangju 501-759, Korea; (S.K.); (S.G.)
| | - Seung Joo Cho
- Department of Biomedical Sciences, College of Medicine, Chosun University, Gwangju 501-759, Korea; (S.K.); (S.G.)
- Department of Cellular Molecular Medicine, College of Medicine, Chosun University, Gwangju 501-759, Korea
- Correspondence: ; Tel.: +82-62-230-7482 or +82-11-479-1010
| |
Collapse
|