Domain mobility plays essential role in the biological function of multidomain systems. The characteristic times of domain motions fall into the interval from nano-to milliseconds, amenable to NMR studies. Proper analysis of NMR relaxation data for these systems in solution has to account for interdomain motions, in addition to the overall tumbling and local intradomain dynamics. Here we propose a model of interdomain mobility in a multi-domain protein, which considers domain reorientations as exchange/interconversion between two distinct conformational states of the molecule, combined with fully anisotropic overall tumbling. Analysis of 15 N-relaxation data for Lys48-linked di-ubiquitin at pH4.5 and 6.8 showed that this model adequately fits the experimental data and allows characterization of both structural and motional properties of di-ubiquitin, thus providing information about the relative orientation of ubiquitin domains in both interconverting states. The analysis revealed that the two domains reorient on a time scale of 9-30 ns, with the amplitudes sufficient for allowing a protein ligand access to the binding sites sequestered at the interface in the closed conformation. The analysis of possible mechanism controlling the equilibrium between the interconverting states in di-ubiquitin points toward protonation of His68, which results in three different charged states of the molecule, with zero, +e, and +2e net charge. Only two of the three states are noticeably populated at pH4.5 or 6.8, which assures applicability of the two-state model to di-ubiquitin at these conditions. We also compare our model with the "extended modelfree" approach and discuss possible future developments of the model.