The isospin, spin and parity dependent potential of a pair of $B$ mesons is computed using Wilson twisted mass lattice QCD with two flavours of degenerate dynamical quarks. The $B$ meson is addressed in the static-light approximation, i.e.\ the $b$ quarks are infinitely heavy. From the results of the $B\,B$ meson-meson potentials, a simple rule can be deduced stating which isospin, spin and parity combinations correspond to attractive and which to repulsive forces. We provide fits to the ground state potentials in the attractive channels and discuss the potentials in the repulsive and excited channels. The attractive channels are most important since they can possibly lead to a bound four-quark state, i.e.\ a $\bar{b}\bar{b}ud$ tetraquark. Using these attractive potentials in the Schr\"odinger equation, we find indication for such a tetraquark state of two static bottom antiquarks and two light $u/d$ quarks with mass extrapolated down to the physical value.