The different ways through which specificity works in orthosteric and allosteric drugs

16R-HETE and 16S-HETE alter human cytochrome P450 1B1 enzyme activity probably through an allosteric mechanism

Cytochrome P450 1B1 (CYP1B1) has been widely associated with the development of cardiac pathologies due to its ability to produce cardiotoxic metabolites like midchain hydroxyeicosatetraenoic acids (HETEs) from arachidonic acid (AA) through an allylic oxidation reaction. 16-HETE is a subterminal HETE that is also produced by CYP-mediated AA metabolism. 19-HETE is another subterminal HETE that was found to inhibit CYP1B1 activity, lower midchain HETEs, and have cardioprotective effects. However, the effect of 16-HETE enantiomers on CYP1B1 has not yet been investigated. We hypothesized that 16(R/S)-HETE could alter the activity of CYP1B1 and other CYP enzymes. Therefore, this study was carried out to investigate the modulatory effect of 16-HETE enantiomers on CYP1B1 enzyme activity, and to examine the mechanisms by which they exert these modulatory effects. To investigate whether these effects are specific to CYP1B1, we also investigated 16-HETE modulatory effects on CYP1A2. Our results showed that 16-HETE enantiomers significantly increased CYP1B1 activity in RL-14 cells, recombinant human CYP1B1, and human liver microsomes, as seen by the significant increase in 7-ethoxyresorufin deethylation rate. On the contrary, 16-HETE enantiomers significantly inhibited CYP1A2 catalytic activity mediated by the recombinant human CYP1A2 and human liver microsomes. 16R-HETE showed stronger effects than 16S-HETE. The sigmoidal binding mode of the enzyme kinetics data demonstrated that CYP1B1 activation and CYP1A2 inhibition occurred through allosteric regulation. In conclusion, our study provides the first evidence that 16R-HETE and 16S-HETE increase CYP1B1 catalytic activity through an allosteric mechanism.

mGlu2 and mGlu3 receptor negative allosteric modulators attenuate the interoceptive effects of alcohol in male and female rats

RATIONALE

The subjective effects of alcohol are associated with alcohol use disorder (AUD) vulnerability and treatment outcomes. The interoceptive effects of alcohol are part of these subjective effects and can be measured in animal models using drug discrimination procedures. The newly developed mGlu2 and mGlu3 negative allosteric modulators (NAMs) are potential therapeutics for AUD and may alter interoceptive sensitivity to alcohol.

OBJECTIVES

To determine the effects of mGlu2 and mGlu3 NAMs on the interoceptive effects of alcohol in rats.

METHODS

Long-Evans rats were trained to discriminate the interoceptive stimulus effects of alcohol (2.0 g/kg, i.g.) from water using both operant (males only) and Pavlovian (male and female) drug discrimination techniques. Following acquisition training, an alcohol dose-response (0, 0.5, 1.0, 2.0 g/kg) experiment was conducted to confirm stimulus control over behavior. Next, to test the involvement of mGlu2 and mGlu3, rats were pretreated with the mGlu2-NAM (VU6001966; 0, 3, 6, 12 mg/kg, i.p.) or the mGlu3-NAM (VU6010572; 0, 3, 6, 12 mg/kg, i.p.) before alcohol administration (2.0 g/kg, i.g.).

RESULTS

In Pavlovian discrimination, male rats showed greater interoceptive sensitivity to 1.0 and 2.0 g/kg alcohol compared to female rats. Both mGlu2-NAM and mGlu3-NAM attenuated the interoceptive effects of alcohol in male and female rats using Pavlovian and operant discrimination. There may be a potential sex difference in response to the mGlu2-NAM at the highest dose tested.

CONCLUSIONS

Male rats may be more sensitive to the interoceptive effects of the 2.0 g/kg alcohol training dose compared to female rats. Both mGlu2-and mGlu3-NAM attenuate the interoceptive effects of alcohol in male and female rats. These drugs may have potential for treatment of AUD in part by blunting the subjective effects of alcohol.

Secondary Sites of the C-type Lectin-Like Fold

C-type lectins are a large superfamily of proteins involved in a multitude of biological processes. In particular, their involvement in immunity and homeostasis has rendered them attractive targets for diverse therapeutic interventions. They share a characteristic C-type lectin-like domain whose adaptability enables them to bind a broad spectrum of ligands beyond the originally defined canonical Ca2+-dependent carbohydrate binding. Together with variable domain architecture and high-level conformational plasticity, this enables C-type lectins to meet diverse functional demands. Secondary sites provide another layer of regulation and are often intricately linked to functional diversity. Located remote from the canonical primary binding site, secondary sites can accommodate ligands with other physicochemical properties and alter protein dynamics, thus enhancing selectivity and enabling fine-tuning of the biological response. In this review, we outline the structural determinants allowing C-type lectins to perform a large variety of tasks and to accommodate the ligands associated with it. Using the six well-characterized Ca2+-dependent and Ca2+-independent C-type lectin receptors DC-SIGN, langerin, MGL, dectin-1, CLEC-2 and NKG2D as examples, we focus on the characteristics of non-canonical interactions and secondary sites and their potential use in drug discovery endeavors.

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.

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

