pith. sign in

arxiv: 2511.16355 · v2 · submitted 2025-11-20 · ⚛️ physics.atm-clus · physics.bio-ph

Contact cluster modeling of allosteric communication in PDZ domains

Emanuel Dorbath , Fabian Rudolf , Adnan Gulzar , Gerhard Stock This is my paper

Pith reviewed 2026-05-17 21:13 UTC · model grok-4.3

classification ⚛️ physics.atm-clus physics.bio-ph
keywords allosteryPDZ domainscontact clustersmolecular dynamicsprotein dynamicsallosteric communicationnonequilibrium simulations
0
0 comments X

The pith

Dynamic contact clusters provide the modular architecture for allosteric communication in PDZ domains.

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

The paper establishes that allosteric communication in PDZ domains occurs through localized contact clusters, which are groups of highly correlated contacts between protein parts. These clusters act as modular units that change cooperatively and communicate with each other through rigid secondary structures in a multistep process. Using extensive molecular dynamics simulations on four photoswitchable PDZ domains totaling about one millisecond, the authors identify recurring clusters like loops connecting beta-sheets and link the protein's response times to individual cluster motions. A sympathetic reader would care because this offers a concrete dynamical picture of how proteins can signal between distant sites without large-scale shape changes, potentially explaining regulation in many biomolecules.

Core claim

The dynamic decomposition of PDZ domains into contact clusters uncovers a modular, dynamics-based architecture that underlies and facilitates long-range allosteric communication. The contact cluster model identifies localized groups of highly correlated contacts that mediate interactions between secondary structure elements, proposing that allostery proceeds through cooperative contact changes within clusters and communication between distant clusters transmitted through rigid secondary structures.

What carries the argument

Contact clusters, defined as localized groups of highly correlated contacts that mediate interactions between secondary structure elements and enable multistep allosteric communication.

If this is right

  • Recurring contact clusters represent shared flexible structural modules such as loops connecting beta-sheets across different PDZ domains.
  • The time scales of the nonequilibrium protein response match the motions of specific contact clusters.
  • Allostery involves multistep processes of cooperative changes within clusters and inter-cluster communication via rigid structures.
  • Variations in domains, ligands, and perturbations alter the contact clusters and their dynamics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the model holds, then targeted disruptions to specific contact clusters could be used to control allosteric pathways in PDZ-containing proteins.
  • This approach might extend to modeling allostery in other protein families by decomposing their dynamics into similar contact-based modules.
  • Experimental validation could involve correlating predicted cluster dynamics with observed relaxation times in spectroscopic measurements.

Load-bearing premise

That the contact clusters identified in the simulations physically transmit the allosteric signal instead of merely reflecting correlated motions from the overall protein dynamics.

What would settle it

Observing that mutations within identified contact clusters fail to selectively impair allosteric communication between distant sites while leaving local dynamics intact would challenge the central claim.

Figures

Figures reproduced from arXiv: 2511.16355 by Adnan Gulzar, Emanuel Dorbath, Fabian Rudolf, Gerhard Stock.

Figure 1
Figure 1. Figure 1: Structure and nonequilibrium response of a photoswitchable PDZ3 domain. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Contact clusters C1 to C7 of PDZ3WT obtained from equilibrium MD sim￾ulations. Compared are contact clusters ob￾tained from (a) shortest heavy-atom inter￾residue distances below 0.45 nm and (b) corre￾sponding Cα-distances below 0.8 nm. atoms, and require this distance to remain be￾low 0.8 nm for at least 10% of the simulation time, see Tab. S5. To compare these “local” Cα-distances with the nearest heavy-a… view at source ↗
Figure 3
Figure 3. Figure 3: Contact clusters (a) and nonequilib￾rium response (b,c) of a photoswitchable PDZ3 domain (PDZ3L6) with a ligand that is one residue longer than the PDZ3 shown in Fig.1. See the caption of Fig.1 for additional informa￾tion. photoswitching of the α3-helix, it is instructive to examine how the nonequilibrium response of PDZ3L6 differs from that of PDZ3. As an ex￾ample, Fig. 3b depicts the time evolution of th… view at source ↗
Figure 4
Figure 4. Figure 4: Contact clusters (a) and nonequi￾librium response (b,c) of a PDZ2 domain (PDZ2S) featuring a photoswitch across its binding pocket. See the caption of Fig.1 for additional information. To explore the nonequilibrium response of PDZ2S following photoswitching of its bind￾ing pocket, we again analyze the timescales of all contact distances and compute the dynam￾ical content across all clusters (see Fig. 4b,c)… view at source ↗
Figure 5
Figure 5. Figure 5: Contact clusters (a) and nonequilib￾rium response (b,c) of a PDZ2 domain (PDZ2L) featuring a photoswitched ligand across its binding pocket. See the caption of Fig.1 for additional information. SAIC analysis, we find six functional clusters as presented in Fig. 5a. While cluster C3, C5, C6, and C7 reappear similar to the other PDZ do￾mains, we in particular obtain clusters C4 and C6’, which are located in … view at source ↗
read the original abstract

