In a series of three papers in the 80’s, Kusuoka and Stroock developed a probabilistic program in order to obtain sharp bounds for the density function of a hypoelliptic SDE driven by a Brownian motion. We aim to investigate how their method can be used to study rough SDEs driven by fractional Brownian motions. In this talk, I will outline Kusuoka and Stroock’s approach and point out where the difficulties are in our current setting.