We consider a discrete particle system of two species coupled through nonlocal interactions driven by the one-dimensional Newtonian potential, with repulsive self-interaction and attractive cross-interaction. After providing a suitable existence theory in a finite-dimensional framework, we explore the behaviour of the particle system in case of collisions and analyse the behaviour of the solutions with initial data featuring particle clusters. Subsequently, we prove that the empirical measure associated to the particle system converges to the unique 2-Wasserstein gradient flow solution of a system of two partial differential equations (PDEs) with nonlocal interaction terms in a proper measure sense. The latter result uses uniform estimates of theĀ Lm-norms of a piecewise constant reconstruction of the density using the particle trajectories.