The program of matrix product states on tensor powers ({mathcal {A}}^{otimes {mathbb {Z}}}) of (C^*)-algebras is carried under the assumption that ({mathcal {A}}) is an arbitrary nuclear C*-algebra. For any shift invariant state (omega ), we demonstrate the existence of an order kernel ideal ({mathcal {K}}_omega ), whose quotient action reduces and factorizes the initial data (({mathcal {A}}^{otimes {mathbb {Z}}}, omega )) to the tuple (({mathcal {A}},{mathcal {B}}_omega = {mathcal {A}}^{otimes {mathbb {N}}^times }/{mathcal {K}}_omega , {mathbb {E}}_omega : text{AA }otimes {mathcal {B}}_omega rightarrow {mathcal {B}}_omega , {bar{omega }}: {mathcal {B}}_omega rightarrow {mathbb {C}})), where ({mathcal {B}}_omega ) is an operator system and ({mathbb {E}}_omega ) and ({bar{omega }}) are unital and completely positive maps. Reciprocally, given a (input) tuple (({mathcal {A}},{mathcal {S}},{mathbb {E}},phi )) that shares similar attributes, we supply an algorithm that produces a shift-invariant state on ({mathcal {A}}^{otimes {mathbb {Z}}}). We give sufficient conditions in which the so constructed states are ergodic and they reduce back to their input data. As examples, we formulate the input data that produces AKLT-type states, this time in the context of infinite dimensional site algebras ({mathcal {A}}), such as the (C^*)-algebras of discrete amenable groups.