Modeling the role of G12V and G13V Ras mutations in the Ras-GAP-catalyzed hydrolysis reaction of guanosine triphosphate.

Cancer-associated point mutations in Ras, in particular, at glycine 12 and glycine 13, affect the normal cycle between inactive GDP-bound and active GTP-bound states. In this work, the role of G12V and G13V replacements in the GAP-stimulated intrinsic GTP hydrolysis reaction in Ras is studied using molecular dynamics (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials. A model molecular system was constructed by motifs of the relevant crystal structure (Protein Data Bank entry 1WQ1 ). QM/MM optimization of geometry parameters in the Ras-GAP-GTP complex and QM/MM-MD simulations were performed with a quantum subsystem comprising a large fraction of the enzyme active site. For the system with wild-type Ras, the conformations fluctuated near the structure ready to be involved in the efficient chemical reaction leading to the cleavage of the phosphorus-oxygen bond in GTP upon approach of the properly aligned catalytic water molecule. Dynamics of the system with the G13V mutant is characterized by an enhanced flexibility in the area occupied by the γ-phosphate group of GTP, catalytic water, and the side chains of Arg789 and Gln61, which should somewhat hinder fast chemical steps. Conformational dynamics of the system with the G12V mutant shows considerable displacement of the Gln61 side chain and catalytic water from their favorable arrangement in the active site that may lead to a marked reduction in the reaction rate. The obtained computational results correlate well with the recent kinetic measurements of the Ras-GAP-catalyzed hydrolysis of GTP.