Molecular Determinants Underlying Binding Specificities of the ABL Kinase Inhibitors: Combining Alanine Scanning of Binding Hot Spots with Network Analysis of Residue Interactions and Coevolution

Quantifying binding specificity and drug resistance of protein kinase inhibitors is of fundamental importance and remains highly challenging due to complex interplay of structural and thermodynamic factors. In this work, molecular simulations and computational alanine scanning are combined with the network-based approaches to characterize molecular determinants underlying binding specificities of the ABL kinase inhibitors. The proposed theoretical framework unveiled a relationship between ligand binding and inhibitor-mediated changes in the residue interaction networks. By using topological parameters, we have described the organization of the residue interaction networks and networks of coevolving residues in the ABL kinase structures. This analysis has shown that functionally critical regulatory residues can simultaneously embody strong coevolutionary signal and high network centrality with a propensity to be energetic hot spots for drug binding. We have found that selective (Nilotinib) and promiscuous (Bosutinib, Dasatinib) kinase inhibitors can use their energetic hot spots to differentially modulate stability of the residue interaction networks, thus inhibiting or promoting conformational equilibrium between inactive and active states. According to our results, Nilotinib binding may induce a significant network-bridging effect and enhance centrality of the hot spot residues that stabilize structural environment favored by the specific kinase form. In contrast, Bosutinib and Dasatinib can incur modest changes in the residue interaction network in which ligand binding is primarily coupled only with the identity of the gate-keeper residue. These factors may promote structural adaptability of the active kinase states in binding with these promiscuous inhibitors. Our results have related ligand-induced changes in the residue interaction networks with drug resistance effects, showing that network robustness may be compromised by targeted mutations of key mediating residues. This study has outlined mechanisms by which inhibitor binding could modulate resilience and efficiency of allosteric interactions in the kinase structures, while preserving structural topology required for catalytic activity and regulation.     

Dynamics-based community analysis and perturbation response scanning of allosteric interaction networks in the TRAP1 chaperone structures dissect molecular linkage between conformational asymmetry and sequential ATP hydrolysis

Allosteric interactions of the Hsp90 chaperones with cochaperones and diverse protein clients can often exhibit distinct asymmetric features that determine regulatory mechanisms and cellular functions in many signaling networks. The recent crystal structures of the mitochondrial Hsp90 isoform TRAP1 in complexes with ATP analogs have provided first evidence of significant asymmetry in the closed dimerized state that triggers independent activity of the chaperone protomers, whereby preferential hydrolysis of the buckled protomer is followed by conformational flipping between protomers and hydrolysis of the second protomer. Despite significant insights in structural characterizations of the TRAP1 chaperone, the atomistic details and mechanics of allosteric interactions that couple sequential ATP hydrolysis with asymmetric conformational switching in the TRAP1 protomers remain largely unknown. In this work, we explored atomistic and coarse-grained simulations of the TRAP1 dimer structures in combination with the ensemble-based network modeling and perturbation response scanning of residue interaction networks to probe salient features underlying allosteric signaling mechanism. This study has revealed that key effector sites that orchestrate allosteric interactions occupy the ATP binding region and N-terminal interface of the buckled protomer, whereas the main sensors of allosteric signals that drive functional conformational changes during ATPase cycle are consolidated near the client binding region of the straight protomer, channeling the energy of ATP hydrolysis for client remodeling. The community decomposition analysis of the interaction networks and reconstruction of allosteric communication pathways in the TRAP1 structures have quantified mechanism of allosteric regulation, revealing control points and interactions that coordinate asymmetric switching during ATP hydrolysis.     

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.     

The energy landscape analysis of cancer mutations in protein kinases

