This article presents an optimal estimator for discrete-time systems disturbed by output white noise, where the proposed algorithm identifies the parameters of a Multiple Input Single Output LPV State Space model. This is an LPV version of a class of algorithms proposed elsewhere for identifying LTI systems. These algorithms use the matchable observable linear identification parameterization that leads to an LTI predictor in a linear regression form, where the ouput prediction is a linear function of the unknown parameters. With a proper choice of the predictor parameters, the optimal prediction error estimator can be approximated. In a previous work, an LPV version of this method, that also used an LTI predictor, was proposed; this LTI predictor was in a linear regression form enablin, in this way, the model estimation to be handled by a Least-Squares Support Vector Machine approach, where the kernel functions had to be filtered by an LTI 2D-system with the predictor dynamics. As a result, it can never approximate an optimal LPV predictor which is essential for an optimal prediction error LPV estimator. In this work, both the unknown parameters and the state-matrix of the output predictor are described as a linear combination of a finite number of basis functions of the scheduling signal; the LPV predictor is derived and it is shown to be also in the regression form, allowing the unknown parameters to be estimated by a simple linear least squares method. Due to the LPV nature of the predictor, a proper choice of its parameters can lead to the formulation of an optimal prediction error LPV estimator. Simulated examples are used to assess the effectiveness of the algorithm. In future work, optimal prediction error estimators will be derived for more general disturbances and the LPV predictor will be used in the Least-Squares Support Vector Machine approach.