Rational design of allosteric modulators: Challenges and successes

Advances of Predicting Allosteric Mechanisms Through Protein Contact in New Technologies and Their Application

Allostery is an intriguing phenomenon wherein the binding activity of a biological macromolecule is modulated via non-canonical binding site, resulting in synchronized functional changes. The mechanics underlying allostery are relatively complex and this review is focused on common methodologies used to study allostery, such as X-ray crystallography, NMR spectroscopy, and HDXMS. Different methodological approaches are used to generate data in different scenarios. For example, X-ray crystallography provides high-resolution structural information, NMR spectroscopy offers dynamic insights into allosteric interactions in solution, and HDXMS provides information on protein dynamics. The residue transition state (RTS) approach has emerged as a critical tool in understanding the energetics and conformational changes associated with allosteric regulation. Allostery has significant implications in drug discovery, gene transcription, disease diagnosis, and enzyme catalysis. Enzymes' catalytic activity can be modulated by allosteric regulation, offering opportunities to develop novel therapeutic alternatives. Understanding allosteric mechanisms associated with infectious organisms like SARS-CoV and bacterial pathogens can aid in the development of new antiviral drugs and antibiotics. Allosteric mechanisms are crucial in the regulation of a variety of signal transduction and cell metabolism pathways, which in turn govern various cellular processes. Despite progress, challenges remain in identifying allosteric sites and characterizing their contribution to a variety of biological processes. Increased understanding of these mechanisms can help develop allosteric systems specifically designed to modulate key biological mechanisms, providing novel opportunities for the development of targeted therapeutics. Therefore, the current review aims to summarize common methodologies that are used to further our understanding of allosteric mechanisms. In conclusion, this review provides insights into the methodologies used for the study of allostery, its applications in in silico modeling, the mechanisms underlying antibody allostery, and the ongoing challenges and prospects in advancing our comprehension of this intriguing phenomenon.

A millisecond coarse-grained simulation approach to decipher allosteric cannabinoid binding at the glycine receptor α1

Glycine receptors (GlyR) are regulated by small-molecule binding at several allosteric sites. Cannabinoids like tetrahydrocannabinol (THC) and N-arachidonyl-ethanol-amide (AEA) potentiate the GlyR response but their mechanism of action is not fully established. By combining millisecond coarse-grained (CG) MD simulations powered by Martini 3 with backmapping to all-atom representations, we have characterized the cannabinoid-binding site(s) at the zebrafish GlyR-α1 active state with atomic resolution. Based on hundreds of thousand ligand-binding events, we find that cannabinoids bind to the transmembrane domain of the receptor at both intrasubunit and intersubunit sites. For THC, the intrasubunit binding mode predicted in simulation is in excellent agreement with recent cryo-EM structures, while intersubunit binding recapitulates in full previous mutagenesis experiments. Intriguingly, AEA is predicted to bind at the same intersubunit site despite the strikingly different chemistry. Statistical analyses of the ligand-receptor interactions highlight potentially relevant residues for GlyR potentiation, offering experimentally testable predictions. The predictions for AEA have been validated by electrophysiology recordings of rationally designed mutants. The results highlight the existence of multiple cannabinoid-binding sites for the allosteric regulation of GlyR and put forward an effective strategy for the identification and structural characterization of allosteric binding sites.

The value of protein allostery in rational anticancer drug design: an update

INTRODUCTION

Allosteric drugs are advantageous. However, they still face hurdles, including identification of allosteric sites that will effectively alter the active site. Current strategies largely focus on identifying pockets away from the active sites into which the allosteric ligand will dock and do not account for exactly how the active site is altered. Favorable allosteric inhibitors dock into sites that are nearby the active sites and follow nature, mimicking diverse allosteric regulation strategies.

AREAS COVERED

The following article underscores the immense significance of allostery in drug design, describes current allosteric strategies, and especially offers a direction going forward. The article concludes with the authors' expert perspectives on the subject.

EXPERT OPINION

To select a productive venue in allosteric inhibitor development, we should learn from nature. Currently, useful strategies follow this route. Consider, for example, the mechanisms exploited in relieving autoinhibition and in harnessing allosteric degraders. Mimicking compensatory, or rescue mutations may also fall into such a thesis, as can molecular glues that capture features of scaffolding proteins. Capturing nature and creatively tailoring its mimicry can continue to innovate allosteric drug discovery.

Prediction of Protein Allosteric Sites with Transfer Entropy and Spatial Neighbor-Based Evolutionary Information Learned by an Ensemble Model

