Fiber formation seen through the high-resolution computational microscope

Supramolecular fibers draw widespread attention due to their role in biological systems and ability to form complex materials exhibiting rich and dynamic behavior. Although the information about the supramolecular structure is encoded in their molecular blocks, a complete understanding of how this information translates into the final structure requires detailed insights into the energy landscape of the process and the possible routes across this landscape. Here, we study the formation of 1,3,5-cyclohexanetricarboxamide fibers by a Markov state model of molecular dynamics simulations with the polarizable CHARMM Drude force-field. We provide insights into all stages of supramolecular fiber formation up to microsecond timescales, starting from primary nucleation, through fiber elongation and secondary nucleation, and finally, the assembly of single fibers into bundles. Our results show that nucleation involves a rapid collapse of dissolved monomers into disordered assemblies, which gradually transform into nuclei and then grow into elongated fibers. Moreover, elongation and secondary nucleation appeared to be competing processes, depending on the density of the monomers adsorbed at the fiber-liquid interface. Finally, bundling involves the initial association of fibers by interactions between surface-exposed groups, followed by stabilization of the bundle by surface reorganization, which allows for favorable interactions between aromatic groups.