zetadin
commented
1 year ago The topology files contain two [ system ] and [ molecules ] each. The first one seems to be inserted into the middle of the dihedral definitions. When parsing of [ molecules ], ligandHybrid doesn't know what to make of the subsequent lines that are supposed to be in [ dihedrals ]. Remove this inserted section, and try again.
34 30 31 36 1 0 0.650844 3 34 30 31 37 1 0 0.650844 3; remove this [ system ] BioSimSpace System
[ molecules ] ;molecule name nr. MOL 1
1 ; end remove 19 18 29 32 9 180 10.46 2 20 8 10 25 1 180 15.167 2 21 3 9 24 1 180 15.167 2
noahharrison64
Hi,
I'm attempting to run non-eq FEP and using some PMX tools to help generate the topologies to do so. Using ligandHybrid fails when I include a certain single ligand (ejm46) in the pair, while working fine with all other ligand pairs. The command I'm using is:
'pmx ligandHybrid -i1 ejm42.pdb -i2 ejm46.pdb -itp1 ejm42.top -itp2 ejm46.top -pairs pairs_a.dat -oA merged.pdb -oitp merged.itp -offitp ffmerged.itp'The resulting error is as follows:
Input files are here. Is there any reason this might be happening? Have checked the top files and the pdb files but can't see any issues. Again, since this works for other pairs including ejm42, but always fails if ejm46 is included in a pair, it's safe to say this issue is specific to ejm46 input data.