The origin of allosteric functional modulation: multiple pre-existing pathways

Role of Protein Dimeric Interface in Allosteric Inhibition of N-Acetyl-Aspartate Hydrolysis by Human Aspartoacylase

The results of molecular modeling suggest a mechanism of allosteric inhibition upon hydrolysis of N-acetyl-aspartate (NAA), one of the most abundant amino acid derivatives in brain, by human aspartoacylase (hAsp). Details of this reaction are important to suggest the practical ways to control the enzyme activity. Search for allosteric sites using the Allosite web server and SiteMap analysis allowed us to identify substrate binding pockets located at the interface between the subunits of the hAsp dimer molecule. Molecular docking of NAA to the pointed areas at the dimer interface predicted a specific site, in which the substrate molecule interacts with the Gly237, Arg233, Glu290, and Lys292 residues. Analysis of multiple long-scaled molecular dynamics trajectories (the total simulation time exceeded 1.5 μs) showed that binding of NAA to the identified allosteric site induced significant rigidity to the protein loops with the amino acid side chains forming gates to the enzyme active site. Application of the protein dynamical network algorithms showed that substantial reorganization of the signal propagation pathways of intersubunit communication in the dimer occurred upon allosteric NAA binding to the remote site. The modeling approaches provide an explanation to the observed decrease of the reaction rate of NAA hydrolysis by hAsp at high substrate concentrations.

How Intrinsic Dynamics Mediates the Allosteric Mechanism in the ABC Transporter Nucleotide Binding Domain Dimer

A protein's architecture facilitates specific motions-intrinsic dynamic modes-that are employed to effect function. Here we used molecular dynamics (MD) simulations to investigate the dynamics of the MJ0796 ABC transporter nucleotide-binding domain (NBD). ABC transporter NBDs form a rotationally symmetric dimer whereby two equivalent active sites are formed at their interface; in complex with a dimer of transmembrane domains they hydrolyze ATP to energize translocation of substrates across cellular membranes. Our data suggest the ABC NBD's ensemble of functional states can be understood predominately in terms of conformational changes between its major subdomains, occurring along two orthogonal dynamic modes. The data show that ligands and oligomeric interactions modulate the equilibrium conformation of the NBD with respect to these motions, suggesting that allostery is achieved by affecting the energetic profile along these two modes. The observed dynamics and allostery integrate consonantly and logically within a mechanistic framework for the ABC NBD dimer, which is supported by a large body of experimental and theoretical data, providing a higher resolution view of the enzyme's dynamic cycle. Our study shows how valuable mechanistic inferences can be derived from accessible short-time scale MD simulations of an enzyme's substructures.

Molecular mechanisms of cooperative binding of transcription factors Runx1-CBFβ-Ets1 on the TCRα gene enhancer

Ets1 is an essential transcription factor (TF) for several important physiological processes, including cell proliferation and differentiation. Its recognition of the enhancer region of the TCRα gene is enhanced by the cooperative binding of the Runx1–CBFβ heterodimer, with the cancelation of phosphorylation-dependent autoinhibition. The detailed mechanism of this interesting cooperativity between Ets1 and the Runx1–CBFβ heterodimer is still largely unclear. Here, we investigated the molecular mechanisms of this cooperativity, by using molecular dynamics simulations. Consequently, we detected high flexibility of the loop region between the HI2 and H1 helices of Ets1. Upon Runx1–CBFβ heterodimer binding, this loop transiently adopts various sub-stable conformations in its interactions with the DNA. In addition, a network analysis suggested an allosteric pathway in the molecular assembly and identified some key residues that coincide with previous experimental studies. Our simulations suggest that the cooperative binding of Ets1 and the Runx1–CBFβ heterodimer alters the DNA conformation and induces sub-stable conformations of the HI2–H1 loop of Ets1. This phenomenon increases the flexibility of the regulatory module, including the HI2 helix, and destabilizes the inhibitory form of this module. Thus, we hypothesize that this effect facilitates Ets1–DNA binding and prevents the phosphorylation-dependent DNA binding autoinhibition.

Conformational diversity analysis reveals three functional mechanisms in proteins

Protein motions are a key feature to understand biological function. Recently, a large-scale analysis of protein conformational diversity showed a positively skewed distribution with a peak at 0.5 Å C-alpha root-mean-square-deviation (RMSD). To understand this distribution in terms of structure-function relationships, we studied a well curated and large dataset of ~5,000 proteins with experimentally determined conformational diversity. We searched for global behaviour patterns studying how structure-based features change among the available conformer population for each protein. This procedure allowed us to describe the RMSD distribution in terms of three main protein classes sharing given properties. The largest of these protein subsets (~60%), which we call "rigid" (average RMSD = 0.83 Å), has no disordered regions, shows low conformational diversity, the largest tunnels and smaller and buried cavities. The two additional subsets contain disordered regions, but with differential sequence composition and behaviour. Partially disordered proteins have on average 67% of their conformers with disordered regions, average RMSD = 1.1 Å, the highest number of hinges and the longest disordered regions. In contrast, malleable proteins have on average only 25% of disordered conformers and average RMSD = 1.3 Å, flexible cavities affected in size by the presence of disordered regions and show the highest diversity of cognate ligands. Proteins in each set are mostly non-homologous to each other, share no given fold class, nor functional similarity but do share features derived from their conformer population. These shared features could represent conformational mechanisms related with biological functions.

Uncoupling of an ammonia channel as a mechanism of allosteric inhibition in anthranilate synthase of Serratia marcescens: dynamic and graph theoretical analysis

Network modeling and molecular dynamic studies reveal the perturbation in communication pathways as a mechanism of allosteric inhibition in anthranilate synthase.

DNA-binding protects p53 from interactions with cofactors involved in transcription-independent functions

