Predicting the Conformational Variability of Oncogenic GTP-bound G12D Mutated KRas-4B Proteins at Cell Membranes

Abstract

KRas proteins are the largest family of mutated Ras isoforms, participating in a wide variety of cancers. Due to their importance, large effort is being carried out on drug development by small-molecule inhibitors. However, understanding protein conformational variability remains a challenge in drug discovery. In the case of the Ras family, their multiple conformational states can affect the binding of potential drug inhibitors. To overcome this challenge, we propose a computational framework based on combined all-atom Molecular Dynamics and Metadynamics simulations able to accurately access conformational variants of the target protein. We tested the methodology using a G12D mutated GTP bound oncogenic KRas-4B protein located at the interface of a DOPC/DOPS/cholesterol model anionic cell membrane. Two main orientations of KRas-4B at the anionic membrane have been obtained and explored. The corresponding angles have been taken as reliable reaction coordinates so that free-energy landscapes have been obtained by well-tempered metadynamics simulations, revealing the local and global minima of KRas-4B binding to the cell membrane, unvealing reactive paths of the system between the two preferential orientations and highlighting opportunities for targeting the unique metastable states through the identification of druggable pockets.