We discuss a general five-dimensional completely anisotropic holographic model with three different spatial scale factors, characterized by a Van der Waals-like phase transition between small and large black holes. A peculiar feature of the model is the relation between anisotropy of the background and anisotropy of the colliding heavy ions geometry. We calculate the holographic entanglement entropy (HEE) of the slab-shaped region, the orientation of which relatively to the beams line and the impact parameter is characterized by the Euler angles. We study the dependences of the HEE and its density on the thermodynamic (temperature, chemical potential) and geometric (parameters of anisotropy, thickness, and orientation of entangled regions) parameters. As a particular case the model with two equal transversal scaling factors is considered. This model is supported by the dilaton and two Maxwell fields. In this case we discuss the HEE and its density in detail: interesting features of this model are jumps of the entanglement entropy and its density near the line of the small/large black hole phase transition. These jumps depend on the anisotropy parameter, chemical potential, and orientation. We also discuss different definitions and behavior of c-functions in this model. The c-function calculated in the Einstein frame decreases while increasing $ell$ for all $ell$ in the isotropic case (in regions of $(mu,T)$-plane far away from the line of the phase transition). We find the non-monotonicity of the c-functions for several anisotropic configurations, which however does not contradict with any of the existing c-theorems since they all base on Lorentz invariance.