Search

Adopting a $300 \, \mu$s-long molecular dynamics (MD) trajectory of the reversible folding of villin headpiece (HP35) published by D. E. Shaw Research, we recently constructed a Markov state model (MSM) of the folding process based on interresidue contacts [J. Chem. Theory Comput. 2023, ${\bf {19}}$, 3391]. The model reproduces the MD folding times of the system and predicts that both the native basin and the unfolded region of the free energy landscape are partitioned into several metastable substates that are structurally well characterized. Recognizing the need to establish well-defined but nontrivial benchmark problems, in this Perspective we study to what extent and in what sense this MSM may be employed as a reference model. To this end, we test the robustness of the MSM by comparing it to models that use alternative combinations of features, dimensionality reduction methods and clustering schemes. The study suggests some main characteristics of the folding of HP35, which should be reproduced by any other competitive model of the system. Moreover, the discussion reveals which parts of the MSM workflow matter most for the considered problem, and illustrates the promises and possible pitfalls of state-based models for the interpretation of biomolecular simulations., Comment: arXiv admin note: text overlap with arXiv:2303.03814

Employing extensive quantum-chemical calculations at the DFT/B3LYP and MP2 level, a quartic force field of isolated N-methylacetamide is constructed. Taking into account 24 vibrational degrees of freedom, the model is employed to perform numerically exact vibrational configuration interaction calculations of the vibrational energy relaxation of the amide I mode. It is found that the energy transfer pathways may sensitively depend on details of the theoretical description. Moreover, the exact reference calculations were used to study the applicability and accuracy of (i) the quasiclassical trajectory method, (ii) time-dependent second-order perturbation theory, and (iii) the instantaneous normal mode description of frequency fluctuations. Based on the results, several strategies to describe vibrational energy relaxation in biomolecular systems are discussed., Comment: 18 pages, 6 figures, submitted to J. Phys. Chem. A

Across many fields, a problem of interest is to predict the transition rates between nodes of a network, given limited stationary state and dynamical information. We give a solution using the principle of Maximum Caliber. We find the transition rate matrix by maximizing the path entropy of a random walker on the network constrained to reproducing a stationary distribution and a few dynamical averages. A main finding here is that when constrained only by the mean jump rate, the rate matrix is given by a square-root dependence of the rate, $\omega_{ab} \propto \sqrt{p_b/p_a}$, on $p_a$ and $p_b$, the stationary state populations at nodes a and b. We give two examples of our approach. First, we show that this method correctly predicts the correlated rates in a biochemical network of two genes, where we know the exact results from prior simulation. Second, we show that it correctly predicts rates of peptide conformational transitions, when compared to molecular dynamics simulations. This method can be used to infer large numbers of rates on known networks where smaller numbers of steady-state node populations are known.

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 is 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 $\sim$1~ns to $\sim$10~$\mu$s, can be rationalize by a remodelling of its rugged free energy landscape, with ver 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., Comment: This unedited earlier version of the manuscript may be downloaded for personal use only. Any other use requires prior permission of the author and the National Academy of Sciences USA. The final manuscript was published in Proceedings of the National Academy of Sciences USA 117, 26031-26039 (2020) and can be found under https://www.pnas.org/content/117/42/26031.short

An azobenzene-derived photoswitch has been covalently cross-linked to two sites of the S-peptide 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.

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 macromo-lecular 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.

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. USA 110 , 11 725–11 730. ( doi:10.1073/pnas.1306323110 )) and non-equilibrium molecular dynamics simulations (Buchenberg et al. 2017 Proc. Natl Acad. Sci. USA 114 , 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’.

A local perturbation of a protein may lead to functional changes at some distal site. An example is the PDZ2 domain of human tyrosine phosphatase 1E, which shows an allosteric transition upon binding to a peptide ligand. Recently Buchli et al. presented a time-resolved study of this transition by covalently linking an azobenzene photoswitch across the binding groove and using a femtosecond laser pulse that triggers the cis-trans photoisomerization of azobenzene. To aid the interpretation of these experiments, in this work seven microsecond runs of all-atom molecular dynamics simulations each for the wild-type PDZ2 in the ligand-bound and -free state, as well as the photoswitchable protein (PDZ2S) in the cis and trans states of the photoswitch, in explicit water were conducted. First the theoretical model is validated by recalculating the available NMR data from the simulations. By comparing the results for PDZ2 and PDZ2S, it is analyzed to what extent the photoswitch indeed mimics the free-bound transition. A detailed description of the conformational rearrangement following the cis-trans photoisomerization of PDZ2S reveals a series of photoinduced structural changes that propagate from the anchor residues of the photoswitch via intermediate secondary structure segments to the C-terminus of PDZ2S. The changes of the conformational distribution of the C-terminal region is considered as the distal response of the isolated allosteric protein.