Allostery is one of the most direct and efficient ways to regulate protein functions. The diverse allosteric sites make it possible to design allosteric modulators of differential selectivity and improved safety compared with those of orthosteric drugs targeting conserved orthosteric sites. Here, we develop an ensemble machine learning method AllosES to predict protein allosteric sites in which the new and effective features are utilized, including the entropy transfer-based dynamic property, secondary structure features, and our previously proposed spatial neighbor-based evolutionary information besides the traditional physicochemical properties. To overcome the class imbalance problem, the multiple grouping strategy is proposed, which is applied to feature selection and model construction. The ensemble model is constructed where multiple submodels are trained on multiple training subsets, respectively, and their results are then integrated to be the final output. AllosES achieves a prediction performance of 0.556 MCC on the independent test set D24, and additionally, AllosES can rank the real allosteric sites in the top three for 83.3/89.3% of allosteric proteins from the test set D24/D28, outperforming the state-of-the-art peer methods. The comprehensive results demonstrate that AllosES is a promising method for protein allosteric site prediction. The source code is available at https://github.com/ChunhuaLab/AllosES.

Navigating the complexity of p53-DNA binding: implications for cancer therapy

Abstract

The tumor suppressor protein p53, a transcription factor playing a key role in cancer prevention, interacts with DNA as its primary means of determining cell fate in the event of DNA damage. When it becomes mutated, it opens damaged cells to the possibility of reproducing unchecked, which can lead to formation of cancerous tumors. Despite its critical role, therapies at the molecular level to restore p53 native function remain elusive, due to its complex nature. Nevertheless, considerable information has been amassed, and new means of investigating the problem have become available.

Objectives

We consider structural, biophysical, and bioinformatic insights and their implications for the role of direct and indirect readout and how they contribute to binding site recognition, particularly those of low consensus. We then pivot to consider advances in computational approaches to drug discovery.

Materials and methods

We have conducted a review of recent literature pertinent to the p53 protein.

Results

Considerable literature corroborates the idea that p53 is a complex allosteric protein that discriminates its binding sites not only via consensus sequence through direct H-bond contacts, but also a complex combination of factors involving the flexibility of the binding site. New computational methods have emerged capable of capturing such information, which can then be utilized as input to machine learning algorithms towards the goal of more intelligent and efficient de novo allosteric drug design.

Conclusions

Recent improvements in machine learning coupled with graph theory and sector analysis hold promise for advances to more intelligently design allosteric effectors that may be able to restore native p53-DNA binding activity to mutant proteins.

Clinical relevance

The ideas brought to light by this review constitute a significant advance that can be applied to ongoing biophysical studies of drugs for p53, paving the way for the continued development of new methodologies for allosteric drugs. Our discoveries hold promise to provide molecular therapeutics which restore p53 native activity, thereby offering new insights for cancer therapies.

Graphical Abstract

Structural representation of the p53 DBD (PDBID 1TUP). DNA consensus sequence is shown in gray, and the protein is shown in blue. Red beads indicate hotspot residue mutations, green beads represent DNA interacting residues, and yellow beads represent both.

STORMM: Structure and topology replica molecular mechanics for chemical simulations

The Structure and TOpology Replica Molecular Mechanics (STORMM) code is a next-generation molecular simulation engine and associated libraries optimized for performance on fast, vectorized central processor units and graphics processing units (GPUs) with independent memory and tens of thousands of threads. STORMM is built to run thousands of independent molecular mechanical calculations on a single GPU with novel implementations that tune numerical precision, mathematical operations, and scarce on-chip memory resources to optimize throughput. The libraries are built around accessible classes with detailed documentation, supporting fine-grained parallelism and algorithm development as well as copying or swapping groups of systems on and off of the GPU. A primary intention of the STORMM libraries is to provide developers of atomic simulation methods with access to a high-performance molecular mechanics engine with extensive facilities to prototype and develop bespoke tools aimed toward drug discovery applications. In its present state, STORMM delivers molecular dynamics simulations of small molecules and small proteins in implicit solvent with tens to hundreds of times the throughput of conventional codes. The engineering paradigm transforms two of the most memory bandwidth-intensive aspects of condensed-phase dynamics, particle-mesh mapping, and valence interactions, into compute-bound problems for several times the scalability of existing programs. Numerical methods for compressing and streamlining the information present in stored coordinates and lookup tables are also presented, delivering improved accuracy over methods implemented in other molecular dynamics engines. The open-source code is released under the MIT license.

Leveraging Bidirectional Nature of Allostery To Inhibit Protein-Protein Interactions (PPIs): A Case Study of PCSK9-LDLR Interaction

