Open rabiakhd opened 4 days ago
This seems to be more an issue the topology. Looking at the data you sent, there are only two chains identified in the topology to begin with.
>>> np.unique(universe.atoms.chainIDs)
array(['A', 'X'], dtype=object)
There is potential to improve the ligand detection though - how would you extract ligand information from this trajectory currently?
The ligand and solvent are grouped together and identified as one chain on MD import Automatic ligand detection from protein-ligand complexes inside trajectories will be awesome.