View source

MSH6 and LDLR explanation case studies

This supplementary experiment asks GPT-6 Astra at max reasoning effort to justify its predictions for one deliberately selected panel: OpenSplice MSH6 exon 7 (E01, ENSE00003677325). All 50 benchmark alleles are retained. Selection was based on mechanistic interpretability after reading the paper and measurements, before inspecting this panel's model explanations. It is not a representative performance sample or a case selected for model failure. On September 17, the user also requested the older LDLR promoter example below, retaining its full 50-allele panel before model-trace inspection.

Complete prompts and responses

The exact final responses, including their original Markdown and all 50 final predictions, are committed with this article. Expand either response to read it here. The downloads also include the exact prompts, baseline answers, and data needed to recompute both correlations. The separate provider-exposed summaries are partial records, not complete or verified accounts of internal reasoning.

The full experiment archive (9.8 MB) preserves all 32 original files, including complete provider payloads, question snapshots, request settings, usage, attempt receipts, and baseline provenance. Its manifest records each file's size and SHA-256 digest; checksums also cover the downloadable archive. It is published separately from official benchmark versions and leaderboard results.

Reference interpretation, fixed before model explanations

The OpenSplice preprint shows the mutational landscape in Fig. 2f and the regulatory map in Fig. 6a, group 6. The Results paragraph beginning “Similar diversity is also observed among highly included exons” interprets MSH6 exon 7 specifically. The authors describe a 90-nt exon with roughly 97% wild-type inclusion and a sparse regulatory architecture: 78.4% neutral and 21.6% enhancer-state sequence across the 185-nt native insert. They locate an upstream enhancer-state segment at insert positions 47–50, a 21-nt acceptor-region segment, and a 15-nt segment across the exon end and donor. The reported states summarize substitution and short-deletion effects; they do not identify bound proteins. The 21-nt deletions below were not used to construct those maps. No direct branchpoint or RNA-binding-protein validation for MSH6's upstream segment was identified in the retrieved main text. Calling that segment a proven branchpoint would exceed the available evidence.

That paragraph records the reference available before model inspection. The expanded literature audit below, added on September 18, identifies earlier functional evidence for the branchpoint motif.

The assay uses a FAS three-exon minigene in HEK293T cells. The supplied sequence is in construct orientation. The 185-nt native insert occupies cassette bases 147–331; the tested exon is 217–306. Add 146 to native-insert coordinates to obtain cassette coordinates. Thus the upstream 47–50 segment maps to 193–196, the acceptor AG to 215–216, and the donor GT to 307–308. Junction recognition is a sequence-grounded interpretation; exact effects and the regulatory-state map require measurements.

Allele Sequence feature in the supplied cassette Measured ΔPSI (percentage points) Reference interpretation
V19, G216A Terminal G of acceptor AG −87.466 Large inclusion loss is consistent with disruption of the canonical acceptor.
V44, G307C First G of donor GT −97.003 Near-total loss is consistent with disruption of the canonical donor.
V41, G306A Final exonic base immediately before donor −26.389 A weaker junction effect than V44; an identical near-null prediction would miss the measured gradation.
V24, deletion 241–261 Exon bases 25–45 +0.297 Substantial interior deletion is tolerated in this reporter.
V28, deletion 261–281 Exon bases 45–65 +0.488 Substantial interior deletion is tolerated in this reporter.
V30, deletion 276–296 Exon bases 60–80 +0.415 Substantial interior deletion is tolerated in this reporter.

Deletion intervals describe removed bases, excluding each normalized VCF anchor. Tiny positive measurements are described as tolerance, not evidence of a specific inclusion-enhancing mechanism. These illustrations were fixed before model inspection; the scored experiment still contains all 50 alleles. Measurements come from the pinned benchmark source panel.

The model can locate the junctions, read alleles and deletion lengths, and use the supplied cell and reporter context. It is not given the gene name, wild-type PSI, measured effects, paper, map labels, or these illustrative candidates. Specific protein binding, branchpoint usage, altered RNA structure, and cryptic splice-product identities would need additional evidence; motif plausibility alone does not establish them. Reporter results also do not establish an effect in endogenous tissues or a clinical consequence.