Allostery, the intriguing phenomenon of long-range communication between distant sites in proteins, plays a central role in biomolecular regulation and signal transduction. While it is commonly attributed to conformational rearrangements, the underlying dynamical mechanisms remain poorly understood. The contact cluster model of allostery [J. Chem. Theory Comput. 2024, 20, 10731-10739] identifies localized groups of highly correlated contacts that mediate interactions between secondary structure elements. This framework proposes that allostery proceeds through a multistep process involving cooperative contact changes within clusters and communication between distant clusters, transmitted through rigid secondary structures. To demonstrate the validity and generality of the model, this Perspective employs extensive molecular dynamics simulations ($\sim1\,$ms total simulation time) of four different photoswitchable PDZ domains and studies how different domains, ligands, and perturbations influence both the contact clusters and their dynamical evolution. These analyses reveal several recurring clusters that represent shared flexible structural modules, such as loops connecting $\beta$-sheets, and show that the characteristic time scales of the nonequilibrium protein response can be directly associated with the motions of individual contact clusters. Thus, the dynamic decomposition of PDZ domains into contact clusters uncovers a modular, dynamics-based architecture that underlies and facilitates long-range allosteric communication.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript applies the authors' prior contact cluster model of allostery to extensive MD simulations (~1 ms total) of four photoswitchable PDZ domains. It identifies recurring contact clusters as shared flexible modules (e.g., loops connecting β-sheets) and associates their dynamical timescales with nonequilibrium response times, proposing a multistep mechanism in which cooperative contact changes within clusters and inter-cluster communication via rigid secondary structures enable long-range allostery.

Significance. If the reported associations prove robust, the work supplies a concrete, dynamics-based modular decomposition of PDZ domains that links local contact fluctuations to functional timescales, offering a testable framework that could generalize to other allosteric systems and complement conformational or energetic models.

major comments (2)
  1. [Results (nonequilibrium response and cluster timescale analysis)] The central claim that contact clusters 'mediate' long-range communication rests on timescale matching and recurrence across systems, yet no section demonstrates causality: there is no test in which contacts inside a cluster are selectively altered (while preserving fold) to measurably change allosteric coupling between distant sites, nor a control showing the cluster decomposition outperforms a null model based on secondary-structure fluctuations alone.
  2. [Methods and cluster identification] Cluster definitions and the multistep communication picture are taken from the authors' earlier J. Chem. Theory Comput. 2024 paper; the present analyses test consistency with new simulations rather than deriving the clusters from first principles independent of that framework, raising the risk that observed associations are partly by construction.
minor comments (2)
  1. [Abstract] The abstract states '~1 ms total simulation time' but provides no per-system breakdown, convergence diagnostics, or error estimates on cluster recurrence or timescale extraction; these details are needed to assess whether post-hoc definitions or incomplete sampling affect the claimed associations.
  2. [Introduction/Results] Notation for contact clusters and their time scales should be defined explicitly on first use, with clear distinction between equilibrium fluctuations and nonequilibrium response observables.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below, clarifying the evidential basis of our claims while acknowledging limitations. Revisions have been made to improve transparency and add supporting analyses.

