To read the full version of this content please select one of the options below:

Global stochastic optimization in hierarchical modeling of ligand/protein binding profiles

Alexander D. Varshavsky (Department of Statistical and Mathematical Sciences, Lilly Research Laboratories, Indianapolis, Indiana, USA)
Duane T. Birnbaum (Department of Biopharmaceutical Development, Eli Lilly and Company, Lilly Research Laboratories, Indianapolis, Indiana, USA,)
John M. Beals (Department of Biopharmaceutical Development, Eli Lilly and Company, Lilly Research Laboratories, Indianapolis, Indiana, USA, and)
Bo E. Saxberg (Department of Statistical and Mathematical Sciences, Lilly Research Laboratories, Indianapolis, Indiana, USA)

Kybernetes

ISSN: 0368-492X

Article publication date: 1 June 2000

Abstract

We present a detailed explanation of a mathematical method and numerical technique applied to solve an irregular non‐linear fitting problem that results from attempts to model the calorimetric profiles generated by the binding of phenolic ligands to the insulin hexamer. The method employed uses a non‐traditional approach to modeling data. Rather than start with a simplified model, we use a hierarchical tree of physical models with different degrees of sophistication. Starting with the model of highest dimension, we work our way to an optimum model which is of a lower dimension and is less complex. The algorithm uses two complementary techniques. First, a sensitivity analysis in the vicinity of the optimal point for each model is used to estimate errors in the parameters; that, in turn, provides the user with insight for model simplification. Second, we utilize the optimized model in the prediction of new experimental curves. The core of the method combines a strategy based on the proper split of the initial global numerical task into three locally independent subtasks, and induces a specific split in the search space. The application of three different optimization techniques (two parametric and one variational) with an alternating objective function defined in corresponding subspaces, in combination with the search along the hierarchical tree of mathematical models, enables us to overcome difficult computational problems, including over‐parametrization. We have obtained very accurate fits to a number of calorimetric curves, resulting in a quantitative description of intrinsic functional (free ligand concentration) and vector (equilibrium coefficients and enthalpies of binding) parameters. These quantitative results can now be used to improve the stability of insulin formulations. We believe that, with small modifications to the model, the method and algorithms presented in this article can be applied to other protein‐ligand systems.

Keywords

Citation

Varshavsky, A.D., Birnbaum, D.T., Beals, J.M. and Saxberg, B.E. (2000), "Global stochastic optimization in hierarchical modeling of ligand/protein binding profiles", Kybernetes, Vol. 29 No. 4, pp. 452-473. https://doi.org/10.1108/03684920010322217

Publisher

:

MCB UP Ltd

Copyright © 2000, MCB UP Limited