LDLR secondary reference, fixed before model explanations

Kircher et al. (2019), Fig. 1d and the LDLR Results section, describe effects clustered at SP1 and SREBP1/SREBP2 sites. The established sterol-responsive sequence ATCACCCCAC occurs at reporter positions 161–170 in the pinned LDLR panel. The primary functional study identifies that sequence in its reporter-construct Methods. The five selected substitutions within it all reduce activity:

Allele Reporter edit Measured log2 activity effect
V27 T162A −2.070
V28 A164C −2.798
V29 A164T −1.797
V30 C165A −2.252
V31 A169G −2.669

These measurements are consistent with disrupting a positive regulatory element; they do not demonstrate altered protein binding separately for every allele. Sequence localization and effect direction are assessable from the model's explanation. Specific TF occupancy, affinity changes, and precise effects require additional evidence. The benchmark supplies the complete insert in construct orientation and an episomal reporter in HepG2 cells, measured after 24 hours with no added treatment. It does not supply gene identity or a sterol-depletion condition. Predictions premised on an unprovided treatment will be identified.

The source is the study's LDLR replicate 2, as recorded in the pinned satMutMPRA source. All 50 alleles remain in the scored panel. Both the 2019 study and the underlying mechanism predate Astra's cutoff, so this secondary case supplies mechanistic interpretability without the primary case's later-release rationale.

Availability evidence

Astra's documented cutoff is April 30, 2026. The publisher version record dates the first preprint to May 23, 2026, before the Figshare initial release on May 24. The September 16 verification retained the publisher XML and version record and traced the study's predecessor-release evidence; see the provenance audit. The selected source uses Figshare v5. Earliest-supplement versus v5 byte identity for these exact 50 measurements has not been established, and the supplementary TSV was not independently checked. A post-cutoff first verified study date is evidence about public availability, not proof of absence from training. General splice-site knowledge predates this assay.

Explanation request and comparison

The experiment replaces only the optional explanation sentence in the original prompt. It requests substantive sequence-grounded biological justifications, candidate IDs, proposed mechanisms, exceptions, alternatives, and uncertainty, without a fixed word limit. It permits established biological knowledge but asks the model not to substitute assumed gene identities or remembered experimental results for a justification grounded in the supplied input. It also explicitly permits saying that the input does not support a mechanistic explanation. This instruction guides the stated justification; it cannot establish whether recall influenced the prediction. The sequence, geometry, VCF, prediction target, and required final JSON remain unchanged. A distinct template version and question digest identify this supplementary prompt. The full question snapshot retains reference measurements for scoring; only its prompt goes to the model.

Each panel used openai/gpt-6-astra, max effort, seed 20260829, and the existing 128,000-token output ceiling including reasoning. Sampling controls are omitted. These settings follow the requested effort and existing benchmark controls; they are not a provider-endorsed biology recipe. The ceiling is a resource limit, not a requested explanation length.

The experiment has a separate $15 user-authorized budget, unrelated to issue #72. The OpenAI route is pinned with fallbacks disabled and all request parameters required. Prompt and completion prices are capped at $10 and $50 per million tokens, with no per-request fee. The recorded preflight bounds both requests together even before an anticipated 50% batch discount. There are no automatic retries. Exact request bodies, price assumptions, and submission receipts are retained beside the result files.

The September 17 batch returned malformed upstream results for both requests, with no usable model output and $0 charged. A September 18 direct attempt was interrupted locally during the first request, without a saved response or cost receipt. The user then clarified that both panels should run directly. The restart reserves the interrupted request's full possible cost and checks the remaining budget before each request. Prompts and inference settings remain unchanged; the existing streaming transport is used to keep the connection alive. All attempts and the transport changes are recorded separately from the frozen selection. A missing response is not treated as a free request.