The growing interest in quantifying the molecular basis of protein kinase activation and allosteric regulation by cancer mutations has fueled computational studies of allosteric signaling in protein kinases. In the present study, we combined computer simulations and the energy landscape analysis of protein kinases to characterize the interplay between oncogenic mutations and locally frustrated sites as important catalysts of allostetric kinase activation. While structurally rigid kinase core constitutes a minimally frustrated hub of the catalytic domain, locally frustrated residue clusters, whose interaction networks are not energetically optimized, are prone to dynamic modulation and could enable allosteric conformational transitions. The results of this study have shown that the energy landscape effect of oncogenic mutations may be allosteric eliciting global changes in the spatial distribution of highly frustrated residues. We have found that mutation-induced allosteric signaling may involve a dynamic coupling between structurally rigid (minimally frustrated) and plastic (locally frustrated) clusters of residues. The presented study has demonstrated that activation cancer mutations may affect the thermodynamic equilibrium between kinase states by allosterically altering the distribution of locally frustrated sites and increasing the local frustration in the inactive form, while eliminating locally frustrated sites and restoring structural rigidity of the active form. The energy landsape analysis of protein kinases and the proposed role of locally frustrated sites in activation mechanisms may have useful implications for bioinformatics-based screening and detection of functional sites critical for allosteric regulation in complex biomolecular systems.     

Exploring Mechanisms of Allosteric Regulation and Communication Switching in the Multiprotein Regulatory Complexes of the Hsp90 Chaperone with Cochaperones and Client Proteins: Atomistic Insights from Integrative Biophysical Modeling and Network Analysis of Conformational Landscapes

Understanding molecular principles underlying Hsp90 chaperone functions and modulation of client activity is fundamental to dissect activation mechanisms of many proteins. In this work, we performed a computational investigation of the Hsp90-Hsp70-Hop-CR client complex to examine allosteric regulatory mechanisms underlying dynamic chaperone interactions and principles of chaperone-dependent client recognition and remodeling. Conformational dynamics analysis using high-resolution coarse-grained simulations and ensemble-based local frustration analysis suggest that the Hsp90 chaperone could recognize and recruit the GR client by invoking reciprocal dynamic exchanges near the intermolecular interfaces with the client. Using mutational scanning of the intermolecular residues in the Hsp90-Hsp70-Hop-GR complex, we identified binding energy hotspots in the regulatory complex. Perturbation-based network analysis and dynamic fluctuations-based modeling of allosteric residue potentials are employed for a detailed analysis of allosteric interaction networks and identification of conformational communication switches. We found that allosteric interactions between the Hsp90, the client-bound Hsp70 and Hop cochaperone can define two allosteric residue clusters that control client recruitment in which the intrinsic Hsp70 allostery is exploited to mediate integration of the Hsp70-bound client into the Hsp90 chaperone system. The results suggest a model of dynamics-driven allostery that enables efficient client recruitment and loading through allosteric couplings between intermolecular interfaces and communication switch centers. This study showed that the Hsp90 interactions with client proteins may operate under dynamic-based allostery in which ensembles of preexisting conformational states and intrinsic allosteric pathways present in the Hsp90 and Hsp70 chaperones can be exploited for recognition and integration of substrate proteins.     

Networks for the allosteric control of protein kinases

The allosteric regulation of protein kinases serves as an efficient strategy for molecular communication, event coupling and interconversion between catalytic states. Recent co-crystal structures have revealed novel ways in which kinases control activity and substrate specificity following phosphorylation, dimerization, or binding to regulatory proteins, substrates and scaffolds. In addition, hydrogen exchange coupled with mass spectrometry is emerging as a complementary strategy to probe the solution behavior of kinases; recent results have shown that allosteric regulation may involve transitions in protein motions as well as structural rearrangements.     

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.     

Role for SUR2A ED domain in allosteric coupling within the K(ATP) channel complex

Allosteric regulation of heteromultimeric ATP-sensitive potassium (K(ATP)) channels is unique among protein systems as it implies transmission of ligand-induced structural adaptation at the regulatory SUR subunit, a member of ATP-binding cassette ABCC family, to the distinct pore-forming K+ (Kir6.x) channel module. Cooperative interaction between nucleotide binding domains (NBDs) of SUR is a prerequisite for K(ATP) channel gating, yet pathways of allosteric intersubunit communication remain uncertain. Here, we analyzed the role of the ED domain, a stretch of 15 negatively charged aspartate/glutamate amino acid residues (948-962) of the SUR2A isoform, in the regulation of cardiac K(ATP) channels. Disruption of the ED domain impeded cooperative NBDs interaction and interrupted the regulation of K(ATP) channel complexes by MgADP, potassium channel openers, and sulfonylurea drugs. Thus, the ED domain is a structural component of the allosteric pathway within the K(ATP) channel complex integrating transduction of diverse nucleotide-dependent states in the regulatory SUR subunit to the open/closed states of the K+-conducting channel pore.     

