Prediction of Protein Allosteric Signalling Pathways and Functional Residues Through Paths of Optimised Propensity

Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. We showcase the approach with three well-studied allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were identified and showed agreement with experimental results, and pivotal signalling residues along the pathway were also revealed, thus providing alternative targets for drug design. By using the computed path scores, we were also able to differentiate the activity of different allosteric modulators.     

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.     

Novel non-active site inhibitor of Cryptosporidium hominis TS-DHFR identified by a virtual screen

The essential enzyme thymidylate synthase-dihydrofolate reductase (TS-DHFR) is a validated drug target for many pathogens, but has been elusive in Cryptosporidium hominis, as active site inhibitors of the enzymes from related parasitic protozoa show decreased potency and lack of species specificity over the human enzymes. As a rational approach to discover novel inhibitors, we conducted a virtual screen of a non-active site pocket in the DHFR linker region. From this screen, we have identified and characterized a noncompetitive inhibitor, flavin mononucleotide (FMN), with micromolar potency that is selective for ChTS-DHFR versus the human enzymes. These results describe a novel allosteric pocket amenable to inhibitor targeting, and a lead compound with which to move towards potent, selective inhibitors of ChTS-DHFR.     

Allosteric inhibitors of Mycobacterium tuberculosis tryptophan synthase

Global dispersion of multidrug resistant bacteria is very common and evolution of antibiotic‐resistance is occurring at an alarming rate, presenting a formidable challenge for humanity. The development of new therapeuthics with novel molecular targets is urgently needed. Current drugs primarily affect protein, nucleic acid, and cell wall synthesis. Metabolic pathways, including those involved in amino acid biosynthesis, have recently sparked interest in the drug discovery community as potential reservoirs of such novel targets. Tryptophan biosynthesis, utilized by bacteria but absent in humans, represents one of the currently studied processes with a therapeutic focus. It has been shown that tryptophan synthase (TrpAB) is required for survival of Mycobacterium tuberculosis in macrophages and for evading host defense, and therefore is a promising drug target. Here we present crystal structures of TrpAB with two allosteric inhibitors of M. tuberculosis tryptophan synthase that belong to sulfolane and indole‐5‐sulfonamide chemical scaffolds. We compare our results with previously reported structural and biochemical studies of another, azetidine‐containing M. tuberculosis tryptophan synthase inhibitor. This work shows how structurally distinct ligands can occupy the same allosteric site and make specific interactions. It also highlights the potential benefit of targeting more variable allosteric sites of important metabolic enzymes.     

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. Proteins 2015; 83:1375–1384. © 2014 Wiley Periodicals, Inc.     

新型冠状病毒相关靶点小分子抑制剂虚拟筛选

为确定治疗新型冠状病毒(SARS-CoV-2)感染的候选药物,开展了针对SARS-CoV-2的药物虚拟筛选研究。以SARS-CoV-2的刺突蛋白(S蛋白)和3CL蛋白酶(主蛋白酶)作为药物靶点,以美国食品药品监督管理局(FDA)批准上市的2 726个小分子药物作为候选,通过分子对接方法,筛选出了3种(阿巴瑞克(Abarelix)、西曲瑞克(Cetrorelix)、鞣酸(Tannic acid))与S蛋白具有较强结合能力的小分子药物,1种(戈舍瑞林(Goserelin))与3CL蛋白酶具有较好结合能力的小分子药物,它们理论上都具有抑制新型冠状病毒复制的效果。将靶向3CL蛋白酶的候选药物与辉瑞公司开发的药物Paxlovid进行比较,发现其作用位点均集中于3CL蛋白酶的第130-200位的残基周围,具有相似的结合位点与相互作用。此外也对候选药物的物理与化学性质及与基因相互作用进行了分析。本研究可为开发新型冠状病毒感染的治疗药物提供参考。     

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BackgroundAccumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.Methods/resultsWe employ two computationally efficient approaches – graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.ConclusionsContact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.General significanceThis study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.     

Discovery of Novel Allosteric Effectors Based on the Predicted Allosteric Sites for Escherichia coli D-3-Phosphoglycerate Dehydrogenase

