Complex Kinetic Models Predict β-Carotene Production and Reveal Flux Limitations in Recombinant <i>Saccharomyces cerevisiae</i> Strains

β-Carotene is a high-value compound with multiple commercial applications as a pigment and due to its antioxidant properties. For its industrial production, precision fermentation using engineered microorganisms has been proposed as an attractive alternative given consumer concerns and technical limitations of traditional production methods such as chemical synthesis and extraction from plants. However, the factors limiting microbial production are complex and remain poorly understood, hindering bioprocess scale-up. To tackle this limitation, we built and evaluated kinetic model ensembles of the native mevalonate and the heterologous β-carotene production pathways in recombinant <i>Saccharomyces cerevisiae</i> strains to identify bottlenecks limiting the production flux. For this task, flux and transcriptomic data from chemostat cultivations were generated and combined with literature information for simulating model structures capturing different degrees of kinetic detail and complexity within the ABC-GRASP framework. Our results showed that detailed kinetic models including both allosteric regulation and complex mechanistic descriptions (e.g., enzyme promiscuity) are necessary to explain the metabolic phenotype of recombinant strains in different conditions. Calculation of flux and concentration response coefficients of the detailed models revealed that the promiscuous CrtYB enzyme exerts the highest control over β-carotene production at different growth rates in the best producer. Simulation of various enzyme and metabolite perturbations confirmed the above result and discarded other seemingly intuitive targets for intervention, e.g., upregulation of ERG10. Overall, this work deepens our understanding about the factors limiting β-carotene production in yeast, providing mechanistic models for <i>in silico</i> metabolic prospection and rational design of genetic interventions.