Human α1,4-galactosyltransferase (A4galt), a Golgi apparatus-resident GT, synthesizes Gb3 glycosphingolipid (GSL) and P1 glycotope on glycoproteins (GPs), which are receptors for Shiga toxin types 1 and 2. Despite the significant role of A4galt in glycosylation processes, the molecular mechanisms underlying its varied acceptor specificities remain poorly understood. Here, we attempted to elucidate A4galt specificity towards GSLs and GPs by exploring its interaction with GTs with various acceptor specificities, GP-specific β1,4-galactosyltransferase 1 (B4galt1) and GSL-specific β1,4-galactosyltransferase isoenzymes 5 and 6 (B4galt5 and B4galt6). Using a novel NanoBiT assay, we found that A4galt can form homodimers and heterodimers with B4galt1 and B4galt5 in two cell lines, human embryonic kidney cells (HEK293T) and Chinese hamster ovary cells (CHO-Lec2). We found that A4galt-B4galts heterodimers preferred N-terminally tagged interactions, while in A4galt homodimers, the favored localization of the fused tag depended on the cell line used. Furthermore, by employing AlphaFold for state-of-the-art structural prediction, we analyzed the interactions and structures of these enzyme complexes. Our analysis highlighted that the A4galt-B4galt5 heterodimer exhibited the highest prediction confidence, indicating a significant role of A4galt heterodimerization in determining enzyme specificity toward GSLs and GPs. These findings enhance our knowledge of A4galt acceptor specificity and may contribute to a better comprehension of pathomechanisms of the Shiga toxin-related diseases.