We devise a supersymmetry-based method for the construction of zero-energy states in graphene. Our method is applied to a two-dimensional massless Dirac equation with a hyperbolic scalar potential. We determine supersymmetric partners of our initial system and derive a reality condition for the transformed potential. The Dirac potentials generated by our method can be used to approximate interactions that are experimentally realizable.