Allosteric Molecules Research Articles

Potential allosteric pockets identification of glucagon receptor based on molecular dynamics simulations

Glucagon receptor (GCGR) is an important target for the treatment of type 2 diabetes mellitus. Although several small molecules with antagonistic activity have been discovered, so far, only one small molecule binding site has been resolved. To discover more novel allosteric pockets and allosteric molecules, we started with the unique full-length inactive conformation of GCGR and applied all-atom molecular dynamics (MD) simulations to obtain extensive dynamic conformations of the GCGR/glucagon complex. For the first time, MDpocket, FTMove and FTMap were used to detect allosteric pockets in simulation trajectories, selecting 4 stable pockets with a total of 14 structures as templates for virtual screening. From the results of virtual screening, 14 compounds were ultimately selected after a series of filtering steps. The cAMP accumulation assay indicated that compound gs6 has antagonistic activity, and MD simulations further revealed the allosteric mechanism of gs6. We are the first to identify new allosteric pockets and allosteric molecules in simulation trajectories of the GCGR/glucagon complex, providing a reference for research on other G-protein-coupled receptors (GPCR). However, there is still considerable room for improvement, such as using more simulation methods to obtain a richer set of dynamic conformations.

Aptamer-Based DNA Allosteric Switch for Regulation of Protein Activity.

Allostery is a fundamental way to regulate the function of biomolecules playing crucial roles in cell metabolism and proliferation and is deemed the second secret of life. Given the limited understanding of the structure of natural allosteric molecules, the development of artificial allosteric molecules brings a huge opportunity to transform the allosteric mechanism into practical applications. In this study, the concept of bionics is introduced into the design of artificial allosteric molecules and an allosteric DNA switch with an activity site and an allosteric site based on two aptamers for selective inhibition of thrombin activity. Compared with the single aptamer, the allosteric switch possesses a significantly enhanced inhibition ability, which can be precisely regulated by converting the switch states. Moreover, the dynamic allosteric switch is further subjected to the control of the DNA threshold circuit for realizing automatic concentration determination and activity inhibition of thrombin. These compelling results confirm that this allosteric switch equipped with self-sensing and information-processing modules puts a new slant on the research of allosteric mechanisms and further application of allosteric tactics in chemical and biomedical fields.

Exploring Binding Pockets in the Conformational States of the SARS-CoV-2 Spike Trimers for the Screening of Allosteric Inhibitors Using Molecular Simulations and Ensemble-Based Ligand Docking.

Understanding mechanisms of allosteric regulation remains elusive for the SARS-CoV-2 spike protein, despite the increasing interest and effort in discovering allosteric inhibitors of the viral activity and interactions with the host receptor ACE2. The challenges of discovering allosteric modulators of the SARS-CoV-2 spike proteins are associated with the diversity of cryptic allosteric sites and complex molecular mechanisms that can be employed by allosteric ligands, including the alteration of the conformational equilibrium of spike protein and preferential stabilization of specific functional states. In the current study, we combine conformational dynamics analysis of distinct forms of the full-length spike protein trimers and machine-learning-based binding pocket detection with the ensemble-based ligand docking and binding free energy analysis to characterize the potential allosteric binding sites and determine structural and energetic determinants of allosteric inhibition for a series of experimentally validated allosteric molecules. The results demonstrate a good agreement between computational and experimental binding affinities, providing support to the predicted binding modes and suggesting key interactions formed by the allosteric ligands to elicit the experimentally observed inhibition. We establish structural and energetic determinants of allosteric binding for the experimentally known allosteric molecules, indicating a potential mechanism of allosteric modulation by targeting the hinges of the inter-protomer movements and blocking conformational changes between the closed and open spike trimer forms. The results of this study demonstrate that combining ensemble-based ligand docking with conformational states of spike protein and rigorous binding energy analysis enables robust characterization of the ligand binding modes, the identification of allosteric binding hotspots, and the prediction of binding affinities for validated allosteric modulators, which is consistent with the experimental data. This study suggested that the conformational adaptability of the protein allosteric sites and the diversity of ligand bound conformations are both in play to enable efficient targeting of allosteric binding sites and interfere with the conformational changes.

