The Madden–Julian oscillation (MJO) typically contains several superclusters and numerous embedded mesoscale convective systems (MCSs). It is hypothesized here that the poorly simulated MJOs in current coarse-resolution global climate models (GCMs) is related to the inadequate treatment of unresolved MCSs. So its parameterization should provide the missing collective effects of MCSs. However, a satisfactory understanding of the upscale impact of MCSs on the MJO is still lacking. A simple two-dimensional multicloud model is used as an idealized GCM with clear deficiencies. Eddy transfer of momentum and temperature by the MCSs, predicted by the mesoscale equatorial synoptic dynamics (MESD) model, is added to this idealized GCM. The upscale impact of westward-moving MCSs promotes eastward propagation of the MJO analog, consistent with the theoretical prediction of the MESD model. Furthermore, the upscale impact of upshear-moving MCSs significantly intensifies the westerly wind burst because of two-way feedback between easterly vertical shear and eddy momentum transfer with low-level eastward momentum forcing. Finally, a basic parameterization of the upscale impact of upshear-moving MCSs modulated by deep heating excess and vertical shear strength significantly improves key features of the MJO analog in the idealized GCM with clear deficiencies. A three-way interaction mechanism between the MJO analog, parameterized upscale impact of MCSs, and background vertical shear is identified.