Galaxy clusters are good targets to search for signals from dark matter (DM)rnannihilation process as they are dominatly composed of DM. From DM annihilationrnprocess neutral and charged pions are beleived to be produced. The neutral pionsrndecay directly into gamma-rays, while the secondary particle e_ produced from therncharged pions are expected to emit gamma-rays through inverse compton scatteringrn(ICS) and non-thermal bremsstrahlung processes. Taking energetic particles thatrnare e_ produced from DM annihilation process for the model of the origin of thernextended secondary particle observed in radio halos, the corresponding gamma-rayrnemission in those regions, we have studied both their morphology and the spectralrnenergy distribution using simulated Coma like galaxy clusters. Using raw galaxyrncluster data from large set of MUltiDark SImulation of galaxy Clusters (MUSIC-rn2) dataset simulation we choose three galaxy clusters for calculating gamma-rayrnemmission from DM annhilation process. For our calculation we have consideredrna neutralino annihilation of bb species and mass of DM model of neutralino massrnof 35 GeV with annihilation cross-section times velocity of = 1_ 10À€€26 cm3rnsð€€€1. Among the three simulated galaxy clusters, simulated galaxy cluster SGC282rnshowes the densest DM distribution almost in all the regions of the map. Therndifferential gamma-ray flux obtained from neutral pion decay process is on averagernan order of magnitude higher than that from non-thermal Bremsstrahlung process,rnfor similar photon energy, for all the three simulated galaxy clusters. For all thernthree gamma-ray emission mechanisims from DM annihilation process, simulatedrngalaxy cluster SGC282 showes dominant differential gamma-ray flux densities.rnOverall we have observed that the gamma-ray emission flux calculated in thesernsimulated galaxy cluster from the DM annihilation process is higher from regionsrnof denser DM halos in the simulated Coma like galaxy cluster as expected.