In this paper, we first prove that the existence of a solution of SDEs under the assumptions that the drift coefficient is of linear growth and path-dependent, and diffusion coefficient is bounded, uniformly elliptic and Hölder continuous. We apply Gaussian upper bound for a probability density function of a solution of SDE without drift coefficient and local Novikov condition, in order to use Maruyama-Girsanov transformation. The aim of this paper is to prove the existence with explicit representations (under linear/super-linear growth condition), Gaussian two-sided bound and Hölder continuity (under sub-linear growth condition) of a probability density function of a solution of SDEs with path-dependent drift coefficient. As an application of explicit representation, we provide the rate of convergence for an Euler-Maruyama (type) approximation, and an unbiased simulation scheme. = σ(t, Y 0,x t )dW t , Y 0,x 0 = x and "local" Novikov condition, in order to use Maruyama-Girsanov transformation. We will also show that the existence of a pdf