Abstract

Interactions between individual members of the B-cell lymphoma 2 (Bcl-2) family of proteins form a regulatory network governing mitochondrial outer membrane permeabilization (MOMP). Bcl-2 family initiated MOMP causes release of the inter-membrane pro-apoptotic proteins to cytosol and creates a cytosolic environment suitable for the executionary phase of apoptosis. We designed the mathematical model of this regulatory network where the synthesis rates of the Bcl-2 family members served as the independent inputs. Using computational simulations, we have then analyzed the response of the model to up-/downregulation of the Bcl-2 proteins. Under several assumptions, and using estimated reaction parameters, a non-linear stimulus-response emerged, whose characteristics are associated with bistability and switch-like behavior. Interestingly, using the principal component analysis (PCA) we have shown that the given model of the Bcl-2 family interactions classifies the random combinations of inputs into two distinct classes, and responds to these by one of the two qualitatively distinct outputs. As we showed, the emergence of this behavior requires specific organization of the interactions between particular Bcl-2 proteins.