Cryo-EM Informed Molecular Dynamics Simulations to Investigate the Disulfide Hydrogel Self-Assembly.

Abstract

Disulfide hydrogels, derived from cysteine-based redox systems, exhibit active self-assembly properties driven by reversible disulfide bond formation, making them a versatile platform for dynamic material design. Detailed cryogenic electron microscopy (cryo-EM) analysis revealed a consistent fiber diameter of 5.4 nm for individual fibers. Using cryo-EM-informed radial positional restraints, all-atom molecular dynamics (MD) simulations were employed to reproduce fibers with dimensions closely matching experimental observations, validated further through simulated cryo-EM images. The MD simulations revealed that the disulfide gelator (CSSC) predominantly adopts an open conformation, with hydrogen bonds emerging as the key intermolecular force stabilizing the fibers. Notably, intermolecular interactions were found to be higher at 70\% conversion to the disulfide gelator compared to 100\%, comparable with past unrestrained simulations. Water molecules and solute-water hydrogen bonds are present throughout the fiber, indicating that the fiber remains hydrated. These findings underscore the potential role of the thiol precursor CSH in stabilizing the transient phase and highlight the importance of CSH-CSSC interplay. This study provides novel insights into molecular mechanisms governing self-assembly and offers strategies for designing tunable materials through controlled assembly conditions.