Theoretical description of molecular interactions remains a challenge for computational chemistry. In particular, systems dominated by static correlation, in which bonds are stretched or twisted, are often beyond capabilities of methods based on a single-electron approximation, being usually a method of choice. For interacting multireference systems, it is necessary to apply either high-level coupled-cluster methods (which however include multireference effects only partially) or base the theoretical description on a multireference wavefunction. Popular multireference methods like CASPT2 (complete active space perturbation theory) do not provide satisfactory results since they may suffer from problems with size consistency and poor accuracy. Recently we have shown that combining a simple multireference wavefunction, perfect-pairing generalized valence bond (GVB) with extended random phase approximation in embedding framework (EERPA) leads to a method EERPA-GVB providing accurate results for challenging multireference systems. In this paper, good performance of EERPA-GVB is confirmed by its application to van der Waals and hydrogen-bonded complexes. In addition, we show that the decomposition of the EERPA-GVB correlation energy into contributions from pairs of geminals can provide useful insight into the investigated interactions.