Mechanical response of silk crystalline units from force-distribution analysis

The outstanding mechanical toughness of silk fibers is thought to be caused by embedded crystalline units acting as cross links of silk proteins in the fiber. Here, we examine the robustness of these highly ordered beta-sheet structures by molecular dynamics simulations and finite element analysis. Structural parameters and stress-strain relationships of four different models, from spider and Bombyx mori silk peptides, in antiparallel and parallel arrangement, were determined and found to be in good agreement with x-ray diffraction data. Rupture forces exceed those of any previously examined globular protein many times over, with spider silk (poly-alanine) slightly outperforming Bombyx mori silk ((Gly-Ala)(n)). All-atom force distribution analysis reveals both intrasheet hydrogen-bonding and intersheet side-chain interactions to contribute to stability to similar extent. In combination with finite element analysis of simplified beta-sheet skeletons, we could ascribe the distinct force distribution pattern of the antiparallel and parallel silk crystalline units to the difference in hydrogen-bond geometry, featuring an in-line or zigzag arrangement, respectively. Hydrogen-bond strength was higher in antiparallel models, and ultimately resulted in higher stiffness of the crystal, compensating the effect of the mechanically disadvantageous in-line hydrogen-bond geometry. Atomistic and coarse-grained force distribution patterns can thus explain differences in mechanical response of silk crystals, opening up the road to predict full fiber mechanics.