🧬 Co-Mutation Resistance Pathway Interaction Simulator

The Primary Driver Mutation and First-Line Targeted Therapy

For over a decade, precision oncology has been organized around the idea of the single "actionable" driver mutation: find the one alteration that a tumor is addicted to, match it to a drug, and watch the tumor respond. EGFR L858R — a point mutation in the kinase domain that constitutively activates the receptor — is the textbook case, matched to third-generation EGFR tyrosine kinase inhibitors (TKIs) like osimertinib. But this single-gene framing, while a triumph of molecular medicine, is also an oversimplification that the rest of this simulation exists to correct.

• ~15–20%: EGFR-mutant NSCLC (global) (higher in East Asian, never-smoker cohorts)
• 18.9 mo: FLAURA trial median PFS (osimertinib 1st-line, NEJM 2018)
• ~80%: Objective response rate (osimertinib 1st-line EGFR-mutant NSCLC)
• ~40% / ~45%: L858R vs Exon19del share (of classical EGFR-activating mutations)

Why EGFR L858R is druggable

EGFR (epidermal growth factor receptor) is a receptor tyrosine kinase that, upon ligand binding, dimerizes and autophosphorylates, triggering RAS-MAPK and PI3K-AKT signaling that drives proliferation and survival. L858R is a leucine-to-arginine substitution in the activation loop of the kinase domain (exon 21) that destabilizes the inactive conformation, locking the receptor in an active, ligand-independent state.

This single amino acid change creates a specific vulnerability: the mutant kinase domain has a different ATP-pocket conformation than wild-type EGFR, and small-molecule TKIs can be designed to bind it with high selectivity. Osimertinib, a third-generation, irreversible EGFR TKI, covalently binds Cys797 in the ATP pocket and is selective for both the activating mutation and for the T790M resistance mutation that emerged against earlier-generation drugs — which is why it largely replaced first-generation TKIs (gefitinib, erlotinib) as first-line therapy after the FLAURA trial.

The FLAURA trial (Soria et al., NEJM 2018) established a median progression-free survival of 18.9 months with first-line osimertinib versus 10.2 months for first-generation TKIs, with an objective response rate near 80%. In the initial weeks of treatment, tumor volume typically falls exponentially — the classic "response curve" that this stage visualizes — reflecting rapid apoptosis of the oncogene-addicted cell population.

Why "single driver" is a simplification, not a lie

The single-driver model is not wrong — it correctly predicts that EGFR-mutant tumors respond to EGFR TKIs, and this prediction has changed the standard of care. But every EGFR-mutant tumor also carries a background of co-occurring alterations picked up on the same next-generation sequencing (NGS) panel that found the driver: TP53 mutations in roughly half of cases, less commonly STK11, KEAP1, RB1, PIK3CA, and cell-cycle gene amplifications.

These co-mutations are not incidental noise. They were present in the tumor genome at diagnosis, before treatment even began, and they measurably shift the prognosis and the resistance trajectory — the subject of the remaining stages. A tumor board that reports only "EGFR L858R positive, start osimertinib" and stops reading the panel is discarding prognostic information that is already sitting in the same sequencing file.

Co-Occurring Mutations — Reading the Whole Panel as an Interaction Network

Modern clinical NGS panels (50 to 500+ genes) do not return a single mutation — they return a list. The task of a molecular tumor board is to interpret that list as an interacting system: which co-occurring alterations modify the primary driver's behavior, which are prognostic on their own, and which change the treatment plan. TP53, STK11, and KEAP1 are the three most clinically consequential co-mutation partners in EGFR- and KRAS-driven lung adenocarcinoma.

• ~45–65%: TP53 co-mutation with EGFR (most common EGFR co-alteration)
• ~15–20%: STK11 co-mutation with KRAS (defines an immune-cold subtype)
• ~15–20%: KEAP1 co-mutation with KRAS (often co-occurs with STK11)
• 50–500+: Typical clinical NGS panel size (genes per assay (e.g. FoundationOne))

Building the mutation interaction network

Treating a tumor genome as a network rather than a list means asking, for each pair of co-occurring alterations, whether they are: (a) independent passengers with no interaction, (b) co-selected because they operate in the same pathway and reinforce each other, or (c) mutually exclusive in most tumors, such that their co-occurrence in one tumor signals an unusual biology.