Allostery: The Good, the Bad, and the Ugly.

With the advent of functional screening, more allosteric molecules are being discovered and developed as possible therapeutic entities. Allosteric proteins are unique because of two specific properties: 1) separate binding sites for allosteric modulators and guests and 2) mandatory alteration of receptor conformation upon binding of allosteric modulators. For G protein-coupled receptors, these properties produce many beneficial effects on pharmacologic systems that are described here. Allosteric discovery campaigns also bring with them added considerations that must be addressed for the endeavor to be successful, and these are described herein as well. SIGNIFICANCE STATEMENT: Recent years have seen the increasing presence of allosteric molecules as possible therapeutic drug candidates. The scientific procedures to characterize these are unique and require special techniques, so it is imperative that scientists understand the new concepts involved in allosteric function. This review examines the reasons why allosteric molecules should be considered as new drug entities and the techniques required to optimize the discovery process for allosteric molecules.

The Role of Sex Hormones in Degenerative Disc Disease.

Narrative review. The purpose of this review is to outline the role of sex hormones, particularly estrogen, in the pathogenesis of degenerative disc disease (DDD). A narrative review of studies discussing sex hormones and intervertebral disc (IVD) degeneration was conducted through a search of bibliographic databases to identify various mechanisms involved in effectuating DDD. Estrogen-deficient states negatively impact various aspects of IVD function. These internal hormone environments reflect routine changes that commonly arise with physiologic aging and can compromise IVD structural integrity through a host of processes. Additionally, allosteric molecules such as micro-RNAs (mi-RNAs) and G protein-coupled estrogen receptors (GPER) antagonists can bind to estrogen receptors and inhibit protective downstream effects with estrogen receptor signaling. Furthermore, cursory studies have observed chondrogenic effects with testosterone supplementation, although the specific mechanism remains unclear. Regulation of sex hormones, namely estrogen and testosterone, significantly impacts the structural integrity and function of IVDs. Uncovering underlying interactions driving these regulatory processes can facilitate development of novel, clinical therapies to treat DDD.

Unmasking allosteric binding sites: Novel targets for GPCR drug discovery

ABSTRACT Introduction Unexpected non-apparent and hidden allosteric binding sites are non-classical and non-apparent allosteric centers in 3-D X-ray protein structures until orthosteric or allosteric ligands bind to them. The orthosteric center of one protomer that modulates binding centers of the other protomers within an oligomer is also an unexpected allosteric site. Furthermore, another partner protein can also produce these effects, acting as an unexpected allosteric modulator. Areas covered This review summarizes both classical and non-classical allosterism. The authors focus on G protein-coupled receptor (GPCR) oligomers as a paradigm of allosteric molecules. Moreover, they show several examples of unexpected allosteric sites such as hidden allosteric sites in a protomer that appear after the interaction with other molecules and the allosterism exerted between orthosteric sites within GPCR oligomer, emphasizing on the allosteric modulations that can occur between binding sites. Expert opinion The study of these new non-classical allosteric sites will expand the diversity of allosteric control on the function of orthosteric sites within proteins, whether GPCRs or other receptors, enzymes or transporters. Moreover, the design of new drugs targeting these hidden allosteric sites or already known orthosteric sites acting as allosteric sites in protein homo- or hetero-oligomers will increase the therapeutic potential of allosterism.

Machine Learning of Allosteric Effects: The Analysis of Ligand-Induced Dynamics to Predict Functional Effects in TRAP1.

