Juvenile myelomonocytic leukemia (JMML), an invasive myeloproliferative neoplasm, is a childhood disease with very high clinical lethality. Somatic mutation E76K in SHP2 is the most commonly identified mutation found in up to 35% of patients with JMML. To investigate the effect of gain-of-function mutation-E76K on SHP2 activity, molecular dynamic simulations on the wild-type SHP2 (SHP2-WT) system and the mutated E76K (SHP2-E76K) system were performed. The evaluation of stability of these two systems indicated that the simulated trajectories were stable after simulation for 3 nanoseconds. The root mean square fluctuation and the per-residue root mean square deviation illustrated that there were two regions (residues Tyr 81-Glu 83 and Glu 258-Leu 261) in the wild-type system and the mutated system, which had large differences. The principal component analysis, dynamic cross correlation maps analysis, as well as secondary structure analysis suggested that the mutated E76K impacted the movement of these two regions in SHP2 protein. Furthermore, residue interaction network analysis, hydrogen bond occupancy, and binding free energies analysis were used to explain how the two regions were specifically affected by the mutant. The results indicated that the primary variances between SHP2-WT and SHP2-E76K were the different interactions between Glu/Lys 76 and Arg 265, Tyr 80 and Leu 77, Leu 77 and Tyr 81, Thr 73 and Glu 258, Ala 75 and Cys 259, Phe 71 and Tyr 81, Ala 75 and Glu 258, and Tyr 73 and Glu/Lys 76. Consequently, these findings here might provide insights into the increased activity in SHP2-E76K.