Given a collection L of n subsets of a finite set R, the Maximum Intersection of k-Subsets problem (kMIS) aims to select a subset \(L' \subseteq L\) with \(|L'| = k\) such that the intersection of the subsets in \(L'\) is maximized. This problem is \(\mathcal{N}\mathcal{P}\) -hard and has applications in data anonymization, bioinformatics, recommendation systems, and various other fields. In this work, we propose an Iterated Local Search (ILS) algorithm for solving the kMIS. The local search phase uses the swap(1,1) neighborhood structure. To enhance efficiency, we introduce an auxiliary data structure which precomputes partial intersections, reducing the complexity of evaluating a neighboring solution from \(O(k^2)\) to O(k). The proposed ILS was evaluated against two state-of-the-art metaheuristics from the literature on 238 benchmark instances. The experimental results demonstrate that our approach efficiently finds high-quality solutions, establishing its competitiveness for solving the kMIS.