A specific modification of Newtonian dynamics known as MOND has been shown to reproduce the dynamics of most astrophysical systems at different scales without invoking non-baryonic dark matter (DM). There is, however, a long-standing unsolved problem when MOND is applied to rich clusters of galaxies in the form of a deficit (by a factor around two) of predicted dynamical mass derived from the virial theorem with respect to observations. In this article we approach the virial theorem using the velocity dispersion of cluster members along the line of sight rather than using the cluster temperature from X-ray data and hydrostatic equilibrium. Analytical calculations of the virial theorem in clusters for Newtonian gravity+DM and MOND are develo...