In terms of L-LTA4 systems, the calculated average rmsd value of K565A system was 0.158 nm which is lower than the rmsd values of WT (0.189 nm), E272Q (0.186 nm), and R563A (0.202 nm) systems ( Figure 4A ). The R563A system has shown the higher rmsd value and thereby indicating the additional effect of this particular mutation. Despite of these small differences in the rmsd values between the L-LTA4 systems, the systems were stabilized well throughout the timescale of simulation. As another method to investigate the stability of the systems, rmsf values of all systems were calculated during the simulation and plotted. From the plot, it was observed that none of the active site residues were fluctuating higher than 0.2 nm and explained the stable nature of the systems over the time scale of simulation. In addition, the number of intramolecular hydrogen bonds were calculated for all the systems and plotted. The average number of hydrogen bonds revealed that R563A system has formed more number of hydrogen bonds (466.9) compared to WT (452.8), E272Q (451.7) and R565A (450.6) systems which displayed the reduced number of hydrogen bonds. This result also has confirmed the stability of the systems despite of small differences between systems ( Figure 4A ).