Massive connections and high data rate services are key players in 5G ecosystem and beyond. To satisfy the requirements of these types of services, non orthogonal multiple Access (NOMA) and 3-dimensional beam-forming (3DBF) can be exploited. In this paper, we devise a novel radio resource allocation and 3D multiple input single output (MISO) BF algorithm in NOMA-based heterogeneous networks (HetNets) at which our main aim is to maximize the profit of mobile virtual network operators (MVNOs). To this end, we consider multiple infrastructure providers (InPs) and MVNOs serving multiple users. Each InP has multiple access points as base stations (BSs) with specified spectrum and multi-beam antenna array in each transmitter that share its spectrum with MVNO's users by employing NOMA. To realize this, we formulate a novel optimization problem at which the main aim is to maximize the revenue of MVNOs, subject to resource limitations and quality of service (QoS) constraints. Since our proposed optimization problem is non-convex and mathematically intractable, we transform it into a convex one by introducing a new optimization variable and converting the variables with adopting successive convex approximation. More importantly, the proposed solution is assessed and compared with the alternative search method and the optimal solution that is obtained with adopting the exhaustive search method. In addition, it is studied from the computational complexity, convergence, and performance perspective. Our simulation results demonstrate that NOMA-3DBF has better performance and increases system throughput and MVNO's revenue compared to orthogonal multiple access with 2DBF by approximately 64%. Especially, by exploiting 3DBF the MVNO's revenue is improved nearly 27% in contrast to 2DBF in high order of antennas.