Golden Gate Overhang Fidelity, Scored Against Real Ligation Data: Three Published Sets vs. Two Naive Designs

What "weakest-link fidelity" and a "risky pair" mean mechanically

In a Golden Gate or MoClo assembly, a Type IIS restriction enzyme (BsaI, BsmBI, etc.) cuts outside its recognition sequence and leaves a 4-nucleotide single-stranded overhang on each fragment end. Two fragments join when their overhangs anneal and get ligated. The enzyme only determines where the cut falls and what overhang sequence is exposed; it is the ligase step (T4 DNA ligase, either as a separate step or as part of a one-pot digestion-ligation with the Type IIS enzyme) that actually decides how specifically that overhang binds and joins to its intended partner versus to some other overhang floating in the same reaction.

Potapov et al. (2018) and Pryor et al. (2020) measured this directly and exhaustively: raw ligation counts for every pairwise combination among all 256 possible four-base overhangs, under real T4 ligase conditions and under real one-pot BsaI-HFv2 conditions respectively. In that data, the count for overhang X ligating to its own perfect complement is the "correct" signal; the counts for X ligating to every other overhang actually present in a given reaction are cross-talk — real, measured mis-ligation events, not a computational estimate of similarity.

SeqBench's scorer (scoreOverhangSet() in the ligase-fidelity module) takes a candidate set of overhangs you actually intend to use together and, for each overhang X in that set, computes fidelity(X) = correctSignal(X) / (correctSignal(X) + crossTalkSignal(X)), where crossTalkSignal(X) sums only the ligation counts between X and the other overhangs that are actually in your set — not all 256 possible overhangs, since only the ones sharing your tube can cross-react with each other. A set's weakest-link fidelity is the minimum of those per-overhang fidelities, not the average; a design is only as reliable as its worst junction, and an average can look fine while one overhang is quietly bad. Mean fidelity is reported alongside for the same reason: to catch cases where the mean and the minimum tell different stories.

A risky pair is an ordered pair (A, B) within your candidate set where B's cross-reaction count against A is at least 5% of A's own correct-signal-plus-that-cross-reaction total (a documented, adjustable default, not a universal biological cutoff). Because the calculation is directional and set-dependent, the same overhang can score cleanly in one candidate set and get flagged in another — cross-talk depends on exactly which other overhangs are competing for the same sticky end in the same reaction, which is why this has to be checked per design, not looked up once as a property of a single overhang.

The honesty point, in full: this does not match NEB's or Potapov's own numbers

SeqBench built this scorer directly from the raw, publicly mirrored 256×256 ligation-count matrices underlying Potapov et al. 2018 and Pryor et al. 2020 — not from NEB's or the papers' own aggregate scoring code, which is not published. An attempt was made to re-derive the exact formula NEB's Ligase Fidelity Tools and the papers themselves use to collapse a whole overhang set down to a single published aggregate fidelity percentage, using these same raw counts. That attempt did not reproduce their published numbers, and the exact aggregation formula is not disclosed anywhere accessible.

The clearest illustration is Pryor et al. 2020's own 20-overhang set (Figure 4, an 11-overhang plant-standard core plus 9 GetSet-added overhangs). The paper itself reports roughly 80% aggregate fidelity for this set using NEB's undisclosed formula. Scoring the identical 20 overhangs against the identical underlying raw counts with SeqBench's transparent per-overhang correct-signal-versus-cross-talk methodology gives 99.89% weakest-link fidelity (generic T4 dataset) and 99.74% (BsaI-HFv2 dataset) — see the next section. Neither number is "wrong"; they answer different questions using different formulas, and only one of those formulas (SeqBench's) is documented well enough for a third party to reproduce from the raw data. Do not compare a weakest-link percentage from this methodology directly against an aggregate percentage from NEB's tool or from a paper's own reported figure as if they were the same statistic — they are not, and no conversion factor between them is published.

What this methodology is useful for is comparing your own candidate overhangs against a documented, reproducible calculation, and against the real published-set benchmarks computed the same way in this article.

