In this contribution, the variational problem for the Kirchhoff plate based on the modified strain gradient theory (MSGT) is derived, and the Euler-Lagrange equations governing the equation of motion are obtained. The Galerkin-type weak form, upon which the finite element method is constructed, is derived from the variational problem. The shape functions which satisfy the governing homogeneous partial differential equation are derived as extensions of Adini-Clough-Melosh (ACM) and Bogner-Fox-Schmit (BFS) plate element formulations by introducing additional curvature degrees of freedom (DOF) on each node. Based on the proposed set of shape functions, 20-, 24-, 28- and 32- DOF modified strain gradient theory-based higher-order Kirchhoff microplate element are proposed. The performance of the elements are demonstrated in terms of various tests and representative boundary value problems. Length scale parameters for gold are also proposed based on experiments reported in literature.