Allosteric drugs: New principles and design approaches

Focusing on an important biomedical implication of allostery - design of allosteric drugs, we describe characteristics of allosteric sites, effectors, and their modes of actions distinguishing them from the orthosteric counterparts and calling for new principles and protocols in the quests for allosteric drugs. We show the importance of considering both binding affinity and allosteric signaling in establishing the structure-activity relationships (SARs) toward design of allosteric effectors, arguing that pairs of allosteric sites and their effector ligands - the site-effector pairs - should be generated and adjusted simultaneously in the framework of what we call directed design protocol. Key ideas and approaches for designing allosteric effectors including reverse perturbation, targeted and agnostic analysis are also discussed here. Several promising computational approaches are highlighted, along with the need for and potential advantages of utilizing generative models to facilitate discovery/design of new allosteric drugs.

A journey from molecule to physiology and in silico tools for drug discovery targeting the transient receptor potential vanilloid type 1 (TRPV1) channel

The heat and capsaicin receptor TRPV1 channel is widely expressed in nerve terminals of dorsal root ganglia (DRGs) and trigeminal ganglia innervating the body and face, respectively, as well as in other tissues and organs including central nervous system. The TRPV1 channel is a versatile receptor that detects harmful heat, pain, and various internal and external ligands. Hence, it operates as a polymodal sensory channel. Many pathological conditions including neuroinflammation, cancer, psychiatric disorders, and pathological pain, are linked to the abnormal functioning of the TRPV1 in peripheral tissues. Intense biomedical research is underway to discover compounds that can modulate the channel and provide pain relief. The molecular mechanisms underlying temperature sensing remain largely unknown, although they are closely linked to pain transduction. Prolonged exposure to capsaicin generates analgesia, hence numerous capsaicin analogs have been developed to discover efficient analgesics for pain relief. The emergence of in silico tools offered significant techniques for molecular modeling and machine learning algorithms to indentify druggable sites in the channel and for repositioning of current drugs aimed at TRPV1. Here we recapitulate the physiological and pathophysiological functions of the TRPV1 channel, including structural models obtained through cryo-EM, pharmacological compounds tested on TRPV1, and the in silico tools for drug discovery and repositioning.

ASD2023: towards the integrating landscapes of allosteric knowledgebase

