Key Allosteric and Active Site Residues of SARS-CoV-2 3CLpro Are Promising Drug Targets.

The main protease of SARS-CoV-2, 3-chymotrypsin-like protease (3CLpro), is a prominent target for antiviral development due to its essential role in the viral life cycle. Research has largely focused on competitive inhibitors of 3CLpro that target the active site. However, allosteric sites distal to the peptide substrate-binding region are also potential targets for the design of reversible noncompetitive inhibitors. Computational analyses have examined the importance of key contacts at allosteric sites of 3CLpro, but these contacts have not been validated experimentally. In this work, four druggable pockets spanning the surface of SARS-CoV-2 3CLpro were predicted: pocket 1 is the active site, whereas pockets 2, 3, and 4 are located away from the active site at the interface of domains II and III. Site-directed alanine mutagenesis of selected residues with important structural interactions revealed that 7 of 13 active site residues (N28, R40, Y54, S147, Y161, D187 and Q192) and 7 of 12 allosteric site residues (T111, R131, N133, D197, N203, D289 and D295) are essential for maintaining catalytically active and thermodynamically stable 3CLpro. Alanine substitution at these key amino acid residues inactivated or reduced the activity of 3CLpro. In addition, the thermodynamic stability of 3CLpro decreased in the presence of some of these mutations. This work provides experimental validation of essential contacts in the active and allosteric sites of 3CLpro that could be targeted with competitive and noncompetitive inhibitors as new therapeutics against COVID-19.

GABAA receptor-mediated seizure liabilities: a mixed-methods screening approach.

GABAA receptors, members of the pentameric ligand-gated ion channel superfamily, are widely expressed in the central nervous system and mediate a broad range of pharmaco-toxicological effects including bidirectional changes to seizure threshold. Thus, detection of GABAA receptor-mediated seizure liabilities is a big, partly unmet need in early preclinical drug development. This is in part due to the plethora of allosteric binding sites that are present on different subtypes of GABAA receptors and the critical lack of screening methods that detect interactions with any of these sites. To improve in silico screening methods, we assembled an inventory of allosteric binding sites based on structural data. Pharmacophore models representing several of the binding sites were constructed. These models from the NeuroDeRisk IL Profiler were used for in silico screening of a compiled collection of drugs with known GABAA receptor interactions to generate testable hypotheses. Amoxapine was one of the hits identified and subjected to an array of in vitro assays to examine molecular and cellular effects on neuronal excitability and in vivo locomotor pattern changes in zebrafish larvae. An additional level of analysis for our compound collection is provided by pharmacovigilance alerts using FAERS data. Inspired by the Adverse Outcome Pathway framework, we postulate several candidate pathways leading from specific binding sites to acute seizure induction. The whole workflow can be utilized for any compound collection and should inform about GABAA receptor-mediated seizure risks more comprehensively compared to standard displacement screens, as it rests chiefly on functional data.

Docking-Based Evidence for the Potential of ImmunoDefender: A Novel Formulated Essential Oil Blend Incorporating Synergistic Antiviral Bioactive Compounds as Promising Mpro Inhibitors against SARS-CoV-2.

Essential oils (Eos) have demonstrated antiviral activity, but their toxicity can hinder their use as therapeutic agents. Recently, some essential oil components have been used within safe levels of acceptable daily intake limits without causing toxicity. The "ImmunoDefender," a novel antiviral compound made from a well-known mixture of essential oils, is considered highly effective in treating SARS-CoV-2 infections. The components and doses were chosen based on existing information about their structure and toxicity. Blocking the main protease (Mpro) of SARS-CoV-2 with high affinity and capacity is critical for inhibiting the virus's pathogenesis and transmission. In silico studies were conducted to examine the molecular interactions between the main essential oil components in "ImmunoDefender" and SARS-CoV-2 Mpro. The screening results showed that six key components of ImmunoDefender formed stable complexes with Mpro via its active catalytic site with binding energies ranging from -8.75 to -10.30 kcal/mol, respectively for Cinnamtannin B1, Cinnamtannin B2, Pavetannin C1, Syzyginin B, Procyanidin C1, and Tenuifolin. Furthermore, three essential oil bioactive inhibitors, Cinnamtannin B1, Cinnamtannin B2, and Pavetannin C, had significant ability to bind to the allosteric site of the main protease with binding energies of -11.12, -10.74, and -10.79 kcal/mol; these results suggest that these essential oil bioactive compounds may play a role in preventing the attachment of the translated polyprotein to Mpro, inhibiting the virus's pathogenesis and transmission. These components also had drug-like characteristics similar to approved and effective drugs, suggesting that further pre-clinical and clinical studies are needed to confirm the generated in silico outcomes.

Molecular basis of drug resistance in smoothened receptor: An in silico study of protein resistivity and specificity.

