1
|
Bernetti M, Bosio S, Bresciani V, Falchi F, Masetti M. Probing allosteric communication with combined molecular dynamics simulations and network analysis. Curr Opin Struct Biol 2024; 86:102820. [PMID: 38688074 DOI: 10.1016/j.sbi.2024.102820] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Revised: 03/28/2024] [Accepted: 04/03/2024] [Indexed: 05/02/2024]
Abstract
Understanding the allosteric mechanisms within biomolecules involved in diseases is of paramount importance for drug discovery. Indeed, characterizing communication pathways and critical hotspots in signal transduction can guide a rational approach to leverage allosteric modulation for therapeutic purposes. While the atomistic signatures of allosteric processes are difficult to determine experimentally, computational methods can be a remarkable resource. Network analysis built on Molecular Dynamics simulation data is particularly suited in this respect and is gradually becoming of routine use. Herein, we collect the recent literature in the field, discussing different aspects and available options for network construction and analysis. We further highlight interesting refinements and extensions, eventually providing our perspective on this topic.
Collapse
Affiliation(s)
- Mattia Bernetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum - University of Bologna, Via Belmeloro 6, 40126 Bologna, Italy; Computational and Chemical Biology, Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy.
| | - Stefano Bosio
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum - University of Bologna, Via Belmeloro 6, 40126 Bologna, Italy; Computational and Chemical Biology, Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy. https://twitter.com/Stefano__Bosio
| | - Veronica Bresciani
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum - University of Bologna, Via Belmeloro 6, 40126 Bologna, Italy; Computational and Chemical Biology, Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy. https://twitter.com/V_Bresciani
| | - Federico Falchi
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum - University of Bologna, Via Belmeloro 6, 40126 Bologna, Italy; Computational and Chemical Biology, Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum - University of Bologna, Via Belmeloro 6, 40126 Bologna, Italy.
| |
Collapse
|
2
|
Qi X, Li G, Liu J, Mou L, Zhang Y, Guo S, Chen X, Li W. Structural and energetic insights into the selective inhibition of PKMYT1 against WEE1. J Biomol Struct Dyn 2024; 42:3010-3018. [PMID: 37345529 DOI: 10.1080/07391102.2023.2225106] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Accepted: 04/30/2023] [Indexed: 06/23/2023]
Abstract
Protein kinase, membrane-associated tyrosine/threonine 1 (PKMYT1), a member of the WEE family and responsible for the regulation of CDK1 phosphorylation, has been considered a promising therapeutic target for cancer therapy. However, the highly structural conservation of the ATP-binding sites of the WEE family poses a challenge to the design of selective inhibitors for PKMYT1. Here, molecular docking, multiple microsecond-length molecular dynamics (MD) simulations and end-point free energy calculations were performed to uncover the molecular mechanism of the binding selectivity of RP-6306 toward PKMYT1 over its highly homologous kinase WEE1. The binding specificity of RP-6306 reported in previous experimental bioassays was clarified by MD simulations and binding free energy calculations. Further, the binding free energy prediction indicated that the binding selectivity of RP-6306 largely derived from the difference in the protein-ligand electrostatic interactions. The per-residue free energy decomposition suggested that the non-conserved gatekeeper residue in the hinge domain of PKMYT1/WEE1, Thr187/Asn376, is the critical factor responsible for the binding selectivity of RP-6306 toward PKMYT1. In addition, a water-mediated hydrogen bond was formed between RP-6306 and Gly191 at the hinge domain in the PKMYT1/RP-6306 complex, which was absent in the WEE1/RP-6306 complex. This study is expected to offer useful information for the design of more potent and selective PKMYT1 inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Xuesen Qi
- Department of Urology, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| | - Guozhen Li
- Department of Urology, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| | - Jiahai Liu
- Department of Gastrointestinal and Anal Diseases Surgery, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| | - Linkai Mou
- Department of Urology, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| | - Yusheng Zhang
- Department of Urology, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| | - Shilin Guo
- Department of Urology, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| | - Xiangyu Chen
- School of Medical Laboratory, Weifang Medical University, Weifang, Shandong, China
| | - Wenxing Li
- Pathology Department, The Affiliated Hospital of Weifang Medical University, Weifang Medical University, Weifang, Shandong, China
| |
Collapse
|
3
|
Li X, Hou C, Yang M, Luo B, Mao N, Chen K, Chen Z, Bai Y. The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin. J Biomol Struct Dyn 2024:1-9. [PMID: 38457488 DOI: 10.1080/07391102.2024.2316763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Accepted: 02/05/2024] [Indexed: 03/10/2024]
Abstract
Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Xiaolong Li
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Canglong Hou
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Mingyuan Yang
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Beier Luo
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Ningfang Mao
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Kai Chen
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Ziqiang Chen
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| | - Yushu Bai
- Department of Orthopedics, Changhai Hospital, Affiliated to Naval Medical University, Shanghai, China
| |
Collapse
|
4
|
Hu F, Wang Y, Zeng J, Deng X, Xia F, Xu X. Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations. J Phys Chem B 2024; 128:1418-1427. [PMID: 38323538 DOI: 10.1021/acs.jpcb.3c07666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2024]
Abstract
In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.
Collapse
Affiliation(s)
- Fangchen Hu
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Yiqiu Wang
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Juan Zeng
- School of Biomedical Engineering, Guangdong Medical University, Dongguan 523808, China
| | - Xianming Deng
- State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Faculty of Medicine and Life Sciences, Xiamen University, Xiamen, Fujian 361102, China
| | - Fei Xia
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai 200433, China
| |
Collapse
|
5
|
Bahar ME, Kim HJ, Kim DR. Targeting the RAS/RAF/MAPK pathway for cancer therapy: from mechanism to clinical studies. Signal Transduct Target Ther 2023; 8:455. [PMID: 38105263 PMCID: PMC10725898 DOI: 10.1038/s41392-023-01705-z] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 11/03/2023] [Accepted: 11/12/2023] [Indexed: 12/19/2023] Open
Abstract
Metastatic dissemination of solid tumors, a leading cause of cancer-related mortality, underscores the urgent need for enhanced insights into the molecular and cellular mechanisms underlying metastasis, chemoresistance, and the mechanistic backgrounds of individuals whose cancers are prone to migration. The most prevalent signaling cascade governed by multi-kinase inhibitors is the mitogen-activated protein kinase (MAPK) pathway, encompassing the RAS-RAF-MAPK kinase (MEK)-extracellular signal-related kinase (ERK) pathway. RAF kinase is a primary mediator of the MAPK pathway, responsible for the sequential activation of downstream targets, such as MEK and the transcription factor ERK, which control numerous cellular and physiological processes, including organism development, cell cycle control, cell proliferation and differentiation, cell survival, and death. Defects in this signaling cascade are associated with diseases such as cancer. RAF inhibitors (RAFi) combined with MEK blockers represent an FDA-approved therapeutic strategy for numerous RAF-mutant cancers, including melanoma, non-small cell lung carcinoma, and thyroid cancer. However, the development of therapy resistance by cancer cells remains an important barrier. Autophagy, an intracellular lysosome-dependent catabolic recycling process, plays a critical role in the development of RAFi resistance in cancer. Thus, targeting RAF and autophagy could be novel treatment strategies for RAF-mutant cancers. In this review, we delve deeper into the mechanistic insights surrounding RAF kinase signaling in tumorigenesis and RAFi-resistance. Furthermore, we explore and discuss the ongoing development of next-generation RAF inhibitors with enhanced therapeutic profiles. Additionally, this review sheds light on the functional interplay between RAF-targeted therapies and autophagy in cancer.
Collapse
Affiliation(s)
- Md Entaz Bahar
- Department of Biochemistry and Convergence Medical Sciences and Institute of Medical Science, Gyeongsang National University, College of Medicine, Jinju, South Korea
| | - Hyun Joon Kim
- Department of Anatomy and Convergence Medical Sciences and Institute of Medical Science, Gyeongsang National University, College of Medicine, Jinju, South Korea
| | - Deok Ryong Kim
- Department of Biochemistry and Convergence Medical Sciences and Institute of Medical Science, Gyeongsang National University, College of Medicine, Jinju, South Korea.
| |
Collapse
|
6
|
Zha J, He J, Wu C, Zhang M, Liu X, Zhang J. Designing drugs and chemical probes with the dualsteric approach. Chem Soc Rev 2023; 52:8651-8677. [PMID: 37990599 DOI: 10.1039/d3cs00650f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
Traditionally, drugs are monovalent, targeting only one site on the protein surface. This includes orthosteric and allosteric drugs, which bind the protein at orthosteric and allosteric sites, respectively. Orthosteric drugs are good in potency, whereas allosteric drugs have better selectivity and are solutions to classically undruggable targets. However, it would be difficult to simultaneously reach high potency and selectivity when targeting only one site. Also, both kinds of monovalent drugs suffer from mutation-caused drug resistance. To overcome these obstacles, dualsteric modulators have been proposed in the past twenty years. Compared to orthosteric or allosteric drugs, dualsteric modulators are bivalent (or bitopic) with two pharmacophores. Each of the two pharmacophores bind the protein at the orthosteric and an allosteric site, which could bring the modulator with special properties beyond monovalent drugs. In this study, we comprehensively review the current development of dualsteric modulators. Our main effort reason and illustrate the aims to apply the dualsteric approach, including a "double win" of potency and selectivity, overcoming mutation-caused drug resistance, developments of function-biased modulators, and design of partial agonists. Moreover, the strengths of the dualsteric technique also led to its application outside pharmacy, including the design of highly sensitive fluorescent tracers and usage as molecular rulers. Besides, we also introduced drug targets, designing strategies, and validation methods of dualsteric modulators. Finally, we detail the conclusions and perspectives.
Collapse
Affiliation(s)
- Jinyin Zha
- College of Pharmacy, Ningxia Medical University, Yinchuan, Ningxia Hui Autonomous Region, China.
- State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Jixiao He
- State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Chengwei Wu
- State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Mingyang Zhang
- State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Xinyi Liu
- College of Pharmacy, Ningxia Medical University, Yinchuan, Ningxia Hui Autonomous Region, China.
- State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Jian Zhang
- College of Pharmacy, Ningxia Medical University, Yinchuan, Ningxia Hui Autonomous Region, China.
- State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China
| |
Collapse
|
7
|
Zhang X, Liang W, Zheng G, Li B. Decoding the deactivation mechanism of R192W mutation of ZAP-70 using molecular dynamics simulations and binding free energy calculations. J Mol Model 2023; 29:371. [PMID: 37953318 DOI: 10.1007/s00894-023-05771-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 10/25/2023] [Indexed: 11/14/2023]
Abstract
CONTEXT ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation. METHODS In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.
Collapse
Affiliation(s)
- Xuehua Zhang
- Suzhou Hospital, Affiliated Hospital of Medical School, Nanjing University, Suzhou, China
| | - Wenqi Liang
- Department of Emergency, Changhai Hospital, The First Affiliated Hospital of Naval Medical University, Shanghai, 200433, China
| | - Guodong Zheng
- Department of VIP Clinic, Changhai Hospital, The First Affiliated Hospital to Naval Medical University, Shanghai, 200433, China.
| | - Bei Li
- Department of VIP Clinic, Changhai Hospital, The First Affiliated Hospital to Naval Medical University, Shanghai, 200433, China.
| |
Collapse
|
8
|
Decoding the Conformational Selective Mechanism of FGFR Isoforms: A Comparative Molecular Dynamics Simulation. Molecules 2023; 28:molecules28062709. [PMID: 36985681 PMCID: PMC10052029 DOI: 10.3390/molecules28062709] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Revised: 03/09/2023] [Accepted: 03/14/2023] [Indexed: 03/19/2023] Open
Abstract
Fibroblast growth factor receptors (FGFRs) play critical roles in the regulation of cell growth, differentiation, and proliferation. Specifically, FGFR2 gene amplification has been implicated in gastric and breast cancer. Pan-FGFR inhibitors often cause large toxic side effects, and the highly conserved ATP-binding pocket in the FGFR1/2/3 isoforms poses an immense challenge in designing selective FGFR2 inhibitors. Recently, an indazole-based inhibitor has been discovered that can selectively target FGFR2. However, the detailed mechanism involved in selective inhibition remains to be clarified. To this end, we performed extensive molecular dynamics simulations of the apo and inhibitor-bound systems along with multiple analyses, including Markov state models, principal component analysis, a cross-correlation matrix, binding free energy calculation, and community network analysis. Our results indicated that inhibitor binding induced the phosphate-binding loop (P-loop) of FGFR2 to switch from the open to the closed conformation. This effect enhanced extensive hydrophobic FGFR2-inhibitor contacts, contributing to inhibitor selectivity. Moreover, the key conformational intermediate states, dynamics, and driving forces of this transformation were uncovered. Overall, these findings not only provided a structural basis for understanding the closed P-loop conformation for therapeutic potential but also shed light on the design of selective inhibitors for treating specific types of cancer.
Collapse
|
9
|
Shi Y, Yasen M, Wang Z, Du T, Ding Y, Li X, Chai Z, Jie C, Ju G, Ji M. The allosteric effect of the upper half of SENP1 contributes to its substrate selectivity for SUMO1 over SUMO2. J Biomol Struct Dyn 2023; 41:12372-12386. [PMID: 36656084 DOI: 10.1080/07391102.2023.2166997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 01/03/2023] [Indexed: 01/20/2023]
Abstract
SUMOylation regulates various cellular process and SENP1 (SUMO-specific protease 1) serves as a SUMO (small ubiquitin-related modifier) specific protease that participates in the SUMO cycle. Given its extensive influences on metabolic activities, SENP1 has gained more and more attentions in clinical treatments. However, there remains a question on why does the SENP1 prefer to process SUMO1 rather than SUMO2. Here, we performed molecular dynamics simulations of SENP1-SUMO1, SENP1-SUMO2, and apo SENP1 systems and observed distinct conformational dynamics in the upper half of the clamp and the three loops in the catalytic center of the SENP1. Principal component analysis revealed that the most prominent canonical variable represented the spatial distribution of the upper half of the clamp, while the openness of the cleft was closely related to the catalytic ability of SENP1. Further analysis of the SENP1-SUMO interactions revealed that the extensive and strong interactions between the SENP1 and SUMO1 were both in the interface of the upper half region and the catalytic center. Dynamic cross-correlation matrix analysis demonstrated that the inter-residue correlations in the SUMO1 system was much stronger, especially in the two essential regions belonging to the upper and lower half of cleft. Based on these observations, we proposed an allosteric propagation model and further testified it using the community analysis. These results revealed the propagation pathway of allosteric communication that contributed to the substrate discrimination of SENP1 upon SUMO1 and SUMO2.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Yuxin Shi
- Department of Medicinal Chemistry and Bioinformatics Center, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
| | - Miersalijiang Yasen
- Department of Orthopedic Surgery, Zhongshan Hospital Xiamen Branch, Fudan University, Xiamen, Fujian, China
| | - Zhijun Wang
- Department of Urology, The Second Affiliated Hospital of Navy Medical University, Shanghai, China
| | - Tingting Du
- Department of VIP Clinic, Changhai Hospital, Navy Medical University, Shanghai, China
| | - Yelei Ding
- Department of Urology, The Second Affiliated Hospital of Navy Medical University, Shanghai, China
| | - Xuefei Li
- Department of Urology, Guixi People's Hospital, Guixi City, China
| | - Zongtao Chai
- Department of Liver Surgery and Transplantation, Key Laboratory of Carcinogenesis and Cancer Invasion of Ministry of Education, Liver Cancer Institute and Zhongshan Hospital, Fudan University, Shanghai, China
- Department of Hepatic Surgery, Shanghai Geriatric Cancer, Shanghai, China
| | - Chen Jie
- Department of Urology, The Second Affiliated Hospital of Navy Medical University, Shanghai, China
| | - Guanqun Ju
- Department of Urology, The Second Affiliated Hospital of Navy Medical University, Shanghai, China
| | - Mingfei Ji
- Department of Urology, The Second Affiliated Hospital of Navy Medical University, Shanghai, China
- Department of Urology, The First Affiliated Hospital of Soochow University, Suzhou, China
| |
Collapse
|
10
|
Cao S, Tan C, Fei A, Hu G, Fu M, Lv J. Insights into pralsetinib resistance to the non-gatekeeper RET kinase G810C mutation through molecular dynamics simulations. J Mol Model 2022; 29:24. [PMID: 36576611 DOI: 10.1007/s00894-022-05429-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Accepted: 12/19/2022] [Indexed: 12/29/2022]
Abstract
OBJECTIVE RET (rearranged during transfection) kinase, as a transmembrane receptor tyrosine kinase, is a therapeutic target for several human cancer such as non-small cell lung cancer (NSCLC) and thyroid cancer. Pralsetinib is a recently approved drug for the treatment of RET-driven NSCLC and thyroid cancers. A single point mutation G810C at the C-lobe of the RET kinase causes pralsetinib resistance to this non-gatekeeper variant. However, the detailed mechanism remains poorly understood. METHODS Here, multiple microsecond molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations, and community network analysis were performed to reveal the mechanism of pralsetinib resistance to the RET G810C mutant. RESULTS The simulations showed that the G810C mutation had a minor effect on the overall conformational dynamics of the RET kinase domain. Energetic analysis suggested that the G810C mutation reduced the binding affinity of pralsetinib to the mutant. Per-residue energy contribution and structural analyses revealed that the hydrogen bonding interactions between pralsetinib and the hinge residues Glu805 and Ala807 were disrupted in the G810C mutant, which were responsible for the decreased binding affinity of pralsetinib to the mutant. CONCLUSIONS The obtained results may provide understanding of the mechanism of pralsetinib resistance to the non-gatekeeper RET G810C mutant.
Collapse
Affiliation(s)
- Shu Cao
- Department of Urology, Ezhou Central Hospital, Hubei, 436000, China
| | - Changbin Tan
- Department of Urology, Ezhou Central Hospital, Hubei, 436000, China
| | - Anhua Fei
- Department of Urology, Ezhou Central Hospital, Hubei, 436000, China
| | - Gangqiang Hu
- Department of Urology, Ezhou Central Hospital, Hubei, 436000, China
| | - Ming Fu
- Department of Urology, Ezhou Central Hospital, Hubei, 436000, China
| | - Jun Lv
- Department of Neurology, Ezhou Central Hospital, Hubei, 436000, China.
| |
Collapse
|
11
|
Li X, Li B, Li J, Yang M, Bai Y, Chen K, Chen Z, Mao N. Mechanistic insights into the role of calcium in the allosteric regulation of the calmodulin-regulated death-associated protein kinase. Front Mol Biosci 2022; 9:1104942. [PMID: 36601586 PMCID: PMC9806222 DOI: 10.3389/fmolb.2022.1104942] [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/22/2022] [Accepted: 12/08/2022] [Indexed: 12/23/2022] Open
Abstract
Calcium (Ca2+) signaling plays an important role in the regulation of many cellular functions. Ca2+-binding protein calmodulin (CaM) serves as a primary effector of calcium function. Ca2+/CaM binds to the death-associated protein kinase 1 (DAPK1) to regulate intracellular signaling pathways. However, the mechanism underlying the influence of Ca2+ on the conformational dynamics of the DAPK1-CaM interactions is still unclear. Here, we performed large-scale molecular dynamics (MD) simulations of the DAPK1-CaM complex in the Ca2+-bound and-unbound states to reveal the importance of Ca2+. MD simulations revealed that removal of Ca2+ increased the anti-correlated inter-domain motions between DAPK1 and CaM, which weakened the DAPK1-CaM interactions. Binding free energy calculations validated the decreased DAPK1-CaM interactions in the Ca2+-unbound state. Structural analysis further revealed that Ca2+ removal caused the significant conformational changes at the DAPK1-CaM interface, especially the helices α1, α2, α4, α6, and α7 from the CaM and the basic loop and the phosphate-binding loop from the DAPK1. These results may be useful to understand the biological role of Ca2+ in physiological processes.
Collapse
Affiliation(s)
- Xiaolong Li
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China
| | - Bo Li
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China
| | - Jun Li
- Department of Orthopedics, Tongji Hospital, School of Medicine, Tongji University, Shanghai, China
| | - Mingyuan Yang
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China
| | - Yushu Bai
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China,*Correspondence: Yushu Bai, ; Kai Chen, ; Ziqiang Chen, ; Ningfang Mao,
| | - Kai Chen
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China,*Correspondence: Yushu Bai, ; Kai Chen, ; Ziqiang Chen, ; Ningfang Mao,
| | - Ziqiang Chen
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China,*Correspondence: Yushu Bai, ; Kai Chen, ; Ziqiang Chen, ; Ningfang Mao,
| | - Ningfang Mao
- Department of Orthopedics, Changhai Hospital, Naval Medical University, Shanghai, China,*Correspondence: Yushu Bai, ; Kai Chen, ; Ziqiang Chen, ; Ningfang Mao,
| |
Collapse
|
12
|
Zhuang H, Ji D, Fan J, Li M, Tao R, Du K, Lu S, Chai Z, Fan X. Mechanistic Insights into the Protection Effect of Argonaute-RNA Complex on the HCV Genome. Biomolecules 2022; 12:1631. [PMID: 36358979 PMCID: PMC9687641 DOI: 10.3390/biom12111631] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 10/25/2022] [Accepted: 10/31/2022] [Indexed: 11/14/2023] Open
Abstract
While host miRNA usually plays an antiviral role, the relentless tides of viral evolution have carved out a mechanism to recruit host miRNA as a viral protector. By complementing miR-122 at the 5' end of the genome, the hepatitis C virus (HCV) gene can form a complex with Argonaute 2 (Ago2) protein to protect the 5' end of HCV RNA from exonucleolytic attacks. Experiments showed that the disruption of the stem-loop 1(SL1) structure and the 9th nucleotide (T9) of HCV site 1 RNA could enhance the affinity of the Ago2 protein to the HCV site 1 RNA (target RNA). However, the underlying mechanism of how the conformation and dynamics of the Ago2: miRNA: target RNA complex is affected by the SL1 and T9 remains unclear. To address this, we performed large-scale molecular dynamics simulations on the AGO2-miRNA complex binding with the WT target, T9-abasic target and SL1-disruption target, respectively. The results revealed that the T9 and SL1 structures could induce the departing motion of the PAZ, PIWI and N domains, propping up the mouth of the central groove which accommodates the target RNA, causing the instability of the target RNA and disrupting the Ago2 binding. The coordinated motion among the PAZ, PIWI and N domains were also weakened by the T9 and SL1 structures. Moreover, we proposed a new model wherein the Ago2 protein could adopt a more constraint conformation with the proximity and more correlated motions of the PAZ, N and PIWI domains to protect the target RNA from dissociation. These findings reveal the mechanism of the Ago2-miRNA complex's protective effect on the HCV genome at the atomic level, which will offer guidance for the design of drugs to confront the protection effect and engineering of Ago2 as a gene-regulation tool.
Collapse
Affiliation(s)
- Haiming Zhuang
- Medicinal Chemistry and Bioinformatics Center, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Dong Ji
- Department of Anesthesiology, Changhai Hospital, Naval Medical University, Shanghai 200433, China
| | - Jigang Fan
- Medicinal Chemistry and Bioinformatics Center, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Mingyu Li
- Medicinal Chemistry and Bioinformatics Center, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Ran Tao
- School of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing 312000, China
| | - Kui Du
- School of Chemistry and Chemical Engineering, Shaoxing University, Shaoxing 312000, China
| | - Shaoyong Lu
- Medicinal Chemistry and Bioinformatics Center, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Zongtao Chai
- Key Laboratory of Carcinogenesis and Cancer Invasion of Ministry of Education, Department of Liver Surgery and Transplantation, Liver Cancer Institute and Zhongshan Hospital, Fudan University, Shanghai 200032, China
- Department of Hepatic Surgery, Shanghai Geriatric Cancer, Shanghai 201104, China
| | - Xiaohua Fan
- Department of Anesthesiology, Changhai Hospital, Naval Medical University, Shanghai 200433, China
| |
Collapse
|
13
|
Wang G, Bai Y, Cui J, Zong Z, Gao Y, Zheng Z. Computer-Aided Drug Design Boosts RAS Inhibitor Discovery. Molecules 2022; 27:molecules27175710. [PMID: 36080477 PMCID: PMC9457765 DOI: 10.3390/molecules27175710] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2022] [Revised: 08/13/2022] [Accepted: 08/26/2022] [Indexed: 11/16/2022] Open
Abstract
The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design.
Collapse
Affiliation(s)
- Ge Wang
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
- College of Stomatology, Shanghai Jiao Tong University, Shanghai 200120, China
| | - Yuhao Bai
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
- College of Stomatology, Shanghai Jiao Tong University, Shanghai 200120, China
| | - Jiarui Cui
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
- College of Stomatology, Shanghai Jiao Tong University, Shanghai 200120, China
| | - Zirui Zong
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
- College of Stomatology, Shanghai Jiao Tong University, Shanghai 200120, China
| | - Yuan Gao
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
- College of Stomatology, Shanghai Jiao Tong University, Shanghai 200120, China
| | - Zhen Zheng
- Medicinal Chemistry and Bioinformatics Center, Shanghai Jiao Tong University School of Medicine, Shanghai 200025, China
- Correspondence:
| |
Collapse
|
14
|
Zhang H, Ni D, Fan J, Li M, Zhang J, Hua C, Nussinov R, Lu S. Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation. J Chem Inf Model 2022; 62:4222-4231. [PMID: 35994329 DOI: 10.1021/acs.jcim.2c00591] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.
Collapse
Affiliation(s)
- Hao Zhang
- Department of Plastic and Reconstructive Surgery, Shanghai Ninth People's Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai 200011, China
| | - Duan Ni
- The Charles Perkins Centre, University of Sydney, Sydney, New South Wales 2006, Australia
| | - Jigang Fan
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Minyu Li
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Jian Zhang
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| | - Chen Hua
- Department of Plastic and Reconstructive Surgery, Shanghai Ninth People's Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai 200011, China
| | - Ruth Nussinov
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research, Cancer Innovation Laboratory, National Cancer Institute, Frederick, Maryland 21702, United States.,Department of Human Genetics and Molecular Medicine, Sackler School of Medicine, Sackler Institute of Molecular Medicine, Tel Aviv University, Tel Aviv 69978, Israel
| | - Shaoyong Lu
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China.,Medicinal Chemistry and Bioinformatics Centre, School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China
| |
Collapse
|
15
|
Berezovsky IN, Nussinov R. Multiscale Allostery: Basic Mechanisms and Versatility in Diagnostics and Drug Design. J Mol Biol 2022; 434:167751. [PMID: 35863488 DOI: 10.1016/j.jmb.2022.167751] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Igor N Berezovsky
- Bioinformatics Institute (BII), Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01, Matrix, Singapore 138671, Singapore; Department of Biological Sciences (DBS), National University of Singapore (NUS), 8 Medical Drive, 117579, Singapore.
| | - Ruth Nussinov
- Computational Structural Biology Section, Frederick National Laboratory for Cancer Research in the Cancer Innovation Laboraory, National Cancer Institute, Frederick, MD 21702, USA; Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Tel Aviv 69978, Israel.
| |
Collapse
|