Allosteric regulation, induced by perturbations at an allosteric site topographically distinct from the orthosteric site, is one of the most direct and efficient ways to fine-tune macromolecular function. The Allosteric Database (ASD; accessible online at http://mdl.shsmu.edu.cn/ASD) has been systematically developed since 2009 to provide comprehensive information on allosteric regulation. In recent years, allostery has seen sustained growth and wide-ranging applications in life sciences, from basic research to new therapeutics development, while also elucidating emerging obstacles across allosteric research stages. To overcome these challenges and maintain high-quality data center services, novel features were curated in the ASD2023 update: (i) 66 589 potential allosteric sites, covering > 80% of the human proteome and constituting the human allosteric pocketome; (ii) 748 allosteric protein-protein interaction (PPI) modulators with clear mechanisms, aiding protein machine studies and PPI-targeted drug discovery; (iii) 'Allosteric Hit-to-Lead,' a pioneering dataset providing panoramic views from 87 well-defined allosteric hits to 6565 leads and (iv) 456 dualsteric modulators for exploring the simultaneous regulation of allosteric and orthosteric sites. Meanwhile, ASD2023 maintains a significant growth of foundational allosteric data. Based on these efforts, the allosteric knowledgebase is progressively evolving towards an integrated landscape, facilitating advancements in allosteric target identification, mechanistic exploration and drug discovery.

Characterization of pH-induced conformational changes in recombinant DENV NS2B-NS3pro

The increasing frequency of Dengue is a cause of severe epidemics and therefore demands strategies for effective prevention, diagnosis, and treatment. DENV-protease is being investigated as a potential therapeutic target. However, due to the flat and highly charged active site of the DENV-protease, designing orthosteric medicines is very difficult. In this study, we have done a thorough analysis of pH-dependent conformational changes in recombinantly expressed DENV protease using various spectroscopic techniques. Our spectroscopic study of DENV protease (NS2B-NS3pro) at different pH conditions gives important insights into the dynamicity of structural conformation. At physiological pH, the DENV-protease exists in a random-coiled state. Lowering the pH promotes the formation of alpha-helical and beta-sheet structures i.e. gain of secondary structure as shown by Far-UV CD. The light scattering and Thioflavin T (ThT)-binding assay proved the aggregation-prone tendency of DENV-protease at pH 4.0. Further, the confocal microscopy image intensity showed the amorphous aggregate formation of DENV protease at pH 4.0. Thus, the DENV protease acquires different conformations with changes in pH conditions. Together, these results have the potential to facilitate the design of a conformation destabilizer-based therapeutic strategy for dengue fever.

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.

Protein conformational ensembles in function: roles and mechanisms

The sequence-structure-function paradigm has dominated twentieth century molecular biology. The paradigm tacitly stipulated that for each sequence there exists a single, well-organized protein structure. Yet, to sustain cell life, function requires (i) that there be more than a single structure, (ii) that there be switching between the structures, and (iii) that the structures be incompletely organized. These fundamental tenets called for an updated sequence-conformational ensemble-function paradigm. The powerful energy landscape idea, which is the foundation of modernized molecular biology, imported the conformational ensemble framework from physics and chemistry. This framework embraces the recognition that proteins are dynamic and are always interconverting between conformational states with varying energies. The more stable the conformation the more populated it is. The changes in the populations of the states are required for cell life. As an example, in vivo, under physiological conditions, wild type kinases commonly populate their more stable "closed", inactive, conformations. However, there are minor populations of the "open", ligand-free states. Upon their stabilization, e.g., by high affinity interactions or mutations, their ensembles shift to occupy the active states. Here we discuss the role of conformational propensities in function. We provide multiple examples of diverse systems, including protein kinases, lipid kinases, and Ras GTPases, discuss diverse conformational mechanisms, and provide a broad outlook on protein ensembles in the cell. We propose that the number of molecules in the active state (inactive for repressors), determine protein function, and that the dynamic, relative conformational propensities, rather than the rigid structures, are the hallmark of cell life.

Advances in adrenergic receptors for the treatment of chronic obstructive pulmonary disease: 2023 update

INTRODUCTION

Strong scientific evidence and large experience support the use of β2-agonists for the symptomatic alleviation of COPD. Therefore, there is considerable effort in discovering highly potent and selective β2-agonists.

AREAS COVERED

Recent research on novel β2-agonists for the treatment of COPD. A detailed literature search was performed in two major databases (PubMed/MEDLINE and Scopus) up to September 2023."

EXPERT OPINION

Compounds that preferentially activate a Gs- or β-arrestin-mediated signaling pathway via β- adrenoceptors (ARs) are more innovative. Pepducins, which target the intracellular region of β2-AR to modulate receptor signaling output, have the most interesting profile from a pharmacological point of view. They stabilize the conformation of the β2-AR and influence its signaling by interacting with the intracellular receptor-G protein interface. New bifunctional drugs called muscarinic antagonist-β2 agonist (MABA), which have both muscarinic receptor (mAChR) antagonism and β2-agonist activity in the same molecule, are a new opportunity. However, all tested compounds have been shown to act predominantly as mAChR antagonists or β2-agonists. An intriguing idea is to utilize allosteric modulators that bind to β2-ARs at sites different than those bound by orthosteric ligands to augment or reduce the signaling transduced by the orthosteric ligand.

AlphaFold, allosteric, and orthosteric drug discovery: Ways forward

Drug discovery is arguably a highly challenging and significant interdisciplinary aim. The stunning success of the artificial intelligence-powered AlphaFold, whose latest version is buttressed by an innovative machine-learning approach that integrates physical and biological knowledge about protein structures, raised drug discovery hopes that unsurprisingly, have not come to bear. Even though accurate, the models are rigid, including the drug pockets. AlphaFold's mixed performance poses the question of how its power can be harnessed in drug discovery. Here we discuss possible ways of going forward wielding its strengths, while bearing in mind what AlphaFold can and cannot do. For kinases and receptors, an input enriched in active (ON) state models can better AlphaFold's chance of rational drug design success.

Protein-mediated interactions in the dynamic regulation of acute inflammation

Protein-mediated interactions are the fundamental mechanism through which cells regulate health and disease. These interactions require physical contact between proteins and their respective targets of interest. These targets include not only other proteins but also nucleic acids and other important molecules as well. These proteins are often involved in multibody complexes that work dynamically to regulate cellular health and function. Various techniques have been adapted to study these important interactions, such as affinity-based assays, mass spectrometry, and fluorescent detection. The application of these techniques has led to a greater understanding of how protein interactions are responsible for both the instigation and resolution of acute inflammatory diseases. These pursuits aim to provide opportunities to target specific protein interactions to alleviate acute inflammation.

Enhancement of Inhibitory Activity by Combining Allosteric Inhibitors Putatively Binding to Different Allosteric Sites on Cathepsin K

BACKGROUND

Cathepsin K, which is involved in bone resorption, is a good target for treating osteoporosis, but no clinically approved medicine has been developed. Recently, allosteric inhibitors with high specificity and few side effects have been attracting attention for use in new medicines.

METHODS

Cathepsin K inhibitors were isolated from the methanol extract of Chamaecrista nomame (Leguminosae) using cathepsin K inhibition activity-assisted multi-step chromatography. Standard kinetic analysis was employed to examine the mechanism of cathepsin K inhibition when an isolated inhibitor and its derivative were used. The allosteric binding of these cathepsin K inhibitors was supported by a docking study using AutoDock vina. Combinations of allosteric cathepsin K inhibitors expected to bind to different allosteric sites were examined by means of cathepsin K inhibition assay.

RESULTS

Two types of cathepsin K inhibitors were identified in the methanol extract of Chamaecrista nomame. One type consisted of cassiaoccidentalin B and torachrysone 8-β-gentiobioside, and inhibited both cathepsin K and B with similar inhibitory potential, while the other type of inhibitor consisted of pheophytin a, and inhibited cathepsin K but not cathepsin B, suggesting that pheophytin a binds to an allosteric site of cathepsin K. Kinetic analysis of inhibitory activity suggested that pheophytin a and its derivative, pheophorbide b, bind allosterically to cathepsin K. This possibility was supported by a docking study on cathepsin K. The cathepsin K inhibitory activity of pheophytin a and pheophorbide b was enhanced by combining them with the allosteric inhibitors NSC 13345 and NSC94914, which bind to other allosteric sites on cathepsin K.

CONCLUSIONS

Different allosteric inhibitors that bind to different sites in combination, as shown in this study, may be useful for designing new allosteric inhibitory drugs with high specificity and few side effects.

Stability Analysis of the Asiatic Acid-COX-2 Complex Using 100 ns Molecular Dynamic Simulations and Its Selectivity against COX-2 as a Potential Anti-Inflammatory Candidate

Asiatic acid, a triterpenoid compound, has been shown to have anti-inflammatory activity through the inhibition of the formation of cyclooxygenase-2 (COX-2) in vitro and in vivo. This study was conducted to determine the binding stability and the inhibitory potential of asiatic acid as an anti-inflammatory candidate. The study involved in vitro testing utilizing a colorimetric kit as well as in silico testing for the pharmacophore modeling and molecular dynamic (MD) simulation of asiatic acid against COX-2 (PDB ID: 3NT1). The MD simulations showed a stable binding of asiatic acid to COX-2 and an RMSD range of 1-1.5 Å with fluctuations at the residues of Phe41, Leu42, Ile45, Arg44, Asp367, Val550, Glu366, His246, and Gly227. The total binding energy of the asiatic acid-COX-2 complex is -7.371 kcal/mol. The anti-inflammatory activity of the asiatic acid inhibition of COX-2 was detected at IC₅₀ values of 120.17 µM. Based on pharmacophore modeling, we discovered that carboxylate and hydroxyl are the two main functional groups that act as hydrogen bond donors and acceptors interacting with the COX-2 enzyme. From the results, it is evident that asiatic acid is a potential anti-inflammatory candidate with high inhibitory activity in relation to the COX-2 enzyme.

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

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

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.

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.

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.

Molecular Pharmacology of P2X Receptors: Exploring Druggable Domains Revealed by Structural Biology

Extracellular ATP is a critical signaling molecule that is found in a wide range of concentrations across cellular environments. The family of nonselective cation channels that sense extracellular ATP, termed P2X receptors (P2XRs), is composed of seven subtypes (P2X₁-P2X₇) that assemble as functional homotrimeric and heterotrimeric ion channels. Each P2XR is activated by a distinct concentration of extracellular ATP, spanning from high nanomolar to low millimolar. P2XRs are implicated in a variety of physiological and pathophysiological processes in the cardiovascular, immune, and central nervous systems, corresponding to the spatiotemporal expression, regulation, and activation of each subtype. The therapeutic potential of P2XRs is an emerging area of research in which structural biology has seemingly exceeded medicinal chemistry, as there are several published P2XR structures but currently no FDA-approved drugs targeting these ion channels. Cryogenic electron microscopy is ideally suited to facilitate structure-based drug design for P2XRs by revealing and characterizing novel ligand-binding sites. This review covers structural elements in P2XRs including the extracellular orthosteric ATP-binding site, extracellular allosteric modulator sites, channel pore, and cytoplasmic substructures, with an emphasis on potential therapeutic ligand development.

Pharmacotherapy Based on ACE2 Targeting and COVID-19 Infection

The new SARS-CoV-2 coronavirus is responsible for the COVID-19 pandemic. A massive vaccination campaign, which is still ongoing, has averted most serious consequences worldwide; however, lines of research are continuing to identify the best drug therapies to treat COVID-19 infection. SARS-CoV-2 penetrates the cells of the host organism through ACE2. The ACE2 protein plays a key role in the renin–angiotensin system (RAS) and undergoes changes in expression during different stages of COVID-19 infection. It appears that an unregulated RAS is responsible for the severe lung damage that occurs in some cases of COVID-19. Pharmacologically modifying the expression of ACE2 could be an interesting line of research to follow in order to avoid the severe complications of COVID-19.

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.

Allosteric regulation of the soluble epoxide hydrolase by nitro fatty acids using a combined experimental and computational approach

The human soluble epoxide hydrolase (hsEH) is a key regulator of epoxy fatty acid (EpFA) metabolism. Inhibition of sEH can maintain endogenous levels of beneficial EpFAs and reduce the levels of their corresponding diol products, thus ameliorating a variety of pathological conditions including cardiovascular, central nervous system and metabolic diseases. The quest for orthosteric drugs that bind directly to the catalytic crevice of hsEH has been prolonged and sustained over the past decades, but the disappointing outcome of clinical trials to date warrants alternative pharmacological approaches. Previously, we have shown that hsEH can be allosterically inhibited by the endogenous electrophilic lipid 15-deoxy-Δ^12,14-Prostaglandin-J₂, via covalent adduction to two cysteines, C423 and C522. In this study, we explore the properties and behaviour of three electrophilic lipids belonging to the class of the nitro fatty acids, namely 9- and 10-nitrooleate and 10-nitrolinoleate. Biochemical and biophysical investigations revealed that, in addition to C423 and C522, nitro fatty acids can covalently bind to additional nucleophilic residues in hsEH C-terminal domain (CTD), two of which predicted in this study to be latent allosteric sites. Systematic mapping of the protein mutational space and evaluation of possible propagation pathways delineated selected residues, both in the allosteric patches and in other regions of the enzyme, envisaged to play a role on allosteric signalling. The responses elicited by the ligands on the covalent adduction sites supports future fragment-based design studies of new allosteric effectors for hsEH with increased efficacy and selectivity.

Outliers in SAR and QSAR: 4. effects of allosteric protein-ligand interactions on the classical quantitative structure-activity relationships

Effects of allosteric interactions on the classical structure-activity relationship (SAR) and quantitative SAR (QSAR) have been investigated. Apprehending the outliers in SAR and QSAR studies can improve the quality, predictability, and use of QSAR in designing unknown compounds in drug discovery research. We explored allosteric protein-ligand interactions as a possible source of outliers in SAR/QSAR. We used glycogen phosphorylase as an example of a protein that has an allosteric site. Examination of the ligand-bound x-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site. The case provided an example of constructive use of QSAR identifying outliers with alternative binding modes. Other allosteric QSARs that captured our attention were the inverted parabola/bilinear QSARs. The x-ray crystal structures and the QSAR analyses indicated that the inverted parabola QSARs could be associated with the conformational changes in the allosteric interactions. Our results showed that the normal parabola, as well as the inverted parabola QSARs, can describe the allosteric interactions. Examination of the ligand-bound X-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site.

Allosteric Site of ACE-2 as a Drug Target for COVID-19

The coronavirus disease 2019 (COVID-19) pandemic has a significant impact on healthcare systems and our lives. Vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) provide protection against SARS-CoV-2. However, mutations in the viral genome are common, raising concerns about the effectiveness of existing vaccines for SARS-CoV-2. The receptor-binding domain (RBD) of SARS-CoV-2 uses angiotensin-converting enzyme-2 (ACE-2) as a gateway to enter host cells. Therefore, the ACE-2-RBD interaction may be targeted by antiviral drugs. In this context, allosteric modulation of ACE-2 may offer a promising approach. It may lead to allosteric inhibition of the interaction between ACE-2 and SARS-CoV-2.

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.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Dynamic Cellular Cartography: Mapping the Local Determinants of Oligodendrocyte Transcription Factor 2 (OLIG2) Function in Live Cells Using Massively Parallel Fluorescence Correlation Spectroscopy Integrated with Fluorescence Lifetime Imaging Microscopy (mpFCS/FLIM)

