scholarly journals Mapping major SARS-CoV-2 drug targets and assessment of druggability using computational fragment screening: Identification of an allosteric small-molecule binding site on the Nsp13 helicase

PLoS ONE ◽  
2021 ◽  
Vol 16 (2) ◽  
pp. e0246181
Author(s):  
Matthew R. Freidel ◽  
Roger S. Armen

The 2019 emergence of, SARS-CoV-2 has tragically taken an immense toll on human life and far reaching impacts on society. There is a need to identify effective antivirals with diverse mechanisms of action in order to accelerate preclinical development. This study focused on five of the most established drug target proteins for direct acting small molecule antivirals: Nsp5 Main Protease, Nsp12 RNA-dependent RNA polymerase, Nsp13 Helicase, Nsp16 2’-O methyltransferase and the S2 subunit of the Spike protein. A workflow of solvent mapping and free energy calculations was used to identify and characterize favorable small-molecule binding sites for an aromatic pharmacophore (benzene). After identifying the most favorable sites, calculated ligand efficiencies were compared utilizing computational fragment screening. The most favorable sites overall were located on Nsp12 and Nsp16, whereas the most favorable sites for Nsp13 and S2 Spike had comparatively lower ligand efficiencies relative to Nsp12 and Nsp16. Utilizing fragment screening on numerous possible sites on Nsp13 helicase, we identified a favorable allosteric site on the N-terminal zinc binding domain (ZBD) that may be amenable to virtual or biophysical fragment screening efforts. Recent structural studies of the Nsp12:Nsp13 replication-transcription complex experimentally corroborates ligand binding at this site, which is revealed to be a functional Nsp8:Nsp13 protein-protein interaction site in the complex. Detailed structural analysis of Nsp13 ZBD conformations show the role of induced-fit flexibility in this ligand binding site and identify which conformational states are associated with efficient ligand binding. We hope that this map of over 200 possible small-molecule binding sites for these drug targets may be of use for ongoing discovery, design, and drug repurposing efforts. This information may be used to prioritize screening efforts or aid in the process of deciphering how a screening hit may bind to a specific target protein.

2019 ◽  
Author(s):  
Michael Olp ◽  
Daniel Sprague ◽  
Stefan Kathman ◽  
Ziyang Xu ◽  
Alexandar Statsyuk ◽  
...  

<p>Brd4, a member of the bromodomain and extraterminal domain (BET) family, has emerged as a promising epigenetic target in cancer and inflammatory disorders. All reported BET family ligands bind within the bromodomain acetyl-lysine binding sites and competitively inhibit BET protein interaction with acetylated chromatin. Alternative chemical probes that act orthogonally to the highly-conserved acetyl-lysine binding sites may exhibit selectivity within the BET family and avoid recently reported toxicity in clinical trials of BET bromodomain inhibitors. Here, we report the first identification of a ligandable site on a bromodomain outside the acetyl-lysine binding site. Inspired by our computational prediction of hotspots adjacent to non-homologous cysteine residues within the <i>C</i>-terminal Brd4 bromodomain (Brd4-BD2), we performed a mid-throughput mass spectrometry screen to identify cysteine-reactive fragments that covalently and selectively modify Brd4. Subsequent mass spectrometry, NMR and computational docking analyses of electrophilic fragment hits revealed a novel ligandable site near Cys356 that is unique to Brd4 among all human bromodomains. This site is orthogonal to the Brd4-BD2 acetyl-lysine binding site as Cys356 modification did not impact binding of the pan-BET bromodomain inhibitor JQ1 in fluorescence polarization assays. Finally, we tethered covalent fragments to JQ1 and performed NanoBRET assays to provide proof of principle that this orthogonal site can be covalently targeted in intact human cells. Overall, we demonstrate the potential of targeting sites orthogonal to bromodomain acetyl-lysine binding sites to develop bivalent and covalent inhibitors that displace Brd4 from chromatin.</p>


2015 ◽  
Vol 112 (52) ◽  
pp. 15910-15915 ◽  
Author(s):  
R. Frederick Ludlow ◽  
Marcel L. Verdonk ◽  
Harpreet K. Saini ◽  
Ian J. Tickle ◽  
Harren Jhoti

Proteins need to be tightly regulated as they control biological processes in most normal cellular functions. The precise mechanisms of regulation are rarely completely understood but can involve binding of endogenous ligands and/or partner proteins at specific locations on a protein that can modulate function. Often, these additional secondary binding sites appear separate to the primary binding site, which, for example for an enzyme, may bind a substrate. In previous work, we have uncovered several examples in which secondary binding sites were discovered on proteins using fragment screening approaches. In each case, we were able to establish that the newly identified secondary binding site was biologically relevant as it was able to modulate function by the binding of a small molecule. In this study, we investigate how often secondary binding sites are located on proteins by analyzing 24 protein targets for which we have performed a fragment screen using X-ray crystallography. Our analysis shows that, surprisingly, the majority of proteins contain secondary binding sites based on their ability to bind fragments. Furthermore, sequence analysis of these previously unknown sites indicate high conservation, which suggests that they may have a biological function, perhaps via an allosteric mechanism. Comparing the physicochemical properties of the secondary sites with known primary ligand binding sites also shows broad similarities indicating that many of the secondary sites may be druggable in nature with small molecules that could provide new opportunities to modulate potential therapeutic targets.


