Structural robustness of biochemical network models-with application to the oscillatory metabolism of activated neutrophils

Sensitivity of biochemical network models to uncertainties in the model structure, with a focus on autonomously oscillating systems, is addressed. Structural robustness, as defined here, concerns the sensitivity of the model predictions with respect to changes in the specific interactions between the network components and encompass, for instance, uncertain kinetic models, neglected intermediate reaction steps and unmodelled transport phenomena. Traditional parametric sensitivity analysis does not address such structural uncertainties and should therefore be combined with analysis of structural robustness. Here a method for quantifying the structural robustness of models for systems displaying sustained oscillations is proposed. The method adopts concepts from robust control theory and is based on adding dynamic perturbations to the network of interacting biochemical components. In addition to providing a measure of the overall robustness, the method is able to identify specific network fragilities. The importance of considering structural robustness is demonstrated through an analysis of a recently proposed model of the oscillatory metabolism in activated neutrophils. The model displays small parametric sensitivities, but is shown to be highly unrobust to small perturbations in some of the network interactions. Identification of specific fragilities reveals that adding a small delay or diffusion term in one of the involved reactions, likely to exist in vivo, completely removes all oscillatory behaviour in the model.