Abstract
We report a new approach that combines molecular dynamics trajectories with time-dependent linear response theory to compute the time evolution of residue fluctuation responses to force perturbations exerted at functional sites. Applying this new approach to TEM-1 beta-lactamase, we observe that the time-resolved response profiles of allosteric sites to perturbations of TEM-1 active sites are distinct from those of non-allosteric residues. Using Fourier transformations, we convert the time domain response profiles to the frequency domain and demonstrate that the frequency space representation of the perturbation response can capture the mutational behavior of each site when applied to deep sequencing mutational data. Furthermore, we show that classification models built on perturbation responses can accurately identify distal positions that regulate antibiotic resistance. These findings provide insights into the contributions of specific residues to resistance-encoded in time-resolved perturbation response behavior and highlight the importance of this new approach in identifying allosteric mutations, opening avenues for the potential characterization of additional allosteric positions without extensive computational simulations.