Binding-induced conformational changes of a protein at regions distant from the binding site may play crucial roles in protein function and regulation. The p53 tumour suppressor is an example of such an allosterically regulated protein. Little is known, however, about how DNA binding can affect distal sites for transcription factors. Furthermore, the molecular details of how a local perturbation is transmitted through a protein structure are generally elusive and occur on timescales hard to explore by simulations. Thus, we employed state-of-the-art enhanced sampling atomistic simulations to unveil DNA-induced effects on p53 structure and dynamics that modulate the recruitment of cofactors and the impact of phosphorylation at Ser215. We show that DNA interaction promotes a conformational change in a region 3 nm away from the DNA binding site. Specifically, binding to DNA increases the population of an occluded minor state at this distal site by more than 4-fold, whereas phosphorylation traps the protein in its major state. In the minor conformation, the interface of p53 that binds biological partners related to p53 transcription-independent functions is not accessible. Significantly, our study reveals a mechanism of DNA-mediated protection of p53 from interactions with partners involved in the p53 transcription-independent signalling. This also suggests that conformational dynamics is tightly related to p53 signalling.

Allosteric therapies for lung cancer

Allostery is a regulation at a distance by conveying information from one site to another and an intrinsic property of dynamic proteins. Allostery plays an essential role in receptor trafficking, signal transmission, controlled catalysis, gene turn on/off, or cell apoptosis. Allosteric mutations are considered as one of causes responsible for cancer development, leading to "allosteric diseases" by stabilizing an active or inactive conformation or changing the dynamic distribution of preexisting propagation pathways. The present article mainly focuses on the potential of allosteric therapies for lung cancer. Allosteric drugs may have several advantages over traditional drugs. The epidermal growth factor receptor mutations and signaling pathways downstream (such as PI3K/AKT/mTOR and RAS/RAF/MEK/ERK pathways) were suggested to play a key role in lung cancer and considered as targets of allosteric therapy. Some allosteric inhibitors for lung cancer-specific targets and a series of preclinical trials of allosteric inhibitors for lung cancer have been developed and reported. We expect that allosteric therapies will gain more attentions to develop combinatorial strategies for lung cancer and metastasis.

Long-Range Signaling in MutS and MSH Homologs via Switching of Dynamic Communication Pathways

Allostery is conformation regulation by propagating a signal from one site to another distal site. This study focuses on the long-range communication in DNA mismatch repair proteins MutS and its homologs where intramolecular signaling has to travel over 70 Å to couple lesion detection to ATPase activity and eventual downstream repair. Using dynamic network analysis based on extensive molecular dynamics simulations, multiple preserved communication pathways were identified that would allow such long-range signaling. The pathways appear to depend on the nucleotides bound to the ATPase domain as well as the type of DNA substrate consistent with previously proposed functional cycles of mismatch recognition and repair initiation by MutS and homologs. A mechanism is proposed where pathways are switched without major conformational rearrangements allowing for efficient long-range signaling and allostery.

Oligomeric interface modulation causes misregulation of purine 5´-nucleotidase in relapsed leukemia

Background

Relapsed acute lymphoblastic leukemia (ALL) is one of the main causes of mortality in childhood malignancies. Previous genetic studies demonstrated that chemoresistant ALL is driven by activating mutations in NT5C2, the gene encoding cytosolic 5´-nucleotidase (cN-II). However, molecular mechanisms underlying this hyperactivation are still unknown. Here, we present kinetic and structural properties of cN-II variants that represent 75 % of mutated alleles in patients who experience relapsed ALL (R367Q, R238W and L375F).

Results

Enzyme kinetics measurements revealed that the mutants are consitutively active without need for allosteric activators. This shows that hyperactivity is not caused by a direct catalytic effect but rather by misregulation of cN-II. X-ray crystallography combined with mass spectrometry-based techniques demonstrated that this misregulation is driven by structural modulation of the oligomeric interface within the cN-II homotetrameric assembly. These specific conformational changes are shared between the studied variants, despite the relatively random spatial distribution of the mutations.

Conclusions

These findings define a common molecular mechanism for cN-II hyperactivity, which provides a solid basis for targeted therapy of leukemia. Our study highlights the cN-II oligomerization interface as an attractive pharmacological target.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0313-y) contains supplementary material, which is available to authorized users.

Investigation of a Catenane with a Responsive Noncovalent Network: Mimicking Long-Range Responses in Proteins

We report a functional synthetic model for studying the noncovalent networks (NCNs) required for complex protein functions. The model [2]-catenane is self-assembled from dipeptide building blocks and contains an extensive network of hydrogen bonds and aromatic interactions. Perturbations to the catenane cause compensating changes in the NCNs structure and dynamics, resulting in long-distance changes reminiscent of a protein. Key findings include the notion that NCNs require regions of negative cooperativity, or "frustrated" noncovalent interactions, as a source of potential energy for driving the response. We refer to this potential energy as latent free energy and describe a mechanistic and energetic model for responsive systems.

Prediction of allosteric sites and mediating interactions through bond-to-bond propensities

Allostery is a fundamental mechanism of biological regulation, in which binding of a molecule at a distant location affects the active site of a protein. Allosteric sites provide targets to fine-tune protein activity, yet we lack computational methodologies to predict them. Here we present an efficient graph-theoretical framework to reveal allosteric interactions (atoms and communication pathways strongly coupled to the active site) without a priori information of their location. Using an atomistic graph with energy-weighted covalent and weak bonds, we define a bond-to-bond propensity quantifying the non-local effect of instantaneous bond fluctuations propagating through the protein. Significant interactions are then identified using quantile regression. We exemplify our method with three biologically important proteins: caspase-1, CheY, and h-Ras, correctly predicting key allosteric interactions, whose significance is additionally confirmed against a reference set of 100 proteins. The almost-linear scaling of our method renders it suitable for high-throughput searches for candidate allosteric sites.

Regulation of integrins in platelets

