The study of intermolecular interactions of proteins has been an important problem for many years. This paper presents an approach to analyze different levels of protein interactions in solutions through a set of the second-and higher-order virial coefficients. The proposed approach is based on the diversified analysis of protein translational collective diffusion and self-diffusion obtained by dynamic light scattering and the pulsed-field gradient NMR (PFG NMR) spectroscopy experimental data. The experimental results were analyzed within the theoretical approach based on Vink's frictional formalism of nonequilibrium thermodynamics and the standard Derjaguin-Landau-Verwey-Overbeekb (DLVO) theory of interactions of colloid particles in electrolyte solutions. The second- and higher-order virial coefficients were obtained to estimate the pairwise and many-body intermolecular interactions in the solutions of globular alpha-chymotrypsin and intrinsically unstructured alpha(s)-casein. The second virial coefficients were calculated from the model of the protein-protein potential of mean force. The description of protein-protein interactions includes a set of interaction potentials: the attractive charge-dipole, dipole-dipole, the dispersion Hamaker, the mean force osmotic-attraction, and the repulsive charge-charge ones. It has been found that the major contribution to the intermolecular alpha(s)-casein interactions is made by the repulsive charge-charge potential, whereas for the case of alpha-chymotrypsin, the contributions from other types of interaction are of importance. It was determined that the model was well suited to describe the interactions of both globular and intrinsically disordered proteins. The suggested combination of Vink's approach and the DLVO theory is novel and holds much promise to make a profound analysis of the processes in systems containing various types of protein molecules.