The completed MSH6 direct response cost $2.381404. Its receipt left insufficient room for LDLR's full standard-price ceiling while retaining the interrupted request's $6.45879 reserve, so LDLR was withheld locally. The subsequent LDLR request pins OpenAI's openai/flex direct endpoint and caps prices at $5/$25 per million input/output tokens. The prompt, max reasoning effort, seed, and 128,000-token ceiling are unchanged. Including that new request's ceiling, the conservative total is $12.065004. The route snapshot and all earlier records remain beside the experiment artifacts.

The September 15 publication's strongest completed configuration was Astra High (overall Spearman 0.598976), as recorded in the original freeze. The September 18 publication includes the user's completed Astra Max baseline, now strongest at 0.610478. This later availability does not change the selected panels. The matching baseline runs are gpt-6-astra-max-opensplice-snv-20260916 and gpt-6-astra-max-satmut-mpra-20260916, both served by OpenAI with max effort, seed 20260829, and a 128,000-token ceiling. The supplementary experiment adds provider-routing price restrictions and uses direct requests. No duplicate baseline is requested; provider and inference-setting differences are retained in the comparison artifacts.

The comparison distinguishes sequence localization and allele/orientation handling, predicted direction and magnitude, proposed causal links, unsupported specificity, contradictions, and unresolved claims. Explanations are stated justifications, not verified descriptions of internal reasoning. Agreement with the paper is not mechanistic proof. Existing deterministic Spearman and Pearson scoring applies to all 50 predictions. A single completion per condition cannot establish a benefit from requesting explanations or separate prompt effects from sampling variability. Supplementary results do not enter the leaderboard.

Primary result: MSH6

The explanation completion returned valid predictions for all 50 alleles without truncation. Spearman was 0.803453, compared with 0.799818 in the matching Astra Max baseline; Pearson was 0.713871, compared with 0.740322. These small, mixed changes from one completion per condition do not establish an improvement from requesting an explanation.

The six illustrations fixed before model inspection show why the overall score does not settle the mechanistic interpretation:

Allele Baseline prediction Explanation prediction Measured ΔPSI
V19, G216A −21 −10 −87.466
V44, G307C −72 −72 −97.003
V41, G306A −45 −25 −26.389
V24, deletion 241–261 −24 −16 +0.297
V28, deletion 261–281 +5 +8 +0.488
V30, deletion 276–296 +4 +6 +0.415

The explanation correctly localizes the canonical acceptor and donor, distinguishes the final exonic donor base from the intronic GT, and describes V44's mutant donor as CT rather than GC. Its graded V41 prediction closely matches the measured effect. It also excludes retained VCF anchors when describing deletion intervals.

For V19 it reconstructs the newly created AG216–217 and proposes a shifted acceptor, predicting only −10 points despite the measured −87.466. This is a plausible sequence hypothesis with an unsuccessful prediction, not evidence that the alternative junction forms. The model explicitly identifies uncertainty about whether shifted products count toward PSI. A post-response check of the paper's Methods, “PSI per barcode”, finds that isoforms were classified by exact matching to expected inclusion or skipping sequences, with expected exonic deletions incorporated into the inclusion sequence. Arbitrary boundary shifts therefore cannot simply be assumed to rescue the scored inclusion target. Such products can be excluded from both inclusion and skipping counts: the reported PSI is inclusion divided by inclusion plus skipping, not inclusion divided by all possible products. A low measured PSI does not prove that no shifted product formed. Our prompt gives the geometry and PSI target but omits this exact-match counting rule. This target ambiguity limits attribution of the discrepancy to biological reasoning alone; the frozen experiment is not changed after seeing the result.

For the interior, the model proposes enhancer loss for V24 and silencer loss for V28/V30, tying these to short SR-associated, RBFOX-like, and PTBP-like motifs. It explicitly calls binding and activity unproven. The measured near-neutral effects do not support the predicted magnitudes or establish those mechanisms. The paper's largely tolerant interior remains the stronger empirical account. Neither these measurements nor the model's caveats establish which proteins bind.

