Abstract
Kinetic models mechanistically link enzyme levels, metabolite concentrations, and allosteric regulation to metabolic reaction fluxes. This coupling allows for the quantitative elucidation of the dynamics of the evolution of metabolite concentrations and metabolic fluxes as a function of time. So far, most large-scale kinetic model parameterizations are carried out using mostly steady-state flux measurements supplemented with metabolomics and/or proteomics data when available. Even though the parameterized kinetic model can trace a temporal evolution of the system, lack of anchoring to temporal data reduces confidence in the dynamics predictions. Notably, the simulation of enzymatic cascade reactions requires a full description of the dynamics of the system as a steady-state is not applicable given that all measured metabolite concentrations vary with time. Here we describe how kinetic parameters fitted to the dynamics of single-enzyme assays remain accurate for the simulation of multi-enzyme cell-free systems. Herein, we demonstrate two extensions for the Kinetic Estimation Tool Capturing Heterogeneous datasets Using Pyomo (KETCHUP) software tool for parameterizing a kinetic model of the cell-free kinetics of formate dehydrogenase (FDH) and 2,3-butanediol dehydrogenase (BDH) through the use of time-course data across various initial conditions. An implemented extension of KETCHUP allowing for the reconciliation of measurement time-lag errors present in datasets was used to parameterize kinetic models using multiple datasets. By combining the kinetic parameters identified by the FDH and BDH assays, accurate simulation of the binary FDH-BDH system was achieved.