The energetic and allosteric landscape for KRAS inhibition.

Thousands of proteins have been validated genetically as therapeutic targets for human diseases1. However, very few have been successfully targeted, and many are considered 'undruggable'. This is particularly true for proteins that function via protein-protein interactions-direct inhibition of binding interfaces is difficult and requires the identification of allosteric sites. However, most proteins have no known allosteric sites, and a comprehensive allosteric map does not exist for any protein. Here we address this shortcoming by charting multiple global atlases of inhibitory allosteric communication in KRAS. We quantified the effects of more than 26,000 mutations on the folding of KRAS and its binding to six interaction partners. Genetic interactions in double mutants enabled us to perform biophysical measurements at scale, inferring more than 22,000 causal free energy changes. These energy landscapes quantify how mutations tune the binding specificity of a signalling protein and map the inhibitory allosteric sites for an important therapeutic target. Allosteric propagation is particularly effective across the central ß-sheet of KRAS, and multiple surface pockets are genetically validated as allosterically active, including a distal pocket in the C-terminal lobe of the protein. Allosteric mutations typically inhibit binding to all tested effectors, but they can also change the binding specificity, revealing the regulatory, evolutionary and therapeutic potential to tune pathway activation. Using the approach described here, it should be possible to rapidly and comprehensively identify allosteric target sites in many proteins.

Genetic switching by the Lac repressor is based on two-state Monod-Wyman-Changeux allostery.

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod-Wyman-Changeux model.

Disease-causing cystathionine ß-synthase linker mutations impair allosteric regulation.

Cystathionine ß-synthase (CBS) catalyzes the committing step in the transsulfuration pathway, which is important for clearing homocysteine and furnishing cysteine. The transsulfuration pathway also generates H2S, a signaling molecule. CBS is a modular protein with a heme and pyridoxal phosphate-binding catalytic core, which is separated by a linker region from the C-terminal regulatory domain that binds S-adenosylmethionine (AdoMet), an allosteric activator. Recent cryo-EM structures reveal that CBS exists in a fibrillar form and undergoes a dramatic architectural rearrangement between the basal and AdoMet-bound states. CBS is the single most common locus of mutations associated with homocystinuria, and, in this study, we have characterized three clinical variants (K384E/N and M391I), which reside in the linker region. The native fibrillar form is destabilized in the variants, and differences in their limited proteolytic fingerprints also reveal conformational alterations. The crystal structure of the truncated K384N variant, lacking the regulatory domain, reveals that the overall fold of the catalytic core is unperturbed. M391I CBS exhibits a modest (1.4-fold) decrease while the K384E/N variants exhibit a significant (∼8-fold) decrease in basal activity, which is either unresponsive to or inhibited by AdoMet. Pre-steady state kinetic analyses reveal that the K384E/N substitutions exhibit pleiotropic effects and that the differences between them are expressed in the second half reaction, that is, homocysteine binding and reaction with the aminoacrylate intermediate. Together, these studies point to an important role for the linker in stabilizing the higher-order oligomeric structure of CBS and enabling AdoMet-dependent regulation.

A critical inter-subunit interaction for the transmission of the allosteric signal in the Agrobacterium tumefaciens ADP-glucose pyrophosphorylase.

ADP-glucose pyrophosphorylase is a key regulatory enzyme involved in starch and glycogen synthesis in plants and bacteria, respectively. It has been hypothesized that inter-subunit communications are important for the allosteric effect in this enzyme. However, no specific interactions have been identified as part of the regulatory signal. The enzyme from Agrobacterium tumefaciens is a homotetramer allosterically regulated by fructose 6-phosphate and pyruvate. Three pairs of distinct subunit-subunit interfaces are present. Here we focus on an interface that features two symmetrical interactions between Arg11 and Asp141 from one subunit with residues Asp141 and Arg11 of the neighbor subunit, respectively. Previously, scanning mutagenesis showed that a mutation at the Arg11 position disrupted the activation of the enzyme. Considering the distance of these residues from the allosteric and catalytic sites, we hypothesized that the interaction between Arg11 and Asp141 is critical for allosteric signaling rather than effector binding. To prove our hypothesis, we mutated those two sites (D141A, D141E, D141N, D141R, R11D, and R11K) and performed kinetic and binding analysis. Mutations that altered the charge affected the regulation the most. To prove that the interaction per se (rather than the presence of specific residues) is critical, we partially rescued the R11D protein by introducing a second mutation (R11D/D141R). This could not restore the activator effect on kcat , but it did rescue the effect on substrate affinity. Our results indicate the critical functional role of Arg11 and Asp141 to relay the allosteric signal in this subunit interface.