Blood platelets prevent bleeding after trauma by forming occlusive aggregates at sites of vascular injury. Platelet aggregation is mediated by the integrin heterodimer αIIbβ3 and occurs when platelet agonists generated at the injury site convert αIIbβ3 from its resting to its active conformation. Active αIIbβ3 is then able to bind macromolecular ligands such as fibrinogen that crosslink adjacent platelets into hemostatic aggregates. Platelets circulate in a plasma milieu containing high concentrations of the principal αIIbβ3 ligand fibrinogen. Thus, αIIbβ3 activity is tightly regulated to prevent the spontaneous formation of platelet aggregates. αIIbβ3 activity is regulated at least three levels. First, intramolecular interactions involving motifs located in the membrane-proximal stalk regions, transmembrane domains, and the membrane-proximal cytosolic tails of αIIb and β3 maintain αIIbβ3 in its inactive conformation. Transmembrane domain interactions appear particularly important because disrupting these interactions causes constitutive αIIbβ3 activation. Second, the agonist-stimulated binding of the cytosolic proteins talin and kindlin-3 to the β3 cytosolic tail rapidly causes αIIbβ3 activation by disrupting the intramolecular interactions constraining αIIbβ3 activity. Third, the strength of ligand binding to active αIIbβ3 seems to be allosterically regulated. Thus, αIIbβ3 exists in a minimum of three interconvertible states: an inactive (resting) state that does not interact with ligands and two active ligand binding states that differ in their affinity for fibrinogen and in the mechanical stability of fibrinogen complexes they form.

NMR Characterization of Information Flow and Allosteric Communities in the MAP Kinase p38γ

The intramolecular network structure of a protein provides valuable insights into allosteric sites and communication pathways. However, a straightforward method to comprehensively map and characterize these pathways is not currently available. Here we present an approach to characterize intramolecular network structure using NMR chemical shift perturbations. We apply the method to the mitogen activated protein kinase (MAPK) p38γ. p38γ contains allosteric sites that are conserved among eukaryotic kinases as well as unique to the MAPK family. How these regulatory sites communicate with catalytic residues is not well understood. Using our method, we observe and characterize for the first time information flux between regulatory sites through a conserved kinase infrastructure. This network is accessed, reinforced, and broken in various states of p38γ, reflecting the functional state of the protein. We demonstrate that the approach detects critical junctions in the network corresponding to biologically significant allosteric sites and pathways.

Pyruvate Occupancy in the Carboxyl Transferase Domain of Pyruvate Carboxylase Facilitates Product Release from the Biotin Carboxylase Domain through an Intermolecular Mechanism

Protein structure, ligand binding, and catalytic turnover contributes to the governance of catalytic events occurring at spatially distinct domains in multifunctional enzymes. Coordination of these catalytic events partially rests on the ability of spatially discrete active sites to communicate with other allosteric and active sites on the same polypeptide chain (intramolecular) or on different polypeptide chains (intermolecular) within the holoenzyme. Often, communication results in long-range effects on substrate binding or product release. For example, pyruvate binding to the carboxyl transferase (CT) domain of pyruvate carboxylase (PC) increases the rate of product release in the biotin carboxylase (BC) domain. In order to address how CT domain ligand occupancy is "sensed" by other domains, we generated functional, mixed hybrid tetramers using the E218A (inactive BC domain) and T882S (low pyruvate binding, low activity) mutant forms of PC. The apparent Ka pyruvate for the pyruvate-stimulated release of Pi catalyzed by the T882S:E218A[1:1] hybrid tetramer was comparable to the wild-type enzyme and nearly 10-fold lower than that for the T882S homotetramer. In addition, the ratio of the rates of oxaloacetate formation to Pi release for the WT:T882S[1:1] and E218A:T882S[1:1] hybrid tetramer-catalyzed reactions was 0.5 and 0.6, respectively, while the T882S homotetramer exhibited a near 1:1 coupling of the two domains, suggesting that the mechanisms coordinating catalytic events is more complicated that we initially assumed. The results presented here are consistent with an intermolecular communication mechanism, where pyruvate binding to the CT domain is "sensed" by domains on a different polypeptide chain within the tetramer.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Allosteric sites can be identified based on the residue-residue interaction energy difference

Allosteric drugs act at a distance to regulate protein functions. They have several advantages over conventional orthosteric drugs, including diverse regulation types and fewer side effects. However, the rational design of allosteric ligands remains a challenge, especially when it comes to the identification allosteric binding sites. As the binding of allosteric ligands may induce changes in the pattern of residue-residue interactions, we calculated the residue-residue interaction energies within the allosteric site based on the molecular mechanics generalized Born surface area energy decomposition scheme. Using a dataset of 17 allosteric proteins with structural data for both the apo and the ligand-bound state available, we used conformational ensembles generated by molecular dynamics simulations to compute the differences in the residue-residue interaction energies in known allosteric sites from both states. For all the known sites, distinct interaction energy differences (>25%) were observed. We then used CAVITY, a binding site detection program to identify novel putative allosteric sites in the same proteins. This yielded a total of 31 "druggable binding sites," of which 21 exhibited >25% difference in residue interaction energies, and were hence predicted as novel allosteric sites. Three of the predicted allosteric sites were supported by recent experimental studies. All the predicted sites may serve as novel allosteric sites for allosteric ligand design. Our study provides a computational method for identifying novel allosteric sites for allosteric drug design.

Rigid Residue Scan Simulations Systematically Reveal Residue Entropic Roles in Protein Allostery

Intra-protein information is transmitted over distances via allosteric processes. This ubiquitous protein process allows for protein function changes due to ligand binding events. Understanding protein allostery is essential to understanding protein functions. In this study, allostery in the second PDZ domain (PDZ2) in the human PTP1E protein is examined as model system to advance a recently developed rigid residue scan method combining with configurational entropy calculation and principal component analysis. The contributions from individual residues to whole-protein dynamics and allostery were systematically assessed via rigid body simulations of both unbound and ligand-bound states of the protein. The entropic contributions of individual residues to whole-protein dynamics were evaluated based on covariance-based correlation analysis of all simulations. The changes of overall protein entropy when individual residues being held rigid support that the rigidity/flexibility equilibrium in protein structure is governed by the La Châtelier's principle of chemical equilibrium. Key residues of PDZ2 allostery were identified with good agreement with NMR studies of the same protein bound to the same peptide. On the other hand, the change of entropic contribution from each residue upon perturbation revealed intrinsic differences among all the residues. The quasi-harmonic and principal component analyses of simulations without rigid residue perturbation showed a coherent allosteric mode from unbound and bound states, respectively. The projection of simulations with rigid residue perturbation onto coherent allosteric modes demonstrated the intrinsic shifting of ensemble distributions supporting the population-shift theory of protein allostery. Overall, the study presented here provides a robust and systematic approach to estimate the contribution of individual residue internal motion to overall protein dynamics and allostery.

