A General Theory for Phenotypic Association in Biological Systems

Giuseppe Battaglia

REVIEW 3 major objections 5 minor 107 references

A General Theory for Phenotypic Association in Biological Systems

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper proves that multivalent binding—the way antibodies, lipoproteins, and T-cell receptors engage surfaces through many weak bonds at once—is exactly the monomer-dimer problem on a bipartite reachability graph, and therefore that…

desk verdict Solid core theorem and an exact three-body term, but the 'never a phase transition' headline is proven only for J=0 and overreaches in the general functional. read the letter →

arxiv 2608.08571 v1 pith:YRR3BP24 submitted 2026-08-09 physics.bio-ph q-bio.MN

Giuseppe Battaglia This is my paper · ORCID

classification physics.bio-phq-bio.MN

keywords multivalencymonomer-dimerproblemmatchingpolynomialsuperselectivitypotentialofmeanforcemany-bodyforceskineticproofreadingbiologicalrecognition

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the recognition between two multivalent objects in water—objects that present several binding groups on a shared scaffold—is governed by a single counting problem: the number of distinct ways to bind equals the number of matchings of a bipartite reachability graph. Identifying the contact partition function with the matching polynomial of that graph imports a classical theorem from statistical mechanics, which forces two structural conclusions: the bound fraction is an analytic function of receptor density, so superselective thresholds are smooth crossovers and never genuine phase transitions, and the equilibrium distribution of the number of simultaneous bonds is always unimodal. From this, three design constraints are derived: a repulsive surface layer is thermodynamically obligatory, bond free energies cannot be treated as independent, and because rates are exponential in free energy, small changes in receptor number can shift binding lifetimes by many orders of magnitude. The same equations then unify antibody avidity, lipoprotein clearance, and T-cell antigen recognition, three systems usually described in separate vocabularies. A sympathetic reader would describe the paper as an attempt to place the whole phenomenology of multivalent recognition on one Hamiltonian, with a rigorous no-phase-transition ceiling.

What carries the argument

The central object is the topology kernel $\Omega(k)=m_k(G)$, the number of $k$-matchings of the bipartite reachability graph; it enters the free energy logarithmically and is the only route through which architecture (ligand spacing, spacer reach, scaffold rigidity, steric blocking) affects thermodynamics. The carrying theorem is the Heilmann–Lieb real-rootedness of the monomer-dimer partition function, which yields analyticity of the bound fraction and log-concavity $\Omega(k)^2\ge\Omega(k-1)\Omega(k+1)$, hence unimodality. The kernel is evaluated exactly for five reference architectures—flexible ligands with mobile receptors, rigid commensurate scaffolds, rigid incommensurate scaffolds, neighbour exclusion, and tethered ligands against dilute receptors—each a different edge set of the same graph. The many-body side is carried by the body-order expansion, with crowder depletion solved exactly: the irreducible three-body term $w^{(3)}=\Pi V_{123}$ is repulsive, so pairwise treatments overestimate depletion attraction by 17–24% at contact.

What would settle it

Track the equilibrium distribution of the number of simultaneous bonds between a single well-defined multivalent unit (for example a DNA-origami construct of valency 8–12) and a receptor-presenting surface while receptor density is swept over orders of magnitude; the theory predicts the distribution is always single-peaked and the bound fraction is smooth. Observing two coexisting bond-number peaks at any density, or a discontinuity in the binding isotherm, would falsify the claim and reveal physics the static-reachability model omits.

Watch

Extended reading notes

Core claim

