A batch sorption method was used to study the removal of few toxic metals onto the Late Cretaceous clays of Aleg formation (Coniacian-Lower Campanian system), Tunisia, in single, binary and multi-component systems. The collected clay samples were used as adsorbents for the removal of Pb(II), Cd(II), Cu(II) and Zn(II) from aqueous solutions. Results show that the natural clay samples were mainly composed of silica, alumina, iron and magnesium oxides. N 2 -adsorption measurements indicated mesoporous materials with modest specific surface area of <71 m 2 /g. Carbonate minerals were the most influencing parameters for heavy metal removal by natural clays in both single and multi-element systems. The affinity sequence was Pb(II)>Cu(II)>Zn(II)>Cd(II) due to the variable physical properties of the studied metals. The maximum adsorption capacity reached 131.58 mg/g in single systems, but decreased to <50.10 mg/g in mixed systems. In single, binary and muti-element systems, the studied clay samples removed substantial amounts of heavy metals, showing better effectiveness than the relevant previous studies. These results suggest that the studied clay samples of the Late Cretaceous clays from Tunisia can be effectively used as natural adsorbents for the removal of toxic heavy metals in aqueous systems.