DeepAlloDriver: a deep learning-based strategy to predict cancer driver mutations.

Driver mutations can contribute to the initial processes of cancer, and their identification is crucial for understanding tumorigenesis as well as for molecular drug discovery and development. Allostery regulates protein function away from the functional regions at an allosteric site. In addition to the known effects of mutations around functional sites, mutations at allosteric sites have been associated with protein structure, dynamics, and energy communication. As a result, identifying driver mutations at allosteric sites will be beneficial for deciphering the mechanisms of cancer and developing allosteric drugs. In this study, we provided a platform called DeepAlloDriver to predict driver mutations using a deep learning method that exhibited >93% accuracy and precision. Using this server, we found that a missense mutation in RRAS2 (Gln72 to Leu) might serve as an allosteric driver of tumorigenesis, revealing the mechanism of the mutation in knock-in mice and cancer patients. Overall, DeepAlloDriver would facilitate the elucidation of the mechanisms underlying cancer progression and help prioritize cancer therapeutic targets. The web server is freely available at: https://mdl.shsmu.edu.cn/DeepAlloDriver.

Coevolution-based prediction of key allosteric residues for protein function regulation.

Allostery is fundamental to many biological processes. Due to the distant regulation nature, how allosteric mutations, modifications, and effector binding impact protein function is difficult to forecast. In protein engineering, remote mutations cannot be rationally designed without large-scale experimental screening. Allosteric drugs have raised much attention due to their high specificity and possibility of overcoming existing drug-resistant mutations. However, optimization of allosteric compounds remains challenging. Here, we developed a novel computational method KeyAlloSite to predict allosteric site and to identify key allosteric residues (allo-residues) based on the evolutionary coupling model. We found that protein allosteric sites are strongly coupled to orthosteric site compared to non-functional sites. We further inferred key allo-residues by pairwise comparing the difference of evolutionary coupling scores of each residue in the allosteric pocket with the functional site. Our predicted key allo-residues are in accordance with previous experimental studies for typical allosteric proteins like BCR-ABL1, Tar, and PDZ3, as well as key cancer mutations. We also showed that KeyAlloSite can be used to predict key allosteric residues distant from the catalytic site that are important for enzyme catalysis. Our study demonstrates that weak coevolutionary couplings contain important information of protein allosteric regulation function. KeyAlloSite can be applied in studying the evolution of protein allosteric regulation, designing and optimizing allosteric drugs, and performing functional protein design and enzyme engineering.

A complete allosteric map of a GTPase switch in its native cellular network.

Allosteric regulation is central to protein function in cellular networks. A fundamental open question is whether cellular regulation of allosteric proteins occurs only at a few defined positions or at many sites distributed throughout the structure. Here, we probe the regulation of GTPases-protein switches that control signaling through regulated conformational cycling-at residue-level resolution by deep mutagenesis in the native biological network. For the GTPase Gsp1/Ran, we find that 28% of the 4,315 assayed mutations show pronounced gain-of-function responses. Twenty of the sixty positions enriched for gain-of-function mutations are outside the canonical GTPase active site switch regions. Kinetic analysis shows that these distal sites are allosterically coupled to the active site. We conclude that the GTPase switch mechanism is broadly sensitive to cellular allosteric regulation. Our systematic discovery of new regulatory sites provides a functional map to interrogate and target GTPases controlling many essential biological processes.

AlloMAPS 2: allosteric fingerprints of the AlphaFold and Pfam-trRosetta predicted structures for engineering and design.

AlloMAPS 2 is an update of the Allosteric Mutation Analysis and Polymorphism of Signalling database, which contains data on allosteric communication obtained for predicted structures in the AlphaFold database (AFDB) and trRosetta-predicted Pfam domains. The data update contains Allosteric Signalling Maps (ASMs) and Allosteric Probing Maps (APMs) quantifying allosteric effects of mutations and of small probe binding, respectively. To ensure quality of the ASMs and APMs, we performed careful and accurate selection of protein sets containing high-quality predicted structures in both databases for each organism/structure, and the data is available for browsing and download. The data for remaining structures are available for download and should be used at user's discretion and responsibility. We believe these massive data can facilitate both diagnostics and drug design within the precision medicine paradigm. Specifically, it can be instrumental in the analysis of allosteric effects of pathological and rescue mutations, providing starting points for fragment-based design of allosteric effectors. The exhaustive character of allosteric signalling and probing fingerprints will be also useful in future developments of corresponding machine learning applications. The database is freely available at: http://allomaps.bii.a-star.edu.sg.

