Abstract
Orientation of ten water molecules bound strongly at the contact surface of the dihydrofolate reductase-methotrexate enzyme-inhibitor complex was determined theoretically. To optimize the orientation of the water molecules, a recent method based on a simple electrostatic model was applied. The electrostatic complementarity in the binary complex was investigated using the lock-and-key model, considering the effect of the water molecules as well. The strongly bound water molecules improve the electrostatic fit in the pteridine region of methotrexate. Their role in the benzoic amide andγ-glutamate region is to decrease the internal energy by creating water bridges among remote polar sites making it possible to form H-bonds. Some modifications in the inhibitor structure were proposed for achieving greater inhibitor potency. The presumably enhanced effect is ascribed to the free energy gain in repelling the water molecules from the contact surface to the bulk of the solvent, and, in other cases, to internal energy decreases due to better electrostatic fit in the enzyme-inhibitor complex.

This publication has 39 references indexed in Scilit: