Designed covalent allosteric modulators: an emerging paradigm in drug discovery

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Landscape of In Silico Tools for Modeling Covalent Modification of Proteins: A Review on Computational Covalent Drug Discovery

Covalent drug discovery has been a challenging research area given the struggle of finding a sweet balance between selectivity and reactivity for these drugs, the lack of which often leads to off-target activities and hence undesirable side effects. However, there has been a resurgence in covalent drug design following the success of several covalent drugs such as boceprevir (2011), ibrutinib (2013), neratinib (2017), dacomitinib (2018), zanubrutinib (2019), and many others. Design of covalent drugs includes many crucial factors, where "evaluation of the binding affinity" and "a detailed mechanistic understanding on covalent inhibition" are at the top of the list. Well-defined experimental techniques are available to elucidate these factors; however, often they are expensive and/or time-consuming and hence not suitable for high throughput screens. Recent developments in in silico methods provide promise in this direction. In this report, we review a set of recent publications that focused on developing and/or implementing novel in silico techniques in "Computational Covalent Drug Discovery (CCDD)". We also discuss the advantages and disadvantages of these approaches along with what improvements are required to make it a great tool in medicinal chemistry in the near future.

Chemoproteomics Reveals the Pan-HER Kinase Inhibitor Neratinib To Target an Arabidopsis Epoxide Hydrolase Related to Phytohormone Signaling

Plant phytohormone pathways are regulated by an intricate network of signaling components and modulators, many of which still remain unknown. Here, we report a forward chemical genetics approach for the identification of functional SA agonists in Arabidopsis thaliana that revealed Neratinib (Ner), a covalent pan-HER kinase inhibitor drug in humans, as a modulator of SA signaling. Instead of a protein kinase, chemoproteomics unveiled that Ner covalently modifies a surface-exposed cysteine residue of Arabidopsis epoxide hydrolase isoform 7 (AtEH7), thereby triggering its allosteric inhibition. Physiologically, the Ner application induces jasmonate metabolism in an AtEH7-dependent manner as an early response. In addition, it modulates PATHOGENESIS RELATED 1 (PR1) expression as a hallmark of SA signaling activation as a later effect. AtEH7, however, is not the exclusive target for this physiological readout induced by Ner. Although the underlying molecular mechanisms of AtEH7-dependent modulation of jasmonate signaling and Ner-induced PR1-dependent activation of SA signaling and thus defense response regulation remain unknown, our present work illustrates the powerful combination of forward chemical genetics and chemical proteomics for identifying novel phytohormone signaling modulatory factors. It also suggests that marginally explored metabolic enzymes such as epoxide hydrolases may have further physiological roles in modulating signaling.

Strategies for Conditional Regulation of Proteins

Design of the next-generation of therapeutics, biosensors, and molecular tools for basic research requires that we bring protein activity under control. Each protein has unique properties, and therefore, it is critical to tailor the current techniques to develop new regulatory methods and regulate new proteins of interest (POIs). This perspective gives an overview of the widely used stimuli and synthetic and natural methods for conditional regulation of proteins.

Mechanistic insights into the clinical Y96D mutation with acquired resistance to AMG510 in the KRASG12C

Special oncogenic mutations in the RAS proteins lead to the aberrant activation of RAS and its downstream signaling pathways. AMG510, the first approval drug for KRAS, covalently binds to the mutated cysteine 12 of KRASG12C protein and has shown promising antitumor activity in clinical trials. Recent studies have reported that the clinically acquired Y96D mutation could severely affect the effectiveness of AMG510. However, the underlying mechanism of the drug-resistance remains unclear. To address this, we performed multiple microsecond molecular dynamics simulations on the KRASG12C−AMG510 and KRASG12C/Y96D−AMG510 complexes at the atomic level. The direct interaction between the residue 96 and AMG510 was impaired owing to the Y96D mutation. Moreover, the mutation yielded higher flexibility and more coupled motion of the switch II and α3-helix, which led to the departing motion of the switch II and α3-helix. The resulting departing motion impaired the interaction between the switch II and α3-helix and subsequently induced the opening and loosening of the AMG510 binding pocket, which further disrupted the interaction between the key residues in the pocket and AMG510 and induced an increased solvent exposure of AMG510. These findings reveal the resistance mechanism of AMG510 to KRASG12C/Y96D, which will help to offer guidance for the development of KRAS targeted drugs to overcome acquired resistance.

A Comprehensive Guide for Assessing Covalent Inhibition in Enzymatic Assays Illustrated with Kinetic Simulations

