Learning Closed-form Equations for Subgrid-scale Closures from High-fidelity Data: Promises and Challenges