[en] Three general cases of dynamical interacting dark energy models ($\mathcal{D}$-class) are investigated in the context of Brans-Dicke cosmology. Some of important cosmological quantities are calculated for every cases as a function of redshift parameter. The most important part of this paper deals with fitting models with two different expansion history: ($\text{SNIa}+{\text{BAO}}_A+{\mathit{O}\mathit{m}\mathit{h}}^2$ and $\text{SNIa}+{\text{BAO}}_A+H(z)$) and with two different sets of data for Hubble parameter. This provides a remarkable feature to could analytically see the effects of each analyses and each data sets on final results. The best fitted values of parameters according to these analyses and data points, ${\mathrm{\chi <!--\; ??\; -->}}_{\mathit{t}\mathit{o}\mathit{t}}^2/\mathit{d}\mathit{o}\mathit{f}$, AIC and BIC are reported. By these diagnostic tools we found that some of these models have no chance against $\mathrm{\Lambda <!--\; ??\; -->}$CDM, even without need to study the structure formation, and could be ruled out. While some (e.g. $\mathcal{B}\mathcal{D}\mathcal{-<!--\; ???\; -->}\mathcal{D}C2$ and $\mathcal{B}\mathcal{D}\mathcal{-<!--\; ???\; -->}\mathcal{D}{A}^{\ast <!--\; ???\; -->}$) render the best fit quality, i.e. the value of AIC and BIC and figures shows that they fit perfectly with overall data and reveals a strong evidence in favor of these two models against $\mathrm{\Lambda <!--\; ??\; -->}$CDM.