Direct Binding Free Energy (DBFE)

• Direct Binding Free Energy (DBFE) is defined as the absolute free energy change when unbound molecules form a complex, central to understanding binding thermodynamics.
• The concept spans diverse methods including rigorous statistical mechanics, end-state non-alchemical approaches, and pathway-based calculations, each with necessary standard-state and restraint corrections.
• Applications include antibody-antigen prediction, protein-ligand binding studies, and ML-driven models, highlighting the importance of accurate solvent treatment and entropy corrections.

Searching arXiv for the cited works to ground the article in published sources. Direct Binding Free Energy (DBFE) denotes the absolute thermodynamic free energy change associated with forming a bound complex from its unbound components. Across the literature, the underlying quantity is usually the standard binding free energy, written as $\Delta G^\circ_{\mathrm{bind}}$ or $\Delta G$, and related to equilibrium association and dissociation constants by standard statistical thermodynamics. The label “DBFE” is not used uniformly: in some works it refers to direct end-state absolute binding free energy estimation without alchemical intermediates, in others to rigorous statistical-mechanical formulations of absolute binding, to path-based physical unbinding free energies, or to machine-learning prediction of experimentally measured $\Delta G$ values (Brocidiacono et al., 12 Mar 2026, Minh, 2012, Yu et al., 27 Aug 2025). This suggests that DBFE is best understood as a family resemblance term: methods differ substantially, but they aim at the same absolute binding thermodynamics rather than the relative mutation quantity $\Delta\Delta G$.

1. Thermodynamic definition and scope

The common thermodynamic object is the standard free energy for association. For a receptor and ligand, one formulation is

$$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$

while for protein-protein binding one paper writes

$$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$

Equivalent equilibrium relations are reported as

$\Delta G = -RT \ln K_a = RT \ln K_d$

and

$\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$

These expressions appear in antibody-antigen prediction, inverse-folding-based protein-protein binding analysis, and solution-phase free-energy treatments (Yu et al., 27 Aug 2025, Jiao et al., 2024, Jensen, 2015).

A central distinction in the literature is between absolute binding free energy and mutational free-energy differences. Relative mutational effects are defined as

$$\Delta\Delta G = \Delta G_{\mathrm{mut}} - \Delta G_{\mathrm{wt}},$$

or, in the protein-protein setting, as the difference between mutant and wild-type binding free energies. Several papers emphasize that DBFE corresponds to the absolute $\Delta G$ of a bound complex, whereas $\Delta G$0 is a derived relative quantity more suitable for mutation scans and affinity optimization (Yu et al., 27 Aug 2025, Jiao et al., 2024).

Standard-state treatment is integral to the definition. Multiple sources write the correction in the form

$\Delta G$1

or equivalent variants involving translational and rotational phase-space factors such as $\Delta G$2 or an effective binding-region volume $\Delta G$3 (Pal et al., 2019, Minh, 2012, Brocidiacono et al., 12 Mar 2026). In path-based formulations, omission of such terms is explicitly identified as the reason a reported free energy should be interpreted as qualitative or path-comparative rather than an absolute standard-state DBFE (Deeks et al., 2023).

2. Rigorous statistical-mechanical formulations

A rigorous formulation is provided by implicit ligand theory (ILT), which separates receptor sampling from ligand sampling while preserving exact expressions for the standard binding free energy. ILT defines an effective interaction energy $\Delta G$4 and a receptor-configuration-dependent binding potential of mean force $\Delta G$5, then expresses the standard binding free energy as

$\Delta G$6

with

$\Delta G$7

This is presented as a direct route to standard binding free energies from docking-like rigid-receptor calculations averaged over a receptor ensemble (Minh, 2012).

The same ILT perspective underlies later FFT-accelerated work, where the binding potential of mean force is estimated efficiently by evaluating rigid-body interaction energies over many ligand poses. For a fixed receptor conformation, the pose-dependent BPMF is written as

$\Delta G$8

and receptor-ensemble averaging yields the final standard binding free energy with the standard-state correction expressed through $\Delta G$9 (Nguyen et al., 2017). In that study, the resulting standard binding free energies for T4 lysozyme ligands are reported to agree well with previous alchemical calculations, with Pearson $\Delta G$0–$\Delta G$1 against a flexible-complex benchmark for 24 systems and $\Delta G$2 against an alternative rigid-receptor method for 141 ligands at fine grid spacing (Nguyen et al., 2017).

