The future digital evolution poses challenges that need to be spectral and energy-efficient, as well as highly reliable and resilient. The non-orthogonal multiple access (NOMA) accomplishes massive connectivity, spectral efficiency, effective bandwidth utilization, and low latency. The proposed work involves the code domain NOMA scheme called Sparse Code Multiple Access (SCMA) which provides shaping gain through multi-dimensional constellation and the best performance in terms of bit error rate (BER). It achieves overloading of users through the non-orthogonal allocation of resources which enhances the spectral efficiency and serves more users. The shaping gain can be further improved by reducing the BER and enhancing the capacity of the channel through constellation shaping. This work employs a probabilistic-shaped (PS) constellation where each symbol is transmitted with different probabilities which achieves a reduction of average symbol power and forward error correction (FEC) through channel coding using polar codes which aid in energy efficiency. The output is two-dimensionally spread over Orthogonal Frequency Code Division Multiplexing (OFCDM) subcarriers to achieve a flexible transmission rate through a variable spreading factor. Computer simulations showed better BER performance under AWGN and Rayleigh channels with remarkable gain in SNR which paves the way for future applications in Fifth Generation (5G) beyond networks.