A computational model for the evaluation of complement system regulation under homeostasis, disease, and drug intervention

The complement system is an intricate defense network that rapidly removes invading pathogens. Although many complement regulators are present to protect host cells under homeostasis, the impairment of Factor H (FH) regulatory mechanism has been associated with several autoimmune and inflammatory diseases. To understand the dynamics involved in the pivotal balance between activation and regulation, we have developed a comprehensive computational model of the alternative and classical pathways of the complement system. The model is composed of 290 ordinary differential equations with 142 kinetic parameters that describe the state of complement system under homeostasis and disorder through FH impairment. We have evaluated the state of the system by generating concentration-time profiles for the biomarkers C3, C3a-desArg, C5, C5a-desArg, Factor B (FB), Ba, Bb, and fC5b-9 that are influenced by complement dysregulation. We show that FH-mediated disorder induces substantial levels of complement activation compared to homeostasis, by generating reduced levels of C3 and FB, and to a lesser extent C5, and elevated levels of C3a-desArg, Ba, Bb, C5a-desArg, and fC5b-9. These trends are consistent with clinically observed biomarkers associated with complement-mediated diseases. Furthermore, we introduced therapy states by modeling known inhibitors of the complement system, a compstatin variant (C3 inhibitor) and eculizumab (C5 inhibitor). Compstatin demonstrates strong restorative effects for early-stage biomarkers, such as C3a-desArg, FB, Ba, and Bb, and milder restorative effects for late-stage biomarkers, such as C5a-desArg and fC5b-9, whereas eculizumab has strong restorative effects on late-stage biomarkers, and negligible effects on early-stage biomarkers. These results highlight the need for patient-tailored therapies that target early complement activation at the C3 level, or late-stage propagation of the terminal cascade at the C5 level, depending on the specific FH-mediated disease and the manifestations of a patient’s genetic profile in complement regulatory function.