Mathematical model of NF-κB regulatory module

The two-feedback-loop regulatory module of nuclear factor kappaB (NF-kappaB) signaling pathway is modeled by means of ordinary differential equations. The constructed model involves two-compartment kinetics of the activators 1kappaB (IKK) and NF-kappaB, the inhibitors A20 and 1kappaBalpha, and their complexes. In resting cells, the unphosphorylated 1kappaBalpha binds to NF-kappaB and sequesters it in an inactive form in the cytoplasm. In response to extracellular signals such as tumor necrosis factor or interleukin-1, IKK is transformed from its neutral form (IKKn) into its active form (IKKa), a form capable of phosphorylating 1kappaBalpha, leading to 1kappaBalpha degradation. Degradation of 1kappaBalpha releases the main activator NF-kappaB, which then enters the nucleus and triggers transcription of the inhibitors and numerous other genes. The newly synthesized 1kappaBalpha leads NF-kappaB out of the nucleus and sequesters it in the cytoplasm, while A20 inhibits IKK converting IKKa into the inactive form (IKKi), a form different from IKKn, no longer capable of phosphorylating 1kappaBalpha. After parameter fitting, the proposed model is able to properly reproduce time behaviour of all variables for which the data are available: NF-kappaB, cytoplasmic 1kappaBalpha, A20 and 1kappaBalpha mRNA transcripts, IKK and IKK catalytic activity in both wild-type and A20-deficient cells. The model allows detailed analysis of kinetics of the involved proteins and their complexes and gives the predictions of the possible responses of whole kinetics to the change in the level of a given activator or inhibitor. (C) 2004 Elsevier Ltd. All rights reserved.