A coculture of Syntrophobacter fumaroxidans and Methanospirillum hungatei was modeled using four biokinetic models, which only differed by the functions used to describe the growth yields (dynamic or constant) and the hydrogen inhibition function (noncompetitive or based on thermodynamics). First, a batch experiment was used to train the model and analyze the fitted parameters. Two fitting procedures were followed by minimizing the error on different indicators. Second, a chemostat experiment was used as a test data set to assess the predictive power of the models. Overall, the four models were able to accurately fit the train data set following both fitting procedures. However, some parameters fitted with the ADM1-like model differed significantly from values reported in the literature and were dependent on the fitting procedure. When applied to the test data set it systematically resulted in positive Gibbs free energy changes values for propionate oxidation, in contradiction with the second law of thermodynamics. On the opposite, the parameters fitted with model including both a thermodynamic-based inhibition function and a dynamic computation of growth yields were more consistent with values reported in the literature and repeatable whatever the fitting procedure. The results highlight the potential of implementing thermodynamic-based functions in biokinetic models.