Allosteric molecules provide a powerful means to modulate protein function. However, the effect of such ligands on distal orthosteric sites cannot be easily described by classical docking methods. Here, we applied machine learning (ML) approaches to expose the links between local dynamic patterns and different degrees of allosteric inhibition of the ATPase function in the molecular chaperone TRAP1. We focused on 11 novel allosteric modulators with similar affinities to the target but with inhibitory efficacy between the 26.3 and 76%. Using a set of experimentally related local descriptors, ML enabled us to connect the molecular dynamics (MD) accessible to ligand-bound (perturbed) and unbound (unperturbed) systems to the degree of ATPase allosteric inhibition. The ML analysis of the comparative perturbed ensembles revealed a redistribution of dynamic states in the inhibitor-bound versus inhibitor-free systems following allosteric binding. Linear regression models were built to quantify the percentage of experimental variance explained by the predicted inhibitor-bound TRAP1 states. Our strategy provides a comparative MD–ML framework to infer allosteric ligand functionality. Alleviating the time scale issues which prevent the routine use of MD, a combination of MD and ML represents a promising strategy to support in silico mechanistic studies and drug design.

Transition in relaxation paths in allosteric molecules: Enzymatic kinetically constrained model

A hierarchy of timescales is ubiquitous in biological systems, where enzymatic reactions play an important role because they can hasten the relaxation to equilibrium. We introduced a statistical physics model of interacting spins that also incorporates enzymatic reactions to extend the classic model for allosteric regulation. Through Monte Carlo simulations, we found that the relaxation dynamics are much slower than the elementary reactions and are logarithmic in time with several plateaus, as is commonly observed for glasses. This is because of the kinetic constraints from the cooperativity via the competition for an enzyme, which has different affinity for molecules with different structures. Our model showed symmetry breaking in the relaxation trajectories that led to inherently kinetic transitions without any correspondence to the equilibrium state. In this paper, we discuss the relevance of these results for diverse responses in biology.

Allostery and Kinetic Proofreading.

Kinetic proofreading is an error correction mechanism present in the processes of the central dogma and beyond and typically requires the free energy of nucleotide hydrolysis for its operation. Though the molecular players of many biological proofreading schemes are known, our understanding of how energy consumption is managed to promote fidelity remains incomplete. In our work, we introduce an alternative conceptual scheme called "the piston model of proofreading" in which enzyme activation through hydrolysis is replaced with allosteric activation achieved through mechanical work performed by a piston on regulatory ligands. Inspired by Feynman's ratchet and pawl mechanism, we consider a mechanical engine designed to drive the piston actions powered by a lowering weight, whose function is analogous to that of ATP synthase in cells. Thanks to its mechanical design, the piston model allows us to tune the "knobs" of the driving engine and probe the graded changes and trade-offs between speed, fidelity, and energy dissipation. It provides an intuitive explanation of the conditions necessary for optimal proofreading and reveals the unexpected capability of allosteric molecules to beat the Hopfield limit of fidelity by leveraging the diversity of states available to them. The framework that we have built for the piston model can also serve as a basis for additional studies of driven biochemical systems.

Selective Allosteric Modulation of N-Terminally Cleaved, but Not Full Length CCL3 in CCR1

Chemokines undergo post-translational modification such as N-terminal truncations. Here, we describe how N-terminal truncation of full length CCL3(1-70) affects its activity at CCR1. Truncated CCL3(5-70) has 10-fold higher potency and enhanced efficacy in β-arrestin recruitment, but less than 2-fold increased potencies in G protein signaling determined by calcium release, cAMP and IP3 formation. Small positive ago-allosteric ligands modulate the two CCL3 variants differently as the metal ion chelator bipyridine in complex with zinc (ZnBip) enhances the binding of truncated, but not full length CCL3, while a size-increase of the chelator to a chloro-substituted terpyridine (ZnClTerp), eliminates its allosteric, but not agonistic action. By employing a series of receptor mutants and in silico modeling we describe residues of importance for chemokine and small molecule binding. Notably, the chemokine receptor-conserved Glu2877.39 interacts with the N-terminal amine of truncated CCL3(5-70) and with Zn2+ of ZnBip, thereby bridging their binding sites and enabling the positive allosteric effect. Our study emphasizes that small allosteric molecules may act differently toward chemokine variants and thus selectively modulate interactions of specific chemokine subsets with their cognate receptors.

Single-Molecule Mechanistic Study of Enzyme Hysteresis.