Disorder transitions and conformational diversity cooperatively modulate biological function in proteins

Structural differences between conformers sustain protein biological function. Here, we studied in a large dataset of 745 intrinsically disordered proteins, how ordered-disordered transitions modulate structural differences between conformers as derived from crystallographic data. We found that almost 50% of the proteins studied show no transitions and have low conformational diversity while the rest show transitions and a higher conformational diversity. In this last subset, 60% of the proteins become more ordered after ligand binding, while 40% more disordered. As protein conformational diversity is inherently connected with protein function our analysis suggests differences in structure-function relationships related to order-disorder transitions.

Emerging Computational Methods for the Rational Discovery of Allosteric Drugs

Allosteric drug development holds promise for delivering medicines that are more selective and less toxic than those that target orthosteric sites. To date, the discovery of allosteric binding sites and lead compounds has been mostly serendipitous, achieved through high-throughput screening. Over the past decade, structural data has become more readily available for larger protein systems and more membrane protein classes (e.g., GPCRs and ion channels), which are common allosteric drug targets. In parallel, improved simulation methods now provide better atomistic understanding of the protein dynamics and cooperative motions that are critical to allosteric mechanisms. As a result of these advances, the field of predictive allosteric drug development is now on the cusp of a new era of rational structure-based computational methods. Here, we review algorithms that predict allosteric sites based on sequence data and molecular dynamics simulations, describe tools that assess the druggability of these pockets, and discuss how Markov state models and topology analyses provide insight into the relationship between protein dynamics and allosteric drug binding. In each section, we first provide an overview of the various method classes before describing relevant algorithms and software packages.

Correlated rigid modes in protein families

A great deal of evolutionarily conserved information is contained in genomes and proteins. Enormous effort has been put into understanding protein structure and developing computational tools for protein folding, and many sophisticated approaches take structure and sequence homology into account. Several groups have applied statistical physics approaches to extracting information about proteins from sequences alone. Here, we develop a new method for sequence analysis based on first principles, in information theory, in statistical physics and in Bayesian analysis. We provide a complete derivation of our approach and we apply it to a variety of systems, to demonstrate its utility and its limitations. We show in some examples that phylogenetic alignments of amino-acid sequences of families of proteins imply the existence of a small number of modes that appear to be associated with correlated global variation. These modes are uncovered efficiently in our approach by computing a non-perturbative effective potential directly from the alignment. We show that this effective potential approaches a limiting form inversely with the logarithm of the number of sequences. Mapping symbol entropy flows along modes to underlying physical structures shows that these modes arise due to correlated compensatory adjustments. In the protein examples, these occur around functional binding pockets.

Molecular Dynamics Simulations Reveal the Mechanisms of Allosteric Activation of Hsp90 by Designed Ligands

Controlling biochemical pathways through chemically designed modulators may provide novel opportunities to develop therapeutic drugs and chemical tools. The underlying challenge is to design new molecular entities able to act as allosteric chemical switches that selectively turn on/off functions by modulating the conformational dynamics of their target protein. We examine the origins of the stimulation of ATPase and closure kinetics in the molecular chaperone Hsp90 by allosteric modulators through atomistic molecular dynamics (MD) simulations and analysis of protein-ligand interactions. In particular, we focus on the cross-talk between allosteric ligands and protein conformations and its effect on the dynamic properties of the chaperone's active state. We examine the impact of different allosteric modulators on the stability, structural and internal dynamics properties of Hsp90 closed state. A critical aspect of this study is the development of a quantitative model that correlates Hsp90 activation to the presence of a certain compound, making use of information on the dynamic adaptation of protein conformations to the presence of the ligand, which allows to capture conformational states relevant in the activation process. We discuss the implications of considering the conformational dialogue between allosteric ligands and protein conformations for the design of new functional modulators.

Dynamic Allostery Mediated by a Conserved Tryptophan in the Tec Family Kinases

Bruton’s tyrosine kinase (Btk) is a Tec family non-receptor tyrosine kinase that plays a critical role in immune signaling and is associated with the immunological disorder X-linked agammaglobulinemia (XLA). Our previous findings showed that the Tec kinases are allosterically activated by the adjacent N-terminal linker. A single tryptophan residue in the N-terminal 17-residue linker mediates allosteric activation, and its mutation to alanine leads to the complete loss of activity. Guided by hydrogen/deuterium exchange mass spectrometry results, we have employed Molecular Dynamics simulations, Principal Component Analysis, Community Analysis and measures of node centrality to understand the details of how a single tryptophan mediates allostery in Btk. A specific tryptophan side chain rotamer promotes the functional dynamic allostery by inducing coordinated motions that spread across the kinase domain. Either a shift in the rotamer population, or a loss of the tryptophan side chain by mutation, drastically changes the coordinated motions and dynamically isolates catalytically important regions of the kinase domain. This work also identifies a new set of residues in the Btk kinase domain with high node centrality values indicating their importance in transmission of dynamics essential for kinase activation. Structurally, these node residues appear in both lobes of the kinase domain. In the N-lobe, high centrality residues wrap around the ATP binding pocket connecting previously described Catalytic-spine residues. In the C-lobe, two high centrality node residues connect the base of the R- and C-spines on the αF-helix. We suggest that the bridging residues that connect the catalytic and regulatory architecture within the kinase domain may be a crucial element in transmitting information about regulatory spine assembly to the catalytic machinery of the catalytic spine and active site.

