How to Run Boltz-2

by Corin Wagen · Jun 6, 2025

Today, a team of researchers from MIT and Recursion released Boltz-2, an open-source protein–ligand co-folding model. Boltz-2 can not only predict the structure of biomolecular complexes from sequences, it also "approaches the accuracy of FEP-based methods" at protein–ligand binding-affinity prediction (source).

In a set of binding-affinity benchmarks, the authors show that Boltz-2 performs almost as well as the industry-standard FEP+ workflow and handily outperforms cheaper physics-based methods like MM/PBSA, although performance is considerably worse on internal targets from Recursion:

Comparison of Boltz-2 to other methods.

Figure 6 from the Boltz-2 paper.

While full assessment of Boltz-2's capabilities will require extensive benchmarking and external verification, it's already possible for scientists to start using Boltz-2 for their own projects. In this post, we provide step-by-step guides on how to run Boltz-2 locally and through Rowan's computational-chemistry platform.

(Curious about how Boltz-2 works? Check out this FAQ to learn more about what the model's trained on, where it can be useful, and where it still has limitations.)

Running Locally

1. Install Boltz-2

Boltz-2 is an open-source model and can be installed from the Python Package Index. You can install this any number of ways; we like using pixi for dependency management.

pixi init
pixi add python=3.12
pixi add --pypi boltz

2. Create a Template .yaml File

Boltz-2 requires a specific input-file syntax. The authors provide several examples in their GitHub repository; here's the example .yaml file for predicting protein–ligand binding affinity.

version: 1  # Optional, defaults to 1
sequences:
  - protein:
      id: A
      sequence: MVTPEGNVSLVDESLLVGVTDEDRAVRSAHQFYERLIGLWAPAVMEAAHELGVFAALAEAPADSGELARRLDCDARAMRVLLDALYAYDVIDRIHDTNGFRYLLSAEARECLLPGTLFSLVGKFMHDINVAWPAWRNLAEVVRHGARDTSGAESPNGIAQEDYESLVGGINFWAPPIVTTLSRKLRASGRSGDATASVLDVGCGTGLYSQLLLREFPRWTATGLDVERIATLANAQALRLGVEERFATRAGDFWRGGWGTGYDLVLFANIFHLQTPASAVRLMRHAAACLAPDGLVAVVDQIVDADREPKTPQDRFALLFAASMTNTGGGDAYTFQEYEEWFTAAGLQRIETLDTPMHRILLARRATEPSAVPEGQASENLYFQ
  - ligand:
      id: B
      smiles: 'N[C@@H](Cc1ccc(O)cc1)C(=O)O'
properties:
  - affinity:
      binder: B

Boltz-2 does not yet support protein–protein binding affinity or predicting binding affinity for multiple ligands.

3. Run Boltz-2 Locally

To run Boltz-2, initialize the environment shell and then use boltz predict to run the model.

pixi shell
boltz predict affinity.yaml --use_msa_server

This call will take a little while to run; make sure your computer has enough disk space to download the model weights! When finished, Boltz-2 will write a bunch of directories and .json files containing predictions. The predictions will be located in output/predictions/[input-file]/affinity-[input-file].json, and will contain predicted IC50 values (in micromolar) and binary probability that the compound is a binder.

Boltz-2 is a complex package and this guide barely scratches the surface. For a full guide to running prediction with Boltz-2, see the authors' documentation.

Running Through Rowan

To quickly use Boltz-2 for binding-affinity prediction, calculations can also be run through Rowan. Creating an account on Rowan is completely free and can be done using any Google-managed email account; create an account here.

1. Choose Workflow

Selecting the protein-ligand co-Folding workflow.

Once you sign in to Rowan, you can select which workflow you want to run. Here, we'll select the "Protein–Ligand Co-Folding" workflow (towards the bottom of the screen).

2. Enter Protein and Ligand

Inputting a sucrose molecule.

Proteins can be specified by sequence; existing protein structures in Rowan won't work, because this is co-folding—we don't want to start with a 3D structure.

Molecules can be loaded into Rowan by name, by SMILES, by input file, or through our provided 2D and 3D editors. Here, we'll input the molecule from the above demo by SMILES.

3. Run Calculation

The finished calculation results.

Once you click "Submit Calculation," we'll allocate a cloud GPU and start running Boltz-2 on your system. The calculations should be done in a few minutes and can be viewed through the browser.

Rowan displays the predicted protein–ligand complex through our 3D viewer, with predicted binding affinity and confidence metrics on the side. The complex can be downloaded as a PDB file for further analysis.

Banner background image

What to Read Next

Studying Scaling in Electron-Affinity Predictions

Studying Scaling in Electron-Affinity Predictions

Testing low-cost computational methods to see if they get the expected scaling effects right.
Sep 10, 2025 · Corin Wagen
Open-Source Projects We Wish Existed

Open-Source Projects We Wish Existed

The lacunæ we've identified in computational chemistry and suggestions for future work.
Sep 9, 2025 · Corin Wagen, Jonathon Vandezande, Ari Wagen, and Eli Mann
How to Make a Great Open-Source Scientific Project

How to Make a Great Open-Source Scientific Project

Guidelines for building great open-source scientific-software projects.
Sep 9, 2025 · Jonathon Vandezande
ML Models for Aqueous Solubility, NNP-Predicted Redox Potentials, and More

ML Models for Aqueous Solubility, NNP-Predicted Redox Potentials, and More

the promise & peril of solubility prediction; our approach and models; pH-dependent solubility; testing NNPs for redox potentials; benchmarking opt. methods + NNPs; an FSM case study; intern farewell
Sep 5, 2025 · Eli Mann, Corin Wagen, and Ari Wagen
Machine-Learning Methods for pH-Dependent Aqueous-Solubility Prediction

Machine-Learning Methods for pH-Dependent Aqueous-Solubility Prediction

Prediction of aqueous solubility for unseen organic molecules remains an outstanding and important challenge in computational drug design.
Sep 5, 2025 · Elias L. Mann, Corin C. Wagen
What Isaiah and Sawyer Learned This Summer

What Isaiah and Sawyer Learned This Summer

Reflections from our other two interns on their time at Rowan and what they learned.
Sep 5, 2025 · Isaiah Sippel and Sawyer VanZanten
Benchmarking OMol25-Trained Models on Experimental Reduction-Potential and Electron-Affinity Data

Benchmarking OMol25-Trained Models on Experimental Reduction-Potential and Electron-Affinity Data

We evaluate the ability of neural network potentials (NNPs) trained on OMol25 to predict experimental reduction-potential and electron-affinity values for a variety of main-group and organometallic species.
Sep 4, 2025 · Sawyer VanZanten, Corin C. Wagen
Which Optimizer Should You Use With NNPs?

Which Optimizer Should You Use With NNPs?

The results of optimizing 25 drug-like molecules with each combination of four optimizers (Sella, geomeTRIC, and ASE's implementations of FIRE and L-BFGS) and four NNPs (OrbMol, OMol25's eSEN Conserving Small, AIMNet2, and Egret-1) & GFN2-xTB.
Sep 4, 2025 · Ari Wagen and Corin Wagen
Double-Ended TS Search and the Invisible Work of Computer-Assisted Drug Design

Double-Ended TS Search and the Invisible Work of Computer-Assisted Drug Design

finding transition states; the freezing-string method; using Rowan to find cool transition states; discussing drug design
Sep 3, 2025 · Jonathon Vandezande, Ari Wagen, Spencer Schneider, and Corin Wagen
The Invisible Work of Computer-Assisted Drug Design

The Invisible Work of Computer-Assisted Drug Design

Everything that happens before the actual designing of drugs, and how Rowan tries to help.
Aug 28, 2025 · Corin Wagen