The production of recombinant proteins is frequently enhanced at the levels of transcription, codon usage, protein folding and secretion. Overproduction of heterologous proteins, however, also directly affects the primary metabolism of the producing cells. By incorporation of the production of a heterologous protein into a genome scale metabolic model of the yeast Pichia pastoris, the effects of overproduction were simulated and gene targets for deletion or overexpression for enhanced productivity were predicted. Overexpression targets were localized in the pentose phosphate pathway and the TCA cycle, while knockout targets were found in several branch points of glycolysis. Five out of 9 tested targets led to an enhanced production of cytosolic human superoxide dismutase (hSOD). Expression of bacterial beta-glucuronidase could be enhanced as well by most of the same genetic modifications. Beneficial mutations were mainly related to reduction of the NADP/H pool and the deletion of fermentative pathways. Overexpression of the hSOD gene itself had a strong impact on intracellular fluxes, most of which changed in the same direction as predicted by the model. In vivo fluxes changed in the same direction as predicted to improve hSOD production. Genome scale metabolic modeling is shown to predict overexpression and deletion mutants which enhance recombinant protein production with high accuracy.