Flow Matching for Optimal Reaction Coordinates of Biomolecular System