The present investigation has dealt with co-biosorption of zinc and copper on the surface of eggshell powder (ESP) in liquid phase. Various process parameters such as pH, temperature, initial concentration of metal ions, agitation rate, contact time, and ESP dose were optimized to obtain the maximum removal of copper (77.36%) and zinc (44.24%). The optimized values of pH, temperature, initial concentration of metal ions, contact time, agitation rate, and ESP dose were 5,308K, 60mg/L, 50min, 180rpm, and 2mg/L, respectively. Various isotherms such as Langmuir, Freundlich, and Temkin model were evaluated. The results indicated that Langmuir model has better suitability in explaining co-biosorption of copper and zinc in terms of higher linear regression (R-2) in range of 0.94-0.99 together with very low values of (2) and sum of square errors from 0.007 to 0.01 and from 0.11 to 0.44, respectively. Contrary to this, the range of linear regression coefficient (R-2) reported for Freundlich and Temkin was quite low 0.74-0.94 and 0.66-0.91, coupled with very high values statistical error function. The preferential order of metal ion biosorption was Cu>Zn. The superior biosorption of Cu over Zn was due to their difference in the molecular weight, atomic radii, and electronegative. The surface characterization revealed the presence of macroporous texture with open void spaces rendering ESP as less efficient biosorbent for the removal of heavy metals such as Cu and Zn.