The protein PCSK9 (proprotein convertase subtilisin/Kexin type 9) negatively regulates the recycling of LDLR (low-density lipoprotein receptor), leading to an elevated plasma level of LDL. Inhibition of PCSK9-LDLR interaction has emerged as a promising therapeutic strategy to manage hypercholesterolemia. However, the large interaction surface area between PCSK9 and LDLR makes it challenging to identify a small molecule competitive inhibitor. An alternative strategy would be to identify distal cryptic sites as targets for allosteric inhibitors that can remotely modulate PCSK9-LDLR interaction. Using several microseconds long molecular dynamics (MD) simulations, we demonstrate that on binding with LDLR, there is a significant conformational change (population shift) in a distal loop (residues 211-222) region of PCSK9. Consistent with the bidirectional nature of allostery, we establish a clear correlation between the loop conformation and the binding affinity with LDLR. Using a thermodynamic argument, we establish that the loop conformations predominantly present in the apo state of PCSK9 would have lower LDLR binding affinity, and they would be potential targets for designing allosteric inhibitors. We elucidate the molecular origin of the allosteric coupling between this loop and the LDLR binding interface in terms of the population shift in a set of salt bridges and hydrogen bonds. Overall, our work provides a general strategy toward identifying allosteric hotspots: compare the conformational ensemble of the receptor between the apo and bound states of the protein and identify distal conformational changes, if any. The inhibitors should be designed to bind and stabilize the apo-specific conformations.

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.

Machine learning approaches in predicting allosteric sites

Allosteric regulation is a fundamental biological mechanism that can control critical cellular processes via allosteric modulator binding to protein distal functional sites. The advantages of allosteric modulators over orthosteric ones have sparked the development of numerous computational approaches, such as the identification of allosteric binding sites, to facilitate allosteric drug discovery. Building on the success of machine learning (ML) models for solving complex problems in biology and chemistry, several ML models for predicting allosteric sites have been developed. In this review, we provide an overview of these models and discuss future perspectives powered by the field of artificial intelligence such as protein language models.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

Modulation of Allostery with Multiple Mechanisms by Hotspot Mutations in TetR

Modulating allosteric coupling offers unique opportunities for biomedical applications. Such efforts can benefit from efficient prediction and evaluation of allostery hotspot residues that dictate the degree of cooperativity between distant sites. We demonstrate that effects of allostery hotspot mutations can be evaluated qualitatively and semiquantitatively by molecular dynamics simulations in a bacterial tetracycline repressor (TetR). The simulations recapitulate the effects of these mutations on abolishing the induction function of TetR and provide a rationale for the different rescuabilities observed to restore allosteric coupling of the hotspot mutations. We demonstrate that the same noninducible phenotype could be the result of perturbations in distinct structural and energetic properties of TetR. Our work underscores the value of explicitly computing the functional free energy landscapes to effectively evaluate and rank hotspot mutations despite the prevalence of compensatory interactions and therefore provides quantitative guidance to allostery modulation for therapeutic and engineering applications.

A Structure-Based Allosteric Modulator Design Paradigm

Importance: Allosteric drugs bound to topologically distal allosteric sites hold a substantial promise in modulating therapeutic targets deemed undruggable at their orthosteric sites. Traditionally, allosteric modulator discovery has predominantly relied on serendipitous high-throughput screening. Nevertheless, the landscape has undergone a transformative shift due to recent advancements in our understanding of allosteric modulation mechanisms, coupled with a significant increase in the accessibility of allosteric structural data. These factors have extensively promoted the development of various computational methodologies, especially for machine-learning approaches, to guide the rational design of structure-based allosteric modulators. Highlights: We here presented a comprehensive structure-based allosteric modulator design paradigm encompassing 3 critical stages: drug target acquisition, allosteric binding site, and modulator discovery. The recent advances in computational methods in each stage are encapsulated. Furthermore, we delve into analyzing the successes and obstacles encountered in the rational design of allosteric modulators. Conclusion: The structure-based allosteric modulator design paradigm holds immense potential for the rational design of allosteric modulators. We hope that this review would heighten awareness of the use of structure-based computational methodologies in advancing the field of allosteric drug discovery.

Modulation of Allostery with Multiple Mechanisms by Hotspot Mutations in TetR

Modulating allosteric coupling offers unique opportunities for biomedical applications. Such efforts can benefit from efficient prediction and evaluation of allostery hotspot residues that dictate the degree of co-operativity between distant sites. We demonstrate that effects of allostery hotspot mutations can be evaluated qualitatively and semi-quantitatively by molecular dynamics simulations in a bacterial tetracycline repressor (TetR). The simulations recapitulate the effects of these mutations on abolishing the induction function of TetR and provide a rationale for the different degrees of rescuability observed to restore allosteric coupling of the hotspot mutations. We demonstrate that the same non-inducible phenotype could be the result of perturbations in distinct structural and energetic properties of TetR. Our work underscore the value of explicitly computing the functional free energy landscapes to effectively evaluate and rank hotspot mutations despite the prevalence of compensatory interactions, and therefore provide quantitative guidance to allostery modulation for therapeutic and engineering applications.

