Thank you Justin I don't understand. Do I use the explicit water molecules or the model? and where do I used the "define = -DFLEXIBLE" option? please explain a bit more!