Eculizumab and pegcetacoplan both inhibit complement — but they act at different nodes. When a patient cohort is activated predominantly at C3, a C5 inhibitor may leave most of the damage unblocked. Transcriptomics lets you look upstream and answer the stratification question before a trial begins.

Why Pathway Position Matters in Complement Therapeutics

The complement system is not a single target; it is a proteolytic cascade with distinct branch points, each susceptible to a different drug class. Eculizumab and ravulizumab (anti-C5 monoclonals) block the cleavage of C5 into C5a and C5b, preventing membrane attack complex (MAC) formation downstream. Pegcetacoplan (anti-C3 peptide) acts further upstream, blocking amplification at C3 convertase and reducing both MAC formation and the generation of C3a and C3b opsonins.

The therapeutic implication is precise: in a cohort where disease pathology is driven by excessive C3 opsonization — as seen in some peripheral neuropathies — anti-C5 therapy leaves the upstream C3 arm intact and may achieve only partial suppression. Conversely, in a MAC-driven membranolytic disease, C5 inhibition is directly on-mechanism. Choosing the right node requires knowing where, in a given tissue and patient population, the cascade is most active.

"In CIDP, IgG4 autoantibodies against nodal and paranodal proteins activate complement at the axonal surface. The dominant damage mechanism — C3d opsonization versus MAC pore formation — varies by antigen and biopsy location."

Bulk RNA-seq from blood or nerve biopsy offers a practical window into this question. Complement system genes — C1QB, C3, C4A, C5, CFB, CFD, CFH, CFI — are robustly transcribed in immune infiltrates, Schwann cells, and macrophages. Their relative expression profiles, interpreted through pathway-level enrichment, can distinguish cohorts where the classical C3 arm dominates from those where late-pathway (C5/MAC) activation is the prevailing signature.

The Three-Step Chained Workflow

BioMate executes this analysis as a sequential chain: each step's output feeds the next, and the final output is a ranked table of complement sub-pathway activation scores with drug-class implications attached.

1
Bulk RNA-seq Differential Expression (DESeq2) Input: raw count matrix from blood or biopsy RNA-seq (e.g. GSE138426 CIDP nerve biopsy, GSE51867 MG PBMC). DESeq2 size-factor normalization, Wald test, shrinkage estimator apeglm. Output: ranked gene list with log2 fold change and adjusted p-value (Benjamini–Hochberg), exported as a .rnk file for downstream enrichment.
2
clusterProfiler / enrichGO Pathway Enrichment Over-representation analysis against GO:0006958 (complement activation), GO:0001906 (cell killing), and MSigDB C2:CP:REACTOME complement gene sets. GSEA using the ranked gene list from Step 1. Complement-specific gene sets: C3/C4/C5/CFB/CFD/CFH modules curated from UniProt complement pathway annotations.
3
Complement Sub-Pathway Ranking Enrichment scores decomposed into C3-convertase signature (C3, C4A, C4B, C2, CFB, CFD, MASP1, MASP2) versus C5-convertase/MAC signature (C5, C6, C7, C8A, C9, CFHR1). Ratio output: C3-dominant / C5-dominant activation index, with 95% bootstrap CI, stratified per cohort. Drug-class recommendation appended based on dominant index.

Reference Datasets: CIDP and Myasthenia Gravis

Two GEO datasets serve as the validation benchmark for this workflow. GSE138426 provides nerve biopsy transcriptomics from CIDP patients versus healthy controls — a tissue-proximal view of complement activation at the site of injury, where macrophage infiltration and Schwann cell reactivity amplify both classical and alternative pathway signals. GSE51867 provides PBMC gene expression from myasthenia gravis patients, capturing the peripheral blood immune milieu including classical pathway activation driven by AChR-specific IgG1 antibodies.