Compartmentalization and integration of molecular processes through diffusion are basic mechanisms through which cells perform biological functions. To characterize these mechanisms in live cells, quantitative and ultrasensitive analytical methods with high spatial and temporal resolution are needed. Here, we present quantitative scanning-free confocal microscopy with single-molecule sensitivity, high temporal resolution (∼10 μs/frame), and fluorescence lifetime imaging capacity, developed by integrating massively parallel fluorescence correlation spectroscopy with fluorescence lifetime imaging microscopy (mpFCS/FLIM); we validate the method, use it to map in live cell location-specific variations in the concentration, diffusion, homodimerization, DNA binding, and local environment of the oligodendrocyte transcription factor 2 fused with the enhanced Green Fluorescent Protein (OLIG2-eGFP), and characterize the effects of an allosteric inhibitor of OLIG2 dimerization on these determinants of OLIG2 function. In particular, we show that cytoplasmic OLIG2-eGFP is largely monomeric and freely diffusing, with the fraction of freely diffusing OLIG2-eGFP molecules being f_D,free^cyt = (0.75 ± 0.10) and the diffusion time τ_D,free^cyt = (0.5 ± 0.3) ms. In contrast, OLIG2-eGFP homodimers are abundant in the cell nucleus, constituting ∼25% of the nuclear pool, some f_D,bound^nuc = (0.65 ± 0.10) of nuclear OLIG2-eGFP is bound to chromatin DNA, whereas freely moving OLIG2-eGFP molecules diffuse at the same rate as those in the cytoplasm, as evident from the lateral diffusion times τ_D,free^nuc = τ_D,free^cyt = (0.5 ± 0.3) ms. OLIG2-eGFP interactions with chromatin DNA, revealed through their influence on the apparent diffusion behavior of OLIG2-eGFP, τ_D,bound^nuc (850 ± 500) ms, are characterized by an apparent dissociation constant K_d,app^OLIG2-DNA = (45 ± 30) nM. The apparent dissociation constant of OLIG2-eGFP homodimers was estimated to be K_d,app^{(OLIG2-eGFP)2} ≈ 560 nM. The allosteric inhibitor of OLIG2 dimerization, compound NSC 50467, neither affects OLIG2-eGFP properties in the cytoplasm nor does it alter the overall cytoplasmic environment. In contrast, it significantly impedes OLIG2-eGFP homodimerization in the cell nucleus, increasing five-fold the apparent dissociation constant, K_{d,app,NSC50467}^{(OLIG2-eGFP)2} ≈ 3 μM, thus reducing homodimer levels to below 7% and effectively abolishing OLIG2-eGFP specific binding to chromatin DNA. The mpFCS/FLIM methodology has a myriad of applications in biomedical research and pharmaceutical industry. For example, it is indispensable for understanding how biological functions emerge through the dynamic integration of location-specific molecular processes and invaluable for drug development, as it allows us to quantitatively characterize the interactions of drugs with drug targets in live cells.

