scholarly journals Hybrid Dynamic Pharmacophore Models as Effective Tools to Identify Novel Chemotypes for Anti-TB Inhibitor Design: A Case Study With Mtb-DapB

2020 ◽  
Vol 8 ◽  
Author(s):  
Chinmayee Choudhury ◽  
Anshu Bhardwaj

Antimicrobial resistance (AMR) is one of the most serious global public health threats as it compromises the successful treatment of deadly infectious diseases like tuberculosis. New therapeutics are constantly needed but it takes a long time and is expensive to explore new biochemical space. One way to address this issue is to repurpose the validated targets and identify novel chemotypes that can simultaneously bind to multiple binding pockets of these targets as a new lead generation strategy. This study reports such a strategy, dynamic hybrid pharmacophore model (DHPM), which represents the combined interaction features of different binding pockets contrary to the conventional approaches, where pharmacophore models are generated from single binding sites. We have considered Mtb-DapB, a validated mycobacterial drug target, as our model system to explore the effectiveness of DHPMs to screen novel unexplored compounds. Mtb-DapB has a cofactor binding site (CBS) and an adjacent substrate binding site (SBS). Four different model systems of Mtb-DapB were designed where, either NADPH/NADH occupies CBS in presence/absence of an inhibitor 2, 6-PDC in the adjacent SBS. Two more model systems were designed, where 2, 6-PDC was linked to NADPH and NADH to form hybrid molecules. The six model systems were subjected to 200 ns molecular dynamics simulations and trajectories were analyzed to identify stable ligand-receptor interaction features. Based on these interactions, conventional pharmacophore models (CPM) were generated from the individual binding sites while DHPMs were created from hybrid-molecules occupying both binding sites. A huge library of 1,563,764 publicly available molecules were screened by CPMs and DHPMs. The screened hits obtained from both types of models were compared based on their Hashed binary molecular fingerprints and 4-point pharmacophore fingerprints using Tanimoto, Cosine, Dice and Tversky similarity matrices. Molecules screened by DHPM exhibited significant structural diversity, better binding strength and drug like properties as compared to the compounds screened by CPMs indicating the efficiency of DHPM to explore new chemical space for anti-TB drug discovery. The idea of DHPM can be applied for a wide range of mycobacterial or other pathogen targets to venture into unexplored chemical space.

2020 ◽  
Author(s):  
Chinmayee Choudhury ◽  
Anshu Bhardwaj

AbstractAntimicrobial resistance (AMR) is one of the most serious global public health threats as it compromises the successful treatment of deadly infectious diseases like tuberculosis. New therapeutics are constantly needed but it takes a long time and is expensive to explore new biochemical space. One way to address this issue is to repurpose the validated targets and identify novel chemotypes that can simultaneously bind to multiple binding pockets of these targets as a new lead generation strategy. This study reports such a strategy, dynamic hybrid pharmacophore model (DHPM), which represents the combined interaction features of different binding pockets contrary to the conventional approaches, where pharmacophore models are generated from single binding sites. We have considered Mtb-DapB, a validated mycobacterial drug target, as our model system to explore the effectiveness of DHPMs to screen novel unexplored compounds. Mtb-DapB has a cofactor binding site (CBS) and an adjacent substrate binding site (SBS). Four different model systems of Mtb-DapB were designed where, either NADPH/NADH occupies CBS in presence/absence of an inhibitor 2, 6-PDC in the adjacent SBS. Two more model systems were designed, where 2, 6-PDC was linked to NADPH and NADH to form hybrid molecules. The six model systems were subjected to 200ns molecular dynamics simulations and trajectories were analyzed to identify stable ligand-receptor interaction features. Based on these interactions, conventional pharmacophore models (CPM) were generated from the individual binding sites while DHPMs were created from hybrid-molecules occupying both binding sites. A huge library of 15, 63,764 publicly available molecules were screened by CPMs and DHPMs. The screened hits obtained from both types of models were compared based on their Hashed binary molecular fingerprints and 4-point pharmacophore fingerprints using Tanimoto, Cosine, Dice and Tversky similarity matrices. Molecules screened by DHPM exhibited significant structural diversity, better binding strength and drug like properties as compared to the compounds screened by CPMs indicating the efficiency of DHPM to explore new chemical space for anti-TB drug discovery. The idea of DHPM can be applied for a wide range of mycobacterial or other pathogen targets to venture into unexplored chemical space.


