Shape corrections to the standard approximate Kohn-Sham exchange-correlation ͑xc͒ potentials are considered with the aim to improve the excitation energies ͑especially for higher excitations͒ calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection ͑GRAC͒ of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends ͑LB͒ are seamlessly connected to the ͑shifted͒ xc potential of Becke and Perdew ͑BP͒ with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N 2 , CO, CH 2 O, C 2 H 4 , C 5 NH 5 , C 6 H 6 , Li 2 , Na 2 , K 2 . The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of ͑model͒ orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine ''local'' density functional with an analytical representation.