Druggability Simulations and X-Ray Crystallography Reveal a Ligand-Binding Site in the GluA3 AMPA Receptor N-Terminal Domain

Allo-targeting of the kinase domain: Insights from in silico studies and comparison with experiments

The eukaryotic protein kinase domain has been a broadly explored target for drug discovery, despite limitations imposed by its high sequence conservation as a shared modular domain and the development of resistance to drugs. One way of addressing those limitations has been to target its potential allosteric sites, shortly called allo-targeting, in conjunction with, or separately from, its conserved catalytic/orthosteric site that has been widely exploited. Allosteric regulation has gained importance as an alternative to overcome the drawbacks associated with the indiscriminate effect of targeting the active site, and it turned out to be particularly useful for these highly promiscuous and broadly shared kinase domains. Yet, allo-targeting often faces challenges as the allosteric sites are not as clearly defined as its orthosteric sites, and the effect on the protein function may not be unambiguously assessed. A robust understanding of the consequence of site-specific allo-targeting on the conformational dynamics of the target protein is essential to design effective allo-targeting strategies. Recent years have seen important advances in in silico identification of druggable sites and distinguishing among them those sites expected to allosterically mediate conformational switches essential to signal transmission. The present opinion underscores the utility of such computational approaches applied to the kinase domain, with the help of comparison between computational predictions and experimental observations.

Diversity of AMPA Receptor Ligands: Chemotypes, Binding Modes, Mechanisms of Action, and Therapeutic Effects

L-Glutamic acid is the main excitatory neurotransmitter in the central nervous system (CNS). Its associated receptors localized on neuronal and non-neuronal cells mediate rapid excitatory synaptic transmission in the CNS and regulate a wide range of processes in the brain, spinal cord, retina, and peripheral nervous system. In particular, the glutamate receptors selective to α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) also play an important role in numerous neurological disorders and attract close attention as targets for the creation of new classes of drugs for the treatment or substantial correction of a number of serious neurodegenerative and neuropsychiatric diseases. For this reason, the search for various types of AMPA receptor ligands and studies of their properties are attracting considerable attention both in academic institutions and in pharmaceutical companies around the world. This review focuses mainly on the advances in this area published since 2017. Particular attention is paid to the structural diversity of new chemotypes of agonists, competitive AMPA receptor antagonists, positive and negative allosteric modulators, transmembrane AMPA regulatory protein (TARP) dependent allosteric modulators, ion channel blockers as well as their binding sites. This review also presents the studies of the mechanisms of action of AMPA receptor ligands that mediate their therapeutic effects.

Protein dynamics developments for the large scale and cryoEM: case study of ProDy 2.0

New computational biophysics pipelines for analysing the global dynamics of structural ensembles and large, dynamic complexes resolved by cryoEM are reviewed.

Cryo-electron microscopy (cryoEM) has become a well established technique with the potential to produce structures of large and dynamic supramolecular complexes that are not amenable to traditional approaches for studying structure and dynamics. The size and low resolution of such molecular systems often make structural modelling and molecular dynamics simulations challenging and computationally expensive. This, together with the growing wealth of structural data arising from cryoEM and other structural biology methods, has driven a trend in the computational biophysics community towards the development of new pipelines for analysing global dynamics using coarse-grained models and methods. At the centre of this trend has been a return to elastic network models, normal mode analysis (NMA) and ensemble analyses such as principal component analysis, and the growth of hybrid simulation methodologies that make use of them. Here, this field is reviewed with a focus on ProDy, the Python application programming interface for protein dynamics, which has been developed over the last decade. Two key developments in this area are highlighted: (i) ensemble NMA towards extracting and comparing the signature dynamics of homologous structures, aided by the recent SignDy pipeline, and (ii) pseudoatom fitting for more efficient global dynamics analyses of large and low-resolution supramolecular assemblies from cryoEM, revisited in the CryoDy pipeline. It is believed that such a renewal and extension of old models and methods in new pipelines will be critical for driving the field forward into the next cryoEM revolution.

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.

AMPA receptor trafficking and LTP: Carboxy-termini, amino-termini and TARPs