Cancer-related Mutations with Local or Long-range Effects on an Allosteric Loop of p53.

The tumor protein 53 (p53) is involved in transcription-dependent and independent processes. Several p53 variants related to cancer have been found to impact protein stability. Other variants, on the contrary, might have little impact on structural stability and have local or long-range effects on the p53 interactome. Our group previously identified a loop in the DNA binding domain (DBD) of p53 (residues 207-213) which can recruit different interactors. Experimental structures of p53 in complex with other proteins strengthen the importance of this interface for protein-protein interactions. We here characterized with structure-based approaches somatic and germline variants of p53 which could have a marginal effect in terms of stability and act locally or allosterically on the region 207-213 with consequences on the cytosolic functions of this protein. To this goal, we studied 1132 variants in the p53 DBD with structure-based approaches, accounting also for protein dynamics. We focused on variants predicted with marginal effects on structural stability. We then investigated each of these variants for their impact on DNA binding, dimerization of the p53 DBD, and intramolecular contacts with the 207-213 region. Furthermore, we identified variants that could modulate long-range the conformation of the region 207-213 using a coarse-grain model for allostery and all-atom molecular dynamics simulations. Our predictions have been further validated using enhanced sampling methods for 15 variants. The methodologies used in this study could be more broadly applied to other p53 variants or cases where conformational changes of loop regions are essential in the function of disease-related proteins.

Positive Allosteric Modulation of α5-GABAA Receptors Reverses Stress-Induced Alterations in Dopamine System Function and Prepulse Inhibition of Startle.

BACKGROUND: Up to 64% of patients diagnosed with posttraumatic stress disorder (PTSD) experience psychosis, likely attributable to aberrant dopamine neuron activity. We have previously demonstrated that positive allosteric modulators of α5-GABAARs can selectively decrease hippocampal activity and reverse psychosis-like physiological and behavioral alterations in a rodent model used to study schizophrenia; however, whether this approach translates to a PTSD model remains to be elucidated. METHODS: We utilized a 2-day inescapable foot shock (IS) procedure to induce stress-related pathophysiology in male Sprague-Dawley rats. We evaluated the effects of intra-ventral hippocampus (vHipp) administration GL-II-73, an α5-GABAAR, or viral overexpression of the α5 subunit, using in vivo electrophysiology and behavioral measures in control and IS-treated rats. RESULTS: IS significantly increased ventral tegmental area dopamine neuron population activity, or the number of dopamine neurons firing spontaneously (n = 6; P = .016), consistent with observation in multiple rodent models used to study psychosis. IS also induced deficits in sensorimotor gating, as measured by reduced prepulse inhibition of startle (n = 12; P = .039). Interestingly, intra-vHipp administration of GL-II-73 completely reversed IS-induced increases in dopamine neuron population activity (n = 6; P = .024) and deficits in prepulse inhibition (n = 8; P = .025), whereas viral overexpression of the α5 subunit in the vHipp was not effective. CONCLUSIONS: Our results demonstrate that pharmacological intervention augmenting α5-GABAAR function, but not α5 overexpression in itself, can reverse stress-induced deficits related to PTSD in a rodent model, providing a potential site of therapeutic intervention to treat comorbid psychosis in PTSD.

Allostery and Missense Mutations as Intermittently Linked Promising Aspects of Modern Computational Drug Discovery.

Drug research and development is a multidisciplinary field with its own successes. Yet, given the complexity of the process, it also faces challenges over the long development stages and even includes those that develop once a drug is marketed, i.e. drug toxicity and drug resistance. Better success can be achieved via well designed criteria in the early drug development stages. Here, we introduce the concepts of allostery and missense mutations, and argue that incorporation of these two intermittently linked biological phenomena into the early computational drug discovery stages would help to reduce the attrition risk in later stages of the process. We discuss the individual or in concert mechanisms of actions of mutations in allostery. Design of allosteric drugs is challenging compared to orthosteric drugs, yet they have been gaining popularity in recent years as alternative systems for the therapeutic regulation of proteins with an action-at-a-distance mode and non-invasive mechanisms. We propose an easy-to-apply computational allosteric drug discovery protocol which considers the mutation effect, and detail it with three case studies focusing on (1) analysis of effect of an allosteric mutation related to isoniazid drug resistance in tuberculosis; (2) identification of a cryptic pocket in the presence of an allosteric mutation of falcipain-2 as a malarial drug target; and (3) deciphering the effects of SARS-CoV-2 evolutionary mutations on a potential allosteric modulator with changes to allosteric communication paths.