Probing molecular mechanisms of the Hsp90 chaperone: biophysical modeling identifies key regulators of functional dynamics

Deciphering functional mechanisms of the Hsp90 chaperone machinery is an important objective in cancer biology aiming to facilitate discovery of targeted anti-cancer therapies. Despite significant advances in understanding structure and function of molecular chaperones, organizing molecular principles that control the relationship between conformational diversity and functional mechanisms of the Hsp90 activity lack a sufficient quantitative characterization. We combined molecular dynamics simulations, principal component analysis, the energy landscape model and structure-functional analysis of Hsp90 regulatory interactions to systematically investigate functional dynamics of the molecular chaperone. This approach has identified a network of conserved regions common to the Hsp90 chaperones that could play a universal role in coordinating functional dynamics, principal collective motions and allosteric signaling of Hsp90. We have found that these functional motifs may be utilized by the molecular chaperone machinery to act collectively as central regulators of Hsp90 dynamics and activity, including the inter-domain communications, control of ATP hydrolysis, and protein client binding. These findings have provided support to a long-standing assertion that allosteric regulation and catalysis may have emerged via common evolutionary routes. The interaction networks regulating functional motions of Hsp90 may be determined by the inherent structural architecture of the molecular chaperone. At the same time, the thermodynamics-based "conformational selection" of functional states is likely to be activated based on the nature of the binding partner. This mechanistic model of Hsp90 dynamics and function is consistent with the notion that allosteric networks orchestrating cooperative protein motions can be formed by evolutionary conserved and sparsely connected residue clusters. Hence, allosteric signaling through a small network of distantly connected residue clusters may be a rather general functional requirement encoded across molecular chaperones. The obtained insights may be useful in guiding discovery of allosteric Hsp90 inhibitors targeting protein interfaces with co-chaperones and protein binding clients.     

Global conformational changes in allosteric proteins. A study of Escherichia coli cAMP receptor protein and muscle pyruvate kinase.

One of the basic features in allosteric regulation involves long range transduction of information. Based on crystallographic data on protein systems that are regulated by allosteric mechanisms, a global conformational change has always been observed. It is, therefore, important and useful to correlate the cooperativity of global structural change with the mode of binding of the regulatory ligand. Two systems were chosen for study, namely Escherichia coli cAMP receptor protein and muscle pyruvate kinase, which show negative and positive cooperativity in the binding of allosteric ligands, respectively. Quantitative titration of the global structural change, monitored by a high precision analytical gel chromatography technique, was conducted as a function of allosteric effector concentration. The results obtained for cAMP receptor protein show that the protein undergoes contraction upon binding of cAMP. The decreases in Stokes radius associated with complex formation are 0.1 +/- 0.1 and 0.7 +/- 0.1 A when one and two cAMP-binding sites are filled, respectively. The results for the pyruvate kinase system show a concerted structural change that quantitatively match the predicted behavior based on equilibrium constants derived from the analysis of steady state kinetic data by a two-state model. Hence, for these two systems, these results show that negative and positive cooperativity are correlated with sequential and concerted modes of structural change, respectively.     

The N-terminal end of the catalytic domain of SRC kinase Hck is a conformational switch implicated in long-range allosteric regulation

Signal transduction in cell growth and proliferation involves regulation of kinases through long-range allostery between remote protein regions. Molecular dynamics free energy calculations are used to clarify the coupling between the catalytic domain of Src kinase Hck and its N-terminal end connecting to the regulatory SH2 and SH3 modules. The N-terminal end is stable in the orientation required for the regulatory modules to remain properly bound only in the inactive catalytic domain. In the active catalytic domain, the N-terminal end prefers a different conformation consistent with dissociation of the regulatory modules. The free energy surface shows that the N-terminal end acts as a reversible two-state conformational switch coupling the catalytic domain to the regulatory modules. Structural analogy with insulin receptor kinase and c-Src suggests that such reversible conformational switching in a critical hinge region could be a common mechanism in long-range allosteric regulation of protein kinase activity.     