Smoothened (SMO) antagonist Vismodegib effectively inhibits the Hedgehog pathway in proliferating cancer cells. In early stage of treatment, Vismodegib exhibited promising outcomes to regress the tumors cells, but ultimately relapsed due to the drug resistive mutations in SMO mostly occurring before (primary mutations G497W) or after (acquired mutations D473H/Y) anti-SMO therapy. This study investigates the unprecedented insights of structural and functional mechanism hindering the binding of Vismodegib with sensitive and resistant mutant variants of SMO (SMOMut ). Along with the basic dynamic understanding of Vismodegib-SMO complexes, network propagation theory based on heat diffusion principles is first time applied here to identify the modules of residues influenced by the individual mutations. The allosteric modulation by GLY497 residue in Vismodegib bound SMO wild-type (SMOWT ) conformation depicts the interconnections of intermediate residues of SMO with the atom of Vismodegib and identify two important motifs (E-X-P-L) and (Q-A-N-V-T-I-G) mediating this allosteric regulation. In this study a novel computational framework based on the heat diffusion principle is also developed, which identify significant residues of allosteric site causing drug resistivity in SMOMut . This framework could also be useful for assessing the potential allosteric sites of different other proteins. Moreover, previously reported novel inhibitor "ZINC12368305," which is proven to make an energetically favorable complex with SMOWT is chosen as a control sample to assess the impact of receptor mutation on its binding and subsequently identify the important factors that govern binding disparity between Vismodegib and ZINC12368305 bound SMOWT/Mut conformations.

Computational insights into the molecular mechanisms of differentiated allosteric modulation at the mu opioid receptor by structurally similar bitopic modulators.

