Druglike Molecules Binding to Large Membrane Proteins: Absolute Binding Free Energy Computation

In this research, we employed the alchemical double-decoupling method alongside restraining potentials, coupled with the FEPMD method, to ascertain the standard binding free energy of a drug-like molecule termed BHQ and three analogous compounds engineered with progressive addition of bulky para-alkyl groups binding to SERCA (Ca²⁺-ATPase of skeletal muscle sarcoplasmic reticulum). Integral transmembrane proteins represent crucial drug targets in numerous therapeutic interventions, presenting computational challenges due to their considerable system sizes. Our approach integrated the generalized born potential method and the spherical solvent boundary potential method, allowing us to explicitly focus on the active binding site while treating the remainder of the system implicitly. We evaluated contributions to the standard binding free energy from distinct interaction potentials: electrostatic, repulsive, dispersive, and restraining potentials, computed separately. The resulting absolute binding free energy of BHQ (11.63 kcal/mol) closely aligns with experimental measurements (10.56 kcal/mol). Notably, an accurate estimation of the absolute binding free energy was achieved for the simplest analog, created with the addition of a single para-methyl group. However, the analog with two para-methyl groups exhibited the highest binding free energy, which disagreed with experimental results. Determining the binding free energy of the BHQ analog engineered with three para-methyl groups presented challenges in convergence and resulted in the lowest free energy among the three.