D-3-phosphoglycerate dehydrogenase (PGDH) from Escherichia coli catalyzes the first critical step in serine biosynthesis, and can be allosterically inhibited by serine. In a previous study, we developed a computational method for allosteric site prediction using a coarse-grained two-state Gō Model and perturbation. Two potential allosteric sites were predicted for E. coli PGDH, one close to the active site and the nucleotide binding site (Site I) and the other near the regulatory domain (Site II). In the present study, we discovered allosteric inhibitors and activators based on site I, using a high-throughput virtual screen, and followed by using surface plasmon resonance (SPR) to eliminate false positives. Compounds 1 and 2 demonstrated a low-concentration activation and high-concentration inhibition phenomenon, with IC₅₀ values of 34.8 and 58.0 µM in enzymatic bioassays, respectively, comparable to that of the endogenous allosteric effector, L-serine. For its activation activity, compound 2 exhibited an AC₅₀ value of 34.7 nM. The novel allosteric site discovered in PGDH was L-serine- and substrate-independent. Enzyme kinetics studies showed that these compounds influenced K_m, k_cat, and k_cat/K_m. We have also performed structure-activity relationship studies to discover high potency allosteric effectors. Compound 2-2, an analog of compound 2, showed the best in vitro activity with an IC₅₀ of 22.3 µM. Compounds targeting this site can be used as new chemical probes to study metabolic regulation in E. coli. Our study not only identified a novel allosteric site and effectors for PGDH, but also provided a general strategy for designing new regulators for metabolic enzymes.     

A bifunctional allosteric site in the dimer interface of procaspase-3

The dimer interface of caspase-3 contains a bifunctional allosteric site in which the enzyme can be activated or inactivated, depending on the context of the protein. In the mature caspase-3, the binding of allosteric inhibitors to the interface results in an order-to-disorder transition in the active site loops. In procaspase-3, by contrast, the binding of allosteric activators to the interface results in a disorder-to-order transition in the active site. We have utilized the allosteric site to identify a small molecule activator of procaspase and to characterize its binding to the protease. The data suggest that an efficient activator must stabilize the active conformer of the zymogen by expelling the intersubunit linker from the interface, and it must interact with active site residues found in the allosteric site. Small molecule activators that fulfill the two requirements should provide scaffolds for drug candidates as a therapeutic strategy for directly promoting procaspase-3 activation in cancer cells.     

Distinct Roles of N-Glycosylation at Different Sites of Corin in Cell Membrane Targeting and Ectodomain Shedding

Corin is a membrane-bound protease essential for activating natriuretic peptides and regulating blood pressure. Human corin has 19 predicted N-glycosylation sites in its extracellular domains. It has been shown that N-glycans are required for corin cell surface expression and zymogen activation. It remains unknown, however, how N-glycans at different sites may regulate corin biosynthesis and processing. In this study, we examined corin mutants, in which each of the 19 predicted N-glycosylation sites was mutated individually. By Western analysis of corin proteins in cell lysate and conditioned medium from transfected HEK293 cells and HL-1 cardiomyocytes, we found that N-glycosylation at Asn-80 inhibited corin shedding in the juxtamembrane domain. Similarly, N-glycosylation at Asn-231 protected corin from autocleavage in the frizzled-1 domain. Moreover, N-glycosylation at Asn-697 in the scavenger receptor domain and at Asn-1022 in the protease domain is important for corin cell surface targeting and zymogen activation. We also found that the location of the N-glycosylation site in the protease domain was not critical. N-Glycosylation at Asn-1022 may be switched to different sites to promote corin zymogen activation. Together, our results show that N-glycans at different sites may play distinct roles in regulating the cell membrane targeting, zymogen activation, and ectodomain shedding of corin.     

Engineering protein allostery: 1.05 A resolution structure and enzymatic properties of a Na+-activated trypsin

Some trypsin-like proteases are endowed with Na⁺-dependent allosteric enhancement of catalytic activity, but this important mechanism has been difficult to engineer in other members of the family. Replacement of 19 amino acids in Streptomyces griseus trypsin targeting the active site and the Na⁺-binding site were found necessary to generate efficient Na⁺ activation. Remarkably, this property was linked to the acquisition of a new substrate selectivity profile similar to that of factor Xa, a Na⁺-activated protease involved in blood coagulation. The X-ray crystal structure of the mutant trypsin solved to 1.05 Å resolution defines the engineered Na⁺ site and active site loops in unprecedented detail. The results demonstrate that trypsin can be engineered into an efficient allosteric protease, and that Na⁺ activation is interwoven with substrate selectivity in the trypsin scaffold.     

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

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

A structurally preserved allosteric site in the MIF superfamily affects enzymatic activity and CD74 activation in D-dopachrome tautomerase

