The addition of boundary localized kinetic and Yukawa terms to the action of a five-dimensional Standard Model would nontrivially modify the Kaluza-Klein spectra and some of the interactions among the Kaluza-Klein excitations compared to the minimal version of this model, in which these boundary terms are not present. In the minimal version of this framework, known as the universal extra dimensional model, special assumptions are made about these unknown, beyond the cutoff contributions to restrict the number of unknown parameters of the theory to be minimum. We estimate the contribution of Kaluza-Klein modes to the branching ratios of Bs(d)→μ+μ- in the framework of the nonminimal universal extra dimensional model, at one-loop level. The results have been compared to the experimental data to constrain the parameters of this model. From the measured decay branching ratio of Bs→μ+μ- (depending on the values of boundary localized parameters), the lower limit on R-1 can be as high as 800 GeV. We have briefly reviewed the bounds on nonminimal universal extra dimensional parameter space coming from electroweak precision observables. The present analysis (Bs→μ+μ-) has ruled out new regions of parameter space in comparison to the analysis of electroweak data. We have revisited the bound on R-1 in the universal extra dimensional model, which came out to be 454 GeV. This limit on R-1 in the universal extra dimensional framework is not as competitive as the limits derived from the consideration of relic density or Standard Model Higgs boson production and decay to W+W-. Unfortunately, the Bd→μ+μ- decay branching ratio would not set any significant limit on R-1 in a minimal or nonminimal universal extra dimensional model. © 2016 American Physical Society.