Large eddy simulation, using a dynamic Smagorinsky sub-grid scale, is used for the prediction of flow structures around two identical spheres fixed side-by-side at a subcritical flow regime with the Reynolds number equal to 5,000. The QUICK discretization method is applied to discretize the convection terms of the Navier-Stokes equation by means of the finite volume approach. This work focuses mainly on the wake structures downstream of the two interactive spheres located at three various dimensionless separation distances between spheres such as G/D = 1.50, 2.00, and 3.00. The obtained results revealed that the interaction between wakes affects the flow structures downstream of spheres. The rate of this interaction is strongly altered as a function of separation distance. On the other hand, some flow data such as mean drag and mean lift coefficients are affected due to the wake interactions. Finally, examination of the Reynolds stress variation along the different lateral axis, L/D, revealed that the nozzle effect does not play a significant role on the turbulence characteristics beyond the G/D = 3.00.