I think the problem is that you are mixing up terminology. The statment about an estimator being BLUE states given the Gauss-Markov assumption, OLS is the best linear unbiased estimator. But as your model is heteroskedastic the Gauss-Markov assumptions no longer hold and the proof of OLS being BLUE is no longer true.
Unfortunalty there is no general result that shows whether a estimater is the "best" in the case of hetroskedacity. In theory if you know the exact functional form of the heteroskedascity you can perfectly correct for the hetroskadicity and you WLS is just as effective as OLS when the Gauss-Markov assumptions hold. But in reality this is never the case, and unless you have a small sample or a very strong argument of why you know the functional form of the hetroskedacity then you are better of using Whites robust standart errors.