Modulated electron radiotherapy (MERT) may potentially be an effective modality for the treatment of shallow tumors, but dose calculation accuracy and delivery efficiency challenges remain. The purpose of this work is to investigate the dose accuracy of MERT delivery using a photon multileaf collimator (pMLC) on a Siemens Primus accelerator. A Monte Carlo (MC)-based inverse treatment planning system was developed for the 3D treatment planning process. Phase space data of 6, 9, 12 and 15 MeV electron beams were commissioned and used as the input source for MC dose calculations. A treatment plan was performed based on the 3D CT data of a heterogeneous 'breast phantom' that mimics a breast cancer patient, and delivered with 22 segments, each associated with a particular energy and Monitor Unit value. Film and ion chamber dosimetry was carefully performed for the conversion from measurement reading to dose, and the results were employed for plan verification using the heterogeneous breast phantom and a solid water phantom. Dose comparisons between measurements and calculations showed agreement within 2% or 1 mm. We conclude that our in-house MC treatment planning system is capable of performing treatment planning and accurate dose calculations for MERT using the pMLC to deliver radiation therapy to the intact breast.