We have calculated ABFEs of trapped waters in protein-ligand complexes using a separation of states approach. To calculate the ABFEs, we simulated a thermodynamic cycle that implemented a harmonic restraint on the trapped water, a solvent repulsion term, and position restraints on the protein's binding site (Figure ). Once the restraints were applied, we decoupled the trapped water using an NES simulation, and then removed the restraints.