1992 ◽  
Vol 285 (2) ◽  
pp. 419-425 ◽  
Author(s):  
U Christensen ◽  
L Mølgaard

The kinetics of a series of Glu-plasminogen ligand-binding processes were investigated at pH 7.8 and 25 degrees C (in 0.1 M-NaCl). The ligands include compounds analogous to C-terminal lysine residues and to normal lysine residues. Changes of the Glu-plasminogen protein fluorescence were measured in a stopped-flow instrument as a function of time after rapid mixing of Glu-plasminogen and ligand at various concentrations. Large positive fluorescence changes (approximately 10%) accompany the ligand-induced conformational changes of Glu-plasminogen resulting from binding at weak lysine-binding sites. Detailed studies of the concentration-dependencies of the equilibrium signals and the rate constants of the process induced by various ligands showed the conformational change to involve two sites in a concerted positive co-operative process with three steps: (i) binding of a ligand at a very weak lysine-binding site that preferentially, but not exclusively, binds C-terminal-type lysine ligands, (ii) the rate-determining actual-conformational-change step and (iii) binding of one more lysine ligand at a second weak lysine-binding site that then binds the ligand more tightly. Further, totally independent initial small negative fluorescence changes (approximately 2-4%) corresponding to binding at the strong lysine-binding site of kringle 1 [Sottrup-Jensen, Claeys, Zajdel, Petersen & Magnusson (1978) Prog. Chem. Fibrinolysis Thrombolysis 3, 191-209] were observed for the C-terminal-type ligands. The finding that the conformational change in Glu-plasminogen involves two weak lysine-binding sites indicates that the effect cannot be assigned to any single kringle and that the problem of whether kringle 4 or kringle 5 is responsible for the process resolves itself. Probably kringle 4 and 5 are both participating. The involvement of two lysine binding-sites further makes the high specificity of Glu-plasminogen effectors more conceivable.


2008 ◽  
Vol 412 (1) ◽  
pp. 103-112 ◽  
Author(s):  
Doreen Thor ◽  
Angela Schulz ◽  
Thomas Hermsdorf ◽  
Torsten Schöneberg

GPCRs (G-protein-coupled receptors) exist in a spontaneous equilibrium between active and inactive conformations that are stabilized by agonists and inverse agonists respectively. Because ligand binding of agonists and inverse agonists often occurs in a competitive manner, one can assume an overlap between both binding sites. Only a few studies report mutations in GPCRs that convert receptor blockers into agonists by unknown mechanisms. Taking advantage of a genetically modified yeast strain, we screened libraries of mutant M3Rs {M3 mAChRs [muscarinic ACh (acetylcholine) receptors)]} and identified 13 mutants which could be activated by atropine (EC50 0.3–10 μM), an inverse agonist on wild-type M3R. Many of the mutations sensitizing M3R to atropine activation were located at the junction of intracellular loop 3 and helix 6, a region known to be involved in G-protein coupling. In addition to atropine, the pharmacological switch was found for other M3R blockers such as scopolamine, pirenzepine and oxybutynine. However, atropine functions as an agonist on the mutant M3R only when expressed in yeast, but not in mammalian COS-7 cells, although high-affinity ligand binding was comparable in both expression systems. Interestingly, we found that atropine still blocks carbachol-induced activation of the M3R mutants in the yeast expression system by binding at the high-affinity-binding site (Ki ∼10 nM). Our results indicate that blocker-to-agonist converting mutations enable atropine to function as both agonist and antagonist by interaction with two functionally distinct binding sites.


Biochemistry ◽  
2017 ◽  
Vol 56 (30) ◽  
pp. 3972-3982 ◽  
Author(s):  
Eline Sijbesma ◽  
Lukasz Skora ◽  
Seppe Leysen ◽  
Luc Brunsveld ◽  
Uwe Koch ◽  
...  

2022 ◽  
Author(s):  
Adam Zemla ◽  
Jonathan E. Allen ◽  
Dan Kirshner ◽  
Felice C. Lightstone

We present a structure-based method for finding and evaluating structural similarities in protein regions relevant to ligand binding. PDBspheres comprises an exhaustive library of protein structure regions (spheres) adjacent to complexed ligands derived from the Protein Data Bank (PDB), along with methods to find and evaluate structural matches between a protein of interest and spheres in the library. Currently, PDBspheres library contains more than 2 million spheres, organized to facilitate searches by sequence and/or structure similarity of protein-ligand binding sites or interfaces between interacting molecules. PDBspheres uses the LGA structure alignment algorithm as the main engine for detecting structure similarities between the protein of interest and library spheres. An all-atom structure similarity metric ensures that sidechain placement is taken into account in the PDBspheres primary assessment of confidence in structural matches. In this paper, we (1) describe the PDBspheres method, (2) demonstrate how PDBspheres can be used to detect and characterize binding sites in protein structures, (3) compare PDBspheres use for binding site prediction with seven other binding site prediction methods using a curated dataset of 2,528 ligand-bound and ligand-free crystal structures, and (4) use PDBspheres to cluster pockets and assess structural similarities among protein binding sites of the 4,876 structures in the refined set of PDBbind 2019 dataset. The PDBspheres library is made publicly available for download at https://proteinmodel.org/AS2TS/PDBspheres