Two real ligation datasets, two different reaction chemistries

Every set below is scored against two independently generated real datasets, because "Golden Gate fidelity" is not one fixed property of a ligase — it depends on the actual reaction conditions.

• generic-t4-37c-1h — T4 DNA ligase alone, 1 hour at 37°C, enzyme-agnostic. Source: Potapov V, Ong JL, Kucera RB, Langhorst BW, Bilotti K, Pryor JM, Cantor EJ, Canton B, Knight TF, Evans TC Jr, Lohman GJS. "Comprehensive Profiling of Four Base Overhang Ligation Fidelity by T4 DNA Ligase and Application to DNA Assembly." ACS Synth Biol. 2018;7(11):2665-2674. doi:10.1021/acssynbio.8b00333. Because the Type IIS enzyme only determines the cut site, not the ligation chemistry, this dataset applies to any Type IIS Golden Gate enzyme paired with a standard T4 ligation step.
• bsai-hfv2 — real one-pot BsaI-HFv2 digestion-ligation, where enzyme and ligase act together in the same reaction rather than as separate steps. Source: Pryor JM, Potapov V, Kucera RB, Bilotti K, Cantor EJ, Lohman GJS. "Enabling one-pot Golden Gate assemblies of unprecedented complexity using data-optimized assembly design." PLOS ONE. 2020;15(9):e0238592. doi:10.1371/journal.pone.0238592. Use this dataset specifically when your actual protocol uses BsaI-HFv2 in a one-pot format; it is not a generic proxy for other enzymes.

Three published overhang sets, scored

These are real, previously published Golden Gate/MoClo overhang sets, scored exactly as specified in their source figures/tables (or, for the CIDAR set, cross-validated against secondary sources where noted).

• 11 standard plant-synbio overhangs — GGAG, TGAC, TCCC, TACT, CCAT, AATG, AGCC, TTCG, GCTT, GGTA, CGCT. Source: Pryor et al. 2020, PLOS ONE 15(9):e0238592, Figure 4. generic-t4-37c-1h: weakest-link fidelity 99.95%, mean fidelity 99.99%, 0 risky pairs. bsai-hfv2: weakest-link fidelity 99.75%, mean fidelity 99.98%, 0 risky pairs.
• 20-overhang set (the 11 above plus 9 GetSet-added overhangs: ACCT, CCGC, ACAA, AACA, GAAA, CAAG, GCAC, TAGA, AAAT) — same source, Figure 4. This is the set the paper itself scores at roughly 80% aggregate fidelity using NEB's undisclosed formula; see the honesty-point section above for why SeqBench's methodology gives a different number by design. generic-t4-37c-1h: weakest-link fidelity 99.89%, mean fidelity 99.98%, 0 risky pairs. bsai-hfv2: weakest-link fidelity 99.74%, mean fidelity 99.95%, 0 risky pairs.
• CIDAR MoClo fusion sites — GGAG, TACT, AATG, AGGT, GCTT, CGCT, TGCC, ACTA. Source: Iverson SV, Haddock TL, Beal J, Densmore DM. "CIDAR MoClo: Improved MoClo Assembly Standard and New E. coli Part Library Enable Rapid Combinatorial Design for Synthetic and Traditional Biology." ACS Synth Biol. 2016;5(1):99-103. PMID: 26479688. This set was cross-validated via secondary sources (moclo.readthedocs.io and an Addgene kit guide); the original paper's own table was not independently re-verified against the primary source by SeqBench. generic-t4-37c-1h: weakest-link fidelity 99.95%, mean fidelity 99.99%, 0 risky pairs. bsai-hfv2: weakest-link fidelity 99.80%, mean fidelity 99.98%, 0 risky pairs.

The baseline this establishes: real published sets cluster at 99.7%+

