The two-dimensional discrete dislocation dynamics (2D DD) method, consisting of parallel straight edge dislocations gliding on independent slip systems in a plane strain model of a crystal, is often used to study complicated boundary value problems in crystal plasticity. However, the absence of truly three dimensional mechanisms such as junction formation means that forest hardening cannot be modeled, unless additional so-called '2.5D' constitutive rules are prescribed for short-range dislocation interactions. Here, results from three dimensional dislocation dynamics (3D DD) simulations in an FCC material are used to define new constitutive rules for short-range interactions and junction formation between dislocations on intersecting slip systems in 2D. The mutual strengthening effect of junctions on preexisting obstacles, such as precipitates or grain boundaries, is also accounted for in the model. The new '2.5D' DD model, with no arbitrary adjustable parameters beyond those obtained from lower scale simulation methods, is shown to predict athermal hardening rates, differences in flow behavior for single and multiple slip, and latent hardening ratios. All these phenomena are well-established in the plasticity of crystals and quantitative results predicted by the model are in good agreement with experimental observations.