A search is presented for a heavy pseudoscalar boson (\text{A}) decaying to a Z boson and a Higgs boson with mass of 125 (\phantom{\rule{0.166667em}{0ex}}\text{GeV}). In the final state considered, the Higgs boson decays to a bottom quark and antiquark, and the Z boson decays either into a pair of electrons, muons, or neutrinos. The analysis is performed using a data sample corresponding to an integrated luminosity of 35.9 (\phantom{\rule{0.166667em}{0ex}}{\text{fb}}^{-1}) collected in 2016 by the CMS experiment at the LHC from proton–proton collisions at a center-of-mass energy of 13 (\phantom{\rule{0.166667em}{0ex}}\text{Te}\phantom{\rule{0.333333em}{0ex}}). The data are found to be consistent with the background expectations. Exclusion limits are set in the context of two-Higgs-doublet models in the (\text{A}) boson mass range between 225 and 1000 (\phantom{\rule{0.166667em}{0ex}}\text{GeV}).