The completed response used 47,597 output tokens, of which the provider reports 42,838 as reasoning tokens, plus 1,527 input tokens. The $2.381404 direct charge is separate from the unresolved interrupted attempt. The baseline's per-panel cost is allocated from a batch receipt, so a cost difference cannot be attributed solely to the explanation instruction.

What the full MSH6 panel reveals

This expansion is a post-response audit, added on September 18. The six examples above remain the original illustrations. Here we examine every allele, with coordinates in the submitted cassette and deletion anchors excluded. The exon spans 217–306. The upstream enhancer-state segment from OpenSplice maps to 193–196; the sequence-defined acceptor AG|G is 215–217 and donor TAG|GTAAGA is 304–312. These short recognition sequences are not claimed to be the full boundaries of the paper's 21-nt and 15-nt enhancer-state segments. An enhancer state inferred from mutagenesis also does not identify an SR protein.

The following groups partition all 50 alleles by affected bases. Acceptor and donor overlaps take priority over interior or intronic location. A deletion can extend beyond its assigned region. These are descriptive means for the selected panel, not estimates of how all mutations in each region behave.

Region Alleles Mean measured ΔPSI Mean baseline Mean explanation
Upstream, without acceptor overlap 9 −12.494 −4.889 −5.889
Acceptor `AG G` overlap 12 −86.382 −29.000
Interior, without either junction overlap 11 −0.048 −3.545 −1.636
Donor `TAG GTAAGA` overlap 17 −68.499 −61.000
Downstream, without donor overlap 1 −0.477 −2.000 −2.000

The model captures the broad concentration of damaging variants at the junctions, which is consistent with its high rank correlation. It substantially compresses the acceptor losses, and the explanation completion moves both junction means closer to zero than the baseline. Interior means conceal opposing predictions of enhancer loss and silencer loss; individual measurements are more informative.

Earlier experiments on this exon

Additional literature changes the evidential status of the proposed branchpoint. Canson et al. (2022), Fig. 3D tested MSH6 exon 7 in pSPL3 minigenes in Ishikawa cells. Within TTCAT, T191G, C192G, and A193G correspond to c.3557−26, −25, and −24; they reported 35%, 1%, and 6% skipping, respectively. Combined T191G/A193G gave 46%. They also discussed alternative branchpoints, including an experimentally inferred site at −28. This supports the motif's functional relevance and possible redundancy, without establishing A193 as the sole branchpoint in OpenSplice's different reporter and cell line. None of these three substitutions is in our 50-allele panel.

OpenSplice cites Mercer et al. (2015) and Taggart et al. (2017) as references 34 and 35. Canson connects those mapping datasets to this exon and supplies the targeted functional experiment; Canson itself is not cited in the retrieved OpenSplice reference list. Thus the literature connection is through the earlier branchpoint studies, rather than a prior MSH6 experiment cited in its exon example.

There is also evidence in endogenous patient RNA. Šetrajčič Dragoš et al. (2022), Table 1 and Fig. 1G reported 6.8% exon-7-skipping reads for c.3646+5G>A, with none detected in controls. Its deleted transcript interval, r.3557_3646, is the same 90-nt exon; the donor substitution maps to cassette G311A. That exact substitution is absent from our panel, although several deletions remove G311. This independently supports donor sensitivity outside the GT dinucleotide. The patient read fraction is not a mutant-only ΔPSI, and should not be substituted for an OpenSplice effect or used to calibrate a different allele. Neither study validates the model's specific SRSF2, RBFOX, or PTBP assignments in the exon interior.

These exon-specific experiments predate Astra's cutoff. The newer OpenSplice measurements do not imply that this exon's regulatory mechanisms were previously unknown, or establish whether earlier knowledge influenced this response.

Upstream losses cannot be assigned to one branchpoint

V07 removes 184–204 and V09 removes 186–206. Both remove the TTCAT motif, including A193, and the entire 193–196 enhancer-state segment. Yet their measured losses differ by 54.747 points: −29.249 versus −83.996. Astra predicts −22 and −30, getting the order but missing most of the separation. The reference poly-T run occupies 201–213; these deletions leave nine versus seven of its Ts. That offers a plausible contributor, alongside changed spacing and deletion junctions. The two alleles do not isolate tract length or a single branchpoint, so their difference cannot establish any one of those mechanisms.