Targeting the mitochondrial chaperone TRAP1: strategies and therapeutic perspectives

TRAP1, the mitochondrial isoform of heat shock protein (Hsp)90 chaperones, is a key regulator of metabolism and organelle homeostasis in diverse pathological states. While selective TRAP1 targeting is an attractive goal, classical active-site-directed strategies have proved difficult, due to high active site conservation among Hsp90 paralogs. Here, we discuss advances in developing TRAP1-directed strategies, from lead modification with mitochondria delivery groups to the computational discovery of allosteric sites and ligands. Specifically, we address the unique opportunities that targeting TRAP1 opens up in tackling fundamental questions on its biology and in unveiling new therapeutic approaches. Finally, we show how crucial to this endeavor is our ability to predict the activities of TRAP1-selective allosteric ligands and to optimize target engagement to avoid side effects.

Inhibition of Nonfunctional Ras

Intuitively, functional states should be targeted; not nonfunctional ones. So why could drugging the inactive K-Ras4BG12Cwork-but drugging the inactive kinase will likely not? The reason is the distinct oncogenic mechanisms. Kinase driver mutations work by stabilizing the active state and/or destabilizing the inactive state. Either way, oncogenic kinases are mostly in the active state. Ras driver mutations work by quelling its deactivation mechanisms, GTP hydrolysis, and nucleotide exchange. Covalent inhibitors that bind to the inactive GDP-bound K-Ras4BG12C conformation can thus work. By contrast, in kinases, allosteric inhibitors work by altering the active-site conformation to favor orthosteric drugs. From the translational standpoint this distinction is vital: it expedites effective pharmaceutical development and extends the drug classification based on the mechanism of action. Collectively, here we postulate that drug action relates to blocking the mechanism of activation, not to whether the protein is in the active or inactive state.

