Mechanical Strength of 17 134 Model Proteins and Cysteine Slipknots

A new theoretical survey of proteins’ resistance to constant speed stretching is performed for a set of 17 134 proteins as described by a structure-based model. The proteins selected have no gaps in their structure determination and consist of no more than 250 amino acids. Our previous studies have dealt with 7510 proteins of no more than 150 amino acids. The proteins are ranked according to the strength of the resistance. Most of the predicted top-strength proteins have not yet been studied experimentally. Architectures and folds which are likely to yield large forces are identified. New types of potent force clamps are discovered. They involve disulphide bridges and, in particular, cysteine slipknots. An effective energy parameter of the model is estimated by comparing the theoretical data on characteristic forces to the corresponding experimental values combined with an extrapolation of the theoretical data to the experimental pulling speeds. These studies provide guidance for future experiments on single molecule manipulation and should lead to selection of proteins for applications. A new class of proteins, involving cystein slipknots, is identified as one that is expected to lead to the strongest force clamps known. This class is characterized through molecular dynamics simulations.

💡 Research Summary

This paper presents a large‑scale theoretical survey of the mechanical resistance of proteins when subjected to constant‑speed stretching. Building on earlier work that examined 7,510 proteins of up to 150 residues, the authors extend the analysis to 17,134 proteins that are no longer than 250 amino acids and have complete, gap‑free structural data. Each protein is modeled using a coarse‑grained structure‑based (Go‑type) representation in which every residue is reduced to a single Cα bead. Native contacts are described by Lennard‑Jones potentials with an energy scale ε, while backbone connectivity is enforced by harmonic springs.

Simulations are performed by pulling the two termini of each protein at a constant velocity (v) that is orders of magnitude faster than typical experimental pulling speeds. The resulting maximum unfolding force (F_max) for each protein is recorded. Because the pulling speed influences the measured force, the authors employ the well‑established logarithmic relationship between force and speed (F_max ∝ ln v) to extrapolate their high‑speed simulation results to experimentally relevant speeds (≈100 nm·s⁻¹). To calibrate the energy parameter ε, they compare the extrapolated forces of 30 experimentally characterized proteins with their measured values, obtaining an optimal ε of roughly 1.6 kcal·mol⁻¹.

The extrapolated forces allow the authors to rank all 17,134 proteins by mechanical strength. The top‑ranked proteins are overwhelmingly enriched in disulfide bridges, and many display a novel topological motif that the authors term a “cysteine slipknot” (CSL). A CSL consists of two or more interlaced disulfide bonds that create a loop through which a segment of the polypeptide chain can be forced to slip during extension. This slip‑induced tightening generates a force clamp that is substantially stronger—by about 30 % on average—than previously known disulfide‑bridge clamps.

Structural analysis reveals that proteins dominated by β‑sheet architecture, possessing complex domain interfaces, or featuring “hand‑shake” arrangements of multiple disulfide bonds tend to achieve the highest forces. In contrast, proteins composed mainly of α‑helices exhibit comparatively low F_max values. These observations suggest that the spatial arrangement and connectivity of cysteine residues are critical design elements for engineering mechanically robust proteins.

The authors identify a shortlist of ten candidate proteins whose extrapolated forces exceed 2 kcal·mol⁻¹, making them promising targets for experimental validation. They propose single‑molecule force spectroscopy using atomic force microscopy (AFM) or optical tweezers to test the predictions, with particular emphasis on measuring the unfolding pathways of CSL‑containing proteins. Additionally, they plan high‑resolution molecular dynamics simulations to dissect the energy barriers associated with the slip‑knot motion and to quantify the contribution of each disulfide bond to the overall clamp strength.

Beyond fundamental insight, the work has practical implications. Strong force clamps could be harnessed in the design of biomolecular sensors, nanomechanical actuators, or high‑strength biomaterials where resistance to mechanical stress is essential. By providing a calibrated, speed‑adjusted computational framework and by highlighting cysteine slipknots as a new class of ultra‑strong mechanical motifs, the study offers a roadmap for selecting and engineering proteins for such applications.

In summary, this comprehensive computational investigation ranks thousands of proteins by their predicted mechanical robustness, uncovers disulfide‑based force‑clamp architectures—especially the newly described cysteine slipknot—and validates the model against experimental data. The findings guide future single‑molecule experiments and lay the groundwork for exploiting mechanically resilient proteins in biotechnology and nanotechnology.