Kinetic response of a photoperturbed allosteric protein

Tuning energy transport in helical protein nanotubes through side-chain modifications

Fibrous proteins are widely used as materials due to their biocompatibility, flexibility, and mechanical properties. With advancements in bioelectronics and flexible materials, there is increasing demand for biocompatible materials with tunable thermal conductivity. Understanding the mechanisms of thermal transport in proteins can facilitate the design of biomaterials with tailored thermal properties. In this study, we use non-equilibrium molecular dynamics (NEMD) to investigate how side-chain mass affects thermal transport in α-helix proteins. We analyze four representative residues - glycine (G), alanine (A), leucine (L), and phenylalanine (F) - and demonstrate that variations in side-chain mass significantly influence thermal conductivity. Results show that heavier side chains hinder heat transport, while lighter side chains enhance it. Phonon analysis reveals that side-chain mass primarily affects the properties of low-frequency acoustic and semi-optical phonons, which are critical for energy transfer. These findings provide insights into the design of protein-based biomaterials with customized thermal properties, offering potential applications in bioelectronics, medical devices, and sustainable materials. STATEMENT OF SIGNIFICANCE: This research explores how side chains in α-helix proteins influence their thermal conductivity through the application of molecular dynamics simulations. By analyzing four types of amino acids with differing side-chain masses, the study demonstrates that lighter side chains enhance heat transport, whereas heavier ones diminish it. This work establishes a direct correlation between protein structural features and their thermal properties, providing the groundwork that could enable the engineering of biomaterials with tailored heat conduction capabilities. The findings have implications for applications in bioelectronics, medical devices, and sustainable materials, where precise thermal management is essential, rendering this research highly relevant to scientists and engineers focused on advancing biocompatible materials with specific thermal characteristics.

Uncovering the Role of Distal Regions in PDK1 Allosteric Activation

Allosteric regulation is a pivotal mechanism governing a wide array of cellular functions. Essential to this process is a flexible biomolecule allowing distant sites to interact through coordinated or sequential conformational shifts. Phosphoinositide-dependent kinase 1 (PDK1) possesses a conserved allosteric binding site, the PIF-pocket, which regulates the kinase's ATP binding, catalytic activity, and substrate interactions. We elucidated the allosteric mechanisms of PDK1 by comparing conformational ensembles of the kinase bound with different small-molecule allosteric modulators in the PIF-pocket with that of the modulator-free kinase. Analysis of over 48 μs of simulations consistently shows that the allosteric modulators predominantly influence the conformational dynamics of specific distal regions from the PIF-pocket, driving allosteric activation. Furthermore, a recently developed advanced difference contact network community analysis is employed to elucidate allosteric communications. This approach integrates multiple conformational ensembles into a single community network, offering a valuable tool for future studies aimed at identifying function-related dynamics in proteins.

Enhanced sampling of protein conformational changes via true reaction coordinates from energy relaxation

The bottleneck in enhanced sampling lies in finding collective variables that effectively accelerate protein conformational changes; true reaction coordinates that accurately predict the committor are the well-recognized optimal choice. However, identifying them requires unbiased natural reactive trajectories, which, paradoxically, require effective enhanced sampling. Using the generalized work functional method, we uncover that true reaction coordinates control both conformational changes and energy relaxation, enabling us to compute them from energy relaxation simulations. Biasing true reaction coordinates accelerates conformational changes and ligand dissociation in PDZ2 domain and HIV-1 protease by 105 to 1015-fold. The resulting trajectories follow natural transition pathways, enabling efficient generation of unbiased reactive trajectories. In contrast, biased trajectories from empirical collective variables display non-physical features. Furthermore, our method uses a single protein structure as input, enabling predictive sampling of conformational changes. These findings unlock access to a broader range of protein functional processes in molecular dynamics simulations.

Log-periodic oscillations as real-time signatures of hierarchical dynamics in proteins

