The vast number of expression hosts available for recombinant protein production have a variety of advantages and disadvantages; none, however, is globally optimal and host selection is frequently a compromise. Strain development requires a holistic approach, which systems biology can supply by delineating experimental data sets with computational modelling. Here, we review recent advances in computational models, in parallel with an expansion of the molecular toolbox, in the pursuit of optimal host strains for industrial protein production.