The role of flow in the self-assembly of dragline spider silk proteins.

Hydrodynamic flow in the spider duct induces conformational changes in dragline spider silk proteins (spidroins) and drives their assembly, but the underlying physical mechanisms are still elusive. Here we address this challenging multiscale problem with a complementary strategy of atomistic and coarse-grained molecular dynamics simulations with uniform flow. The conformational changes at the molecular level were analyzed for single-tethered spider silk peptides. Uniform flow leads to coiled-to-stretch transitions and pushes alanine residues into ß sheet and poly-proline II conformations. Coarse-grained simulations of the assembly process of multiple semi-flexible block copolymers using multi-particle collision dynamics reveal that the spidroins aggregate faster but into low-order assemblies when they are less extended. At medium-to-large peptide extensions (50%-80%), assembly slows down and becomes reversible with frequent association and dissociation events, whereas spidroin alignment increases and alanine repeats form ordered regions. Our work highlights the role of flow in guiding silk self-assembly into tough fibers by enhancing alignment and kinetic reversibility, a mechanism likely relevant also for other proteins whose function depends on hydrodynamic flow.

Molecular Dynamics Simulations of Molecules in Uniform Flow.

Flow at the molecular level induces shear-induced unfolding of single proteins and can drive their assembly, the mechanisms of which are not completely understood. To be able to analyze the role of flow on molecules, we present uniform-flow molecular dynamics simulations at atomic level. The pull module of the GRoningen MAchine for Chemical Simulations package was extended to be able to force-group atoms within a defined layer of the simulation box. Application of this external enforcement to explicit water molecules, together with the coupling to a thermostat, led to a uniform terminal velocity of the solvent water molecules. We monitored the density of the whole system to establish the conditions under which the simulated flow is well-behaved. A maximal velocity of 1.3 m/s can be generated if a pull slice of 8 nm is used, and high velocities would require larger pull slices to still maintain a stable density. As expected, the target velocity increases linearly with the total external force applied. Finally, we suggest an appropriate setup to stretch a protein by uniform flow, in which protein extensions depend on the flow conditions. Our implementation provides an efficient computational tool to investigate the effect of the flow at the molecular level.