In this paper, we investigate a (3+1)-dimensional modified Zakharov-Kuznetsov equation, which describes the nonlinear plasma-acoustic waves in a multicomponent magnetised plasma. With the aid of the Hirota method and symbolic computation, bilinear forms and one-, two- and three-soliton solutions are derived. The characteristics and interaction of the solitons are discussed graphically. We present the effects on the soliton’s amplitude by the nonlinear coefficients which are related to the ratio of the positive-ion mass to negative-ion mass, number densities, initial densities of the lower- and higher-temperature electrons and ratio of the lower temperature to the higher temperature for electrons, as well as by the dispersion coefficient, which is related to the ratio of the positive-ion mass to the negative-ion mass and number densities. Moreover, using the Lie symmetry group theory, we derive the Lie point symmetry generators and the corresponding symmetry reductions, through which certain analytic solutions are obtained via the power series expansion method and the ( G ′/ G ) expansion method. We demonstrate that such an equation is strictly self-adjoint, and the conservation laws associated with the Lie point symmetry generators are derived.