We calculate the membrane-induced interaction between inclusions, in terms of the membrane stretching and bending moduli and the spontaneous curvature. We find that the membrane-induced interaction between inclusions varies nonmonotonically as a function of the inclusion spacing. The location of the energy minimum depends on the spontaneous curvature and the membrane perturbation decay length, where the latter is set by the membrane moduli. The membrane perturbation energy increases with the inclusion radius. The Ornstein-Zernike theory, with the Percus-Yevick closure, is used to calculate the radial distribution function of inclusions. We find that when the spontaneous curvature is zero, the interaction between inclusions due to the membrane deformation is qualitatively similar to the hard-core interaction. However, in the case of finite spontaneous curvature, the effective interaction is dramatically modified.