The fusion of two Mirkovic-Vilonen cycles is a degeneration of their product, defined using the Beilinson-Drinfeld Grassmannian. In this paper, we put in place a conceptually elementary approach to computing this product in type $A$. We do so by transferring the problem to a fusion of generalized orbital varieties using the Mirkovic-Vybornov isomorphism. As an application, we explicitly compute all cluster exchange relations in the coordinate ring of the upper-triangular subgroup of $GL_4$, confirming that all the cluster variables are contained in the Mirkovic-Vilonen basis.