An azobenzene-derived photoswitch has been covalently cross-linked to two sites of the Speptide in the RNase S complex in a manner that the alpha-helical content of the S-peptide reduces upon cis-to-trans isomerization of the photoswitch. Three complementary experimental techniques have been employed, isothermal titration calorimetry, circular dichroism spectroscopy and intrinsic tyrosine fluorescence quenching, to determine the binding affinity of the S-peptide to the S-protein in the two states of the photoswitch. Five mutants with the photoswitch attached to different sites of the S-peptide have been explored, with the goal to maximize the change in binding affinity upon photoswitching, and to identify the mechanisms that determine the binding affinity. With regard to the first goal, one mutant has been identified, which binds with reasonable affinity in the one state of the photoswitch, while specific binding is completely switched off in the other state. With regard to the second goal, accompanying molecular dynamics simulations combined with a quantitative structure activity relationship revealed that the alpha-helicity of the S-peptide in the binding pocket correlates surprisingly well with measured dissociation constants. Moreover, the simulations show that both configurations of all S-peptides exhibit quite well-defined structures, even in apparently disordered states.

While allostery is of paramount importance for protein regulation, the underlying dynamical process of ligand (un)binding at one site, resulting time evolution of the protein structure, and change of the binding affinity at a remote site are not well understood. Here the ligand-induced conformational transition in a widely studied model system of allostery, the PDZ2 domain, is investigated by transient infrared spectroscopy accompanied by molecular dynamics simulations. To this end, an azobenzene-derived photoswitch is linked to a peptide ligand in a way that its binding affinity to the PDZ2 domain changes upon switching, thus initiating an allosteric transition in the PDZ2 domain protein. The subsequent response of the protein, covering four decades of time, ranging from ∼1 ns to ∼μs, can be rationalized by a remodeling of its rugged free-energy landscape, with very subtle shifts in the populations of a small number of structurally well-defined states. It is proposed that structurally and dynamically driven allostery, often discussed as limiting scenarios of allosteric communication, actually go hand-in-hand, allowing the protein to adapt its free-energy landscape to incoming signals.

Recent time-resolved experiments and accompanying molecular dynamics simulations allow us to monitor the flow of vibrational energy in biomolecules. As a simple means to describe these experimental and simulated data, Buchenberg et al. [J. Phys. Chem. Lett. 7, 25 (2016)] suggested a master equation model that accounts for the energy transport from an initially excited residue to some target residue. The transfer rates of the model were obtained from two scaling rules, which account for the energy transport through the backbone and via tertiary contacts, respectively, and were parameterized using simulation data of a small α-helical protein at low temperatures. To extend the applicability of the model to general proteins at room temperature, here a new parameterization is presented, which is based on extensive nonequilibrium molecular dynamics simulations of a number of model systems. With typical transfer times of 0.5–1 ps between adjacent residues, backbone transport represents the fastest channel of energy flow. It is well described by a diffusive-type scaling rule, which requires only an overall backbone diffusion coefficient and interatom distances as input. Contact transport, e.g., via hydrogen bonds, is considerably slower (6–30 ps) at room temperature. A new scaling rule depending on the inverse square contact distance is suggested, which is shown to successfully describe the energy transport in the allosteric protein PDZ3. Since both scaling rules require only the structure of the considered system, the model provides a simple and general means to predict energy transport in proteins. To identify the pathways of energy transport, Monte Carlo Markov chain simulations are performed, which highlight the competition between backbone and contact transport channels.

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.

The noncanonical amino acid azidohomoalanine (Aha) is known to be an environment-sensitive infrared probe for the site-specific investigation of protein structure and dynamics. Here, the capability of that label is explored to detect protein ligand interactions by incorporating it in the vicinity of the binding groove of a PDZ2 domain. Circular dichroism and isothermal titration calorimetry measurements reveal that the perturbation of the protein system by mutation is negligible, with minimal influence on protein stability and binding affinity. Two-dimensional infrared spectra exhibit small (1-3 cm(-1)) but clearly measurable red shifts of the Aha vibrational frequency upon binding of two different peptide ligands, while accompanying molecular dynamics simulations suggest that these red shifts are induced by polar contacts with side chains of the peptide ligands. Hence, Aha is a versatile and minimally invasive vibrational label that is not only able to report on large structural changes during, e.g., protein folding, but also on very subtle changes of the electrostatic environment upon ligand binding.