A generalization of the Keller–Segel model for chemotactic systems is studied. In this model there are several populations interacting via several sensitivity agents in a two-dimensional domain. The dynamics of the population is determined by a Fokker–Planck system of equations, coupled with a system of diffusion equations for the chemical agents. Conditions for global existence of solutions and equilibria are discussed, as well as the possible existence of time-periodic attractors. The analysis is based on a variational functional associated with the system.