Unmasking an Allosteric Binding Site of the Papain-like Protease in SARS-CoV-2: Molecular Dynamics Simulations of Corticosteroids

To date, mechanistic insights into many clinical drugs against COVID-19 remain unexplored. Dexamethasone, a corticosteroid, is one of them. While treating the entire corticosteroid database, including vitamins D2 and D3, with cutting-edge computational techniques, several intriguing results are unfolded. From the top-notch candidates, dexamethasone is likely to inhibit the viral main protease (Mpro), with vitamin D3 exhibiting multitarget [Mpro, papain-like protease (PLpro), and nucleocapsid protein (N-pro)] roles and ciclesonide's dynamic flipping disinterring a cryptic allosteric site in the PLpro enzyme. The results rationalize why these drugs improve the health of COVID-19 patients. Understanding an enzyme's secret binding site is essential to understanding how the enzyme works and how to inhibit its function. Ciclesonide's allosteric inhibition could not only jeopardize PLpro's catalytic role in polyprotein processing but also make it less vulnerable to the host body's defense machinery. Hotspot residues in the identified allosteric site could be considered for effective therapeutic designs against PLpro.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Discovery of Hit Compounds Targeting the P4 Allosteric Site of K-RAS, Identified through Ensemble-Based Virtual Screening

Mutants of Ras are oncogenic drivers of a large number of human tumors. Despite being recognized as an attractive target for the treatment of cancer, the high affinity for its substrate tagged the protein as undruggable for a few years. The identification of cryptic pockets on the protein surface gave the opportunity to identify molecules capable of acting as allosteric modulators. Several molecules were disclosed in recent years, with sotorasib and adagrasib already approved for clinical use. The present study makes use of computational methods to characterize eight prospective allosteric pockets (P1-P8) in K-Ras, four of which (P1-P4) were previously characterized in the literature. The present study also describes the results of a virtual screening study focused on the discovery of hit compounds, binders of the P4 site that can be considered as peptidomimetics of a fragment of the SOS αI helix, a guanine exchange factor of Ras. After a detailed description of the computational procedure followed, we disclose five hit compounds, prospective binders of the P4 allosteric site that exhibit an inhibitory capability higher than 30% in a cell proliferation assay at 50 μM.

CryoEM analysis of the essential native UDP-glucose pyrophosphorylase from Aspergillus nidulans reveals key conformations for activity regulation and function

Invasive aspergillosis is one of the most serious clinical invasive fungal infections, resulting in a high case fatality rate among immunocompromised patients. The disease is caused by saprophytic molds in the genus Aspergillus, including Aspergillus fumigatus, the most significant pathogenic species. The fungal cell wall, an essential structure mainly composed of glucan, chitin, galactomannan, and galactosaminogalactan, represents an important target for the development of antifungal drugs. UDP (uridine diphosphate)-glucose pyrophosphorylase (UGP) is a central enzyme in the metabolism of carbohydrates that catalyzes the biosynthesis of UDP-glucose, a key precursor of fungal cell wall polysaccharides. Here, we demonstrate that the function of UGP is vital for Aspergillus nidulans (AnUGP). To understand the molecular basis of AnUGP function, we describe a cryoEM structure (global resolution of 3.5 Å for the locally refined subunit and 4 Å for the octameric complex) of a native AnUGP. The structure reveals an octameric architecture with each subunit comprising an N-terminal α-helical domain, a central catalytic glycosyltransferase A-like (GT-A-like) domain, and a C-terminal (CT) left-handed β-helix oligomerization domain. AnUGP displays unprecedented conformational variability between the CT oligomerization domain and the central GT-A-like catalytic domain. In combination with activity measurements and bioinformatics analysis, we unveil the molecular mechanism of substrate recognition and specificity for AnUGP. Altogether, our study not only contributes to understanding the molecular mechanism of catalysis/regulation of an important class of enzymes but also provides the genetic, biochemical, and structural groundwork for the future exploitation of UGP as a potential antifungal target. IMPORTANCE Fungi cause diverse diseases in humans, ranging from allergic syndromes to life-threatening invasive diseases, together affecting more than a billion people worldwide. Increasing drug resistance in Aspergillus species represents an emerging global health threat, making the design of antifungals with novel mechanisms of action a worldwide priority. The cryoEM structure of UDP (uridine diphosphate)-glucose pyrophosphorylase (UGP) from the filamentous fungus Aspergillus nidulans reveals an octameric architecture displaying unprecedented conformational variability between the C-terminal oligomerization domain and the central glycosyltransferase A-like catalytic domain in the individual protomers. While the active site and oligomerization interfaces are more highly conserved, these dynamic interfaces include motifs restricted to specific clades of filamentous fungi. Functional study of these motifs could lead to the definition of new targets for antifungals inhibiting UGP activity and, thus, the architecture of the cell wall of filamentous fungal pathogens.