Author(s):  
David Slochower ◽  
Niel Henriksen ◽  
Lee-Ping Wang ◽  
John Chodera ◽  
David Mobley ◽  
...  

<div><div><div><p>Designing ligands that bind their target biomolecules with high affinity and specificity is a key step in small- molecule drug discovery, but accurately predicting protein-ligand binding free energies remains challenging. Key sources of errors in the calculations include inadequate sampling of conformational space, ambiguous protonation states, and errors in force fields. Noncovalent complexes between a host molecule with a binding cavity and a drug-like guest molecules have emerged as powerful model systems. As model systems, host-guest complexes reduce many of the errors in more complex protein-ligand binding systems, as their small size greatly facilitates conformational sampling, and one can choose systems that avoid ambiguities in protonation states. These features, combined with their ease of experimental characterization, make host-guest systems ideal model systems to test and ultimately optimize force fields in the context of binding thermodynamics calculations.</p><p><br></p><p>The Open Force Field Initiative aims to create a modern, open software infrastructure for automatically generating and assessing force fields using data sets. The first force field to arise out of this effort, named SMIRNOFF99Frosst, has approximately one tenth the number of parameters, in version 1.0.5, compared to typical general small molecule force fields, such as GAFF. Here, we evaluate the accuracy of this initial force field, using free energy calculations of 43 α and β-cyclodextrin host-guest pairs for which experimental thermodynamic data are available, and compare with matched calculations using two versions of GAFF. For all three force fields, we used TIP3P water and AM1-BCC charges. The calculations are performed using the attach-pull-release (APR) method as implemented in the open source package, pAPRika. For binding free energies, the root mean square error of the SMIRNOFF99Frosst calculations relative to experiment is 0.9 [0.7, 1.1] kcal/mol, while the corresponding results for GAFF 1.7 and GAFF 2.1 are 0.9 [0.7, 1.1] kcal/mol and 1.7 [1.5, 1.9] kcal/mol, respectively, with 95% confidence ranges in brackets. These results suggest that SMIRNOFF99Frosst performs competitively with existing small molecule force fields and is a parsimonious starting point for optimization.</p></div></div></div>


2021 ◽  
Author(s):  
Rishal Aggarwal ◽  
Akash Gupta ◽  
Vineeth Chelur ◽  
C. V. Jawahar ◽  
U. Deva Priyakumar

<div> A structure-based drug design pipeline involves the development of potential drug molecules or ligands that form stable complexes with a given receptor at its binding site. A prerequisite to this is finding druggable and functionally relevant binding sites on the 3D structure of the protein. Although several methods for detecting binding sites have been developed beforehand, a majority of them surprisingly fail in the identification and ranking of binding sites accurately. The rapid adoption and success of deep learning algorithms in various sections of structural biology beckons the usage of such algorithms for accurate binding site detection. As a combination of geometry based software and deep learning, we report a novel framework, DeepPocket that utilises 3D convolutional neural networks for the rescoring of pockets identified by Fpocket and further segments these identified cavities on the protein surface. Apart from this, we also propose another dataset SC6K containing protein structures submitted in the Protein Data Bank (PDB) from January 2018 till February 2020 for ligand binding site (LBS) detection. DeepPocket's results on various binding site datasets and SC6K highlights its better performance over current state-of-the-art methods and good generalization ability over novel structures. </div><div><br></div>


2019 ◽  
Vol 47 (W1) ◽  
pp. W345-W349 ◽  
Author(s):  
Lukas Jendele ◽  
Radoslav Krivak ◽  
Petr Skoda ◽  
Marian Novotny ◽  
David Hoksza

AbstractPrankWeb is an online resource providing an interface to P2Rank, a state-of-the-art method for ligand binding site prediction. P2Rank is a template-free machine learning method based on the prediction of local chemical neighborhood ligandability centered on points placed on a solvent-accessible protein surface. Points with a high ligandability score are then clustered to form the resulting ligand binding sites. In addition, PrankWeb provides a web interface enabling users to easily carry out the prediction and visually inspect the predicted binding sites via an integrated sequence-structure view. Moreover, PrankWeb can determine sequence conservation for the input molecule and use this in both the prediction and result visualization steps. Alongside its online visualization options, PrankWeb also offers the possibility of exporting the results as a PyMOL script for offline visualization. The web frontend communicates with the server side via a REST API. In high-throughput scenarios, therefore, users can utilize the server API directly, bypassing the need for a web-based frontend or installation of the P2Rank application. PrankWeb is available at http://prankweb.cz/, while the web application source code and the P2Rank method can be accessed at https://github.com/jendelel/PrankWebApp and https://github.com/rdk/p2rank, respectively.


Sign in / Sign up

Export Citation Format

Share Document