Specialized proteins serve as scaffolds sculpting strongly curved membranes of intracellular organelles. Effective membrane shaping requires segregation of these proteins into domains and is, therefore, critically dependent on the protein-protein interaction. Interactions mediated by membrane elastic deformations have been extensively analyzed within approximations of large inter-protein distances, small extents of the protein-mediated membrane bending and small deviations of the protein shapes from isotropic spherical segments. At the same time, important classes of the realistic membrane-shaping proteins have strongly elongated shapes with large and highly anisotropic curvature. Here we investigated, computationally, the membrane mediated interaction between proteins or protein oligomers representing membrane scaffolds with strongly anisotropic curvature, and addressed, quantitatively, a specific case of the scaffold geometrical parameters characterizing BAR domains, which are crucial for membrane shaping in endocytosis. In addition to the previously analyzed contributions to the interaction, we considered a repulsive force stemming from the entropy of the scaffold orientation. We computed this interaction to be of the same order of magnitude as the well-known attractive force related to the entropy of membrane undulations. We demonstrated the scaffold shape anisotropy to cause a mutual aligning of the scaffolds and to generate a strong attractive interaction bringing the scaffolds close to each other to equilibrium distances much smaller than the scaffold size. We computed the energy of interaction between scaffolds of a realistic geometry to constitute tens of kBT, which guarantees a robust segregation of the scaffolds into domains. 