The rapidly growing market of monoclonal antibodies (mAbs) within the biopharmaceutical industry has incentivised numerous works on the design of more efficient production processes. Protein A affinity chromatography is regarded as one of the best processes for the capture of mAbs. Although the screening of Protein A resins has been previously examined, process flexibility has not been considered to date. Examining performance alongside flexibility is crucial for the design of processes that can handle disturbances arising from the feed stream. In this work, we present a model-based approach for the identification of design spaces, enhanced by machine learning. We demonstrate its capabilities on the design of a Protein A chromatography unit, screening five industrially relevant resins. The computational results favourably compare to experimental data and a resin performance comparison is presented. An improvement on the computational time by a factor of 300,000 is achieved using the machine learning aided methodology. This allowed for the identification of 5,120 different design spaces in only 19 h.