Let \(n\ge 3\) and let \(\mathscr {M}_n(\mathbb {C})\) denote the algebra of all \(n\times n\) complex matrices. For any matrix \(A\in \mathscr {M}_n(\mathbb {C})\) , let \(A^*\) and \(A^{\dagger }\) denote the adjoint and Moore–Penrose inverse of A, respectively. For a positive scalar \(\alpha \) , the Bourhim–Mbekhta transformation \(\Theta _\alpha \) is defined on \(\mathscr {M}_n(\mathbb {C})\) by \( \Theta _\alpha (A):= \frac{1}{\alpha +1}\left( \alpha A + A^{*\dagger }\right) , \quad (A\in \mathscr {M}_n(\mathbb {C})). \) In this paper, we show that if \(\phi \) is a map on \(\mathscr {M}_n(\mathbb {C})\) , then \(\Theta _\alpha \!\left( \phi (A)\phi (B)\right) \) and \(\Theta _\alpha (AB)\) are unitarily similar for all \(A,B\in \mathscr {M}_n(\mathbb {C})\) if and only if there exist a unitary matrix \(U\in \mathscr {M}_n(\mathbb {C})\) and a scalar \(\epsilon \) with \(\epsilon ^2=1\) , such that \(\phi (A)=\epsilon UAU^*\) for all \(A\in \mathscr {M}_n(\mathbb {C})\) . We also characterize all maps \(\phi \) on \(\mathscr {M}_n(\mathbb {C})\) for which \(\Theta _\alpha (ABA)\) and \(\Theta _\alpha (\phi (A)\phi (B)\phi (A))\) are unitarily similar for every \(A,B\in \mathscr {M}_n(\mathbb {C})\) .