A conceptually different but still rigorous direct treatment appears in Monte Carlo integration. There, the gas-phase DBFE is obtained from a Boltzmann-weighted interaction enthalpy plus entropic penalties inferred from insertion statistics and standard-state volume restriction:

$\Delta G$3

The method is described as direct because it does not require alchemical transformations, restraints, or pose predefinition, although the reported implementation uses a rigid receptor and approximate solvation corrections (Clark et al., 2017).

These statistical-mechanical treatments share a common principle: DBFE can be written as a free-energy difference between physical bound and unbound states, but practical estimators differ in how they generate overlap, define the binding region, and treat solvent and configurational entropy.

3. Direct end-state versus alchemical DBFE

One major use of the term DBFE refers specifically to end-state absolute binding free energy methods that avoid alchemical intermediates. “Binding Free Energies without Alchemy” presents a method explicitly named Direct Binding Free Energy (DBFE), formulated in implicit solvent and based on only three simulations per ligand: apo receptor, free ligand, and holo complex (Brocidiacono et al., 12 Mar 2026). The method constructs a thermodynamic cycle with a decoupled unrestrained state, a decoupled restrained state, a clash-filtered decoupled state, and the coupled bound state:

$\Delta G$4

Its restraining potential combines a Gaussian translational term and a Bingham rotational term,

$\Delta G$5

with an analytic standard-state correction and an entropic clash-filtering term

$\Delta G$6

The coupled-versus-clash-filtered leg is estimated by MBAR (Brocidiacono et al., 12 Mar 2026).

In the host-guest benchmark reported there, DBFE in OBC2 implicit solvent achieves RMSE $\Delta G$7 kcal/mol, Pearson $\Delta G$8, and Spearman $\Delta G$9, outperforming OBC2 double decoupling in correlation and strongly outperforming OBC2 MM/GBSA, while TIP3P double decoupling remains best overall. On a protein-ligand benchmark, DBFE gives RMSE $\Delta\Delta G$0 kcal/mol, $\Delta\Delta G$1, and $\Delta\Delta G$2, which is better in RMSE than OBC2 double decoupling but below TIP3P explicit-solvent performance (Brocidiacono et al., 12 Mar 2026).

A related non-alchemical framework based on the energy representation (ER) theory of solution also computes binding free energies directly from endpoint ensembles. There the binding free energy is written as

$\Delta\Delta G$3

where the bound and dissociated solvation free energies are obtained from energy-distribution functionals rather than alchemical interpolation (Okita et al., 2024). For aspirin binding to $\Delta\Delta G$4-cyclodextrin, the ER-OR method reports $\Delta\Delta G$5 kcal mol$\Delta\Delta G$6 for the primary complex and $\Delta\Delta G$7 kcal mol$\Delta\Delta G$8 for the secondary complex, close to BAR values of $\Delta\Delta G$9 and $$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$0 kcal mol$$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$1, respectively (Okita et al., 2024).

By contrast, alchemical DBFE methods estimate the same absolute quantity through nonphysical intermediate Hamiltonians. The dopamine D3 receptor study uses a single-decoupling method with a $$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$2-dependent Hamiltonian

$$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$3

plus a standard-state correction

$$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$4

Its central finding is that enclosed hydration effects must be embedded into the solvation model: without HSA-informed hydration terms, RMSE is $$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$5 kcal/mol and correlation is $$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$6, whereas inclusion of enclosed hydration yields RMSE $$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$7 kcal/mol and $$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$8 (Pal et al., 2019).

This contrast clarifies a recurring ambiguity. In some publications, DBFE means “absolute binding free energy” irrespective of whether the estimator is direct or alchemical; in others, especially (Brocidiacono et al., 12 Mar 2026), it means an explicitly non-alchemical end-state method. The distinction is methodological, not thermodynamic.

4. Pathway-based and restraint-based direct calculations

A separate class of direct methods computes DBFE from physical binding or unbinding pathways. In the iMD-VR-FE framework, interactive molecular dynamics in virtual reality is used to generate ligand unbinding paths, which are then projected into a six-dimensional collective-variable space and profiled by umbrella sampling with WHAM (Deeks et al., 2023). The path collective variables are

$$\Delta G^\circ_{\mathrm{bind}} = G^\circ_{\mathrm{aq}}(RL) - G^\circ_{\mathrm{aq}}(R) - G^\circ_{\mathrm{aq}}(L),$$9

and

$$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$0