Structure-Activity Relationship Study of the High-Affinity Neuropeptide Y4 Receptor Positive Allosteric Modulator VU0506013

Positive allosteric modulators targeting the Y₄ receptor (Y₄R), a G protein-coupled receptor (GPCR) involved in the regulation of satiety, offer great potential in anti-obesity research. In this study, we selected 603 compounds by using quantitative structure-activity relationship (QSAR) models and tested them in high-throughput screening (HTS). Here, the novel positive allosteric modulator (PAM) VU0506013 was identified, which exhibits nanomolar affinity and pronounced selectivity toward the Y₄R in engineered cell lines and mouse descending colon mucosa natively expressing the Y₄R. Based on this lead structure, we conducted a systematic SAR study in two regions of the scaffold and presented a series of 27 analogues with modifications in the N- and C-terminal heterocycles of the molecule to obtain insight into functionally relevant positions. By mutagenesis and computational docking, we present a potential binding mode of VU0506013 in the transmembrane core of the Y₄R. VU0506013 presents a promising scaffold for developing in vivo tools to move toward anti-obesity drug research focused on the Y₄R.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K⁺, providing evidence for synergism between lysine/K⁺ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K⁺. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K⁺ binding, revealing that riboswitch folding (k_fold) and unfolding (k_unfold) rate constants increase and decrease dramatically with K⁺, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K⁺ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K⁺ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na⁺ and Cs⁺ studies to demonstrate that K⁺ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K⁺ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

How protein topology controls allosteric regulations

Allostery is an important regulatory mechanism of protein functions. Among allosteric proteins, certain protein structure types are more observed. However, how allosteric regulation depends on protein topology remains elusive. In this study, we extracted protein topology graphs at the fold level and found that known allosteric proteins mainly contain multiple domains or subunits and allosteric sites reside more often between two or more domains of the same fold type. Only a small fraction of fold-fold combinations are observed in allosteric proteins, and homo-fold-fold combinations dominate. These analyses imply that the locations of allosteric sites including cryptic ones depend on protein topology. We further developed TopoAlloSite, a novel method that uses the kernel support vector machine to predict the location of allosteric sites on the overall protein topology based on the subgraph-matching kernel. TopoAlloSite successfully predicted known cryptic allosteric sites in several allosteric proteins like phosphopantothenoylcysteine synthetase, spermidine synthase, and sirtuin 6, demonstrating its power in identifying cryptic allosteric sites without performing long molecular dynamics simulations or large-scale experimental screening. Our study demonstrates that protein topology largely determines how its function can be allosterically regulated, which can be used to find new druggable targets and locate potential binding sites for rational allosteric drug design.

Potential Self-Peptide Inhibitors of the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in viral replication, cleaving viral polyproteins into functional proteins. This makes Mpro an important drug target. Mpro consists of an N-terminal catalytic domain and a C-terminal α-helical domain (MproC). Previous studies have shown that peptides derived from a given protein sequence (self-peptides) can affect the folding and, in turn, the function of that protein. Since the SARS-CoV-1 MproC is known to stabilize its Mpro and regulate its function, we hypothesized that SARS-CoV-2 MproC-derived self-peptides may modulate the folding and the function of SARS-CoV-2 Mpro. To test this, we studied the folding of MproC in the presence of various self-peptides using coarse-grained structure-based models and molecular dynamics simulations. In these simulations of MproC and one self-peptide, we found that two self-peptides, the α1-helix and the loop between α4 and α5 (loop4), could replace the equivalent native sequences in the MproC structure. Replacement of either sequence in full-length Mpro should, in principle, be able to perturb Mpro function albeit through different mechanisms. Some general principles for the rational design of self-peptide inhibitors emerge: The simulations show that prefolded self-peptides are more likely to replace native sequences than those which do not possess structure. Additionally, the α1-helix self-peptide is kinetically stable and once inserted rarely exchanges with the native α1-helix, while the loop4 self-peptide is easily replaced by the native loop4, making it less useful for modulating function. In summary, a prefolded α1-derived peptide should be able to inhibit SARS-CoV-2 Mpro function.

Development of a cell-free screening assay for the identification of direct PERK activators

