Based on the finite deformation equation of protein bubble with viscoelastic film, and taking into account the Bjerknes effect that one oscillating bubble acts on the ather in Bingham liquid, a nonlinear equation describing the vibration of two equal-radius protein bubbles in Bingham fluid is built. The numerical simulation is used to solve the above nonlinear equation. The results show that, increasing the plastic viscosity of Bingham liquid the protein bubble wall will vibrate with higher decrement velocity of amplitude, and with lower frequency; shortening the distance between bubbles, the bubble wall will vibrate with a higher frequency and a higher decrement velocity of amplitude. Furthermore, the smaller the bubble size is, the higher the increment of frequency and decrement velocity of amplitude are. The two bubbles in Bingham fluid will vibrate with higher frequency and higher decrement velocity of amplitude than that of a single bubble.