Bruton’s tyrosine kinase (Btk) belongs to the Tec family of protein tyrosine kinases, and plays a crucial role in the signaling pathway in B-cells. Alteration of Btk activity results in the serious immunological disorder, X-linked agammaglobulinemia. Btk is a multi-domain protein and the activity of the kinase domain is regulated by the adjacent non-catalytic domains, which mediate their effect by means of a conserved tryptophan residue. In this work, we have investigated the mechanism of regulation by this tryptophan residue, W395, in the linker preceding the Btk kinase domain. Using hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations we identify structural elements within the kinase domain that are required for function by transmitting the allosteric effects of W395. Molecular Dynamics simulations further guided us to delineate the kinase domain into dynamically correlated sets of residues using community analysis, thereby identifying the important communication nodes that connect the various elements of the kinase domain required for function. The analyses performed indicate clearly how the W395A mutant changes the communication pathway required for function.

Function changing mutations in glucocorticoid receptor evolution correlate with their relevance to mode coupling

Nonlinear effects in protein dynamics are expected to play role in function, particularly of allosteric nature, by facilitating energy transfer between vibrational modes. A recently proposed method focusing on the non-Gaussian shape of the configurational population near equilibrium projects this information onto real space in order to identify the aminoacids relevant to function. We here apply this method to three ancestral proteins in glucocorticoid receptor (GR) family and show that the mutations that restrict functional activity during GR evolution correlate significantly with locations that are highlighted by the nonlinear contribution to the near-native configurational distribution. Our findings demonstrate that the analysis of nonlinear effects in protein dynamics can be harnessed into a predictive tool for functional site determination.

Structure-Based Statistical Mechanical Model Accounts for the Causality and Energetics of Allosteric Communication

Allostery is one of the pervasive mechanisms through which proteins in living systems carry out enzymatic activity, cell signaling, and metabolism control. Effective modeling of the protein function regulation requires a synthesis of the thermodynamic and structural views of allostery. We present here a structure-based statistical mechanical model of allostery, allowing one to observe causality of communication between regulatory and functional sites, and to estimate per residue free energy changes. Based on the consideration of ligand free and ligand bound systems in the context of a harmonic model, corresponding sets of characteristic normal modes are obtained and used as inputs for an allosteric potential. This potential quantifies the mean work exerted on a residue due to the local motion of its neighbors. Subsequently, in a statistical mechanical framework the entropic contribution to allosteric free energy of a residue is directly calculated from the comparison of conformational ensembles in the ligand free and ligand bound systems. As a result, this method provides a systematic approach for analyzing the energetics of allosteric communication based on a single structure. The feasibility of the approach was tested on a variety of allosteric proteins, heterogeneous in terms of size, topology and degree of oligomerization. The allosteric free energy calculations show the diversity of ways and complexity of scenarios existing in the phenomenology of allosteric causality and communication. The presented model is a step forward in developing the computational techniques aimed at detecting allosteric sites and obtaining the discriminative power between agonistic and antagonistic effectors, which are among the major goals in allosteric drug design.

The 50th anniversary of Monod-Changeux-Jacob seminal paper “Allosteric proteins and cellular control systems” became the hallmark of a new wave in the allostery studies and the turning point in our vision of allostery and its implications in protein engineering and drug design. Recent experimental and theoretical works clearly show relevance of allosteric phenomenon to drug design, unraveling advantages of allosteric drugs in comparison to traditional orthosteric compounds. Remarkable simplicity of allosteric effectors and, at the same time, their potentially high specificity is one of the most important traits. The non conserved nature of allosteric ligands is a basis for avoiding drug resistance, and existence of latent regulatory sites make them attractive drug targets. The model presented in this work provides a theoretical framework for the quantification of the causality and energetics of allosteric regulation, which is a prerequisite for design of effector molecules with required characteristics. The synthesis between the thermodynamics of allostery and the intrinsic atomic nature of proteins and their interactions with the allosteric effectors accomplished in this work is a small initial step in the long endeavor towards future allosteric drugs.

Hydrophobic Core Variations Provide a Structural Framework for Tyrosine Kinase Evolution and Functional Specialization

Protein tyrosine kinases (PTKs) are a group of closely related enzymes that have evolutionarily diverged from serine/threonine kinases (STKs) to regulate pathways associated with multi-cellularity. Evolutionary divergence of PTKs from STKs has occurred through accumulation of mutations in the active site as well as in the commonly conserved hydrophobic core. While the functional significance of active site variations is well understood, relatively little is known about how hydrophobic core variations contribute to PTK evolutionary divergence. Here, using a combination of statistical sequence comparisons, molecular dynamics simulations, mutational analysis and in vitro thermostability and kinase assays, we investigate the structural and functional significance of key PTK-specific variations in the kinase core. We find that the nature of residues and interactions in the hydrophobic core of PTKs is strikingly different from other protein kinases, and PTK-specific variations in the core contribute to functional divergence by altering the stability and dynamics of the kinase domain. In particular, a functionally critical STK-conserved histidine that stabilizes the regulatory spine in STKs is selectively mutated to an alanine, serine or glutamate in PTKs, and this loss-of-function mutation is accommodated, in part, through compensatory PTK-specific interactions in the core. In particular, a PTK-conserved phenylalanine in the I-helix appears to structurally and functionally compensate for the loss of STK-histidine by interacting with the regulatory spine, which has far-reaching effects on enzyme activity, inhibitor sensing, and stability. We propose that hydrophobic core variations provide a selective advantage during PTK evolution by increasing the conformational flexibility, and therefore the allosteric potential of the kinase domain. Our studies also suggest that Tyrosine Kinase Like kinases such as RAF are intermediates in PTK evolutionary divergence inasmuch as they share features of both PTKs and STKs in the core. Finally, our studies provide an evolutionary framework for identifying and characterizing disease and drug resistance mutations in the kinase core.

