Plasmodium vivax-induced malaria is one of the leading causes of morbidity and mortality in sub-tropical and tropical regions and known to infect 2.85 billion people globally. The continual rise and propagation of resistance against anti-malarial drugs is a prerequisite to identify a possible vaccine candidate for Plasmodium vivax (P. vivax). Circumsporozoite protein (CSP) is an important immunogen of malaria parasite that has conserved the CSP structure as an immune dominant B-cell epitope. In the current study, we focused on designing multi-epitope vaccines (MEVs) using various immunoinformatics tools against Pakistani based allelic variants VK210 and VK247 of P. vivax CSP (PvCSP) gene. Antigenicity, allergic potential and physicochemical parameters of both PvCSP variants were assessed for the designed MEVs and are within acceptable range suitable for post experimental investigations. The three-dimensional structures of both MEV shave been predicted ab initio, optimized, and validated by using different online servers. Structure and from residues perspectives, the MEVs are stable and are free from aggregation-prone regions. The stability of both MEVs has been improved by a disulfide engineering approach. To estimate the binding energy and stability of the MEVs, molecular docking simulation and binding free energy calculations with TLR-4 immune receptor have been conducted. The expression of both MEVs produced in Escherichia coli K12 expression system by in silico cloning was significant. Immune simulation revealed that the proposed MEVs induce strong humoral and cellular immunological responses, in addition to significant production of interleukins and cytokines. In conclusions, we believed that the MEVs proposed in current research, using combine approach of immunoinformatics, structural biology and biophysical approaches, could induce protective and effective immune responses against P. vivax and the experimental validation of our findings could contribute to the development of potential malaria vaccine.