Hybrid Dynamical Simulation Reveals Apparent Stiffening of Flexible Protein Lattices Driving Membrane Bending
Abstract
Membrane-deforming protein lattices play a central role in essential and pathogenic remodeling processes, including clathrin-mediated endocytosis and viral budding. Simulating these systems at biologically relevant length and time scales requires mesoscale approaches that preserve structural detail. We present a hybrid simulation framework that couples a particle-based flexible protein lattice to a continuum membrane model, enabling systematic investigation of how lattice geometry and rigidity influence membrane remodeling. We validate the coupled model by comparing simulation results with theoretical predictions for membranes under increasing tension. Using buckling-based deformations of pre-assembled clathrin lattices, we quantify the lattice flexural rigidity and establish a direct relationship between the force constants in the coarse-grained energy function and the emergent lattice mechanics. We then compare this flexural rigidity, κl,flex, to an effective rigidity, κl,eff, commonly used in continuum descriptions of sphere-forming protein assemblies. Although κl,flex is set solely by the energy function, κl,eff depends on lattice size and connectivity, with the two measures converging for weakly-connected lattices. As a result, with lattice growth, the effective rigidity relevant for spherical bud formation increases and the relaxation timescale decreases. This size-dependent stiffening highlights the importance of structural details in interpreting lattice mechanics and cautions against assuming a single constant stiffness. We demonstrate the generality of the method by applying it to pre-assembled viral lattices generated with reaction-diffusion simulations. This work provides a validated framework for simulating deformable protein assemblies of diverse geometry coupled to membrane dynamics and remodeling.
Related
- Dynamical superfluid and Bose-insulator phases in quantized microcavity lattices
- Interference-induced dynamical crossovers in quantum walks on disordered square lattices
- diff-soup: GPU-native Differentiable Simulation and Inference for Large-scale Stiff Dynamical Systems
- Optimizing quantum encodings for analog simulation through dynamical algebra reachability
- Extraordinary stiffening and reciprocal control in a DNA condensate/actin composite
BibTeX
@misc{foley2026hybrid,
title = {{Hybrid Dynamical Simulation Reveals Apparent Stiffening of Flexible Protein Lattices Driving Membrane Bending}},
author = {Samuel L. Foley and Margaret E. Johnson},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2607.06378},
eprint = {2607.06378},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.06378},
}
RIS
TY - UNPB
TI - Hybrid Dynamical Simulation Reveals Apparent Stiffening of Flexible Protein Lattices Driving Membrane Bending
AU - Samuel L. Foley
AU - Margaret E. Johnson
PY - 2026
DA - 2026/10/27
PB - arXiv
DO - 10.48550/arxiv.2607.06378
UR - https://arxiv.org/abs/2607.06378
AB - Membrane-deforming protein lattices play a central role in essential and pathogenic remodeling processes, including clathrin-mediated endocytosis and viral budding. Simulating these systems at biologically relevant length and time scales requires mesoscale approaches that preserve structural detail. We present a hybrid simulation framework that couples a particle-based flexible protein lattice to a continuum membrane model, enabling systematic investigation of how lattice geometry and rigidity influence membrane remodeling. We validate the coupled model by comparing simulation results with theoretical predictions for membranes under increasing tension. Using buckling-based deformations of pre-assembled clathrin lattices, we quantify the lattice flexural rigidity and establish a direct relationship between the force constants in the coarse-grained energy function and the emergent lattice mechanics. We then compare this flexural rigidity, κ_l,flex, to an effective rigidity, κ_l,eff, commonly used in continuum descriptions of sphere-forming protein assemblies. Although κ_l,flex is set solely by the energy function, κ_l,eff depends on lattice size and connectivity, with the two measures converging for weakly-connected lattices. As a result, with lattice growth, the effective rigidity relevant for spherical bud formation increases and the relaxation timescale decreases. This size-dependent stiffening highlights the importance of structural details in interpreting lattice mechanics and cautions against assuming a single constant stiffness. We demonstrate the generality of the method by applying it to pre-assembled viral lattices generated with reaction-diffusion simulations. This work provides a validated framework for simulating deformable protein assemblies of diverse geometry coupled to membrane dynamics and remodeling.
ER -