ABSTRACT
Elastic proteins and derived biomaterials contain numerous tandemly repeated peptides along their sequences, ranging from a few copies to hundreds. These repetitions are responsible for their biochemical, biological and biomechanical properties. These sequences are considered to be intrinsically disordered, and the variations in their behavior are actually mainly due to their high flexibility and lack of stable secondary structures originating from their unique amino acid sequences. Consequently, the simulation of elastic proteins and large elastomeric biomaterials using classical molecular dynamics is an important challenge. Here, we propose a novel approach that allows the application of the DURABIN protocol to repeated elastin-like peptides (r-ELPs) in a simple way. Four large r-ELPs were studied to evaluate our method, which was developed for simulating extracellular matrix proteins at the mesoscopic scale. After structure clustering applied on molecular dynamic trajectories of constitutive peptides (5-mers and 6-mers), the main conformations were used as starting points to define the corresponding primitives, further used as rigid body fragments in our program. Contributions derived from electrostatic and molecular hydrophobicity potentials were tested to evaluate their influence on the interactions during simple mesoscopic simulations. The CHLORAINE approach, despite the thinner granularity due to the size of the patterns used, was included in the DURABIN protocol and emerges as a promising way to simulate elastic macromolecular systems.