2020 ◽  
Author(s):  
Chinmayee Choudhury ◽  
Anshu Bhardwaj

Abstract Antimicrobial resistance (AMR) is one of the most serious global public health threats as it compromises the successful treatment of deadly infectious diseases like tuberculosis. New therapeutics are constantly needed but it takes a long time and is expensive to explore new biochemical space. One way to address this issue is to repurpose the validated targets and identify novel chemotypes that can simultaneously bind to multiple binding pockets of these targets as a new lead generation strategy. This study reports such a strategy, dynamic hybrid pharmacophore model (DHPM), which represents the combined interaction features of different binding pockets contrary to the conventional approaches, where pharmacophore models are generated from single binding sites. We have considered Mtb-DapB, a validated mycobacterial drug target, as our model system to explore the effectiveness of DHPMs to screen novel unexplored compounds. Mtb-DapB has a cofactor binding site (CBS) and an adjacent substrate binding site (SBS). Four different model systems of Mtb-DapB were designed where, either NADPH/NADH occupies CBS in presence/absence of an inhibitor 2, 6-PDC in the adjacent SBS. Two more model systems were designed, where 2, 6-PDC was linked to NADPH and NADH to form hybrid molecules. The six model systems were subjected to 200ns molecular dynamics simulations and trajectories were analysed to identify stable ligand-receptor interaction features. Based on these interactions, Conventional pharmacophore models (CPM) were generated from the individual binding sites while DHPMs were created from hybrid-molecules occupying both binding sites. A huge library of 15, 63,764 publically available molecules were screened by CPMs and DHPMs. The screened hits obtained from both types of models were compared based on their Hashed binary molecular fingerprints and 4 point pharmacophore fingerprints using Tanimoto, Cosine, Dice and Tversky similarity matrices. Molecules screened by DHPM exhibited significant structural diversity, better binding strength and drug like properties as compared to the compounds screened by CPMs and the reported anti-mycobacterial molecules indicating the efficiency of DHPM to explore new chemical space for anti-TB drug discovery. The idea of DHPM can be applied for a wide range of mycobacterial or other pathogen targets to venture into unexplored chemical space.


2016 ◽  
Vol 113 (5) ◽  
pp. E644-E653 ◽  
Author(s):  
Nurit Degani-Katzav ◽  
Revital Gortler ◽  
Lilach Gorodetzki ◽  
Yoav Paas

The invertebrate glutamate-gated chloride-selective receptors (GluClRs) are ion channels serving as targets for ivermectin (IVM), a broad-spectrum anthelmintic drug used to treat human parasitic diseases like river blindness and lymphatic filariasis. The native GluClR is a heteropentamer consisting of α and β subunit types, with yet unknown subunit stoichiometry and arrangement. Based on the recent crystal structure of a homomeric GluClαR, we introduced mutations at the intersubunit interfaces where Glu (the neurotransmitter) binds. By electrophysiological characterization of these mutants, we found heteromeric assemblies with two equivalent Glu-binding sites at β/α intersubunit interfaces, where the GluClβ and GluClα subunits, respectively, contribute the “principal” and “complementary” components of the putative Glu-binding pockets. We identified a mutation in the IVM-binding site (far away from the Glu-binding sites), which significantly increased the sensitivity of the heteromeric mutant receptor to both Glu and IVM, and improved the receptor subunits’ cooperativity. We further characterized this heteromeric GluClR mutant as a receptor having a third Glu-binding site at an α/α intersubunit interface. Altogether, our data unveil heteromeric GluClR assemblies having three α and two β subunits arranged in a counterclockwise β-α-β-α-α fashion, as viewed from the extracellular side, with either two or three Glu-binding site interfaces.


1997 ◽  
Vol 324 (2) ◽  
pp. 547-553 ◽  
Author(s):  
Hyungtae KIM ◽  
William D. PENNIE ◽  
Yi SUN ◽  
Nancy H. COLBURN

