Exploring KRas Protein Dynamics: An Integrated Molecular Dynamics Analysis of KRas Wild and Mutant Variants

Abstract

This study employs a comprehensive approach combining protein retrieval, sequence alignment, and molecular dynamics simulations to investigate the structural dynamics and stability of wild-type KRas and its mutated variants (G12C, G12D, G12V, and G13D). The selected protein structures were retrieved from the Protein Data Bank (PDB) and prepared by using visual molecular dynamics (VMD) software. Sequence alignment using Clustal Omega provided a detailed comparison of the amino acid sequences, focusing on key mutation sites. Molecular dynamics simulations, performed with Gromacs, revealed distinct conformational changes and stability patterns in the wild-type and mutated KRas proteins over 100 ns. Clustering analysis identified higher conformational changes in the second α-helix of the mutated variants. The root-mean-square deviation (RMSD) distribution analysis showed variant-specific conformational dynamics, with G12V and G12D exhibiting slightly higher average RMSD values. Furthermore, clustering and RMSD analyses of specific amino acid residues (12, 13, 51, and 118) highlighted their roles in maintaining overall stability and influencing structural dynamics. The results indicate that mutations at positions 12 and 13 disrupt normal cycling between wild and mutated variants, leading to the persistent activation of KRas. Additionally, principal component analysis (PCA) elucidated unique conformational dynamics in mutated variants. Free energy landscape (FEL) analysis revealed alterations in the thermodynamic stability of mutated variants compared with the wild type. Overall, this study provides a detailed understanding of the structural changes associated with oncogenic mutations in KRas, offering insights crucial for targeted therapeutic strategies in KRas-driven cancers.