1
|
Zhao Y, Tian Y, Pang X, Li G, Shi S, Yan A. Classification of FLT3 inhibitors and SAR analysis by machine learning methods. Mol Divers 2023:10.1007/s11030-023-10640-8. [PMID: 37142889 DOI: 10.1007/s11030-023-10640-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Accepted: 03/17/2023] [Indexed: 05/06/2023]
Abstract
FMS-like tyrosine kinase 3 (FLT3) is a type III receptor tyrosine kinase, which is an important target for anti-cancer therapy. In this work, we conducted a structure-activity relationship (SAR) study on 3867 FLT3 inhibitors we collected. MACCS fingerprints, ECFP4 fingerprints, and TT fingerprints were used to represent the inhibitors in the dataset. A total of 36 classification models were built based on support vector machine (SVM), random forest (RF), eXtreme Gradient Boosting (XGBoost), and deep neural networks (DNN) algorithms. Model 3D_3 built by deep neural networks (DNN) and TT fingerprints performed best on the test set with the highest prediction accuracy of 85.83% and Matthews correlation coefficient (MCC) of 0.72 and also performed well on the external test set. In addition, we clustered 3867 inhibitors into 11 subsets by the K-Means algorithm to figure out the structural characteristics of the reported FLT3 inhibitors. Finally, we analyzed the SAR of FLT3 inhibitors by RF algorithm based on ECFP4 fingerprints. The results showed that 2-aminopyrimidine, 1-ethylpiperidine,2,4-bis(methylamino)pyrimidine, amino-aromatic heterocycle, [(2E)-but-2-enyl]dimethylamine, but-2-enyl, and alkynyl were typical fragments among highly active inhibitors. Besides, three scaffolds in Subset_A (Subset 4), Subset_B, and Subset_C showed a significant relationship to inhibition activity targeting FLT3.
Collapse
Affiliation(s)
- Yunyang Zhao
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, Beijing University of Chemical Technology, 15 BeiSanHuan East Road, P.O. Box 53, Beijing, 100029, People's Republic of China
| | - Yujia Tian
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, Beijing University of Chemical Technology, 15 BeiSanHuan East Road, P.O. Box 53, Beijing, 100029, People's Republic of China
| | - Xiaoyang Pang
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, Beijing University of Chemical Technology, 15 BeiSanHuan East Road, P.O. Box 53, Beijing, 100029, People's Republic of China
| | - Guo Li
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, Beijing University of Chemical Technology, 15 BeiSanHuan East Road, P.O. Box 53, Beijing, 100029, People's Republic of China
| | - Shenghui Shi
- College of Information Science and Technology, Beijing University of Chemical Technology, 15 BeiSanHuan East Road, Beijing, 100029, People's Republic of China.
| | - Aixia Yan
- State Key Laboratory of Chemical Resource Engineering, Department of Pharmaceutical Engineering, Beijing University of Chemical Technology, 15 BeiSanHuan East Road, P.O. Box 53, Beijing, 100029, People's Republic of China.
| |
Collapse
|
2
|
Three-Dimensional-QSAR and Relative Binding Affinity Estimation of Focal Adhesion Kinase Inhibitors. Molecules 2023; 28:molecules28031464. [PMID: 36771129 PMCID: PMC9919860 DOI: 10.3390/molecules28031464] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 01/17/2023] [Accepted: 01/31/2023] [Indexed: 02/05/2023] Open
Abstract
Precise binding affinity predictions are essential for structure-based drug discovery (SBDD). Focal adhesion kinase (FAK) is a member of the tyrosine kinase protein family and is overexpressed in a variety of human malignancies. Inhibition of FAK using small molecules is a promising therapeutic option for several types of cancer. Here, we conducted computational modeling of FAK-targeting inhibitors using three-dimensional structure-activity relationship (3D-QSAR), molecular dynamics (MD), and hybrid topology-based free energy perturbation (FEP) methods. The structure-activity relationship (SAR) studies between the physicochemical descriptors and inhibitory activities of the chemical compounds were performed with reasonable statistical accuracy using CoMFA and CoMSIA. These are two well-known 3D-QSAR methods based on the principle of supervised machine learning (ML). Essential information regarding residue-specific binding interactions was determined using MD and MM-PB/GBSA methods. Finally, physics-based relative binding free energy (ΔΔGRBFEA→B) terms of analogous ligands were estimated using alchemical FEP simulation. An acceptable agreement was observed between the experimental and computed relative binding free energies. Overall, the results suggested that using ML and physics-based hybrid approaches could be useful in synergy for the rational optimization of accessible lead compounds with similar scaffolds targeting the FAK receptor.
Collapse
|
3
|
Ghosh S, Cho SJ. Comparative binding affinity analysis of dual
CDK2
/
FLT3
inhibitors. B KOREAN CHEM SOC 2022. [DOI: 10.1002/bkcs.12625] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Suparna Ghosh
- Department of Biomedical Sciences College of Medicine, Chosun University Gwangju Republic of Korea
| | - Seung Joo Cho
- Department of Biomedical Sciences College of Medicine, Chosun University Gwangju Republic of Korea
- Department of Cellular Molecular Medicine College of Medicine, Chosun University Gwangju Republic of Korea
| |
Collapse
|
4
|
Binding Studies and Lead Generation of Pteridin-7(8H)-one Derivatives Targeting FLT3. Int J Mol Sci 2022; 23:ijms23147696. [PMID: 35887060 PMCID: PMC9319409 DOI: 10.3390/ijms23147696] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2022] [Revised: 07/05/2022] [Accepted: 07/05/2022] [Indexed: 02/01/2023] Open
Abstract
Ligand modification by substituting chemical groups within the binding pocket is a popular strategy for kinase drug development. In this study, a series of pteridin-7(8H)-one derivatives targeting wild-type FMS-like tyrosine kinase-3 (FLT3) and its D835Y mutant (FL3D835Y) were studied using a combination of molecular modeling techniques, such as docking, molecular dynamics (MD), binding energy calculation, and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. We determined the protein–ligand binding affinity by employing molecular mechanics Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA), fast pulling ligand (FPL) simulation, linear interaction energy (LIE), umbrella sampling (US), and free energy perturbation (FEP) scoring functions. The structure–activity relationship (SAR) study was conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), and the results were emphasized as a SAR scheme. In both the CoMFA and CoMSIA models, satisfactory correlation statistics were obtained between the observed and predicted inhibitory activity. The MD and SAR models were co-utilized to design several new compounds, and their inhibitory activities were anticipated using the CoMSIA model. The designed compounds with higher predicted pIC50 values than the most active compound were carried out for binding free energy evaluation to wild-type and mutant receptors using MM-PB/GBSA, LIE, and FEP methods.
Collapse
|
5
|
Ghosh S, Cho SJ. Structure–activity
relationship and
in silico
development of
c‐Met
kinase inhibitors. B KOREAN CHEM SOC 2022. [DOI: 10.1002/bkcs.12537] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Suparna Ghosh
- Department of Biomedical Sciences College of Medicine, Chosun University Gwangju Republic of Korea
| | - Seung Joo Cho
- Department of Biomedical Sciences College of Medicine, Chosun University Gwangju Republic of Korea
- Department of Cellular Molecular Medicine College of Medicine, Chosun University Gwangju Republic of Korea
| |
Collapse
|
6
|
Structural Insights from Molecular Modeling of Isoindolin-1-One Derivatives as PI3Kγ Inhibitors against Gastric Carcinoma. Biomedicines 2022; 10:biomedicines10040813. [PMID: 35453562 PMCID: PMC9030798 DOI: 10.3390/biomedicines10040813] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Revised: 03/24/2022] [Accepted: 03/29/2022] [Indexed: 01/15/2023] Open
Abstract
The upregulation of phosphoinositol-3-kinase γ (PI3Kγ) is deemed to be positively correlated with tumor-associated-macrophage (TAM)-mediated gastric carcinoma (GC). PI3Kγ suppresses tumor necrosis factor-alpha (TNF-α) and interleukin-12 (IL-12) through activation of the AKT/mTOR pathway, which promotes the immunosuppressant phenotype of TAM. Unlike α and β isoforms, δ and γ isoforms are primarily distributed in leucocytes and macrophages. Dual inhibitors against PI3Kδ and PI3Kγ have been proven to have merits in targeting solid tumors. Furthermore, it has been found that PI3Kδ is activated by cytokines, while PI3Kγ is activated by G-protein-coupled receptors (GPCRs). This facilitates determining the functional difference between these two isoforms. For this goal, selective inhibitors would be immensely helpful. In the current manuscript, we conducted various molecular modeling studies with a series of isoindolin-1-one derivatives as potent PI3Kγ inhibitors by combining molecular docking, molecular dynamics (MD), molecular mechanics, Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA) binding free energy calculation, and three-dimensional structure–activity relationship (3D-QSAR) study. To evaluate the selectivity of γ isoform over δ, the molecular modeling studies of idelalisib analogs reported as PI3Kδ inhibitors were also investigated. The contour polyhedrons were generated from the comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) around the ligand-bound active site for both isoforms, which could emphasize plausible explanations for the physicochemical factors that affect selective ligand recognition. The binding modalities of the two isoforms using CoMFA and MD models were compared, which suggested some key differences in the molecular interactions with the ligands and could be summarized as three subsites (one affinity subsite near the C-helix and DFG and two hydrophobic subsites). In the context of the structure–activity relationship (SAR), several new compounds were designed using a fragment-substitution strategy with the aim of selectively targeting PI3Kγ. The pIC50 values of the designed compounds were predicted by the 3D-QSAR models, followed by the MM-PB/GBSA binding energy estimation. The overall findings suggest that the designed compounds have the potential to be used as PI3Kγ inhibitors with a higher binding affinity and selectivity.
Collapse
|