The central claim is that the equilibrium of two multivalent units is exactly the monomer-dimer problem on the reachability graph $G=(L\cup R,E)$, where an edge $(i,j)$ exists precisely when ligand $i$ can physically reach receptor $j$. The contact partition function $\Xi(x)=\sum_k \Omega(k)x^k$ is the matching generating polynomial of $G$, with $\Omega(k)=m_k(G)$ the number of $k$-matchings. Because the Heilmann–Lieb theorem guarantees that this polynomial has only real, non-positive roots for every graph and every set of non-negative activities, the bound fraction is an analytic, strictly monotonic function of receptor density with no singularity; a superselective threshold is a sharp crossover, never a phase transition. The same theorem implies the coefficient sequence is log-concave, so the equilibrium bond-number distribution is unimodal—there is always a single most-probable engagement $k^\star$, never two competing ones. The paper then assembles the general functional $W_{AB}(h)=W_{\mathrm{rep}}(h)-k_BT\ln\sum_k \Omega(k;g,\xi)\exp[-\beta(\sum_\zeta k_\zeta\varepsilon_\zeta(h)+\sum_{\zeta\le\zeta'}J_{\zeta\zeta'}k_\zeta k_{\zeta'})]$, whose limits recover classical monovalent affinity, the standard combinatorial-entropy theory of multivalency, and a many-body regime where engagement can self-limit.

Load-bearing premise

The theory assumes a binding encounter is fully described by a static reachability graph whose edges say whether a ligand can physically reach a receptor, with all thermodynamics and kinetics flowing from the matching counts of that graph; if real encounters involve time-dependent conformations, rebinding, induced fit, or interactions that cannot be represented as a binary edge, the monomer-dimer mapping and its no-phase-transition conclusion do not apply to the real system.

Editorial extensions

If this is right

Any multivalent construct needs a repulsive steric layer; without one, surfaces in physiological salt fall into a primary minimum about $21\,k_BT$ deep and bind irreversibly.
Bond additivity is quantitatively wrong: scaffold connectivity, corona compression and crowder depletion generate irreducible many-body couplings that can either sharpen or blunt a threshold, and pairwise models systematically bias computed binding constants.
Kinetic selectivity can far exceed equilibrium selectivity, because lifetimes grow exponentially with well depth: in the paper's worked example, a tenfold increase in receptor number moves a contact's residence time from roughly 10 ms to months.
The antibody isotype hierarchy, the ApoE2 lipoprotein clearance defect, and the T-cell memory threshold are one physics: how a fixed total avidity is distributed over bonds, with antigen size gating usable valency and kinetic proofreading manufacturing specificity beyond the equilibrium ceiling.
Discrimination beyond what free energies allow is possible only by breaking detailed balance: T-cell antigen recognition is modelled as kinetic proofreading, whose specificity grows as the $N$-th power of the dwell-time ratio while sensitivity falls.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The matching-polynomial identification implies a strong equivalence principle the paper leaves implicit: two chemically different constructs whose reachability graphs have the same matching polynomial are thermodynamically indistinguishable at equilibrium, so 'matching-polynomial equivalence' could serve as a design criterion for multivalent therapeutics.
Because the theory states that architecture changes only the per-bond affinity required, not the achievable sharpness, it predicts that sharply different scaffolds (DNA origami, flexible polymers, colloids) should display the same maximal selectivity exponent once per-bond free energies are matched—a quantitative test that existing superselectivity data could be re-analysed to check.
The ageing section's saddle-node bifurcation suggests a general experimental signature: physiological decline should be preceded by critical slowing down in the kinetics of multivalent contacts, making lifetime fluctuations a candidate early-warning observable for tissue dysfunction.
The valency-amplification principle—per-bond defects are multiplied by valency—generalises beyond ApoE2: any high-valency interface, engineered or evolved, is a point of fragility where small mutations produce outsized functional effects, which could inform the interpretation of disease-associated variants at multivalent binding sites.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

Read Battaglia's 'General Theory for Phenotypic Association'. The paper is better than the title suggests, and the title is the problem.

What is actually new: a clean transposition of the Heilmann-Lieb theorem to multivalent binding, showing that for the standard superselectivity partition function Ξ(x)=Σ Ω(k)x^k the bond-number distribution is log-concave and there is no phase transition; an exact three-body depletion term (eq. 29) with a surprising sign and a sharp geometric onset at q*≈0.155; and a four-architecture comparison showing topology costs per-bond affinity rather than sharpness. These are real contributions, clearly written, with the limitations mostly acknowledged.

The soft spot is the one the abstract hides. The no-phase-transition and log-concavity results are proven for the additive J=0 kernel. But Section 8 argues additivity is generically false, and the general functional (eq. 32) includes a many-body coupling J entering as Jk². That is not a monomer-dimer model; Heilmann-Lieb does not apply. For cooperative J<0 the effective weights Ω(k)e^{β|J|k²} can violate their own inequality (18) and produce a bimodal bond-number distribution. The stress-test example (ℓ=2, weights 1,20,10000) is correct. So the abstract's unconditional 'never a phase transition' and 'always a single most probable engagement' are not consequences of the general functional. This is a load-bearing overclaim, because Part I advertises it as the general result.

Also, J itself is not derived. Table 5 catalogs mechanisms and signs, but only depletion is computed exactly. The biological applications (isotypes, ApoE, T cells) are illustrative; the authors say so honestly, but 'general theory' overpromises. Minor: reference [39] is incomplete and there are typos.

The core theorem is correct and worth having, and the exact three-body term is a genuine calculation. With serious reviewing the path is clear: qualify the phase-transition and unimodality claims to the J=0 case, or supply a new proof for J>0, and present the biology as worked examples rather than validation.

I would send it to peer review with heavy revision requested. It deserves the referee time.

Referee Report

3 major / 5 minor

Summary. This manuscript proposes a statistical-mechanical theory of multivalent biological recognition. It defines a multivalent unit as a scaffold carrying binding moieties and a repulsive corona, represents a binding encounter as a matching problem on a bipartite reachability graph, and identifies the topology kernel Ω(k) with the number of k-matchings of that graph. The contact partition function is then the matching polynomial, so Heilmann-Lieb real-rootedness is invoked to conclude that multivalent binding is always a smooth crossover, never a phase transition, and that the bond-number distribution is log-concave and unimodal. The paper argues that additivity of bond free energies generically fails, computes an exact three-body depletion correction, introduces a many-body coupling matrix J, and writes a general functional (Eq. 32). Part II adds dynamics via a master equation and Smoluchowski equation, Kramers lifetime amplification, kinetic proofreading, and ageing as a slow drift ending in a saddle-node bifurcation. Part III applies the framework to immunoglobulins, lipoproteins, and T cell recognition.

Significance. The monomer-dimer identification and its additive-limit consequences are rigorous and valuable: real-rootedness gives a parameter-free bound on threshold sharpness, and the exact three-body depletion calculation is a concrete, checkable result with a clear physical mechanism. The paper synthesizes a large literature and makes falsifiable quantitative claims in the additive and anti-cooperative regimes. Its breadth across antibodies, lipoproteins, and T cells is genuinely stimulating. However, the central advertised conclusion is asserted for the general functional (Eq. 32), while the proof covers only the J=0 additive limit that Section 8 itself identifies as generically false; for cooperative J<0 the log-concavity and no-phase-transition conclusions are not established and can fail. The paper therefore needs substantive qualification or additional proof before its headline claims can be accepted.

major comments (3)

[§6.3, Eqs. (17), (31), (32)] The headline claim—'never a phase transition' and 'always a single most probable engagement'—is proven only for the matching partition function Ξ(x)=Σ_k Ω(k)x^k with x=e^{-βε}>0. The general functional (32) replaces x^k by exp[-β(εk+Σ_{ζ≤ζ'}J_{ζζ'}k_ζk_ζ')]; for a single receptor type this is Σ_k Ω(k)e^{-β(εk+Jk^2)}, which is not a monomer-dimer partition function, and Heilmann-Lieb does not apply. For cooperative J<0 (Table 5, Fig. 4), the coefficient sequence Ω(k)e^{β|J|k^2} can violate inequality (18): with ℓ=2, Ω=(1,2,1), and βJ=-ln10, the weights are (1,20,10^4), and 20^2 < 1·10^4. Thus the unconditional statements in the abstract and Section 6.3 are not consequences of the general functional; they hold only in the J=0 additive limit that Section 8 itself identifies as generically false. The abstract and Sections 6.3 and 9 must be qualified, or a new proof supplied for J≠0.
[§9 and Fig. 4] Figure 4 and the surrounding text compute θ and k⋆ for the 'synergistic' case J=-0.15 kBT and assert that attractive coupling sharpens the threshold beyond the additive result. Since for J<0 the log-concavity guarantee of Eq. (18) is lost, the existence and uniqueness of k⋆ are not established for this case; the caption's statement that k⋆ is 'the most probable engagement' presumes unimodality, which is exactly what J can destroy. Please either prove log-concavity/unimodality for J<0 or restrict the quantitative claims in this section to J=0 and J>0, where the product of log-concave sequences preserves the property.
[§6.1, Eq. (12), and §21] The entire no-phase-transition and log-concavity result rests on the assumption that a binding encounter is represented by a static bipartite reachability graph G whose edge set is fixed by geometry. The manuscript does not test whether time-dependent conformations, rebinding, induced fit, or receptor redistribution can be absorbed into a static G; if they cannot, the theorem does not apply to the advertised biological systems (antibodies, lipoproteins, T cells). Section 21 acknowledges that data are missing, but the abstract's scope claim ('One set of equations then covers antibodies, lipoproteins, and T cell recognition') goes beyond what is demonstrated. Please state explicitly the modeling assumptions under which the theorem applies and, ideally, test the framework against at least one experimental dataset.

minor comments (5)

[§3.2, Eq. (5)] The integral notation in Eq. (5) is garbled: the numerator should integrate over the polar angle θ from 0 to φ0 (with dθ), and the denominator over the full sphere; as printed, 'sin θ dϕ' is dimensionally inconsistent.
[§4 and §7] There are typographical errors: 'capatable' in Section 4 should be 'capable', and 'convenintely' in Section 7 should be 'conveniently'.
[References] Reference [39] is incomplete: it ends with '[citation to be completed]' and cannot be verified. Please complete this reference before resubmission.
[Appendix A and Fig. 5] The Kramers prefactor τ0=1 ns is an estimate and shifts all absolute residence times in Fig. 5 uniformly; this caveat is stated only in the appendix and should appear with Fig. 5 so readers do not take the absolute milliseconds-to-months values literally.
[§8.2] The text states that the three-body depletion expression is validated against Monte Carlo and implemented in released code, but no code URL or Monte Carlo details are provided; please include them or give an explicit reference.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the matching-polynomial theorem is external and parameter-free; the J-coupling gap is a correctness caveat, not a circular reduction.

full rationale

The central derivation is self-contained: Eq. 12 identifies the topology kernel with matching counts, Eq. 17 is the matching-generating polynomial, and the real-rootedness/log-concavity statement is imported from the external Heilmann-Lieb theorem [35]. The paper does not fit parameters and then re-predict the same fitted quantities; the three-body depletion term (Eqs. 28-29) is computed exactly from geometry and validated by Monte Carlo, and the immunoglobulin, lipoprotein, and T-cell numbers follow from stated valencies, literature affinities, and the stated model. Self-citations [25,26] support the existence of density-reading carriers but are not load-bearing for the theorem. There is a genuine overreach: Section 6.3's Eq. 18 guarantee also backs Section 12's uniqueness of k*, yet the general functional of Eq. 32 includes the J coupling; sufficiently cooperative J<0 can violate Eq. 18. That is a correctness/scope caveat, not a circular step. Reference [39] is a placeholder citation and is missing support, but again not circular. No derived prediction is definitionally identical to an input or to a fitted quantity, so the circularity score is 0.

Assumptions & free parameters 6 free parameters · 5 assumptions · 0 invented entities

The central theorem depends only on the standard Heilmann-Lieb result and the static graph model. The broader theory adds hand-set parameters (epsilon, J, Wrep, tau0, kp, delta) and assumes a quadratic many-body truncation. No new particles or forces are postulated. The ledger shows the framework is not fully self-contained: the many-body couplings are asserted to have microscopic origins but are not computed for the biological systems.

free parameters (6)

per-bond free energy epsilon = -22.8, -7.9, -8, -9.7 kBT in fig 4; -4 kBT in fig 5; -4.5 to -7.5 kBT in fig 8
Chosen so curves cross a common threshold in fig 4 or to represent isoform differences in fig 8; central to the shape of the predicted curves.
many-body coupling J = 0, -0.15, +0.60 kBT in fig 4; 1 kBT in fig 5
Hand-set; no independent measurement or derivation for the systems treated. Controls self-limiting engagement and threshold sharpness.
corona toll Wrep = 9 kBT in fig 4; 12 kBT in fig 5
Chosen to illustrate thresholds; sets the price of repulsion.
Kramers attempt time tau0 = 1 ns
Estimate from diffusivity; the paper acknowledges it shifts all residence times by a common factor.
proofreading rate kp = 0.1 s^-1 in fig 9
Representative value for T cell signalling; not derived.
contact layer thickness delta = 2 nm in fig 8
Defines receptor availability area Ac = 2 pi R delta for lipoprotein calculations; not justified.

assumptions (5)

standard math Heilmann-Lieb theorem on monomer-dimer systems (real roots, no phase transition)
Invoked in section 6.3 to conclude smooth crossover and log-concavity; treated as an external rigorous result.
domain assumption Potential of mean force (eq 1) is a valid free energy and fast coordinates equilibrate quickly
Section 2.2 assumes the solvent is always at equilibrium with respect to the units, based on timescale separation from table 1.
domain assumption Reachability graph is a binary relation determined by geometry (eq 14)
Section 6.1 models physical reach as a yes/no edge; all thermodynamics and kinetics flow from matching counts of this graph.
ad hoc to paper All non-additive interactions captured by a quadratic coupling matrix J (eq 31)
Section 8.4 truncates the body-order expansion at second order with no estimate of higher-order terms; J is posited rather than derived.
domain assumption Biological parameters in Table 6 (valencies, per-bond energies) are taken from literature as inputs
Part III uses reported Fab valencies and receptor affinities (e.g., Bruhns 2009) to compute thresholds; these are not derived in the paper.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A General Theory for Phenotypic Association in Biological Systems." pith.science (2026). https://pith.science/paper/YRR3BP24

@misc{pith2026260808571,
  author       = {Pith},
  title        = {Pith review of: A General Theory for Phenotypic Association in Biological Systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YRR3BP24}},
  note         = {Machine review of arXiv:2608.08571}
}

read the original abstract

Biological recognition rarely rests on one strong bond. It works by forming many weak ones at once, between crowded, deformable surfaces in water. This review develops that process as a problem in statistical mechanics. Counting the ways two multivalent objects can bind proves to be the classical monomer-dimer problem on a graph, with a rigorous consequence: the apparent switching of multivalent binding is always a smooth crossover, never a phase transition. Three constraints follow. A repulsive surface layer is obligatory rather than a design choice; bonds do not act independently; and since free energies enter rates exponentially, small changes in receptor number shift binding lifetimes by orders of magnitude. One set of equations then covers antibodies, lipoproteins, and T cell recognition. In each, what decides the outcome is not the strength of any single bond but how a fixed total is spread over many: affinity is a property of a molecule, selectivity a property of an assembly.

Figures

Figures reproduced from arXiv: 2608.08571 by the authors.

**Figure 1.** The elementary contributions to the pair potential of mean force in water, all in units of kBT. (a) The specific, directional terms are contact-ranged; the non-specific ones are weak or area-driven. (b) Repulsive terms: a hard core, a double layer screened to 0.79 nm, hydration oscillating with the width of a water molecule, and a polymer corona that dominates out to twice its height. (c) The net PMF between two uni… view at source ↗

**Figure 2.** The topology kernel is a k-matching count of the bipartite reachability graph, eq. (12). (a) Reachability graphs for four architectures with ℓ = 8 ligands (circles, on the scaffold) and R = 14 receptors (squares); an edge means the ligand can reach the receptor. (b) The resulting kernels, computed exactly, spanning six orders of magnitude at fixed valency. (c) Selectivity exponent α when the per-bond free energy is … view at source ↗

**Figure 3.** Pairwise additivity failing, exactly. (a) Three units in crowders with size ratio q = 1: the pairwise sum overestimates the attraction relative to the exact three-body result, the difference being the irreducible term w (3) = +ΠV123 of eq. (29), which is repulsive. (b) Nonadditivity at contact versus crowder size, showing the geometric onset at q ∗ ≃ 0.155 of eq. (30) and growth to about 24% [PITH_FULL_IMAGE:figur… view at source ↗

Figures from the paper (6 more)

**Figure 4.** Figure 4: The three limits of the general functional, eq. (32). (a–c) Schematics: a monovalent unit, for which the topology kernel is trivial and the coupling matrix has nothing to couple; an additive multivalent unit, whose bonds ignore one another; and a many-body unit, in whi…

**Figure 5.** Figure 5: Dynamics of a multivalent contact, computed from eq. (31) with ℓ = 12, ε = −4 kBT, corona toll Wrep = 12 kBT and self-coupling J = 1 kBT. (a) The engagement landscape. Approach must first pay the corona toll at k = 0, which acts as the transition state; the bound minim…

**Figure 6.** Figure 6: (a) The antigen face of immunoglobulins, from eq. (21), across the full isotype valency series: IgG monomer (ℓ = 2), secretory IgA dimer (ℓ = 4) and tetramer (ℓ = 8), IgM pentamer (ℓ = 10) and hexamer (ℓ = 12). Valencies count Fab arms per assembled molecule; IgE and I…

**Figure 7.** Figure 7: The effector face. Engagement probability against the number of presented Fc stems, computed from eq. (21) with monomeric affinities from Bruhns et al. [29]. High-affinity FcεRI engages at a single stem; low-affinity FcγRIIIa requires tens, and therefore an immune comp…

**Figure 8.** Figure 8: Lipoproteins as multivalent units, computed from eq. (31) with receptor availability set by the contact area Ac = 2πRδ and δ = 2 nm. (a) Avidity against receptor density for four classes. Multivalent ApoE-bearing particles reach far higher avidity than the single-conta…

**Figure 9.** Figure 9: T cell recognition as a driven, non-equilibrium process. (a) Kinetic proofreading, eq. (40), with kp = 0.1 s −1 . Increasing the number of proofreading steps sharpens discrimination between a short-lived self complex and a long-lived agonist complex, at the cost of ove…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

107 extracted references · 79 canonical work pages

[39]

[citation to be completed]

Machine-learned coarse-grained many-body potentials for ligand-stabilised nanoparticles (2023). [citation to be completed]

work page 2023
[1]

J. N. Israelachvili, Intermolecular and Surface Forces , 3rd ed., Academic Press (2011)

2011
[2]

P. I. Kitov and D. R. Bundle, On the nature of the multivalency eﬀect: a thermodynamic model , J. Am. Chem. Soc. 125, 16271 (2003). 35

2003
[3]

F. J. Martinez-Veracoechea and D. Frenkel, Designing super selectivity in multivalent nano-particle binding , Proc. Natl. Acad. Sci. USA 108, 10963 (2011)

2011
[4]

Asakura and F

S. Asakura and F. Oosawa, On interaction between two bodies immersed in a solution of macromolecules, J. Chem. Phys. 22, 1255 (1954)

1954
[5]

B. M. Axilrod and E. Teller, Interaction of the van der Waals type between three atoms , J. Chem. Phys. 11, 299 (1943)

1943
[6]

C. P. Brangwynne et al. , Germline P granules are liquid droplets that localize by controlled dissolution/condensation, Science 324, 1729 (2009)

2009
[7]

T. J. Nott et al. , Phase transition of a disordered nuage protein generates environmentally responsive membraneless organelles , Mol. Cell 57, 936 (2015)

2015

Show all 107 references

[8]

Nicolis and I

G. Nicolis and I. Prigogine, Self-Organization in Nonequilibrium Systems , Wiley, New York (1977)

1977
[9]

J. L. England, Dissipative adaptation in driven self-assembly , Nat. Nanotechnol. 10, 919 (2015)

2015
[10]

P. Bak, C. Tang and K. Wiesenfeld, Self-organized criticality: an explanation of the 1/ f noise, Phys. Rev. Lett. 59, 381 (1987)

1987
[11]

E. N. Lorenz, Deterministic nonperiodic ﬂow , J. Atmos. Sci. 20, 130 (1963)

1963
[12]

S. H. Strogatz, Nonlinear Dynamics and Chaos , 2nd ed., Westview Press (2015)

2015
[13]

Goulian, R

M. Goulian, R. Bruinsma and P. Pincus, Long-range forces in heterogeneous ﬂuid membranes , Europhys. Lett. 22, 145 (1993)

1993
[14]

Azadbakht, T

A. Azadbakht, T. R. Weikl and D. J. Kraft, Many-body eﬀects in membrane-mediated interactions between colloidal particles , arXiv:2403.05253 (2024)

2024 arXiv
[15]

J. C. Pérez-Ángel et al. , Quantifying many-body contributions to depletion forces in hard-sphere mixtures, arXiv:2509.03342 (2025)

2025 arXiv
[16]

Chandler, Interfaces and the driving force of hydrophobic assembly , Nature 437, 640 (2005)

D. Chandler, Interfaces and the driving force of hydrophobic assembly , Nature 437, 640 (2005)

2005
[17]

H.-X. Zhou, G. Rivas and A. P. Minton, Macromolecular crowding and conﬁnement , Annu. Rev. Biophys. 37, 375 (2008)

2008
[18]

Alexander, Adsorption of chain molecules with a polar head: a scaling description , J

S. Alexander, Adsorption of chain molecules with a polar head: a scaling description , J. Phys. (Paris) 38, 983 (1977)

1977
[19]

de Gennes, Polymers at an interface: a simpliﬁed view , Adv

P.-G. de Gennes, Polymers at an interface: a simpliﬁed view , Adv. Colloid Interface Sci. 27, 189 (1987)

1987
[20]

S. T. Milner, T. A. Witten and M. E. Cates, Theory of the grafted polymer brush , Macromolecules 21, 2610 (1988)

1988
[21]

Helfrich, Elastic properties of lipid bilayers: theory and possible experiments , Z

W. Helfrich, Elastic properties of lipid bilayers: theory and possible experiments , Z. Naturforsch. C 28, 693 (1973)

1973
[22]

T. Curk, S. Angioletti-Uberti, J. Dobnikar, D. Frenkel et al. , Receptor clustering tunes and sharpens the selectivity of multivalent binding , Proc. Natl. Acad. Sci. USA 122, e2417159122 (2025). 36

2025
[23]

G. V. Dubacheva, T. Curk and R. P. Richter, Determinants of superselectivity: practical concepts for application in biology and medicine , Acc. Chem. Res. 56, 729 (2023)

2023
[24]

Ravnik, M

J. Ravnik, M. Chabaud, U. Bren, G. V. Dubacheva and T. Curk, Kinetic superselectivity in multivalent binding , arXiv:2605.27019 (2026)

2026 arXiv
[25]

X. Tian, S. Angioletti-Uberti and G. Battaglia, On the design of precision nanomedicines , Sci. Adv. 6, eaat0919 (2020)

2020
[26]

Liu et al

Q. Liu et al. , Combinatorial entropy behaviour leads to range selective binding in ligand-receptor interactions, Nat. Commun. 11, 4836 (2020)

2020
[27]

Linne, D

C. Linne, D. Visco, S. Angioletti-Uberti, L. Laan and D. J. Kraft, Direct visualization of superselective colloid-surface binding mediated by multivalent interactions , Proc. Natl. Acad. Sci. USA 118, e2106036118 (2021)

2021
[28]

H. Bila, K. Paloja, V. Caroprese, A. Kononenko and M. M. C. Bastings, Multivalent pattern recognition through control of nano-spacing in low-valency super-selective materials , J. Am. Chem. Soc. 144, 21576 (2022)

2022
[29]

Bruhns et al

P. Bruhns et al. , Speciﬁcity and aﬃnity of human Fc γ receptors and their polymorphic variants for human IgG subclasses , Blood 113, 3716 (2009)

2009
[30]

J. J. Hopﬁeld, Kinetic proofreading: a new mechanism for reducing errors in biosynthetic processes requiring high speciﬁcity , Proc. Natl. Acad. Sci. USA 71, 4135 (1974)

1974
[31]

Ninio, Kinetic ampliﬁcation of enzyme discrimination , Biochimie 57, 587 (1975)

J. Ninio, Kinetic ampliﬁcation of enzyme discrimination , Biochimie 57, 587 (1975)

1975
[32]

T. W. McKeithan, Kinetic proofreading in T-cell receptor signal transduction , Proc. Natl. Acad. Sci. USA 92, 5042 (1995)

1995
[33]

L. G. Valiant, The complexity of computing the permanent , Theor. Comput. Sci. 8, 189 (1979)

1979
[34]

P. W. Kasteleyn, Graph theory and crystal physics , in Graph Theory and Theoretical Physics , ed. F. Harary, Academic Press, London, pp. 43–110 (1967)

1967
[35]

O. J. Heilmann and E. H. Lieb, Theory of monomer-dimer systems , Commun. Math. Phys. 25, 190 (1972)

1972
[36]

Mézard and A

M. Mézard and A. Montanari, Information, Physics, and Computation , Oxford University Press, Oxford (2009)

2009
[37]

Battiston et al

F. Battiston et al. , The physics of higher-order interactions in complex systems , Nat. Phys. 17, 1093 (2021)

2021
[38]

Drautz, Atomic cluster expansion for accurate and transferable interatomic potentials , Phys

R. Drautz, Atomic cluster expansion for accurate and transferable interatomic potentials , Phys. Rev. B 99, 014104 (2019)

2019
[40]

S. B. Zimmerman and S. O. Trach, Estimation of macromolecule concentrations and excluded volume eﬀects for the cytoplasm of Escherichia coli , J. Mol. Biol. 222, 599 (1991)

1991
[41]

R. J. Ellis, Macromolecular crowding: obvious but underappreciated , Trends Biochem. Sci. 26, 597 (2001). 37

2001
[42]

O. G. Berg, R. B. Winter and P. H. von Hippel, Diﬀusion-driven mechanisms of protein translocation on nucleic acids , Biochemistry 20, 6929 (1981)

1981
[43]

S. F. Banani, H. O. Lee, A. A. Hyman and M. K. Rosen, Biomolecular condensates: organizers of cellular biochemistry , Nat. Rev. Mol. Cell Biol. 18, 285 (2017)

2017
[44]

Chiti and C

F. Chiti and C. M. Dobson, Protein misfolding, functional amyloid, and human disease , Annu. Rev. Biochem. 75, 333 (2006)

2006
[45]

Chiti and C

F. Chiti and C. M. Dobson, Protein misfolding, amyloid formation, and human disease: a summary of progress over the last decade , Annu. Rev. Biochem. 86, 27 (2017)

2017
[46]

T. P. J. Knowles, M. Vendruscolo and C. M. Dobson, The amyloid state and its association with protein misfolding diseases , Nat. Rev. Mol. Cell Biol. 15, 384 (2014)

2014
[47]

Weinbaum, J

S. Weinbaum, J. M. Tarbell and E. R. Damiano, The structure and function of the endothelial glycocalyx layer, Annu. Rev. Biomed. Eng. 9, 121 (2007)

2007
[48]

Möckl, The emerging role of the mammalian glycocalyx in functional membrane organization and immune system regulation , Front

L. Möckl, The emerging role of the mammalian glycocalyx in functional membrane organization and immune system regulation , Front. Cell Dev. Biol. 8, 253 (2020)

2020
[49]

López-Otín, M

C. López-Otín, M. A. Blasco, L. Partridge, M. Serrano and G. Kroemer, The hallmarks of aging , Cell 153, 1194 (2013)

2013
[50]

N. T. Southall, K. A. Dill and A. D. J. Haymet, A view of the hydrophobic eﬀect , J. Phys. Chem. B 106, 521 (2002)

2002
[51]

J. N. Israelachvili and H. Wennerström, Role of hydration and water structure in biological and colloidal interactions , Nature 379, 219 (1996)

1996
[52]

Leikin, V

S. Leikin, V. A. Parsegian, D. C. Rau and R. P. Rand, Hydration forces, Annu. Rev. Phys. Chem. 44, 369 (1993)

1993
[53]

J. G. Kirkwood, Statistical mechanics of ﬂuid mixtures , J. Chem. Phys. 3, 300 (1935)

1935
[54]

Roux, The calculation of the potential of mean force using computer simulations , Comput

B. Roux, The calculation of the potential of mean force using computer simulations , Comput. Phys. Commun. 91, 275 (1995)

1995
[55]

H. C. Hamaker, The London–van der Waals attraction between spherical particles , Physica 4, 1058 (1937)

1937
[56]

V. A. Parsegian, Van der Waals Forces , Cambridge University Press (2006)

2006
[57]

Derjaguin and L

B. Derjaguin and L. Landau, Theory of the stability of strongly charged lyophobic sols , Acta Physicochim. URSS 14, 633 (1941)

1941
[58]

E. J. W. Verwey and J. T. G. Overbeek, Theory of the Stability of Lyophobic Colloids , Elsevier, Amsterdam (1948)

1948
[59]

J. V. Jokerst, T. Lobovkina, R. N. Zare and S. S. Gambhir, Nanoparticle PEGylation for imaging and therapy , Nanomedicine 6, 715 (2011)

2011
[60]

J. S. Suk, Q. Xu, N. Kim, J. Hanes and L. M. Ensign, PEGylation as a strategy for improving nanoparticle-based drug and gene delivery , Adv. Drug Deliv. Rev. 99, 28 (2016). 38

2016
[61]

M. P. Monopoli, C. Åberg, A. Salvati and K. A. Dawson, Biomolecular coronas provide the biological identity of nanosized materials , Nat. Nanotechnol. 7, 779 (2012)

2012
[62]

C. D. Walkey and W. C. W. Chan, Understanding and controlling the interaction of nanomaterials with proteins in a physiological environment , Chem. Soc. Rev. 41, 2780 (2012)

2012
[63]

Salvati et al

A. Salvati et al. , Transferrin-functionalized nanoparticles lose their targeting capabilities when a biomolecule corona adsorbs on the surface , Nat. Nanotechnol. 8, 137 (2013)

2013
[64]

P. W. K. Rothemund, Folding DNA to create nanoscale shapes and patterns , Nature 440, 297 (2006)

2006
[65]

Vrij, Polymers at interfaces and the interactions in colloidal dispersions , Pure Appl

A. Vrij, Polymers at interfaces and the interactions in colloidal dispersions , Pure Appl. Chem. 48, 471 (1976)

1976
[66]

H. N. W. Lekkerkerker and R. Tuinier, Colloids and the Depletion Interaction , Springer (2011)

2011
[67]

Muto, Force between nonpolar molecules, Proc

Y. Muto, Force between nonpolar molecules, Proc. Phys. Math. Soc. Japan 17, 629 (1943)

1943
[68]

B. J. Reynwar et al. , Aggregation and vesiculation of membrane proteins by curvature-mediated interactions, Nature 447, 461 (2007)

2007
[69]

T. R. Weikl, Membrane-mediated cooperativity of proteins, Annu. Rev. Phys. Chem. 69, 521 (2018)

2018
[70]

von Smoluchowski, Versuch einer mathematischen Theorie der Koagulationskinetik kolloider Lösungen, Z

M. von Smoluchowski, Versuch einer mathematischen Theorie der Koagulationskinetik kolloider Lösungen, Z. Phys. Chem. 92, 129 (1917)

1917
[71]

Risken, The Fokker–Planck Equation , 2nd ed., Springer (1989)

H. Risken, The Fokker–Planck Equation , 2nd ed., Springer (1989)

1989
[72]

H. A. Kramers, Brownian motion in a ﬁeld of force and the diﬀusion model of chemical reactions, Physica 7, 284 (1940)

1940
[73]

Hänggi, P

P. Hänggi, P. Talkner and M. Borkovec, Reaction-rate theory: ﬁfty years after Kramers , Rev. Mod. Phys. 62, 251 (1990)

1990
[74]

Scheﬀer et al

M. Scheﬀer et al. , Early-warning signals for critical transitions , Nature 461, 53 (2009)

2009
[75]

François, G

P. François, G. Voisinne, E. D. Siggia, G. Altan-Bonnet and M. Vergassola, Phenotypic model for early T-cell activation displaying sensitivity, speciﬁcity, and antagonism , Proc. Natl. Acad. Sci. USA 110, E888 (2013)

2013
[76]

D. M. Czajkowsky and Z. Shao, The human IgM pentamer is a mushroom-shaped molecule with a ﬂexural bias , Proc. Natl. Acad. Sci. USA 106, 14960 (2009)

2009
[77]

G. D. Victora and M. C. Nussenzweig, Germinal centers , Annu. Rev. Immunol. 30, 429 (2012)

2012
[78]

Kumar Bharathkar, B

S. Kumar Bharathkar, B. W. Parker, A. G. Malyutin, N. Haloi, K. E. Huey-Tubman, E. Tajkhorshid and B. M. Stadtmueller, The structures of secretory and dimeric immunoglobulin A, eLife 9, e56098 (2020)

2020
[79]

Kumar, C

N. Kumar, C. P. Arthur, C. Ciferri and M. L. Matsumoto, Structure of the secretory immunoglobulin A core , Science 367, 1008 (2020). 39

2020
[80]

Oskam, P

N. Oskam, P. Ooijevaar-de Heer, N. I. L. Derksen, S. Kruithof, S. W. de Taeye, G. Vidarsson, S. Reijm, T. Kissel, R. E. M. Toes and T. Rispens, At critically low antigen densities, IgM hexamers outcompete both IgM pentamers and IgG1 for human complement deposition and compleme...

2022
[81]

Boes, Role of natural and immune IgM antibodies in immune responses , Mol

M. Boes, Role of natural and immune IgM antibodies in immune responses , Mol. Immunol. 37, 1141 (2000)

2000
[82]

M. R. Ehrenstein and C. A. Notley, The importance of natural IgM: scavenger, protector and regulator, Nat. Rev. Immunol. 10, 778 (2010)

2010
[83]

Bruhns and F

P. Bruhns and F. Jönsson, Mouse and human FcR eﬀector functions , Immunol. Rev. 268, 25 (2015)

2015
[84]

Nimmerjahn and J

F. Nimmerjahn and J. V. Ravetch, Fcγ receptors as regulators of immune responses , Nat. Rev. Immunol. 8, 34 (2008)

2008
[85]

Kinet, The high-aﬃnity IgE receptor (Fc εRI): from physiology to pathology , Annu

J.-P. Kinet, The high-aﬃnity IgE receptor (Fc εRI): from physiology to pathology , Annu. Rev. Immunol. 17, 931 (1999)

1999
[86]

H. J. Gould and B. J. Sutton, IgE in allergy and asthma today , Nat. Rev. Immunol. 8, 205 (2008)

2008
[87]

C. A. Diebolder et al. , Complement is activated by IgG hexamers assembled at the cell surface , Science 343, 1260 (2014)

2014
[88]

M. S. Brown and J. L. Goldstein, A receptor-mediated pathway for cholesterol homeostasis , Science 232, 34 (1986)

1986
[89]

J. L. Goldstein and M. S. Brown, The LDL receptor , Arterioscler. Thromb. Vasc. Biol. 29, 431 (2009)

2009
[90]

J. P. Segrest et al. , Structure of apolipoprotein B-100 in low density lipoproteins , J. Lipid Res. 42, 1346 (2001)

2001
[91]

R. W. Mahley, Apolipoprotein E: cholesterol transport protein with expanding role in cell biology , Science 240, 622 (1988)

1988
[92]

R. W. Mahley and S. C. Rall, Apolipoprotein E: far more than a lipid transport protein , Annu. Rev. Genomics Hum. Genet. 1, 507 (2000)

2000
[93]

K. H. Weisgraber, T. L. Innerarity and R. W. Mahley, Abnormal lipoprotein receptor-binding activity of the human E apoprotein due to cysteine-arginine interchange at a single site , J. Biol. Chem. 257, 2518 (1982)

1982
[94]

R. W. Mahley, Y. Huang and S. C. Rall, Pathogenesis of type III hyperlipoproteinemia (dysbetalipoproteinemia), J. Lipid Res. 40, 1933 (1999)

1999
[95]

K. J. Williams, Molecular processes that handle and mishandle dietary lipids , J. Clin. Invest. 118, 3247 (2008)

2008
[96]

B. D. Evavold and P. M. Allen, Separation of IL-4 production from Th cell proliferation by an altered T cell receptor ligand , Science 252, 1308 (1991). 40

1991
[97]

M. M. Davis et al. , Ligand recognition by αβ T cell receptors , Annu. Rev. Immunol. 16, 523 (1998)

1998
[98]

Altan-Bonnet and R

G. Altan-Bonnet and R. N. Germain, Modeling T cell antigen discrimination based on feedback control of digital ERK responses , PLoS Biol. 3, e356 (2005)

2005
[99]

Aleksic et al

M. Aleksic et al. , Dependence of T cell antigen recognition on T cell receptor–peptide MHC conﬁnement time , Immunity 32, 163 (2010)

2010
[100]

C. C. Govern, M. K. Paczosa, A. K. Chakraborty and E. S. Huseby, Fast on-rates allow short dwell time ligands to activate T cells , Proc. Natl. Acad. Sci. USA 107, 8724 (2010)

2010
[101]

D. J. Irvine, M. A. Purbhoo, M. Krogsgaard and M. M. Davis, Direct observation of ligand recognition by T cells , Nature 419, 845 (2002)

2002
[102]

M. A. Purbhoo, D. J. Irvine, J. B. Huppa and M. M. Davis, T cell killing does not require the formation of a stable mature immunological synapse , Nat. Immunol. 5, 524 (2004)

2004
[103]

Grakoui et al

A. Grakoui et al. , The immunological synapse: a molecular machine controlling T cell activation , Science 285, 221 (1999)

1999
[104]

M. L. Dustin and A. S. Depoil, New insights into the T cell synapse from single molecule techniques, Nat. Rev. Immunol. 11, 672 (2011)

2011
[105]

Sallusto, D

F. Sallusto, D. Lenig, R. Förster, M. Lipp and A. Lanzavecchia, Two subsets of memory T lymphocytes with distinct homing potentials and eﬀector functions , Nature 401, 708 (1999)

1999
[106]

D. L. Farber, N. A. Yudanin and N. P. Restifo, Human memory T cells: generation, compartmentalization and homeostasis , Nat. Rev. Immunol. 14, 24 (2014)

2014
[107]

C. W. Gardiner, Stochastic Methods: A Handbook for the Natural and Social Sciences , 4th ed., Springer (2009). 41

2009

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [39]

[citation to be completed]

Machine-learned coarse-grained many-body potentials for ligand-stabilised nanoparticles (2023). [citation to be completed]

work page 2023

[2] [1]

J. N. Israelachvili, Intermolecular and Surface Forces , 3rd ed., Academic Press (2011)

2011

[3] [2]

P. I. Kitov and D. R. Bundle, On the nature of the multivalency eﬀect: a thermodynamic model , J. Am. Chem. Soc. 125, 16271 (2003). 35

2003

[4] [3]

F. J. Martinez-Veracoechea and D. Frenkel, Designing super selectivity in multivalent nano-particle binding , Proc. Natl. Acad. Sci. USA 108, 10963 (2011)

2011

[5] [4]

Asakura and F

S. Asakura and F. Oosawa, On interaction between two bodies immersed in a solution of macromolecules, J. Chem. Phys. 22, 1255 (1954)

1954

[6] [5]

B. M. Axilrod and E. Teller, Interaction of the van der Waals type between three atoms , J. Chem. Phys. 11, 299 (1943)

1943

[7] [6]

C. P. Brangwynne et al. , Germline P granules are liquid droplets that localize by controlled dissolution/condensation, Science 324, 1729 (2009)

2009

[8] [7]

T. J. Nott et al. , Phase transition of a disordered nuage protein generates environmentally responsive membraneless organelles , Mol. Cell 57, 936 (2015)

2015