Cubic phases were self-assembled from random starting configurations of fixed MO/water ratios. The required number of MO and water molecules were randomly inserted into a 15 x 15 x 15 nm box and energy was minimized using the steepest descent algorithm for 1000 steps. The system was then equilibrated for 10ns with all nonbonded interactions set to the equivalent of the water self-interaction for the Martini 3 force field, following previously described approaches for correctly optimizing the box volume 173 . After equilibration, the system nonbonded interactions were reset to their default values to allow self-assembly to take place in a 1 μs simulation.