Journal of Shanghai Jiao Tong University (Medical Science) ›› 2025, Vol. 45 ›› Issue (6): 705-716.doi: 10.3969/j.issn.1674-8115.2025.06.005

• Basic research • Previous Articles     Next Articles

Study on the mechanism of KRAS R68G secondary mutation-induced resistance to KRASG12D-targeted inhibitor MRTX1133

WANG Gaoming1,2,3, CUI Ran4, LI Yanjing2,3(), LIU Yingbin1,2,3()   

  1. 1.Department of Biliary-Pancreatic Surgery, Renji Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai 200127, China
    2.State Key Laboratory of Systems Medicine for Cancer, Shanghai Cancer Institute, Renji Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai 200127, China
    3.Shanghai Key Laboratory for Cancer Systems Regulation and Clinical Translation (CSRCT-SHANGHAI), Renji Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai 200127, China
    4.Department of Hepatopancreatobiliary Surgery, East Hospital Affiliated Tongji University, Shanghai 200120, China
  • Received:2024-12-18 Accepted:2025-03-07 Online:2025-06-23 Published:2025-06-23
  • Contact: LI Yanjing, LIU Yingbin E-mail:liyanjing@sjtu.edu.cn;laoniulyb@163.com
  • Supported by:
    State Key Laboratory of Systems Medicine for Cancer(ZZ-RCPY-23-25);National Natural Science Foundation of China(32471528)

Abstract:

Objective ·To explore the mechanism at the atomic level by which the KRASG12D/R68G mutation induces tumor cell resistance to MRTX1133. Methods ·The crystal structure data of the KRASG12D-MRTX1133 complex were obtained from the RCSB Protein Data Bank (PDB). PyMOL software was used to mutate arginine at position 68 of KRAS to glycine (R68G), constructing the initial conformations of the KRASG12D-MRTX1133 and KRASG12D/R68G-MRTX1133 complexes. The LEaP module was used to build simulation systems under periodic boundary conditions. The ff19SB force field was applied to standard amino acids in KRAS, GAFF (general AMBER force field) to MRTX1133, and TIP3P (intermolecular potential three point) to water molecules. Energy minimization was performed using the Amber software suite. The systems were then heated to 300 K, followed by NVT (constant volume and temperature) equilibration and NPT (constant pressure and temperature) production. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), principal component analysis (PCA) and solvent-accessible surface area (SASA) of MRTX1133 and GDP were analyzed using cpptraj. The number of hydrogen bonds between regions and the dynamic cross-correlation matrix (DCCM) of amino acid movements were also calculated. Results ·RMSD analysis showed greater structural variation in KRAS in the KRASG12D/R68G system compared to the KRASG12D system. RMSF analysis revealed significantly higher fluctuations in the Switch Ⅰ and Switch Ⅱ regions of the KRASG12D/R68G system. PCA indicated that Switch Ⅰ and Switch Ⅱ in the KRASG12D/R68G system were more frequently in an open conformation. The distances between Switch Ⅰ and the P-loop, and between Switch Ⅱ and the P-loop, were larger in the KRASG12D/R68G system, indicating an expanded binding pocket for GDP and MRTX1133 compared to the KRASG12D system. SASA analysis indicated that both GDP and MRTX1133 had increased solvent exposure in the KRASG12D/R68G system. DCCM analysis revealed more decoupled movements among the Switch Ⅰ, Switch Ⅱ and P-loop regions in the KRASG12D/R68G system. Conclusion ·The KRASG12D/R68G mutation disrupts the interactions between the Switch Ⅰ and Switch Ⅱ regions, leading to their separation and the opening of the MRTX1133 binding pocket. This increases the solvent exposure of MRTX1133, accelerates its dissociation, and ultimately results in KRASG12D/R68G resistance to MRTX1133.

Key words: MRTX1133, KRAS secondary mutation, molecular dynamics simulation, resistance mechanism

CLC Number: