The control of spoilage yeasts in wines is crucial to avoid organoleptic deviations in wine production. Traditionally, sulfur dioxide (SO2) was used to control them; nevertheless, SO2 influence on human health and its use is criticized. Biocontrol emerges as an alternative in wine pre-fermentation, but there is limited development in its applicability. Managing kinetics is relevant in the microbial interaction process. pH was identified as a factor affecting the interaction kinetics of Wickerhamomyces anomalus killer biocontrol on Zygosaccharomyces rouxii. Mathematical modeling allows insight into offline parameters and the influence of physicochemical factors in the environment. Incorporating submodels that explain manipulable factors (pH), the process can be optimized to achieve the best-desired outcomes. The aim of this study was to model and optimize, using a constant and a variable pH profile, the interaction of killer biocontrol W. anomalus vs. Z. rouxii to reduce the spoilage population in pre-fermentation. The evaluated biocontrol was W. anomalus against the spoilage yeast Z. rouxii in wines. The kinetic interactions of yeasts were studied at different pH levels maintained constant over time. The improved Ramón-Portugal model was adopted using the AMIGO2 toolbox for Matlab. A static optimization of a constant pH profile was performed using the Monte Carlo method, and a dynamic optimization was carried out using a method based on Fourier series and orthogonal polynomials. The model fit with an adjusted R2 of 0.76. Parametric analyses were consistent with the model behavior. Variable vs. constant optimization achieved a lower initial spoilage population peak (99% less) and reached a lower final population (99% less) in a reduced time (100 vs. 140 h). These findings reveal that control with a variable profile would allow an early sequential inoculation of S. cerevisiae. The models explained parameters that are difficult to quantify, such as general inhibitor concentration and toxin concentration. Also, the models indicate higher biocontrol efficiency parameters, such as toxin emission or sensitivity to it, and lower fitness of the contaminant, at pH levels above 3.7 during biocontrol. From a technological standpoint, the study highlights the importance of handling variable profiles in the controller associated with the pH management actuators in the process without incurring additional costs.