Deciphering the T790M/L858R-Selective Inhibition Mechanism of an Allosteric Inhibitor of EGFR: Insights from Molecular Simulations

Allosteric inhibitors have lately received great attention because of their unique advantages, representing a more suitable choice for combinatory therapeutics targeting resistance-relevant signaling cascades. Among the various inhibitors, an allosteric small-molecule inhibitor, JBJ-04-125-02, has been proven to be effective against EGFR^T790M/L858R mutant in vivo and in vitro. Herein, an in silico approach was adopted to shed light on the deep understanding of the higher selectivity of JBJ-04-125-02 against EGFR^T790M/L858R mutant than wild-type EGFR. Our results indicate that JBJ-04-125-02 prefers to bind with the EGFR^T790M/L858R mutant, stabilizes the inactive conformation, and further allosterically affects the conformations and dynamics of the interlobe cleft, including both the allosteric site and the ATP-binding site. Furthermore, docking results confirm that the binding of JBJ-04-125-02 at the allosteric site decreases the binding affinity of ANP (an ATP analogue) at the orthosteric site, especially for the Mut-holo one, which might further inhibit the function of EGFR. The present work provides a clear picture of the mutant-selective inhibition mechanism of an allosteric inhibitor of EGFR. The findings might pave the way for designing allosteric drugs targeting EGFR mutant lung cancer patients, which also takes a step forward in terms of drug resistance caused by protein mutations.

The Discovery of a Putative Allosteric Site in the SARS-CoV-2 Spike Protein Using an Integrated Structural/Dynamic Approach

