Comparative effects of oncogenic mutations G12C, G12V, G13D, and Q61H on local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated protein in human cancers. However, until very recently, its oncogenic mutants were viewed as undruggable. To develop inhibitors that directly target oncogenic K-Ras mutants, we need to understand both their mutant-specific and pan-mutant dynamics and conformations. Recently, we have investigated how the most frequently observed K-Ras mutation in cancer patients, G12D, changes its local dynamics and conformations (Vatansever et al., 2019). Here, we extend our analysis to study and compare the local effects of other frequently observed oncogenic mutations, G12C, G12V, G13D and Q61H. For this purpose, we have performed Molecular Dynamics (MD) simulations of each mutant when active (GTP-bound) and inactive (GDP-bound), analyzed their trajectories, and compared how each mutant changes local residue conformations, inter-protein distance distributions, local flexibility and residue pair correlated motions. Our results reveal that in the four active oncogenic mutants we have studied, the α2 helix moves closer to the C-terminal of the α3 helix. However, P-loop mutations cause α3 helix to move away from Loop7, and only G12 mutations change the local conformational state populations of the protein. Furthermore, the motions of coupled residues are mutant-specific: G12 mutations lead to new negative correlations between residue motions, while Q61H destroys them. Overall, our findings on the local conformational states and protein dynamics of oncogenic K-Ras mutants can provide insights for both mutant-selective and pan-mutant targeted inhibition efforts.