Identification of Unintuitive Features of Sumoylation through Mathematical Modeling.

Sumoylation is a multistep, multienzymatic post-translational modification in which a small ubiquitin-like modifier protein (SUMO) is attached to the target. We present the first mathematical model for sumoylation including enzyme mechanism details such as autosumoylation of E2 and multifunctional nature of SENP. Simulations and analysis reveal three nonobvious properties for the long term response, modeled as an open system: (i) the steady state sumoylation level is robust to variation in several enzyme properties; (ii) even when autosumoylation of E2 results in equal or higher activity, the target sumoylation levels are lower; and (iii) there is an optimal SENP concentration at which steady state target sumoylation level is maximum. These results are qualitatively different for a short term response modeled as a closed system, where e.g. sumoylation always decreases with increasing SENP levels. Simulations with multiple targets suggest that the available SUMO is limiting, indicating a possible explanation for the experimentally observed low fractional sumoylation. We predict qualitative differences in system responses at short post-translational and longer transcriptional time scales. We thus use this mechanism-based model to explain system properties and generate testable hypotheses for existence and mechanism of unexpected responses.

Method for identification of sensitive nodes in Boolean models of biological networks.

Biological systems are often represented as Boolean networks and analysed to identify sensitive nodes which on perturbation disproportionately change a predefined output. There exist different kinds of perturbation methods: perturbation of function, perturbation of state and perturbation in update scheme. Nodes may have defects in interpretation of the inputs from other nodes and calculation of the node output. To simulate these defects and systematically assess their effect on the system output, two new function perturbations, referred to as 'not of function' and 'function of not', are introduced. In the former, the inputs are assumed to be correctly interpreted but the output of the update rule is perturbed; and in the latter, each input is perturbed but the correct update rule is applied. These and previously used perturbation methods were applied to two existing Boolean models, namely the human melanogenesis signalling network and the fly segment polarity network. Through mathematical simulations, it was found that these methods successfully identified nodes earlier found to be sensitive using other methods, and were also able to identify sensitive nodes which were previously unreported.