Identification of core allosteric sites through temperature- and nucleus-invariant chemical shift covariance

Allosteric regulation is essential to control biological function. In addition, allosteric sites offer a promising venue for selective drug targeting. However, accurate mapping of allosteric sites remains challenging since allostery relies on often subtle, yet functionally relevant, structural and dynamical changes. A viable approach proposed to overcome such challenge is chemical shift covariance analysis (CHESCA). Although CHESCA offers an exhaustive map of allosteric networks, it is critical to define the core allosteric sites to be prioritized in subsequent functional studies or in the design of allosteric drugs. Here, we propose two new CHESCA-based methodologies, called temperature CHESCA (T-CHESCA) and CLASS-CHESCA, aimed at narrowing down allosteric maps to the core allosteric residues. Both T- and CLASS-CHESCAs rely on the invariance of core inter-residue correlations to changes in the chemical shifts of the active and inactive conformations interconverting in fast exchange. In T-CHESCA the chemical shifts of such states are modulated through temperature changes, while in CLASS-CHESCA through variations in the spin-active nuclei involved in pairwise correlations. T- and CLASS-CHESCAs, as well as complete-linkage CHESCA, were applied to the cAMP-binding domain of the exchange protein directly activated by cAMP (EPAC), which serves as a prototypical allosteric switch. Residues consistently identified by the three CHESCA methods were found in previously identified EPAC allosteric core sites. Hence, T-, CLASS-, and CL-CHESCA provide a toolset to establish allosteric site hierarchy and triage allosteric sites to be further analyzed by mutations and functional assays. Furthermore, the core allosteric networks selectively revealed through T- and CLASS-CHESCA are expected to facilitate the mechanistic understanding of disease-related mutations and the design of selective allosteric modulators.     

Structure-Based Network Analysis of Activation Mechanisms in the ErbB Family of Receptor Tyrosine Kinases: The Regulatory Spine Residues Are Global Mediators of Structural Stability and Allosteric Interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.     

Correlating allostery with rigidity

Allosteric proteins demonstrate the phenomenon of a ligand binding to a protein at a regulatory or effector site and thereby changing the chemical affinity of the catalytic site. As such, allostery is extremely important biologically as a regulatory mechanism for molecular concentrations in many cellular processes. One particularly interesting feature of allostery is that often the catalytic and effector sites are separated by a large distance. Structural comparisons of allosteric proteins resolved in both inactive and active states indicate that a variety of structural rearrangement and changes in motions may contribute to general allosteric behavior. In general it is expected that the coupling of catalytic and regulatory sites is responsible for allosteric behavior. We utilize a novel examination of allostery using rigidity analysis of the underlying graph of the protein structures. Our results indicate a general global change in rigidity associated with allosteric transitions where the R state is more rigid than the T state. A set of allosteric proteins with heterotropic interactions is used to test the hypothesis that catalytic and effector sites are structurally coupled. Observation of a rigid path connecting the effector and catalytic sites in 68.75% of the structures points to rigidity as a means by which the distal sites communicate with each other and so contribute to allosteric regulation. Thus structural rigidity is shown to be a fundamental underlying property that promotes cooperativity and non-locality seen in allostery.     

Amino acid biosynthesis: new architectures in allosteric enzymes.

This review focuses on the allosteric controls in the Aspartate-derived and the branched-chain amino acid biosynthetic pathways examined both from kinetic and structural points of view. The objective is to show the differences that exist among the plant and microbial worlds concerning the allosteric regulation of these pathways and to unveil the structural bases of this diversity. Indeed, crystallographic structures of enzymes from these pathways have been determined in bacteria, fungi and plants, providing a wonderful opportunity to obtain insight into the acquisition and modulation of allosteric controls in the course of evolution. This will be examined using two enzymes, threonine synthase and the ACT domain containing enzyme aspartate kinase. In a last part, as many enzymes in these pathways display regulatory domains containing the conserved ACT module, the organization of ACT domains in this kind of allosteric enzymes will be reviewed, providing explanations for the variety of allosteric effectors and type of controls observed.     

