Currently, the dye industry is increasing its production as a consequence of the growing need for their products in different manufacturing sectors, such as textiles, plastics, food, paper, etc... Thereafter, these industries generate very large volumes of effluents contaminated by these dyes, which require proper removal treatment before final discharge of the effluents into the environment. In this study, artichoke leaves were used as an economical and eco-friendly bio-adsorbent for Bengal Rose (BR) dye removal. Bio-adsorbent obtained from artichoke leaves was ground to powder size. The resulting powder was characterized by different methods, such as Brunauer-Emmett-Teller (BET) surface area analysis, scanning electron microscopy(SEM), X-ray Diffraction (XRD), Fourier transfer infrared (FTIR), pH at point of zero charge (pHpzc), equilibrium pH, iodine number, methylene blue number, phenol number, density, Energy dispersive X-ray spectroscopy (EDX) and Thermo-gravimetric analysis (TGA). Thereafter, the bio-adsorbent was used to study its capability for removing BR dye by testing contact time, initial concentration of dye and temperature. The results show that the saturation of bio-sorbent was reached after 40 min and the removal rate of BR dye by artichoke leaves powder (ALP) was 4.07 mg/g, which corresponds to a removal efficiency of 80.1%. A design of experiences (DOE) based on a two-level full factorial design (23) was used to study the effects of different parameters, such as pH, temperature and bio-adsorbent dosage on BR dye removal efficiency. The obtained results show that the highest removal efficiency was 86.5% for the optimized values of pH (4), temperature (80 °C) and bio-adsorbent dosage (8 g/L). Furthermore, a satisfying accordance between experimental and predicted data was observed. The kinetic and isotherm studies show that the pseudo-second order model simulated adequately the obtained data and it was found that Langmuir and Temkin isotherm models are liable and suitable for evaluating the adsorption process performance. Free energy change of adsorption (ΔG°), enthalpy change (ΔH°) and entropy change (ΔS°) were furthermore calculated to predict the nature of the adsorption process.