AMPA receptors (AMPARs) are fundamental elements in excitatory synaptic transmission and synaptic plasticity in the CNS. Long term potentiation (LTP), a form of synaptic plasticity which contributes to learning and memory formation, relies on the accumulation of AMPARs at the postsynapse. This phenomenon requires the coordinated recruitment of different elements in the AMPAR complex. Based on recent research reviewed herein, we propose an updated AMPAR trafficking and LTP model which incorporates both extracellular as well as intracellular mechanisms. This article is part of the special Issue on 'Glutamate Receptors - AMPA receptors'.

[Ligands of the AMPA-subtype glutamate receptors: mechanisms of action and novel chemotypes]

Ionotropic glutamate receptors of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) subtype play a key role in synaptic plasticity representing one of the mechanisms for learning and memory formation. They can also serve as targets for the development of novel classes of pharmaceuticals for the treatment or substantive correction of many serious neurodegenerative and psychoneurological disorders. The search and studies of various types of AMPA receptor ligands attract considerable attention from academic organizations and pharmaceutical companies all over the world. This review mainly focuses on recent advances in this field. The architecture and operational mechanism of the receptor as well as its major binding sites and ligand types are considered. Special attention is paid to the studies of mechanisms of action and novel chemotypes of AMPA receptor agonists and competitive antagonists, positive and negative allosteric modulators, auxiliary protein-dependent allosteric modulators, and ion channel blockers.

Network-Based Target Prioritization and Drug Candidate Identification for Multiple Sclerosis: From Analyzing "Omics Data" to Druggability Simulations

Multiple sclerosis (MS) is the most common chronic inflammatory demyelinating disease of the central nervous system. While the drugs currently available for MS provide symptomatic benefit, there is no curative treatment. The emergence of large-scale multiomics data and network theory provide new opportunities for drug discovery in MS, as these are promising strategies for developing novel drugs. In this study, we proposed a computational framework that combined biomolecular network modeling and structural dynamics analysis to facilitate the discovery of new drugs with potential activity in MS. First, we developed a new shortest path-based algorithm that prioritized differentially expressed genes using a newly topological and functional exploration of protein-protein interaction network. Then, pathway enrichment analysis and an assessment of target druggability suggested that TNF-α-induced protein 3 (TNFAIP3), which is involved in NF-κ B signaling, could be a potential therapeutic target for MS. Finally, druggability simulations and mutation enrichment analysis of the TNFAIP3 dimer presented two druggable sites. Follow-up pharmacophore model-based virtual screening of the two sites yielded 30 hit compounds with low energy scores. In summary, this novel method based on analyzing "omics data" and performing druggability simulations, is a systematic approach that unravels disease mechanisms and links them to the chemical space to develop treatments and can be applied to other complex diseases.

A comprehensive binding study illustrates ligand recognition in the periplasmic binding protein PotF

Periplasmic binding proteins (PBPs) are ubiquitous receptors in gram-negative bacteria. They sense solutes and play key roles in nutrient uptake. Escherichia coli's putrescine receptor PotF has been reported to bind putrescine and spermidine. We reveal that several similar biogenic polyamines are recognized by PotF. Using isothermal titration calorimetry paired with X-ray crystallography of the different complexes, we unveil PotF's binding modes in detail. The binding site for PBPs is located between two lobes that undergo a large conformational change upon ligand recognition. Hence, analyzing the influence of ligands on complex formation is crucial. Therefore, we solved crystal structures of an open and closed apo state and used them as a basis for molecular dynamics simulations. In addition, we accessed structural behavior in solution for all complexes by ¹H-¹⁵N HSQC NMR spectroscopy. This combined analysis provides a robust framework for understanding ligand binding for future developments in drug design and protein engineering.

Essential site scanning analysis: A new approach for detecting sites that modulate the dispersion of protein global motions