The time-dependent relaxation of a dynamical system may exhibit a power-law behavior that is superimposed by log-periodic oscillations. D. Sornette [Phys. Rep. 297, 239 (1998)] showed that this behavior can be explained by a discrete scale invariance of the system, which is associated with discrete and equidistant timescales on a logarithmic scale. Examples include such diverse fields as financial crashes, random diffusion, and quantum topological materials. Recent time-resolved experiments and molecular dynamics simulations suggest that discrete scale invariance may also apply to hierarchical dynamics in proteins, where several fast local conformational changes are a prerequisite for a slow global transition to occur. Employing entropy-based timescale analysis and Markov state modeling to a simple one-dimensional hierarchical model and biomolecular simulation data, it is found that hierarchical systems quite generally give rise to logarithmically spaced discrete timescales. By introducing a one-dimensional reaction coordinate that collectively accounts for the hierarchically coupled degrees of freedom, the free energy landscape exhibits a characteristic staircase shape with two metastable end states, which causes the log-periodic time evolution of the system. The period of the log-oscillations reflects the effective roughness of the energy landscape and can, in simple cases, be interpreted in terms of the barriers of the staircase landscape.

Biomolecular infrared spectroscopy: making time for dynamics

Time resolved infrared spectroscopy of biological molecules has provided a wealth of information relating to structural dynamics, conformational changes, solvation and intermolecular interactions. Challenges still exist however arising from the wide range of timescales over which biological processes occur, stretching from picoseconds to minutes or hours. Experimental methods are often limited by vibrational lifetimes of probe groups, which are typically on the order of picoseconds, while measuring an evolving system continuously over some 18 orders of magnitude in time presents a raft of technological hurdles. In this Perspective, a series of recent advances which allow biological molecules and processes to be studied over an increasing range of timescales, while maintaining ultrafast time resolution, will be reviewed, showing that the potential for real-time observation of biomolecular function draws ever closer, while offering a new set of challenges to be overcome.

Distinguishing between concerted, sequential and barrierless conformational changes: Folding versus allostery

Characterization of transition and intermediate states of reactions provides insights into their mechanisms and is often achieved through analysis of linear free energy relationships. Such an approach has been used extensively in protein folding studies but less so for analyzing allosteric transitions. Here, we point out analogies in ways to characterize pathways and intermediates in folding and allosteric transitions. Achieving an understanding of the mechanisms by which proteins undergo allosteric switching is important in many cases for obtaining insights into how they function.

Versatile Femtosecond Laser Synchronization for Multiple-Timescale Transient Infrared Spectroscopy

Several ways to electronically synchronize different types of amplified femtosecond laser systems are presented based on a single freely programmable electronics hardware: arbitrary-detuning asynchronous optical sampling (ADASOPS), as well as actively locking two femtosecond laser oscillators, albeit not necessarily to the same round-trip frequency. They allow us to rapidly probe a very wide range of timescales, from picoseconds to potentially seconds, in a single transient absorption experiment without the need to move any delay stage. Experiments become possible that address a largely unexplored aspect of many photochemical reactions, in particular in the context of photo-catalysis as well as photoactive proteins, where an initial femtosecond trigger very often initiates a long-lasting cascade of follow-up processes. The approach is very versatile and allows us to synchronize very different lasers, such as a Ti:Sa amplifier and a 100 kHz Yb-laser system. The jitter of the synchronization, and therewith the time-resolution in the transient experiment, lies in the range from 1 to 3 ps, depending on the method. For illustration, transient IR measurements of the excited state solvation and decay of a metal carbonyl complex as well as the full reaction cycle of bacteriorhodopsin are shown. The pros and cons of the various methods are discussed, with regard to the scientific question one might want to address, and also with regard to the laser systems that might be already existent in a laser lab.

MCL-1 promiscuity and the structural resilience of its binding partners

The allosteric protein MCL-1 and its natural inhibitors, the BH3-only proteins PUMA, BIM, and NOXA regulate apoptosis by interacting promiscuously within an entangled binding network. Little is known about the transient processes and dynamic conformational fluctuations that are the basis for the formation and stability of the MCL-1/BH3-only complex. In this study, we designed photoswitchable versions of MCL-1/PUMA and MCL-1/NOXA, and investigated the protein response after an ultrafast photo-perturbation with transient infrared spectroscopy. We observed partial α-helical unfolding in all cases, albeit on strongly varying timescales (1.6 ns for PUMA, 9.7 ns for the previously studied BIM, and 85 ns for NOXA). These differences are interpreted as a BH3-only-specific "structural resilience" to defy the perturbation while remaining in MCL-1's binding pocket. Thus, the presented insights could help to better understand the differences between PUMA, BIM, and NOXA, the promiscuity of MCL-1, in general, and the role of the proteins in the apoptotic network.