Reverse genetics of Escherichia coli glycerol kinase allosteric regulation and glucose control of glycerol utilization in vivo

Reverse genetics is used to evaluate the roles in vivo of allosteric regulation of Escherichia coli glycerol kinase by the glucose-specific phosphocarrier of the phosphoenolpyruvate:glycose phosphotransferase system, IIA(Glc) (formerly known as III(glc)), and by fructose 1,6-bisphosphate. Roles have been postulated for these allosteric effectors in glucose control of both glycerol utilization and expression of the glpK gene. Genetics methods based on homologous recombination are used to place glpK alleles with known specific mutations into the chromosomal context of the glpK gene in three different genetic backgrounds. The alleles encode glycerol kinases with normal catalytic properties and specific alterations of allosteric regulatory properties, as determined by in vitro characterization of the purified enzymes. The E. coli strains with these alleles display the glycerol kinase regulatory phenotypes that are expected on the basis of the in vitro characterizations. Strains with different glpR alleles are used to assess the relationships between allosteric regulation of glycerol kinase and specific repression in glucose control of the expression of the glpK gene. Results of these studies show that glucose control of glycerol utilization and glycerol kinase expression is not affected by the loss of IIA(Glc) inhibition of glycerol kinase. In contrast, fructose 1,6-bisphosphate inhibition of glycerol kinase is the dominant allosteric control mechanism, and glucose is unable to control glycerol utilization in its absence. Specific repression is not required for glucose control of glycerol utilization, and the relative roles of various mechanisms for glucose control (catabolite repression, specific repression, and inducer exclusion) are different for glycerol utilization than for lactose utilization.     

An overview of structure,function, and regulation of pyruvate kinases

In the last step of glycolysis Pyruvate kinase catalyzes the irreversible conversion of ADP and phosphoenolpyruvate to ATP and pyruvic acid, both crucial for cellular metabolism. Thus pyruvate kinase plays a key role in controlling the metabolic flux and ATP production. The hallmark of the activity of different pyruvate kinases is their tight modulation by a variety of mechanisms including the use of a large number of physiological allosteric effectors in addition to their homotropic regulation by phosphoenolpyruvate. Binding of effectors signals precise and orchestrated movements in selected areas of the protein structure that alter the catalytic action of these evolutionarily conserved enzymes with remarkably conserved architecture and sequences. While the diverse nature of the allosteric effectors has been discussed in the literature, the structural basis of their regulatory effects is still not well understood because of the lack of data representing conformations in various activation states. Results of recent studies on pyruvate kinases of different families suggest that members of evolutionarily related families follow somewhat conserved allosteric strategies but evolutionarily distant members adopt different strategies. Here we review the structure and allosteric properties of pyruvate kinases of different families for which structural data are available.     

Distinct conformational dynamics and allosteric networks in alpha tryptophan synthase during active catalysis

Experimental observations of enzymes under active turnover conditions have brought new insight into the role of protein motions and allosteric networks in catalysis. Many of these studies characterize enzymes under dynamic chemical equilibrium conditions, in which the enzyme is actively catalyzing both the forward and reverse reactions during data acquisition. We have previously analyzed conformational dynamics and allosteric networks of the alpha subunit of tryptophan synthase under such conditions using NMR. We have proposed that this working state represents a four to one ratio of the enzyme bound with the indole‐3‐glycerol phosphate substrate (E:IGP) to the enzyme bound with the products indole and glyceraldehyde‐3‐phosphate (E:indole:G3P). Here, we analyze the inactive D60N variant to deconvolute the contributions of the substrate‐ and products‐bound states to the working state. While the D60N substitution itself induces small structural and dynamic changes, the D60N E:IGP and E:indole:G3P states cannot entirely account for the conformational dynamics and allosteric networks present in the working state. The act of chemical bond breakage and/or formation, or possibly the generation of an intermediate, may alter the structure and dynamics present in the working state. As the enzyme transitions from the substrate‐bound to the products‐bound state, millisecond conformational exchange processes are quenched and new allosteric connections are made between the alpha active site and the surface which interfaces with the beta subunit. The structural ordering of the enzyme and these new allosteric connections may be important in coordinating the channeling of the indole product into the beta subunit.