read point-by-point responses

  1. Referee: [Results (nonequilibrium response and cluster timescale analysis)] The central claim that contact clusters 'mediate' long-range communication rests on timescale matching and recurrence across systems, yet no section demonstrates causality: there is no test in which contacts inside a cluster are selectively altered (while preserving fold) to measurably change allosteric coupling between distant sites, nor a control showing the cluster decomposition outperforms a null model based on secondary-structure fluctuations alone.

    Authors: We agree that direct causal tests via selective contact perturbation would strengthen the mediation claim. Such interventions are technically demanding in large-scale MD while preserving fold and were outside the scope of this Perspective, which instead tests the model's predictive power through recurrence and timescale associations across four independent PDZ systems. We have added a dedicated paragraph in the revised Discussion acknowledging this correlative nature of the evidence and proposing future targeted simulations (e.g., contact-restraint or mutation protocols). We have also included a new supplementary figure and text comparing cluster-derived communication pathways against a secondary-structure fluctuation null model, demonstrating that contact clusters capture additional dynamical correlations not explained by secondary structure alone. revision: partial

  2. Referee: [Methods and cluster identification] Cluster definitions and the multistep communication picture are taken from the authors' earlier J. Chem. Theory Comput. 2024 paper; the present analyses test consistency with new simulations rather than deriving the clusters from first principles independent of that framework, raising the risk that observed associations are partly by construction.

    Authors: The algorithmic definition of contact clusters (groups of contacts exhibiting high mutual correlation in their time series) is indeed taken from the prior framework, as the present work applies and tests that model for generality. However, clusters are identified independently and de novo in each of the new ~1 ms trajectories using the same correlation threshold and clustering procedure, without seeding or biasing toward previously observed PDZ clusters. The recurrence of analogous modules across different domains, ligands, and photoswitch perturbations therefore constitutes an out-of-sample test. We have revised the Methods section to explicitly state that cluster detection is performed independently on the current dataset and have added clarifying language in the Introduction and Results that the multistep mechanism is a hypothesis whose consistency is being evaluated rather than presupposed. revision: yes

Circularity Check

1 steps flagged

Contact cluster model and multistep picture imported via self-citation; new simulations test consistency rather than derive independently

specific steps
  1. self citation load bearing [Abstract]
    "The contact cluster model of allostery [J. Chem. Theory Comput. 2024, 20, 10731-10739] identifies localized groups of highly correlated contacts that mediate interactions between secondary structure elements. This framework proposes that allostery proceeds through a multistep process involving cooperative contact changes within clusters and communication between distant clusters, transmitted through rigid secondary structures."

    The load-bearing definition of contact clusters, their role as mediators, and the multistep communication mechanism are taken directly from the cited prior work (same author group) rather than obtained by first-principles analysis or external validation in the present study. The new simulations then map observed dynamics onto this pre-existing framework, so the claimed 'modular, dynamics-based architecture' reduces to consistency testing of an imported ansatz.

full rationale

The paper's core architecture (localized contact clusters mediating multistep allosteric communication via cooperative changes and rigid secondary structures) is explicitly attributed to a 2024 JCTC reference. This paper applies the framework to four new PDZ systems with ~1 ms of MD, reports recurring clusters and timescale associations, but does not re-derive the cluster definition or uniqueness of the multistep mechanism from the present trajectories or external controls. The central claim therefore rests on the prior self-publication for its conceptual foundation while adding empirical consistency checks. This matches a moderate self-citation load that does not fully collapse the derivation but makes the model itself non-independent within the present manuscript.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on the contact cluster definition and multistep allostery mechanism from the prior publication, plus standard assumptions of molecular dynamics (force field accuracy, sufficient sampling, and that photoswitch perturbations mimic natural signals). No new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • domain assumption Contact clusters identified by correlation analysis in MD trajectories correspond to the physical units that transmit allosteric signals.
    Invoked when the authors associate cluster motions with nonequilibrium response timescales.
  • domain assumption The multistep process of cooperative contact changes within clusters and communication between clusters is general across PDZ domains.
    Stated as the framework being demonstrated for validity and generality.

