Recent experimental studies of the multiplicity dependence of heavy quark (HQ) production in proton-proton collisions at 7 TeV showed a strong non-linear increase of the HQ multiplicity as a function of the charged particle multiplicity. We try to understand this behavior using the EPOS3 approach. Two issues play an important role: multiple scattering, in particular its impact on multiplicity fluctuations, and the collective hydrodynamic expansion. The results are very robust with respect to many details of the modeling, which means that these data contain valuable information about very basic features of the reaction mechanism in proton-proton collisions.