SARS-CoV-2 has caused the largest pandemic of the twenty-first century (COVID-19), threatening the life and economy of all countries in the world. The identification of novel therapies and vaccines that can mitigate or control this global health threat is among the most important challenges facing biomedical sciences. To construct a long-term strategy to fight both SARS-CoV-2 and other possible future threats from coronaviruses, it is critical to understand the molecular mechanisms underlying the virus action. The viral entry and associated infectivity stems from the formation of the SARS-CoV-2 spike protein complex with angiotensin-converting enzyme 2 (ACE2). The detection of putative allosteric sites on the viral spike protein molecule can be used to elucidate the molecular pathways that can be targeted with allosteric drugs to weaken the spike-ACE2 interaction and, thus, reduce viral infectivity. In this study, we present the results of the application of different computational methods aimed at detecting allosteric sites on the SARS-CoV-2 spike protein. The adopted tools consisted of the protein contact networks (PCNs), SEPAS (Affinity by Flexibility), and perturbation response scanning (PRS) based on elastic network modes. All of these methods were applied to the ACE2 complex with both the SARS-CoV2 and SARS-CoV spike proteins. All of the adopted analyses converged toward a specific region (allosteric modulation region [AMR]), present in both complexes and predicted to act as an allosteric site modulating the binding of the spike protein with ACE2. Preliminary results on hepcidin (a molecule with strong structural and sequence with AMR) indicated an inhibitory effect on the binding affinity of the spike protein toward the ACE2 protein.

G-Protein coupled receptors: structure and function in drug discovery

The G-protein coupled receptors (GPCRs) superfamily comprise similar proteins arranged into families or classes thus making it one of the largest in the mammalian genome. GPCRs take part in many vital physiological functions making them targets for numerous novel drugs. GPCRs share some distinctive features, such as the seven transmembrane domains, they also differ in the number of conserved residues in their transmembrane domain. Here we provide an introductory and accessible review detailing the computational advances in GPCR pharmacology and drug discovery. An overview is provided on family A-C GPCRs; their structural differences, GPCR signalling, allosteric binding and cooperativity. The dielectric constant (relative permittivity) of proteins is also discussed in the context of site-specific environmental effects.

Effect of sodium butyrate on HDAC8 mRNA expression in colorectal cancer cell lines and molecular docking study of LHX1 - sodium butyrate interaction

Colorectal cancer (CRC) is the third most common type of cancer and the fourth leading cause of cancer related deaths worldwide. The Histone Deacetylase 8 (HDAC8) gene is a gene with unique features which can be used as a potential target for drug design. The LHX1 transcription factor is an important transcription factor for this gene. The aim of this study was to investigate the effect of sodium butyrate (NaB) as a histone deacetylase inhibitor (HDACi) on the expression of the HDAC8 gene in the colorectal cancer cell line, and the molecular docking of the LHX1 transcription factor with NaB. For this purpose, HCT-116 and HT-29 cell lines were treated with different concentrations of NaB (6.25 mM to 150 mM) at 24, 48 and 72 hours. Subsequently, RNA was extracted from the treated and untreated cells and cDNA was synthesized. Quantitative Real-Time-PCR was done to investigate the mRNA expression of HDAC8. Molecular docking was also performed to investigate the interaction between NaB and LHX1. Based on Real-time-PCR results, the concentration of 150 mM of NaB after 24 hours in HT-29 and HCT-116 cell lines caused a significant reduction in mRNA expression of HDAC8 (P<0.05). After 48 hours of treatment, there was a significant decrease in the mRNA expression of HDAC8 at all concentrations (P<0.05). The docking results showed that LHX1 and NaB interacted best at the lowest energy levels. Our results also showed that NaB bonded strongly to LHX1. In addition, our results demonstrated that NaB bound to the LHX1 transcription factor and inhibited the function of this factor and consequently decreased the transcription from the HDAC8 gene which resulted in cell death. Future studies are needed to assess the likely molecular mechanisms of NaB action on gene expression.

Regulation of Thrombin Activity with Ultrasmall Nanoparticles: Effects of Surface Chemistry

Nanomaterials displaying well-tailored sizes and surface chemistries can provide novel ways with which to modulate the structure and function of enzymes. Recently, we showed that gold nanoparticles (AuNPs) in the ultrasmall size regime could perform as allosteric effectors inducing partial inhibition of thrombin activity. We now find that the nature of the AuNP surface chemistry controls the interactions to the anion-binding exosites 1 and 2 on the surface of thrombin, the allosterically induced changes to the active-site conformation, and, by extension, the enzymatic activity. Ultrasmall AuNPs passivated with p-mercaptobenzoic acid ligands (AuMBA) and a peptide-based (Ac-ECYN) biomimetic coat (AuECYN) were utilized in our investigations. Remarkably, we found that while AuMBA binds to exosites 1 and 2, AuECYN interacts primarily with exosite 2. It was further established that AuMBA behaves as a "mild denaturant" of thrombin leading to catalytic dysfunction over time. Conversely, AuECYN resembles a proper allosteric effector leading to partial and reversible inhibition of the activity. Collectively, our findings reveal how the distinct binding modes of different AuNP types may uniquely influence thrombin structure and catalysis. The present study further contributes to our understanding of how synthetic nanomaterials could be exploited in the allosteric regulation of enzymes.

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Distinctive communication networks in inactive states of β2 -adrenergic receptor: Mutual information and entropy transfer analysis

