Finite and infinite matrix product states for Gutzwiller projected mean-field wavefunctions and their applications to quantum spin systems