Ventilation-induced tumour motion remains a challenge for the accuracy of proton therapy treatments in lung patients. We investigated the feasibility of using a 4D virtual CT (4D-vCT) approach based on deformable image registration (DIR) and motion-aware 4D CBCT reconstruction (MA-ROOSTER) to enable accurate daily proton dose calculation using a gantry-mounted CBCT scanner tailored to proton therapy. Methods: Ventilation correlated data of 10 breathing phases were acquired from a porcine exvivo functional lung phantom using CT and CBCT. 4D-vCTs were generated by (1) DIR of the mid-position 4D-CT to the mid-position 4D-CBCT (reconstructed with the MA-ROOSTER) using a diffeomorphic Morphons algorithm and (2) subsequent propagation of the obtained mid-position vCT to the individual 4D-CBCT phases. Proton therapy treatment planning was performed to evaluate dose calculation accuracy of the 4D-vCTs. A robust treatment plan delivering a nominal dose of 60 Gy was generated on the average intensity image of the 4D-CT for an approximated internal target volume (ITV). Dose distributions were then recalculated on individual phases of the 4D-CT and the 4D-vCT based on the optimized plan. Dose accumulation was performed for 4D-vCT and 4D-CT using DIR of each phase to the mid position, which was chosen as reference. Dose based on the 4D-vCT was then evaluated against the dose calculated on 4D-CT both, phase-by-phase as well as accumulated, by comparing dose volume histogram (DVH) values (Dmean, D2%, D98%, D95%) for the ITV, and by a 3D-gamma index analysis (global, 3% / 3mm, 5 Gy, 20 Gy and 30 Gy dose thresholds). Results: Good agreement was found between the 4D-CT and 4D-vCT-based ITV-DVH curves. The relative differences ((CT-vCT)/CT) between accumulated values of ITV Dmean, D2%, D95% and D98% for the 4D-CT and 4D-vCT-based dose distributions were-0.2%, 0.0%,-0.1% and-0.1%, respectively. Phase specific values varied between-0.5%-0.2%,-0.2%-0.5%,-3.5%-1.5% and-5.7%-2.3%. The relative difference of accumulated Dmean over the lungs was 2.3% and Dmean for the phases varied between-5.4% and 5.8%. The gamma pass-rates with 5 Gy, 20 Gy and 30 Gy thresholds for the accumulated doses were 96.7.%, 99.6% and 99.9%, respectively. Phase-by-phase comparison yielded pass-rates between 86%-97%, 88%-98% and 94%-100%. 2 Conclusions: Feasibility of the suggested 4D-vCT workflow using proton therapy specific imaging equipment was shown. Results indicate the potential of the method to be applied for daily 4D proton dose estimation.