Arabinogalactan proteins are very abundant, heavily glycosylated plant cell wall proteins. They are intensively studied because of their crucial role in plant development as well as their function in plant defence. Research of these biomacromolecules is complicated by the lack of tools for their analysis and characterisation due to their extreme heterogeneity. One of the few available tools for detection, isolation, characterisation, and functional studies of arabinogalactan proteins is Yariv reagents. Yariv reagent is a synthetic aromatic glycoconjugate originally prepared as an antigen for immunization. Later, it was found that this compound can precipitate arabinogalactan proteins, namely, their ß-D-(1→3)-galactan structures. Even though this compound has been intensively used for decades, the structural basis of arabinogalactan protein precipitation by Yariv is not known. Multiple biophysical studies have been published, but none of them attempted to elucidate the three-dimensional structure of the Yariv-galactan complex. Here we use a series of molecular dynamics simulations of systems containing one or multiple molecules of ß-D-galactosyl Yariv reagent with or without oligo ß-D-(1→3)-galactan to predict the structure of the complex. According to our model of Yariv-galactan complexes, Yariv reagent forms stacked oligomers stabilized by π-π and CH/π interactions. These oligomers may contain irregularities. Galactan structures crosslink these Yariv oligomers. The results were compared with studies in literature.