Janssen Ligand tutorial: why no -H2Hpolar flag in atoms_to_morph.py?

Dear all,

I am doing the Janssen ligand tutorial: Ligand modification free energy calculations. During the atom mapping step using atoms_to_morph.py the following command is used:

$ python atoms_to_morph.py -i1 lig_1a/lig_1a.pdb -i2 lig_5/lig_5.pdb -o pairs.dat -alignment -H2H

I however noticed that this command introduces dummy hydrogens also on the nitrogens by not pairing the polar hydrogens. Upon testing a few flags for the atoms_to_morph.py script, I found that these polar hydrogen dummy atoms can be removed using the -H2Hpolar flag for the atoms_to_morph.py script. Why does the tutorial propose to add dummy atoms for polar hydrogens? Is there a specific reason for this?

Thank you in advance!

Kind regards,


I usually let polar hydrogens go via dummy transition, because these atoms might be involved in forming hydrogen bonds, i.e. I try to avoid hydrogen positioning in one molecule based on the coordinates of another. It should have minimal effect on the perturbation size, as there is another hydrogen occupying the same space in another molecule.