The "violin model": Looking at community networks for dynamic allostery

Role of the αC-β4 loop in protein kinase structure and dynamics

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

Illuminating Protein Allostery by Chemically Accurate Contact Response Analysis (ChACRA)

Decoding allostery at the atomic level is essential for understanding the relationship between a protein's sequence, structure, and dynamics. Recently, we have shown that decomposing temperature responses of inter-residue contacts can reveal allosteric couplings and provide useful insight into the functional dynamics of proteins. The details of this Chemically Accurate Contact Response Analysis (ChACRA) are presented here along with its application to two well-known allosteric proteins. The first protein, IGPS, is a model of ensemble allostery that lacks clear structural differences between the active and inactive states. We show that the application of ChACRA reveals the experimentally identified allosteric coupling between effector and active sites of IGPS. The second protein, ATCase, is a classic example of allostery with distinct active and inactive structural states. Using ChACRA, we directly identify the most significant residue level interactions underlying the enzyme's cooperative behavior. Both test cases demonstrate the utility of ChACRA's unsupervised machine learning approach for dissecting allostery at the residue level.

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus

Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85β subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85β with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85β. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85β-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.

A Fresh Look at the Normal Mode Analysis of Proteins: Introducing Allosteric Co-Vibrational Modes

We propose a new way of utilizing normal modes to study protein conformational transitions. Instead of considering individual modes independently, we show that a weighted mixture of low-frequency vibrational modes can reveal dynamic information about the conformational mechanism in more detail than any single mode can. The weights in the mixed mode, termed the allosteric covibrational mode, are determined using a simple model where the conformational transition is viewed as a perturbation of the coupled harmonic oscillator associated with either of the two conformations. We demonstrate our theory in a biologically relevant example of high pharmaceutical interest involving the V617F mutation of Janus 2 tyrosine kinase (JAK2).

Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution

SRC kinase associated phosphoprotein 1 (SKAP1), an adaptor for protein assembly, plays an important role in the immune system such as stabilizing immune synapses. Understanding how these functions are controlled at the level of the protein-protein interactions is necessary to describe these processes and to develop therapeutics. Here, we dissected the SKAP1 modular organization to recognize SRC kinases and compared it to that of its paralog SRC kinase associated phosphoprotein 2 (SKAP2). Different conserved motifs common to either both proteins or specific to SKAP2 were found using this comparison. Two modules harboring different binding properties between SKAP1 and SKAP2 were identified: one composed of two conserved motifs located in the second interdomain interacting at least with the SH2 domain of SRC kinases and a second one composed of the DIM domain modulated by the SH3 domain and the activation of SRC kinases. This work suggests a convergent evolution of the binding properties of some SRC kinases interacting specifically with either SKAP1 or SKAP2.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study, we employ this paradigm to answer a basic question: in enzyme superfamilies, where the catalytic mechanism, active sites, and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as representatives of the conserved protein tyrosine phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of the catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

Protein dynamics underlying allosteric regulation

Allostery is the mechanism by which information and control are propagated in biomolecules. It regulates ligand binding, chemical reactions, and conformational changes. An increasing level of experimental resolution and control over allosteric mechanisms promises a deeper understanding of the molecular basis for life and powerful new therapeutics. In this review, we survey the literature for an up-to-date biological and theoretical understanding of protein allostery. By delineating five ways in which the energy landscape or the kinetics of a system may change to give rise to allostery, we aim to help the reader grasp its physical origins. To illustrate this framework, we examine three systems that display these forms of allostery: allosteric inhibitors of beta-lactamases, thermosensation of TRP channels, and the role of kinetic allostery in the function of kinases. Finally, we summarize the growing power of computational tools available to investigate the different forms of allostery presented in this review.

Drive-specific selection in multistable mechanical networks

Systems with many stable configurations abound in nature, both in living and inanimate matter, encoding a rich variety of behaviors. In equilibrium, a multistable system is more likely to be found in configurations with lower energy, but the presence of an external drive can alter the relative stability of different configurations in unexpected ways. Living systems are examples par excellence of metastable nonequilibrium attractors whose structure and stability are highly dependent on the specific form and pattern of the energy flow sustaining them. Taking this distinctively lifelike behavior as inspiration, we sought to investigate the more general physical phenomenon of drive-specific selection in nonequilibrium dynamics. To do so, we numerically studied driven disordered mechanical networks of bistable springs possessing a vast number of stable configurations arising from the two stable rest lengths of each spring, thereby capturing the essential physical properties of a broad class of multistable systems. We found that there exists a range of forcing amplitudes for which the attractor states of driven disordered multistable mechanical networks are fine-tuned with respect to the pattern of external forcing to have low energy absorption from it. Additionally, we found that these drive-specific attractor states are further stabilized by precise matching between the multidimensional shape of their orbit and that of the potential energy well they inhabit. Lastly, we showed evidence of drive-specific selection in an experimental system and proposed a general method to estimate the range of drive amplitudes for drive-specific selection.

Protein Kinase Structure and Dynamics: Role of the αC-β4 Loop

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, Local Spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.