Hysteresis is an important feature of enzyme-catalyzed reactions, as it reflects the influence of enzyme regulation in the presence of ligands such as substrates or allosteric molecules. In typical kinetic studies of enzyme activity, hysteretic behavior is observed as a “lag” or “burst” in the time course of the catalyzed reaction. These lags and bursts are due to the relatively slow transition from one state to another state of the enzyme molecule, with different states having different kinetic properties. However, it is difficult to understand the underlying mechanism of hysteresis by observing bulk reactions because the different enzyme molecules in the population behave stochastically. In this work, we studied the hysteretic behavior of mutant β-glucuronidase (GUS) using a high-throughput single-molecule array platform and investigated the effect of thermal treatment on the hysteresis.

The pharmacologic characterization of allosteric molecules: Gq protein activation

Context: Drugs such as positive allosteric modulators (PAMs) produce complex behaviors when acting on tissues in different physiological contexts in vivo.Objective: This study describes the use of functional assays of varying receptor sensitivity to unveil the various behaviors of PAMs and thus quantify allosteric effect through system independent scales.Materials and methods: Muscarinic receptor activation with acetylcholine (ACh) was used to the demonstrate activity of the PAM agonist 1–(4-methoxybenzyl)-4-oxo-1,4-dihydroquinoline-3-carboxylic acid, Benzyl quinolone carboxylic acid (BQCA) in terms of direct agonism, potentiation of ACh affinity, and ACh efficacy. Concentration–response curves were fit to the functional allosteric model to yield indices of agonism (τB), effects on affinity (α cooperativity), and efficacy (β cooperativity).Results: It is shown that a highly sensitive functional assay revealed the direct efficacy of BQCA as an agonist and relatively insensitive cells (produced by chemical alkylation of muscarinic receptor with phenoxybenzamine) revealed a positive allosteric effect of BQCA on ACh efficacy. A wide range of functional assay sensitivities produced a complex pattern of behavior for BQCA all of which was accurately quantified through the system-independent parameters of the functional allosteric model.Conclusions: The study of complex allosteric molecules in a range of functional assays of varying sensitivity allows the measurement of the complete array of activities of these molecules on receptors and also better predicts which will be seen with these in vivo where a range of tissue sensitivities is encountered.

Combinatorial Control through Allostery.

Many instances of cellular signaling and transcriptional regulation involve switch-like molecular responses to the presence or absence of input ligands. To understand how these responses come about and how they can be harnessed, we develop a statistical mechanical model to characterize the types of Boolean logic that can arise from allosteric molecules following the Monod-Wyman-Changeux (MWC) model. Building upon previous work, we show how an allosteric molecule regulated by two inputs can elicit AND, OR, NAND, and NOR responses but is unable to realize XOR or XNOR gates. Next, we demonstrate the ability of an MWC molecule to perform ratiometric sensing-a response behavior where activity depends monotonically on the ratio of ligand concentrations. We then extend our analysis to more general schemes of combinatorial control involving either additional binding sites for the two ligands or an additional third ligand and show how these additions can cause a switch in the logic behavior of the molecule. Overall, our results demonstrate the wide variety of control schemes that biological systems can implement using simple mechanisms.

Progress in Allosteric Database.

