Intervertebral disc (IVD) degeneration (IDD) involves imbalance between the anabolic and the catabolic processes that regulate the extracellular matrix of its tissues. These processes are complex, and improved integration of knowledge is needed. Accordingly, we present a nucleus pulposus cell (NPC) regulatory network model (RNM) that integrates critical biochemical interactions in IVD regulation and can replicate experimental results. The RNM was built from a curated corpus of 130 specialized journal articles. Proteins were represented as nodes that interact through activation and inhibition edges. Semi-quantitative steady states (SS) of node activations were calculated. Then, a full factorial sensitivity analysis (SA) identified which out of the RNM 15 cytokines, and 4 growth factors affected most the structural proteins and degrading enzymes. The RNM was further evaluated against metabolic events measured in non-healthy human NP explant cultures, after 2 days of 1ng/ml IL-1B catabolic induction. The RNM represented successfully an anabolic basal SS, as expected in normal IVD. IL-1B was able to increase catabolic markers and angiogenic factors and decrease matrix proteins. Such activity was confirmed by the explant culture measurements. The SA identified TGF-β and IL1RA as the two most powerful rescue mediators. Accordingly, TGFβ signaling-based IDD treatments have been proposed and IL-1RA gene therapy diminished the expression of proteases. It resulted challenging to simulate rescue strategies by IL-10, but interestingly, IL-1B could not induce IL-10 expression in the explant cultures. Our RNM was confronted to independent in vitro measurements and stands for a unique model, to integrate soluble protein signaling and explore IDD.
