Modeling the co-evolution of supermassive black holes and galaxies: I. Black hole scaling relations and the Active Galactive Nuclei luminosity function