Lignin is a heteropolymer that is thought to form in the cell wall by combinatorial radical coupling of monolignols. Here, we present a simulation model of in vitro lignin polymerization, based on the combinatorial coupling theory, which allows us to predict the reaction conditions controlling the primary structure of lignin polymers. Our model predicts two controlling factors for the b-O-4 content of syringyl-guaiacyl lignins: the supply rate of monolignols and the relative amount of supplied sinapyl alcohol monomers. We have analyzed the in silico degradability of the resulting lignin polymers by cutting the resulting lignin polymers at b-O-4 bonds. These are cleaved in analytical methods used to study lignin composition, namely thioacidolysis and derivatization followed by reductive cleavage, under pulping conditions, and in some lignocellulosic biomass pretreatments.