In this paper, we discuss estimates of transition densities of subordinate Brownian motions in open subsets of Euclidean space. When $D$ is a $C^{1,1}$ domain, we establish sharp two-sided estimates for the transition densities of a large class of subordinate Brownian motions in $D$ whose scaling order is not necessarily strictly below $2$. Our estimates are explicit and written in terms of the dimension, the Euclidean distance between two points, the distance to the boundary and Laplace exponent of the corresponding subordinator only.