TP53 loss-of-function mutations are the most frequent co-alteration across nearly all solid tumor driver contexts, including EGFR-mutant lung adenocarcinoma. TP53 encodes the "guardian of the genome" — a transcription factor that triggers cell-cycle arrest, senescence, or apoptosis in response to DNA damage or oncogenic stress. Loss of TP53 does not itself drive proliferation, but it removes the brake that would otherwise force an EGFR-activated cell into apoptosis, and it increases genomic instability, accelerating the acquisition of further resistance mutations.

STK11 (which encodes the kinase LKB1) and KEAP1 frequently co-occur with each other and with KRAS mutations. STK11/LKB1 is a master regulator of cellular energy metabolism via AMPK signaling and also suppresses inflammatory signaling; KEAP1 loss constitutively activates the NRF2 oxidative-stress response, giving the tumor a metabolic and antioxidant advantage. Tumors with both KRAS and STK11/KEAP1 alterations characteristically show a markedly "cold" tumor immune microenvironment — low T-cell infiltration, low PD-L1 expression — which becomes clinically decisive in Stage 4 of this simulation.

From gene list to clinical panel report

A molecular tumor board reviewing an NGS report works through several layers: (1) identify actionable driver alterations with an approved or trial-available matched therapy; (2) flag co-occurring alterations with established prognostic or resistance associations (TP53, STK11, KEAP1, RB1); (3) note variant allele frequency (VAF) for each mutation, which can hint at clonal versus subclonal origin — a co-mutation at similar VAF to the driver is likely present in the founding clone, while a much lower VAF may indicate a smaller resistant subclone; (4) cross-reference with liquid biopsy (ctDNA) trends if serial samples are available.

This is why the simulation renders co-mutations as nodes connected to the driver by interaction edges rather than as an independent list: STK11 and KEAP1, for instance, are drawn close together because they so often co-occur and compound each other's immune-evasion effect, while TP53 sits nearer the driver because it directly modifies the apoptotic response to EGFR blockade.

Resistance Mutation Emergence — T790M and C797S at the Drug-Binding Site

No targeted therapy cures a tumor outright at the genomic level; it imposes an enormous selective pressure on a genetically heterogeneous population of cancer cells. Rare pre-existing or newly acquired subclones carrying a mutation that disrupts drug binding are selected for and eventually dominate the relapsing tumor. For EGFR-mutant lung cancer, this is the single most consistent resistance story in precision oncology: T790M against first/second-generation TKIs, and C797S against third-generation osimertinib.

• ~50–60%: T790M in 1st/2nd-gen TKI resistance (of acquired-resistance biopsies)
• ~7–10%: C797S in osimertinib resistance (of resistance cases (rises post-osimertinib))
• ~18–19 mo: Median time to resistance (osimertinib) (consistent with FLAURA PFS)
• untreatable by any single EGFR TKI: C797S in cis with T790M (requires 1st+3rd gen combination)

T790M — the classic gatekeeper mutation

T790M is a threonine-to-methionine substitution at the "gatekeeper" residue of the EGFR ATP-binding pocket. The bulkier methionine side chain sterically hinders binding of first- and second-generation reversible TKIs (gefitinib, erlotinib, afatinib) while also restoring the kinase's affinity for ATP — a double mechanism of resistance that both blocks the drug and re-activates signaling. T790M is found in roughly 50–60% of tumors biopsied at the time of acquired resistance to first/second-generation TKIs, which is precisely why third-generation osimertinib (designed to remain active against T790M-mutant EGFR) was developed and why it moved to first-line use.

On osimertinib, EGFR mutation status must be reinterpreted: because osimertinib already covers T790M, the resistance mechanisms that emerge on osimertinib are different, and a repeat biopsy or liquid biopsy at progression is essential rather than assuming the same T790M pathway is at play.

C797S — resistance to the resistance drug

C797S mutates the cysteine residue (Cys797) that osimertinib covalently binds to form its irreversible bond with the kinase. Losing this cysteine abolishes covalent engagement, restoring ATP binding and kinase activity even in the presence of drug. C797S is detected in roughly 7–10% of osimertinib-resistant biopsies, though the reported frequency varies by treatment line and detection method (tissue vs. ctDNA), and it becomes proportionally more common the earlier osimertinib is used in the treatment sequence.

