We established a predictive kinetic metabolic-flux model for the 21 enzymes and 24 metabolites of the monolignol biosynthetic pathway using Populus trichocarpa secondary differentiating xylem. To establish this model, a comprehensive study was performed to obtain the reaction and inhibition kinetic parameters of all 21 enzymes based on functional recombinant proteins. A total of 104 Michaelis-Menten kinetic parameters and 85 inhibition kinetic parameters were derived from these enzymes. Through mass spectrometry, we obtained the absolute quantities of all 21 pathway enzymes in the secondary differentiating xylem. This extensive experimental data set, generated from a single tissue specialized in wood formation, was used to construct the predictive kinetic metabolic-flux model to provide a comprehensive mathematical description of the monolignol biosynthetic pathway. The model was validated using experimental data from transgenic P. trichocarpa plants. The model predicts how pathway enzymes affect lignin content and composition, explains a longstanding paradox regarding the regulation of monolignol subunit ratios in lignin, and reveals novel mechanisms involved in the regulation of lignin biosynthesis. This model provides an explanation of the effects of genetic and transgenic perturbations of the monolignol biosynthetic pathway in flowering plants.