with a metric-corrected distance in CV space. The unbinding and binding free energies are then estimated from PMF end-state differences,

$$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$1

For benzamidine–trypsin, the reported path-based estimate is $$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$2 kcal mol$$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$3, substantially more favorable than the experimental value of approximately $$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$4 kcal mol$$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$5 because the implementation used implicit solvent, protein backbone restraints, and no explicit standard-state or restraint-release corrections (Deeks et al., 2023).

The paper is explicit that rigorous absolute DBFE from such PMF protocols requires standard-state and restraint corrections, including terms of the form

$$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$6

and, when orientational restraints are used, analytic correction formulas akin to the Boresch treatment (Deeks et al., 2023). This makes clear that pathway methods can produce either qualitative mechanistic free-energy profiles or quantitative DBFEs, depending on whether the full thermodynamic bookkeeping is included.

Attach–pull–release, umbrella sampling, funnel metadynamics, and related physical-coordinate approaches are described in comparative context in the SAMPL9 cyclodextrin study. That work notes that direct absolute routes can become difficult when binding involves large host reorganization and multiple competing poses, circumstances under which the authors chose the Alchemical Transfer Method instead of a direct PMF-style absolute route (Khuttan et al., 2023). This suggests that “directness” in DBFE is often limited by overlap and coordinate-design issues rather than by thermodynamic principle.

5. Machine-learning uses of DBFE

In machine-learning papers, DBFE is used as a prediction target rather than a simulated free energy. TopoBind models antibody-antigen binding free energy as a scalar thermodynamic quantity $$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$7 in kcal/mol and treats this as the direct, absolute binding free energy of the complex (Yu et al., 27 Aug 2025). Its curated dataset begins with 1705 antibody-antigen complexes from 472 unique PDB IDs and, after modality filtering, retains 303 PDB IDs corresponding to 1398 usable instances, split into 978/209/211 train/validation/test instances with no PDB ID overlap (Yu et al., 27 Aug 2025). The method fuses pooled ESM-2 embeddings with a 100-dimensional topological descriptor built from contact map metrics, interface contact statistics, distance map statistics, and persistent homology invariants up to dimension 2. Bidirectional cross-attention aligns sequence and topology, and the main model uses a Lasso regressor on the fused embedding (Yu et al., 27 Aug 2025).

The reported test performance for TopoBind with Lasso is MSE $$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$8, RMSE $$\Delta G_{\mathrm{bind}} = G_{AB}^{\mathrm{bound}} - \left(G_A^{\mathrm{unbound}} + G_B^{\mathrm{unbound}}\right).$$9, MAE $\Delta G = -RT \ln K_a = RT \ln K_d$0, Pearson $\Delta G = -RT \ln K_a = RT \ln K_d$1, $\Delta G = -RT \ln K_a = RT \ln K_d$2, and classification accuracy $\Delta G = -RT \ln K_a = RT \ln K_d$3 for strong binders defined by $\Delta G = -RT \ln K_a = RT \ln K_d$4 kcal/mol. The best topology hyperparameters are a contact threshold near $\Delta G = -RT \ln K_a = RT \ln K_d$5 Å and persistent-homology top-$\Delta G = -RT \ln K_a = RT \ln K_d$6, with performance described as robust near those choices (Yu et al., 27 Aug 2025).

A different ML framing appears in the Boltzmann-Aligned inverse folding model for protein-protein interactions. There, DBFE is approximated from inverse-folding log-likelihood differences between the bound complex and the decoupled monomers:

$\Delta G = -RT \ln K_a = RT \ln K_d$7

The explicit treatment of the unbound state is presented as the key physical inductive bias. On antibody-antigen binding energies in the SAbDab test, the resulting DBFE estimator achieves a Spearman correlation of $\Delta G = -RT \ln K_a = RT \ln K_d$8, outperforming DSMBind and other baselines according to the paper (Jiao et al., 2024).

These ML uses broaden the meaning of DBFE. The target remains the absolute $\Delta G = -RT \ln K_a = RT \ln K_d$9 of binding, but the estimator may be a supervised regressor trained on experimental labels or a Boltzmann-aligned proxy derived from sequence-conditioned structural likelihoods. A plausible implication is that, in current ML usage, DBFE often functions as the experimentally defined endpoint quantity against which surrogate models are calibrated rather than as a directly simulated thermodynamic observable.

6. Practical difficulties, corrections, and recurring controversies

