A search for the production of a Higgs boson in association with a single top quark is presented, focusing on leptonic signatures provided by the , , and decay modes. Due to strong interference of the two main leading-order diagrams, the production cross section of this process is highly sensitive to the relative sign of the top-Higgs coupling modifier, , and the coupling modifier of vector bosons to the Higgs, . The analysis exploits signatures with two same-sign leptons or three leptons in the final state, and uses the 2016 data sample collected with the CMS detector at the LHC at a center of mass energy of 13~\mathrm{TeV}, which corresponds to an integrated luminosity of 35.9~\mathrm{fb}^{-1}. Multivariate techniques are used to discriminate the signal from the dominant backgrounds. The analysis yields a confidence level (C.L.) upper limit on the combined production cross section times branching ratio of 0.64\,\mathrm{pb}, with an expected limit of 0.32\,\mathrm{pb}, for a scenario with and . Values of outside the range of to are excluded at C.L., assuming .