Despite the wealth of methods developed for exploring the molecular basis of allostery in biomolecular systems, there is still a need for structure-based predictive tools that can efficiently detect susceptible sites for triggering allosteric responses. Toward this goal, we introduce here an elastic network model (ENM)-based method, Essential Site Scanning Analysis (ESSA). Essential sites are here defined as residues that would significantly alter the protein's global dynamics if bound to a ligand. To mimic the crowding induced upon substrate binding, the heavy atoms of each residue are incorporated as additional network nodes into the α-carbon-based ENM, and the resulting shifts in soft mode frequencies are used as a metric for evaluating the essentiality of each residue. Results on a dataset of monomeric proteins indicate the enrichment of allosteric and orthosteric binding sites, as well as global hinge regions among essential residues, highlighting the significant role of these sites in controlling the overall structural dynamics. Further integration of ESSA with information on predicted pockets and their local hydrophobicity density enables successful predictions of allosteric pockets for both ligand-bound and -unbound structures. ESSA can be efficiently applied to large multimeric systems. Three case studies, namely (i) G-protein binding to a GPCR, (ii) heterotrimeric assembly of the Ser/Thr protein phosphatase PP2A, and (iii) allo-targeting of AMPA receptor, demonstrate the utility of ESSA for identifying essential sites and narrowing down target allosteric sites identified by druggability simulations.

Pharmmaker: Pharmacophore modeling and hit identification based on druggability simulations

Recent years have seen progress in druggability simulations, that is, molecular dynamics simulations of target proteins in solutions containing drug‐like probe molecules to characterize their drug‐binding abilities, if any. An important consecutive step is to analyze the trajectories to construct pharmacophore models (PMs) to use for virtual screening of libraries of small molecules. While considerable success has been observed in this type of computer‐aided drug discovery, a systematic tool encompassing multiple steps from druggability simulations to pharmacophore modeling, to identifying hits by virtual screening of libraries of compounds, has been lacking. We address this need here by developing a new tool, Pharmmaker, building on the DruGUI module of our ProDy application programming interface. Pharmmaker is composed of a suite of steps: (Step 1) identification of high affinity residues for each probe molecule type; (Step 2) selecting high affinity residues and hot spots in the vicinity of sites identified by DruGUI; (Step 3) ranking of the interactions between high affinity residues and specific probes; (Step 4) obtaining probe binding poses and corresponding protein conformations by collecting top‐ranked snapshots; and (Step 5) using those snapshots for constructing PMs. The PMs are then used as filters for identifying hits in structure‐based virtual screening. Pharmmaker, accessible online at http://prody.csb.pitt.edu/pharmmaker/, can be used in conjunction with other tools available in ProDy.

Structural biology of glutamate receptor ion channels: towards an understanding of mechanism

Ionotropic glutamate receptors (iGluRs) are tetrameric ion channels that mediate signal transmission at neuronal synapses, where they contribute centrally to the postsynaptic plasticity that underlies learning and memory. Receptor activation by l-glutamate triggers complex allosteric cascades that are transmitted through the layered and highly flexible receptor assembly culminating in opening a cation-selective pore. This process is shaped by the arrangement of the four core subunits as well as the presence of various auxiliary subunits, and is subject to regulation by an array of small molecule modulators targeting a number of sites throughout the complex. Here, we discuss recent structures of iGluR homomers and heteromers illuminating the organization and subunit arrangement of the core tetramer, co-assembled with auxiliary subunits and in complex with allosteric modulators.

Molecular Dynamics Simulations of the Allosteric Modulation of the Adenosine A2A Receptor by a Mini-G Protein

Through their coupling to G proteins, G Protein-Coupled Receptors (GPCRs) trigger cellular responses to various signals. Some recent experiments have interestingly demonstrated that the G protein can also act on the receptor by favoring a closed conformation of its orthosteric site, even in the absence of a bound agonist. In this work, we explored such an allosteric modulation by performing extensive molecular dynamics simulations on the adenosine A2 receptor (A2aR) coupled to the Mini-Gs protein. In the presence of the Mini-Gs, we confirmed a restriction of the receptor’s agonist binding site that can be explained by a modulation of the intrinsic network of contacts of the receptor. Of interest, we observed similar effects with the C-terminal helix of the Mini-Gs, showing that the observed effect on the binding pocket results from direct local contacts with the bound protein partner that cause a rewiring of the whole receptor’s interaction network.