Proteins evolve new functions through accumulation of mutations in the primary sequence. Understanding how naturally occurring mutations shape protein function can provide insights into how non-natural mutations contribute to disease. Here, we identify sequence variants associated with the functional specialization of tyrosine kinases, a large and medically important class of proteins associated with the evolution of complex multicellular functions and diseases such as cancer. We find that mutations distal from the active site contribute to functional specialization by altering the stability, activity and dynamics of the kinase core. Our findings have implications for understanding the evolution of allosteric regulation in tyrosine kinases, and in predicting the structural and functional impact of disease and drug resistance mutations in the kinase core.

Allosteric Mechanisms of Molecular Machines at the Membrane: Transport by Sodium-Coupled Symporters

Solute transport across cell membranes is ubiquitous in biology as an essential physiological process. Secondary active transporters couple the unfavorable process of solute transport against its concentration gradient to the energetically favorable transport of one or several ions. The study of such transporters over several decades indicates that their function involves complex allosteric mechanisms that are progressively being revealed in atomistic detail. We focus on two well-characterized sodium-coupled symporters: the bacterial amino acid transporter LeuT, which is the prototype for the "gated pore" mechanism in the mammalian synaptic monoamine transporters, and the archaeal GltPh, which is the prototype for the "elevator" mechanism in the mammalian excitatory amino acid transporters. We present the evidence for the role of allostery in the context of a quantitative formalism that can reconcile biochemical and biophysical data and thereby connects directly to recent insights into the molecular structure and dynamics of these proteins. We demonstrate that, while the structures and mechanisms of these transporters are very different, the available data suggest a common role of specific models of allostery in their functions. We argue that such allosteric mechanisms appear essential not only for sodium-coupled symport in general but also for the function of other types of molecular machines in the membrane.

Controlling Allosteric Networks in Proteins

Allosteric transition, defined as conformational changes induced by ligand binding, is one of the fundamental properties of proteins. Allostery has been observed and characterized in many proteins, and has been recently utilized to control protein function via regulation of protein activity. Here, we review the physical and evolutionary origin of protein allostery, as well as its importance to protein regulation, drug discovery, and biological processes in living systems. We describe recently developed approaches to identify allosteric pathways, connected sets of pairwise interactions that are responsible for propagation of conformational change from the ligand-binding site to a distal functional site. We then present experimental and computational protein engineering approaches for control of protein function by modulation of allosteric sites. As an example of application of these approaches, we describe a synergistic computational and experimental approach to rescue the cystic-fibrosis-associated protein cystic fibrosis transmembrane conductance regulator, which upon deletion of a single residue misfolds and causes disease. This example demonstrates the power of allosteric manipulation in proteins to both elucidate mechanisms of molecular function and to develop therapeutic strategies that rescue those functions. Allosteric control of proteins provides a tool to shine a light on the complex cascades of cellular processes and facilitate unprecedented interrogation of biological systems.

Conformational Dynamics and Allostery in Pyruvate Kinase

Pyruvate kinase catalyzes the final step in glycolysis and is allosterically regulated to control flux through the pathway. Two models are proposed to explain how Escherichia coli pyruvate kinase type 1 is allosterically regulated: the "domain rotation model" suggests that both the domains within the monomer and the monomers within the tetramer reorient with respect to one another; the "rigid body reorientation model" proposes only a reorientation of the monomers within the tetramer causing rigidification of the active site. To test these hypotheses and elucidate the conformational and dynamic changes that drive allostery, we performed time-resolved electrospray ionization mass spectrometry coupled to hydrogen-deuterium exchange studies followed by mutagenic analysis to test the activation mechanism. Global exchange experiments, supported by thermostability studies, demonstrate that fructose 1,6-bisphosphate binding to the allosteric domain causes a shift toward a globally more dynamic ensemble of conformations. Mapping deuterium exchange to peptides within the enzyme highlight site-specific regions with altered conformational dynamics, many of which increase in conformational flexibility. Based upon these and mutagenic studies, we propose an allosteric mechanism whereby the binding of fructose 1,6-bisphosphate destabilizes an α-helix that bridges the allosteric and active site domains within the monomeric unit. This destabilizes the β-strands within the (β/α)8-barrel domain and the linked active site loops that are responsible for substrate binding. Our data are consistent with the domain rotation model but inconsistent with the rigid body reorientation model given the increased flexibility at the interdomain interface, and we can for the first time explain how fructose 1,6-bisphosphate affects the active site.

Using Local States To Drive the Sampling of Global Conformations in Proteins

Conformational changes associated with protein function often occur beyond the time scale currently accessible to unbiased molecular dynamics (MD) simulations, so that different approaches have been developed to accelerate their sampling. Here we investigate how the knowledge of backbone conformations preferentially adopted by protein fragments, as contained in precalculated libraries known as structural alphabets (SA), can be used to explore the landscape of protein conformations in MD simulations. We find that (a) enhancing the sampling of native local states in both metadynamics and steered MD simulations allows the recovery of global folded states in small proteins; (b) folded states can still be recovered when the amount of information on the native local states is reduced by using a low-resolution version of the SA, where states are clustered into macrostates; and (c) sequences of SA states derived from collections of structural motifs can be used to sample alternative conformations of preselected protein regions. The present findings have potential impact on several applications, ranging from protein model refinement to protein folding and design.

Mismatch binding, ADP-ATP exchange and intramolecular signaling during mismatch repair

The focus of this article is on the DNA binding and ATPase activities of the mismatch repair (MMR) protein, MutS-our current understanding of how this protein uses ATP to fuel its actions on DNA and initiate repair via interactions with MutL, the next protein in the pathway. Structure-function and kinetic studies have yielded detailed views of the MutS mechanism of action in MMR. How MutS and MutL work together after mismatch recognition to enable strand-specific nicking, which leads to strand excision and synthesis, is less clear and remains an active area of investigation.

Mechanisms of Differential Allosteric Modulation in Homologous Proteins: Insights from the Analysis of Internal Dynamics and Energetics of PDZ Domains

