How the ligand-induced reorganization of protein internal energies is coupled to conformational events

Here, we introduce a novel computational method to identify the protein substructures most likely to support the functionally oriented structural deformations that occur upon ligand-binding. To this aim, we study the modulation of protein energetics along the trajectory of a molecular dynamics simulation of different proteins in the presence and in the absence of their respective ligands, namely, human FGF, human second PDZ from human PTP1E/PTPL1, and the N terminal domain of human Hsp90. The method is based on the idea that a subset of protein residues (hotspots) may initiate the global response via the disassembly and reassembly of interactions, which is reflected in the modulation of the overall protein energetics. To identify structural hotspots and dynamic states linked to the onset of functionally relevant conformational transitions, we define an energy profile to monitor the protein energetics, based on a previously introduced approach that highlights the essential nonbonded couplings among all res...