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

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. Ca2+ 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.

Diastolic chamber properties of the left ventricle assessed by global fitting of pressure-volume data: improving the gold standard of diastolic function

“…Global fitting methods have demonstrated increased robustness in a number of biological processes (7,17), but had never been used to analyze LV chamber properties. With this method, we have been able to assess for the first time the full active and passive mechanical properties in the intact heart (1, 18, 24) using conventional PV data.…”

“…and plotted for each data set for the first beat in each run, adding the term for P0, ‫ץ‬P/(‫ץ‬P0/P0). Resolution matrices were also computed to rule out interdependencies among different parameters (17). Average processing time for a 10-beat data set was 1 min using a 2-GHz, quad-core i5 Intel workstation.…”

Thin Filament Structure and the Steric Blocking Model

“…Thus, local troponin-and myosininduced shifts of semirigid tropomyosin on the actin filament can be expected to be propagated along the tropomyosin cable with limited decrement, producing expected filament cooperative effects. The integrated behavior of thin filament components has been modeled by treating overlapping tropomyosins as continuous flexible chains or other representations (e.g., see 59,139,150,151,194,195,218).…”

Two-Layer Elastographic 3-D Traction Force Microscopy

“…To solve the overdetermined problem, we followed a least-squares approach using an iterative Levenberg-Marquardt algorithm that minimized a global cost function to calculate σ . This iterative algorithm is similar to those employed in our previous studies303132 and is summarized in Fig. 1b.…”