Allostery is a general phenomenon in proteins whereby a perturbation at one site reverberates into a functional change at another one, through modulation of its conformational dynamics. Herein, we address the problem of how the molecular signal encoded by a ligand is differentially transmitted through the structures of two homologous PDZ proteins: PDZ2, which responds to binding with structural and dynamical changes in regions distal from the ligand site, and PDZ3, which is characterized by less-intense dynamical variations. We use novel methods of analysis of MD simulations in the unbound and bound states to investigate the determinants of the differential allosteric behavior of the two proteins. The analysis of the correlations between the redistribution of stabilization energy and local fluctuation patterns highlights the nucleus of residues responsible for the stabilization of the 3D fold, the stability core, as the substructure that defines the difference in the allosteric response: in PDZ2, it undergoes a consistent dynamic and energetic reorganization, whereas in PDZ3, it remains largely unperturbed. Specifically, we observe for PDZ2 a significant anticorrelation between the motions of distal loops and residues of the stability core and differences in the correlation patterns between the bound and unbound states. Such variation is not observed in PDZ3, indicating that its energetics and internal dynamics are less affected by the presence/absence of the ligand. Finally, we propose a model with a direct link between the modulation of the structural, energetic and dynamic properties of a protein, and its allosteric response to a perturbation.

Dancing through Life: Molecular Dynamics Simulations and Network-Centric Modeling of Allosteric Mechanisms in Hsp70 and Hsp110 Chaperone Proteins

Hsp70 and Hsp110 chaperones play an important role in regulating cellular processes that involve protein folding and stabilization, which are essential for the integrity of signaling networks. Although many aspects of allosteric regulatory mechanisms in Hsp70 and Hsp110 chaperones have been extensively studied and significantly advanced in recent experimental studies, the atomistic picture of signal propagation and energetics of dynamics-based communication still remain unresolved. In this work, we have combined molecular dynamics simulations and protein stability analysis of the chaperone structures with the network modeling of residue interaction networks to characterize molecular determinants of allosteric mechanisms. We have shown that allosteric mechanisms of Hsp70 and Hsp110 chaperones may be primarily determined by nucleotide-induced redistribution of local conformational ensembles in the inter-domain regions and the substrate binding domain. Conformational dynamics and energetics of the peptide substrate binding with the Hsp70 structures has been analyzed using free energy calculations, revealing allosteric hotspots that control negative cooperativity between regulatory sites. The results have indicated that cooperative interactions may promote a population-shift mechanism in Hsp70, in which functional residues are organized in a broad and robust allosteric network that can link the nucleotide-binding site and the substrate-binding regions. A smaller allosteric network in Hsp110 structures may elicit an entropy-driven allostery that occurs in the absence of global structural changes. We have found that global mediating residues with high network centrality may be organized in stable local communities that are indispensable for structural stability and efficient allosteric communications. The network-centric analysis of allosteric interactions has also established that centrality of functional residues could correlate with their sensitivity to mutations across diverse chaperone functions. This study reconciles a wide spectrum of structural and functional experiments by demonstrating how integration of molecular simulations and network-centric modeling may explain thermodynamic and mechanistic aspects of allosteric regulation in chaperones.

Small Molecule Targeting of Protein-Protein Interactions through Allosteric Modulation of Dynamics

The protein–protein interaction (PPI) target class is particularly challenging, but offers potential for “first in class” therapies. Most known PPI small molecules are orthosteric inhibitors but many PPI sites may be fundamentally intractable to this approach. One potential alternative is to consider more attractive, remote small molecule pockets; however, on the whole, allostery is poorly understood and difficult to discover and develop. Here we review the literature in order to understand the basis for allostery, especially as it can apply to PPIs. We suggest that the upfront generation of sophisticated and experimentally validated dynamic models of target proteins can aid in target choice and strategy for allosteric intervention to produce the required functional effect.

Cyclophilin40 isomerase activity is regulated by a temperature-dependent allosteric interaction with Hsp90

Cyclophilin 40 (Cyp40) comprises an N-terminal cyclophilin domain with peptidyl-prolyl isomerase (PPIase) activity and a C-terminal tetratricopeptide repeat (TPR) domain that binds to the C-terminal-EEVD sequence common to both heat shock protein 70 (Hsp70) and Hsp90. We show in the present study that binding of peptides containing the MEEVD motif reduces the PPIase activity by ∼30%. CD and fluorescence assays show that the TPR domain is less stable than the cyclophilin domain and is stabilized by peptide binding. Isothermal titration calorimetry (ITC) shows that the affinity for the-MEEVD peptide is temperature sensitive in the physiological temperature range. Results from these biophysical studies fit with the MD simulations of the apo and holo (peptide-bound) structures which show a significant reduction in root mean square (RMS) fluctuation in both TPR and cyclophilin domains when-MEEVD is bound. The MD simulations of the apo-protein also highlight strong anti-correlated motions between residues around the PPIase-active site and a band of residues running across four of the seven helices in the TPR domain. Peptide binding leads to a distortion in the shape of the active site and a significant reduction in these strongly anti-correlated motions, providing an explanation for the allosteric effect of ligand binding and loss of PPIase activity. Together the experimental and MD results suggest that on heat shock, dissociation of Cyp40 from complexes mediated by the TPR domain leads to an increased pool of free Cyp40 capable of acting as an isomerase/chaperone in conditions of cellular stress.

The Mechanism of ATP-Dependent Allosteric Protection of Akt Kinase Phosphorylation

Kinases use ATP to phosphorylate substrates; recent findings underscore the additional regulatory roles of ATP. Here, we propose a mechanism for allosteric regulation of Akt1 kinase phosphorylation by ATP. Our 4.7-μs molecular dynamics simulations of Akt1 and its mutants in the ATP/ADP bound/unbound states revealed that ATP occupancy of the ATP-binding site stabilizes the closed conformation, allosterically protecting pT308 by restraining phosphatase access and key interconnected residues on the ATP→pT308 allosteric pathway. Following ATP→ADP hydrolysis, pT308 is exposed and readily dephosphorylated. Site-directed mutagenesis validated these predictions and indicated that the mutations do not impair PDK1 and PP2A phosphatase recruitment. We further probed the function of residues around pT308 at the atomic level, and predicted and experimentally confirmed that Akt1(H194R/R273H) double mutant rescues pathology-related Akt1(R273H). Analysis of classical Akt homologs suggests that this mechanism can provide a general model of allosteric kinase regulation by ATP; as such, it offers a potential avenue for allosteric drug discovery.

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.