The activation of the unfolded protein response, particularly via the PERK pathway, has been suggested as a promising therapeutic approach in tauopathies, a group of neurodegenerative disorders characterized by the abnormal phosphorylation and aggregation of tau protein. So far, a shortage of available direct PERK activators has been limiting the progresses in this field. Our study aimed at the development of a cell-free screening assay enabling the detection of novel direct PERK activators. By applying the catalytic domain of recombinant human PERK, we initially determined ideal conditions of the kinase assay reaction, including parameters such as optimal kinase concentration, temperature, and reaction time. Instead of using PERK's natural substrate proteins, eIF2α and NRF2, we applied SMAD3 as phosphorylation-accepting protein and successfully detected cell-free PERK activation and inhibition by selected modulators (e.g., calcineurin-B, GSK2606414). The developed assay revealed to be sufficiently stable and robust to assess an activating EC50-value. Additionally, our results suggested that PERK activation may take place independent of the active site which can be blocked by a kinase inhibitor. Finally, we confirmed the applicability of the assay by measuring PERK activation by MK-28, a recently described PERK activator. Overall, our data show that a cell-free luciferase-based assay with the recombinant human PERK kinase domain and SMAD3 as substrate protein is capable of detecting PERK activation, which enables to screen large compound libraries for direct PERK activators, in a high-throughput-based approach. These activators will be useful for deepening our understanding of the PERK signaling pathway, and may also lead to the identification of new therapeutic drug candidates for neurodegenerative tauopathies.

DREAMM: a web-based server for drugging protein-membrane interfaces as a novel workflow for targeted drug design

SUMMARY

The allosteric modulation of peripheral membrane proteins (PMPs) by targeting protein-membrane interactions with drug-like molecules represents a new promising therapeutic strategy for proteins currently considered undruggable. However, the accessibility of protein-membrane interfaces by small molecules has been so far unexplored, possibly due to the complexity of the interface, the limited protein-membrane structural information and the lack of computational workflows to study it. Herein, we present a pipeline for drugging protein-membrane interfaces using the DREAMM (Drugging pRotein mEmbrAne Machine learning Method) web server. DREAMM works in the back end with a fast and robust ensemble machine learning algorithm for identifying protein-membrane interfaces of PMPs. Additionally, DREAMM also identifies binding pockets in the vicinity of the predicted membrane-penetrating amino acids in protein conformational ensembles provided by the user or generated within DREAMM.

AVAILABILITY AND IMPLEMENTATION

DREAMM web server is accessible via https://dreamm.ni4os.eu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploring the druggability of the binding site of aurovertin, an exogenous allosteric inhibitor of FOF1-ATP synthase

In addition to playing a central role in the mitochondria as the main producer of ATP, F_OF₁-ATP synthase performs diverse key regulatory functions in the cell membrane. Its malfunction has been linked to a growing number of human diseases, including hypertension, atherosclerosis, cancer, and some neurodegenerative, autoimmune, and aging diseases. Furthermore, inhibition of this enzyme jeopardizes the survival of several bacterial pathogens of public health concern. Therefore, F_OF₁-ATP synthase has emerged as a novel drug target both to treat human diseases and to combat antibiotic resistance. In this work, we carried out a computational characterization of the binding sites of the fungal antibiotic aurovertin in the bovine F₁ subcomplex, which shares a large identity with the human enzyme. Molecular dynamics simulations showed that although the binding sites can be described as preformed, the inhibitor hinders inter-subunit communications and exerts long-range effects on the dynamics of the catalytic site residues. End-point binding free energy calculations revealed hot spot residues for aurovertin recognition. These residues were also relevant to stabilize solvent sites determined from mixed-solvent molecular dynamics, which mimic the interaction between aurovertin and the enzyme, and could be used as pharmacophore constraints in virtual screening campaigns. To explore the possibility of finding species-specific inhibitors targeting the aurovertin binding site, we performed free energy calculations for two bacterial enzymes with experimentally solved 3D structures. Finally, an analysis of bacterial sequences was carried out to determine conservation of the aurovertin binding site. Taken together, our results constitute a first step in paving the way for structure-based development of new allosteric drugs targeting F_OF₁-ATP synthase sites of exogenous inhibitors.

The orchestrated signaling by PI3Kα and PTEN at the membrane interface

The oncogene PI3Kα and the tumor suppressor PTEN represent two antagonistic enzymatic activities that regulate the interconversion of the phosphoinositide lipids PI(4,5)P2 and PI(3,4,5)P3 in membranes. As such, they are defining components of phosphoinositide-based cellular signaling and membrane trafficking pathways that regulate cell survival, growth, and proliferation, and are often deregulated in cancer. In this review, we highlight aspects of PI3Kα and PTEN interplay at the intersection of signaling and membrane trafficking. We also discuss the mechanisms of PI3Kα- and PTEN- membrane interaction and catalytic activation, which are fundamental for our understanding of the structural and allosteric implications on signaling at the membrane interface and may aid current efforts in pharmacological targeting of these proteins.

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.

Identification and Inhibition of the Druggable Allosteric Site of SARS-CoV-2 NSP10/NSP16 Methyltransferase through Computational Approaches

