1
|
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
|
2
|
Zeng J, Chen J, Xia F, Cui Q, Deng X, Xu X. Identification of functional substates of KRas during GTP hydrolysis with enhanced sampling simulations. Phys Chem Chem Phys 2022; 24:7653-7665. [PMID: 35297922 PMCID: PMC8972078 DOI: 10.1039/d2cp00274d] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
As the hub of major signaling pathways, Ras proteins are implicated in 19% of tumor-caused cancers due to perturbations in their conformational and/or catalytic properties. Despite numerous studies, the functions of the conformational substates for the most important isoform, KRas, remain elusive. In this work, we perform an extensive simulation analysis on the conformational landscape of KRas in its various chemical states during the GTP hydrolysis cycle: the reactant state KRasGTP·Mg2+, the intermediate state KRasGDP·Pi·Mg2+ and the product state KRasGDP·Mg2+. The results from enhanced sampling simulations reveal that State 1 of KRasGTP·Mg2+ has multiple stable substates in solution, one of which might account for interacting with GEFs. State 2 of KRasGTP·Mg2+ features two substates "Tyr32in" and "Tyr32out", which are poised to interact with effectors and GAPs, respectively. For the intermediate state KRasGDP·Pi·Mg2+, Gln61 and Pi are found to assume a broad set of conformations, which might account for the weak oncogenic effect of Gln61 mutations in KRas in contrast to the situation in HRas and NRas. Finally, the product state KRasGDP·Mg2+ has more than two stable substates in solution, pointing to a conformation-selection mechanism for complexation with GEFs. Based on these results, some specific inhibition strategies for targeting the binding sites of the high-energy substates of KRas during GTP hydrolysis are discussed.
Collapse
Affiliation(s)
- Juan Zeng
- School of Biomedical Engineering, Guangdong Medical University, Dongguan 523808, China
| | - Jian Chen
- School of Chemistry and Molecular Engineering, NYU-ECNU Center for Computational Chemistry at NYU Shanghai, East China Normal University, Shanghai 200062, China.
| | - Fei Xia
- School of Chemistry and Molecular Engineering, NYU-ECNU Center for Computational Chemistry at NYU Shanghai, East China Normal University, Shanghai 200062, China.
| | - Qiang Cui
- Departments of Chemistry, Physics and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA
| | - Xianming Deng
- State Key Laboratory of Cellular Stress Biology Innovation Center for Cell Signaling Network, School of Life Sciences, Xiamen University, Fujian 361101, 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
|
3
|
Zhang J, Liu Z, Zhao W, Yin X, Zheng X, Liu C, Wang J, Wang E. Discovery of Small Molecule NSC290956 as a Therapeutic Agent for KRas Mutant Non-Small-Cell Lung Cancer. Front Pharmacol 2022; 12:797821. [PMID: 35069209 PMCID: PMC8766838 DOI: 10.3389/fphar.2021.797821] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Accepted: 11/22/2021] [Indexed: 11/13/2022] Open
Abstract
HRas-GTP has a transient intermediate state with a “non-signaling open conformation” in GTP hydrolysis and nucleotide exchange. Due to the same hydrolysis process and the structural homology, it can be speculated that the active KRas adopts the same characteristics with the “open conformation.” This implies that agents locking this “open conformation” may theoretically block KRas-dependent signaling. Applying our specificity-affinity drug screening approach, NSC290956 was chosen by high affinity and specificity interaction with the “open conformation” structure HRasG60A-GppNp. In mutant KRas-driven non-small-cell lung cancer (NSCLC) model system, NSC290956 effectively suppresses the KRas-GTP state and gives pharmacological KRas inhibition with concomitant blockages of both the MAPK-ERK and AKT-mTOR pathways. The dual inhibitory effects lead to the metabolic phenotype switching from glycolysis to mitochondrial metabolism, which promotes the cancer cell death. In the xenograft model, NSC290956 significantly reduces H358 tumor growth in nude mice by mechanisms similar to those observed in the cells. Our work indicates that NSC290956 can be a promising agent for the mutant KRas-driven NSCLC therapy.
Collapse
Affiliation(s)
- Jiaxin Zhang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China.,Department of Chemistry, University of Science and Technology of China, Hefei, China
| | - Zuojia Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China
| | - Wenjing Zhao
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China
| | - Xunzhe Yin
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China
| | - Xiliang Zheng
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China
| | - Chuanbo Liu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China
| | - Jin Wang
- Department of Chemistry and Physics, State University of New York, Stony Brook, NY, United States
| | - Erkang Wang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, China.,Department of Chemistry, University of Science and Technology of China, Hefei, China
| |
Collapse
|
4
|
Tanaka S, Tokutomi S, Hatada R, Okuwaki K, Akisawa K, Fukuzawa K, Komeiji Y, Okiyama Y, Mochizuki Y. Dynamic Cooperativity of Ligand-Residue Interactions Evaluated with the Fragment Molecular Orbital Method. J Phys Chem B 2021; 125:6501-6512. [PMID: 34124906 DOI: 10.1021/acs.jpcb.1c03043] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
By the splendid advance in computation power realized with the Fugaku supercomputer, it has become possible to perform ab initio fragment molecular orbital (FMO) calculations for thousands of dynamic structures of protein-ligand complexes in a parallel way. We thus carried out electron-correlated FMO calculations for a complex of the 3C-like (3CL) main protease (Mpro) of the new coronavirus (SARS-CoV-2) and its inhibitor N3 incorporating the structural fluctuations sampled by classical molecular dynamics (MD) simulation in hydrated conditions. Along with a statistical evaluation of the interfragment interaction energies (IFIEs) between the N3 ligand and the surrounding amino-acid residues for 1000 dynamic structure samples, in this study we applied a novel approach based on principal component analysis (PCA) and singular value decomposition (SVD) to the analysis of IFIE data in order to extract the dynamically cooperative interactions between the ligand and the residues. We found that the relative importance of each residue is modified via the structural fluctuations and that the ligand is bound in the pharmacophore in a dynamic manner through collective interactions formed by multiple residues, thus providing new insight into structure-based drug discovery.
Collapse
Affiliation(s)
- Shigenori Tanaka
- Graduate School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Shusuke Tokutomi
- Graduate School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Ryo Hatada
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Koji Okuwaki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Kazuki Akisawa
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Kaori Fukuzawa
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan.,Department of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, 6-6-11 Aoba, Aramaki, Aoba-ku, Sendai 980-8579, Japan.,Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Yuto Komeiji
- Biomedical Research Institute, AIST, Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
| | - Yoshio Okiyama
- Division of Medicinal Safety Science, National Institute of Health Sciences, 3-25-26 Tonomachi, Kawasaki-ku, Kawasaki, Kanagawa 201-9501, Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan.,Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
5
|
Pálfy G, Menyhárd DK, Perczel A. Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery. Cancer Metastasis Rev 2020; 39:1075-1089. [PMID: 32815102 PMCID: PMC7680338 DOI: 10.1007/s10555-020-09917-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/22/2020] [Accepted: 06/22/2020] [Indexed: 12/11/2022]
Abstract
Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.
Collapse
Affiliation(s)
- Gyula Pálfy
- Laboratory of Structural Chemistry and Biology, Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117, Budapest, Hungary
- Protein Modeling Group HAS-ELTE, Institute of Chemistry, Eötvös Loránd University, P.O.B. 32, Budapest, 1538, Hungary
| | - Dóra K Menyhárd
- Laboratory of Structural Chemistry and Biology, Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117, Budapest, Hungary.
- Protein Modeling Group HAS-ELTE, Institute of Chemistry, Eötvös Loránd University, P.O.B. 32, Budapest, 1538, Hungary.
| | - András Perczel
- Laboratory of Structural Chemistry and Biology, Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, 1117, Budapest, Hungary.
- Protein Modeling Group HAS-ELTE, Institute of Chemistry, Eötvös Loránd University, P.O.B. 32, Budapest, 1538, Hungary.
| |
Collapse
|
6
|
Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics. Sci Rep 2018; 8:12976. [PMID: 30154518 PMCID: PMC6113251 DOI: 10.1038/s41598-018-31234-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Accepted: 08/13/2018] [Indexed: 12/14/2022] Open
Abstract
The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.
Collapse
|