Amber-MD / cpptraj

Biomolecular simulation trajectory/data analysis.
Other
138 stars 64 forks source link

Improve atom mapping between chiral centers #1111

Closed drroe closed 1 month ago

drroe commented 1 month ago

Version 6.29.5.

Previously when mapping atoms bonded to chiral centers, atommap would use a simple cutoff scheme comparing the improper dihedral angle formed by (M-C-M-U), where M is a mapped atom, C is the chiral center, and U is an unmapped atom. This cutoff had been intentionally chosen to be tight (10 degrees) to prevent bad mapping, but if one structure was distorted with respect to the other this could cause the mapping to fail. Now in the most common chiral case (chiral center, 2 mapped atoms, 2 unmapped atoms) the sign of each angle is used instead of a delta, i.e. the negative angle in reference is mapped to the negative angle in target and vice versa. There is still a warning printed if the delta between the angles is outside the cutoff indicating the structures may be distorted.

Also added a keyword to atommap, chiralimpcut <cut> that can be used to increase the default cutoff to beyond 10 degrees. Update the manual.