We present an isogeometric collocation method for solving the biharmonic equation over planar bilinearly parameterized multi-patch domains. The developed approach is based on the use of the globally C4-smooth isogeometric spline space (Kapl and Vitrih, 2021) to approximate the solution of the considered partial differential equation, and proposes as collocation points two different choices, namely on the one hand the Greville points and on the other hand the so-called superconvergent points. Several examples demonstrate the potential of our collocation method for solving the biharmonic equation over planar multi-patch domains, and numerically study the convergence behavior of the two types of collocation points with respect to the L2-norm as well as to equivalents of the Hs-seminorms for 1≤s≤4. In the studied case of spline degree p=9, the numerical results indicate in case of the Greville points a convergence of order O(hp−3) independent of the considered (semi)norm, and show in case of the superconvergent points an improved convergence of order O(hp−2) for all (semi)norms except for the equivalent of the H4-seminorm, where the order O(hp−3) is anyway optimal.