A model of independent pair parton interactions is proposed, according to which, hadron interactions are represented by a set of independent binary parton collisions. The final multiplicity distribution is described by a convolution of the negative binomial distributions in each of the partonic collisions. As a result, it is given by a weighted sum of negative binomial distributions with parameters multiplied by the number of active pairs. Its shape and moments are considered. Experimental data on multiplicity distributions in high energy pp processes are well fitted by these distributions. Predictions for the CERN Large Hadron Collider and higher energies are presented. The difference between e + e − and pp processes is discussed.