Tissue inhibitor of metalloproteinases-3 (TIMP-3) is an extracellular-matrix-associated protein that suppresses tumorigenicity or invasion in several model systems. We have identified, by in vitro footprinting, six AP-1 (activator protein-1) or AP-1-like binding sites in the mouse TIMP-3 promoter that bind purified c-Jun homodimers. Electrophoretic mobility shift assays revealed that the non-consensus fifth AP-1 binding site (AP-720; nt -720 to -714) had the strongest binding activity for recombinant c-Jun protein, and that the fourth binding site (AP-763; nt -763 to -754) and AP-720 showed strong binding activity for cellular nuclear proteins. Antibody supershift and blocking experiments suggest that AP-720, but not AP-763, binds authentic AP-1 components. Transient transfection reporter assays of deletion constructs showed that the region spanning AP-720 has the highest transcriptional activity, and that sequences 5′ to this region (nt -2846 to -747) may contain negative regulatory elements. The deletion construct containing about 500 nt 5′ to the transcriptional start, but no AP-1 sites, showed lower but significant activity, suggesting both AP-1-dependent and -independent regulation of the mouse TIMP-3 promoter. Mutational inactivation of AP-720 abolished the activity increment that distinguished the reporter construct containing both AP-720 and sixth AP-1 binding site (AP-617; nt -617 to -611) from that containing only AP-617. In summary, we report here that both AP-1 and non-AP-1 elements contribute to activity, with the non-consensus AP-1 site at -720 showing the greatest functional significance among the AP-1 sites.


2020 ◽  
Vol 95 (1) ◽  
Author(s):  
Tetyana Lukash ◽  
Tatiana Agback ◽  
Francisco Dominguez ◽  
Nikita Shiliaev ◽  
Chetan Meshram ◽  
...  

ABSTRACT Decades of insufficient control have resulted in unprecedented spread of chikungunya virus (CHIKV) around the globe, and millions have suffered from the highly debilitating disease. Nevertheless, the current understanding of CHIKV-host interactions and adaptability of the virus to replication in mosquitoes and mammalian hosts is still elusive. Our new study shows that four-and-a-half LIM domain protein (FHL1) is one of the host factors that interact with the hypervariable domain (HVD) of CHIKV nsP3. Unlike G3BPs, FHL1 is not a prerequisite of CHIKV replication, and many commonly used cell lines do not express FHL1. However, its expression has a detectable stimulatory effect(s) on CHIKV replication, and Fhl1 knockout (KO) cell lines demonstrate slower infection spread. Nuclear magnetic resonance (NMR)-based studies revealed that the binding site of FHL1 in CHIKV nsP3 HVD overlaps that of another proviral host factor, CD2AP. The structural data also demonstrated that FHL1-HVD interaction is mostly determined by the LIM1 domain of FHL1. However, it does not mirror binding of the entire protein, suggesting that other LIM domains are involved. In agreement with previously published data, our biological experiments showed that interactions of CHIKV HVD with CD2AP and FHL1 have additive effects on the efficiency of CHIKV replication. This study shows that CHIKV mutants with extensive modifications of FHL1- or both FHL1- and CD2AP-binding sites remain viable and develop spreading infection in multiple cell types. Our study also demonstrated that other members of the FHL family can bind to CHIKV HVD and thus may be involved in viral replication. IMPORTANCE Replication of chikungunya virus (CHIKV) is determined by a wide range of host factors. Previously, we have demonstrated that the hypervariable domain (HVD) of CHIKV nsP3 contains linear motifs that recruit defined families of host proteins into formation of functional viral replication complexes. Now, using NMR-based structural and biological approaches, we have characterized the binding site of the cellular FHL1 protein in CHIKV HVD and defined the biological significance of this interaction. In contrast to previously described binding of G3BP to CHIKV HVD, the FHL1-HVD interaction was found to not be a prerequisite of viral replication. However, the presence of FHL1 has a stimulatory effect on CHIKV infectivity and, subsequently, the infection spread. FHL1 and CD2AP proteins were found to have overlapping binding sites in CHIKV HVD and additive proviral functions. Elimination of the FHL1-binding site in the nsP3 HVD can be used for the development of stable, attenuated vaccine candidates.


2000 ◽  
Vol 20 (1) ◽  
pp. 389-401 ◽  
Author(s):  
Elisabetta Soldaini ◽  
Susan John ◽  
Stefano Moro ◽  
Julie Bollenbacher ◽  
Ulrike Schindler ◽  
...  

