Abstract

In this work, we have employed docking and atomistic molecular dynamics (MD) simulations supported by complementary experiments using atomic force microscopy, rheology and spectroscopy to investigate the self-assembled structure of β-sitosterol and γ-oryzanol molecules into cylindrical tubules in a non-aqueous solvent. Docking models of several phytosterols, including sitosterol, with oryzanol and other sterol-esters demonstrate that for systems to form tubules, the phytosterol sterane group must be stacked in a wedge shape with the esters sterane group, and a hydrogen bond must form between the hydroxyl group of the phytosterol and the carbonyl group of the ester. Molecular dynamics of the self-assembled structure were initiated with the molecules in a roughly cylindrical configuration, as suggested from previous experimental studies, and the configurations were found to be stable during 50 ns simulations. We performed MD simulations of two tubules in proximity to better understand the aggregation of these fibrils and how the fibrils interact in order to stick together. We found that an interfibril network of non-covalent bonds, in particular van der Waals and π-π contacts, which is formed between the ferulic acid groups of oryzanol through the hydroxyl, methoxy and aromatic groups, is responsible for the surface-to-surface interactions between fibrils; an observation supported by molecular spectroscopy. We believe these interactions are of primary importance in creating a strong organogel network.