Covalent inhibition has become more accepted in the past two decades, as illustrated by the clinical approval of several irreversible inhibitors designed to covalently modify their target. Elucidation of the structure-activity relationship and potency of such inhibitors requires a detailed kinetic evaluation. Here, we elucidate the relationship between the experimental read-out and the underlying inhibitor binding kinetics. Interactive kinetic simulation scripts are employed to highlight the effects of in vitro enzyme activity assay conditions and inhibitor binding mode, thereby showcasing which assumptions and corrections are crucial. Four stepwise protocols to assess the biochemical potency of (ir)reversible covalent enzyme inhibitors targeting a nucleophilic active site residue are included, with accompanying data analysis tailored to the covalent binding mode. Together, this will serve as a guide to make an educated decision regarding the most suitable method to assess covalent inhibition potency. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol I: Progress curve analysis of substrate association competition Basic Data Analysis Protocol 1A: Two-step irreversible covalent inhibition Basic Data Analysis Protocol 1B: One-step irreversible covalent inhibition Basic Data Analysis Protocol 1C: Two-step reversible covalent inhibition Basic Data Analysis Protocol 1D: Two-step irreversible covalent inhibition with substrate depletion Basic Protocol II: Incubation time-dependent potency IC₅₀ (t) Basic Data Analysis Protocol 2: Two-step irreversible covalent inhibition Basic Protocol III: Preincubation time-dependent inhibition without dilution Basic Data Analysis Protocol 3: Preincubation time-dependent inhibition without dilution Basic Data Analysis Protocol 3Ai: Two-step irreversible covalent inhibition Alternative Data Analysis Protocol 3Aii: Two-step irreversible covalent inhibition Basic Data Analysis Protocol 3Bi: One-step irreversible covalent inhibition Alternative Data Analysis Protocol 3Bii: One-step irreversible covalent inhibition Basic Data Analysis Protocol 3C: Two-step reversible covalent inhibition Basic Protocol IV: Preincubation time-dependent inhibition with dilution/competition Basic Data Analysis Protocol 4: Preincubation time-dependent inhibition with dilution Basic Data Analysis Protocol 4Ai: Two-step irreversible covalent inhibition Alternative Data Analysis Protocol 4Aii: Two-step irreversible covalent inhibition Basic Data Analysis Protocol 4Bi: One-step irreversible covalent inhibition Alternative Data Analysis Protocol 4Bii: One-step irreversible covalent inhibition.

Understanding the P-Loop Conformation in the Determination of Inhibitor Selectivity Toward the Hepatocellular Carcinoma-Associated Dark Kinase STK17B

As a member of the death-associated protein kinase family of serine/threonine kinases, the STK17B has been associated with diverse diseases such as hepatocellular carcinoma. However, the conformational dynamics of the phosphate-binding loop (P-loop) in the determination of inhibitor selectivity profile to the STK17B are less understood. Here, a multi-microsecond length molecular dynamics (MD) simulation of STK17B in the three different states (ligand-free, ADP-bound, and ligand-bound states) was carried out to uncover the conformational plasticity of the P-loop. Together with the analyses of principal component analysis, cross-correlation and generalized correlation motions, secondary structural analysis, and community network analysis, the conformational dynamics of the P-loop in the different states were revealed, in which the P-loop flipped into the ADP-binding site upon the inhibitor binding and interacted with the inhibitor and the C-lobe, strengthened the communication between the N- and C-lobes. These resulting interactions contributed to inhibitor selectivity profile to the STK17B. Our results may advance our understanding of kinase inhibitor selectivity and offer possible implications for the design of highly selective inhibitors for other protein kinases.

The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs

G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.

Selection Method of DNA-Encoded Chemical Library for Irreversible Covalent Binders

Compounds acting through irreversible covalent interactions with therapeutic targets represent an important mechanism of drug action. Here I describe a selection method for DNA-encoded libraries to discover irreversible covalent binders to target proteins. This method offers an enabling tool in drug discovery for therapeutic targets that may be undruggable for reversible inhibition.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.

Myeloid cell leukemia 1 (MCL-1): Structural characteristics and application in cancer therapy

Apoptosis, a major hallmark of cancer cells, regulates cellular fate and homeostasis. BCL-2 (B-cell CLL/Lymphoma 2) protein family is popularly known to mediate the intrinsic mode of apoptosis, of which MCL-1 is a crucial member. Myeloid cell leukemia 1 (MCL-1) is an anti-apoptotic oncoprotein and one of the most investigated members of the BCL-2 family. It is commonly known to be genetically altered, aberrantly overexpressed, and primarily associated with drug resistance in various human cancers. Recent advancements in the development of selective MCL-1 inhibitors and evaluating their effectiveness in cancer treatment establish its popularity as a molecular target. The overall aim is the selective induction of apoptosis in cancer cells by using a single or combination of BCL-2 family inhibitors. Delineating the precise molecular mechanisms associated with MCL-1-mediated cancer progression will certainly improve the efficacy of clinical interventions aimed at MCL-1 and hence patient survival. This review is structured to highlight the structural characteristics of MCL-1, its specific interactions with NOXA, MCL-1-regulatory microRNAs, and at the same time focus on the emerging therapeutic strategies targeting our protein of interest (MCL-1), alone or in combination with other treatments.