New AG sequences do not demonstrate scored rescue

Reconstructing each mutant confirms four nearby AGs discussed by the model. Positions below refer to original cassette bases even when deletion brings them together. Sequence creation is directly checkable; splice usage is a hypothesis.

Allele Proposed AG bases Explanation ΔPSI Measured ΔPSI
V16, A215G 214 and 215 +2 −27.564
V17, deletion G216 215 and 217 −7 −90.182
V18, deletion 216–236 215 and 237 −6 −91.585
V19, G216A 216 and 217 −10 −87.466

The shared failure is treating a plausible nearby site as sufficient rescue of the assay's expected inclusion product. The counting limitation described above prevents deciding from these ΔPSIs whether shifted products actually formed. V21 provides a complementary example: G217T preserves the canonical intronic AG but yields −60.673, versus Astra's −6. Acceptor recognition extends beyond the dinucleotide, and correctly finding the AG does not quantify its strength.

The interior tolerates the model's proposed motif losses

All 11 interior alleles have measured effects between −1.120 and +0.624. V24, V28, and V30 each remove 21 bases, leaving a 69-nt exon, yet measure +0.297, +0.488, and +0.415. Astra instead predicts −16, +8, and +6, invoking an SR-associated enhancer (CCAG, 253–256), an RBFOX-like motif (TGCATG, 265–270), and a PTBP-like pyrimidine motif (TCTCT, 282–286). These net effects do not support the predicted magnitudes. They do not prove absence of binding, or exclude compensating effects of a long deletion. The 21-nt deletions were excluded from fitting the paper's SpliceMap, but still come from the same assay; they are not independent biological replication.

The paper's approximately 97% wild-type PSI also leaves only about 3 points of possible inclusion gain. Predictions of +8 and +6 exceed that available range. Wild-type PSI was withheld from the prompt, making this a limitation in effect-size calibration with missing information. It does not by itself negate the ranking score or establish why the model chose those values.

Mutant donor reconstruction is informative, but incomplete

Some of the strongest reasoning is about the actual sequence left after deletion. Here the bar marks the original exon–intron boundary, retained through the edit.

Allele Reconstructed donor context Explanation ΔPSI Measured ΔPSI
V39, deletion 304–306 AAT|GTAAGA −19 −22.460
V40, deletion G306 TTA|GTAAGA −43 −68.912
V41, G306A TAA|GTAAGA −25 −26.389
V44, G307C TAG|CTAAGA −72 −97.003
V45, deletion 309–311 TAG|GTACAT −39 −44.307
V46, deletion 309–314 TAG|GTTTAA −63 −97.047
V47, A309C TAG|GTCAGA −24 −77.079
V49, deletion 310–330 TAG|GTAACA −35 −29.863

Astra correctly orders the one-base V40 deletion as worse than the three-base V39 deletion, and the six-base V46 deletion as worse than the 21-base V49 deletion. In V49, donor position +5 is now the C at original cassette position 332, in the fixed FAS intron. This is a concrete reason to supply the complete reporter context: the deletion draws on sequence outside the native insert. The sizeable V47 miss, despite retaining GT, shows the limit of assigning a moderate penalty from a recognizable donor alone. The broad spatial pattern is convincing; several detailed rescue and motif-strength predictions remain wrong.

Secondary result: LDLR

The direct Flex completion returned valid predictions for all 50 alleles without truncation. Spearman was 0.498629, compared with 0.472678 in the matching Astra Max baseline; Pearson was 0.450955, compared with 0.607670. The explanation therefore gives a slightly better ordering but weaker correspondence between the numerical effect sizes. One completion per condition does not establish a prompt effect, and the service routes differ.

The model correctly identifies ATCACCCCAC at 161–170 and predicts decreased activity for all five preselected substitutions within it. This agrees with the direction of the measurements and the paper's SREBP-site interpretation. It does not establish altered occupancy for any individual allele.