Across three independent, peer-reviewed overhang sets, ranging from 8 to 20 overhangs, scored against two chemically distinct ligation datasets — six scores in total — every single one falls between 99.74% and 99.95% weakest-link fidelity, mean fidelities sit between 99.95% and 99.99%, and not one risky pair is flagged anywhere. That consistency, by itself, is useful information: it is an empirical baseline for what a competently designed, functionally validated Golden Gate/MoClo overhang set actually looks like under this methodology, independent of the specific ligase chemistry.

The practical rule of thumb that follows: if your own candidate overhang set scores meaningfully below roughly 99% weakest-link fidelity — and especially if it drops into the 90s or lower, or has any flagged risky pairs — treat that as a real, mechanistically grounded warning rather than noise. It is not a false-alarm-prone heuristic; it is the same measurement that puts every real published set tested here at 99.7% or above.

Two naive sets, worked in full

These two sets are constructed illustrative examples, not published or field-used sets — they exist purely to show how a naive, "looks fine to the eye" overhang choice degrades under real ligation data.

• Naive near-identical set: AAAA, AAAT, AATA, ATAA. Constructed as a single-base-shift pattern (each overhang looks like a one-position slide of the last) over a low-complexity, all-A/T background. generic-t4-37c-1h: weakest-link fidelity 96.39%, mean fidelity 98.78%, 0 risky pairs. bsai-hfv2: weakest-link fidelity 95.49%, mean fidelity 98.12%, 0 risky pairs.
• Naive single-base-swap set: GGAG, GGAT, GGAC, GGAA. Constructed by holding the first three bases fixed ("GGA") and varying only the last base — a pattern that looks systematic and readable on paper but gives the ligase very little distinguishing sequence to work with. generic-t4-37c-1h: weakest-link fidelity 87.53%, mean fidelity 95.50%, 1 risky pair flagged — GGAC↔GGAT, cross-reaction count 536. bsai-hfv2: weakest-link fidelity 85.44%, mean fidelity 88.27%, 5 risky pairs flagged — GGAC→GGAT (count 64), GGAT→GGAC (count 45), GGAT→GGAA (count 39), GGAA→GGAG (count 39), GGAG→GGAA (count 22).

What a flagged risky pair means for your actual assembly

A risky pair is not a similarity score or a computational guess — it is a real, measured ligation event count from an actual T4 ligase or BsaI-HFv2 reaction. Under bsai-hfv2 conditions, for example, GGAC's own row in the underlying published matrix shows 64 measured ligation events to a GGAT-presenting end; that count alone was large enough relative to GGAC's own correct-partner-plus-cross-reaction total to clear the 5% risk threshold, which is why the pair is flagged. GGAT→GGAC, GGAT→GGAA, GGAA→GGAG, and GGAG→GGAA each independently clear the same threshold in the same set — five distinct flagged directions among just four overhangs.

Mechanically, in a real one-pot Golden Gate reaction built on this four-overhang set, that means a fraction of the fragment ends presenting one of these overhangs will ligate to the wrong neighbor rather than failing to ligate at all. The resulting molecule still gets built and can still transform a cell; it is simply assembled at the wrong junction. That produces colonies that pass a crude "did anything grow" check but carry a misassembled construct, which then surfaces later as a failed functional screen or a sequence-verification mismatch rather than as an obvious ligation failure. This is the concrete, mechanical reason a low-diversity overhang choice is a design risk and not a stylistic nitpick: it changes measured molecular behavior, not just how the sequences look side by side.

It is also worth noting what the two naive examples show by contrast with each other: the near-identical AAAA-family set (differing from its neighbors by a one-base shift at varying positions) actually scores better than the GGAG-family set (differing from its neighbors only at the fixed last position), despite looking like the more repetitive, riskier-looking design on paper. Real ligation specificity is not a simple function of overall string similarity or edit distance — which position varies, and what the flanking bases are, changes measured outcomes in ways that are not obvious from eyeballing a list of four-letter codes. That is precisely why this needs to be checked against real ligation data rather than judged by inspection.

A workflow for checking your own overhangs before ordering fragments