ABSTRACT We have defined the optimal binding sites for Stat5a and Stat5b homodimers and found that they share similar core TTC(T/C)N(G/A)GAA interferon gamma-activated sequence (GAS) motifs. Stat5a tetramers can bind to tandemly linked GAS motifs, but the binding site selection revealed that tetrameric binding also can be seen with a wide range of nonconsensus motifs, which in many cases did not allow Stat5a binding as a dimer. This indicates a greater degree of flexibility in the DNA sequences that allow binding of Stat5a tetramers than dimers. Indeed, in an oligonucleotide that could bind both dimers and tetramers, it was possible to design mutants that affected dimer binding without affecting tetramer binding. A spacing of 6 bp between the GAS sites was most frequently selected, demonstrating that this distance is favorable for Stat5a tetramer binding. These data provide insights into tetramer formation by Stat5a and indicate that the repertoire of potential binding sites for this transcription factor is broader than expected.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7818 ◽  
Author(s):  
Hua Yang ◽  
Yan-Lin Liu ◽  
Yuan-Yuan Tao ◽  
Wei Yang ◽  
Chun-Ping Yang ◽  
...  

The bamboo snout beetle Cyrtotrachelus buqueti is a widely distributed wood-boring pest found in China, and its larvae cause significant economic losses because this beetle targets a wide range of host plants. A potential pest management measure of this beetle involves regulating olfactory chemoreceptors. In the process of olfactory recognition, pheromone-binding proteins (PBPs) play an important role. Homology modeling and molecular docking were conducted in this study for the interaction between CbuqPBP1 and dibutyl phthalate to better understand the relationship between PBP structures and their ligands. Site-directed mutagenesis and binding experiments were combined to identify the binding sites of CbuqPBP1 and to explore its ligand-binding mechanism. The 3D structural model of CbuqPBP1 has six a-helices. Five of these a-helices adopt an antiparallel arrangement to form an internal ligand-binding pocket. When docking dibutyl phthalate within the active site of CbuqPBP1, a CH-π interaction between the benzene ring of dibutyl phthalate and Phe69 was observed, and a weak hydrogen bond formed between the ester carbonyl oxygen and His53. Thus, Phe69 and His53 are predicted to be important residues of CbuqPBP1 involved in ligand recognition. Site-directed mutagenesis and fluorescence assays with a His53Ala CbuqPBP1 mutant showed no affinity toward ligands. Mutation of Phe69 only affected binding of CbuqPBP1 to cedar camphor. Thus, His53 (Between α2 and α3) of CbuqPBP1 appears to be a key binding site residue, and Phe69 (Located at α3) is a very important binding site for particular ligand interactions.


2020 ◽  
Author(s):  
Tetyana Lukash ◽  
Tatiana Agback ◽  
Francisco Dominguez ◽  
Nikita Shiliaev ◽  
Chetan Meshram ◽  
...  

ABSTRACTDecades of insufficient control resulted in unprecedented spread of chikungunya virus (CHIKV) around the globe and millions already suffered from the highly debilitating disease. Nevertheless, the current understanding of CHIKV-host interactions and adaptability of the virus to replication in mosquitoes and mammalian hosts is still elusive. Our new study shows that four-and-a-half LIM domain protein (FHL1) is one of the host factors that interact with hypervariable domain (HVD) of CHIKV nsP3. Unlike G3BPs, FHL1 is not a pre-requisite of CHIKV replication, and many commonly used cell lines do not express FHL1. However, its expression has detectable stimulatory effect(s) on CHIKV replication, and the Fhl1 KO cell lines demonstrate slower infection spread. The NMR-based studies revealed that the binding site of FHL1 in CHIKV nsP3 HVD overlaps with that of another pro-viral host factor, CD2AP. The structural data also demonstrated that FHL1-HVD interaction is mostly determined by LIM1 domain of FHL1. However, it does not mirror binding of the entire protein, suggesting that other LIM domains are involved. In agreement with previously published data, our biological experiments showed that interactions of CHIKV HVD with CD2AP and FHL1 have additive positive effects on the efficiency of CHIKV replication. This study shows that CHIKV mutants with extensive modifications of FHL1- or both FHL1- and CD2AP- binding sites remain viable and develop spreading infection in multiple cell types. Thus, such modifications of HVD may improve live CHIKV vaccine candidates in terms of their safety and stability of the attenuated phenotype.IMPORTANCEReplication of chikungunya virus (CHIKV) is determined by a wide range of host factors. Previously, we have demonstrated that the hypervariable domain (HVD) of CHIKV nsP3 protein contains linear motifs that recruit defined families of host proteins into formation of functional viral replication complexes. Now, using NMR-based structural and biological approaches, we have characterized the binding site of cellular FHL1 protein in CHIKV HVD and defined the biological significance of this interaction. In contrast to previously described binding of G3BP to CHIKV HVD, the FHL1-HVD interaction was found to not be a prerequisite of viral replication. However, the presence of FHL1 has a stimulatory effect on CHIKV infectivity and subsequently, the infection spread. FHL1 and CD2AP proteins were found to have overlapping binding sites in CHIKV HVD and additive pro-viral functions. Elimination of FHL1-binding site in nsP3 HVD can be used for the development of stable, live attenuated vaccine candidates.


