The molecular states in coupled channel system of $B_s^{0}N-B\Lambda -B^{*}\Lambda -B\Sigma -B^{*}\Sigma$ are investigated within one-boson-exchange model that includes $S$-$D$ wave mixing and a tunable short-range $\delta (r)$ term. Bound, resonant, and virtual states are searched by analytically continuing the $S$ matrix in the complex energy plane. In the single-channel analysis, the $ B$ and $ B^$ systems with $J^P=1/2^{-}$ and $3/2^{-}$ are found to be attractive and form three bound states with reasonable cutoff region, corresponding to $1/2^{-}(B\Sigma )$, $1/2^{-}(B^{*}\Sigma )$, and $3/2^{-}(B^{*}\Sigma )$ . When coupled-channel dynamics is included, these states evolve into near-threshold poles that produce enhancements in the $B_s^{0}N$, $ B$, and $ B^$ invariant mass spectra. The pole below the $ B^$ threshold in the $J^P=1/2^{-}$ system depends strongly on the treatment of the short-range $\delta (r)$ term and becomes a virtual state when this contribution is removed, while the other two resonant poles and their line shapes are only moderately affected. The pole coupling and partial-width analyses indicate that the invariant mass spectrum of the $B_s^{0}N$ and $ B^{()}$ channels are the best place to hunt the predicted states. These results support the existence of near-threshold molecular pentaquark $P_{\bar{b}s}$ and provide useful guidance in future experimental searches.