This exact check is exposed directly as the golden_gate_fidelity tool on SeqBench's REST API and MCP server, so it can be called programmatically from a design script against a candidate overhang list rather than typed in by hand each time. A lighter, standalone check for basic overhang uniqueness and palindrome collisions is also available in the Golden Gate mode of SeqBench's Cloning & Assembly Simulator, useful as a quick sanity pass before running the fuller ligation-fidelity comparison described here.

1. Identify the real reaction chemistry you will actually use — a standard T4 ligation step after Type IIS digestion (score against generic-t4-37c-1h), or a genuine one-pot BsaI-HFv2 digestion-ligation (score against bsai-hfv2). Do not default to whichever dataset is listed first.
2. List the exact overhangs as you plan to order/synthesize them for one reaction — every junction, no accidental duplicates (a duplicate overhang used at two junctions in the same reaction is a design error the scorer will refuse to score, since it can't tell which fragment a repeated overhang belongs to).
3. Score the set and read the weakest-link fidelity first, not the mean. Compare it against the ~99.7%+ baseline established by the three published sets above.
4. If any risky pairs are flagged, treat the flagged list as a prioritized set of specific overhangs to reconsider, not a single pass/fail verdict on the whole design. Swap only the flagged overhangs where possible.
5. Re-score after each change. A set that clears both the weakest-link threshold and has zero flagged risky pairs against the chemistry you're actually using is a reasonable candidate to move forward with synthesizing fragments for.
6. If you must keep an overhang that shows some flagged cross-talk (for example because it is fixed by a required promoter or part-standard convention), make that a documented, deliberate trade-off rather than an unnoticed default.

Common mistakes and what this scoring does not tell you

• It does not reproduce NEB's or Potapov's own published aggregate fidelity percentages for named sets, as covered above — do not directly compare a weakest-link percentage from this methodology against a vendor-reported aggregate percentage for the same set.
• It only reflects the specific conditions actually measured in the underlying papers — T4 ligase alone at 1 hour/37°C, or real one-pot BsaI-HFv2 digestion-ligation. A different Type IIS enzyme's real one-pot chemistry, a different incubation time or temperature, a different buffer, or a different ligase entirely was not measured and is not modeled by either dataset.
• It only scores pairwise cross-talk among the overhangs you actually list. It does not model full combinatorial behavior in a shared reagent pool against all 256 possible overhang sequences, and it does not capture higher-order, multi-way ligation effects beyond the pairwise counts in the underlying matrices.
• It does not check anything else about your fragments: no scanning for internal Type IIS recognition sites, no GC-content or secondary-structure assessment of the overhang region itself, and no prediction of overall assembly efficiency, which depends heavily on fragment stoichiometry, backbone linearization/dropout efficiency, transformation efficiency, and colony-screening depth — none of which this scorer touches.
• The 5% risk threshold is a documented, adjustable default, not a universal biological cutoff. A ten-fragment assembly amplifies the practical cost of even one borderline overhang far more than a two-fragment assembly does, so tightening the threshold for larger assemblies is a reasonable choice, not a workaround.
• Zero flagged risky pairs does not mean zero mis-assembly risk. It means no single pairwise cross-reaction in the underlying published data cleared the chosen threshold; smaller cross-reaction counts below that threshold still exist in the raw matrices and are already folded into the reported mean and weakest-link fidelity numbers, just not called out individually as "risky."
• This is in vitro biochemical ligation data. It says nothing about how a correctly assembled construct behaves after transformation — plasmid toxicity, low-copy origin effects, or mutation during propagation are entirely outside its scope, and a perfect fidelity score is not a guarantee that your construct will work in vivo.

Related tools

Related guides

• Gibson Assembly vs. Golden Gate vs. Restriction Cloning: Choosing a Strategy →
• Cloning a Gene From Scratch: A Complete Workflow From Sequence to Verified Construct →
• Type IIS Enzymes for Golden Gate and MoClo Assembly: BsaI, BsmBI, BbsI, and SapI Compared →