Hi All,
I'm trying to predict whether heteroatoms (S or O) in my ligand can form non-bonding interactions (See JACS, 2013, 135, 10475) with amino acids such as a Methionine or Cysteine in my protein active site as a possible explanation for selectivity and potency. From the docked poses it does appear that the heteroatom comes reasonably close to the sulfurs in the amino acids (around 3.5 Angstrom). Is there a modeling software that can estimate the strength of the interaction between these heteroatoms in the complex?
Thank you so much for your time and comments!