Can molecular dynamics simulations predict the mechanical behavior of protein complexes? Can simulations decipher the role of protein domains of unknown function in large macromolecular complexes? Here, we employ a wide-sampling computational approach to demonstrate that molecular dynamics simulations, when carefully performed and combined with single-molecule atomic force spectroscopy experiments, can predict and explain the behavior of highly mechanostable protein complexes. As a test case, we studied a previously unreported homologue from Ruminococcus flavefaciens called X-module-Dockerin (XDoc) bound to its partner Cohesin (Coh). By performing dozens of short simulation replicas near the rupture event, and analyzing dynamic network fluctuations, we were able to generate large simulation statistics and directly compare them with experiments to uncover the mechanisms involved in mechanical stabilization. Our single-molecule force spectroscopy experiments show that the XDoc-Coh homologue complex withstands forces up to 1 nN at loading rates of 105 pN/s. Our simulation results reveal that this remarkable mechanical stability is achieved by a protein architecture that directs molecular deformation along paths that run perpendicular to the pulling axis. The X-module was found to play a crucial role in shielding the adjacent protein complex from mechanical rupture. These mechanisms of protein mechanical stabilization have potential applications in biotechnology for the development of systems exhibiting shear enhanced adhesion or tunable mechanics.

Bacterial colonization of the human intestine requires firm adhesion of bacteria to insoluble substrates under hydrodynamic flow. Here we report the molecular mechanism behind an ultrastable protein complex responsible for resisting shear forces and adhering bacteria to cellulose fibers in the human gut. Using single-molecule force spectroscopy (SMFS), singlemolecule FRET (smFRET), and molecular dynamics (MD) simulations, we resolve two binding modes and three unbinding reaction pathways of a mechanically ultrastable R. champanellensis (Rc) Dockerin:Cohesin (Doc:Coh) complex. The complex assembles in two discrete binding modes with significantly different mechanical properties, with one breaking at 500 pN and the other at~200 pN at loading rates from 1-100 nN s −1. A neighboring X-module domain allosterically regulates the binding interaction and inhibits one of the lowforce pathways at high loading rates, giving rise to a catch bonding mechanism that manifests under force ramp protocols. Multi-state Monte Carlo simulations show strong agreement with experimental results, validating the proposed kinetic scheme. These results explain mechanistically how gut microbes regulate cell adhesion strength at high shear stress through intricate molecular mechanisms including dual-binding modes, mechanical allostery and catch bonds.

We used single-molecule AFM force spectroscopy (AFM-SMFS) in combination with click chemistry to mechanically dissociate anticalin, a non-antibody protein binding scaffold, from its target (CTLA-4), by pulling from eight different anchor residues. We found that pulling on the anticalin from residue 60 or 87 resulted in significantly higher rupture forces and a decrease in k off by 2–3 orders of magnitude over a force range of 50–200 pN. Five of the six internal anchor points gave rise to complexes significantly more stable than N- or C-terminal anchor points, rupturing at up to 250 pN at loading rates of 0.1–10 nN s–1. Anisotropic network modeling and molecular dynamics simulations helped to explain the geometric dependency of mechanostability. These results demonstrate that optimization of attachment residue position on therapeutic binding scaffolds can provide large improvements in binding strength, allowing for mechanical affinity maturation under shear stress without mutation of binding interface residues.