Powerful technique as it is to study the collective fluctuations and structural transitions in under-damped conservative systems, normal mode analysis cannot be applied to the optical matter (OM) system, an overdamped system subject to non-conservative forces, which consists of (nano-)particle constituents in solution that can self-organize into ordered arrays bound by electrodynamic interactions. We propose a data-driven approach based on principal component analysis (PCA) of deviations from a reference structure to determine the soft collective modes of non-conservative overdamped systems, such as OM structures, and harmonic linear discriminant analysis (HLDA) of time trajectories to estimate the reaction coordinate for structural transitions, which is demonstrated for structural isomerization of a six-particle OM system. The reaction coordinate we discover is in accord with committor analysis and the identified mechanism is in agreement with experiment. The PCA-HLDA approach to data-driven discovery of reaction coordinates aids in the understanding of non-conservative and overdamped systems.