Allele Baseline prediction Explanation prediction Measured log2 effect
V27, T162A −0.90 −0.95 −2.070
V28, A164C −1.60 −1.65 −2.798
V29, A164T −1.50 −1.90 −1.797
V30, C165A −1.55 −2.10 −2.252
V31, A169G −1.50 −1.60 −2.669

The model acknowledges that motif-core recognition supports effect direction more strongly than fine ordering. The ordering of V28 and V29 is nevertheless reversed: it predicts A164T to be more damaging than A164C, whereas A164C has the larger measured loss. The C/T-rich flanking elements are also correctly localized, with Sp-family binding presented as a hypothesis requiring evidence.

What the established LDLR elements reveal

This expanded audit was added after inspecting the responses, on September 18. It does not change the five preselected illustrations or the 50 scored variants. The reference is now more specific than a search for plausible short motifs: satMutMPRA Fig. 1d and Supplementary Table 12 connect the activity landscape to previously studied promoter elements and alleles. The GPN-Star paper, Fig. 4D labels four blocks FP1, SREBP1, SREBP2, and SP1 and examines dependencies within and between them.

The nomenclature needs care. The satMutMPRA supplement, Table 12, pp. 92–93 defines the numbered labels as SRE repeats (with a duplicated “SREBP1” in the caption). They should not be treated as proof that SREBP-1 and SREBP-2 bind different sites. Classic functional studies identify a distal Sp1 site in repeat 1, the SRE within repeat 2, and a proximal Sp1 site in repeat 3. For example, Peeters et al. found that a repeat-1 deletion abolished Sp1 binding in a footprint assay and reduced reporter activity to about 10% of normal. Astra's Sp-family hypothesis for the upstream C/T-rich sequence is therefore compatible with established biology, despite the figure's SREBP1 label.

The following mapping uses the submitted insert's 1-based, inclusive coordinates. Genomic GRCh38 position is reporter position + 11,089,230; the corresponding NM_000527.4 upstream coordinate is position − 319. The latter is relative to translation, not a transcription-start coordinate. These offsets are checked against every selected source allele; V34, for example, is reporter C177T = chr19:11089407 C>T = c.-142C>T. FP1 uses the experimentally footprinted sequence. The other three intervals are the precise cores discussed by Astra, within the known regulatory modules; they do not define complete binding footprints or classify their flanks as inert.

Literature element Reporter sequence interval Selected variants What Astra recognized
FP1 81–100, GAGCTTCACGGGTTAAAAAG V07–V12 Does not identify the footprint; treats G81A/T86A as without a supported direction and the A/T-rich tract as nearly neutral.
Repeat 1, “SREBP1” figure label 126–136, CTCCTCCTCTT V18–V24 Locates the C/T-rich core and proposes Sp-family cooperation, but assigns much smaller losses than measured.
Repeat 2 / SRE, “SREBP2” figure label 161–170, ATCACCCCAC V27–V31 Correctly recognizes the canonical SRE and shared negative direction.
Repeat 3 / SP1 177–187, CTCCTCCCCCT V34–V38 Locates the core but underweights its importance, especially when proposing compensating ETS recruitment.

FP1 is a missed regulatory element, not merely a missing name. Mehta et al. footprinted this exact 20-base sequence, tested substitutions in reporters, and observed sequence-specific nuclear-protein binding. Their induction experiments involved sterol depletion; we do not transfer that treatment or its effect size to this untreated MPRA, nor infer the identity of the bound protein. The independent benchmark measurements nevertheless show losses for all six selected FP1 variants. V07 (G81A) is the most damaging allele in our entire 50-variant panel at −3.784, but both baseline (−0.03) and explanation (0.00) effectively miss it. V09–V12, the four substitutions in the tract Astra discusses at 93–99, measure −1.692 to −2.083 despite predictions of only −0.04 to −0.10. Declining to call this a canonical TATA box is reasonable; treating that lack of a TATA match as insufficient evidence for a substantial effect misses FP1.