Integrating atomistic molecular dynamics simulations, experiments, and network analysis to study protein dynamics: strength in unity

In the last years, we have been observing remarkable improvements in the field of protein dynamics. Indeed, we can now study protein dynamics in atomistic details over several timescales with a rich portfolio of experimental and computational techniques. On one side, this provides us with the possibility to validate simulation methods and physical models against a broad range of experimental observables. On the other side, it also allows a complementary and comprehensive view on protein structure and dynamics. What is needed now is a better understanding of the link between the dynamic properties that we observe and the functional properties of these important cellular machines. To make progresses in this direction, we need to improve the physical models used to describe proteins and solvent in molecular dynamics, as well as to strengthen the integration of experiments and simulations to overcome their own limitations. Moreover, now that we have the means to study protein dynamics in great details, we need new tools to understand the information embedded in the protein ensembles and in their dynamic signature. With this aim in mind, we should enrich the current tools for analysis of biomolecular simulations with attention to the effects that can be propagated over long distances and are often associated to important biological functions. In this context, approaches inspired by network analysis can make an important contribution to the analysis of molecular dynamics simulations.

AIM for Allostery: Using the Ising Model to Understand Information Processing and Transmission in Allosteric Biomolecular Systems

In performing their biological functions, molecular machines must process and transmit information with high fidelity. Information transmission requires dynamic coupling between the conformations of discrete structural components within the protein positioned far from one another on the molecular scale. This type of biomolecular “action at a distance” is termed allostery. Although allostery is ubiquitous in biological regulation and signal transduction, its treatment in theoretical models has mostly eschewed quantitative descriptions involving the system's underlying structural components and their interactions. Here, we show how Ising models can be used to formulate an approach to allostery in a structural context of interactions between the constitutive components by building simple allosteric constructs we termed Allosteric Ising Models (AIMs). We introduce the use of AIMs in analytical and numerical calculations that relate thermodynamic descriptions of allostery to the structural context, and then show that many fundamental properties of allostery, such as the multiplicative property of parallel allosteric channels, are revealed from the analysis of such models. The power of exploring mechanistic structural models of allosteric function in more complex systems by using AIMs is demonstrated by building a model of allosteric signaling for an experimentally well-characterized asymmetric homodimer of the dopamine D2 receptor.

Oligomerization and nanocluster organization render specificity

Nanoclusters are anchored to membranes, either within them or in the cytoplasm latched onto the cytoskeleton, whose reorganization can regulate their activity. Nanoclusters have been viewed in terms of cooperativity and activation; here we perceive nanocluster organization from a conformational standpoint. This leads us to suggest that while single molecules encode activity, nanoclusters induce specificity, and that this is their main evolutionary aim. Distinct, isoform-specific nanocluster organization can drive the preferred effector (and ligand) interactions and thereby designate signalling pathways. The absence of detailed structural information across the nanocluster, due to size and dynamics, hinders an in-depth grasp of its mechanistic features; however, available data already capture some of the principles and their functional 'raison d'être'. Collectively, clustering lends stability and reduces the likelihood of proteolytic cleavage; it also increases the effective local concentration and enables efficient cooperative activation. However, clustering does not determine the ability of the single molecule to function. Drugs targeting nanoclusters can attenuate activity by hampering cooperativity; however, this may not perturb activation and signalling, which originate from the molecules themselves, and as such, are likely to endure. What then is the major role of nanoclustering? Assuming that single molecules evolved first, with a subsequent increase in cellular complexity and emergence of highly similar isoform variants, evolution faced the threat of signalling promiscuity. We reason that this potential risk was thwarted by oligomerization and clustering; clustering confers higher specificity, and a concomitant extra layer of cellular control. In our Ras example, signalling will be more accurate as a dimer than as a monomer, where its isomer specificity could be compromised.

Mechanism of the Association between Na+ Binding and Conformations at the Intracellular Gate in Neurotransmitter:Sodium Symporters

Neurotransmitter:sodium symporters (NSSs) terminate neurotransmission by Na⁺-dependent reuptake of released neurotransmitters. Previous studies suggested that Na⁺-binding reconfigures dynamically coupled structural elements in an allosteric interaction network (AIN) responsible for function-related conformational changes, but the intramolecular pathway of this mechanism has remained uncharted. We describe a new approach for the modeling and analysis of intramolecular dynamics in the bacterial NSS homolog LeuT. From microsecond-scale molecular dynamics simulations and cognate experimental verifications in both LeuT and human dopamine transporter (hDAT), we apply the novel method to identify the composition and the dynamic properties of their conserved AIN. In LeuT, two different perturbations disrupting Na⁺ binding and transport (i.e. replacing Na⁺ with Li⁺ or the Y268A mutation at the intracellular gate) affect the AIN in strikingly similar ways. In contrast, other mutations that affect the intracellular gate (i.e. R5A and D369A) do not significantly impair Na⁺ cooperativity and transport. Our analysis shows these perturbations to have much lesser effects on the AIN, underscoring the sensitivity of this novel method to the mechanistic nature of the perturbation. Notably, this set of observations holds as well for hDAT, where the aligned Y335A, R60A, and D436A mutations also produce different impacts on Na⁺ dependence. Thus, the detailed AIN generated from our method is shown to connect Na⁺ binding with global conformational changes that are critical for the transport mechanism. That the AIN between the Na⁺ binding sites and the intracellular gate in bacterial LeuT resembles that in eukaryotic hDAT highlights the conservation of allosteric pathways underlying NSS function.