The macrophage migration inhibitory factor (MIF) family of cytokines contains multiple ligand-binding sites and mediates immunomodulatory processes through an undefined mechanism(s). Previously, we reported a dynamic relay connecting the MIF catalytic site to an allosteric site at its solvent channel. Despite structural and functional similarity, the MIF homolog D-dopachrome tautomerase (also called MIF-2) has low sequence identity (35%), prompting the question of whether this dynamic regulatory network is conserved. Here, we establish the structural basis of an allosteric site in MIF-2, showing with solution NMR that dynamic communication is preserved in MIF-2 despite differences in the primary sequence. X-ray crystallography and NMR detail the structural consequences of perturbing residues in this pathway, which include conformational changes surrounding the allosteric site, despite global preservation of the MIF-2 fold. Molecular simulations reveal MIF-2 to contain a comparable hydrogen bond network to that of MIF, which was previously hypothesized to influence catalytic activity by modulating the strength of allosteric coupling. Disruption of the allosteric relay by mutagenesis also attenuates MIF-2 enzymatic activity in vitro and the activation of the cluster of differentiation 74 receptor in vivo, highlighting a conserved point of control for nonoverlapping functions in the MIF superfamily.     

Enzyme–substrate complexes of allosteric citrate synthase: Evidence for a novel intermediate in substrate binding

The citrate synthase (CS) of Escherichia coli is an allosteric hexameric enzyme specifically inhibited by NADH. The crystal structure of wild type (WT) E. coli CS, determined by us previously, has no substrates bound, and part of the active site is in a highly mobile region that is shifted from the position needed for catalysis. The CS of Acetobacter aceti has a similar structure, but has been successfully crystallized with bound substrates: both oxaloacetic acid (OAA) and an analog of acetyl coenzyme A (AcCoA). We engineered a variant of E. coli CS wherein five amino acids in the mobile region have been replaced by those in the A. aceti sequence. The purified enzyme shows unusual kinetics with a low affinity for both substrates. Although the crystal structure without ligands is very similar to that of the WT enzyme (except in the mutated region), complexes are formed with both substrates and the allosteric inhibitor NADH. The complex with OAA in the active site identifies a novel OAA-binding residue, Arg306, which has no functional counterpart in other known CS-OAA complexes. This structure may represent an intermediate in a multi-step substrate binding process where Arg306 changes roles from OAA binding to AcCoA binding. The second complex has the substrate analog, S-carboxymethyl-coenzyme A, in the allosteric NADH-binding site and the AcCoA site is not formed. Additional CS variants unable to bind adenylates at the allosteric site show that this second complex is not a factor in positive allosteric activation of AcCoA binding.     

Learning About Allosteric Drugs and Ways to Design Them

While the accelerating quest for precision medicine requires new individually targeting and selective drugs, and the ability to work with so-called undruggable targets, the realm of allosteric drugs meeting this need remains largely uncharted. Generalizing the observations on two major drug targets with widely observed inherent allostery, GPCRs and kinases, we describe and discuss basic allosteric modes of action that are universally applicable in all types of structures and functions. Using examples of Class A GPCRs and CMGC protein kinases, we show how Allosteric Signalling and Probing Fingerprints can be used to identify potential allosteric sites and reveal effector-leads that may serve as a starting point for the development of allosteric drugs targeting these regulatory sites. A set of distinct characteristics of allosteric ligands was established, which highlights the versatility of their design and make them advantageous before their orthosteric counterparts in personalized medicine. We argue that rational design of allosteric drugs should begin with the search for latent sites or design of non-natural binding sites followed by fragment-based design of allosteric ligands and by the mutual adjustment of the site-ligand pair in order to achieve required drug efficacy. On the basis of the perturbative nature and reversibility of allosteric communication, we propose a generic protocol for computational design of allosteric effectors, enabling also the allosteric tuning of biologics, in obtaining allosteric control over protein functions.     

Analysis of binding site hot spots on the surface of Ras GTPase

We have recently discovered an allosteric switch in Ras, bringing an additional level of complexity to this GTPase whose mutants are involved in nearly 30% of cancers. Upon activation of the allosteric switch, there is a shift in helix 3/loop 7 associated with a disorder to order transition in the active site. Here, we use a combination of multiple solvent crystal structures and computational solvent mapping (FTMap) to determine binding site hot spots in the “off” and “on” allosteric states of the GTP-bound form of H-Ras. Thirteen sites are revealed, expanding possible target sites for ligand binding well beyond the active site. Comparison of FTMaps for the H and K isoforms reveals essentially identical hot spots. Furthermore, using NMR measurements of spin relaxation, we determined that K-Ras exhibits global conformational dynamics very similar to those we previously reported for H-Ras. We thus hypothesize that the global conformational rearrangement serves as a mechanism for allosteric coupling between the effector interface and remote hot spots in all Ras isoforms. At least with respect to the binding sites involving the G domain, H-Ras is an excellent model for K-Ras and probably N-Ras as well. Ras has so far been elusive as a target for drug design. The present work identifies various unexplored hot spots throughout the entire surface of Ras, extending the focus from the disordered active site to well-ordered locations that should be easier to target.     

