Systems Biology Modeling of the Complement System Under Immune Susceptible Pathogens
Nehemiah T. Zewde
Rohaine V. Hsu
Dimitrios Morikis
Giulia Palermo- 1Department of Bioengineering, University of California, Riverside, Riverside, CA, United States
- 2Department of Chemistry, University of California, Riverside, Riverside, CA, United States
The complement system is assembled from a network of proteins that function to bring about the first line of defense of the body against invading pathogens. However, complement deficiencies or invasive pathogens can hijack complement to subsequently increase susceptibility of the body to infections. Moreover, invasive pathogens are increasingly becoming resistant to the currently available therapies. Hence, it is important to gain insights into the highly dynamic interaction between complement and invading microbes in the frontlines of immunity. Here, we developed a mathematical model of the complement system composed of 670 ordinary differential equations with 328 kinetic parameters, which describes all three complement pathways (alternative, classical, and lectin) and includes description of mannose-binding lectin, collectins, ficolins, factor H-related proteins, immunoglobulin M, and pentraxins. Additionally, we incorporate two pathogens: (type 1) complement susceptible pathogen and (type 2) Neisseria meningitidis located in either nasopharynx or bloodstream. In both cases, we generate time profiles of the pathogen surface occupied by complement components and the membrane attack complex (MAC). Our model shows both pathogen types in bloodstream are saturated by complement proteins, whereas MACs occupy <<1.0% of the pathogen surface. Conversely, the MAC production in nasopharynx occupies about 1.5–10% of the total N. meningitidis surface, thus making nasal MAC levels at least about eight orders of magnitude higher. Altogether, we predict complement-imbalance, favoring overactivation, is associated with nasopharynx homeostasis. Conversely, orientating toward complement-balance may cause disruption to the nasopharynx homeostasis. Thus, for sporadic meningococcal disease, our model predicts rising nasal levels of complement regulators as early infection biomarkers.
Introduction
The human respiratory tract is a complex system whose primary function is to exchange oxygen and carbon dioxide to meet the metabolic requirements of an organism. However, this physiological need also creates vulnerability through constant exposure to toxic particles and pathogens [1–3]. To counteract this susceptibility, defense mechanism is deployed in the respiratory tract to clear entrapped particles and target invading pathogens [4]. Part of the innate immunity, the complement system also forms a defensive unit in the respiratory tract to identify and remove infectious agents [4–7]. The complement system is composed of three pathways known as alternative (AP), classical (CP), and lectin (LP) that coordinate an immune response to directly engage invading pathogens. Once activated, the complement system achieves host immunity by generating anaphylatoxins, opsonins, and the membrane attack complex (MAC) [8–19]. Additionally, complement forms a bridge between innate and adaptive immunity. Activation of the complement system assists in lymphocyte stimulation, regulation, and antigen presentation. For instance, complement protein C3d, which tags the pathogen surface, has acted as a molecular adjuvant by lowering the threshold for B-cell activation within the lymphoid follicle. Complement-tagged immune complexes also promote the development of germinal center, which subsequently is essential for maintenance of memory B-cell repertoires.
To counteract the deadly effects of complement activation, microorganisms evolved several strategies to evade the immune system and potentially cause infections [20–24]. These counterstrategies include recruitment of complement system regulators to the surface of pathogen, secretion of proteins that inhibit complement activation, and proteases that cleave complement effector molecules. For instance, human-specific Neisseria meningitidis is an important cause of meningitis and septicemia; it is associated with mortality rates ranging from 10 to 40% [25, 26]. Of the many counterstrategies used by N. meningitidis, one method involves evading the complement system by expressing a surface lipoprotein known as factor H binding protein (fHbp) [27–30]. Presence of this lipoprotein promotes survival of N. meningitidis by recruiting the complement regulators factor H (FH) and factor H-like protein 1 (FHL-1) to its surface. In addition to recruiting AP regulators (FH and FHL-1), N. meningitidis uses additional membrane proteins such as porin A (porA) and Meningococcal surface fibril (Msf) to recruit complement regulators C4b-binding protein (C4BP) and vitronectin (Vn), respectively [20]. By utilizing these strategies, N. meningitidis inhibits all three pathways of the complement system (AP, CP, and LP) and effectively interferes with the bactericidal function of complement.
The importance of the complement system against invading microbes is reflected in studies that show how genetic deficiencies in complement lead to recurrent infections [31]. For instance, deficiencies in the classical pathway are associated with Streptococcus pneumoniae bacteremia, and Haemophilus influenzae meningitis [31–33]. In addition to the classical pathway, deficiencies in the lectin pathway are associated with an increased risk of meningococcal and pneumococcal infections [31, 32, 34]. Lastly, deficiencies in the alternative pathway are also associated with invasive diseases. Individuals with inherited and acquired deficiency in C3 develop recurrent and severe infections, such as meningitis, bacteremia, pneumonia, and otitis media [31, 32, 35, 36].
Deficiencies in the MAC components are also associated with an increased predisposition to invasive meningococcal disease and disseminated gonococcal infection [32]. Individuals with defective MAC components have 7,000- to 10,000-fold higher risk of developing meningococcal disease [37]. Furthermore, 40–50% of MAC-deficient individuals will suffer from recurrent infections at a rate ~100–150 times greater than the normal population. Moreover, using complement inhibitors to interfere with MAC assembly might also increase the susceptibility of acquiring bacterial infections [38–40].
Altogether, a competent complement system that spans all three pathways forms a pivotal junction in the fight against pathogens. Knowing this importance, computational models have been assembled to shed light on the dynamics of complement activation and propagation [41–46]. However, to our knowledge, a comprehensive model that covers all three pathways of the complement system has not been developed. Here, we have implemented a computational predictive model of the complement system that contains all three pathways: alternative, classical, and lectin. Building on the frameworks of our previous models [47, 48], the presented approach contains immunoglobulin M (IgM), mannose-binding lectin (MBL), collectins, and ficolins, factor H-related proteins 1–5 (FHR1–5), and pentraxins (C-reactive protein, CRP; serum-amyloid P component, SAP; and long pentraxin 3, PTX3).
Additionally, since complement activation leads to indirect killing through opsonophagocytosis and direct killing through MAC production, we used our model to determine the complement response in the presence of two types of pathogens: (type 1) complement-susceptible pathogen incapable of immune evasion, (type 2) model of N. meningitidis that recruits complement regulators (C4BP, Vn, FH, and FHL-1), and also an activator (FHR-3). Since bactericidal kinetics of MAC may take 3 h [19], our calculations were performed within a 3-h timespan. A summary of the complement cascade used in generating our model is shown in Figure 1. By comparing both pathogen types, we can use the model to examine how the complement activation varies in response to invasion by different pathogen groups. Moreover, comparing the immune response on different pathogen types may also identify potential strategies to boost bactericidal activities against invasive- or therapy-resistant pathogens. These insights can also be used to gain an understanding of which strategies to evade the complement system that can lead to infections. Lastly, incorporating both pathogen types can be used to elucidate how environmental changes, such as the bloodstream or nasal cavity, may alter the response of complement on the level of opsonization or MAC deposition.
Figure 1. Summary of the complement system biochemical network. Activation and propagation of the alternative pathway are shown in blue, the classical pathway is shown in yellow, and the lectin pathway in conjunction with different segments of humoral immunity is shown in magenta. The interconnection between all pathways is shown in green. The continued propagation terminates in the formation of either fC5b-9 (soluble MAC) or surface-bound MAC. LP is used as an umbrella term that equates to pattern recognition molecules, such as MBL trimers and tetramers, collectin liver 1 (CL-L1), collectin kidney 1 (CL-K1), heteromeric CL-L1/K1 known as CL-LK, ficolin-1 (M-ficolin), ficolin-2 (L-ficolin), and ficolin-3 (H-ficolin), MBL-associated serine proteases (MASPs)-1, MASP-2, MASP-3, MBL-associated proteins (MAp19 and MAp44). Complement proteins that make up the regulatory checkpoints such FH, FHL-1, and Vn are shown in red. Activated C1 (C1*) will use C1s* to cleave C4 into C4a and nascent C4b (nC4b). Similar to nC3b, nC4b can also covalently attach to nearby cells and initiate complement.
Results
To initiate our modeling efforts, we focused on complement targeting of type 1 pathogens. We compared the different modes of complement activations: (i) fluid phase (FP) of alternative and classical pathways, (ii) CP, (iii) LP, and (iv) FHR1–5 in generating complement components, such as C3b, proconvertases, C3/C5 convertases, and MAC pores that occupy the pathogen surface within 180 min (Figure 2A). Under FP, about 99.0% of the total pathogen surface remains unoccupied (Figure 2A blue), whereas, under CP (Figure 2A green), the pathogen surface is fully occupied as shown in Figure 2A and Supplementary Material 1 (green). Moreover, under LP (Figure 2A and Supplementary Material 1, magenta) and FHR1–5 (Figure 2A and Supplementary Material 1, yellow), complement components nearly saturate the pathogen surface. The CP-based surface occupation however took longer (~10 min) compared to LP/FHR-based surface occupation (<1 min) as shown in Figure 2A and Supplementary Material 1. Next, we generated time profiles of MAC by using the same modes of the above-mentioned complement activation. As shown in the insets of Figure 2A and Supplementary Material 1, MAC production under all four modes occupied <<1.0% of the total pathogen surface. However, MAC production was the highest under FP activation (Figure 2A inset) and subsequently followed in decreasing order by FHR1–5 (Supplementary Figure 1 inset, yellow), LP (Supplementary Figure 1 inset, magenta), and CP (Supplementary Figure 1 inset, green). Overall, these results show under CP, LP, and FHR1–5 complement-mediated surface occupation is the primary response against type 1 pathogens. And despite occupying <1.0% of the total pathogen surface, FP produced the highest MAC level in comparison to the different modes of complement activation.
Figure 2. Percent of surface occupation (type 1 pathogens) under the different modes of complement activation of (i) FP in blue (ii) CP in green, (iii) LP in magenta, and (iv) FHR1–5 in yellow, and complement system with pentraxins (CSP). (A) FP occupied <1.0% of the pathogen surface, whereas all other modes of complement activation rapidly saturate the pathogen surface. The inset shows membrane attack complex (MAC) production is the highest under FP, but all four modes occupied <<1.0% of the pathogen surface. (B) Pathogen surface is rapidly saturated by CSP components as shown in blue. The inset shows MAC production is similar to (A) inset where MACs occupied <<1.0% of the pathogen surface.
We then combined FP, CP, LP, FHR1–5, and pentraxins (PTX3, SAP, and CRP) to generate a combined mode containing both the complement system and pentraxins (CSP). We compared the effects of CSP on pathogen surface occupation and MAC production. As shown in Figure 2B, type 1 pathogens (blue) are rapidly occupied by CSP in <1 min. Furthermore, as shown in Figure 2B inset (yellow), CSP produce MAC levels that occupy <<1.0% of the pathogen surface in 180 min.
After generating complement profiles on type 1 pathogens, we continued our modeling efforts of the complement system against N. meningitidis in the bloodstream and nasopharynx. N. meningitidis evades the complement system by recruiting complement regulators C4BP, Vn, FH, and FHL-1 to its surface. C4BP and Vn are recruited by using membrane proteins known as porA and Msf, respectively. Additionally, alternative pathways regulators (FH and FHL-1) are also recruited to the surface using a membrane lipoprotein known as fHbp. However, fHbp also recruits a promotor of complement activation known as FHR-3. Hence survival of N. meningitidis may depend on the dynamics between complement regulators (C4BP, Vn, FH, and FHL-1) and promotor (FHR-3).
To better understand this dynamic, we first picked the complement activation mode that favored MAC production (not favoring opsonophagocytosis) since individuals with MAC deficiency have 7,000- to 10,000-fold higher risk of developing meningococcal disease [37] and complement-based opsonophagocytosis is still achievable in the absence of MAC deficiency. And then we generated three conditions: (i) FP activation and recruitment of FHR-3 on N. meningitidis (ii) FP activation and recruitment of C4BP, Vn, FH, and FHL-1 on N. meningitidis, and (iii) FP activation and recruitment of C4BP, Vn, FH, FHL-1, and FHR-3 on N. meningitidis.
To generate percent surface occupation levels and MAC production on N. meningitidis, we started our studies by using nasal complement levels (Supplementary Table 1). As shown in Figure 3A (blue) for condition (i), we observe a surface occupation of about 40.0% within 15 min. This level of surface occupation is maintained until about 80 min, and there after there is another increase in the occupation level to reach a final value of 52.8% in 180 min. Under condition (ii) as shown in Figure 3A (green), we observe a more rapid occupation of N. meningitidis (<5 min). However, this phase of rapid occupation is followed by a transient increase in the available surface and then a subsequent decrease to reach a final occupation level of 81.1%. Lastly, under condition (iii) as shown in Figure 3A (red), we generate a time profile similar to condition (ii). However, under condition (iii), the final occupation level settles to 83.1% in 180 min. After generating complement occupation levels on N. meningitidis, we continued our calculations for MAC levels (Figure 3B). Condition (i) generates the highest MAC level (Figure 3B, blue) with 6.1% of the pathogen surface being occupied by MAC pores. Condition (ii) generates the second highest MAC level (Figure 3B, green) with 2.3% of the pathogen surface being occupied by MAC pores. And lastly, the lowest level of MAC pores is generated under condition (iii) where MACs occupy (Figure 3B, red) 2.1% of the pathogen surface. Altogether, condition (i) produced the highest MAC level, followed in order by conditions (ii) and (iii).
Figure 3. Nasal complement profiles on N. meningitidis that occupy the pathogen surface under three conditions: (i) FP and FHR-3 recruitment in blue; (ii) FP and C4BP, Vn, FH, and FHL-1 recruitment in green; and (iii) FP and FHR-3, C4BP, Vn, FH, and FHL-1 recruitment in red. (A) Conditions (ii) and (iii) rapidly occupy the surface of N. meningitidis and account for 81.1 and 83.1% of its surface, respectively. Condition (i) occupies a lower percent of the bacterial surface (52.8%). (B) Conversely, condition (i) produces the highest membrane attack complex (MAC) level, occupying 6.1% N. meningitidis surface, whereas, conditions (ii) and (iii) led to lower MAC production that occupied 2.3 and 2.1% of the bacterial surface, respectively.
Our model shows a majority of the recruited proteins are not involved in biochemical cascades that propagate/inhibit complement but as entities in complex with the pathogen surface (e.g., pathogen: FHR-3) as shown in Supplementary Figures 2A–C for conditions (i), (ii), and (iii), respectively. For instance, for condition (ii) shown as a green time profile in Supplementary Figure 2B, out of the total 81.1% surface occupation by complement proteins (180 min), Vn (magenta) in complex with just the pathogen surface accounts for 41.1 of 81.1%, FH and FHL-1 (cyan) account for 21.8%, and C4BP (brown) accounts for 9.8%. Thus, a total of 72.7% out of 81.1% of N. meningitidis surface is occupied by just recruiting complement proteins (C4BP, Vn, FH, and FHL-1) as entities with pathogen surface. Moreover, for condition (iii) shown as a red time profile in Supplementary Figure 2C, out of the total pathogen surface occupation of 83.1%, Vn (magenta) in complex with just the pathogen surface accounts for 36.6% of 83.1%, FH and FHL-1 (cyan) account for 19.4%, FHR-3 (yellow) accounts for 11.0%, and C4BP (brown) accounts for 8.7%. Thus, a total of 75.7% out of 83.1% of N. meningitidis surface is occupied by just recruiting complement proteins in complex with pathogen surface.
We next generated time profiles for surface occupation and MAC production in the bloodstream by using the same three above-mentioned conditions. As shown in Figure 4, the surface of N. meningitidis is rapidly occupied under all three conditions. However, under condition (i) shown in blue, complement components occupy 84.5% of the pathogen surface and do not fully occupy the surface as conditions (ii) and (iii), shown in green and red, respectively. The result for condition (ii) is under the time profile of condition (iii). Later, we generated the MAC production profiles under the same conditions. The MAC production under all the three conditions occupied <<1.0% of the pathogen surface as shown in Figure 4 inset. Condition (i), however, generated the highest MACs followed by conditions (ii) and (iii). MAC levels for conditions (ii) and (iii) are overlapping as shown in Figure 4 inset. Overall, these results show MAC levels on N. meningitidis are similar to type 1 pathogens in the bloodstream where MAC levels in both pathogen types occupy <<1.0% of the total surface. However, the MAC production under condition (i) in nasopharynx is about eight orders of magnitude higher to that generated under condition (i) in the bloodstream. Moreover, under conditions (ii) and (iii) in the nasopharynx, MAC levels are about 10 orders of magnitude higher compared to those generated under conditions (ii) and (iii) in the bloodstream.
Figure 4. Bloodstream time profiles on N. meningitidis under three conditions: (i) FP and FHR-3 recruitment in blue; (ii) FP and C4BP, Vn, FH, and FHL-1 recruitment in green; and (iii) FP and FHR-3, C4BP, Vn, FH, and FHL-1 recruitment in red. All three conditions rapidly occupy the bacterial surface. However, condition (i) occupied a lower amount of 84.5% of the bacterial surface. Time profile of condition (ii) is under the profile of condition (iii). Moreover, membrane attack complex (MAC) profiles in the inset show condition (i) produced the highest amount, whereas under conditions (ii) and (iii) MAC levels are lower and overlapping.
Similar to the above-mentioned time profiles, we also generated profiles for recruited complement proteins on the surface of the pathogen to determine if these proteins are mostly propagating/inhibiting complement cascade or bound to the pathogen surface as free entities. The results are shown in Supplementary Figures 3A–C for conditions (i), (ii), and (iii), respectively. In the bloodstream, our model shows almost all of the occupied surface on the N. meningitidis is due to the complement proteins in complex with the pathogen surface as entities, such as pathogen: FHR-3 or pathogen:Vn. For instance, for conditions (ii) shown as a green time profile in Supplementary Figure 3B, out of the total 99.9736% surface occupation by complement proteins (180 min), Vn (magenta) in complex with just the pathogen surface accounts for 57.32%, FH and FHL-1 (cyan) account for 30.44%, and C4BP (brown) accounts for 12.21%. In total, almost all of 99.9736% of N. meningitidis surface is occupied by just recruiting complement proteins C4BP, Vn, FH, and FHL-1, and subsequently leaving little room for other complement components, such as MACs. Moreover, although a similar response is also seen for condition (iii) shown as a red time profile in Supplementary Figure 3C, FHR-3 (yellow) accounts for only 0.14% of the occupied surface, whereas complement regulators (recruited) account for almost all of the remaining percent of surface occupation on N. meningitidis.
We continued our modeling efforts where we examined if FHR-3 was absent in the nasopharynx can enhance complement activation (FP) on N. meningitidis offset the protective effects of recruiting C4BP, Vn, FH, and FHL-1. For this, we first divide our initial concentration into three modules: AP module [C3, factor B (FB), factor D (FD), and properdin], CP module [C1, C1q, (C1rC1s)2, C2, and C4], and terminal module (C5, C6, C7, C8, and C9). We then increased the concentration of a single protein to 20.0% of its serum level, while keeping the concentration of others constant at 12.0% of their serum levels (Supplementary Table 1). This choice is based on the assumption that nasal complement proteins, like that of measured nasal C3 and immunoglobulin G (IgG) [49], can be as high as 20.0% of their serum levels. All modules were compared to the MAC production on N. meningitidis that recruits C4BP, Vn, FH, and FHL-1 and has no complement enhancement in the nasopharynx as mentioned above in Figure 3B for condition (ii). In Figure 5A, we present our results for the terminal module, whereas AP and CP modules are shown in Supplementary Figures 4A,B, respectively. Compared to MAC levels on N. meningitidis without terminal module enhancement (2.3% surface occupation, red), increasing the concentration of C9 had the largest effect in the MAC production with MAC levels occupying 4.0% of the pathogen surface in 180 min (yellow). Furthermore, increasing the concentration of C5 or C6 also led to a higher MAC level. As shown in Figure 5A, increasing the level of C5 increased the MAC production to occupy 3.6% (green) and increasing C6 lead to a MAC level occupying 3.3% of the N. meningitidis surface in 180 min. Conversely, increasing the concentration of C7 (cyan) or C8 (black) had minor effects on the MAC production. Lastly, increasing AP or CP modules also had minor effects on the MAC production (Supplementary Figures 4A,B). Overall, our results show increasing the terminal components C5, C6, and C9 to 20.0% of their serum levels led to an increase in MAC levels despite the absence of FHR-3 and the presence of the protective effects in recruiting C4BP, Vn, FH, and FHL-1.
Figure 5. Nasal complement profiles on N. meningitidis under complement enhancement, varying recruitment capabilities, and absence of FHR-3. (A) Complement proteins are increased to 20.0% of their serum values. Compared to the membrane attack complex (MAC) deposition under N. meningitidis that recruits C4BP, Vn, FH, and FHL-1 (red, 2.3% occupation), the terminal module enhancement of proteins C5 (green, 3.3% occupation), C6 (blue, 3.6% occupation), and C9 (yellow, 4.0% occupation) led to higher MAC levels. However, enhancing the concentrations of C7 (cyan) and C8 (black) had minor effects. Time profile for N. meningitidis that recruits C4BP, Vn, FH, and FHL-1 (red) was increased to show the source of comparison. (B) Compared to the MAC deposition under N. meningitidis that recruits C4BP, Vn, FH, and FHL-1 (red), removing the ability of N. meningitidis to recruit Vn (yellow) had the largest effect with MAC level increasing by 2-fold. This was followed by an increase in the MAC deposition through the removal of FH and FHL-1 (green) recruitment capabilities. Lastly, removing the ability of N. meningitidis to recruit C4BP (blue) also increased the MAC deposition, but this had the smallest effect in comparison to Vn, FH, and FHL-1.
Since enhancing terminal components (C5, C6, and C9) led to an increase in the MAC deposition on N. meningitidis, we next examined which of the recruited complement regulators (C4BP, Vn, FH, and FHL-1) had a greater effect on the MAC modulation (absent FHR-3). For this, we first carried out calculations by using nasal complement levels without enhancement and removing (making association rate constant zero) the ability of N. meningitidis to recruit C4BP, while still having the ability to recruit Vn, FH, and FHL-1. Subsequently, we performed a similar calculation by removing the ability of N. meningitidis to recruit Vn (making association rate constant zero), while still having the ability to recruit C4BP, FH, and FHL-1. We then completed our calculations by performing the same procedures for the recruitment of FH and FHL-1. All permutations are compared to N. meningitidis that recruits all four regulators and have no enhancement in complement activation as shown in Figure 5B (red, 2.3% surface occupation). Our model shows removing the ability of N. meningitidis to recruit Vn led to the largest effect with MAC level occupying 4.2% of the pathogen surface in 180 min (Figure 5B, yellow). This equates to about a 2-fold increase in the MAC deposition compared to N. meningitidis that recruits all four regulators and has no enhancement (Figure 5B, red). This is followed by the loss in the ability to recruit FH and FHL-1, where MACs occupy 3.0% of the surface in 180 min (Figure 5B, green). Lastly, removing the ability of N. meningitidis to recruit C4BP led to the smallest increase in the MAC occupation (2.5%) in 180 min (Figure 5B, blue). Altogether, our model shows removing the ability to recruit Vn had the largest effect on the MAC production, followed by FH and FHL-1, and then C4BP.
We concluded our studies by exploring how the board range in nasal regulatory levels (0.2–0.0002% of their serum levels) affects MAC production on N. meningitidis (absent FHR-3). This range is based on the assumption that the nasal complement regulatory levels can vary similarly to that of mucosal FH [50]. For this, we reduced the concentrations of C1 inhibitor (C1-INH), C4BP, Vn, clusterin (Cn), FH, and FHL-1 between 0.2 and 0.0002% of their respective serum levels (Supplementary Table 1). As shown in Figure 6, as we reduce the concentrations of complement regulators from 0.2 to 0.0002% of their serum levels, we see an increase in the MAC deposition. The MAC production increases from occupying 1.3% of the surface under 0.2% (red) to an occupation level of 6.4% under nasal regulatory levels of 0.02% of their serum levels (blue). Further reducing the nasal regulatory levels to 0.002% (green) or 0.0002% (yellow) of their serum levels has increased the MAC occupation level to about 10% of the pathogen surface. In contrast, increasing the concentration of complement regulators to 2.0% of their serum levels reduced the MAC production to occupy 0.008% of the pathogen surface as shown in Figure 6 in magenta. From this, our model predicts if complement proteins are maintained at 12.0% of their serum levels, and nasal C1-INH, C4BP, Vn, Cn, FH, and FHL-1 are maintained at 0.02% of their serum levels or lower, MAC pores can go from occupying 6.4% (0.02%) to about 10% (under 0.002% or 0.0002%) of the total N. meningitidis surface. Conversely, increasing the concentration of complement regulators to 0.2% led to a decrease in the MAC occupation of 1.3%. Moreover, further increasing the concentration of regulators to 2.0% led to reduced MAC occupation levels of 0.008%. Altogether, our model shows mounting higher levels of MAC may require a concentration barrier between complement proteins and regulators (C1-INH, C4BP, Vn, Cn, FH, and FHL-1). However, as this concentration barrier is reduced, the MAC deposition is potently regulated on the surface of N. meningitidis.
Figure 6. Nasal complement profiles on N. meningitidis with varying regulatory levels and absence of FHR-3. Concentrations of C1-INH, C4BP, Vn, Cn, FH, and FHL-1 are varied between 2.0 and 0.0002% of their respective serum levels. The highest membrane attack complex (MAC) levels with the surface occupation of about 10.0% are produced when the concentration of the complement regulators is reduced to either 0.002% (green) or 0.0002% (yellow) of their serum levels. This is followed by 0.02% (blue) of serum levels where MACs occupied 6.4% of the bacterial surface. Lastly, the MAC production under 0.2% (red) of serum occupies 1.3% of the bacterial surface, and subsequently increasing the complement regulators to 2.0% of their serum levels led to MAC pores occupying 0.008% of the bacterial surface.
Discussion
In this study, we developed a computational model of the complement system that includes all three pathways and the crosstalk with different subsets of the immune system. Our comprehensive model includes the alternative, classical, and lectin pathways of the complement system and components of the humoral immunity such as IgG, IgM, CRP, SAP, and PTX3. Biomolecular details and kinetic parameters of the complement system are assembled from the available literature data, whereas unknown parameters are estimated/assumed as mentioned in section “Methods.” Subsequently, we compare the different phases of complement activation by generating time profiles of the pathogen surface that is occupied by the complement components and MAC pores (C5b-918).
Our computational model shows that pathogens (type 1 and type 2) are rapidly occupied by complement components in the bloodstream. However, in both pathogen types, our model shows MACs occupied <<1.0% of the total pathogen surface. Thus, our model predicts if a pathogen enters the bloodstream, phagocytosis, following opsonization, may serve as the primary mechanism of eliminating pathogens. To highlight this point, in the absence of MACs, specific antibodies in conjunction with complement activation can elicit bacterial killing of N. meningitidis through opsonophagocytosis [51–53]. Furthermore, opsonophagocytosis may also compensate for the MAC deficiency as impairment in assembling MACs increases the risk for Neisseria and disseminated gonococcal infections, whereas not predisposing individuals for contracting infections caused by other Gram-negative pathogens [54]. Altogether, pathogen elimination through opsonization followed by phagocytosis may form a potent immune response in the absence of MAC-based directed response.
In contrast to the terminal cascade in the bloodstream, MAC levels in the nasopharynx were much higher. Production of MAC under condition (i) occupied 6.1% of the total pathogen surface (N. meningitidis), followed by condition (ii) with 2.3% surface occupation, and then condition (iii) with 2.1% surface occupation, as shown in Figure 3B. This difference in nasal MAC production equates to about eight orders of magnitude higher in the nasopharynx, compared to MAC levels on pathogens in the bloodstream. And this elevated MAC production in nasopharynx stems mainly from the differences between levels of complement activators and regulators. Nasopharyngeal secretions of C3 are estimated to be between 4.2 and 20.2% of serum C3 levels [49], whereas levels of mucosal FH are below 0.2% (in a range between 0.2 and 0.0002%) of plasma FH levels [50]. These concentration differences in nasopharynx will favor overactivation of the complement system as nasal C3 levels are 50- to 233-fold greater than nasal FH levels. Conversely, complement activation in the bloodstream is regulated as the mean plasma C3 (7.1 μM) is about 2-fold higher than plasma FH (3.2 μM). It should be noted, unlike other regulators of the alternative pathway, FH is the main regulator of AP that targets complement proteins such as C3(H2O) and C3b or complement complexes such as C3/C5 convertases C3bBb.
Altogether, N. meningitidis may be susceptible to MAC-mediated lysis due to overactivation of the complement system in the nasopharynx. Our model shows that the MAC deposition on N. meningitidis is dependent on the recruitment of complement regulators (Figure 5B). Removing the ability of N. meningitidis to recruit complement regulators led to enhanced MAC deposition. The role of complement regulators is shown in recent genome-wide association studies where variants on the complement regulator factor H gene are associated with susceptibility to meningococcal disease, caused by N. meningitidis [55, 56]. It should be noted that FHL-1 is an alternative splice variant of FH and that the earlier study [56] did also show genetic variation within FHR-3. Although these studies did not find any association with genes of other complement regulators to the susceptibility of meningococcal disease, our model predicts recruitment of Vn as an important factor used to evade the complement system. And this may explain why the majority of the sequenced meningococcal isolates contain Msf (used to recruit Vn) [57].
Increasing the concentration of complement proteins to 20% of their serum values showed the terminal module (C5, C6, and C9) enhanced MAC production and subsequently counteracted the effects of recruiting complement regulators by porA, Msf, and fHbp (Figure 5A). In contrast, increasing the concentrations AP and CP modules had minor effects on the MAC production as shown in Supplementary Figures 4A,B, respectively. To our surprise, increasing the concentration of properdin in the AP module did not enhance the MAC deposition. In our current model, properdin is taken as a stabilizer of the convertases and not as a recognition molecule. This choice in properdin function comes from recent studies showing properdin binding to zymosan, Escherichia coli, human umbilical vein endothelial cells, and N. meningitidis is solely dependent on the initial C3b deposition [58, 59]. Although past studies have shown properdin as a recognition molecule that initiates AP [60–64], caution needs to be taken when interpreting these results depending on the serum conditions of intact C3 and molecular structure of properdin [58, 65]. However, a recent study showed generating a highly polymerized artificial form of properdin by recombinant techniques was able to directly bind to N. meningitidis and elicit an enhanced complement activation [66]. High-order oligomers of properdin have additional functions that include direct binding to pathogens and serving as a platform for the complement activation. And this function may potentially be manipulated for therapeutic interventions to target invasive pathogens, such as N. meningitidis.
Our findings show that levels of complement regulators C1-INH, C4BP, Vn, Cn, FH, and FHL-1 may have to be reassessed in conjunction with their local host microenvironments to determine if they function as protective or not. Balanced complement activation and regulation are essential to avoid inflammatory disorders; however, an imbalanced complement system (favoring overactivation) may signify complement homeostasis in regions, such as the nasopharynx that experiences constant insults from toxins and pathogens [67–69]. And if genetic and environmental factors, which affect complement regulatory (FH) plasma levels [70], also affect nasal levels, i.e., increasing their levels, pathogenic or commensal strains may use the elevated levels of complement regulatory proteins to modulate the MAC deposition and potentially develop into invasive disease. Altogether, if the complement activation is balanced in the nasopharynx as shown in the bloodstream under normal conditions, this new state may promote increased survival of pathogens in the nasopharynx. And this subsequently makes elevated nasal levels of C1-INH, C4BP, Vn, Cn, FH, and FHL-1 as early biomarkers for an individual potentially developing meningococcal disease. A summary of our model prediction is shown in Figure 7.
Figure 7. Model predictions for the susceptibility of developing meningococcal disease. States associated with a dysfunctional complement system due to an overactivation may have protective effects in the nasal cavity against invasive pathogens.
A major limitation of this study is that the findings are exclusively computational. Experimental studies were not conducted to either support or refute the findings. Future studies will need to be conducted to examine the predictions of the model. Additionally, our model does not account for the flow-mediated transport of complement proteins and complexes. This will not account for multiscale dynamics of complement reactions on the pathogen surface that occurs in a different time and spatial scales. Moreover, the presence of flow-mediated transport will affect both the assembly and deposition of MAC pores on the pathogen surface. Subsequently, altering the lysis rate of pathogens. Thus, by not including the transport of MAC proteins, further modeling efforts should be taken to better compare the dynamics of MAC (and complement activation) in the nasal cavity.
Methods
Computational Model
Our model is based on the manual construction of biochemical reactions of the complement system. Figure 1 summarizes the biochemical reactions that form a network of proteins, fragments, and complexes that either activate or regulate the complement system. This network of interactions is also coupled with other segments of the humoral immunity that includes IgM, CRP, SAP, and PTX3. The formulated biochemical model entails the different phases of complement activation, propagation, and regulation in both fluid and cellular surfaces. Human erythrocytes were used as a host cell model, whereas Escherichia coli was used as a pathogen cell model [47]. We incorporated two types of pathogens: (type 1) complement susceptible pathogen, and (type 2) model of N. meningitidis located in either nasopharynx or bloodstream. Moreover, using our biochemical network, we assembled a system of 670 ordinary differential equations (ODEs) with 328 kinetic parameters. Our system was organized into modules of equations that describe the biochemical cascade for initiation and propagation, host cell, pathogen, regulation, proteins, and fragments as shown in Supplementary Equation 1. The equations describe the reaction rates of complement proteins, fragments, and complexes in conjunctions with immunoglobins and pentraxins. Enzymatic reactions are based on Michaelis–Menten kinetics with substrate competitions taken into account for the complement complexes, such as C3/C5 convertases [71]. Model parameters of concentration and kinetic parameters, when available, are acquired from the published literature data or experimental measurements (Supplementary Table 2). We used the ode15s solver of MATLAB (MathWorks, Natick, MA, USA) to solve our ODEs.
For unknown kinetic rate constants in our model, we assumed those parameters to be the same for proteins that are structurally or functionally homologous with experimentally determined parameter values. For instance, the decay rate of C3bBb by FH is known [72], whereas the decay rate of C3(H2O)Bb by the regulatory function of FH is not known. However, we assumed that these two rates to be the same for the following reasons. First, C3(H2O) is known as a C3b-like molecule that has a homologous function as C3b by forming a proconvertase with FB. Second, proconvertase C3(H2O)B, just like C3bB, is activated by the serine protease factor D to form convertases [C3(H2O)Bb and C3bBb]. Lastly, without the FH involvement, both convertases of C3(H2O)Bb and C3bBb have very similar decay rates of 9.0 × 10−3 s−1 and 7.7 × 10−3 s−1, respectively [73]. Similarly, we assumed complement regulators: complement receptor 1 (CR1), decay-accelerating factor (DAF), and C4BP accelerate the decay of C3/C5 convertases with the same kinetic rate constant as FH. We made this assumption because all three regulators are functionally homologous to FH by either inhibiting or deactivating C3/C5 convertases.
We performed global sensitivity analysis by using multiparametric sensitivity analysis (MPSA) as described in our previous paper [47] to identify complement parameters for MAC production under FP activation and recruitment of C4BP, Vn, FH, FHL-1, and FHR-3 on N. meningitidis. We emphasized in the nasopharynx as MAC production was more pronounced unlike the MAC production in bloodstream that occupied <<1% of the pathogen surface. Briefly, we first selected all the kinetic parameters in our biochemical network of the complement system. We then set broad ranges for each parameter value to cover all feasible variations and subsequently used the Latin hypercube sampling method to generate our random parameter vectors. We simulated our model for each parameter set to calculate the objective function value for 1,585 kinetic variations. The objective function is the sum of squared errors between the observed and perturbed system output values. Subsequently, the given set of parameters were determined to be acceptable or unacceptable by comparing the objective function value to a threshold of 50% divisions of the 1,585 sorted objective functions. Lastly, we calculated the cumulative frequency for each parameter with both acceptable and unacceptable cases. The sensitivity of each parameter was evaluated using Kolmogorov–Smirnov statistic. The parameter ranges used in MPSA can be found in Supplementary Table 3 of Supporting Information and the result is shown in Supplementary Figure 5.
Biochemical Model
Alternative Pathway
The alternative pathway auto-activates in the fluid phase by the spontaneous hydrolysis of the complement protein C3. Once activated, C3 will generate a C3b-like molecule known as C3(H2O) [74]. This activated molecule is targeted by the complement protease FB to form a proconvertase known as C3(H2O)B [75]. Subsequently, this complex is targeted by the serine protease FD that cleaves the bound FB complex into Ba and Bb, whereas Ba is released into the fluid, and Bb remains attached to C3(H2O) to form the initial C3 convertase of the alternative pathway, C3(H2O)Bb [73, 74]. This short-lived enzyme catalyzes the cleavage of C3 into a smaller fragment C3a and a larger molecule called nascent C3b (nC3b) [73, 74]. The first product, C3a plays a major role in mediating inflammation by inducing histamine production by mast and basophil cells, resulting in vasodilatation [13, 76, 77], whereas in other cases, C3a plays an anti-inflammatory role by preventing neutrophil migration and degranulation [13, 78]. Nascent fluid C3b (nfC3b) is highly reactive due to its exposed internal thioester bond and is capable of indiscriminately attaching to cell surfaces via covalent linkage [79–81]. However, the newly formed nfC3b has a short half-life of 60 μs before it loses the ability to covalently attach to nearby cell surfaces and remain solely in the fluid phase to form fluid C3b (fC3b) [82]. Similar to the cascade of reactions initiated by C3(H2O), fC3b or surface-bound nC3b can associate with FB, followed by FD cleavage to form the C3/C5 convertase C3bBb [73, 83]. This reaction scheme is also followed by nC3b if it interacts with IgG already in complex with nC3b to form IgGC3bC3b and subsequently interacts with FB and FD to form the C3 convertase IgGC3bC3bBb [84].
Formation of surface-bound convertases will then recruit the complement protein properdin to form the C3/C5 convertase, such as C3bBbP [85, 86]. The recruitment of properdin stabilizes C3/C5 convertases by extending their half-lives by 10-fold and hence initiating the amplification loop. The presence of stabilized convertases will cleave more C3 (into C3a and nC3b) and C5 (into C5a and C5b) than convertases without properdin. The smaller fragments C3a and C5a are mediators of inflammation, but C5a is more potent [13]. However, the larger fragments, nC3b, will initiate the same cascade of reaction by interacting with FB and FD to form C3bBb, whereas C5b initiates the terminal step by first interacting with C6 to form C5b6. This step in complex formation between C5b and C6 is affected by a 5-min short half-life of C5b before irreversibly losing its binding ability to C6 [87]. However, if C5b6 complex is formed, C7 and C8 will sequentially interact C5b6 to form C5b-8. This complex (C5b-8) forms a small pore that later expands by the insertion of multiple C9s to form the MAC [16, 88, 89].
Classical Pathway
Similar to the spontaneous activation of C3 in the alternative pathway, the classical pathway is also spontaneously activated in the fluid phase through its first component of complement, C1 [9–11]. C1 is composed of two complexes of C1q and the heterotetramer (C1rC1s) containing two copies of C1r and C1s serine proteases. Activation of C1 proceeds by the autocatalytic cleavage of C1r, which subsequently activates/cleaves C1s subunits (C1s*). However, nC4b has a short half-life for its exposed thioester and hence becomes rapidly hydrolyzed in solution to generate the fluid C4b (fC4b) [90]. This molecule cannot covalently attach to cellular surfaces but can associate with complement C2 to form the classical pathway proconvertase C4bC2. Soon after, C1* will cleave C2 to form the C3/C5 convertase of the classical pathway, C4b2a. The same cascade of reactions is followed by surface-bound C4b to form C4b2a. On the one hand, the formation of C4b2a leads to cleavage of C3 into C3a and nC3b to propagate the alternative pathway. On the other hand, nC3b or nC4b molecules can also form covalent linkage with each other to form homo- or heterodimers, such as C3bC3b/C4bC4b and C3bC4b that can interact with FB or C2 to form C3/C5 convertases (C3bC3bBb/C4bC4bC2a and C3bC4bBb/C3bC4bC2a) [91–93].
In addition to fluid phase activation, the classical pathway can directly act on cells through the presence or absence of antigen–antibody complexes [94]. In the presence of antigen–antibody immune complexes, C1q (as part of C1) will bind to fragment crystallizable (Fc) portion region of IgG clusters and IgM [94–98]. However, in the absence of antigen–antibody immune complexes, C1q recognizes different targets, such as the short pentraxin CRP and the long pentraxin PTX3. Once C1 binds to its different targets, it follows the same cascade of reaction as in the fluid phase by the auto-activation of C1r, which subsequently activates C1s (C1*). This will lead to cleavage of C4 and C2 to form the convertase C4b2a. Overall, the classical pathway can recognize various targets through the recognition molecule C1q. But once bound, C1 goes through the same cascade of reactions to generate C1*.
Lectin Pathway
Unlike the alternative and classical pathways, the lectin pathway contains numerous pattern recognition molecules (PRMs) in complex with either three types of serine proteases or with two types of non-enzymatic protein molecules. The PRMs of the lectin pathway include MBL, CL-L1, CL-K1, CL-LK, M-ficolin, L-ficolin, and H-ficolin [99, 100]. These molecules have a common feature where they can associate with serine proteases known as MASP-1, MASP-2, and MASP-3. In addition to associating with MASPs, the pattern recognition molecules can also interact with MAp19 and MAp44 that do not contain a serine protease domain [101, 102]. The serine proteases (MASP-1,−2, and−3) are in their zymogen forms and follow subsequent mechanisms to become activated with various cross-activation steps. For instance, MASP-1 in complex with PRM auto-activates and cleaves zymogen MASP-2 and MASP-3 [103, 104]. Furthermore, activated MASP-1 cleaves zymogen MASP-1, whereas activated MASP-2 cross-activates zymogen MASP-1 [103, 105]. However, once PRM–MASP-1 or PRM–MASP-2 are activated, they both can propagate the lectin pathway by cleaving complement proteins C2 and C4 to generate the C3/C5 convertase C4b2a [106, 107]. Once this convertase is formed, it can generate C3a and nC3b or C5a and C5b. And through this mechanism, the lectin pathway can initiate either the alternative pathway or the terminal step to form MACs.
Factor H-Related Proteins
Factor H is the main regulator of an alternative pathway that deactivates C3b and accelerates the decay of C3/C5 convertases. However, factor H belongs to a family protein that includes factor H-like protein 1 (FHL-1) and five factor H-related proteins (FHR-1, FHR-2, FHR-3, FHR-4, and FHR-5) [108–110]. The genes of FHRs are located downstream of the FH gene and each gene codes for complement control proteins (CCP) that are homologous to CCP domains 6–9 and 18–20 of FH [108–110]. However, FHRs lack the regulatory domains of CCP1–4 of FH/FHL-1 and hence implies FHRs may not inhibit complement. Recent evidence suggests that factor H-related proteins may function as positive regulators of the complement system that compete against FH and further enhance complement activation and propagation. For instance, FHR-1 was shown to enhance complement in several ways, such as binding to C3b and allowing convertase formation, not accelerating the decay of convertases, and not inhibiting convertase decay activity of FH [111]. Furthermore, FHR-4 in complex with CRP was able to activate complement and result in C3 and C4 deposition [112]. Lastly, FHR-5 was shown to generate the C3/C5 convertase when incubated with C3b, FB, FD, and properdin [113].
Regulation
Due to the potency of complement activation and propagation, multiple regulatory checkpoints are present to ensure host cell protection. In the alternative pathway, FH, followed by FHL-1, are the two main regulators present in both fluid and cellular surfaces. Both FH and FHL-1 bind to C3b and act as cofactors for Factor I (FI) cleavage of C3b into inactivated C3b (iC3b) [114–118]. Additionally, FH and FHL-1 also target convertases, such as C3bBb and accelerate the rate of decay into its respective components of C3b and Bb [114–116]. Similarly, the classical and lectin pathways also contain two main regulatory proteins known as C1-INH and C4BP that act in both fluid and surface phases. Activated C1, MASP-1, and MASP-2 in complex with PRMs are all regulated by C1-INH [119, 120], whereas C4b is targeted by C4BP and subsequently deactivated into C4d with the help of the serum protease FI [121]. Furthermore, C4BP targets convertases, such as C4b2a and accelerates their decay into their respective components [121].
As a secondary checkpoint, host cells also contain membrane-bound regulators, such as CR1 that regulates activation of the complement system. CR1 regulates in a similar manner as FH/C4BP by binding to either C3b/C4b and using the protease FI to deactivate the complement proteins into iC3b/C4d [122]. Furthermore, iC3b is further cleaved into C3dg by the action of CR1 in collaboration with FI. CR1 can also induce decay acceleration of C3/C5 convertases into their respective parts [123]. In addition, another membrane-bound regulator, DAF, also targets C3/C5 convertases and rapidly dissociates the components [124, 125], whereas DAF only targets convertases, another cell-bound regulator, membrane cofactor protein (MCP), deactivates C3b or C4b with the aid of FI. Since our system uses erythrocytes as a cellular model [47], MCP is not included because it is not expressed on erythrocytes. Lastly, the terminal cascade is regulated through the action of vitronectin or clusterin that inhibits surface binding of C5b-7, whereas CD59 (protectin) inhibits C9 polymerization to regulate MAC pore formation [116, 126].
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
NZ performed the computations, analyzed data, and wrote the manuscript. RH analyzed data. DM conceived the original research. GP supervised research, performed computations, analyzed data, and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This material is based upon work supported by the National Institute of Health through the Grant R01 EY027440. This work was also possible thanks to funding provided to the Palermo lab by the National Science Foundation under Grant no. CHE-1905374 and by the National Institute Of General Medical Sciences of the National Institutes of Health under Award Number R01GM141329 (to GP). This work also used supercomputing resources awarded by the Extreme Science and Engineering Discovery Environment (XSEDE) under allocation award TG-MCB160059 (to GP). Computational resources were also made available by the National Energy Research Scientific Computing Center (NERSC) through the award M3807 (to GP).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor is currently organizing a Research Topic with one of the authors DM.
Acknowledgments
The content of this manuscript has been published as a part of the Ph.D. thesis of Zewde [127]. In honor and memory of our beloved mentor Professor Dimitrios Morikis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2021.603704/full#supplementary-material
The Supporting Information (SI) is provided in separate SI Appendix PDF files, including: (i) SI main text, (ii) Supplementary Table 1, (iii) Supplementary Table 2, (iv) Supplementary Table 3, and (v) Supplementary Equations 1.
References
1. Bosch AATM, Biesbroek G, Trzcinski K, Sanders EAM, Bogaert D. Viral and bacterial interactions in the upper respiratory tract. PLoS Pathog. (2013) 9:1003057. doi: 10.1371/journal.ppat.1003057
2. Gorguner M, Akgun M. Acute inhalation injury. Eurasian J Med. (2010) 42:28–35. doi: 10.5152/eajm.2010.09
3. Dickson RP, Erb-Downward JR, Martinez FJ, Huffnagle GB. The microbiome and the respiratory tract. Annu Rev Physiol. (2016) 78:481–504. doi: 10.1146/annurev-physiol-021115-105238
4. Iwasaki A, Foxman EF, Molony RD. Early local immune defenses in the respiratory tract. Nat Rev Immunol. (2017) 17:7–20. doi: 10.1038/nri.2016.117
5. Nicod LP. Lung defences: an overview. Eur Respir Rev. (2005) 14:45–50. doi: 10.1183/09059180.05.00009501
6. Kulkarni HS, Liszewski MK, Brody SL, Atkinson JP. The complement system in the airway epithelium: an overlooked host defense mechanism and therapeutic target? J Allergy Clin Immunol. (2018) 141:1582–6.e1. doi: 10.1016/j.jaci.2017.11.046
7. Pandya PH, Wilkes DS. Complement system in lung disease. Am J Respir Cell Mol Biol. (2014) 51:467–73. doi: 10.1165/rcmb.2013-0485TR
8. Ricklin D, Reis ES, Mastellos DC, Gros P, Lambris JD. Complement component C3 - The “Swiss Army Knife” of innate immunity and host defense. Immunol Rev. (2016) 274:33–58. doi: 10.1111/imr.12500
9. Ziccardi RJ. Spontaneous activation of the first component of human complement (C1) by an intramolecular autocatalytic mechanism. J Immunol Baltim Md 1950. (1982) 128:2500–4.
10. Tseng Y, Poon PH, Zavodszky P, Schumaker VN. Spontaneous activation of serum C1 in vitro. Role of C1 inhibitor. J Immunol. (1991) 147:1884–90.
11. Bianchino AC, Poon PH, Schumaker VN. A mechanism for the spontaneous activation of the first component of complement, C1, and its regulation by C1-inhibitor. J Immunol. (1988) 141:3930–6.
12. Schatz-Jakobsen JA, Pedersen DV, Andersen GR. Structural insight into proteolytic activation and regulation of the complement system. Immunol Rev. (2016) 274:59–73. doi: 10.1111/imr.12465
13. Merle NS, Noe R, Halbwachs-Mecarelli L, Fremeaux-Bacchi V, Roumenina LT. Complement system part II: role in immunity. Front Immunol. (2015) 6:257. doi: 10.3389/fimmu.2015.00257
14. Reis ES, Mastellos DC, Hajishengallis G, Lambris JD. New insights into the immune functions of complement. Nat Rev Immunol. (2019) 19:503–4. doi: 10.1038/s41577-019-0168-x
15. Hajishengallis G, Reis ES, Mastellos DC, Ricklin D, Lambris JD. Novel mechanisms and functions of complement. Nat Immunol. (2017) 18:1288–98. doi: 10.1038/ni.3858
16. Sharp TH, Koster AJ, Gros P. Heterogeneous MAC initiator and pore structures in a lipid bilayer by phase-plate cryo-electron tomography. Cell Rep. (2016) 15:1–8. doi: 10.1016/j.celrep.2016.03.002
17. Serna M, Giles JL, Morgan BP, Bubeck D. Structural basis of complement membrane attack complex formation. Nat Commun. (2016) 7:ncomms10587. doi: 10.1038/ncomms10587
18. Menny A, Serna M, Boyd CM, Gardner S, Joseph AP, Morgan BP, et al. CryoEM reveals how the complement membrane attack complex ruptures lipid bilayers. Nat Commun. (2018) 9:5316. doi: 10.1038/s41467-018-07653-5
19. Berends ETM, Mohan S, Miellet WR, Ruyken M, Rooijakkers SHM. Contribution of the complement Membrane Attack Complex to the bactericidal activity of human serum. Mol Immunol. (2015) 65:328–35. doi: 10.1016/j.molimm.2015.01.020
20. Hovingh ES, van den Broek B, Jongerius I. Hijacking complement regulatory proteins for bacterial immune evasion. Front Microbiol. (2016) 7:2004. doi: 10.3389/fmicb.2016.02004
21. Bennett KM, Rooijakkers SHM, Gorham RD. Let's tie the knot: marriage of complement and adaptive immunity in pathogen evasion, for better or worse. Front Microbiol. (2017) 8:89. doi: 10.3389/fmicb.2017.00089
22. Zipfel PF, Hallström T, Riesbeck K. Human complement control and complement evasion by pathogenic microbes–tipping the balance. Mol Immunol. (2013) 56:152–60. doi: 10.1016/j.molimm.2013.05.222
23. Lambris JD, Ricklin D, Geisbrecht BV. Complement evasion by human pathogens. Nat Rev Microbiol. (2008) 6:132–42. doi: 10.1038/nrmicro1824
24. Garcia BL, Zwarthoff SA, Rooijakkers SHM, Geisbrecht BV. Novel evasion mechanisms of the classical complement pathway. J Immunol Baltim Md 1950. (2016) 197:2051–60. doi: 10.4049/jimmunol.1600863
25. Stephens DS, Greenwood B, Brandtzaeg P. Epidemic meningitis, meningococcaemia, and Neisseria meningitidis. Lancet Lond Engl. (2007) 369:2196–210. doi: 10.1016/S0140-6736(07)61016-2
26. Martinón-Torres F. Deciphering the burden of meningococcal disease: conventional and under-recognized elements. J Adolesc Health Off Publ Soc Adolesc Med. (2016) 59:S12–20. doi: 10.1016/j.jadohealth.2016.03.041
27. Schneider MC, Prosser BE, Caesar JJE, Kugelberg E, Li S, Zhang Q, et al. Neisseria meningitidis recruits factor H using protein mimicry of host carbohydrates. Nature. (2009) 458:890–3. doi: 10.1038/nature07769
28. Caesar JJE, Lavender H, Ward PN, Exley RM, Eaton J, Chittock E, et al. Competition between antagonistic complement factors for a single protein on N. meningitidis rules disease susceptibility. eLife. (2014) 3:4008. doi: 10.7554/eLife.04008
29. Schneider MC, Exley RM, Chan H, Feavers I, Kang Y-H, Sim RB, et al. Functional significance of factor H binding to Neisseria meningitidis. J Immunol Baltim Md 1950. (2006) 176:7566–75. doi: 10.4049/jimmunol.176.12.7566
30. Granoff DM, Welsch JA, Ram S. Binding of complement factor H (fH) to Neisseria meningitidis is specific for human fH and inhibits complement activation by rat and rabbit sera. Infect Immun. (2009) 77:764–9. doi: 10.1128/IAI.01191-08
31. Skattum L, van Deuren M, van der Poll T, Truedsson L. Complement deficiency states and associated infections. Mol Immunol. (2011) 48:1643–55. doi: 10.1016/j.molimm.2011.05.001
32. Ram S, Lewis LA, Rice PA. Infections of people with complement deficiencies and patients who have undergone splenectomy. Clin Microbiol Rev. (2010) 23:740–80. doi: 10.1128/CMR.00048-09
33. Jönsson G, Truedsson L, Sturfelt G, Oxelius V-A, Braconier JH, Sjöholm AG. Hereditary C2 deficiency in Sweden: frequent occurrence of invasive infection, atherosclerosis, and rheumatic disease. Medicine. (2005) 84:23–34. doi: 10.1097/01.md.0000152371.22747.1e
34. Munthe-Fog L, Hummelshøj T, Honoré C, Madsen HO, Permin H, Garred P. Immunodeficiency associated with FCN3 mutation and ficolin-3 deficiency. N Engl J Med. (2009) 360:2637–44. doi: 10.1056/NEJMoa0900381
35. Wright V, Hibberd M, Levin M. Genetic polymorphisms in host response to meningococcal infection: the role of susceptibility and severity genes. Vaccine. (2009) 27(Suppl.2):B90–102. doi: 10.1016/j.vaccine.2009.05.002
36. S Reis E, Falcão DA, Isaac L. Clinical aspects and molecular basis of primary deficiencies of complement component C3 and its regulatory proteins factor I and factor H. Scand J Immunol. (2006) 63:155–68. doi: 10.1111/j.1365-3083.2006.01729.x
37. Figueroa JE, Densen P. Infectious diseases associated with complement deficiencies. Clin Microbiol Rev. (1991) 4:359–95. doi: 10.1128/CMR.4.3.359
38. Heesterbeek DAC, Angelier ML, Harrison RA, Rooijakkers SHM. Complement and bacterial infections: from molecular mechanisms to therapeutic applications. J Innate Immun. (2018) 10:455–64. doi: 10.1159/000491439
39. Kawakami T, Nakazawa H, Kurasawa Y, Sakai H, Nishina S, Senoo N, et al. Severe infection of Pseudomonas aeruginosa during eculizumab therapy for paroxysmal nocturnal hemoglobinuria. Intern Med Tokyo Jpn. (2018) 57:127–30. doi: 10.2169/internalmedicine.9151-17
40. Vellanki VS, Bargman JM. Aspergillus Niger peritonitis in a peritoneal dialysis patient treated with eculizumab. Ren Fail. (2014) 36:631–3. doi: 10.3109/0886022X.2014.882712
41. Sagar A, Dai W, Minot M, LeCover R, Varner JD. Reduced order modeling and analysis of the human complement system. PLoS ONE. (2017) 12:187373. doi: 10.1371/journal.pone.0187373
42. Lang SN, Germerodt S, Glock C, Skerka C, Zipfel PF, Schuster S. Molecular crypsis by pathogenic fungi using human factor H. A numerical model. PLoS ONE. (2019) 14:e0212187. doi: 10.1371/journal.pone.0212187
43. Liu B, Zhang J, Tan PY, Hsu D, Blom AM, Leong B, et al. A computational and experimental study of the regulatory mechanisms of the complement system. PLoS Comput Biol. (2011) 7:e1001059. doi: 10.1371/journal.pcbi.1001059
44. Hirayama H, Yoshii K, Ojima H, Kawai N, Gotoh S, Fukuyama Y. Linear systems analysis of activating processes of complement system as a defense mechanism. Biosystems. (1996) 39:173–85. doi: 10.1016/0303-2647(96)01617-6
45. Korotaevskiy AA, Hanin LG, Khanin MA. Non-linear dynamics of the complement system activation. Math Biosci. (2009) 222:127–43. doi: 10.1016/j.mbs.2009.10.003
46. Bakshi S, Cunningham F, Nichols E-M, Biedzka-Sarek M, Neisen J, Petit-Frere S, et al. Mathematical modelling of alternative pathway of complement system. Bull Math Biol. (2020) 82:33. doi: 10.1007/s11538-020-00708-z
47. Zewde N Jr RDG, Dorado A, Morikis D. Quantitative modeling of the alternative pathway of the complement system. PLoS ONE. (2016) 11:e0152337. doi: 10.1371/journal.pone.0152337
48. Zewde N, Morikis D. A computational model for the evaluation of complement system regulation under homeostasis, disease, and drug intervention. PLoS ONE. (2018) 13:e0198644. doi: 10.1371/journal.pone.0198644
49. Prellner K, Stjernquist-Desatnik A. A simple method for collecting nasopharyngeal secretions, using cross-linked dextran. Clin Otolaryngol Allied Sci. (1985) 10:85–7. doi: 10.1111/j.1365-2273.1985.tb01172.x
50. Glennie S, Gritzfeld JF, Pennington SH, Garner-Jones M, Coombes N, Hopkins MJ, et al. Modulation of nasopharyngeal innate defenses by viral coinfection predisposes individuals to experimental pneumococcal carriage. Mucosal Immunol. (2016) 9:56–67. doi: 10.1038/mi.2015.35
51. Granoff DM. Relative importance of complement-mediated bactericidal and opsonic activity for protection against meningococcal disease. Vaccine. (2009) 27:B117–25. doi: 10.1016/j.vaccine.2009.04.066
52. Forthal DN. Functions of antibodies. Microbiol Spectr. (2014) 2:1–17. doi: 10.1128/microbiolspec.AID-0019-2014
53. Fijen CA, Bredius RG, Kuijper EJ, Out TA, De Haas M, De Wit AP, et al. The role of Fcgamma receptor polymorphisms and C3 in the immune defence against Neisseria meningitidis in complement-deficient individuals. Clin Exp Immunol. (2000) 120:338–45. doi: 10.1046/j.1365-2249.2000.01208.x
54. Lewis LA, Ram S. Meningococcal disease and the complement system. Virulence. (2014) 5:98–126. doi: 10.4161/viru.26515
55. Martinón-Torres F, Png E, Khor CC, Davila S, Wright VJ, Sim KS, et al. Natural resistance to meningococcal disease related to CFH loci: meta-analysis of genome-wide association studies. Sci Rep. (2016) 6:35842. doi: 10.1038/srep35842
56. Davila S, Wright VJ, Khor CC, Sim KS, Binder A, Breunis WB, et al. Genome-wide association study identifies variants in the CFH region associated with host susceptibility to meningococcal disease. Nat Genet. (2010) 42:772–6. doi: 10.1038/ng.640
57. Andreae CA, Sessions RB, Virji M, Hill DJ. Bioinformatic analysis of meningococcal Msf and Opc to inform vaccine antigen design. PLoS ONE. (2018) 13:e0193940. doi: 10.1371/journal.pone.0193940
58. Harboe M, Johnson C, Nymo S, Ekholt K, Schjalm C, Lindstad JK, et al. Properdin binding to complement activating surfaces depends on initial C3b deposition. Proc Natl Acad Sci USA. (2017) 114:E534–9. doi: 10.1073/pnas.1612385114
59. Harboe M, Garred P, Lindstad JK, Pharo A, Müller F, Stahl GL, et al. The role of properdin in zymosan- and Escherichia coli-induced complement activation. J Immunol Baltim Md 1950. (2012) 189:2606–13. doi: 10.4049/jimmunol.1200269
60. Pillemer L, Blum L, Lepow IH, Ross OA, Todd EW, Wardlaw AC. The properdin system and immunity. I. Demonstration and isolation of a new serum protein, properdin, and its role in immune phenomena. Science. (1954) 120:279–85. doi: 10.1126/science.120.3112.279
61. Cortes C, Ferreira VP, Pangburn MK. Native properdin binds to Chlamydia pneumoniae and promotes complement activation. Infect Immun. (2011) 79:724–31. doi: 10.1128/IAI.00980-10
62. Spitzer D, Mitchell LM, Atkinson JP, Hourcade DE. Properdin can initiate complement activation by binding specific target surfaces and providing a platform for de novo convertase assembly. J Immunol Baltim Md 1950. (2007) 179:2600–8. doi: 10.4049/jimmunol.179.4.2600
63. Saggu G, Cortes C, Emch HN, Ramirez G, Worth RG, Ferreira VP. Identification of a novel mode of complement activation on stimulated platelets mediated by properdin and C3(H2O). J Immunol Baltim Md 1950. (2013) 190:6457–67. doi: 10.4049/jimmunol.1300610
64. Kemper C, Mitchell LM, Zhang L, Hourcade DE. The complement protein properdin binds apoptotic T cells and promotes complement activation and phagocytosis. Proc Natl Acad Sci USA. (2008) 105:9023–8. doi: 10.1073/pnas.0801015105
65. Chen JY, Cortes C, Ferreira VP. Properdin: a multifaceted molecule involved in inflammation and diseases. Mol Immunol. (2018) 102:58–72. doi: 10.1016/j.molimm.2018.05.018
66. Ali YM, Hayat A, Saeed BM, Haleem KS, Alshamrani S, Kenawy HI, et al. Low-dose recombinant properdin provides substantial protection against Streptococcus pneumoniae and Neisseria meningitidis infection. Proc Natl Acad Sci USA. (2014) 111:5301–6. doi: 10.1073/pnas.1401011111
67. Caswell JL. Failure of respiratory defenses in the pathogenesis of bacterial pneumonia of cattle. Vet Pathol. (2014) 51:393–409. doi: 10.1177/0300985813502821
68. Rawlings BA, Higgins TS, Han JK. Bacterial pathogens in the nasopharynx, nasal cavity, and osteomeatal complex during wellness and viral infection. Am J Rhinol Allergy. (2013) 27:39–42. doi: 10.2500/ajra.2013.27.3835
69. Man WH, de Steenhuijsen Piters WAA, Bogaert D. The microbiota of the respiratory tract: gatekeeper to respiratory health. Nat Rev Microbiol. (2017) 15:259–70. doi: 10.1038/nrmicro.2017.14
70. Esparza-Gordillo J, Soria JM, Buil A, Almasy L, Blangero J, Fontcuberta J, et al. Genetic and environmental factors influencing the human factor H plasma levels. Immunogenetics. (2004) 56:77–82. doi: 10.1007/s00251-004-0660-7
71. Schäuble S, Stavrum AK, Puntervoll P, Schuster S, Heiland I. Effect of substrate competition in kinetic models of metabolic networks. FEBS Lett. (2013) 587:2818–24. doi: 10.1016/j.febslet.2013.06.025
72. Wu J, Wu Y-Q, Ricklin D, Janssen BJC, Lambris JD, Gros P. Structure of complement fragment C3b-factor H and implications for host protection by complement regulators. Nat Immunol. (2009) 10:728–33. doi: 10.1038/ni.1755
73. Pangburn MK, Müller-Eberhard HJ. The C3 convertase of the alternative pathway of human complement. Enzymic properties of the bimolecular proteinase. Biochem J. (1986) 235:723–30. doi: 10.1042/bj2350723
74. Pangburn MK, Schreiber RD, Müller-Eberhard HJ. Formation of the initial C3 convertase of the alternative complement pathway. Acquisition of C3b-like activities by spontaneous hydrolysis of the putative thioester in native C3. J Exp Med. (1981) 154:856–67. doi: 10.1084/jem.154.3.856
75. Laich A, Sim RB. Complement C4bC2 complex formation: an investigation by surface plasmon resonance. Biochim Biophys Acta. (2001) 1544:96–112. doi: 10.1016/S0167-4838(00)00208-9
76. Kretzschmar T, Jeromin A, Gietz C, Bautsch W, Klos A, Köhl J, et al. Chronic myelogenous leukemia-derived basophilic granulocytes express a functional active receptor for the anaphylatoxin C3a. Eur J Immunol. (1993) 23:558–61. doi: 10.1002/eji.1830230239
77. el-Lati SG, Dahinden CA, Church MK. Complement peptides C3a- and C5a-induced mediator release from dissociated human skin mast cells. J Invest Dermatol. (1994) 102:803–6. doi: 10.1111/1523-1747.ep12378589
78. Coulthard LG, Woodruff TM. Is the complement activation product C3a a pro-inflammatory molecule? Re-evaluating the evidence and the myth. J Immunol Baltim Md 1950. (2015) 194:3542–8. doi: 10.4049/jimmunol.1403068
79. Law SK, Lichtenberg NA, Levine RP. Evidence for an ester linkage between the labile binding site of C3b and receptive surfaces. J Immunol Baltim Md 1950. (1979) 123:1388–94.
80. Tack BF, Harrison RA, Janatova J, Thomas ML, Prahl JW. Evidence for presence of an internal thiolester bond in third component of human complement. Proc Natl Acad Sci USA. (1980) 77:5764–8. doi: 10.1073/pnas.77.10.5764
81. Law SK, Levine RP. Interaction between the third complement protein and cell surface macromolecules. Proc Natl Acad Sci USA. (1977) 74:2701–5. doi: 10.1073/pnas.74.7.2701
82. Sim RB, Twose TM, Paterson DS, Sim E. The covalent-binding reaction of complement component C3. Biochem J. (1981) 193:115–27. doi: 10.1042/bj1930115
83. Rawal N, Pangburn MK. C5 convertase of the alternative pathway of complement. Kinetic analysis of the free and surface-bound forms of the enzyme. J Biol Chem. (1998) 273:16828–35. doi: 10.1074/jbc.273.27.16828
84. Lutz HU, Jelezarova E. Complement amplification revisited. Mol Immunol. (2006) 43:2–12. doi: 10.1016/j.molimm.2005.06.020
85. Pedersen DV, Roumenina L, Jensen RK, Gadeberg TA, Marinozzi C, Picard C, et al. Functional and structural insight into properdin control of complement alternative pathway amplification. EMBO J. (2017) 36:1084–99. doi: 10.15252/embj.201696173
86. Hourcade DE. The role of properdin in the assembly of the alternative pathway C3 convertases of complement. J Biol Chem. (2006) 281:2128–32. doi: 10.1074/jbc.M508928200
87. Cooper NR, Müller-Eberhard HJ. The reaction mechanism of human C5 in immune hemolysis. J Exp Med. (1970) 132:775–93. doi: 10.1084/jem.132.4.775
88. Bayly-Jones C, Bubeck D, Dunstone MA. The mystery behind membrane insertion: a review of the complement membrane attack complex. Philos Trans R Soc Lond B Biol Sci. (2017) 372:221. doi: 10.1098/rstb.2016.0221
89. Zalman LS, Müller-Eberhard HJ. Comparison of channels formed by poly C9, C5b-8 and the membrane attack complex of complement. Mol Immunol. (1990) 27:533–7. doi: 10.1016/0161-5890(90)90072-8
90. Sepp A, Dodds AW, Anderson MJ, Campbell RD, Willis AC, Law SK. Covalent binding properties of the human complement protein C4 and hydrolysis rate of the internal thioester upon activation. Protein Sci Publ Protein Soc. (1993) 2:706–16. doi: 10.1002/pro.5560020502
91. Meri S, Pangburn MK. A mechanism of activation of the alternative complement pathway by the classical pathway: protection of C3b from inactivation by covalent attachment to C4b. Eur J Immunol. (1990) 20:2555–61. doi: 10.1002/eji.1830201205
92. Rawal N, Pangburn M. Formation of high-affinity C5 convertases of the alternative pathway of complement. J Immunol Baltim Md 1950. (2001) 166:2635–42. doi: 10.4049/jimmunol.166.4.2635
93. Rawal N, Pangburn MK. Formation of high affinity C5 convertase of the classical pathway of complement. J Biol Chem. (2003) 278:38476–83. doi: 10.1074/jbc.M307017200
94. Kishore U, Ghai R, Greenhough TJ, Shrive AK, Bonifati DM, Gadjeva MG, et al. Structural and functional anatomy of the globular domain of complement protein C1q. Immunol Lett. (2004) 95:113–28. doi: 10.1016/j.imlet.2004.06.015
95. Diebolder CA, Beurskens FJ, de Jong RN, Koning RI, Strumane K, Lindorfer MA, et al. Complement is activated by IgG hexamers assembled at the cell surface. Science. (2014) 343:1260–3. doi: 10.1126/science.1248943
96. Bally I, Inforzato A, Dalonneau F, Stravalaci M, Bottazzi B, Gaboriaud C, et al. Interaction of C1q with pentraxin 3 and IgM revisited: mutational studies with recombinant C1q variants. Front Immunol. (2019) 10:461. doi: 10.3389/fimmu.2019.00461
97. Wang G, de Jong RN, van den Bremer ETJ, Beurskens FJ, Labrijn AF, Ugurlar D, et al. Molecular basis of assembly and activation of complement component C1 in complex with immunoglobulin G1 and antigen. Mol Cell. (2016) 63:135–45. doi: 10.1016/j.molcel.2016.05.016
98. Ugurlar D, Howes SC, de Kreuk B-J, Koning RI, de Jong RN, Beurskens FJ, et al. Structures of C1-IgG1 provide insights into how danger pattern recognition activates complement. Science. (2018) 359:794–7. doi: 10.1126/science.aao4988
99. Garred P, Genster N, Pilely K, Bayarri-Olmos R, Rosbjerg A, Ma YJ, et al. A journey through the lectin pathway of complement-MBL and beyond. Immunol Rev. (2016) 274:74–97. doi: 10.1111/imr.12468
100. Degn SE, Thiel S. Humoral pattern recognition and the complement system. Scand J Immunol. (2013) 78:181–93. doi: 10.1111/sji.12070
101. Degn SE, Thiel S, Nielsen O, Hansen AG, Steffensen R, Jensenius JC. MAp19, the alternative splice product of the MASP2 gene. J Immunol Methods. (2011) 373:89–101. doi: 10.1016/j.jim.2011.08.006
102. Degn SE, Hansen AG, Steffensen R, Jacobsen C, Jensenius JC, Thiel S. MAp44, a human protein associated with pattern recognition molecules of the complement system and regulating the lectin pathway of complement activation. J Immunol Baltim Md 1950. (2009) 183:7371–8. doi: 10.4049/jimmunol.0902388
103. Megyeri M, Harmat V, Major B, Végh Á, Balczer J, Héja D, et al. Quantitative characterization of the activation steps of mannan-binding lectin (MBL)-associated serine proteases (MASPs) points to the central role of MASP-1 in the initiation of the complement lectin pathway. J Biol Chem. (2013) 288:8922–34. doi: 10.1074/jbc.M112.446500
104. Degn SE, Kjaer TR, Kidmose RT, Jensen L, Hansen AG, Tekin M, et al. Complement activation by ligand-driven juxtaposition of discrete pattern recognition complexes. Proc Natl Acad Sci USA. (2014) 111:13445–50. doi: 10.1073/pnas.1406849111
105. Héja D, Kocsis A, Dobó J, Szilágyi K, Szász R, Závodszky P, et al. Revised mechanism of complement lectin-pathway activation revealing the role of serine protease MASP-1 as the exclusive activator of MASP-2. Proc Natl Acad Sci USA. (2012) 109:10498–503. doi: 10.1073/pnas.1202588109
106. Ambrus G, Gál P, Kojima M, Szilágyi K, Balczer J, Antal J, et al. Natural substrates and inhibitors of mannan-binding lectin-associated serine protease-1 and−2: a study on recombinant catalytic fragments. J Immunol. (2003) 170:1374–82. doi: 10.4049/jimmunol.170.3.1374
107. Rossi V, Teillet F, Thielens NM, Bally I, Arlaud GJ. Functional characterization of complement proteases C1s/Mannan-binding lectin-associated serine protease-2 (MASP-2) chimeras reveals the higher C4 recognition efficacy of the MASP-2 complement control protein modules. J Biol Chem. (2005) 280:41811–8. doi: 10.1074/jbc.M503813200
108. Józsi M, Tortajada A, Uzonyi B, Goicoechea de Jorge E, Rodríguez de Córdoba S. Factor H-related proteins determine complement-activating surfaces. Trends Immunol. (2015) 36:374–84. doi: 10.1016/j.it.2015.04.008
109. Skerka C, Chen Q, Fremeaux-Bacchi V, Roumenina LT. Complement factor H related proteins (CFHRs). Mol Immunol. (2013) 56:170–80. doi: 10.1016/j.molimm.2013.06.001
110. Medjeral-Thomas N, Pickering MC. The complement factor H-related proteins. Immunol Rev. (2016) 274:191–201. doi: 10.1111/imr.12477
111. Csincsi ÁI, Szabó Z, Bánlaki Z, Uzonyi B, Cserhalmi M, Kárpáti É, et al. FHR-1 binds to C-reactive protein and enhances rather than inhibits complement activation. J Immunol Baltim Md 1950. (2017) 199:292–303. doi: 10.4049/jimmunol.1600483
112. Hebecker M, Okemefuna AI, Perkins SJ, Mihlan M, Huber-Lang M, Józsi M. Molecular basis of C-reactive protein binding and modulation of complement activation by factor H-related protein 4. Mol Immunol. (2010) 47:1347–55. doi: 10.1016/j.molimm.2009.12.005
113. Csincsi ÁI, Kopp A, Zöldi M, Bánlaki Z, Uzonyi B, Hebecker M, et al. Factor H-related protein 5 interacts with pentraxin 3 and the extracellular matrix and modulates complement activation. J Immunol Baltim Md 1950. (2015) 194:4963–73. doi: 10.4049/jimmunol.1403121
114. Dopler A, Guntau L, Harder MJ, Palmer A, Höchsmann B, Schrezenmeier H, et al. Self vs. non-self discrimination by the soluble complement regulators factor H and FHL-1. J Immunol Baltim Md 1950. (2019) 202:2082–94. doi: 10.4049/jimmunol.1801545
115. Zipfel PF, Skerka C. FHL-1/reconectin: a human complement and immune regulator with cell-adhesive function. Immunol Today. (1999) 20:135–40. doi: 10.1016/S0167-5699(98)01432-7
116. Schmidt CQ, Lambris JD, Ricklin D. Protection of host cells by complement regulators. Immunol Rev. (2016) 274:152–71. doi: 10.1111/imr.12475
117. Narkhede YB, Gautam AK, Hsu RV, Rodriguez W, Zewde NT, Harrison RES, et al. Role of electrostatic hotspots in the selectivity of complement control proteins toward human and bovine complement inhibition. Front Mol Biosci. (2021) 8:618068. doi: 10.3389/fmolb.2021.618068
118. Harrison RES, Zewde NT, Narkhede YB, Hsu RV, Morikis D, Vullev VI, et al. Factor H-inspired design of peptide biomarkers of the complement C3d protein. ACS Med Chem Lett. (2020) 11:1054–9. doi: 10.1021/acsmedchemlett.9b00663
119. Rossi V, Bally I, Ancelet S, Xu Y, Frémeaux-Bacchi V, Vivès RR, et al. Functional characterization of the recombinant human C1 inhibitor serpin domain: insights into heparin binding. J Immunol Baltim Md 1950. (2010) 184:4982–9. doi: 10.4049/jimmunol.0902016
120. Kerr FK, Thomas AR, Wijeyewickrema LC, Whisstock JC, Boyd SE, Kaiserman D, et al. Elucidation of the substrate specificity of the MASP-2 protease of the lectin complement pathway and identification of the enzyme as a major physiological target of the serpin, C1-inhibitor. Mol Immunol. (2008) 45:670–7. doi: 10.1016/j.molimm.2007.07.008
121. Ermert D, Blom AM. C4b-binding protein: the good, the bad and the deadly. Novel functions of an old friend. Immunol Lett. (2016) 169:82–92. doi: 10.1016/j.imlet.2015.11.014
122. Krych-Goldberg M, Atkinson JP. Structure-function relationships of complement receptor type 1. Immunol Rev. (2001) 180:112–22. doi: 10.1034/j.1600-065X.2001.1800110.x
123. Krych-Goldberg M, Hauhart RE, Subramanian VB, Yurcisin BM, Crimmins DL, Hourcade DE, et al. Decay accelerating activity of complement receptor type 1 (CD35). Two active sites are required for dissociating C5 convertases. J Biol Chem. (1999) 274:31160–8. doi: 10.1074/jbc.274.44.31160
124. Harris CL, Abbott RJM, Smith RA, Morgan BP, Lea SM. Molecular dissection of interactions between components of the alternative pathway of complement and decay accelerating factor (CD55). J Biol Chem. (2005) 280:2569–78. doi: 10.1074/jbc.M410179200
125. Claire L, Harris DMP. Decay-accelerating factor must bind both components of the complement alternative pathway C3 convertase to mediate efficient decay. J Immunol Baltim Md 1950. (2007) 178:352–9. doi: 10.4049/jimmunol.178.1.352
126. Morgan BP, Boyd C, Bubeck D. Molecular cell biology of complement membrane attack. Semin Cell Dev Biol. (2017) 72:124–32. doi: 10.1016/j.semcdb.2017.06.009
127. Zewde NT. Computational Modeling and Mechanistic Elucidation of the Complement System Under Homeostasis, Disease, and Therapeutic Interventions. (Electronic Thesis and Dissertation). University of California Riverside, Riverside, CA, United States (2019). Available online at: https://escholarship.org/uc/item/5gw6w7g9 (accessed September 14, 2019).
Keywords: predictive modeling, computational systems biology, protein network, predictive network biology, complement system, Neisseria meningitidis
Citation: Zewde NT, Hsu RV, Morikis D and Palermo G (2021) Systems Biology Modeling of the Complement System Under Immune Susceptible Pathogens. Front. Phys. 9:603704. doi: 10.3389/fphy.2021.603704
Received: 07 September 2020; Accepted: 31 March 2021;
Published: 29 April 2021.
Edited by:
Jorge Bernardino De La Serna, Imperial College London, United KingdomReviewed by:
Eva M. Sevick, University of Texas Health Science Center, United StatesGennady Verkhivker, Chapman University, United States
Copyright © 2021 Zewde, Hsu, Morikis and Palermo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giulia Palermo, giulia.palermo@ucr.edu; Dimitrios Morikis, dmorikis@ucr.edu