This research investigates the intricate dynamics of dissipative non-slip magnetohydrodynamic (MHD) nanofluid flow, characterized by variable viscosity and thermal conductivity, under the influence of an Arrhenius chemical reaction. The inclusion of the Arrhenius chemical reaction adds complexity through heat generation or absorption, impacting temperature and concentration gradients. The study is motivated by the extensive applications of nanofluids in engineering and industrial processes, where precise control of heat and mass transfer is critical. We develop a comprehensive mathematical model that incorporates the variable properties of the nanofluid, the effects of the Lorentz force due to the applied magnetic field, and the temperature-dependent reaction rates dictated by the Arrhenius equation. The formulated governing equations were non-dimensionalised to identify the flow governing parameters. Finite Element Method (FEM), grid generation, solution algorithms, and post-processing to analyse velocity, temperature, and concentration distributions were used to obtain the numerical methods to solve fluid flow problems based on the Navier-Stokes equations, involving concepts of discretization. pdsolve subpackage in Maple 2023 was used to numerically solve PDEs with specific initial and boundary conditions, incorporating the plot and display commands for graphical analysis, and the results are presented and discussed. The findings reveal that the interplay between parameters like Hartmann number, Darcy parameter, and heat generation or absorption profoundly influences flow behaviour and thermal characteristics. The reactivity parameter is crucial, dictating the rate of chemical reactions and affecting system dynamics. This research enhances understanding of the interdependencies among fluid properties, chemical reactions, and external parameters in nanofluid flows.