The construction of AdS/QCD models for heavy quarkonia is a confusing problem because phenomenologically it is not clear which pattern of spectrum should they describe in the first approximation. The second problem lies in the fact that there is no universal way for incorporating the quark masses into such models. We propose a solution to the second problem. One can consider the static limit with very large quark masses when the Coulomb interaction dominates and thus determines in the first approximation the spectrum of quarkonia. We make use of this limit for identification of a relevant bottom-up holographic model. It is shown that such a model has a simple structure: it represents the standard soft-wall holographic model, in which, however, the quadratic dilaton background is replaced by the linear one. The given replacement converts a Regge-like spectrum into a Hydrogen-like one. The mass scale introduced by the linear dilaton matches the quark mass in a natural way. The resulting model is analyzed for the scalar, vector and tensor cases. We also demonstrate that the spectrum of considered holographic model shares several common features with the recently found spectrum of open string scattering amplitude for strings ending on a D-brane in the AdS space.