This paper deals with the numerical simulation of transient magnetohydrodynamics natural convective boundary layer flow of a nanofluid over an inclined plate. In the modeling of nanofluids, dynamic effects including the Brownian motion and thermophoresis are taken into account. Numerical solutions have been computed via the Galerkin-finite element method. The effects of angle of inclination, buoyancy-ratio parameter, Brownian motion, thermophoresis and magnetic field are taken into account and controlled by non-dimensional parameters. To compute the rate of convergence and error of the computed numerical solution, the double mesh principle is used. Similarity solutions are calculated and presented graphically for non-dimensional velocity, temperature, local rate of heat and mass transfer with pertinent parameters. The modified Nusselt number decreases with increasing inclination angle, buoyancy-ratio parameter, Brownian motion and thermophoresis parameter, whereas it increases with increasing Prandtl number. Validation of the results is achieved with earlier results for forced convective flow and non-magnetic studies. Such problems have several applications in engineering and petroleum industries such as electroplating, chemical processing of heavy metals and solar water heaters. External magnetic fields play an important role in electrical power generation, inclination/acceleration sensors, fine-tuning of the final materials to industrial specification because of their controlling behaviour on the flow characteristics of nanofluids.