The clinical significance of C797S depends critically on its allelic relationship to T790M when both are present: if C797S occurs in cis (same DNA strand/allele) with T790M, no single existing EGFR TKI generation can bind the triple-mutant kinase, and a combination of a first-generation reversible TKI (which does not require Cys797) plus a third-generation TKI is sometimes used to cover both mutant and non-T790M alleles. If C797S occurs in trans (different allele) from T790M, first-generation and third-generation TKIs can be combined to independently suppress each resistant clone.

Beyond on-target mutations, roughly 30–40% of osimertinib resistance is driven by off-target/bypass mechanisms rather than a new EGFR mutation at all — MET amplification (~15%), HER2 amplification, PIK3CA mutation, BRAF fusion, or histologic transformation to small-cell lung cancer (~3–10%) — underscoring why comprehensive re-biopsy, not a single-gene resistance assay, is standard at progression.

Roughly half of tumors progressing on osimertinib show no identifiable resistance mutation in EGFR at all — the escape route runs through a bypass pathway (MET amplification being the most common) or through histologic transformation, not through the drug-binding pocket. This is a central argument for full-panel re-biopsy at every progression event rather than testing only for the "expected" resistance mutation.

Co-Mutation Impact on Progression-Free Survival and Drug Sensitivity

Co-occurring mutations are not just resistance bystanders — they independently and measurably shift outcomes. TP53 co-mutation shortens progression-free survival on EGFR TKIs. STK11 and KEAP1 co-mutation with KRAS blunts response to immune checkpoint inhibitors, converting an otherwise reasonable immunotherapy candidate into a poor one. Reading the co-mutation panel changes the survival curve a patient is actually on, not just the theoretical curve for their single driver mutation.

• ~11 mo: EGFR+TP53 median PFS (vs ~22 mo TP53-wildtype, TKI-treated)
• ~7–11%: KRAS+STK11/KEAP1 ORR to ICI (vs ~30–35% KRAS-only (checkpoint inhibitor))
• often high: KEAP1-mutant tumor mutational burden (yet still immunotherapy-resistant)
• worst subtype: TP53 exon 8 (DNA-binding domain) (vs. other TP53 disruptive mutations)

TP53 co-mutation and EGFR TKI outcomes

Multiple cohort studies (including analyses from MSKCC and Chinese multi-center registries) have found that EGFR-mutant, TP53 co-mutated tumors have shorter progression-free and overall survival on EGFR TKIs than TP53-wildtype tumors — commonly reported medians around 11 months versus roughly 22 months. The proposed mechanism connects directly to Stage 1's biology: TP53 loss removes the apoptotic checkpoint that would otherwise eliminate cells under oncogenic and drug-induced stress, and it increases genomic instability, which accelerates the acquisition of on-target and off-target resistance mutations such as T790M and MET amplification.

Not all TP53 mutations are equal: disruptive mutations affecting the DNA-binding domain (particularly exon 8) are associated with worse outcomes than non-disruptive missense variants elsewhere in the gene, illustrating that co-mutation interpretation benefits from mutation-level, not just gene-level, granularity.

STK11/KEAP1 co-mutation and immunotherapy resistance

In KRAS-mutant lung adenocarcinoma, STK11 (LKB1) and KEAP1 co-mutations define a molecular subtype with a distinctly "cold" tumor immune microenvironment: reduced CD8+ T-cell infiltration, low PD-L1 expression, and elevated expression of immunosuppressive cytokines (IL-6, CXCL7). Multiple large retrospective and prospective analyses (including Skoulidis et al., Cancer Discovery 2018, and subsequent validation cohorts) found that objective response rates to PD-1/PD-L1 checkpoint inhibitors drop from roughly 30–35% in KRAS-mutant/STK11-wildtype tumors to roughly 7–11% in KRAS+STK11 co-mutated tumors — despite these tumors often carrying a high tumor mutational burden, a biomarker that would otherwise predict good immunotherapy response.

This is a clinically decisive example of why co-mutation status must inform first-line treatment selection: for a KRAS-mutant, STK11/KEAP1 co-mutated tumor, chemo-immunotherapy combinations or non-immunotherapy strategies may be preferred over single-agent checkpoint inhibition, even though the KRAS driver alone would not predict this.

