Reaction coordinate flows for model reduction of molecular kinetics