The organization of Golgi-resident proteins is crucial for sorting molecules within the secretory pathway and regulating posttranslational modifications. However, evaluating changes to Golgi organization can be challenging, often requiring extensive experimental investigations. Here, we propose a systems biology approach in which changes to Golgi-resident protein sorting and localization can be deduced using cellular N-glycan profiles as the only experimental input.The approach detailed here utilizes the influence of Golgi organization on N-glycan biosynthesis to investigate the mechanisms involved in establishing and maintaining Golgi organization. While N-glycosylation is carried out in a non-template-driven manner, the distribution of N-glycan biosynthetic enzymes within the Golgi ensures this process is not completely random. Therefore, changes to N-glycan profiles provide clues into how altered cell phenotypes affect the sorting and localization of Golgi-resident proteins. Here, we generate a stochastic simulation of N-glycan biosynthesis to produce a simulated glycan profile similar to that obtained experimentally and then combine this with Bayesian fitting to enable inference of changes in enzyme amounts and localizations. Alterations to Golgi organization are evaluated by calculating how the fitted enzyme parameters shift when moving from simulating the glycan profile of one cellular state (e.g., a wild type) to an altered cellular state (e.g., a mutant). Our approach illustrates how an iterative combination of mathematical systems biology and minimal experimental cell biology can be utilized to maximally integrate biological knowledge to gain insightful knowledge of the underlying mechanisms in a manner inaccessible to either alone.