In this article, we investigate the contribution of the higher twist Feynman diagrams to the large-$p_T$ inclusive pion production cross section in proton-proton collisions and present the general formulae for the higher twist differential cross sections in the case of the running coupling and frozen coupling approaches. The structure of infrared renormalon singularities of the higher twist subprocess cross section and the resummed expression (the Borel sum) for it are found. We compared the resummed higher twist cross sections with the ones obtained in the framework of the frozen coupling approximation and leading twist cross section. We obtain, that ratio $R$ for all values of the transverse momentum $p_{T}$ of the pion identical equivalent to ratio $r$. It is shown that the resummed result depends on the choice of the meson wave functions used in calculation. Phenomenological effects of the obtained results are discussed.