Since its emergence in early 2019, the respiratory infectious virus, SARS-CoV-2, has ravaged the health of millions of people globally and has affected almost every sphere of life. Many efforts are being made to combat the COVID-19 pandemic’s emerging and recurrent waves caused by its evolving and more infectious variants. As a result, novel and unexpected targets for SARS-CoV-2 have been considered for drug discovery. 2′-O-Methyltransferase (nsp10/nsp16) is a significant and appealing target in the SARS-CoV-2 life cycle because it protects viral RNA from the host degradative enzymes via a cap formation process. In this work, we propose prospective allosteric inhibitors that target the allosteric site, SARS-CoV-2 MTase. Four drug libraries containing ~119,483 compounds were screened against the allosteric site of SARS-CoV-2 MTase identified in our research. The identified best compounds exhibited robust molecular interactions and alloscore-score rankings with the allosteric site of SARS-CoV-2 MTase. Moreover, to further assess the dynamic stability of these compounds (CHEMBL2229121, ZINC000009464451, SPECS AK-91811684151, NCI-ID = 715319), a 100 ns molecular dynamics simulation, along with its holo-form, was performed to provide insights on the dynamic nature of these allosteric inhibitors at the allosteric site of the SARS-CoV-2 MTase. Additionally, investigations of MM-GBSA binding free energies revealed a good perspective for these allosteric inhibitor–enzyme complexes, indicating their robust antagonistic action on SARS-CoV-2 (nsp10/nsp16) methyltransferase. We conclude that these allosteric repressive agents should be further evaluated through investigational assessments in order to combat the proliferation of SARS-CoV-2.

Affinity Assays for Cannabinoids Detection: Are They Amenable to On-Site Screening?

Roadside testing of illicit drugs such as tetrahydrocannabinol (THC) requires simple, rapid, and cost-effective methods. The need for non-invasive detection tools has led to the development of selective and sensitive platforms, able to detect phyto- and synthetic cannabinoids by means of their main metabolites in breath, saliva, and urine samples. One may estimate the time passed from drug exposure and the frequency of use by corroborating the detection results with pharmacokinetic data. In this review, we report on the current detection methods of cannabinoids in biofluids. Fluorescent, electrochemical, colorimetric, and magnetoresistive biosensors will be briefly overviewed, putting emphasis on the affinity formats amenable to on-site screening, with possible applications in roadside testing and anti-doping control.

Molecular Dynamics Simulations Establish the Molecular Basis for the Broad Allostery Hotspot Distributions in the Tetracycline Repressor

It is imperative to identify the network of residues essential to the allosteric coupling for the purpose of rationally engineering allostery in proteins. Deep mutational scanning analysis has emerged as a function-centric approach for identifying such allostery hotspots in a comprehensive and unbiased fashion, leading to observations that challenge our understanding of allostery at the molecular level. Specifically, a recent deep mutational scanning study of the tetracycline repressor (TetR) revealed an unexpectedly broad distribution of allostery hotspots throughout the protein structure. Using extensive molecular dynamics simulations (up to 50 μs) and free energy computations, we establish the molecular and energetic basis for the strong anticooperativity between the ligand and DNA binding sites. The computed free energy landscapes in different ligation states illustrate that allostery in TetR is well described by a conformational selection model, in which the apo state samples a broad set of conformations, and specific ones are selectively stabilized by either ligand or DNA binding. By examining a range of structural and dynamic properties of residues at both local and global scales, we observe that various analyses capture different subsets of experimentally identified hotspots, suggesting that these residues modulate allostery in distinct ways. These results motivate the development of a thermodynamic model that qualitatively explains the broad distribution of hotspot residues and their distinct features in molecular dynamics simulations. The multifaceted strategy that we establish here for hotspot evaluations and our insights into their mechanistic contributions are useful for modulating protein allostery in mechanistic and engineering studies.

Opioid ligands addressing unconventional binding sites and more than one opioid receptor subtype

Opioid receptors (ORs) represent one of the most significant groups of G-protein coupled receptor (GPCR) drug targets and also act as prototypical models for GPCR function. In a constant effort to develop drugs with less side effects, and tools to explore the ORs nature and function, various (poly)pharmacological ligand design approaches have been performed. That is, besides classical ligands, a great number of bivalent ligands (i.e. aiming on two distinct OR subtypes), univalent heteromer-selective ligands and bitopic and allosteric ligands have been synthesized for the ORs. The scope of our review is to present the most important of the aforementioned ligands, highlight their properties and exhibit the current state-of-the-art pallet of promising drug candidates or useful molecular tools for the ORs.

Computational Design of Inhibitors Targeting the Catalytic β Subunit of Escherichia coli FOF1-ATP Synthase

