Presented in this current study is the numerical analysis of magnetohydrodynamics Williamson nanofluid flow over an exponentially stretching surface. The most important aspect of the investigation is that the effects of the magnetic field, chemical reaction and thermal radiation in the fluid flow are taken into account. The partial differential equations governing the present Williamson nanofluid flow problem were observed to be highly nonlinear and coupled. Suitable similarity transformations were used to transmute the coupled system of nonlinear partial differential equations governing the fluid flow into a linear system. The linear system was solved numerically using the spectral quasi-linearization method. The MATLAB bvp4c numerical technique and a comparison with existing results for the skin friction coefficient were used to confirm the appropriateness of the method in solving the current problem. The influence of some pertinent physical parameters on the fluid’s velocity, temperature and concentration profiles were displayed graphically. The effects of all the physical parameters on the skin friction coefficient, Nusselt number and Sherwood number were portrayed in a tabular form. It was noted that enhancing the thermal radiation parameter reduces the fluid’s temperature, Nusselt number and the skin friction coefficient, while the Sherwood number is improved.