Rosmarinic acid prevents refractory bacterial pneumonia through regulating Keap1/Nrf2-mediated autophagic pathway and mitochondrial oxidative stress

Methicillin-resistant Staphylococcus aureus (MRSA) is the leading cause of bacterial pneumonia, featured with exuberant inflammatory cytokine production, extensive oxidative stress and tissue injury. The Keap1/Nrf2 system is the major apparatus essential for host defense against oxidative and electrophilic stresses of both exogenous and endogenous origins, representing a logical target for host-directed strategy to treat severe inflammatory diseases including MRSA-induced pneumonia. In an effort to search therapeutics for bacterial pneumonia, we identify rosmarinic acid (RA) as a covalent modifier of Keap1 and hence an activator of Nrf2. Specifically, RA forms a covalent bond with the cysteine 151 of Keap1 in BTB domain, and blocks its association with Nrf2 for proteasome-mediated degradation. Consequently, RA treatment caused the increased Nrf2 nuclear translocation to initiate antioxidant and mitochondrial biogenic programs, as well as macrophage bactericidal activity through inducing autophagic pathway, which eventually led to expedited bacterial eradication, inflammation resolution, and disease recovery. Collectively, our findings establish RA as a specific inducer of Nrf2 and show its potential to prevent MRSA pneumonia.

Drug design targeting active posttranslational modification protein isoforms

Modern drug design aims to discover novel lead compounds with attractable chemical profiles to enable further exploration of the intersection of chemical space and biological space. Identification of small molecules with good ligand efficiency, high activity, and selectivity is crucial toward developing effective and safe drugs. However, the intersection is one of the most challenging tasks in the pharmaceutical industry, as chemical space is almost infinity and continuous, whereas the biological space is very limited and discrete. This bottleneck potentially limits the discovery of molecules with desirable properties for lead optimization. Herein, we present a new direction leveraging posttranslational modification (PTM) protein isoforms target space to inspire drug design termed as "Post-translational Modification Inspired Drug Design (PTMI-DD)." PTMI-DD aims to extend the intersections of chemical space and biological space. We further rationalized and highlighted the importance of PTM protein isoforms and their roles in various diseases and biological functions. We then laid out a few directions to elaborate the PTMI-DD in drug design including discovering covalent binding inhibitors mimicking PTMs, targeting PTM protein isoforms with distinctive binding sites from that of wild-type counterpart, targeting protein-protein interactions involving PTMs, and hijacking protein degeneration by ubiquitination for PTM protein isoforms. These directions will lead to a significant expansion of the biological space and/or increase the tractability of compounds, primarily due to precisely targeting PTM protein isoforms or complexes which are highly relevant to biological functions. Importantly, this new avenue will further enrich the personalized treatment opportunity through precision medicine targeting PTM isoforms.

Covalent allosteric modulation: An emerging strategy for GPCRs drug discovery

Designing covalent allosteric modulators brings new opportunities to the field of drug discovery towards G-protein-coupled receptors (GPCRs). Targeting an allosteric binding pocket can allow a modulator to have protein subtype selectivity and low drug resistance. Utilizing covalent warheads further enables the modulator to increase the binding potency and extend the duration of action. This review starts with GPCR allosteric modulation to discuss the structural biology of allosteric binding pockets, the different types of allosteric modulators, as well as the advantages of employing allosteric modulation. This is followed by a discussion on covalent modulators to clarify how covalent ligands can benefit the receptor modulation and to illustrate moieties that can commonly be used as covalent warheads. Finally, case studies are presented on designing class A, B, and C GPCR covalent allosteric modulators to demonstrate successful stories on combining allosteric modulation and covalent binding. Limitations and future perspectives are also covered.

Impact of Allosteric Modulation in Drug Discovery: Innovation in Emerging Chemical Modalities

Recent years have seen an unprecedented level of innovation in allosteric drug discovery and development, with multiple drug candidates advancing into clinical studies. From early examples of allosteric drugs like GABA_A receptor modulators (benzodiazepines) in the 1960s to more recent GPCR negative allosteric modulators of CCR5 (maraviroc) approved in 2007, the opportunities for interrogating allosteric sites in drug discovery have expanded to other target classes such as protein-protein interactions, kinases, and nuclear hormone receptors. In this Innovation Letter, the authors highlight the latest advances of allosteric drug discovery from different target classes and novel emerging chemical modalities beyond small molecules.