Mapping the energetic and allosteric landscapes of protein binding domains.

Allosteric communication between distant sites in proteins is central to biological regulation but still poorly characterized, limiting understanding, engineering and drug development1-6. An important reason for this is the lack of methods to comprehensively quantify allostery in diverse proteins. Here we address this shortcoming and present a method that uses deep mutational scanning to globally map allostery. The approach uses an efficient experimental design to infer en masse the causal biophysical effects of mutations by quantifying multiple molecular phenotypes-here we examine binding and protein abundance-in multiple genetic backgrounds and fitting thermodynamic models using neural networks. We apply the approach to two of the most common protein interaction domains found in humans, an SH3 domain and a PDZ domain, to produce comprehensive atlases of allosteric communication. Allosteric mutations are abundant, with a large mutational target space of network-altering 'edgetic' variants. Mutations are more likely to be allosteric closer to binding interfaces, at glycine residues and at specific residues connecting to an opposite surface within the PDZ domain. This general approach of quantifying mutational effects for multiple molecular phenotypes and in multiple genetic backgrounds should enable the energetic and allosteric landscapes of many proteins to be rapidly and comprehensively mapped.

Allostery: Allosteric Cancer Drivers and Innovative Allosteric Drugs.

Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.

Engineering and exploiting synthetic allostery of NanoLuc luciferase.

Allostery enables proteins to interconvert different biochemical signals and form complex metabolic and signaling networks. We hypothesize that circular permutation of proteins increases the probability of functional coupling of new N- and C- termini with the protein's active center through increased local structural disorder. To test this we construct a synthetically allosteric version of circular permutated NanoLuc luciferase that can be activated through ligand-induced intramolecular non-covalent cyclisation. This switch module is tolerant of the structure of binding domains and their ligands, and can be used to create biosensors of proteins and small molecules. The developed biosensors covers a range of emission wavelengths and displays sensitivity as low as 50pM and dynamic range as high as 16-fold and could quantify their cognate ligand in human fluids. We apply hydrogen exchange kinetic mass spectroscopy to analyze time resolved structural changes in the developed biosensors and observe ligand-mediated folding of newly created termini.

Known allosteric proteins have central roles in genetic disease.

Allostery is a form of protein regulation, where ligands that bind sites located apart from the active site can modify the activity of the protein. The molecular mechanisms of allostery have been extensively studied, because allosteric sites are less conserved than active sites, and drugs targeting them are more specific than drugs binding the active sites. Here we quantify the importance of allostery in genetic disease. We show that 1) known allosteric proteins are central in disease networks, contribute to genetic disease and comorbidities much more than non-allosteric proteins, and there is an association between being allosteric and involvement in disease; 2) they are enriched in many major disease types like hematopoietic diseases, cardiovascular diseases, cancers, diabetes, or diseases of the central nervous system; 3) variants from cancer genome-wide association studies are enriched near allosteric proteins, indicating their importance to polygenic traits; and 4) the importance of allosteric proteins in disease is due, at least partly, to their central positions in protein-protein interaction networks, and less due to their dynamical properties.

Structural insights into the activation of chemokine receptor CXCR2.

The C-X-C motif chemokine CXCL8 (interleukin-8, IL-8) and its receptor chemokine receptor 2 (CXCR2) mediate neutrophil migration during cell development and inflammatory responses and thus are related to numerous inflammatory diseases and cancers. We have determined the cryo-electron microscopy structure of CXCL8 bound CXCR2 coupled to Gi protein, as well as the crystal structure of inactive CXCR2 in complex with a designed allosteric antagonist. These results reveal the binding modes between CXCL8 and CXCR2, CXCR2 and G protein, and the detailed binding pattern of the allosteric antagonist, 00767013. Further structural analysis of the inactive- and active- states of CXCR2 reveals the unique shallow-pocket activation mechanism of C-X-C chemokine receptors and promotes our understanding on how a G protein-coupled receptor (GPCR) is activated by an endogenous protein molecule. In addition, the cholesterol molecule is observed in the activated CXCR2 structure, providing the structural basis of the potential allosteric modulation role of cholesterol in chemokine receptors.