Dynamics and Allosteric Information Pathways of Unphosphorylated c-Cbl

Human c-Cbl is a RING-type ligase and plays a central role in the protein degradation cascade. To elucidate its conformational changes related to substrate binding, we performed molecular dynamics simulations of different variants/states of c-Cbl for a cumulative time of 68 μs. Our simulations demonstrate that before the substrate binds, the RING domain samples a broad set of conformational states at a biologically relevant salt concentration, including the closed, partially open, and fully open states, whereas substrate binding leads to a restricted conformational sampling. Phe378 and the C-terminal region play an essential role in stabilizing the partially open state. To visualize the allosteric signal transmission pathways from the substrate-binding site to the 40 Å apart RING domain and identify the critical residues for allostery, we have created a subgraph from the optimal and suboptimal paths. Redundant paths are seen in the SH2 domain where the substrate binds, while the major bottlenecks are found at the junction between the SH2 domain and the linker helix region as well as that between the SH2 domain and the 4H bundle. These bottlenecks separate the paths into two overall routes. The nodes/residues at the bottlenecks on the subgraph are considered allosteric hot spots. This subgraph approach provides a general tool for network visualization and determination of critical residues for allostery. The structurally and allosterically critical residues identified in our work are testable and would provide valuable insights into the emerging strategies for drug discovery, such as targeted protein degradation.

Thermal tuning of protein hydration in a hyperthermophilic enzyme

