Accounts for the interaction of both water and small molecule probes with a protein’s surface. MixMD recognizes conserved and displaceable water sites. It employs an occupancy-based metric to find sites which are consistently occupied by water even in the presence of probe molecules. This method is able to determine which functional groups are capable of displacing which water sites. It stores the displacement of water sites by common functional groups.