  In this paper, we study behaviour of Wilson loops in the boost-invariant
nonequilibrium anisotropic quark-gluon plasma produced in heavy-ion collisions
within the holographic approach. We describe the thermalization studying the
evolution of the Vaidya metric in the boost-invariant and spatially anisotropic
background. To probe the system during this process we calculate rectangular
Wilson loops oriented in different spatial directions. We find that anisotropic
effects are more visible for the Wilson loops lying in the transversal plane
unlike the Wilson loops with partially longitudinal orientation. In particular,
we observe that the Wilson loops can thermalizes first unlike to the order of
the isotropic model. We see that Wilson loops on transversal contours have the
shortest thermalization time. We also calculate the string tension and the
pseudopotential at different temperatures for the static quark-gluon plasma. We
show that the pseudopotential related to the configuration on the transversal
plane has the screened Cornell form. We also show that the jet-quenching
parameter related with the average of the light-like Wilson loop exhibits the
dependence on orientations.