Mutual information and entropy transfer analysis employed on two inactive states of human beta-2 adrenergic receptor (β₂ -AR) unraveled distinct communication pathways. Previously, a so-called "highly" inactive state of the receptor was observed during 1.5 microsecond long molecular dynamics simulation where the largest intracellular loop (ICL3) was swiftly packed onto the G-protein binding cavity, becoming entirely inaccessible. Mutual information quantifying the degree of correspondence between backbone-C_α fluctuations was mostly shared between intra- and extra-cellular loop regions in the original inactive state, but shifted to entirely different regions in this latest inactive state. Interestingly, the largest amount of mutual information was always shared among the mobile regions. Irrespective of the conformational state, polar residues always contributed more to mutual information than hydrophobic residues, and also the number of polar-polar residue pairs shared the highest degree of mutual information compared to those incorporating hydrophobic residues. Entropy transfer, quantifying the correspondence between backbone-C_α fluctuations at different timesteps, revealed a distinctive pathway directed from the extracellular site toward intracellular portions in this recently exposed inactive state for which the direction of information flow was the reverse of that observed in the original inactive state where the mobile ICL3 and its intracellular surroundings drove the future fluctuations of extracellular regions.

Unraveling allosteric landscapes of allosterome with ASD

Allosteric regulation is one of the most direct and efficient ways to fine-tune protein function; it is induced by the binding of a ligand at an allosteric site that is topographically distinct from an orthosteric site. The Allosteric Database (ASD, available online at http://mdl.shsmu.edu.cn/ASD) was developed ten years ago to provide comprehensive information related to allosteric regulation. In recent years, allosteric regulation has received great attention in biological research, bioengineering, and drug discovery, leading to the emergence of entire allosteric landscapes as allosteromes. To facilitate research from the perspective of the allosterome, in ASD 2019, novel features were curated as follows: (i) >10 000 potential allosteric sites of human proteins were deposited for allosteric drug discovery; (ii) 7 human allosterome maps, including protease and ion channel maps, were built to reveal allosteric evolution within families; (iii) 1312 somatic missense mutations at allosteric sites were collected from patient samples from 33 cancer types and (iv) 1493 pharmacophores extracted from allosteric sites were provided for modulator screening. Over the past ten years, the ASD has become a central resource for studying allosteric regulation and will play more important roles in both target identification and allosteric drug discovery in the future.

A kainate receptor-selective RNA aptamer

Kainate and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors are two major, closely related receptor subtypes in the glutamate ion channel family. Excessive activities of these receptors have been implicated in a number of central nervous system diseases. Designing potent and selective antagonists of these receptors, especially of kainate receptors, is useful for developing potential treatment strategies for these neurological diseases. Here, we report on two RNA aptamers designed to individually inhibit kainate and AMPA receptors. To improve the biostability of these aptamers, we also chemically modified these aptamers by substituting their 2'-OH group with 2'-fluorine. These 2'-fluoro aptamers, FB9s-b and FB9s-r, were markedly resistant to RNase-catalyzed degradation, with a half-life of ∼5 days in rat cerebrospinal fluid or serum-containing medium. Furthermore, FB9s-r blocked AMPA receptor activity. Aptamer FB9s-b selectively inhibited GluK1 and GluK2 kainate receptor subunits, and also GluK1/GluK5 and GluK2/GluK5 heteromeric kainate receptors with equal potency. This inhibitory profile makes FB9s-b a powerful template for developing tool molecules and drug candidates for treatment of neurological diseases involving excessive activities of the GluK1 and GluK2 subunits.

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.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

Allosteric inhibition of α-thrombin enzymatic activity with ultrasmall gold nanoparticles

The catalytic activity of enzymes can be regulated by interactions with synthetic nanoparticles (NPs) in a number of ways. To date, however, the potential use of NPs as allosteric effectors has not been investigated in detail. Importantly, targeting allosteric (distal) sites on the enzyme surface could afford unique ways to modulate the activity, allowing for either enzyme activation, partial or full inhibition. Using p-mercaptobenzoic acid-coated ultrasmall gold NPs (AuMBA) and human α-thrombin as a model system, here we experimentally tested the hypothesis that enzyme activity could be regulated through ultrasmall NP interactions at allosteric sites. We show that AuMBA interacted selectively and reversibly around two positively charged regions of the thrombin surface (exosites 1 and 2) and away from the active site. NP complexation at the exosites transmitted long-range structural changes over to the active site, altering both substrate binding affinity and catalysis. Significantly, thrombin activity was partially reduced - but not completely inhibited - by interactions with AuMBA. These findings indicate that interactions of proteins with ultrasmall NPs may mimic a typical biomolecular complexation event, and suggest the prospect of using ultrasmall particles as synthetic receptors to allosterically regulate protein function.