Abstract:The estimation of chemical kinetic rate constants for any non-trivial model is complex due to the nonlinear effects of second order chemical reactions. We developed an algorithm to accomplish this goal based on the Damped Least Squares (DLS) inversion method and then tested the effectiveness of this method on the McKillop-Geeves (MG) model of thin filament regulation. The kinetics of MG model is defined by a set of nonlinear ordinary differential equations (ODEs) that predict the evolution of troponin-tropomyosin-actin and actin-myosin states. The values of the rate constants are estimated by integrating these ODEs numerically and fitting them to a series of stopped-flow pyrene fluorescence transients of myosin-S1 fragment binding to regulated actin in solution. The accuracy and robustness of the estimated rate constants are evaluated for DLS and two other methods, namely quasi-Newton (QN) and simulated annealing (SA). The comparison of these methods revealed that SA provides the best estimates of the model parameters because of its global optimization scheme. However it converges slowly and does quantify the uniqueness of the estimated parameters. On the other hand the QN method converges rapidly but only if the initial guess of the parameters is close to the optimum values, otherwise it diverges. Overall, the DLS method proves to be the most convenient method. It converges fast and was able to provide excellent estimates of kinetic parameters. Furthermore, DLS provides the model resolution matrix, which quantifies the interdependence of model parameters thereby evaluating the uniqueness of their estimated values. This property is essential for estimating of the dependence of the model parameters on experimental conditions (e.g. Ca(2+) concentration) when it is assessed from noisy experimental data such as pyrene fluorescence from stopped-flow transients. The advantages of the DLS method observed in this study should be further examined in other physicochemical systems to firmly establish the observed effectiveness of DSL vs. the other parameter estimation methods.

On the reliable estimation of sequential Monod kinetic parameters

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

Bayesian Inference for Kinetic Models of Biotransformation Using A Generalized Rate Equation

Comparison of various estimation methods for the parameters of Michaelis-Menten equation based on in vitro elimination kinetic simulation data

Parameter Estimation in Biokinetic Degradation Models in Wastewater Treatment - A Novel Approach Relevant for Micro-pollutant Removal

Data-driven Bayesian estimation of Monod kinetics

An approach for estimating microbial growth and chlorinated ethenes biotransformation in groundwater based on discrete observations

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Parameter Estimation for Reactive Chromatography Model by Bayesian Inference and Parallel Sequential Monte Carlo

Frequency Adaptive Metadynamics for the Calculation of Rare-Event Kinetics.

Studies in Parameter Estimation for Analytical Transport Equations

Estimation of Kinetic Parameters for Chemical Reactions Based on Monte Carlo Simulation

Uncertainty Quantification of Linear Scaling, Machine Learning, and DFT Derived Thermodynamics for the Catalytic Partial Oxidation of Methane on Rhodium

Determination of growth kinetics of microorganisms linked with 1,4-dioxane degradation in a consortium based on two improved methods

Acceleration and sensitivity analysis of lattice kinetic Monte Carlo simulations using parallel processing and rate constant rescaling.

Resolution and uniqueness of estimated parameters of a model of thin filament regulation in solution

Analytical methods for fitting integrated rate equations. A discontinuous assay

Towards Kinetic Modeling of Global Metabolic Networks with Incomplete Experimental Input on Kinetic Parameters

Omics data for sampling thermodynamically feasible kinetic models

Improvement of Numerical Simulation of Bimolecular Reactive Transport Using Time-Dependent Parameters

Kinetic profiling of metabolic specialists demonstrates stability and consistency of in vivo enzyme turnover numbers