Maximal Operators Associated to Hilbert Transforms along Variable Parabolas in $\mathbb{R}^2$