Water at the protein surface is an active biological molecule that plays a critical role in many functional processes. Using NMR-restrained MD simulations, we here addressed how protein hydration is tuned at high biological temperatures by analysing homologous acylphosphatase enzymes (AcP) possessing similar structure and dynamics under very different thermal conditions. We found that the hyperthermophilic Sso AcP at 80°C interacts with a lower number of structured waters in the first hydration shell than its human homologous mt AcP at 37°C. Overall, the structural and dynamical properties of waters at the surface of the two enzymes resulted similar in the first hydration shell, including solvent molecules residing in the active site. By contrast the dynamical content of water molecules in the second hydration shell was found to diverge, with higher mobility observed in Sso AcP at 80°C. Taken together the results delineate the subtle differences in the hydration properties of mt AcP and Sso AcP, and indicate that the concept of corresponding states with equivalent dynamics in homologous mesophilic and hyperthermophylic proteins should be extended to the first hydration shell.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β₂ adrenergic receptor (β₂AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β₂AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β₂AR, it seems not particularly suited to explain allosteric phenomena.

Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication

While allostery is of paramount importance for protein signaling and regulation, the underlying dynamical process of allosteric communication is not well understood. The PDZ3 domain represents a prime example of an allosteric single-domain protein, as it features a well-established long-range coupling between the C-terminal α₃-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α₃-helix to initiate a conformational change of PDZ3 that propagates from the C-terminus to the bound ligand within 200 ns. Performing extensive nonequilibrium molecular dynamics simulations, the modeling of the experiment reproduces the measured time scales and reveals a detailed picture of the allosteric communication in PDZ3. In particular, a correlation analysis identifies a network of contacts connecting the α₃-helix and the core of the protein, which move in a concerted manner. Representing a one-step process and involving direct α₃-ligand contacts, this cooperative transition is considered as the elementary step in the propagation of conformational change.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

Information Flow and Allosteric Communication in Proteins

Based on Schreiber's work on transfer entropy, a molecular theory of nonlinear information transfer in proteins is developed. The joint distribution function for residue fluctuations is expressed in terms of tensor Hermite polynomials which conveniently separate harmonic and nonlinear contributions to information transfer. The harmonic part of information transfer is expressed as the difference between time dependent and independent mutual information. Third order nonlinearities are discussed in detail. Amount and speed of information transfer between residues, important for understanding allosteric activity in proteins, are discussed. While mutual information shows the maximum amount of information that may be transferred between two residues, it does not explain the actual amount of transfer nor the transfer rate of information. For this, dynamic equations of the system are needed. The solution of the Langevin equation and molecular dynamics trajectories are used in the present work for this purpose. Allosteric communication in Human NAD-dependent isocitrate dehydrogenase is studied as an example. Calculations show that several paths contribute collectively to information transfer. Important residues on these paths are identified. Time resolved information transfer between these residues, their amplitudes and transfer rates, which are in agreement with time resolved ultraviolet resonance Raman measurements in general, are estimated. Estimated transfer rates are in the order of 1-20 megabits per second. Information transfer from third order contributions are one to two orders of magnitude smaller than the harmonic terms, showing that harmonic analysis is a good approximation to information transfer.

Using azobenzene photocontrol to set proteins in motion

Controlling the activity of proteins with azobenzene photoswitches is a potent tool for manipulating their biological function. With the help of light, it is possible to change binding affinities, control allostery or manipulate complex biological processes, for example. Additionally, owing to their intrinsically fast photoisomerization, azobenzene photoswitches can serve as triggers that initiate out-of-equilibrium processes. Such switching of the activity initiates a cascade of conformational events that can be accessed with time-resolved methods. In this Review, we show how the potency of azobenzene photoswitching can be combined with transient spectroscopic techniques to disclose the order of events and experimentally observe biomolecular interactions in real time. This strategy will further our understanding of how a protein can accommodate, adapt and readjust its structure to answer an incoming signal, revealing more of the dynamical character of proteins.

Combining steady state and temperature jump IR spectroscopy to investigate the allosteric effects of ligand binding to dsDNA

Changes in the structural dynamics of double stranded (ds)DNA upon ligand binding have been linked to the mechanism of allostery without conformational change, but direct experimental evidence remains elusive. To address this, a combination of steady state infrared (IR) absorption spectroscopy and ultrafast temperature jump IR absorption measurements has been used to quantify the extent of fast (∼100 ns) fluctuations in (ds)DNA·Hoechst 33258 complexes at a range of temperatures. Exploiting the direct link between vibrational band intensities and base stacking shows that the absolute magnitude of the change in absorbance caused by fast structural fluctuations following the temperature jump is only weakly dependent on the starting temperature of the sample. The observed fast dynamics are some two orders of magnitude faster than strand separation and associated with all points along the 10-base pair duplex d(GCATATATCC). Binding the Hoechst 33258 ligand causes a small but consistent reduction in the extent of these fast fluctuations of base pairs located outside of the ligand binding region. These observations point to a ligand-induced reduction in the flexibility of the dsDNA near the binding site, consistent with an estimated allosteric propagation length of 15 Å, about 5 base pairs, which agrees well with both molecular simulation and coarse-grained statistical mechanics models of allostery leading to cooperative ligand binding.

Sequence of Events during Peptide Unbinding from RNase S: A Complete Experimental Description

The phototriggered unbinding of the intrinsically disordered S-peptide from the RNase S complex is studied with the help of transient IR spectroscopy, covering a wide range of time scales from 100 ps to 10 ms. To that end, an azobenzene moiety has been linked to the S-peptide in a way that its helicity is disrupted by light, thereby initiating its complete unbinding. The full sequence of events is observed, starting from unfolding of the helical structure of the S-peptide on a 20 ns time scale while still being in the binding pocket of the S-protein, S-peptide unbinding after 300 μs, and the structural response of the S-protein after 3 ms. With regard to the S-peptide dynamics, the binding mechanism can be classified as an induced fit, while the structural response of the S-protein is better described as conformational selection.

Intrinsic Dynamics of Protein-Peptide Unbinding

The dynamics of peptide-protein binding and unbinding of a variant of the RNase S system has been investigated. To initiate the process, a photoswitchable azobenzene moiety has been covalently linked to the S-peptide, thereby switching its binding affinity to the S-protein. Transient fluorescence quenching was measured with the help of a time-resolved fluorometer, which has been specifically designed for these experiments and is based on inexpensive light-emitting diodes and laser diodes only. One mutant shows on-off behavior with no specific binding detectable in one of the states of the photoswitch. Unbinding is faster by at least 2 orders of magnitude, compared to that of other variants of the RNase S system. We conclude that unbinding is essentially barrier-less in that case, revealing the intrinsic dynamics of the unbinding event, which occurs on a time scale of a few hundred microseconds in a strongly stretched-exponential manner.

Dynamics underlying hydroxylation selectivity of cytochrome P450cam

Structural heterogeneity and the dynamics of the complexes of enzymes with substrates can determine the selectivity of catalysis; however, fully characterizing how remains challenging as heterogeneity and dynamics can vary at the spatial level of an amino acid residue and involve rapid timescales. We demonstrate the nascent approach of site-specific two-dimensional infrared (IR) spectroscopy to investigate the archetypical cytochrome P450, P450cam, to better delineate the mechanism of the lower regioselectivity of hydroxylation of the substrate norcamphor in comparison to the native substrate camphor. Specific locations are targeted throughout the enzyme by selectively introducing cyano groups that have frequencies in a spectrally isolated region of the protein IR spectrum as local vibrational probes. Linear and two-dimensional IR spectroscopy were applied to measure the heterogeneity and dynamics at each probe and investigate how they differentiate camphor and norcamphor recognition. The IR data indicate that the norcamphor complex does not fully induce a large-scale conformational change to a closed state of the enzyme adopted in the camphor complex. Additionally, a probe directed at the bound substrate experiences rapidly interconverting states in the norcamphor complex that explain the hydroxylation product distribution. Altogether, the study reveals large- and small-scale structural heterogeneity and dynamics that could contribute to selectivity of a cytochrome P450 and illustrates the approach of site-selective IR spectroscopy to elucidate protein dynamics.

Hierarchical dynamics in allostery following ATP hydrolysis monitored by single molecule FRET measurements and MD simulations

We report on a study that combines advanced fluorescence methods with molecular dynamics (MD) simulations to cover timescales from nanoseconds to milliseconds for a large protein. This allows us to delineate how ATP hydrolysis in a protein causes allosteric changes at a distant protein binding site, using the chaperone Hsp90 as test system. The allosteric process occurs via hierarchical dynamics involving timescales from nano- to milliseconds and length scales from Ångstroms to several nanometers. We find that hydrolysis of one ATP is coupled to a conformational change of Arg380, which in turn passes structural information via the large M-domain α-helix to the whole protein. The resulting structural asymmetry in Hsp90 leads to the collapse of a central folding substrate binding site, causing the formation of a novel collapsed state (closed state B) that we characterise structurally. We presume that similar hierarchical mechanisms are fundamental for information transfer induced by ATP hydrolysis through many other proteins.

Sensing the allosteric force

Allosteric regulation is an innate control in most metabolic and signalling cascades that enables living organisms to adapt to the changing environment by tuning the affinity and regulating the activity of target proteins. For a microscopic understanding of this process, a protein system has been designed in such a way that allosteric communication between the binding and allosteric site can be observed in both directions. To that end, an azobenzene-derived photoswitch has been linked to the α3-helix of the PDZ3 domain, arguably the smallest allosteric protein with a clearly identifiable binding and allosteric site. Photo-induced trans-to-cis isomerisation of the photoswitch increases the binding affinity of a small peptide ligand to the protein up to 120-fold, depending on temperature. At the same time, ligand binding speeds up the thermal cis-to-trans back-isomerisation rate of the photoswitch. Based on the energetics of the four states of the system (cis vs trans and ligand-bound vs free), the concept of an allosteric force is introduced, which can be used to drive chemical reactions.

Protonation-Induced Dynamic Allostery in PDZ Domain: Evidence of Perturbation-Independent Universal Response Network

Dynamic allostery is a relatively new paradigm where certain external perturbations may lead to modulation of conformational dynamics at a distant part of a protein without significant changes in the overall structure. While most well-characterized examples of dynamic allostery involve binding with other entities like small molecules, peptides, or nucleic acids, in this work we demonstrate that chemical modifications like protonation may lead to significant dynamical allosteric response in a PDZ domain protein. Tuning the protonation states of two histidine residues (H317 and H372), we identify the allosteric pathways responsible for the dynamic response. Interestingly, the same set of residues that constitute the allosteric response network upon ligand binding seem to be responsible for protonation-induced dynamic allostery. Thus, we propose the existence of an inherent universal response network in signaling proteins, where the same set of residues can respond to varying types of external perturbations in terms of rearrangement of hydrogen-bonded network and redistribution of electrostatic interaction energies.

Real-time observation of ligand-induced allosteric transitions in a PDZ domain

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.

A Time-Dependent Quantum Approach to Allostery and a Comparison With Light-Harvesting in Photosynthetic Phenomenon

The allosteric effect is one of the most important processes in regulating the function of proteins, and the elucidation of this phenomenon plays a significant role in understanding emergent behaviors in biological regulation. In this process, a perturbation, generated by a ligand in a part of the macromolecule (the allosteric site), moves along this system and reaches a specific (active) site, dozens of Ångströms away, with a great efficiency. The dynamics of this perturbation in the macromolecule can model precisely the allosteric process. In this article, we will be studying the general characteristics of allostery, using a time-dependent quantum approach to obtain rules that apply to this kind of process. Considering the perturbation as a wave that moves within the molecular system, we will characterize the allosteric process with three of the properties of this wave in the active site: (1) t_a, the characteristic time for reaching that site, (2) A_a, the amplitude of the wave in this site, and (3) B_a, its corresponding spectral broadening. These three parameters, together with the process mechanism and the perturbation efficiency in the process, can describe the phenomenon. One of the main purposes of this paper is to link the parameters t_a, A_a, and B_a and the perturbation efficiency to the characteristics of the system. There is another fundamental process for life that has some characteristics similar to allostery: the light-harvesting (LH) process in photosynthesis. Here, as in allostery, two distant macromolecular sites are involved—two sites dozens of Ångströms away. In both processes, it is particularly important that the perturbation is distributed efficiently without dissipating in the infinite degrees of freedom within the macromolecule. The importance of considering quantum effects in the LH process is well documented in literature, and the quantum coherences are experimentally proven by time-dependent spectroscopic techniques. Given the existing similarities between these two processes in macromolecules, in this work, we suggest using Quantum Mechanics (QM) to study allostery.

Unfolding Dynamics of a Photoswitchable Helical Peptide

We present an atomistic force field for the azo-moiety of the photoswitchable FK-11-X peptide. We use the parameters to study the unfolding of the peptide through molecular dynamics simulations. The unfolded ensemble contains many different structures, ranging from a partially unfolded peptide to a fully unfolded structure. The averaged computed far-ultraviolet circular dichroism (CD) spectrum of the set of structures, which was simulated using the newly developed force field, agrees well with experiment. The rate of the simulated unfolding process was estimated to have a time constant of 5.80 ± 0.03 ns from the time evolution of the CD spectrum of the peptide, computed from the backbone conformations sampled over 40 simulated trajectories. Our estimated time constant is faster than, but not inconsistent with, previous experimental estimates from time-resolved infrared and optical rotatory dispersion spectroscopy.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Unexpected large impact of small charges on surface frictions with similar wetting properties

Generally, the interface friction on solid surfaces is regarded as consistent with wetting behaviors, characterized by the contact angles. Here using molecular dynamics simulations, we find that even a small charge difference (≤0.36 e) causes a change in the friction coefficient of over an order of magnitude on two-dimensional material and lipid surfaces, despite similar contact angles. This large difference is confirmed by experimentally measuring interfacial friction of graphite and MoS₂ contacting on water, using atomic force microscopy. The large variation in the friction coefficient is attributed to the different fluctuations of localized potential energy under inhomogeneous charge distribution. Our results help to understand the dynamics of two-dimensional materials and biomolecules, generally formed by atoms with small charge, including nanomaterials, such as nitrogen-doped graphene, hydrogen-terminated graphene, or MoS₂, and molecular transport through cell membranes.

Single-molecule FRET methods to study the dynamics of proteins at work

Feynman commented that "Everything that living things do can be understood in terms of the jiggling and wiggling of atoms". Proteins can jiggle and wiggle large structural elements such as domains and subunits as part of their functional cycles. Single-molecule fluorescence resonance energy transfer (smFRET) is an excellent tool to study conformational dynamics and decipher coordinated large-scale motions within proteins. smFRET methods introduced in recent years are geared toward understanding the time scales and amplitudes of function-related motions. This review discusses the methodology for obtaining and analyzing smFRET temporal trajectories that provide direct dynamic information on transitions between conformational states. It also introduces correlation methods that are useful for characterizing intramolecular motions. This arsenal of techniques has been used to study multiple molecular systems, from membrane proteins through molecular chaperones, and we examine some of these studies here. Recent exciting methodological novelties permit revealing very fast, submillisecond dynamics, whose relevance to protein function is yet to be fully grasped.

Modulation of Structural Heterogeneity Controls Phytochrome Photoswitching

Phytochromes sense red/far-red light and control many biological processes in plants, fungi, and bacteria. Although the crystal structures of dark- and light-adapted states have been determined, the molecular mechanisms underlying photoactivation remain elusive. Here, we demonstrate that the conserved tongue region of the PHY domain of a 57-kDa photosensory module of Deinococcus radiodurans phytochrome changes from a structurally heterogeneous dark state to an ordered, light-activated state. The results were obtained in solution by utilizing a laser-triggered activation approach detected on the atomic level with high-resolution protein NMR spectroscopy. The data suggest that photosignaling of phytochromes relies on careful modulation of structural heterogeneity of the PHY tongue.

Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study

To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Functional Role of Solvent Entropy and Conformational Entropy of Metal Binding in a Dynamically Driven Allosteric System

Allostery is a regulatory phenomenon whereby ligand binding to one site influences the binding of the same or a different ligand to another site on a macromolecule. The physical origins of allosteric regulation remain under intense investigation. In general terms, ligand-induced structural changes, perturbations of residue-specific dynamics, and surrounding solvent molecules all potentially contribute to the global energetics of allostery. While the role of solvent is generally well understood in regulatory events associated with major protein structural rearrangements, the degree to which protein dynamics impact solvent degrees of freedom is unclear, particularly in cases of dynamically driven allostery. With the aid of new crystal structures, extensive calorimetric and residue-specific dynamics studies over a range of time scales and temperatures, we dissect for the first time the relative degree to which changes in solvent entropy and residue-specific dynamics impact dynamically driven, allosteric inhibition of DNA binding by Zn in the zinc efflux repressor, CzrA (chromosomal zinc-regulated repressor). We show that non-native residue-specific dynamics in allosterically impaired CzrA mutants are accompanied by significant perturbations in solvent entropy that cannot be predicted from crystal structures. We conclude that functional dynamics are not necessarily restricted to protein residues but involve surface water molecules that may be responding to ligand (Zn)-mediated perturbations in protein internal motions that define the conformational ensemble, rather than major structural rearrangements.

A non-equilibrium approach to allosteric communication

While the theory of protein folding is well developed, including concepts such as rugged energy landscape, folding funnel, etc., the same degree of understanding has not been reached for the description of the dynamics of allosteric transitions in proteins. This is not only due to the small size of the structural change upon ligand binding to an allosteric site, but also due to challenges in designing experiments that directly observe such an allosteric transition. On the basis of recent pump-probe-type experiments (Buchli et al. 2013 Proc. Natl Acad. Sci. USA110, 11 725-11 730. (doi:10.1073/pnas.1306323110)) and non-equilibrium molecular dynamics simulations (Buchenberg et al. 2017 Proc. Natl Acad. Sci. USA114, E6804-E6811. (doi:10.1073/pnas.1707694114)) studying an photoswitchable PDZ2 domain as model for an allosteric transition, we outline in this perspective how such a description of allosteric communication might look. That is, calculating the dynamical content of both experiment and simulation (which agree remarkably well with each other), we find that allosteric communication shares some properties with downhill folding, except that it is an 'order-order' transition. Discussing the multiscale and hierarchical features of the dynamics, the validity of linear response theory as well as the meaning of 'allosteric pathways', we conclude that non-equilibrium experiments and simulations are a promising way to study dynamical aspects of allostery.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

Azobenzene photocontrol of peptides and proteins

The last few years have witnessed significant advances in the use of light as a stimulus to control biomolecular interactions. Great efforts have been devoted to the development of genetically encoded optobiological and small photochromic switches. Newly discovered small molecules now allow researchers to build molecular systems that are sensitive to a wider range of wavelengths of light than ever before with improved switching fidelities and increased lifetimes of the photoactivated states. Because these molecules are relatively small and adopt predictable conformations they are well suited as tools to interrogate cellular function in a spatially and temporally contolled fashion and for applications in photopharmacology.

Azobenzene-based small molecular photoswitches for protein modulation

This review highlights the design strategies of azobenzene photoswitches as well as their applications in the manipulation of biological systems.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.