Recognition of the SRE is better than recognition of the whole promoter. Every selected variant in each of these four intervals reduces activity. The means below include every selected allele in each interval, including the C128 deletion in repeat 1; the deletion's retained VCF anchor is excluded when assigning its position. These are descriptive summaries of this deliberately selected panel, not new benchmark scores or evidence of a prompt effect.

Interval Alleles Mean measured Mean baseline Mean explanation
FP1 6 −2.078 −0.032 −0.050
Repeat 1 core 7 −2.408 −0.820 −0.364
SRE core 5 −2.317 −1.410 −1.640
SP1 / repeat 3 core 5 −2.691 −0.570 −0.480
All other sampled positions 27 −0.229 −0.097 −0.071

The explanation singles out the SRE as its strongest loss region, although the sampled FP1 and Sp1-associated variants can be equally or more damaging. It also moves the repeat-1 and repeat-3 means closer to zero than the baseline. That pattern helps explain the weaker Pearson correlation without claiming a causal effect of the explanation instruction. “Other” includes regulatory flanks and does not mean biologically unimportant sequence.

Cooperation changes the interpretation of motif gains

Sanchez et al. experimentally showed that SREBP and Sp1 cooperate, that Sp1 orientation matters, and that other tested TF sites could not simply replace the adjacent Sp1 site. An intact SRE therefore does not imply preserved promoter activity after damage to a cooperating element. This is a more constrained reference model than adding independent benefits for every newly created short motif.

The following post-response illustrations examine the complete group for which the explanation proposes beneficial ETS-site changes. They were not part of the preselected five motif substitutions above, and do not change the scored 50-allele panel:

Allele Proposed favorable ETS change Explanation prediction Measured log2 effect
V18, C126T New reverse-strand core +0.15 −1.709
V34, C177T New reverse-strand core +0.15 −3.090
V36, C182T New core offsets C-rich-site loss −0.10 −3.518
V40, A191G New forward-strand core +0.25 +0.504
V43, deletion of G243 Improved flank of an existing core +0.35 +0.755

The new short cores are present in the reconstructed sequences, but their presence does not establish productive factor binding. V18/V34 have the wrong predicted sign, and the proposed compensation in V36 misses a large measured loss. V40/V43 agree in direction; that agreement does not demonstrate the proposed ETS mechanism. The explanation itself cautions that binding, cooperation, and competition require experimental evidence. These are uncertain mechanistic assignments and unsuccessful predictions, rather than fabricated sequence motifs.

V34 is particularly informative because the literature tests the same allele. Table 12 identifies c.-142C>T in the SP1 region and reports −3.09 for LDLR.2, matching our pinned value. In the older transcription-start numbering this is −49C>T: Mozas et al. measured an 80% reduction in HepG2 reporter activity and severely impaired Sp1 binding in gel-shift assays. Those older measurements use a different construct and are not interchangeable with our numerical target, but they independently support loss of an established activator site. Astra notices the new TTCC, yet its proposed +0.15 net effect fails to account for that documented loss mechanism. There is no evidence here that ETS binds the new sequence or rescues activity.

What GPN-Star adds, and what it cannot establish

GPN-Star's Fig. 4D and Supplementary Fig. 14 provide a useful complementary analysis: the model's nucleotide-dependency map resolves the four named regions and connections between them, including the SREBP2-labelled region and SP1. It changes the question from “does this short motif look plausible?” to “has the model captured the interacting regulatory architecture?” Astra's final answer mentions cooperation but misses FP1 and gives unsupported compensation substantial weight at known activator elements. This is our interpretation of its stated explanation, not an observation of its internal computation.

The GPN-Star map is model-derived, not a new binding or double-mutant assay. Its MPRA track also reuses the Kircher measurements: the authors filter by barcode support and significance, average absolute allele effects by position, fill missing positions with zero, and smooth over five positions. The archived LDLR notebook preserves that construction. The track therefore neither supplies independent experimental replication nor tests signed predictions for our 50 alleles. We use it to motivate checking elements and interactions, not to claim a quantitative GPN-Star-versus-Astra performance comparison or experimental epistasis.

