Constrained molecuar dynamics calculations