Identification of the New Covalent Allosteric Binding Site of Fructose-1,6-bisphosphatase with Disulfiram Derivatives toward Glucose Reduction

Fructose 1,6-bisphosphatase (FBPase) has attracted substantial interest as a target associated with cancer and type 2 diabetes. Herein, we found that disulfiram and its derivatives can potently inhibit FBPase by covalently binding to a new C128 allosteric site distinct from the original C128 site in APO FBPase. Further identification of the allosteric inhibition mechanism reveals that the covalent binding of a fragment of 214 will result in the movement of C128 and the dissociation of helix H4 (123-128), which in turn allows S123 to more easily form new hydrogen bonds with K71 and D74 in helix H3 (69-72), thereby inhibiting FBPase activity. Notably, both disulfiram and 212 might moderately reduce blood glucose output in vivo. Therefore, our current findings not only identify a new covalent allosteric site of FBPase but also establish a structural foundation and provide a promising way for the design of covalent allosteric drugs for glucose reduction.

Thiazolidine-2,4-dione-based irreversible allosteric IKK-β kinase inhibitors: Optimization into in vivo active anti-inflammatory agents

Selective kinase inhibitors development is a cumbersome task because of ATP binding sites similarities across kinases. On contrast, irreversible allosteric covalent inhibition offers opportunity to develop novel selective kinase inhibitors. Previously, we reported thiazolidine-2,4-dione lead compounds eliciting in vitro irreversible allosteric inhibition of IKK-β. Herein, we address optimization into in vivo active anti-inflammatory agents. We successfully developed potent IKK-β inhibitors with the most potent compound eliciting IC₅₀ = 0.20 μM. Cellular assay of a set of active compounds using bacterial endotoxin lipopolysaccharide (LPS)-stimulated macrophages elucidated significant in vitro anti-inflammatory activity. In vitro evaluation of microsomal and plasma stabilities showed that the promising compound 7a is more stable than compound 7p. Finally, in vivo evaluation of 7a, which has been conducted in a model of LPS-induced septic shock in mice, showed its ability to protect mice against septic shock induced mortality. Accordingly, this study presents compound 7a as a novel potential irreversible allosteric covalent inhibitor of IKK-β with verified in vitro and in vivo anti-inflammatory activity.

The chemical diversity and structure-based discovery of allosteric modulators for the PIF-pocket of protein kinase PDK1

Phosphoinositide-dependent protein kinase-1 (PDK1) is an important protein in mediating the PI3K-AKT pathway and is thus identified as a promising target. The catalytic activity of PDK1 is tightly regulated by allosteric modulators, which bind to the PDK1 Interacting Fragment (PIF) pocket of the kinase domain that is topographically distinct from the orthosteric, ATP binding site. Allosteric modulators by attaching to the less conserved PIF-pocket have remarkable advantages such as higher selectivity, less side effect, and lower toxicity. Targeting allosteric PIF-pocket of PDK1 has become the focus of recent attention. In this review, we summarise the current advances in the structure-based discovery of PDK1 allosteric modulators. We will first present the three-dimensional structure of PDK1 and illustrate the allosteric regulatory mechanism of PDK1 through the modulation of the PIF-pocket. Then, the recent advances of PDK1 allosteric modulators targeting the PIF-pocket will be recapitulated detailly according to the structural similarity of allosteric modulators.

Does Ras Activate Raf and PI3K Allosterically?

The mechanism through which oncogenic Ras activates its effectors is vastly important to resolve. If allostery is at play, then targeting allosteric pathways could help in quelling activation of MAPK (Raf/MEK/ERK) and PI3K (PI3K/Akt/mTOR) cell proliferation pathways. On the face of it, allosteric activation is reasonable: Ras binding perturbs the conformational ensembles of its effectors. Here, however, we suggest that at least for Raf, PI3K, and NORE1A (RASSF5), that is unlikely. Raf's long disordered linker dampens effective allosteric activation. Instead, we suggest that the high-affinity Ras–Raf binding relieves Raf's autoinhibition, shifting Raf's ensemble from the inactive to the nanocluster-mediated dimerized active state, as Ras also does for NORE1A. PI3K is recruited and allosterically activated by RTK (e.g., EGFR) at the membrane. Ras restrains PI3K's distribution and active site orientation. It stabilizes and facilitates PIP₂ binding at the active site and increases the PI3K residence time at the membrane. Thus, RTKs allosterically activate PI3Kα; however, merging their action with Ras accomplishes full activation. Here we review their activation mechanisms in this light and draw attention to implications for their pharmacology.

