Adsorption on activated carbons of five pharmaceutical molecules (ibuprofen, diclofenac, naproxen, paracetamol, and amoxicillin) in aqueous mixtures has been investigated by molecular simulations using the Grand Canonical Monte Carlo (GCMC) method. A virtual nanoporous carbon model based on polyaromatic units with defects and polar-oxygenated sites was used for this purpose. The simulation results show excellent agreement with available experimental data. The adsorption capacities of the carbons for the five drugs were quite different and were linked, essentially, to their molecular dimensions and atom affinities. The uptake behavior follows the trend PRM > DCF, NPX > IBP > AMX in all the studied structures. This work is a further step in order to describe macroscopic adsorption performance of activated carbons in drug removal applications.