Experimental determination of electron density distribution in crystals by means of high-resolution X-ray diffraction allows, among others, for studying the details of intra- and inter-molecular interactions. In case of co-crystals, this method may help in finding the conditions of creating such species. The results of such analysis for two co-crystals containing betaines, namely trigonelline (TRG: nicotinic acid N-methylbetaine, IUPAC name: 1-methylpyridinium-3-carboxylate) and N-methylpiperidine betaine (MPB: 1-methylpiperidinium-1-yl-carboxylate) with p-hydroxybenzoic acid (HBA) are reported. TRG-HBA crystallizes as a hydrate. For both of the co-crystals, high-quality diffraction data were collected up to sinθ/λ = 1.13 Å−1. Hansen-Coppens multipolar model was then applied for modelling the electron density distribution and Atoms-In-Molecules approach was used for detailed analysis of interactions in crystals. A number of intermolecular interactions was identified, ranging from strong O-H···O hydrogen bonds through C-H···O to C-H···π and π···π interactions. Correlations between the geometrical characteristics of the contacts and the features of their critical points were analyzed in detail. Atomic charges show that in zwitterionic species there are regions of opposite charges, rather than charges that are localized on certain atoms. In case of MPB-HBA, a significant charge transfer between the components of co-crystal (0.5 e) was found, as opposed to TRG-HBA, where all of the components are almost neutral.