Dynamical systems analysis of mitochondrial BAK activation kinetics predicts resistance to BH3 domains

The molecular mechanism underlying mitochondrial BAK activation during apoptosis remains highly controversial. Two seemingly conflicting models have been proposed. In one, BAK requires so-called activating BH3 only proteins (aBH3) to initiate its conformation change. In the other, displacement from inhibitory pro-survival BCL-2 proteins (PBPs) and monomerization of BAK by PBP selective dissociator BH3-only proteins (dBH3) is sufficient.

To better understand the kinetic implications of these conflicting but highly evidence-based models, we have conducted a deterministic, dynamical systems analysis to explore the kinetics underlying the first step of BAK activation, as a non-linear reaction system. We show that dBH3 induced BAK activation is efficient, even in the absence of aBH3s, provided constitutive interaction of PBPs with open conformation BAK occurs in an adenoviral E1B 19K-like manner. The pattern of PBP expression robustly predicts the efficacy of dBH3s.

Our findings accommodate the prevailing BAK activation models as potentially coexisting mechanisms capable of initiating BAK activation, and supports a model based approach for predicting resistance to therapeutically relevant small molecule BH3 mimetics.