From Matrix Interpolation to Tensorized Simulation of High-Dimensional Random Variables: with Applications to Rare Event Estimation