The development of supersonic combustion engines (scramjet-type) presents several challenges. Numerical simulations of scramjet engine combustion have become an attractive alternative to experimental investigations. However, supersonic combustion simulations still have several challenges: such as adequate modelling of shock/boundary layer and turbulence/chemistry interactions. The present work presents two large-eddy simulation probability density function (LES-PDF) novel formulations developed for high-speed applications. One is a conservative joint-scalar approach, where the joint-probability for reactive scalars and energy is solved, while the other is a joint velocity-scalar. The LES-PDF transport equations are solved using the Eulerian stochastic fields technique implemented in a density-based compressible solver. The performance of the models is verified through the simulations of two and three-dimensional supersonic reacting mixing layers and compared against DNS data from the literature. The results show that the joint-scalar formulation is accurate and robust, while the joint velocity-scalar closures require further development.