Several recurring error sources cut across DBFE methodologies. Standard-state handling is one of the most universal. ILT, SDM, end-state DBFE, ER theory, and pathway approaches all require a correction converting a confined or restrained binding region into the $\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$0 M reference state (Minh, 2012, Pal et al., 2019, Brocidiacono et al., 12 Mar 2026, Okita et al., 2024). When such terms are omitted, as in the qualitative iMD-VR-FE implementation, the resulting free energy is explicitly not an absolute standard-state DBFE (Deeks et al., 2023).

Hydration and solvent representation are another major issue. The D3 receptor study shows that confined, energetically frustrated waters can dominate binding thermodynamics and that neglecting them can produce gross underestimation of affinity in implicit-solvent alchemical DBFE (Pal et al., 2019). The 2015 solution-phase review similarly identifies explicit solvent and ion effects, wrong ionization states, errors in ionic solvation free energies, molecular symmetry, anharmonicity, and insufficient conformational sampling as factors that can each contribute roughly $\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$1–$\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$2 kcal/mol errors at $\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$3 K (Jensen, 2015).

Multiple poses and chemical species complicate any absolute DBFE definition when the experimentally observed affinity is an aggregate over substates. The SAMPL9 cyclodextrin study treats pose-specific and species-specific binding constants through a combination formula,

$\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$4

equivalently,

$\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$5

For mixtures of species, the aggregation extends to

$\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$6

The paper’s conclusion is that apparently simple host-guest systems may require explicit treatment of multiple poses, enantiomers, protonation states, and slow conformational modes for reliable free-energy prediction (Khuttan et al., 2023).

Interpretability also varies strongly by method. Classical statistical-mechanical and end-state approaches can identify physically meaningful corrections such as clash entropy, hydration displacement, confinement volume, or reorganization free energy (Brocidiacono et al., 12 Mar 2026, Pal et al., 2019). ML approaches can be predictive but are only indirectly thermodynamic: sparse Lasso weights improve interpretability in TopoBind, yet cross-attention and topology encoders remain black-box components according to the paper (Yu et al., 27 Aug 2025).

A common misconception is that “direct” automatically means more physically rigorous. The literature does not support that equivalence. Some direct methods are fully rigorous when all corrections are included, as in ILT-based and ER-based formulations (Minh, 2012, Okita et al., 2024). Others are explicitly approximate or qualitative because of solvent models, rigid receptors, missing corrections, or proxy-based inference (Deeks et al., 2023, Clark et al., 2017, Jiao et al., 2024).

7. Comparative perspective

DBFE now spans several methodological regimes with a shared target but different computational philosophies.

Regime	Core idea	Representative paper
Rigorous direct statistical mechanics	Receptor/ligand sampling separated or insertion-based free-energy averaging	(Minh, 2012, Nguyen et al., 2017, Clark et al., 2017)
End-state non-alchemical ABFE	Bound and decoupled end states connected without $\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$7 intermediates	(Brocidiacono et al., 12 Mar 2026, Okita et al., 2024)
Alchemical absolute binding free energy	$\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$8-dependent Hamiltonians with standard-state and restraint corrections	(Pal et al., 2019, Khuttan et al., 2023)
Pathway/PMF direct methods	Physical unbinding pathways profiled along CVs	(Deeks et al., 2023)
Machine-learning DBFE prediction	Experimental or proxy $\Delta G^\circ_{\mathrm{bind}} = -k_B T \ln K_a.$9 predicted from structure/sequence representations	(Yu et al., 27 Aug 2025, Jiao et al., 2024)

Taken together, the literature shows that DBFE is less a single algorithm than a unifying thermodynamic objective. In the narrowest modern sense, especially in “Binding Free Energies without Alchemy,” DBFE denotes an end-state absolute binding free energy method that avoids alchemical intermediates (Brocidiacono et al., 12 Mar 2026). In the broader sense used across statistical mechanics, simulation, and ML, it denotes the absolute binding free energy $$\Delta\Delta G = \Delta G_{\mathrm{mut}} - \Delta G_{\mathrm{wt}},$$0 itself—distinguished from $$\Delta\Delta G = \Delta G_{\mathrm{mut}} - \Delta G_{\mathrm{wt}},$$1 and from qualitative binding scores—regardless of whether it is obtained by rigorous ensemble averages, physical-pathway PMFs, alchemical cycles, or predictive surrogates (Yu et al., 27 Aug 2025, Jiao et al., 2024, Minh, 2012).