REVIEW 2 major objections 1 minor
Approximate calculation of the binding energy between 17$\beta$-estradiol and human estrogen receptor alpha
T0 review · 2 major / 1 minor · reviewed 2026-05-24 · grok-4.3
Pith's one-line read Multilayer fragment molecular orbital calculations estimate the binding energy of 17β-estradiol to estrogen receptor alpha at -67.2 kcal/mol.
desk verdict This paper reports one MFMO binding energy of -67.2 kcal/mol for estradiol-ER but the MD cluster sampling and water-shell subtraction lack the checks needed to trust the difference. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Multilayer fragment molecular orbital (MFMO) method applied to MD cluster structures of bound and free states each surrounded by equal-sized explicit water shells.
What would settle it
A direct experimental measurement of the absolute binding free energy or a higher-level ab initio calculation performed on the identical cluster structures that differs substantially from -67.2 kcal/mol.
Extended reading notes
Core claim
Using FMO2-RHF/STO-3G:MP2/6-31G(d) on MD-derived clusters with equal-sized water shells around bound and free states, the weighted binding energy between 17β-estradiol and the ligand-binding domain of human estrogen receptor alpha equals -67.2 kcal/mol, arising principally from electrostatic contacts at Glu 353, His 524 and Arg 394 plus dispersion contacts at hydrophobic sites.
Load-bearing premise
The chosen cluster structures drawn from the molecular dynamics trajectories represent the thermodynamic ensemble of bound and free states without large systematic energy errors.
Editorial extensions
If this is right
- Dispersion forces contribute attractively between 17β-estradiol and hydrophobic binding-site residues.
- Strong electrostatic interactions occur between 17β-estradiol and the charged or polarized residues Glu 353, His 524 and Arg 394.
- The same model-system construction can be reused to rank binding affinities of other candidate endocrine-disrupting chemicals.
- Solvation is treated explicitly by placing identical numbers of water molecules around each state before subtracting the energies.
Reading between the lines
- - Extending the same matched-water-shell protocol to additional ligands could generate a ranked list of binding strengths for screening purposes.
- - Discrepancies between this value and measured dissociation constants would point to the size of the remaining entropic or dynamic contributions omitted by the static-cluster approach.
- - The residue-level interaction map could guide targeted mutagenesis experiments that test which contacts dominate the measured affinity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports an approximate calculation of the binding energy between 17β-estradiol and human estrogen receptor alpha via multilayer fragment molecular orbital (MFMO/FMO2) calculations at the RHF/STO-3G:MP2/6-31G(d) level. Structures are taken from MD trajectories via cluster analysis; the model subtracts total energies of the bound complex, free receptor, free ligand, and a reference water cluster, each solvated by equal-sized water shells. The reported weighted binding energy is -67.2 kcal/mol, accompanied by qualitative observations of dispersion contacts with hydrophobic residues and electrostatic contacts with Glu353, His524, and Arg394.
Significance. If the central numerical result holds after validation, the work supplies a concrete, parameter-free estimate obtained from a fragment-based quantum method on an explicitly solvated biomolecular system. This could serve as a template for similar EDC-binding calculations. The approach avoids empirical fitting to the target binding value and uses a consistent water-shell cancellation scheme, both of which are methodological strengths.
major comments (2)
- [Abstract (model-systems paragraph)] Abstract (model-systems paragraph): the binding energy is obtained directly as a difference of four total energies evaluated on cluster-analysis snapshots; no convergence diagnostics are supplied for cluster count, MD trajectory length, or reweighting, nor is any comparison shown between the selected structures and a larger ensemble or an alternative sampling protocol. Because the reported -67.2 kcal/mol is a single weighted difference, systematic bias in the chosen snapshots directly affects the central numerical claim.
- [Abstract] Abstract: the final result is stated as a single scalar (-67.2 kcal/mol) with neither reported standard error from the weighting procedure nor any comparison to experimental binding free energies. This omission leaves the precision and physical plausibility of the difference unquantified.
minor comments (1)
- [Abstract] The abstract refers to a 'weighted binding energy' without stating the weighting scheme, the number of clusters retained, or the criterion used to select them.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on sampling validation and uncertainty quantification. We respond to each major comment below.
read point-by-point responses
-
Referee: [Abstract (model-systems paragraph)] Abstract (model-systems paragraph): the binding energy is obtained directly as a difference of four total energies evaluated on cluster-analysis snapshots; no convergence diagnostics are supplied for cluster count, MD trajectory length, or reweighting, nor is any comparison shown between the selected structures and a larger ensemble or an alternative sampling protocol. Because the reported -67.2 kcal/mol is a single weighted difference, systematic bias in the chosen snapshots directly affects the central numerical claim.
Authors: The referee is correct that the original manuscript supplies no explicit convergence diagnostics or comparisons against alternative sampling protocols. The structures were obtained via standard cluster analysis of the MD trajectory, with the reported value computed as a population-weighted average. We will revise the manuscript to include the number of clusters retained, their populations, and a short statement on the clustering protocol used. revision: yes
-
Referee: [Abstract] Abstract: the final result is stated as a single scalar (-67.2 kcal/mol) with neither reported standard error from the weighting procedure nor any comparison to experimental binding free energies. This omission leaves the precision and physical plausibility of the difference unquantified.
Authors: We agree that the absence of a reported uncertainty measure and the lack of any experimental comparison leave the result incompletely contextualized. The weighting derives directly from the cluster populations, but no standard error was computed. Moreover, the quantity we calculate is an approximate electronic binding energy in an explicit-solvent model and is not equivalent to an experimental binding free energy. In revision we will add a clarifying paragraph on these distinctions and the associated limitations. revision: partial
Circularity Check
No circularity: binding energy is direct difference of four computed total energies
full rationale
The paper obtains the reported -67.2 kcal/mol value strictly as a weighted difference of FMO2-RHF/STO-3G:MP2/6-31G(d) total energies for the bound complex, free receptor, ligand, and reference water cluster (each with equal-sized water shells). Structures are taken from MD + cluster analysis, but the energy difference itself contains no fitted parameter, no self-citation load-bearing step, and no reduction of the output to the input by construction. The central numerical claim is therefore independent of the target value and does not match any of the enumerated circularity patterns.
Assumptions & free parameters
assumptions (2)
- domain assumption Fragment molecular orbital method with the stated layering accurately recovers the total energy of the full system
- domain assumption Molecular dynamics force field and cluster analysis produce representative bound and unbound conformations
Cite this review
Pith. "Pith review of Approximate calculation of the binding energy between 17$\beta$-estradiol and human estrogen receptor alpha." pith.science (2026). https://pith.science/paper/ZOBOYZCU
@misc{pith2026190710808,
author = {Pith},
title = {Pith review of: Approximate calculation of the binding energy between 17$\beta$-estradiol and human estrogen receptor alpha},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZOBOYZCU}},
note = {Machine review of arXiv:1907.10808}
}
abstract
Estrogen receptors (ERs) are a group of proteins activated by 17$\beta$-estradiol. The endocrine-disrupting chemicals (EDCs) mimic estrogen action by bind directly to the ligand binding domain of ER. From this perspective, ER represent a good model for identifying and assessing the health risk of potential EDCs. This ability is best reflected by the ligand-ER binding energy. Multilayer fragment molecular orbital (MFMO) calculations were performed which allowed us to obtain the binding energy using a calculation scheme that considers the molecular interactions that occur on the following model systems: the bound and free receptor, 17$\beta$-estradiol and a water cluster. The bound and free receptor and 17$\beta$-estradiol were surrounded by a water shell containing the same number of molecules as the water cluster. The structures required for MFMO calculations were obtained from molecular dynamics simulations and cluster analysis. Attractive dispersion interactions were observed between 17$\beta$-estradiol and the binding site hydrophobic residues. In addition, strong electrostatic interactions were found between 17$\beta$-estradiol and the following charged/polarized residues: Glu 353, His 524 and Arg 394. The FMO2-RHF/STO-3G:MP2/6-31G(d) weighted binding energy was of -67.2 kcal/mol. We hope that the model developed in this study can be useful for identifying and assessing the health risk of potential EDCs.
Reviewed May 24, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.