MG with anti-AChR antibodies is the best-characterized complement-mediated neuromuscular disease: IgG1 autoantibodies fix C1q at the neuromuscular junction, activating the classical pathway and driving MAC-mediated destruction of post-synaptic folds. This pathway architecture — classical, C5-convergent — is why eculizumab achieved phase III efficacy in refractory generalized MG (REGAIN trial, Hehir et al., 2017). CIDP presents a more heterogeneous picture: nodal IgG4 antibodies do not fix complement, but sural nerve biopsies show C3d deposition on Schwann cells and axons, implicating macrophage-mediated alternative pathway amplification upstream of C5.

GO PATHWAY ENRICHMENT — CIDP vs CONTROL GO:0005602 complement activation REACTOME: Complement cascade GO:0072562 blood microparticle REACTOME: Initial triggering GO:0006956 complement activation, classical REACTOME: Terminal pathway GO:0035458 cellular response to interferon 0.05 0.10 0.15 0.20 0.25 0.30 Gene Ratio 38 44 28 22 31 18 56 adj. p-value: < 1×10⁻⁷ < 0.001 < 0.05 Gene count
Figure 1. Simulated GO pathway enrichment dot plot from DESeq2 output on GSE138426 (CIDP nerve biopsy vs. control). Bubble size = gene count in pathway; color = adjusted p-value (BH). Classical complement activation (GO:0006956) and the full cascade (GO:0005602) dominate the top-ranked terms, with terminal pathway (REACTOME) scoring lower — consistent with C3-arm predominance in CIDP.

Heatmap: Complement Gene Expression Across Patient Samples

Before sub-pathway ranking, it is useful to visualize the raw normalized expression of key complement genes across individual samples. The heatmap below represents variance-stabilizing-transformed (VST) counts for eight complement-related genes across ten CIDP samples and five controls, revealing inter-sample heterogeneity that the aggregate enrichment score would otherwise collapse.

COMPLEMENT GENE EXPRESSION — VST COUNTS (CIDP NERVE BIOPSY) CFB C3 C5 C4A CFH CFI C1QB C9 CIDP (n=10) CTL (n=5) C1 C2 C3 C4 C5 C6 C7 C8 C9 C10 N1 N2 N3 N4 N5 Low High VST expression CFH: inverse pattern (regulatory — higher in controls)
Figure 2. Simulated heatmap of eight complement genes across 10 CIDP and 5 control nerve biopsy samples (VST-normalized counts). CFB, C3, C4A, and C1QB are strongly upregulated in CIDP (red), confirming classical and alternative pathway activation upstream of C5. CFH is inversely patterned, consistent with loss of complement regulation in disease tissue. C9 (MAC component) shows only moderate upregulation — supporting C3-arm predominance over terminal pathway.

Extending to GBS: Radar Chart of Cohort-Level Sub-Pathway Scores

Adding Guillain-Barré syndrome (GBS) as a third comparator cohort illustrates why the sub-pathway decomposition is clinically valuable. GBS — in particular the axonal variants AMAN and AMSAN — features anti-ganglioside antibodies that fix complement at nodes of Ranvier and drive MAC-mediated axonal injury. Published nerve biopsy data (Hafer-Macko et al.) demonstrate intense C3d and C5b-9 co-deposition in acute GBS axons. The cohort therefore provides a contrast case to CIDP where C5-convertase and terminal pathway scores should be disproportionately elevated.

COMPLEMENT SUB-PATHWAY ACTIVATION SCORES (NORMALIZED) C3 Convertase Classical MAC / Terminal Alternative Regulatory Loss 0.25 0.50 0.75 1.00 MG (anti-AChR) CIDP (nerve biopsy) GBS (AMAN variant) Scores normalized 0–1 per axis from GSEA NES; bootstrapped 95% CI omitted for clarity.
Figure 3. Simulated radar chart comparing complement sub-pathway activation scores across three neuroimmune disease cohorts. MG (red) and GBS/AMAN (green, dashed) show high MAC/Terminal scores — consistent with C5 inhibitor mechanism on-target. CIDP (blue) shows disproportionate C3-convertase and alternative pathway scores with lower MAC enrichment, suggesting that C3-level blockade (pegcetacoplan) may capture a larger fraction of the driving pathology.

