Electrospray ionization (ESI) is commonly used in macromolecular mass spectrometry, yet the dynamics of macromolecules in ESI droplets are not well understood. In this study, a Monte Carlo based model was developed, which can predict the efficiency of electrospray ionization for macromolecules, i.e., the number of macromolecular ions produced per macromolecules electrosprayed. The model takes into account ESI droplet evaporation, macromolecular diffusion within the droplet, droplet fissions, and the statistical nature of the ESI process. Two idealized representations of macromolecular analytes were developed, describing cluster prone, droplet surface inactive macromolecules and droplet surface active macromolecules, respectively. It was found that surface active macromolecules are preferentially ionized over surface inactive cluster prone macromolecules when the initial droplet size is large and the analyte concentration in solution is high. Simulations showed that ESI efficiency decreases with increasing initial droplet size and analyte molecular weight, and is influenced by analyte surface activity, the properties of the solvent, and the variance of the droplet size distribution. Model predictions are qualitatively supported by experimental measurements of macromolecular electrospray ionization made previously. Overall, this study demonstrates the potential capabilities of Monte Carlo based ESI models. Future developments in such models will allow for more accurate predictions of macromolecular ESI intensity.