In latent interactions of dislocations, junction formation is one of the most important phenomena that contribute to the evolution of strength. In this work, the latent hardening coefficients for pure aluminum are estimated using 3D multiscale dislocation dynamics program (MDDP). Three well-known junction configurations, namely, the Hirth lock, the glissile junction, and the Lomer lock, are studied using 3D discrete dislocation dynamics simulations. The evolution of strength is discussed as a function of the resolved shear stress (RSS) and the number of junctions for the three junctions investigated. Hirth lock and Lomer lock are found to be the weakest and strongest junctions, respectively. Collinear reaction of dislocations does not form a junction but causes a higher strength than a Lomer lock. Quantitative and qualitative results are compared with those found in the literature.