Assay handling is generally careful: the model retains deletion anchors, recognizes that no separate minimal promoter is present, and does not assume sterol depletion. It declines to assert TATA-box, upstream-ORF, or methylation mechanisms from insufficient context. We checked its 17 stated reference-sequence features against the supplied insert, as well as the five ETS-related mutant contexts above. These checks establish sequence consistency, not factor occupancy.

The final answer is grounded in sequence positions, but the provider-exposed summary explicitly refers to the “classic LDL receptor promoter.” The prompt's mechanistic framing therefore does not demonstrate removal of recognition or recall. This partial provider summary is not a verified account of internal reasoning, and cannot establish how much recognition affected the predictions. The study and underlying mechanism predate the model's cutoff.

The response used 41,624 output tokens, including 37,639 reported reasoning tokens, plus 1,076 input tokens. Its reported cost was $1.04732125. Together, the two completed explanation responses cost $3.42872525; the original failed batch cost $0, and the interrupted direct request's charge remains unresolved. Keeping its full $6.45879 reserve gives a conservative total of $9.88751525, below the authorized $15. No additional completions are planned.

Reproduction

From the repository root, with the evaluator installed:

uv run --locked python projects/blog-analysis/scripts/mechanism_case.py prepare-pair \
  --output .vepbench/experiments/issue97-msh6-ldlr-20260917

Preparation is offline and refuses to overwrite the experiment directory. It verifies the selected question against the frozen digest and saves both question snapshots, the exact proposed request, and the selection record before inference. Production records use the existing evaluator and retain complete provider payloads, usage, reasoning, failures, and costs outside Git. The small, exact readable responses and comparison extracts are committed beside this article. Preparing these files does not submit a request.

The explicit run --experiment PATH subcommand submits the two-request batch once, reading only OPENROUTER_API_KEY for authentication. Use the existing vepbench batch status and vepbench batch collect commands with the saved batch-state.json to retrieve it. An existing execution plan prevents repeat submission even if no final result was collected.

The recover-direct --experiment PATH subcommand is reserved for the reviewed unbilled batch failure. It checks the receipt and preserved failures, reproduces the original request bodies, and reserves a single direct attempt before sending anything. An existing recovery plan blocks repetition even if no response was saved; unknown or billed outcomes require a separate budget review.

The reviewed September 18 restart uses resume-interrupted-direct. It preserves the interrupted attempt, records each new request before contacting the provider, and bounds total spending including the interrupted request's unresolved charge. It refuses to repeat an existing restart plan.

After both result files are available, reproduce the deterministic comparison:

uv run --locked python projects/blog-analysis/scripts/mechanism_case.py compare \
  --experiment .vepbench/experiments/issue97-msh6-ldlr-20260917 \
  --panel msh6 \
  --baseline PATH_TO_BASELINE_RESULTS.jsonl \
  --results PATH_TO_EXPLANATION_RESULTS.jsonl \
  --output .vepbench/experiments/issue97-msh6-ldlr-20260917/msh6-comparison.json

Repeat with --panel ldlr and a distinct output path for the secondary case.

The committed evidence is generated from those saved records, without model calls:

uv run --locked python projects/blog-analysis/scripts/mechanism_case.py export-evidence \
  --experiment .vepbench/experiments/issue97-msh6-ldlr-20260917 \
  --baseline .vepbench/experiments/issue97-msh6-ldlr-20260917/baseline-20260918/results.jsonl \
  --output /tmp/vepbench-mechanism-evidence

Use a fresh output directory; the exporter refuses replacement. The committed evidence is sufficient for offline score and element-audit replay via uv run --locked pytest projects/blog-analysis/tests/test_mechanism_case.py. This checks the exact prompts, file hashes, all predictions, both correlations, allele mappings and interval summaries without access to private local records. The annotation intervals are defined in projects/blog-analysis/config/ldlr-elements.json and projects/blog-analysis/config/msh6-elements.json; no literature labels were added to the model prompt after inference. The linked HF archive preserves the full provider payloads and attempt receipts separately from the readable extracts.