Novel Allosteric Modulators of G Protein-coupled Receptors

G protein-coupled receptors (GPCRs) are allosteric proteins, because their signal transduction relies on interactions between topographically distinct, yet conformationally linked, domains. Much of the focus on GPCR allostery in the new millennium, however, has been on modes of targeting GPCR allosteric sites with chemical probes due to the potential for novel therapeutics. It is now apparent that some GPCRs possess more than one targetable allosteric site, in addition to a growing list of putative endogenous modulators. Advances in structural biology are also shedding new insights into mechanisms of allostery, although the complexities of candidate allosteric drugs necessitate rigorous biological characterization.     

Probing the Activity Modification Space of the Cysteine Peptidase Cathepsin K with Novel Allosteric Modifiers

Targeting allosteric sites is gaining increasing recognition as a strategy for modulating the activity of enzymes, especially in drug design. Here we investigate the mechanisms of allosteric regulation of cathepsin K as a representative of cysteine cathepsins and a promising drug target for the treatment of osteoporosis. Eight novel modifiers are identified by computational targeting of predicted allosteric sites on the surface of the enzyme. All act via hyperbolic kinetic mechanisms in presence of low molecular mass substrates, as expected for allosteric effectors. Two compounds have sizable effects on enzyme activity using interstitial collagen as a natural substrate of cathepsin K and four compounds show a significantly stabilizing effect on cathepsin K. The concept of activity modification space is introduced to obtain a global perspective of the effects elicited by the modifiers. Analysis of the activity modification space reveals that the activity of cathepsin K is regulated via multiple, different allosteric mechanisms.     

Tricyclic GyrB/ParE (TriBE) Inhibitors: A New Class of Broad-Spectrum Dual-Targeting Antibacterial Agents

Increasing resistance to every major class of antibiotics and a dearth of novel classes of antibacterial agents in development pipelines has created a dwindling reservoir of treatment options for serious bacterial infections. The bacterial type IIA topoisomerases, DNA gyrase and topoisomerase IV, are validated antibacterial drug targets with multiple prospective drug binding sites, including the catalytic site targeted by the fluoroquinolone antibiotics. However, growing resistance to fluoroquinolones, frequently mediated by mutations in the drug-binding site, is increasingly limiting the utility of this antibiotic class, prompting the search for other inhibitor classes that target different sites on the topoisomerase complexes. The highly conserved ATP-binding subunits of DNA gyrase (GyrB) and topoisomerase IV (ParE) have long been recognized as excellent candidates for the development of dual-targeting antibacterial agents with broad-spectrum potential. However, to date, no natural product or small molecule inhibitors targeting these sites have succeeded in the clinic, and no inhibitors of these enzymes have yet been reported with broad-spectrum antibacterial activity encompassing the majority of Gram-negative pathogens. Using structure-based drug design (SBDD), we have created a novel dual-targeting pyrimidoindole inhibitor series with exquisite potency against GyrB and ParE enzymes from a broad range of clinically important pathogens. Inhibitors from this series demonstrate potent, broad-spectrum antibacterial activity against Gram-positive and Gram-negative pathogens of clinical importance, including fluoroquinolone resistant and multidrug resistant strains. Lead compounds have been discovered with clinical potential; they are well tolerated in animals, and efficacious in Gram-negative infection models.     

The flavivirus protease as a target for drug discovery

Many flaviviruses are significant human pathogens causing considerable disease burdens, including encephalitis and hemorrhagic fever, in the regions in which they are endemic. A paucity of treatments for flaviviral infections has driven interest in drug development targeting proteins essential to flavivirus replication, such as the viral protease. During viral replication, the flavivirus genome is translated as a single polyprotein precursor, which must be cleaved into individual proteins by a complex of the viral protease, NS3, and its cofactor, NS2B. Because this cleavage is an obligate step of the viral life-cycle, the flavivirus protease is an attractive target for antiviral drug development. In this review, we will survey recent drug development studies targeting the NS3 active site, as well as studies targeting an NS2B/NS3 interaction site determined from flavivirus protease crystal structures.