2018 ◽  
Author(s):  
Rajiv Gandhi Govindaraj ◽  
Michal Brylinski

AbstractBackgroundDetecting similar ligand-binding sites in globally unrelated proteins has a wide range of applications in modern drug discovery, including drug repurposing, the prediction of side effects, and drug-target interactions. Although a number of techniques to compare binding pockets have been developed, this problem still poses significant challenges.ResultsWe evaluate the performance of three algorithms to calculate similarities between ligand-binding sites, APoc, SiteEngine, and G-LoSA. Our assessment considers not only the capabilities to identify similar pockets and to construct accurate local alignments, but also the dependence of these alignments on the sequence order. We point out certain drawbacks of previously compiled datasets, such as the inclusion of structurally similar proteins, leading to an overestimated performance. To address these issues, a rigorous procedure to prepare unbiased, high-quality benchmarking sets is proposed. Further, we conduct a comparative assessment of techniques directly aligning binding pockets to indirect strategies employing structure-based virtual screening with AutoDock Vina and rDock.ConclusionsThorough benchmarks reveal that G-LoSA offers a fairly robust overall performance, whereas the accuracy of APoc and SiteEngine is satisfactory only against easy datasets. Moreover, combining various algorithms into a meta-predictor improves the performance of existing methods to detect similar binding sites in unrelated proteins by 5-10%. All data reported in this paper are freely available at https://osf.io/6ngbs/.


2021 ◽  
Vol 648 ◽  
pp. A84
Author(s):  
G. Molpeceres ◽  
V. Zaverkin ◽  
N. Watanabe ◽  
J. Kästner

Context. Molecular hydrogen (H2) is the most abundant interstellar molecule and plays an important role in the chemistry and physics of the interstellar medium. The interaction of H2 with interstellar ices is relevant for several processes (e.g., nuclear spin conversion and chemical reactions on the surface of the ice). To model surface processes, quantities such as binding energies and sticking coefficients are required. Aims. We provide sticking coefficients and binding energies for the H2/CO system. These data are absent in the literature so far and could help modelers and experimentalists to draw conclusions on the H2/CO interaction in cold molecular clouds. Methods. Ab initio molecular dynamics simulations, in combination with neural network potentials, were employed in our simulations. Atomistic neural networks were trained against density functional theory calculations on model systems. We sampled a wide range of H2 internal energies and three surface temperatures. Results. Our results show that the binding energy for the H2/CO system is low on average, − 157 K for amorphous CO and −266 K for crystalline CO. This carries several implications for the rest of the work. H2 binding to crystalline CO is stronger by 109 K than to amorphous CO, while amorphous CO shows a wider H2 binding energy distribution. Sticking coefficients are never unity and vary strongly with surface temperature, but less so with ice phase, with values between 0.95 and 0.17. With the values of this study, between 17 and 25% of a beam of H2 molecules at room temperature would stick to the surface, depending on the temperature of the surface and the ice phase. Residence times vary by several orders of magnitude between crystalline and amorphous CO, with the latter showing residence times on the order of seconds at 5 K. H2 may diffuse before desorption in amorphous ices, which might help to accommodate it in deeper binding sites. Conclusions. Based on our results, a significant fraction of H2 molecules will stick on CO ice under experimental conditions, even more so under the harsh conditions of prestellar cores. However, with the low H2–CO binding energies, residence times of H2 on CO ice before desorption are too short to consider a significant population of H2 molecules on pure CO ices. Diffusion is possible in a time window before desorption, which might help accommodate H2 on deeper binding sites, which would increase residence times on the surface.


Sign in / Sign up

Export Citation Format

Share Document