{"title":"基于配体的化学信息学和自由能量启发的分子模拟,优先选择和优化多发性骨髓瘤治疗中g蛋白偶联受体激酶6 (GRK6)抑制剂。","authors":"Arnab Bhattacharjee, Supratik Kar, Probir Kumar Ojha","doi":"10.1016/j.compbiolchem.2025.108347","DOIUrl":null,"url":null,"abstract":"<p><p>Multiple myeloma (MM) is the second most frequently diagnosed hematological malignancy, presenting limited treatment options with no curative potential and significant drug resistance. Recent studies involving genetic knockdown established the crucial role of GRK6 in upholding the viability of MM cells, emphasizing the need to identify potential inhibitors. Computational exploration of GRK6 inhibitors has not been attempted previously. Herein, the present study reports a multilayered lead prioritization and optimization framework using chemometrics and molecular simulations. 2D QSAR studies revealed that hydrogen bonding and polar interactions enhanced GRK6 inhibitory activity, while increased electron accessibility posed a risk of off-target effects. The pharmacophore hypothesis (DDHRRR_1) featured two hydrogen bond donors, one hydrophobic region, and three aromatic rings, laying the foundation for the 3D QSAR models. Hydrophobic groups, such as pyridine and pyrazole, were shown to enhance inhibition, while smaller groups, like ethyl and hydroxyl, reduced activity. 12,557 DrugBank compounds were screened using the developed chemometric models and molecular docking in tandem, which led to the identification of 7 potential parent leads for subsequent QSAR-guided structural optimizations. 350 lead analogs were generated and the top 4 were further analyzed using molecular docking, ADMET, molecular dynamics, and metadynamics analysis based on Principal Component Analysis (PCA), Probability Density Function (PDF), and Free Energy Landscapes (FEL). Upon cumulative retrospection, we propose a novel analog of DB07168 (DB07168-A13) (docking score: -11.2 kcal/mol, MM-GBSA binding energy: -55.2 kcal/mol) as the most promising GRK6 inhibitor, warranting further in vitro validation, for addressing prospective therapeutic intervention in MM.</p>","PeriodicalId":93952,"journal":{"name":"Computational biology and chemistry","volume":"115 ","pages":"108347"},"PeriodicalIF":0.0000,"publicationDate":"2025-01-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Ligand-based cheminformatics and free energy-inspired molecular simulations for prioritizing and optimizing G-protein coupled receptor kinase-6 (GRK6) inhibitors in multiple myeloma treatment.\",\"authors\":\"Arnab Bhattacharjee, Supratik Kar, Probir Kumar Ojha\",\"doi\":\"10.1016/j.compbiolchem.2025.108347\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<p><p>Multiple myeloma (MM) is the second most frequently diagnosed hematological malignancy, presenting limited treatment options with no curative potential and significant drug resistance. Recent studies involving genetic knockdown established the crucial role of GRK6 in upholding the viability of MM cells, emphasizing the need to identify potential inhibitors. Computational exploration of GRK6 inhibitors has not been attempted previously. Herein, the present study reports a multilayered lead prioritization and optimization framework using chemometrics and molecular simulations. 2D QSAR studies revealed that hydrogen bonding and polar interactions enhanced GRK6 inhibitory activity, while increased electron accessibility posed a risk of off-target effects. The pharmacophore hypothesis (DDHRRR_1) featured two hydrogen bond donors, one hydrophobic region, and three aromatic rings, laying the foundation for the 3D QSAR models. Hydrophobic groups, such as pyridine and pyrazole, were shown to enhance inhibition, while smaller groups, like ethyl and hydroxyl, reduced activity. 12,557 DrugBank compounds were screened using the developed chemometric models and molecular docking in tandem, which led to the identification of 7 potential parent leads for subsequent QSAR-guided structural optimizations. 350 lead analogs were generated and the top 4 were further analyzed using molecular docking, ADMET, molecular dynamics, and metadynamics analysis based on Principal Component Analysis (PCA), Probability Density Function (PDF), and Free Energy Landscapes (FEL). Upon cumulative retrospection, we propose a novel analog of DB07168 (DB07168-A13) (docking score: -11.2 kcal/mol, MM-GBSA binding energy: -55.2 kcal/mol) as the most promising GRK6 inhibitor, warranting further in vitro validation, for addressing prospective therapeutic intervention in MM.</p>\",\"PeriodicalId\":93952,\"journal\":{\"name\":\"Computational biology and chemistry\",\"volume\":\"115 \",\"pages\":\"108347\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2025-01-13\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Computational biology and chemistry\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/10.1016/j.compbiolchem.2025.108347\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computational biology and chemistry","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1016/j.compbiolchem.2025.108347","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Ligand-based cheminformatics and free energy-inspired molecular simulations for prioritizing and optimizing G-protein coupled receptor kinase-6 (GRK6) inhibitors in multiple myeloma treatment.
Multiple myeloma (MM) is the second most frequently diagnosed hematological malignancy, presenting limited treatment options with no curative potential and significant drug resistance. Recent studies involving genetic knockdown established the crucial role of GRK6 in upholding the viability of MM cells, emphasizing the need to identify potential inhibitors. Computational exploration of GRK6 inhibitors has not been attempted previously. Herein, the present study reports a multilayered lead prioritization and optimization framework using chemometrics and molecular simulations. 2D QSAR studies revealed that hydrogen bonding and polar interactions enhanced GRK6 inhibitory activity, while increased electron accessibility posed a risk of off-target effects. The pharmacophore hypothesis (DDHRRR_1) featured two hydrogen bond donors, one hydrophobic region, and three aromatic rings, laying the foundation for the 3D QSAR models. Hydrophobic groups, such as pyridine and pyrazole, were shown to enhance inhibition, while smaller groups, like ethyl and hydroxyl, reduced activity. 12,557 DrugBank compounds were screened using the developed chemometric models and molecular docking in tandem, which led to the identification of 7 potential parent leads for subsequent QSAR-guided structural optimizations. 350 lead analogs were generated and the top 4 were further analyzed using molecular docking, ADMET, molecular dynamics, and metadynamics analysis based on Principal Component Analysis (PCA), Probability Density Function (PDF), and Free Energy Landscapes (FEL). Upon cumulative retrospection, we propose a novel analog of DB07168 (DB07168-A13) (docking score: -11.2 kcal/mol, MM-GBSA binding energy: -55.2 kcal/mol) as the most promising GRK6 inhibitor, warranting further in vitro validation, for addressing prospective therapeutic intervention in MM.