Boltz-2 predicts the 3D structure of biomolecular complexes and, for protein–small-molecule systems, estimates whether the ligand binds and how strong that binding may be. You can run it online without installing Python, downloading model weights, or configuring a GPU.
To use Boltz-2 online:
- Open the ProteinIQ Boltz-2 webserver.
- Add the protein, ligand, DNA, or RNA chains in your complex.
- Optionally provide an MSA, structural template, or experimental constraints.
- Keep the default sampling settings for your first run.
- Submit the job.
- Inspect the predicted structure before interpreting confidence or affinity scores.
- Download the structure and supporting data for further analysis.
ProteinIQ accepts sequences, molecular structures, database identifiers, and common ligand formats, then displays the completed prediction in an interactive structure viewer alongside confidence, affinity, and downloadable output files.[1]
Boltz-2 can answer several different questions, but its outputs are not interchangeable. Use the predicted coordinates to study how a complex might assemble, the binary affinity score to separate likely binders from decoys, and the continuous affinity value to compare related active compounds. None of these outputs should be treated as an experimental binding measurement.[2]
What is Boltz-2?
Boltz-2 is an open-source biomolecular structure prediction model developed by researchers at MIT, Recursion, Valence Labs, and collaborating institutions. It was introduced in 2025 as the successor to Boltz-1. Like other all-atom cofolding models, it can predict complexes containing proteins, small molecules, DNA, and RNA.[3]
Its distinguishing feature is a separate affinity model for protein–small-molecule interactions. A single run can therefore produce both a proposed complex structure and two affinity outputs:
| Output | What it answers |
|---|---|
| Binding probability | Is this molecule more likely to behave like a binder or a decoy? |
| Affinity value | How does this molecule compare with other active compounds for the same target and assay context? |
The original paper reported an average Pearson correlation of 0.66 on a four-target subset of an FEP benchmark, with inference more than 1,000 times faster than the compared free-energy calculations. Results varied substantially outside that benchmark. On eight blinded internal assays, correlations exceeded 0.55 on three assays and were lower on the remaining five. The paper is a preprint, so these numbers should be treated as reported model results rather than a universal estimate of real-world accuracy.[3]
Boltz-2 is therefore most useful as a hypothesis-generation and prioritization model. It can suggest a pose, highlight uncertain regions, and help rank compounds before more expensive calculations or experiments.
How does Boltz-2 work?
Boltz-2 combines a shared molecular representation, an all-atom diffusion model, and separate confidence and affinity modules.
The trunk represents the complete complex
Boltz-2 first converts the input into molecular tokens. Protein residues, nucleotides, ligands, templates, multiple sequence alignments, and constraints are represented together.
A neural network called the trunk repeatedly updates two kinds of information:
- A representation of each individual token
- A pairwise representation describing how every token may relate to every other token
This lets the model reason about relationships within a protein, between protein chains, and across protein–ligand or protein–nucleic-acid interfaces. The trunk is based on PairFormer-style blocks and can incorporate evolutionary information from an MSA when one is provided.[3]
A diffusion model generates atomic coordinates
The structure module begins with noisy atomic coordinates and gradually denoises them into a complete three-dimensional structure. This process resembles diffusion-based image generation, except the model is generating atom positions rather than pixels.
Boltz-2 normally performs 200 diffusion sampling steps and three recycling iterations. Recycling sends the model's intermediate representation through the network again so it can refine the predicted geometry. These defaults are sensible for most first runs.[2]
Boltz-2 can also apply inference-time potentials. These potentials push sampled structures toward user-defined pocket or contact constraints. They can be useful when an experiment, homologous structure, or mutational study already identifies the likely interface.
The confidence module scores the prediction
After generating a structure, Boltz-2 estimates several forms of uncertainty:
pLDDTmeasures confidence in local atomic or residue-level geometry.pTMestimates confidence in the overall fold.ipTMfocuses on the arrangement of separate chains or molecular components.PAEestimates uncertainty in the relative positions of pairs of residues or tokens.PDEestimates pairwise distance error in ångströms.
Most confidence scores run from 0 to 1, with higher values indicating greater model confidence. PDE is measured as an error, so lower values are better. Confidence is not the same as correctness, but it helps identify which regions and interfaces deserve the most scrutiny.[2]
The affinity module evaluates protein–ligand binding
The affinity module receives the predicted structure and pairwise molecular representation. It then uses two prediction heads:
- A binary head trained to distinguish binders from decoys
- A continuous head trained to predict relative binding strength
The continuous model was trained using several experimental endpoint types, including IC50, Ki, Kd, AC50, and related measurements, standardized onto a common logarithmic scale. Its result should therefore be read as an IC50-like estimate of binding strength, not a literal prediction of the exact endpoint that would be measured in a particular assay.[3]
How to use Boltz-2 online
1. Open the Boltz-2 webserver
Go to Boltz-2 online.
The form lets you assemble a complex one molecule at a time. Each component receives a chain ID such as A, B, or C. Keep an eye on these IDs because templates, covalent bonds, pocket definitions, and contact constraints refer to them.
A job must contain at least one protein, DNA, or RNA chain. A ligand by itself is not a foldable system.
2. Add the molecules in your complex
Choose the input type that matches each component:
| Component | Accepted input | Current ProteinIQ limit |
|---|---|---|
| Protein | FASTA, .fa, PDB, or an RCSB chain | 10 chains |
| Ligand | SMILES, SDF, MOL, MOL2, or PubChem | 10 ligands |
| Ligand from the PDB Chemical Component Dictionary | CCD code such as ATP, NAD, HEM, or SAH | 10 ligands |
| DNA | FASTA or plain nucleotide sequence | 10 chains |
| RNA | FASTA or plain nucleotide sequence | 10 chains |
| Template | PDB, mmCIF, or an RCSB structure | 5 templates |
| Precomputed protein MSA | A3M | 10 alignments |
ProteinIQ currently accepts complexes containing up to 5,000 residues in total. Memory use grows quickly with complex size, so removing irrelevant chains or long disordered tails can make large jobs more manageable.[1]
Paste protein sequences using the one-letter amino-acid alphabet. Headers are optional when entering a single sequence, but valid FASTA is easier to reuse and audit. The TXT to FASTA converter can clean up unformatted sequence lists before submission.
When starting from a PDB entry, specify the exact biological chain rather than assuming every chain in the deposited structure belongs in the prediction.
3. Add the ligand correctly
For a custom compound, paste a SMILES string or upload an SDF, MOL, or MOL2 file. You can also retrieve a compound from PubChem.
Use a CCD code when the ligand is a standardized PDB component such as ATP, NAD, HEM, or an existing covalent inhibitor. CCD input is especially important for covalent-bond constraints because Boltz needs standardized atom names from the component definition.
For affinity prediction, keep the intended binder unambiguous. The official Boltz-2 implementation supports one small-molecule affinity binder per prediction, and the affinity model is designed for a ligand binding to a protein. It does not predict protein–protein, protein–DNA, protein–RNA, or ligand–RNA binding affinity.[2]
Ligand size also matters. The official implementation rejects affinity ligands with more than 128 non-hydrogen atoms and advises against using molecules much larger than 56 atoms, which was the upper limit represented during affinity training. Larger molecules may still be structurally included in some workflows, but their affinity scores are outside the model's recommended range.[2]
For peptide ligands, add the peptide as a protein chain rather than encoding it as a large small molecule.
4. Decide whether to generate an MSA
A multiple sequence alignment gives the model evolutionary information about which residues tend to vary together.
ProteinIQ can generate an MSA through ColabFold. The available search depths are:
| Setting | Maximum sequences | When to use it |
|---|---|---|
| Shallow | 2,048 | Fast exploratory runs |
| Normal | 8,192 | Default choice for most proteins |
| Deep | 16,384 | Difficult targets where added search time is acceptable |
You can also provide a precomputed A3M alignment. When several protein chains are present, alignments are assigned in protein-chain submission order.
For a natural protein with known homologs, an MSA is usually worth testing when prediction quality matters. Single-sequence mode is faster, but the Boltz documentation warns that omitting an MSA can reduce accuracy. For designed proteins, orphan proteins, or a quick input check, an MSA may add little value.[2]
A practical workflow is to run one inexpensive single-sequence prediction first. Once the input and chain arrangement look correct, repeat the job with a normal-depth MSA and compare the result.
5. Add templates only when they provide useful information
A structural template can guide Boltz-2 toward a known fold or conformational state. Upload a PDB or mmCIF file, or retrieve a structure from RCSB.
Templates are useful when:
- A close homolog has already been solved
- You need a particular domain orientation
- The protein has several known conformational states
- You want to preserve an experimentally supported backbone arrangement
Map each template to the intended query chain. Do not assume that a template automatically applies to every chain in the job.
A poor template can bias the result in the wrong direction. Run an untemplated prediction as a control whenever the template's relevance is uncertain.
6. Add pocket, contact, or covalent constraints when justified
Constraints let you provide information that Boltz-2 would otherwise have to infer.
A pocket constraint identifies a ligand and the residues expected to form its binding site:
C|A:45,A:46|6.0|trueThis example asks ligand chain C to remain near residues 45 and 46 of protein chain A, with a maximum distance of 6 Å.
A residue contact constraint connects two positions:
A:10,B:5|8.0|trueA covalent-bond constraint identifies the exact atoms to connect:
A:12:SG,B:1:C22This example connects the sulfur atom of cysteine 12 in protein chain A to atom C22 in CCD ligand B.
Pocket and contact entries marked force=true become strict steering restraints only when Use constraint potentials is enabled. Without the potential, the model receives the constraint but may not satisfy it closely. Even forced constraints must be checked in the completed structure rather than assumed to have worked.[2]
Custom SMILES, SDF, MOL, and MOL2 ligands cannot currently be used for covalent atom constraints in the ProteinIQ interface. Use a CCD component whose standardized atom names match the RCSB component definition.
7. Choose the sampling settings
For the first meaningful run, change as little as possible.
| Setting | Recommended starting value | When to change it |
|---|---|---|
| Number of samples | 1 for input testing, then 3–5 | Increase when poses or interfaces are uncertain |
| Recycling steps | 3 | Increase only when testing whether extra refinement helps |
| Sampling steps | 200 | Increase for a deliberate quality-versus-runtime comparison |
| Step scale | 1.5 | Lower values produce more diverse samples |
| Affinity sampling steps | 200 | Keep at default unless benchmarking |
| Output format | CIF | Use PDB only for software that cannot read CIF |
| Save PAE matrix | Off | Enable for domain and interface uncertainty analysis |
| Method conditioning | None | Set only when targeting a specific experimental structural style |
The official defaults are three recycling steps, 200 structure sampling steps, a Boltz-2 step scale of 1.5, 200 affinity sampling steps, and five internal affinity samples. CIF is the default structure format.[2]
More samples are usually more useful than immediately increasing the number of diffusion steps. Multiple samples let you see whether the same fold, interface, or ligand pose appears consistently.
Leave confidence filtering disabled for the first run. A low-confidence sample can still reveal alternative conformations or explain why the model is uncertain.
8. Submit the job
Give the run a descriptive name that identifies the target, ligand, and important settings. For example:
KRAS_G12C_U4U_MSA_5samplesReview the chain order and ligand identity one last time, then submit the job.
The completed result opens with three main views:
- Viewer: Interactive inspection of predicted structures
- Data: Confidence and affinity results in tabular form
- Files: Structures, JSON data, matrices, and supporting files
All result files can be downloaded for downstream analysis.[1]
How to interpret Boltz-2 results
Inspect the structure first
Start with the highest-ranked sample, but also compare the alternatives.
Check:
- Whether the expected domains are folded
- Whether separate chains form a plausible interface
- Whether the ligand occupies a chemically sensible pocket
- Whether covalent or contact constraints were satisfied
- Whether the ligand has the expected stereochemistry and geometry
- Whether different samples converge on the same pose
A high confidence score cannot rescue a chemically impossible ligand pose. Visual inspection and basic structural validation should come before affinity interpretation.
For ligand complexes, you can extract or convert the ligand into a supported molecular format and run PoseBusters to check bond geometry, clashes, stereochemistry, and other physical plausibility issues.
Read confidence at the right level
pLDDT is mainly a local confidence measure. It can identify uncertain loops or flexible termini, but it does not by itself prove that two chains are arranged correctly.
For multicomponent complexes, pay particular attention to ipTM and the cross-chain regions of the PAE matrix. A protein can have confident local folds while the relative placement of its chains or ligand remains uncertain.
The combined confidence_score is useful for ranking samples from the same job. It should not be treated as a calibrated probability that the entire structure is correct.[2]
Use affinity probability for hit discovery
affinity_probability_binary estimates how binder-like the ligand appears relative to decoys.
Use it when screening diverse molecules and asking which compounds are worth investigating further. Higher values indicate greater binder-like probability, but the official documentation does not define one universal cutoff that works across all targets.
Calibrate the score using known binders and presumed non-binders for your own target. A threshold that works for one protein family or assay may perform poorly for another.[2]
Use affinity value for related active compounds
affinity_pred_value is intended primarily for comparing compounds that are already active, particularly within a chemically related series.
The output is:
log10(IC50 in μM)Lower values indicate stronger predicted binding. Convert the output back to micromolar units with:
IC50 (μM) = 10 ^ affinity_pred_value| Boltz-2 value | Converted concentration |
|---|---|
| -3 | 0.001 μM = 1 nM |
| -2 | 0.01 μM = 10 nM |
| -1 | 0.1 μM = 100 nM |
| 0 | 1 μM |
| 1 | 10 μM |
| 2 | 100 μM |
For example, a predicted value of -1.5 corresponds to approximately 0.032 μM, or 32 nM.
This conversion does not make the result an experimental IC50. Because the training data combines several measurement types and assay conditions, use the score for relative prioritization within a consistent target and compound series. The official documentation specifically discourages using the continuous value to rank mixtures of active and inactive molecules.[2]
Boltz-2 examples
KRAS G12C with a covalently attached ligand
This example predicts human KRAS G12C with the U4U ligand from the PDB Chemical Component Dictionary. It demonstrates a protein–ligand run with an explicit covalent bond between the mutated cysteine and a named ligand atom.
- Inputs: one 168-residue KRAS G12C protein chain and U4U as a CCD ligand
- Non-default settings:
Number of samples=5to compare repeated predictions;Generate MSA= on withNormaldepth to add evolutionary context;Step scale=1.638to preserve the submitted sampling configuration;Use constraint potentials= on;Covalent bonds=A:12:SG,B:1:C22
The Structure view keeps the U4U ligand and the surrounding KRAS fold visible while retaining all five ranked predictions. The covalent-bond definition tells Boltz-2 which atoms should be connected; the image does not establish that the predicted geometry or binding mode is experimentally correct.
All five predictions return an Affinity Pred Value of −0.401 log10(IC50) μM and an Affinity Probability Binary of 0.950. Their average Complex Plddt is 0.960, with a standard deviation of 0.002. This consistency is useful for inspecting the run, but predicted affinity and confidence are model outputs rather than experimental measurements of potency or covalent engagement.
Transcription-factor dimer on a DNA duplex
This four-chain example combines two complementary 15-nucleotide DNA strands with two copies of the same transcription-factor construct. It demonstrates joint prediction of a protein dimer and a sequence-specific DNA duplex.
- Inputs: DNA strands
TGGGTCACGTGTTCCandAGGAACACGTGACCC, plus two identical 89-residue protein chains - Non-default settings:
Number of samples=5to assess prediction consistency;Generate MSA= on withNormaldepth for the protein chains;Step scale=1.638to preserve the submitted sampling configuration
The Structure view shows the two protein chains arranged around the predicted DNA duplex. It is useful for examining the overall assembly and whether the nucleic acid remains duplex-like, but it does not prove sequence-specific recognition or identify experimentally validated contacts.
The top-ranked prediction has pTM = 0.912 and Iptm = 0.913. Across all five predictions, average Complex Plddt is 0.966 and average Complex Iplddt is 0.979, indicating consistent confidence in the returned assembly. Because this job has no ligand, the displayed affinity value 0.000 and probability 0.500 are not protein–ligand binding results and should not be interpreted as affinity measurements.
Human U1A protein bound to a U1 snRNA hairpin
This example uses the 97-residue human U1A RNA-recognition domain and the 21-nucleotide U1 snRNA hairpin from PDB 1URN. It adds a compact protein–RNA workflow to the protein–ligand and protein–DNA examples above.
- Inputs: U1A protein sequence from PDB 1URN and RNA sequence
AAUCCAUUGCACUCCGGAUUU - Non-default settings: none; the job uses one structure sample, no generated MSA, no constraints, and CIF output
The Structure view shows a compact protein–RNA assembly for the single returned prediction. It demonstrates that Boltz-2 accepts RNA as a first-class polymer input rather than treating it as a small-molecule ligand.
The result reports Confidence Score = 0.934, pTM = 0.935, Iptm = 0.881, Complex Plddt = 0.947, and Complex Iplddt = 0.953. These values support inspection of the predicted overall fold and interface, but they do not validate the exact base contacts, binding specificity, or agreement with the experimental 1URN structure.
What can you use Boltz-2 for?
Predicting protein–ligand complexes from sequence
Boltz-2 can begin with a protein sequence and ligand rather than a prepared receptor structure. This makes it useful when no experimental structure is available.
The result can provide an initial binding-pose hypothesis for later docking, molecular dynamics, free-energy calculations, mutagenesis planning, or experimental structure determination.
Ranking a lead series
For a group of related active molecules tested against the same target, the continuous affinity output can help prioritize compounds before more expensive calculations.
Keep the protein sequence, templates, constraints, and model settings constant across the series. Changing the structural context between compounds introduces variation that can obscure the comparison.
Screening for likely binders
The binary affinity probability is better suited to broad hit discovery than the continuous value. Use known positive and negative controls to assess whether the score separates compounds on your target before expanding to a larger virtual screening library.
Predicting protein–protein and protein–nucleic-acid structures
Boltz-2 can predict complexes containing multiple proteins, DNA, or RNA. This includes transcription factor–DNA complexes, protein–RNA complexes, and assemblies containing several biomolecule types.
These jobs produce structures and confidence estimates, but not protein–protein or nucleic-acid binding affinity.[2]
Modeling covalent complexes
When the ligand exists in the PDB Chemical Component Dictionary, Boltz-2 can model a specified covalent connection between a protein residue and ligand atom.
This is useful for covalent inhibitors, cofactors, and other systems where the bond is already known. It does not establish that a covalent reaction will occur. The bond is an input assumption.
What are the limitations of Boltz-2?
Boltz-2 is fast enough to test many hypotheses, but its affinity output has a narrower scope than its structure model.
The main limitations are:
- Affinity prediction currently supports one small molecule binding to a protein.
- Ligands much larger than 56 non-hydrogen atoms are outside the recommended training range.
- Affinity prediction depends on the model generating the correct pocket, pose, protonation context, and protein state.
- The affinity module does not explicitly model all effects from waters, ions, cofactors, multimeric partners, or assay conditions.
- Large ligand-induced conformational changes can be missed.
- A confident structure does not guarantee an accurate affinity estimate.
- The original Boltz-2 paper remains a preprint.[3]
Later evaluations reinforce the need for target-specific calibration. A 2026 peer-reviewed benchmark covering 356 targets reported a mean absolute error of roughly 0.9 on its standardized affinity scale and characterized performance as useful but moderate. A separate target-specific preprint found weak-to-moderate global agreement and poor overlap among the top-ranked compounds for two large datasets. These studies used different datasets and reference methods, so they should not be collapsed into one accuracy number.[4]
The practical rule is simple: test Boltz-2 against compounds with known outcomes on your target before trusting its ranking on unknown compounds.
Boltz-2 alternatives
The best alternative depends on what information you already have.
| Tool | Choose it when |
|---|---|
| OpenFold-3 | You want an independent all-atom complex structure prediction |
| Chai-1 | You want to cross-check a biomolecular complex with another cofolding model |
| DiffDock-L | You already have a receptor structure and want blind ligand-pose prediction |
| GNINA | You know the binding region and want docking with convolutional-neural-network scoring |
| AutoDock Vina | You want a conventional, configurable docking search within a defined box |
| LMI4Boltz | Your main concern is running Boltz-style predictions on large complexes with lower memory requirements |
Boltz-2 differs from conventional docking because it can generate the receptor and complex together from sequence. This is useful without an experimental structure, but it also means that errors in the predicted protein conformation can propagate into the ligand pose and affinity score.
When a high-quality receptor structure and known pocket are available, AutoDock Vina or GNINA provides more direct control over the search region. DiffDock-L is useful when the pocket is unknown but the receptor structure is already prepared.[5]
For important decisions, agreement between structurally different methods is more informative than a single high score. A ligand pose supported by Boltz-2, a docking model, structural constraints, and experimental structure–activity relationships is more credible than one supported by Boltz-2 alone.
Frequently asked questions
Can I use Boltz-2 without coding?
Yes. The ProteinIQ Boltz-2 webserver exposes sequence, ligand, MSA, template, constraint, sampling, and output settings through a browser interface. You do not need to install the command-line package or manage GPU infrastructure.
Does Boltz-2 need a protein structure?
No. You can provide a protein sequence and let Boltz-2 predict the protein and ligand complex together.
A known structure can still be useful as a template. Conventional docking tools generally require a prepared receptor structure, while Boltz-2 can begin from sequence.
Does Boltz-2 need an MSA?
No, but the official documentation warns that single-sequence mode can reduce accuracy. Start without an MSA when checking an input or working with a designed protein, then compare the result with an MSA-enabled run when suitable homologs exist.[2]
Can Boltz-2 predict protein–protein binding affinity?
No. It can predict the structure of a protein–protein complex, but its affinity module is designed for a small molecule binding to a protein.[2]
Is Boltz-2 better than AlphaFold 3?
Not across every structure prediction task. In the Boltz-2 paper, the model was competitive with other cofolding systems but remained behind AlphaFold 3 overall on the authors' structure benchmark, including antibody–antigen prediction.
Boltz-2's main distinction is its open implementation and integrated protein–small-molecule affinity module, not universal structural superiority.[3]
How many samples should I generate?
Use one sample to verify that the input runs correctly. For a real prediction, three to five samples provide a basic view of structural consistency without making the job unnecessarily large.
Increase to 10 or 20 when alternative poses, flexible interfaces, or constraint satisfaction are central to the question. Compare whether the samples converge rather than selecting the highest score without inspection.
A practical first-run recipe
For a first protein–ligand prediction:
- Add one protein sequence.
- Add one ligand as SMILES or SDF.
- Leave templates and constraints empty.
- Run one sample without an MSA to verify the input.
- Repeat with a normal-depth MSA and five samples.
- Keep three recycling steps and 200 sampling steps.
- Download CIF output and the confidence data.
- Inspect the ligand pose and interface.
- Use affinity probability for binder classification.
- Use affinity value only when comparing related active molecules.
Run the first job with Boltz-2 online, then add complexity only when the result or scientific question gives you a reason to do so.
Sources▼
- Use Boltz-2 Online ProteinIQ · August 16, 2026. https://proteiniq.io/app/boltz-2
- boltz/docs/prediction.md GitHub (jwohlwend/boltz) · August 16, 2026. https://github.com/jwohlwend/boltz/blob/main/docs/prediction.md
- Boltz-2: Towards Accurate and Efficient Binding Affinity Prediction bioRxiv · 2025. https://doi.org/10.1101/2025.06.14.659707
- ChEMBL-Derived Benchmark Dataset and Computational Results of Boltz-2-Based Binding Affinity Prediction Chem-Bio Informatics Journal · 2026. https://doi.org/10.1273/cbij.26.11
- Use AutoDock Vina Online ProteinIQ · August 16, 2026. https://proteiniq.io/app/autodock-vina

Founder and computational chemist, ProteinIQ
Dr. Matic Broz is the founder of ProteinIQ and a computational chemist. He completed a PhD focused on protein structure, molecular dynamics, and neural networks, and writes about structural biology and scientific software.