Liquid biopsy monitoring across the treatment course

Because co-mutation burden and resistance-mutation emergence both evolve over time, serial circulating tumor DNA (ctDNA) monitoring via liquid biopsy has become a core surveillance tool. ctDNA can detect a rising resistance-mutation variant allele frequency (VAF) weeks to months before radiographic progression is visible on CT, allowing earlier treatment adaptation. It also avoids the risks and delays of repeat tissue biopsy, particularly important for co-mutation panels that need updating at every line of therapy. Limitations remain: ctDNA shedding varies by tumor burden and site (CNS metastases shed poorly into peripheral blood), and a negative liquid biopsy at progression does not rule out a resistance mechanism — tissue re-biopsy is still recommended when liquid biopsy is uninformative.

Sequential and Combination Therapy — Adapting the Treatment Line to the Evolving Panel

Because both the co-mutation background and the resistance mutation landscape evolve under treatment pressure, modern management of EGFR- and KRAS-driven lung cancer is not a single prescription but a sequence: start with the best available first-line agent, monitor for resistance with imaging and liquid biopsy, re-biopsy at progression, and select the next line based on what the tumor has actually become — not what it was at diagnosis.

• Osimertinib: Standard 1st-line EGFR+ NSCLC (3rd-gen TKI, covers T790M pre-emptively)
• Osimertinib + MET-TKI: MET-amplified resistance option (e.g. savolitinib, capmatinib combination)
• Sotorasib/adagrasib ± combo: KRAS G12C + STK11/KEAP1 option (chemo-IO often preferred over IO alone)
• ~3–10%: Small-cell transformation on TKI (requires switch to platinum-etoposide regimen)

Building the sequential-therapy timeline

A representative modern sequence for EGFR-mutant NSCLC: (1) first-line osimertinib, chosen because it is active against both the primary activating mutation and pre-emptively against T790M; (2) serial imaging plus periodic or symptom-triggered liquid biopsy to watch for rising ctDNA VAF; (3) upon radiographic or molecular progression, comprehensive re-biopsy (tissue preferred, liquid biopsy if tissue is not feasible) to identify the actual resistance mechanism rather than assuming it; (4a) if an on-target EGFR mutation such as C797S is found, therapy selection depends on its allelic relationship to any co-occurring T790M — trans-configured mutations may respond to a rational combination of first- and third-generation TKIs, while cis-configured triple mutants generally require investigational agents or a switch to chemotherapy-based regimens; (4b) if a bypass mechanism such as MET amplification is found, adding a MET inhibitor (savolitinib, capmatinib) to continued EGFR blockade is an evidence-supported strategy; (4c) if transformation to small-cell histology is found, treatment switches entirely to a platinum-etoposide small-cell regimen, since EGFR-targeted therapy no longer has a molecular target to act on.

Combination strategy for co-mutated KRAS tumors

For KRAS G12C-mutant tumors, the newer KRAS inhibitors (sotorasib, adagrasib) provide direct driver-targeted therapy, but their durability and the choice of combination partner is shaped by the co-mutation panel. In tumors co-mutated with STK11 and/or KEAP1 — the immune-cold subtype described in Stage 4 — checkpoint inhibitor monotherapy underperforms, so chemo-immunotherapy combinations or KRAS-inhibitor-based combinations (e.g., with a MEK inhibitor or SHP2 inhibitor to blunt adaptive resistance signaling) are often favored over an immunotherapy-first approach. Co-occurring TP53 mutation, by contrast, does not by itself argue against a particular drug class but reinforces the expectation of a shorter response duration and the need for closer monitoring.

The unifying principle across all these branches is the same one that motivates this entire simulation: a mutation panel is read as an interacting system across the full arc of treatment, not as a single actionable line item read once at diagnosis. The tumor board's job at each decision point is to ask what the co-mutation and resistance-mutation landscape looks like right now, and choose the next line of therapy accordingly.

The National Comprehensive Cancer Network (NCCN) and IASLC guidelines now explicitly recommend comprehensive genomic re-profiling (tissue or liquid biopsy) at every progression event in EGFR- and KRAS-driven NSCLC, rather than empiric switching to the "next" drug in a fixed sequence — reflecting the field-wide shift from single-driver thinking to co-mutation-pathway thinking that this simulation traces from Stage 1 to Stage 5.