With the uncontrolled growth of multidrug-resistant bacteria, there is an urgent need to search for new therapeutic targets, to develop drugs with novel modes of bactericidal action. FoF1-ATP synthase plays a crucial role in bacterial bioenergetic processes, and it has emerged as an attractive antimicrobial target, validated by the pharmaceutical approval of an inhibitor to treat multidrug-resistant tuberculosis. In this work, we aimed to design, through two types of in silico strategies, new allosteric inhibitors of the ATP synthase, by targeting the catalytic β subunit, a centerpiece in communication between rotor subunits and catalytic sites, to drive the rotary mechanism. As a model system, we used the F1 sector of Escherichia coli, a bacterium included in the priority list of multidrug-resistant pathogens. Drug-like molecules and an IF1-derived peptide, designed through molecular dynamics simulations and sequence mining approaches, respectively, exhibited in vitro micromolar inhibitor potency against F1. An analysis of bacterial and Mammalia sequences of the key structural helix-turn-turn motif of the C-terminal domain of the β subunit revealed highly and moderately conserved positions that could be exploited for the development of new species-specific allosteric inhibitors. To our knowledge, these inhibitors are the first binders computationally designed against the catalytic subunit of FOF1-ATP synthase.

Targeting key proteins involved in transcriptional regulation for cancer therapy: Current strategies and future prospective

The key proteins involved in transcriptional regulation play convergent roles in cellular homeostasis, and their dysfunction mediates aberrant gene expressions that underline the hallmarks of tumorigenesis. As tumor progression is dependent on such abnormal regulation of transcription, it is important to discover novel chemical entities as antitumor drugs that target key tumor-associated proteins involved in transcriptional regulation. Despite most key proteins (especially transcription factors) involved in transcriptional regulation are historically recognized as undruggable targets, multiple targeting approaches at diverse levels of transcriptional regulation, such as epigenetic intervention, inhibition of DNA-binding of transcriptional factors, and inhibition of the protein-protein interactions (PPIs), have been established in preclinically or clinically studies. In addition, several new approaches have recently been described, such as targeting proteasomal degradation and eliciting synthetic lethality. This review will emphasize on accentuating these developing therapeutic approaches and provide a thorough conspectus of the drug development to target key proteins involved in transcriptional regulation and their impact on future oncotherapy.

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.

Design and engineering of allosteric communications in proteins

Allostery in proteins plays an important role in regulating protein activities and influencing many biological processes such as gene expression, enzyme catalysis, and cell signaling. The process of allostery takes place when a signal detected at a site on a protein is transmitted via a mechanical pathway to a functional site and, thus, influences its activity. The pathway of allosteric communication consists of amino acids that form a network with covalent and non-covalent bonds. By mutating residues in this allosteric network, protein engineers have successfully established novel allosteric pathways to achieve desired properties in the target protein. In this review, we highlight the most recent and state-of-the-art techniques for allosteric communication engineering. We also discuss the challenges that need to be overcome and future directions for engineering protein allostery.

Predicting protein–membrane interfaces of peripheral membrane proteins using ensemble machine learning

Abnormal protein–membrane attachment is involved in deregulated cellular pathways and in disease. Therefore, the possibility to modulate protein–membrane interactions represents a new promising therapeutic strategy for peripheral membrane proteins that have been considered so far undruggable. A major obstacle in this drug design strategy is that the membrane-binding domains of peripheral membrane proteins are usually unknown. The development of fast and efficient algorithms predicting the protein–membrane interface would shed light into the accessibility of membrane–protein interfaces by drug-like molecules. Herein, we describe an ensemble machine learning methodology and algorithm for predicting membrane-penetrating amino acids. We utilize available experimental data from the literature for training 21 machine learning classifiers and meta-classifiers. Evaluation of the best ensemble classifier model accuracy yields a macro-averaged F1 score = 0.92 and a Matthews correlation coefficient = 0.84 for predicting correctly membrane-penetrating amino acids on unknown proteins of a validation set. The python code for predicting protein–membrane interfaces of peripheral membrane proteins is available at https://github.com/zoecournia/DREAMM.

Wandering beyond small molecules: peptides as allosteric protein modulators

Recent years have seen the rise of allosteric modulation as an innovative approach for drug design and discovery, efforts which culminated in the development of several clinical candidates. Allosteric modulation of many drug targets, including mainly membrane-embedded receptors, have been vastly explored through small molecule screening campaigns, but much less attention has been paid to peptide-based allosteric modulators. However, peptides have a significant impact on the pharmaceutical industry due to the typically higher potency and selectivity for their targets, as compared with small molecule therapeutics. Therefore, peptides represent one of the most promising classes of molecules that can modulate key biological pathways. Here, we report on the allosteric modulation of proteins (ranging from G protein-coupled receptors to specific protein-protein interactions) by peptides for applications in drug discovery.