pith-pipeline@v0.9.0 · 5532 in / 1342 out tokens · 33663 ms · 2026-05-17T21:13:11.683058+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

  1. [1]

    In- trinsic dynamics of enzymes in the un- bound state and relation to allosteric reg- ulation.Curr

    (3) Bahar, I.; Chennubhotla, C.; Tobi, D. In- trinsic dynamics of enzymes in the un- bound state and relation to allosteric reg- ulation.Curr. Opin. Struct. Biol.2007, 17, 633–640. (4) Cui, Q.; Karplus, M. Allostery and coop- erativity revisited.Protein Science2008, 17, 1295–1307. (5) Changeux, J.-P. Allostery and the Monod- Wyman-Changeux Model After 50 ...

    work page 2007

  2. [2]

    J.; Nussinov, R

    (8) Tsai, C. J.; Nussinov, R. A Unified View of ”How Allostery Works”.PLoS Comput. Biol.2014,10. (9) Thirumalai, D.; Hyeon, C.; Zhu- ravlev, P. I.; Lorimer, G. H. Symme- try, Rigidity, and Allosteric Signaling: From Monomeric Proteins to Molecu- lar Machines.Chem. Rev.2019,119, 6788–6821. (10) Br¨ uschweiler, S.; Schanda, P.; Kloiber, K.; Brutscher, B.; K...

    work page 2014

  3. [3]

    Time- resolved observation of protein allosteric 11 communication.Proc

    (15) Buchenberg, S.; Sittel, F.; Stock, G. Time- resolved observation of protein allosteric 11 communication.Proc. Natl. Acad. Sci. USA2017,114, E6804–E6811. (16) Zheng, Y.; Cui, Q. Multiple Pathways and Time Scales for Conformational Transi- tions in apo-Adenylate Kinase.J. Chem. Theory Comput.2018,14, 1716–1726. (17) Ayaz, P.; Lyczek, A.; Paung, Y.; Min...

    work page 2018

  4. [4]

    Protein Allostery and Conformational Dynamics.Chem

    (18) Guo, J.; Zhou, H.-X. Protein Allostery and Conformational Dynamics.Chem. Rev.2016,116, 6503–6515. (19) Dokholyan, N. V. Controlling Allosteric Networks in Proteins.Chem. Rev.2016, 116, 6463–6487. (20) Fuentes, E.; Der, C.; Lee, A. Ligand- dependent dynamics and intramolecular signaling in a PDZ domain.J. Mol. Biol. 2004,335, 1105–1115. (21) Petit, C....

    work page 2016

  5. [5]

    J.; Domingo, J.; Schmiedel, J

    (27) Faure, A. J.; Domingo, J.; Schmiedel, J. M.; Hidalgo-Carcedo, C.; Diss, G.; Lehner, B. Mapping the ener- getic and allosteric landscapes of protein binding domains.Nature2022,604, 175–183. (28) Bozovic, O.; Ruf, J.; Zanobini, C.; Jankovic, B.; Buhrke, D.; Johnson, P. J. M.; Hamm, P. The Speed of Allosteric Signaling Within a Single-Domain Pro- tein.J...

    work page 2021

  6. [6]

    A.; Walser, R.; Donten, M

    (29) Buchli, B.; Waldauer, S. A.; Walser, R.; Donten, M. L.; Pfister, R.; Bloechliger, N.; Steiner, S.; Caflisch, A.; Zerbe, O.; Hamm, P. Kinetic response of a pho- toperturbed allosteric protein.Proc. Natl. Acad. Sci. USA2013,110, 11725–11730. (30) Bozovic, O.; Zanobini, C.; Gulzar, A.; Jankovic, B.; Buhrke, D.; Post, M.; Wolf, S.; Stock, G.; Hamm, P. Re...

    work page 2022

  7. [7]

    Perspective: Iden- tification of Collective Coordinates and Metastable States of Protein Dynamics.J

    (33) Sittel, F.; Stock, G. Perspective: Iden- tification of Collective Coordinates and Metastable States of Protein Dynamics.J. Chem. Phys.2018,149, 150901. 12 (34) Sittel, F.; Filk, T.; Stock, G. Principal component analysis on a torus: Theory and application to protein dynamics.J. Chem. Phys.2017,147, 244101. (35) Bowman, G. R., Pande, V. S., No´ e, F.,...

    work page 2018

  8. [8]

    Constructing Markov State Models to elucidate the functional conformational changes of complex biomolecules.WIREs Comp

    (36) Wang, W.; Cao, S.; Zhu, L.; Huang, X. Constructing Markov State Models to elucidate the functional conformational changes of complex biomolecules.WIREs Comp. Mol. Sci.2018,8, e1343. (37) Scherer, M. K.; Husic, B. E.; Hoff- mann, M.; Paul, F.; Wu, H.; No´ e, F. Vari- ational selection of features for molecu- lar kinetics.J. Chem. Phys.2019,150, 194108...

    work page 2018

  9. [9]

    (39) Latzer, J.; Shen, T.; Wolynes, P. G. Con- formational Switching upon Phosphoryla- tion: A Predictive Framework Based on Energy Landscape Principles.Biochem. 2008,47, 2110–2122. (40) Di Paola, L.; De Ruvo, M.; Paci, P.; San- toni, D.; Giuliani, A. Protein Contact Net- works: An Emerging Paradigm in Chem- istry.Chem. Rev.2013,113, 1598–1613. (41) Ernst...

    work page 2008

  10. [10]

    Establishing a Framework of Using Residue–Residue Interactions in Protein Difference Network Analysis.J

    (43) Yao, X.-Q.; Momin, M.; Hamelberg, D. Establishing a Framework of Using Residue–Residue Interactions in Protein Difference Network Analysis.J. Chem. Inf. Model.2019,59, 3222–3228. (44) Diez, G.; Nagel, D.; Stock, G. Correlation- based feature selection to identify func- tional dynamics in proteins.J. Chem. Theory Comput.2022,18, 5079 –

    work page 2019

  11. [11]

    (45) Ali, A. A. A. I.; Dorbath, E.; Stock, G. Al- losteric communication mediated by pro- tein contact clusters: A dynamical model. J. Chem. Theory Comput.2024,20, 10731 – 10739. (46) Stock, G.; Hamm, P. A Nonequilibrium Approach to Allosteric Communication. Phil. Trans. B2018,373, 20170187. (47) L´ orenz-Fonfr´ ıa, V. A.; Kandori, H. Trans- formation of ...

    work page 2024

  12. [12]

    Distribution of Re- ciprocal of Interatomic Distances: A Fast Structural Metric.J

    (50) Zhou, T.; Caflisch, A. Distribution of Re- ciprocal of Interatomic Distances: A Fast Structural Metric.J. Chem. Theory Com- put.2012,8, 2930 –

    work page 2012

  13. [13]

    Toward a Benchmark for Markov State Models: The Folding of HP35.J

    (51) Nagel, D.; Sartore, S.; Stock, G. Toward a Benchmark for Markov State Models: The Folding of HP35.J. Phys. Chem. Lett. 2023,14, 6956–6967. (52) Peterson, F. C.; Penkert, R. R.; Volk- man, B. F.; Prehoda, K. E. Cdc42 Reg- ulates the Par-6 PDZ Domain through 13 an Allosteric CRIB-PDZ Transition.Mol. Cell2004,13, 665–676. (53) Li, H.; Ma, A. Enhanced sa...

    work page 2023

  14. [14]

    Diffusion in rough potentials

    (54) Zwanzig, R. Diffusion in rough potentials. Proc. Natl. Acad. Sci. (USA)1988,85,

    work page 1988

  15. [15]

    M.; Smith, J

    (55) Neusius, T.; Daidone, I.; Sokolov, I. M.; Smith, J. C. Subdiffusion in Peptides Originates from the Fractal-Like Structure of Configuration Space.Phys. Rev. Lett. 2008,100, 188103. (56) Best, R. B.; Hummer, G. Coordinate- dependent diffusion in protein folding. Proc. Natl. Acad. Sci. USA2010,107, 1088 –

    work page 2008

  16. [16]

    Estimat- ing ruggedness of free-energy landscapes of small globular proteins from principal component analysis of molecular dynam- ics trajectories.Phys

    (57) Volkhardt, A.; Grubm¨ uller, H. Estimat- ing ruggedness of free-energy landscapes of small globular proteins from principal component analysis of molecular dynam- ics trajectories.Phys. Rev. E2022,105, 044404. (58) Janke, P.; Dorbath, E.; Stock, G.; Hamm, P. Perspective: Universal Relax- ation Behaviour of Photoactive Proteins. J. Chem. Phy.2025, (59...

    work page 2025

  17. [17]

    Enhanced Sampling with Ma- chine Learning.Annu

    (62) Mehdi, S.; Smith, Z.; Herron, L.; Zou, Z.; Tiwary, P. Enhanced Sampling with Ma- chine Learning.Annu. Rev. Phys. Chem. 2024,75, 347–370. (63) Abraham, M. J.; Murtola, T.; Schulz, R.; Pall, S.; Smith, J. C.; Hess, B.; Lin- dahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomput- ers.Software...

    work page 2024

  18. [18]

    L.; Dror, R

    (66) Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J. L.; Dror, R. O.; Shaw, D. E. Improved side-chain torsion potentials for the Amber ff99SB protein force field.Proteins2010,78, 1950–1958. (67) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liqui...

    work page 1950

  19. [19]

    Log- periodic oscillations as real-time signa- tures of hierarchical dynamics in proteins

    14 (71) Dorbath, E.; Gulzar, A.; Stock, G. Log- periodic oscillations as real-time signa- tures of hierarchical dynamics in proteins. J. Chem. Phys.2024,160, 074103. 15 TOC Graphic 16

    work page 2024