In this article, we obtain a new set of nuclear Parton distribution functions (nuclear PDFs) at next-to-leading order and next-to-next-to-leading order accuracy in perturbative QCD. The common nuclear deep-inelastic scattering (DIS) data analyzed in our study are complemented by the available charged-current neutrino DIS data with nuclear targets and data from Drell-Yan cross-section measurements for several nuclear targets. In addition, the most recent DIS data from the Jefferson Lab CLAS and Hall C experiments are also added to our data sample. For these specific datasets, we consider the impact of target mass corrections and higher twist effects which are expected to be important in the region of large x and intermediate-to-low Q2. Our analysis is based on a publicly available open-source tool, APFEL, which has been modified to be applicable for our analysis of nuclear PDFs. Heavy-quark contributions to nuclear DIS are considered within the framework of the FONLL general-mass variable-flavor-number scheme. The most recent CT18 PDFs are used as baseline proton PDFs. The uncertainties of nuclear PDFs are determined using the standard Hessian approach. The results of our global QCD analysis are compared with existing nuclear PDF sets and with the fitted cross-sections, for which our set of nuclear PDFs provides a very good description.