We study the subelliptic heat kernels of the CR three-dimensional solvable Lie groups. We first classify all left-invariant sub-Riemannian structures on three-dimensional solvable Lie groups and obtain representations of these groups. We give expressions for the heat kernels on these groups and obtain heat semigroup gradient bounds using a new type of curvature-dimension inequality.