Emergence of allosteric drug-resistance mutations: new challenges for allosteric drug discovery

Allosteric drugs have several significant advantages over traditional orthosteric drugs, encompassing higher selectivity and lower toxicity. Although allosteric drugs have potential advantages as therapeutic agents to treat human diseases, allosteric drug-resistance mutations still occur, rendering these drugs ineffective. Here, we review the emergence of allosteric drug-resistance mutations with an emphasis on examples covering clinically important therapeutic targets, including Breakpoint cluster region-Abelson tyrosine kinase (Bcr-Abl), Akt kinase [also called Protein Kinase B (PKB)], isocitrate dehydrogenase (IDH), MAPK/ERK kinase (MEK), and SRC homology 2 domain-containing phosphatase 2 (SHP2). We also discuss challenges associated with tackling allosteric drug resistance and the possible strategies to overcome this issue.

Optimization study towards more potent thiazolidine-2,4-dione IKK-β modulator: Synthesis, biological evaluation and in silico docking simulation

Inhibition of IKK-β (inhibitor of nuclear factor kappa-B kinase subunit beta) has been broadly documentedas a promising approach for treatment of acute and chronic inflammatory diseases, cancer, and autoimmune diseases. Recently, we have identified a novel class of thiazolidine-2,4-diones as structurally novel modulators for IKK-β. Herein, we report a hit optimization study via analog synthesis strategy aiming to acquire more potent derivative(s), probe the structure activity relationship (SAR), and get reasonable explanations for the elicited IKK-β inhibitory activities though an in silico docking simulation study. Accordingly, a new series of eighteen thiazolidine-2,4-dione derivatives was rationally designed, synthesized, identified with different spectroscopic techniques and biologically evaluated as noteworthy IKK-β potential modulators. Successfully, new IKK-β potent modulators were obtained, including the most potent analog up-to-date 7m with IC₅₀ value of 260 nM. A detailed structure activity relationship (SAR) was discussed and a mechanistic study for 7m was carried out indicating its irreversible inhibition mode with IKK-β (K_inact value = 0.01 (min^-1). Furthermore, the conducted in silico simulation study provided new insights for the binding modes of this novel class of modulators with IKK-β.

Drugging K-RasG12C through covalent inhibitors: Mission possible?

Ras, whose mutants are present in approximately 30% of human tumours, is one of the most important oncogenes. Drugging Ras is thus regarded as the quest for the Holy Grail in cancer therapeutics development. Despite more than three decades of efforts, drug discovery targeting Ras constantly fails, rendering Ras undruggable, due to its smooth surface and picomolar affinity towards guanosine substrates. The most frequently mutated isoform of Ras is K-Ras, accounting for >85% of Ras-driven cancers, and one majority of them is the G12C mutation. Recent advances in structural biology shed light on drugging Ras, and one of the cutting-edge breakthroughs is the design of covalent G12C-specific inhibitors targeting the mutated cysteine. This type of inhibitor can be classified into substrate-competitive orthosteric inhibitors and non-competitive allosteric inhibitors. They display improved selectivity and enhanced potency due to their G12-specific and irreversible covalent binding nature. Thus, they represent a new hope for revolutionizing the conventional characterization of Ras as "undruggable" and pave a promising avenue for further drug discovery. Here, we provide comprehensive structural and medicinal chemical insights into K-Ras covalent inhibitors specific for the G12C mutant. We first present an in-depth analysis of the conformations of the inhibitor binding pockets. Then, all the latest covalent ligands selectively inhibiting K-Ras^G12C are reviewed. Finally, we examine the current challenges faced by this new class of anti-Ras inhibitors.

The structural basis of the autoinhibition mechanism of glycogen synthase kinase 3β (GSK3β): molecular modeling and molecular dynamics simulation studies

The autoinhibition phenomenon has been frequently observed in enzymes and represents an important regulatory strategy to fine-tune enzyme activity. Evolution has exploited this mechanism to reduce enzymatic activity. Glycogen synthase kinase 3β (GSK3β) undergoes autoinhibition via the phosphorylation of Ser9 at the N-terminus of the kinase, which, acting as a pseudosubstrate, occupies the catalytic domain of GSK3β and subsequently blocks primed substrates from having access to the catalytic domain. The detailed structural basis of the autoinhibition mechanism of GSK3β by the pseudosubstrate, however, has not yet been fully resolved. Here, a three-dimensional model of the binary GSK3β-pseudosubstrate complex was built via the molecular modeling method. Based on the constructed model, extensive molecular dynamics (MD) simulations and subsequent molecular mechanics generalized Born/surface area (MM_GBSA) calculations were performed on the wild-type GSK3β-pseudosubstrate complex and three mutated systems (R4A, R6A, and S9A). Analyses of MD simulations and binding free energies revealed that the phosphorylation of Ser9 is the prerequisite for the autoinhibition of GSK3β, and both mutations of Arg4 and Arg6 to alanine markedly reduced the binding affinities of the mutated pseudosubstrate to the GSK3β catalytic domain, thereby disrupting the autoinhibition of the kinase. This study highlights the importance of Ser9, Arg6, and Arg4 in modulating the autoinhibition mechanism of GSK3β, contributing to a deeper understanding of GSK3β biology.Communicated by Ramaswamy H. Sarma.

Emerging roles of allosteric modulators in the regulation of protein-protein interactions (PPIs): A new paradigm for PPI drug discovery

Protein-protein interactions (PPIs) are closely implicated in various types of cellular activities and are thus pivotal to health and disease states. Given their fundamental roles in a wide range of biological processes, the modulation of PPIs has enormous potential in drug discovery. However, owing to the general properties of large, flat, and featureless interfaces of PPIs, previous attempts have demonstrated that the generation of therapeutic agents targeting PPI interfaces is challenging, rendering them almost "undruggable" for decades. To date, rapid progress in chemical and structural biology techniques has promoted the exploitation of allostery as a novel approach in drug discovery. By attaching to allosteric sites that are topologically and spatially distinct from PPI interfaces, allosteric modulators can achieve improved physiochemical properties. Thus, allosteric modulators may represent an alternative strategy to target intractable PPIs and have attracted intense pharmaceutical interest. In this review, we first briefly introduce the characteristics of PPIs and then present different approaches for investigating PPIs, as well as the latest methods for modulating PPIs. Importantly, we comprehensively review the recent progress in the development of allosteric modulators to inhibit or stabilize PPIs. Finally, we conclude with future perspectives on the discovery of allosteric PPI modulators, especially the application of computational methods to aid in allosteric PPI drug discovery.

Allosteric Modulation of Ionotropic Glutamate Receptors: An Outlook on New Therapeutic Approaches To Treat Central Nervous System Disorders

The allosteric targeting of ionotropic glutamate receptors (iGluRs) is a valuable approach for treating various central nervous system (CNS) disorders. In this frame, this Innovations provides a summary of the state-of-the art in the development of allosteric modulators for iGluRs and offers an outlook regarding innovative strategies for treating neurological diseases.

Allosteric Modalities for Membrane-Bound Receptors: Insights from Drug Hunting for Brain Diseases

Medicinal chemists are accountable for embedding the appropriate drug target profile into the molecular architecture of a clinical candidate. An accurate characterization of the functional effects following binding of a drug to its biological target is a fundamental step in the discovery of new medicines, informing the translation of preclinical efficacy and safety observations into human trials. Membrane-bound proteins, particularly ion channels and G protein-coupled receptors (GPCRs), are biological targets prone to allosteric modulation. Investigations using allosteric drug candidates and chemical tools suggest that their functional effects may be tailored with a high degree of translational alignment, making them molecular tools to correct pathophysiological functional tone and enable personalized medicine when a causative target-to-disease link is known. We present select examples of functional molecular fine-tuning of allosterism and discuss consequences relevant to drug design.

Allosteric Methods and Their Applications: Facilitating the Discovery of Allosteric Drugs and the Investigation of Allosteric Mechanisms

Allostery, or allosteric regulation, is the phenomenon in which protein functional activity is altered by the binding of an effector at an allosteric site that is topographically distinct from the orthosteric, active site. As one of the most direct and efficient ways to regulate protein function, allostery has played a fundamental role in innumerable biological processes of all living organisms, including enzyme catalysis, signal transduction, cell metabolism, and gene transcription. It is thus considered as "the second secret of life". The abnormality of allosteric communication networks between allosteric and orthosteric sites is associated with the pathogenesis of human diseases. Allosteric modulators, by attaching to structurally diverse allosteric sites, offer the potential for differential selectivity and improved safety compared with orthosteric drugs that bind to conserved orthosteric sites. Harnessing allostery has thus been regarded as a novel strategy for drug discovery. Despite much progress having been made in the repertoire of allostery since the turn of the millennium, the identification of allosteric drugs for therapeutic targets and the elucidation of allosteric mechanisms still present substantial challenges. These challenges are derived from the difficulties in the identification of allosteric sites and mutations, the assessment of allosteric protein-modulator interactions, the screening of allosteric modulators, and the elucidation of allosteric mechanisms in biological systems. To address these issues, we have developed a panel of allosteric services for specific allosteric applications over the past decade, including (i) the creation of the Allosteric Database, with the aim of providing comprehensive allosteric information such as allosteric proteins, modulators, sites, pathways, etc., (ii) the construction of the ASBench benchmark of high-quality allosteric sites for the development of computational methods for predicting allosteric sites, (iii) the development of Allosite and AllositePro for the prediction of the location of allosteric sites in proteins, (iv) the development of the Alloscore scoring function for the evaluation of allosteric protein-modulator interactions, (v) the development of Allosterome for evolutionary analysis of query allosteric sites/modulators within the human proteome, (vi) the development of AlloDriver for the prediction of allosteric mutagenesis, and (vii) the development of AlloFinder for the virtual screening of allosteric modulators and the investigation of allosteric mechanisms. Importantly, we have validated computationally predicted allosteric sites, mutations, and modulators in the real cases of sirtuin 6, casein kinase 2α, phosphodiesterase 10A, and signal transduction and activation of transcription 3. Furthermore, our developed allosteric methods have been widely exploited by other users around the world for allosteric research. Therefore, these allosteric services are expected to expedite the discovery of allosteric drugs and the investigation of allosteric mechanisms.

Are the Apo Proteins Suitable for the Rational Discovery of Allosteric Drugs?

Allosteric modulators, by targeting the less-conserved allosteric sites, represent an innovative strategy in drug discovery. These modulators have a distinctive advantage over orthosteric ligands that attach to the conserved, functional orthosteric sites. However, in structure-based drug design, it remains unclear whether allosteric protein structures determined without orthosteric ligand binding are suitable for allosteric drug screening. In this study, we performed large-scale conformational samplings of six representative allosteric proteins uncomplexed ( apo) and complexed ( holo) with orthosteric ligands to explore the effect of orthosteric site binding on the conformational dynamics of allosteric sites. The results, coupled with the redocking evaluation of allosteric modulators to their apo and holo proteins using their MD trajectories, indicated that orthosteric site binding had an effect on the dynamics of the allosteric sites and allosteric modulators preferentially bound to their holo proteins. According to the analysis data, we constructed a new correlation model for quantifying the allosteric site change driven by substrate binding to the orthosteric site. These results highlight the strong demand to select holo allosteric proteins as initial inputs in structure-based allosteric drug screening when the distance between orthosteric and allosteric sites in the protein is below 5 Å, which is expected to contribute to allosteric drug discovery.

Molecular dynamics simulations provide insights into the origin of gleevec's selectivity toward human tyrosine kinases

Protein kinases are critical drug targets against cancer. Since the discovery of Gleevec, a specific inhibitor of Abl kinase, the capability of this drug to distinguish between Abl and other tyrosine kinases, such as Src, has been intensely investigated but the origin of Gleevec's selectivity to Abl against Src is less studied. Here, we performed molecular dynamics (MD) simulations, dynamical cross-correlation matrices (DCCM), dynamical network analysis, and binding free energy calculations to explore Gleevec's selectivity based on the crystal structures of Abl, Src, and their common ancestors (ANC-AS) and the two constructed mutation systems (AS→Abl and AS→Src). MD simulations revealed that the conformation of the phosphate-binding loop (P-loop) was altered significantly in the AS→Abl system. DCCM results unraveled that mutations increased anticorrelated motions in the AS→Abl system. Community network analysis suggested that the P-loop established special contacts in the AS→Abl system that are devoid in the AS→Src system. The binding free energy calculations unveiled that the affinity of Gleevec to AS→Abl increased to near the Abl level, whereas its affinity to AS→Src decreased to near the Src level. Analysis of individual residue contributions showed that the differences were located mainly at the P-loop. This study is valuable for understanding the sensitivity of Gleevec to human tyrosine kinases. Communicated by Ramaswamy H. Sarma.

The Taxonomy of Covalent Inhibitors

Covalent enzyme inhibitors are widely applied as biochemical tools and therapeutic agents. As a complement to categorization of these inhibitors by reactive group or modification site, we present a categorization by mechanism, which highlights common advantages and disadvantages inherent to each approach. Established categories for reversible and irreversible covalent inhibition are reviewed with representative examples given for each class, including covalent reversible inhibitors, slow substrates, residue-specific reagents, affinity labels (classical, quiescent, and photoaffinity), and mechanism-based inactivators. The relationships of these categories to proteomic profiling probes (activity-based and reactivity-based) as well as complementary approaches such as prodrug and soft drug design are also discussed. A wide variety of strategies are used to balance reactivity and selectivity in the design of covalent enzyme inhibitors. Use of a shared terminology is encouraged to clearly convey these mechanisms, to relate them to prior use of covalent inhibitors in enzymology, and to facilitate the development of more effective covalent inhibitors.

Covalent Allosteric Probe for the Metabotropic Glutamate Receptor 2: Design, Synthesis, and Pharmacological Characterization

Covalent labeling of G protein-coupled receptors (GPCRs) by small molecules is a powerful approach to understand binding modes, mechanism of action, pharmacology, and even facilitate structure elucidation. We report the first covalent positive allosteric modulator (PAM) for a class C GPCR, the mGlu₂ receptor. Three putatively covalent mGlu₂ PAMs were designed and synthesized. Pharmacological characterization identified 2 to bind the receptor covalently. Computational modeling combined with receptor mutagenesis revealed T791^7.29×30 as the likely position of covalent interaction. We show how this covalent ligand can be used to characterize the PAM binding mode and that it is a valuable tool compound in studying receptor function and binding kinetics. Our findings advance the understanding of the mGlu₂ PAM interaction and suggest that 2 is a valuable probe for further structural and chemical biology approaches.

Engineered control of enzyme structural dynamics and function

Enzymes undergo a range of internal motions from local, active site fluctuations to large-scale, global conformational changes. These motions are often important for enzyme function, including in ligand binding and dissociation and even preparing the active site for chemical catalysis. Protein engineering efforts have been directed towards manipulating enzyme structural dynamics and conformational changes, including targeting specific amino acid interactions and creation of chimeric enzymes with new regulatory functions. Post-translational covalent modification can provide an additional level of enzyme control. These studies have not only provided insights into the functional role of protein motions, but they offer opportunities to create stimulus-responsive enzymes. These enzymes can be engineered to respond to a number of external stimuli, including light, pH, and the presence of novel allosteric modulators. Altogether, the ability to engineer and control enzyme structural dynamics can provide new tools for biotechnology and medicine.

Computational Insights into the Interactions between Calmodulin and the c/nSH2 Domains of p85α Regulatory Subunit of PI3Kα: Implication for PI3Kα Activation by Calmodulin

Calmodulin (CaM) and phosphatidylinositide-3 kinase (PI3Kα) are well known for their multiple roles in a series of intracellular signaling pathways and in the progression of several human cancers. Crosstalk between CaM and PI3Kα has been an area of intensive research. Recent experiments have shown that in adenocarcinoma, K-Ras4B is involved in the CaM-PI3Kα crosstalk. Based on experimental results, we have recently put forward a hypothesis that the coordination of CaM and PI3Kα with K-Ras4B forms a CaM-PI3Kα-K-Ras4B ternary complex, which leads to the formation of pancreatic ductal adenocarcinoma. However, the mechanism for the CaM-PI3Kα crosstalk is unresolved. Based on molecular modeling and molecular dynamics simulations, here we explored the potential interactions between CaM and the c/nSH2 domains of p85α subunit of PI3Kα. We demonstrated that CaM can interact with the c/nSH2 domains and the interaction details were unraveled. Moreover, the possible modes for the CaM-cSH2 and CaM-nSH2 interactions were uncovered and we used them to construct a complete CaM-PI3Kα complex model. The structural model of CaM-PI3Kα interaction not only offers a support for our previous ternary complex hypothesis, but also is useful for drug design targeted at CaM-PI3Kα protein-protein interactions.

Molecular Dynamics Simulations and Dynamic Network Analysis Reveal the Allosteric Unbinding of Monobody to H-Ras Triggered by R135K Mutation

Ras proteins, as small GTPases, mediate cell proliferation, survival and differentiation. Ras mutations have been associated with a broad spectrum of human cancers and thus targeting Ras represents a potential way forward for cancer therapy. A recently reported monobody NS1 allosterically disrupts the Ras-mediated signaling pathway, but its efficacy is reduced by R135K mutation in H-Ras. However, the detailed mechanism is unresolved. Here, using molecular dynamics (MD) simulations and dynamic network analysis, we explored the molecular mechanism for the unbinding of NS1 to H-Ras and shed light on the underlying allosteric network in H-Ras. MD simulations revealed that the overall structures of the two complexes did not change significantly, but the H-Ras–NS1 interface underwent significant conformational alteration in the mutant Binding free energy analysis showed that NS1 binding was unfavored after R135K mutation, which resulted in the unfavorable binding of NS1. Furthermore, the critical residues on H-Ras responsible for the loss of binding of NS1 were identified. Importantly, the allosteric networks for these important residues were revealed, which yielded a novel insight into the allosteric regulatory mechanism of H-Ras.