But if you could simulate the environment just around the ligand, then you could calculate H-bond and angles energies in a quantum perspective. The limitation is the lack of integration time, but that can be avoided using a frame that corresponds to the average system pose. Or am I too wrong? -- Marcelo Depólo