Drug-Class Implications from the Dominant Index

The final output of the BioMate chain is a per-cohort dominant-pathway index with drug-class annotations. The table below summarizes the expected stratification across the three diseases, with therapeutics mapped to their mechanism of action and the enrichment signature that predicts on-target activity.

Cohort Dominant Signature C3-Index C5-Index Drug Class Example Agent
MG (anti-AChR+) Classical / MAC 0.55 0.92 C5 Inhibitor Eculizumab, ravulizumab
CIDP (nerve biopsy) C3-Arm / Alternative 0.82 0.45 C3 Inhibitor Pegcetacoplan
GBS (AMAN) MAC / Terminal 0.65 0.90 C5 Inhibitor Eculizumab + cemdisiran (siRNA)

Cemdisiran: siRNA Silencing as a Next-Generation C5 Strategy

Cemdisiran (Alnylam Pharmaceuticals) is a GalNAc-conjugated siRNA that silences hepatic C5 gene expression, reducing circulating C5 protein and thereby suppressing MAC formation without occupying the C5 protein target. Unlike monoclonal antibodies, which require dosing in proportion to circulating C5 concentration, cemdisiran ablates the supply of the substrate protein itself. This makes it particularly attractive for diseases where C5 production is amplified — including acute-phase inflammatory states — and for combination strategies with eculizumab, where cemdisiran pre-dosing reduces the stoichiometric burden on the antibody.

The transcriptomic stratification workflow described here is directly applicable to next-indication selection for cemdisiran: if a patient cohort presents a C5-dominant enrichment profile (high MAC/Terminal GSEA score, elevated C5 and C6–C9 expression relative to CFB and C3), the mechanistic rationale for C5 silencing is clear. If the profile is C3-dominant, the upstream protein supply node matters less than the convertase activity itself, and the siRNA strategy may be less effective on its own.

What BioMate automates in this chain

The three-step workflow — DESeq2 → clusterProfiler enrichGO → sub-pathway ranking — runs end-to-end on BioMate from a raw count matrix. The platform handles size-factor normalization, BH multiple-testing correction, GO and MSigDB annotation download, GSEA permutation testing, and complement module decomposition. Output is a structured findings report with the dominant-pathway index table, heatmap, and drug-class annotation, ready for inclusion in a target selection memo or IND justification.


References

  1. 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
  2. Ricklin D, Hajishengallis G, Yang K, Lambris JD. Complement: a key system for immune surveillance and homeostasis. Nat Immunol. 2010;11(9):785–797. doi:10.1038/ni.1923
  3. Howard JF Jr, Utsugisawa K, Benatar M, et al. Safety and efficacy of eculizumab in anti-acetylcholine receptor antibody-positive refractory generalised myasthenia gravis (REGAIN): a phase 3, randomised, double-blind, placebo-controlled, multicentre study. Lancet Neurol. 2017;16(12):976–986. doi:10.1016/S1474-4422(17)30369-1
  4. Zimmermann TS, et al. GalNAc-siRNA conjugates: robust and reversible inhibition of complement C5 with improved potency and durability. Mol Ther. 2017;25(1):71–78. doi:10.1016/j.ymthe.2016.10.014. For the cemdisiran clinical program in PNH (phase 1/2 combination with eculizumab), see ClinicalTrials.gov NCT03682705.
  5. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–287. doi:10.1089/omi.2011.0118
  6. Hafer-Macko CE, Sheikh KA, Li CY, et al. Immune attack on the Schwann cell surface in acute inflammatory demyelinating polyneuropathy. Ann Neurol. 1996;39(5):625–635.
  7. Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database hallmark gene set collection. Cell Syst. 2015;1(6):417–425. (MSigDB; C2:CP:REACTOME complement pathway)
  8. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. doi:10.1186/s13059-014-0550-8