This study presents a numerical procedure to optimize the shape of staggered dimple surface to enhance turbulent heat transfer in a rectangular channel. The RBNN method is used as an optimization technique with Reynolds-averaged Navier-Stokes analysis of fluid flow and heat transfer with shear stress transport (SST) turbulence model. The dimple depth-to-dimple print diameter (d/D), channel height-to-dimple print diameter ratio (H/D), and dimple print diameter-to-pitch ratio (D/S) are chosen as design variables. The objective function is defined as a linear combination of heat transfer related term and friction loss related term with a weighting factor. Latin Hypercube Sampling (LHS) is used to determine the training points as a mean of the design of experiment. The optimum shape shows remarkable performance in comparison with a reference shape.