An allosteric mechanism refers to the biological regulation process wherein macromolecules propagate the effect of ligand binding at one site to a spatially distant orthosteric locus, thus affecting activity. The theory has remained a trending topic in biology research for over 50years, since the understanding of allostery is fundamental for gleaning numerous biological processes and developing new drug therapies. In the past two decades, the allosteric paradigm has evolved into more descriptive models, with ever-expanding amounts of experimental data pertaining to newly identified allosteric molecules. The AlloSteric Database (ASD, accessible at http://mdl.shsmu.edu.cn/ASD ), which is a comprehensive knowledge repository, has provided the public with integrated information encompassing allosteric proteins, modulators, sites, pathways, and networks to investigate allostery since 2009. In this chapter, we introduce the history and usage of the ASD and give attention to specific applications that have benefited from the ASD.

Discovery of a novel conformational equilibrium in urokinase-type plasminogen activator

Although trypsin-like serine proteases have flexible surface-exposed loops and are known to adopt higher and lower activity conformations, structural determinants for the different conformations have remained largely obscure. The trypsin-like serine protease, urokinase-type plasminogen activator (uPA), is central in tissue remodeling processes and also strongly implicated in tumor metastasis. We solved five X-ray crystal structures of murine uPA (muPA) in the absence and presence of allosteric molecules and/or substrate-like molecules. The structure of unbound muPA revealed an unsuspected non-chymotrypsin-like protease conformation in which two β-strands in the core of the protease domain undergoes a major antiparallel-to-parallel conformational transition. We next isolated two anti-muPA nanobodies; an active-site binding nanobody and an allosteric nanobody. Crystal structures of the muPA:nanobody complexes and hydrogen-deuterium exchange mass spectrometry revealed molecular insights about molecular factors controlling the antiparallel-to-parallel equilibrium in muPA. Together with muPA activity assays, the data provide valuable insights into regulatory mechanisms and conformational flexibility of uPA and trypsin-like serine proteases in general.

Mechanistic Models Fit to Variable Temperature Calorimetric Data Provide Insights into Cooperativity

Allostery pervades macromolecular function and drives cooperative binding of ligands to macromolecules. To decipher the mechanisms of cooperative ligand binding it is necessary to define at a microscopic level the structural and thermodynamic consequences of binding of each ligand to its allosterically coupled site(s). However, dynamic sampling of alternative conformations (microstates) in allosteric molecules complicates interpretation of both structural and thermodynamic data. Isothermal titration calorimetry has the potential to directly quantify the thermodynamics of allosteric interactions, but usually falls short of enabling mechanistic insight. This is because 1) its measurements reflect the sum of overlapping caloric processes involving binding-linked population shifts within and between microstates, and 2) data are generally fit with phenomenological binding polynomials that are underdetermined. Nevertheless, temperature-dependent binding data have the potential to resolve overlapping thermodynamic processes, while mechanistically constrained models enable hypothesis testing and identification of informative parameters. We globally fit temperature-dependent isothermal titration calorimetry data for binding of 11 tryptophan ligands to the homo-undecameric trp RNA-binding Attenuation Protein from Bacillus stearothermophilus using nearest-neighbor statistical thermodynamic models. This approach allowed us to distinguish alternative nearest-neighbor interaction models, and quantifies the thermodynamic contribution of neighboring ligands to individual binding sites. We also perform conventional Hill equation modeling and illustrate how comparatively limited it is in quantitative or mechanistic value. This work illustrates the potential of mechanistically constrained global fitting of binding data to yield the microscopic thermodynamic parameters essential for deciphering mechanisms of cooperativity in a wide range of ligand-regulated homo-oligomeric assemblies.

Effect of allosteric molecules on structure and drug affinity of HIV-1 protease by molecular dynamics simulations

Recent experiments show that small molecules can bind onto the allosteric sites of HIV-1 protease (PR), which provides a starting point for developing allosteric inhibitors. However, the knowledge of the effect of such binding on the structural dynamics and binding free energy of the active site inhibitor and PR is still lacking. Here, we report 200ns long molecular dynamics simulation results to gain insight into the influences of two allosteric molecules (1H-indole-6-carboxylic acid, 1F1 and 2-methylcyclohexano, 4D9). The simulations demonstrate that both allosteric molecules change the PR conformation and stabilize the structures of PR and the inhibitor; the residues of the flaps are sensitive to the allosteric molecules and the flexibility of the residues is pronouncedly suppressed; the additions of the small molecules to the allosteric sites strengthen the binding affinities of 3TL-PR by about 12–15kal/mol in the binding free energy, which mainly arises from electrostatic term. Interestingly, it is found that the action mechanisms of 1F1 and 4D9 are different, the former behaviors like a doorman that keeps the inhibitor from escape and makes the flaps (door) partially open; the latter is like a wedge that expands the allosteric space and meanwhile closes the flaps. Our data provide a theoretical support for designing the allosteric inhibitor.