Targeting the mu opioid receptor (MOR) by applying orthosteric ligands is the most frequently employed method to treat opioid use disorder (OUD). Unfortunately, most of MOR orthosteric ligands produce severe side effects, mainly due to their low selectivity over other opioid receptors. In contrast, some G protein-coupled receptor allosteric modulators have been reported to exhibit high subtype selectivity and can effectively modulate the potency and/or efficacy of orthosteric ligands. Recently, NAQ and its analog NCQ were identified as novel MOR bitopic modulators. Interestingly, NAQ and NCQ were similar in structure but exhibited different efficacy profiles to the MOR. NAQ exhibited an antagonism activity to the MOR while NCQ showed a partial agonism activity to the MOR. In the present study, molecular modeling methods were applied to explore the putative molecular mechanisms of their different functional profiles to the MOR. When NAQ binding with the inactive MOR, the 'address' portion of NAQ interacted with the MOR allosteric site but showed no significant allosteric modulation of the efficacy of the 'message' portion of NAQ. However, when NCQ binding with the inactive and active MOR, the 'address' portion of NCQ seemed to be able to positively modulate the efficacy of the 'message' portion of NCQ at varying levels. Evidentially, the substituents at the 1'- and 4'-positions of the isoquinoline ring of NCQ seemed to play a critical role in the modulatory function of the 'address' portion of NCQ. These findings will be invaluable to develop our next generation of MOR bitopic modulators with high affinity and subtype selectivity to potentially treat OUD.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors.

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the ß2 adrenergic (ß2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

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.

Advances in NMR Methods to Identify Allosteric Sites and Allosteric Ligands.

NMR allows assessment of protein structure in solution. Unlike conventional X-ray crystallography that provides snapshots of protein conformations, all conformational states are simultaneously accessible to analysis by NMR. This is a significant advantage for discovery and characterization of allosteric effects. These effects are observed when binding at one site of the protein affects another distinct site through conformational transitions. Allosteric regulation of proteins has been observed in multiple physiological processes in health and disease, providing an opportunity for the development of allosteric inhibitors. These compounds do not directly interact with the orthosteric site of the protein but influence its structure and function. In this book chapter, we provide an overview on how NMR methods are utilized to identify allosteric sites and to discover novel inhibitors, highlighting examples from the field. We also describe how NMR has contributed to understanding of allosteric mechanisms and propose that it is likely to play an important role in clarification and further development of key concepts of allostery.

Recent Insights from Molecular Dynamics Simulations for G Protein-Coupled Receptor Drug Discovery.

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. Some ligands that bind to GPCRs activate different downstream signaling pathways. G protein activation, or -arrestin biased signaling, involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. -arrestin biased signaling has become a hot target for structure-based drug discovery. However, challenges include that there are few crystal structures available in the Protein Data Bank and that GPCRs are highly dynamic. Hence, molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions with receptors and ligands. Here, we highlight recent MD simulation studies and enhanced sampling methods used to study biased G protein-coupled receptor signaling and their conformational dynamics as well as applications to drug discovery.

Covalent-Allosteric Inhibitors to Achieve Akt Isoform-Selectivity.

Isoforms of protein kinase Akt are involved in essential processes including cell proliferation, survival, and metabolism. However, their individual roles in health and disease have not been thoroughly evaluated. Thus, there is an urgent need for perturbation studies, preferably mediated by highly selective bioactive small molecules. Herein, we present a structure-guided approach for the design of structurally diverse and pharmacologically beneficial covalent-allosteric modifiers, which enabled an investigation of the isoform-specific preferences and the important residues within the allosteric site of the different isoforms. The biochemical, cellular, and structural evaluations revealed interactions responsible for the selective binding profiles. The isoform-selective covalent-allosteric Akt inhibitors that emerged from this approach showed a conclusive structure-activity relationship and broke ground in the development of selective probes to delineate the isoform-specific functions of Akt kinases.

Identification of new allosteric sites and modulators of AChE through computational and experimental tools.

Allosteric sites on proteins are targeted for designing more selective inhibitors of enzyme activity and to discover new functions. Acetylcholinesterase (AChE), which is most widely known for the hydrolysis of the neurotransmitter acetylcholine, has a peripheral allosteric subsite responsible for amyloidosis in Alzheimer's disease through interaction with amyloid ß-peptide. However, AChE plays other non-hydrolytic functions. Here, we identify and characterise using computational tools two new allosteric sites in AChE, which have allowed us to identify allosteric inhibitors by virtual screening guided by structure-based and fragment hotspot strategies. The identified compounds were also screened for in vitro inhibition of AChE and three were observed to be active. Further experimental (kinetic) and computational (molecular dynamics) studies have been performed to verify the allosteric activity. These new compounds may be valuable pharmacological tools in the study of non-cholinergic functions of AChE.

Nanosecond Dynamics Regulate the MIF-Induced Activity of CD74.

Macrophage migration inhibitory factor (MIF) activates CD74, which leads to severe disorders including inflammation, autoimmune diseases and cancer under pathological conditions. Molecular dynamics (MD) simulations up to one microsecond revealed dynamical correlation between a residue located at the opening of one end of the MIF solvent channel, previously thought to be a consequence of homotrimerization, and residues in a distal region responsible for CD74 activation. Experiments verified the allosteric regulatory site and identified a pathway to this site via the MIF ß-strands. The reported findings provide fundamental insights on a dynamic mechanism that controls the MIF-induced activation of CD74.

Binding hotspots on K-ras: consensus ligand binding sites and other reactive regions from probe-based molecular dynamics analysis.

We have used probe-based molecular dynamics (pMD) simulations to search for interaction hotspots on the surface of the therapeutically highly relevant oncogenic K-Ras G12D. Combining the probe-based query with an ensemble-based pocket identification scheme and an analysis of existing Ras-ligand complexes, we show that (i) pMD is a robust and cost-effective strategy for binding site identification, (ii) all four of the previously reported ligand binding sites are suitable for structure-based ligand design, and (iii) in some cases probe binding and expanded sampling of configurational space enable pocket expansion and increase the likelihood of site identification. Furthermore, by comparing the distribution of hotspots in nonpocket-like regions with known protein- and membrane-interacting interfaces, we propose that pMD has the potential to predict surface patches responsible for protein-biomolecule interactions. These observations have important implications for future drug design efforts and will facilitate the search for potential interfaces responsible for the proposed transient oligomerization or interaction of Ras with other biomolecules in the cellular milieu.

Harnessing allostery: a novel approach to drug discovery.

Allostery is the most direct and efficient way for regulation of biological macromolecule function, ranging from the control of metabolic mechanisms to signal transduction pathways. Allosteric modulators target to allosteric sites, offering distinct advantages compared to orthosteric ligands that target to active sites, such as greater specificity, reduced side effects, and lower toxicity. Allosteric modulators have therefore drawn increasing attention as potential therapeutic drugs in the design and development of new drugs. In recent years, advancements in our understanding of the fundamental principles underlying allostery, coupled with the exploitation of powerful techniques and methods in the field of allostery, provide unprecedented opportunities to discover allosteric proteins, detect and characterize allosteric sites, design and develop novel efficient allosteric drugs, and recapitulate the universal features of allosteric proteins and allosteric modulators. In the present review, we summarize the recent advances in the repertoire of allostery, with a particular focus on the aforementioned allosteric compounds.

Learning About Allosteric Drugs and Ways to Design Them.

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

Allosteric Binding Sites of the SARS-CoV-2 Main Protease: Potential Targets for Broad-Spectrum Anti-Coronavirus Agents.

The current pandemic caused by the COVID-19 disease has reached everywhere in the world and has affected every aspect of our lives. As of the current data, the World Health Organization (WHO) has reported more than 300 million confirmed COVID-19 cases worldwide and more than 5 million deaths. Mpro is an enzyme that plays a key role in the life cycle of the SARS-CoV-2 virus, and it is vital for the disease progression. The Mpro enzyme seems to have several allosteric sites that can hinder the enzyme catalytic activity. Furthermore, some of these allosteric sites are located at or nearby the dimerization interface which is essential for the overall Mpro activity. In this review paper, we investigate the potential of the Mpro allosteric site to act as a drug target, especially since they interestingly appear to be resistant to mutation. The work is illustrated through three subsequent sections: First, the two main categories of Mpro allosteric sites have been explained and discussed. Second, a total of six pockets have been studied and evaluated for their druggability and cavity characteristics. Third, the experimental and computational attempts for the discovery of new allosteric inhibitors have been illustrated and discussed. To sum up, this review paper gives a detailed insight into the feasibility of developing new Mpro inhibitors to act as a potential treatment for the COVID-19 disease.

Mapping the binding sites of challenging drug targets.

An increasing number of medically important proteins are challenging drug targets because their binding sites are too shallow or too polar, are cryptic and thus not detectable without a bound ligand or located in a protein-protein interface. While such proteins may not bind druglike small molecules with sufficiently high affinity, they are frequently druggable using novel therapeutic modalities. The need for such modalities can be determined by experimental or computational fragment based methods. Computational mapping by mixed solvent molecular dynamics simulations or the FTMap server can be used to determine binding hot spots. The strength and location of the hot spots provide very useful information for selecting potentially successful approaches to drug discovery.

Mechanism of allosteric activation of SIRT6 revealed by the action of rationally designed activators.

The recent discovery of activator compounds binding to an allosteric site on the NAD+-dependent protein lysine deacetylase, sirtuin 6 (SIRT6) has attracted interest and presents a pharmaceutical target for aging-related and cancer diseases. However, the mechanism underlying allosteric activation of SIRT6 by the activator MDL-801 remains largely elusive because no major conformational changes are observed upon activator binding. By combining molecular dynamics simulations with biochemical and kinetic analyses of wild-type SIRT6 and its variant M136A, we show that conformational rotation of 2-methyl-4-fluoro-5-bromo substituent on the right phenyl ring (R-ring) of MDL-801, which uncovers previously unseen hydrophobic interactions, contributes to increased activating deacetylation activity of SIRT6. This hypothesis is further supported by the two newly synthesized MDL-801 derivatives through the removal of the 5-Br atom on the R-ring (MDL-801-D1) or the restraint of the rotation of the R-ring (MDL-801-D2). We further propose that the 5-Br atom serves as an allosteric driver that controls the ligand allosteric efficacy. Our study highlights the effect of allosteric enzyme catalytic activity by activator binding and provides a rational approach for enhancing deacetylation activity.

Big data and artificial intelligence discover novel drugs targeting proteins without 3D structure and overcome the undruggable targets.

The discovery of targeted drugs heavily relies on three-dimensional (3D) structures of target proteins. When the 3D structure of a protein target is unknown, it is very difficult to design its corresponding targeted drugs. Although the 3D structures of some proteins (the so-called undruggable targets) are known, their targeted drugs are still absent. As increasing crystal/cryogenic electron microscopy structures are deposited in Protein Data Bank, it is much more possible to discover the targeted drugs. Moreover, it is also highly probable to turn previous undruggable targets into druggable ones when we identify their hidden allosteric sites. In this review, we focus on the currently available advanced methods for the discovery of novel compounds targeting proteins without 3D structure and how to turn undruggable targets into druggable ones.

Both intra and inter-domain interactions define the intrinsic dynamics and allosteric mechanism in DNMT1s.

DNA methyltransferase 1 (DNMT1), a large multidomain enzyme, is believed to be involved in the passive transmission of genomic methylation patterns via methylation maintenance. Yet, the molecular mechanism of interaction networks underlying DNMT1 structures, dynamics, and its biological significance has yet to be fully characterized. In this work, we used an integrated computational strategy that combined coarse-grained and atomistic simulations with coevolution information and network modeling of the residue interactions for the systematic investigation of allosteric dynamics in DNMT1. The elastic network modeling has proposed that the high plasticity of RFTS has strengthened the correlated behaviors of DNMT1 structures through the hinge sites located at the RFTS-CD interface, which mediate the collective motions between domains. The perturbation response scanning (PRS) analysis combined with the enrichment analysis of disease mutations have further highlighted the allosteric potential of the RFTS domain. Furthermore, the long-range paths connect the intra-domain interactions through the TRD interface and catalytic interface, emphasizing some key inter-domain interactions as the bridges in the global allosteric regulation of DNMT1. The observed interplay between conserved intra-domain networks and dynamical plasticity encoded by inter-domain interactions provides insights into the intrinsic dynamics and functional evolution, as well as the design of allosteric modulators of DNMT1 based on the TRD interface.