Allosteric mechanism of signal transduction in the two-component system histidine kinase PhoQ.

Transmembrane signaling proteins couple extracytosolic sensors to cytosolic effectors. Here, we examine how binding of Mg2+ to the sensor domain of an E. coli two component histidine kinase (HK), PhoQ, modulates its cytoplasmic kinase domain. We use cysteine-crosslinking and reporter-gene assays to simultaneously and independently probe the signaling state of PhoQ's sensor and autokinase domains in a set of over 30 mutants. Strikingly, conservative single-site mutations distant from the sensor or catalytic site strongly influence PhoQ's ligand-sensitivity as well as the magnitude and direction of the signal. Data from 35 mutants are explained by a semi-empirical three-domain model in which the sensor, intervening HAMP, and catalytic domains can adopt kinase-promoting or inhibiting conformations that are in allosteric communication. The catalytic and sensor domains intrinsically favor a constitutively 'kinase-on' conformation, while the HAMP domain favors the 'off' state; when coupled, they create a bistable system responsive to physiological concentrations of Mg2+. Mutations alter signaling by locally modulating domain intrinsic equilibrium constants and interdomain couplings. Our model suggests signals transmit via interdomain allostery rather than propagation of a single concerted conformational change, explaining the diversity of signaling structural transitions observed in individual HK domains.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9.

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

Biological signal processing filters via engineering allosteric transcription factors.

Signal processing is critical to a myriad of biological phenomena (natural and engineered) that involve gene regulation. Biological signal processing can be achieved by way of allosteric transcription factors. In canonical regulatory systems (e.g., the lactose repressor), an INPUT signal results in the induction of a given transcription factor and objectively switches gene expression from an OFF state to an ON state. In such biological systems, to revert the gene expression back to the OFF state requires the aggressive dilution of the input signal, which can take 1 or more d to achieve in a typical biotic system. In this study, we present a class of engineered allosteric transcription factors capable of processing two-signal INPUTS, such that a sequence of INPUTS can rapidly transition gene expression between alternating OFF and ON states. Here, we present two fundamental biological signal processing filters, BANDPASS and BANDSTOP, that are regulated by D-fucose and isopropyl-ß-D-1-thiogalactopyranoside. BANDPASS signal processing filters facilitate OFF-ON-OFF gene regulation. Whereas, BANDSTOP filters facilitate the antithetical gene regulation, ON-OFF-ON. Engineered signal processing filters can be directed to seven orthogonal promoters via adaptive modular DNA binding design. This collection of signal processing filters can be used in collaboration with our established transcriptional programming structure. Kinetic studies show that our collection of signal processing filters can switch between states of gene expression within a few minutes with minimal metabolic burden-representing a paradigm shift in general gene regulation.

Systems-level effects of allosteric perturbations to a model molecular switch.

Molecular switch proteins whose cycling between states is controlled by opposing regulators1,2 are central to biological signal transduction. As switch proteins function within highly connected interaction networks3, the fundamental question arises of how functional specificity is achieved when different processes share common regulators. Here we show that functional specificity of the small GTPase switch protein Gsp1 in Saccharomyces cerevisiae (the homologue of the human protein RAN)4 is linked to differential sensitivity of biological processes to different kinetics of the Gsp1 (RAN) switch cycle. We make 55 targeted point mutations to individual protein interaction interfaces of Gsp1 (RAN) and show through quantitative genetic5 and physical interaction mapping that Gsp1 (RAN) interface perturbations have widespread cellular consequences. Contrary to expectation, the cellular effects of the interface mutations group by their biophysical effects on kinetic parameters of the GTPase switch cycle and not by the targeted interfaces. Instead, we show that interface mutations allosterically tune the GTPase cycle kinetics. These results suggest a model in which protein partner binding, or post-translational modifications at distal sites, could act as allosteric regulators of GTPase switching. Similar mechanisms may underlie regulation by other GTPases, and other biological switches. Furthermore, our integrative platform to determine the quantitative consequences of molecular perturbations may help to explain the effects of disease mutations that target central molecular switches.