We present the first application of the angle-dependent 3-Point Correlation Function (3PCF) to the density fields magnetohydrodynamic (MHD) turbulence simulations intended to model interstellar (ISM) turbulence. Previous work has demonstrated that the angle-averaged bispectrum, the 3PCF's Fourier-space analog, is sensitive to the sonic and Alfv\'enic Mach numbers of turbulence. Here we show that introducing angular information via multipole moments with respect to the triangle opening angle offers considerable additional discriminatory power on these parameters. We exploit a fast, order $N_{\rm g} \log N_{\rm g}$ ($N_{\rm g}$ the number of grid cells used for a Fourier Transform) 3PCF algorithm to study a suite of MHD turbulence simulations with 10 different combinations of sonic and Alfv\'enic Mach numbers over a range from sub to super-sonic and sub to super-Alfv\'{e}nic. The 3PCF algorithm's speed for the first time enables full quantification of the time-variation of our signal: we study 9 timeslices for each condition, demonstrating that the 3PCF is sufficiently time-stable to be used as an ISM diagnostic. In future, applying this framework to 3-D dust maps will enable better treatment of dust as a cosmological foreground as well as reveal conditions in the ISM that shape star formation.