On coupled oscillator dynamics and incident behaviour patterns in slime mould Physarum polycephalum: emergence of wave packets, global streaming clock frequencies and anticipation of periodic stimuli

Slime mould Physarum polycephalum is a single cell which physically oscillates via contraction of actomyosin in order to achieve motility. Several of its apparently ‘intelligent’ behaviour patterns such as anticipatory responses to periodic stimuli have recently been attributed as functions of the coupling between the oscillating intracellular reactions which drive its rhythmic muscular contraction, but the mechanisms that underlie these phenomena have not yet been experimentally verified. Through laboratory investigations in which we entrain the P. polycephalum plasmodium via periodic ultraviolet light exposure we find that this phenomenon is likely to result from biasing its various oscillating life processes through altering local concentration profiles of various allosteric molecules and their effectors. This temporarily overwrites the global streaming clock frequency and eradicates the wave packets usually observed in slime mould biomechanical oscillation. This response is likened to an intracellular chemical memory. We proceed to present a multi-agent model in which we demonstrate that travelling waves and oscillatory clock frequencies may emerge in the virtual organism’s biomechanical oscillator, although anticipatory responses cannot be replicated by simple mechanical interactions. We conclude by arguing that these phenomena are best characterised as analogue computation and discuss practical applications therein.

Structural Studies of the HIV-1 Integrase Protein: Compound Screening and Characterization of a DNA-Binding Inhibitor.

Understanding the HIV integrase protein and mechanisms of resistance to HIV integrase inhibitors is complicated by the lack of a full length HIV integrase crystal structure. Moreover, a lentiviral integrase structure with co-crystallised DNA has not been described. For these reasons, we have developed a structural method that utilizes free software to create quaternary HIV integrase homology models, based partially on available full-length prototype foamy virus integrase structures as well as several structures of truncated HIV integrase. We have tested the utility of these models in screening of small anti-integrase compounds using randomly selected molecules from the ZINC database as well as a well characterized IN:DNA binding inhibitor, FZ41, and a putative IN:DNA binding inhibitor, HDS1. Docking studies showed that the ZINC compounds that had the best binding energies bound at the IN:IN dimer interface and that the FZ41 and HDS1 compounds docked at approximately the same location in integrase, i.e. behind the DNA binding domain, although there is some overlap with the IN:IN dimer interface to which the ZINC compounds bind. Thus, we have revealed two possible locations in integrase that could potentially be targeted by allosteric integrase inhibitors, that are distinct from the binding sites of other allosteric molecules such as LEDGF inhibitors. Virological and biochemical studies confirmed that HDS1 and FZ41 share a similar activity profile and that both can inhibit each of integrase and reverse transcriptase activities. The inhibitory mechanism of HDS1 for HIV integrase seems to be at the DNA binding step and not at either of the strand transfer or 3' processing steps of the integrase reaction. Furthermore, HDS1 does not directly interact with DNA. The modeling and docking methodology described here will be useful for future screening of integrase inhibitors as well as for the generation of models for the study of integrase drug resistance.

The Role of a Sodium Ion Binding Site in the Allosteric Modulation of the A2A Adenosine G Protein-Coupled Receptor

The function of G protein-coupled receptors (GPCRs) can be modulated by a number of endogenous allosteric molecules. In this study, we used molecular dynamics, radioligand binding, and thermostability experiments to elucidate the role of the recently discovered sodium ion binding site in the allosteric modulation of the human A(2A) adenosine receptor, conserved among class A GPCRs. While the binding of antagonists and sodium ions to the receptor was noncompetitive in nature, the binding of agonists and sodium ions appears to require mutually exclusive conformational states of the receptor. Amiloride analogs can also bind to the sodium binding pocket, showing distinct patterns of agonist and antagonist modulation. These findings suggest that physiological concentrations of sodium ions affect functionally relevant conformational states of GPCRs and can help to design novel synthetic allosteric modulators or bitopic ligands exploiting the sodium ion binding pocket.