Elastic deformation and failure in protein filament bundles: Atomistic simulations and coarse-grained modeling

The synthetic peptide RAD16-II has shown promise in tissue engineering and drug delivery. It has been studied as a vehicle for cell delivery and controlled release of IGF-1 to repair infarcted cardiac tissue, and as a scaffold to promote capillary formation for an in vitro model of angiogenesis. The structure of RAD16-II is hierarchical, with monomers forming long β-sheets that pair together to form filaments; filaments form bundles approximately 30–60 nm in diameter; branching networks of filament bundles form macroscopic gels. We investigate the mechanics of shearing between the two β-sheets constituting one filament, and between cohered filaments of RAD16-II. This shear loading is found in filament bundle bending or in tensile loading of fibers composed of partial-length filaments. Molecular dynamics simulations show that time to failure is a stochastic function of applied shear stress, and that for a given loading time behavior is elastic for sufficiently small shear loads. We propose a coarse-grained model based on Langevin dynamics that matches molecular dynamics results and facilities extending simulations in space and time. The model treats a filament as an elastic string of particles, each having potential energy that is a periodic function of its position relative to the neighboring filament. With insight from these simulations, we discuss strategies for strengthening RAD16-II and similar materials.