Thanks for the advice. Yes, I am working on a not very common case. In my lab, we are designing protein molecules that can assemble into cages (
Science 2012, 336(6085) p1129). This method was used successfully to generate tetrahedral cages and now I am trying to make icosahedral cages. At the design stage, we fuse in silico a natural C2 protein dimer (taken from PDB) to a natural C5 protein pentamer (also from PDB) with a long helix linker. In this way a fusion molecule with both 5-fold and 2-fold symmetries can be generated. However, such model won't be oriented in any specific ways that can be used with the sym command. One way to fix this is to write a python script to align the 5-fold axis of the fusion molecule to the z-axis and then align the 2-fold axis (not of the molecule, need to be derived) to the y-axis. Then the sym command can be used by specifying n25 as the orientation. I am just wondering if there's any quick way of doing this.