Search

Supplementary Table S1 and Fig S1 and S2. Supplementary Table 1: qPCR primers used in this study; Supplementary Figure S1: Structures of compounds screened in Figure 1A and 1B; Supplementary Figure S2: Proximity ligation assay (PLA) analyzing SIN3A-MAD interactions

Supplementary Table S2. List of genes with expression regulated by selamectin treatment

Triple-negative breast cancers (TNBC) lacking estrogen, progesterone, and HER2 receptors account for 10% to 20% of breast cancer and are indicative of poor prognosis. The development of effective treatment strategies therefore represents a pressing unmet clinical need. We previously identified a molecularly targeted approach to target aberrant epigenetics of TNBC using a peptide corresponding to the SIN3 interaction domain (SID) of MAD. SID peptide selectively blocked binding of SID-containing proteins to the paired α-helix (PAH2) domain of SIN3, resulting in epigenetic and transcriptional modulation of genes associated with epithelial–mesenchymal transition (EMT). To find small molecule inhibitor (SMI) mimetics of SID peptide, we performed an in silico screen for PAH2 domain–binding compounds. This led to the identification of the avermectin macrocyclic lactone derivatives selamectin and ivermectin (Mectizan) as candidate compounds. Both selamectin and ivermectin phenocopied the effects of SID peptide to block SIN3–PAH2 interaction with MAD, induce expression of CDH1 and ESR1, and restore tamoxifen sensitivity in MDA-MB-231 human and MMTV-Myc mouse TNBC cells in vitro. Treatment with selamectin or ivermectin led to transcriptional modulation of genes associated with EMT and maintenance of a cancer stem cell phenotype in TNBC cells. This resulted in impairment of clonogenic self-renewal in vitro and inhibition of tumor growth and metastasis in vivo. Underlining the potential of avermectins in TNBC, pathway analysis revealed that selamectin also modulated the expression of therapeutically targetable genes. Consistent with this, an unbiased drug screen in TNBC cells identified selamectin-induced sensitization to a number of drugs, including those targeting modulated genes. Mol Cancer Ther; 14(8); 1824–36. ©2015 AACR.

Biophysical studies have established that the thyrotropin (TSH) receptor (TSHR) undergoes posttranslational modifications including dimerization. Following our earlier simulation of a TSHR–transmembrane domain (TMD) monomer (called TSHR-TMD-TRIO) we have now proceeded with a molecular dynamics simulation (MD) of TSHR-TMD dimerization using this improved membrane-embedded model. The starting structure was the TMD protein with all extracellular and intracellular loops and internal waters, which was placed in the relative orientation of the model originally generated with Brownian dynamics. Furthermore, this model was embedded in a DPPC lipid bilayer further solvated with water and added salt. Data from the MD simulation studies showed that the dimeric subunits stayed in the same relative orientation and distance during the 1000 ns of study. Comparison of representative conformations of the individual monomers when dimerized with the conformations from the monomer simulation showed subtle differences as represented by the backbone root mean square deviations. Differences in the conformations of the ligand-binding sites, suggesting variable affinities for these “hot spots,” were also revealed by comparing the docking scores of 46 small-molecule ligands that included known TSHR agonists and antagonists as well as their derivatives. These data add further insight into the tendency of the TSHR-TMD to form dimeric and oligomeric structures and show that the differing conformations influence small-molecule binding sites within the TMD.

The virtual high-throughput screening (vHTS) approach has been widely used for large database screening to identify potential lead compounds for drug discovery. Due to its high computational demands, docking that allows receptor flexibility has been a challenging problem for virtual screening. Therefore, the selection of protein target conformations is crucial to produce useful vHTS results. Since only a single protein structure is used to screen large databases in most vHTS studies, the main challenge is to reduce false negative rates in selecting compounds for in vitro tests. False negatives are most likely to occur when using apo structures or homology models of protein targets due to the small volume of the binding pocket formed by incorrect side-chain conformations. Even holo protein structures can exhibit high false negative rates due to ligand-induced fit effects, since the shape of the binding pocket highly depends on its bound ligand. To reduce false negative rates and improve success rates for vHTS in drug discovery, we have developed a new Monte Carlo-based approach that optimizes the binding pocket of protein targets. This newly developed Monte Carlo pocket optimization (MCPO) approach was assessed on several datasets showing promising results. The binding pocket optimization approach could be a useful tool for vHTS-based drug discovery, especially in cases when only apo structures or homology models are available.

For the characterization of various aspects of protein structures, four useful concepts are discussed: chameleon sequences, circular variance, mutual proximity, and a subsequence-based foldability score. These concepts were used in estimating foldability of globular, intrinsically disordered and fold-switching proteins, properties of protein-protein interfaces, quantifying sphericity, helping to improve protein-protein docking scores, and estimating the effect of mutations on stability. A conjecture about the Achilles' heel of proteins is presented as well.

Anion Exchanger 1 (AE1, SLC4A1) is the primary bicarbonate (HCO3-) transporter expressed in erythrocyte membranes where it mediates transport of CO2between lungs and other tissues via import/export of bicarbonate. It is also a key regulator of erythrocyte structure and antigenic recognition. Previous biochemical studies, and a low-resolution crystal structure of the transmembrane domain have provided initial insight into AE1 structure and function. However, key questions remain regarding substrate binding and transport as well as the mechanism of inhibition. The orientation of the intracellular domain as well as the localization of lipid and sterol binding sites also remain enigmatic. We herein present seven novel high resolution cryo-EM structures of the full length human transporter in the apo, bicarbonate-bound, and several inhibitor-bound states combined with uptake- and computational studies. To our knowledge, these studies represent the first full length human, and substrate bound, SLC4 transporter structure. Our results reveal important molecular details about substrate binding and transport, as well as the diverse mechanisms of AE1 inhibition by both research chemicals and prescription drugs. We also provide novel insights into the full-length transporter architecture, identify the conformational space of the Diego blood antigen system and elucidate multiple lipid and sterol binding sites.

Semi-empirical calculations on the torsion potential for dihedral angle of rivastigmine linking NAP and Carbamyle moieties consistently show regions of several discontinuities and a cusp indicating molecular instability and eventual break-up of rivastigmine observed in the X-ray structure. The phenomena can be explained both by definition of large classical force or quantum nature of chemical bond break-up. Also, to better understand the molecular properties and quantum energetics of the inhibitor molecules, we have performed several ab initio based calculations on all four inhibitors at equilibrium geometry, in ground state and gas phase using the density functional theory level wB97X/6-31G* and HF/6-31G*. A number of properties like computational vibrational (IR), Raman and nuclear magnetic resonance (NMR) spectra as well as HOMO and LUMO orbital energies at optimized geometries have been computed by SPARTAN16 and Gaussian16 utilities. Also, the thermodynamic and QSAR properties of the inhibitors have been assessed and compared by a number of different semi-quantum, Hartree-Fock and density functional methods. The theoretical NMR and IR spectra have been benchmarked against experimental spectrum to compare and assess suitability of the computational methodologies and basis set levels for the calculations.Communicated by Ramaswamy H. Sarma.

The TSH receptor (TSHR) is the major autoantigen in Graves’ disease (GD). Bioinformatic analyses predict the existence of several human TSHR isoforms from alternative splicing, which can lead to the coexpression of multiple receptor forms. The most abundant of these is TSHRv1.3. In silico modeling of TSHRv1.3 demonstrated the structural integrity of this truncated receptor isoform and its potential binding of TSH. Tissue profiling revealed wide expression of TSHRv1.3, with a predominant presence in thyroid, bone marrow, thymus, and adipose tissue. To gain insight into the role of this v1.3 receptor isoform in thyroid pathophysiology, we cloned the entire open reading frame into a mammalian expression vector. Immunoprecipitation studies demonstrated that both TSHR-stimulating antibody and human TSH could bind v1.3. Furthermore, TSHRv1.3 inhibited the stimulatory effect of TSH and TSHR-Ab MS-1 antibody on TSHR-induced cAMP generation in a dose-dependent manner. To confirm the antigenicity of v1.3, we used a peptide ELISA against two different epitopes. Of 13 GD samples, 11 (84.6%) were positive for a carboxy terminal peptide and 10 (76.9%) were positive with a junction region peptide. To demonstrate that intracellular v1.3 could serve as an autoantigen and modulate disease, we used double-transfected Chinese hamster ovary cells that expressed both green fluorescent protein (GFP)-tagged TSHRv1.3 and full-length TSHR. We then induced cell stress and apoptosis using a TSHR monoclonal antibody and observed the culture supernatant contained v1.3-GFP protein, demonstrating the release of the intracellular receptor variant by this mechanism.

We describe the discovery of an agonist of the nuclear receptor NR2F1 that specifically activates dormancy programs in malignant cells. The agonist led to a self-regulated increase in NR2F1 mRNA and protein and downstream transcription of a novel dormancy program. This program led to growth arrest of an HNSCC PDX line, human cell lines, and patient-derived organoids in 3D cultures and in vivo. This effect was lost when NR2F1 was knocked out by CRISPR-Cas9. RNA sequencing revealed that agonist treatment induces transcriptional changes associated with inhibition of cell cycle progression and mTOR signaling, metastasis suppression, and induction of a neural crest lineage program. In mice, agonist treatment resulted in inhibition of lung HNSCC metastasis, even after cessation of the treatment, where disseminated tumor cells displayed an NR2F1hi/p27hi/Ki-67lo/p-S6lo phenotype and remained in a dormant single-cell state. Our work provides proof of principle supporting the use of NR2F1 agonists to induce dormancy as a therapeutic strategy to prevent metastasis.

The thyroid-stimulating hormone receptor (TSHR) is a G-protein-coupled receptor group A family member with 7 transmembrane helices. We generated 3 new models of its entire transmembrane region using a 600 ns molecular simulation. The simulation started from our previously published model, which we have now revised by also modeling the intracellular loops and the C-terminal tail, adding internal waters and embedding it into a lipid bilayer with a water layer and with ions added to complete the system. We have named this model TSHR-TMD-TRIO since 3 representative dominant structures were then extracted from the simulation trajectory and compared with the original model. These structures each showed small but significant changes in the relative positions of the helices. The 3 models were also used as targets to dock a set of small molecules that are known active compounds including a new TSHR antagonist (BT362), which confirmed the appropriateness of the model with some small molecules showing significant preference for one or other of the structures.

Disseminated tumor cells (DTCs) in secondary organs often remain dormant for a long period of time before re-awakening and growing into overt metastases. We have previously identified NR2F1/COUP-TF1, an orphan nuclear receptor, as a master regulator of tumor cell dormancy in head and neck squamous cell carcinoma (HNSCC) and other cancer types. Here we describe the identification and function of a novel NR2F1 agonist herein referred to as compound 26 (C26). C26 was found to specifically activate NR2F1 in HNSCC cells, leading to increased NR2F1 transcription and nuclear protein accumulation. C26-mediated activation of NR2F1 induced growth arrest of HNSCC PDX line and cell lines in 3D cultures in vitro and on chicken embryo chorioallantoic membrane (CAM) in vivo. The effect of C26 on growth arrest was lost when NR2F1 was knocked out by CRISPR/Cas9. C26-induced growth arrest was mediated by activation of an NR2F1-regulated dormancy program, including upregulation of cyclin-dependent kinase (CDK) inhibitor p27 and the transcription factors retinoic acid receptor β (RARβ) and Sox9. In mice bearing HNSCC PDX tumors, combined adjuvant and neo-adjuvant treatment with C26 resulted in complete inhibition of lethal lung metastasis. Mechanistic analysis showed that lung DTCs in C26-treated mice displayed an NR2F1hi/p27hi/Ki-67lo phenotype, which kept them dormant in a single-cell state preventing their outgrowth into overt metastases. Our work reveals a novel NR2F1 agonist and provides a proof of principle strategy supporting that inducing DTC dormancy using NR2F1 agonists could be used as a therapeutic strategy to prevent metastasis.

Molecular simulations from small molecules to large bio-macromolecules and polymer systems are routinely used to simulate thermodynamics properties of interests by molecular mechanics-based potentials. In a recent paper, via three different semi-empirical methods, we reported quantum singularities in molecular mechanics torsion potentials as signature of chemical bond break-up process revealed under experimental X-ray as broken chemical moieties. In this present work, applying first principle methods of Hartree-Fock, Density Functional as well as Moller-Plesset techniques we have reconfirmed the previous general predictions of singularities in the torsion potential for the case of water dimer that connects two water monomers by weak hydrogen bond. Due to quantum nature of chemical bond breaking process leading to break-point conditions in otherwise connected molecular topology aided by molecular mechanics based potentials, the singularities are also suggestive of large forces as onset in the bond-breaking process. We have presented the details of these novel interesting findings in this paper. These results of quantum singularities can have significant impacts to improve current force fields and can open up new areas we define as “Fracture Molecular Mechanics” or “Fracture Force Field” in overlap regions of molecular and quantum mechanics based approaches to explore and account for chemical bond-breaking mechanisms in molecular simulation techniques.

To gain further insight into the binding of the normal and variant human TSHβ subunits (TSHβ and TSHβv), we modeled the 2 monomeric proteins and studied their interaction with the TSH receptor ectodomain (TSHR-ECD) using molecular dynamics simulation Furthermore, analyzed their bioactivity in vitro using recombinant proteins to confirm that such binding was physiologically relevant. Examining the interaction of TSHβ and TSHβv with the TSHR-ECD model using molecular dynamic simulation revealed strong binding of these proteins to the receptor ECD. The specificity of TSHβ and TSHβv binding to the TSHR-ECD was examined by analyzing the hydrogen-bonding residues of these subunits to the FSH receptor ECD, indicating the inability of these molecules to bind to the FSH receptors. Furthermore, the modelling suggests that TSHβ and TSHβv proteins clasped the concave surface of the leucine rich region of the TSHR ECD in a similar way to the native TSH using dynamic hydrogen bonding. These mutually exclusive stable interactions between the subunits and ECD residues included some high-affinity contact sites corresponding to binding models of native TSH. Furthermore, we cloned TSHβ and TSHβv proteins using the entire coding ORF and purified the flag-tagged proteins. The expressed TSHβ subunit proteins retained bioactivity both in a coculture system as well as with immune-purified proteins. In summary, we showed that such interactions can result in a functional outcome and may exert physiological or pathophysiological effects in immune cells.

It has been shown recently by Porter & Looger that a significant number of proteins exist that can form more than one stable fold. This note examines the sequences of these fold-switching proteins by (a) calculating their foldability scores recently introduced by the present author and (b) comparing the propensity of chameleon sequences in fold switchers and in non-fold switchers. It has been found that the average foldability score of the fold switchers indicates weaker foldability. As for the propensity of chameleon sequences of length 5-7 it was found, somewhat surprisingly, that there is only a very small difference between the fold switchers and the non-fold switchers. Furthermore, when looking at amino acid propensities, for several amino acids there was even an opposing trend in the deviation of their propensities from the overall amino acid propensities. This article is protected by copyright. All rights reserved.

Several properties of amino acid sequences corresponding to proteins that are known to fold are compared to those of randomly generated sequences and to sequences of intrinsically disordered proteins in order to find properties that distinguish folding sequences from the rest. The properties studied included helix and sheet propensities from secondary structure prediction, adjacency correlations, directionality correlations, as well as propensities of all possible triplets and quadruplets. Small differences between known folded and random sequences were observed for the adjacency and directional correlations, and significant differences were seen on the triplet and especially on the quadruplet propensities. Based on the differences in the adjacency, triplet or quadruplet propensities folding scores were defined and used to test the accuracy of foldability prediction based on these statistics. The best predictions were obtained from the quadruplet propensities.

We have previously shown that TSH has modulatory influences on osteoblasts and osteoclasts and can act as a bone protection ligand. We have also shown that a novel TSH-β variant protein is produced locally by murine and human bone marrow macrophages (Endocrinology 154:4919-26, 2013 and Endocrinology 157:3658-67, 2016). Furthermore, we have also shown by molecular modelling that both TSH-β and TSH-βv can bind to the native TSH receptor. In order to characterize the biological activity of this TSH-β variant we have inserted the gene sequence together with a FLAG tag into the P-secTag2 vector which was then cloned into mammalian HEK-293 cells. Western blot analysis identified 17KD and a 12KD proteins corresponding to the wild type TSH-β and the TSH-βv nucleotide sequence in cell culture lysates and, less so in supernatants, by using a FLAG monoclonal antibody and a TSH-β polyclonal antibody. Protein bands of ~34 KD and 24 KD were observed only in the absence of DTT suggesting dimeric forms of TSH-β and TSH-βv respectively. Detection by Coomassie staining followed by sequence analysis using mass spectrometry confirmed the identity of these proteins after column purification. Furthermore, cyclic AMP responses were generated by both supernatant preparations with an 8-fold increase in luciferase activity with TSH-β and a 4-fold increase with TSH-βv. These data demonstrate the inherent bioactivity of single TSH-β chains which was mirrored by TSH-βv. Local paracrine action by macrophage derived TSH-βv on osteoblasts and osteoclasts is, therefore, most likely and raises the possibility of this variant form having significant osteoprotective activity.

The recently developed statistical measure for the type of residue-residue contact at protein complex interfaces, based on a parameter-free definition of contact, has been used to define a contact score that is correlated with the likelihood of correctness of a proposed complex structure. Comparing the proposed contact scores on the native structure and on a set of model structures the proposed measure was shown to generally favor the native structure but in itself was not able to reliably score the native structure to be the best. Adjusting the scores of redocking experiments with the contact score showed that the adjusted score was able to move up the ranking of the native-like structure among the proposed complexes when the native-like was not ranked the best by the respective program. Tests on docking of unbound proteins compared the contact scores of the complexes with the contact score of the crystal structure again showing the tendency of the contact score to favor native-like conformations. The possibility of using the contact score to improve the determination of biological dimers in a crystal structure was also explored. Proteins 2017; 85:235-241. © 2016 Wiley Periodicals, Inc.

The objective of this study was to develop an in silico screening model for characterization of potential novel ligands from commercial drug libraries able to functionally activate certain olfactory receptors (ORs), which are members of the class A rhodopsin-like family of G protein couple receptors (GPCRs), in the brain of murine models of concussion. We previously found that concussions may significantly influence expression of certain ORs, for example, OR4M1 in subjects with a history of concussion/traumatic brain injury (TBI). In this study, we built a 3-D OR4M1 model and used it in in silico screening of potential novel ligands from commercial drug libraries. We report that in vitro activation of OR4M1 with the commercially available ZINC library compound 10915775 led to a significant attenuation of abnormal tau phosphorylation in embryonic cortico-hippocampal neuronal cultures derived from NSE-OR4M1 transgenic mice, possibly through modulation of the JNK signaling pathway. The attenuation of abnormal tau phosphorylation was rather selective since ZINC10915775 significantly decreased tau phosphorylation on tau Ser202/T205 (AT8 epitope) and tau Thr212/Ser214 (AT100 epitope), but not on tau Ser396/404 (PHF-1 epitope). Moreover, no response of ZINC10915775 was found in control hippocampal neuronal cultures derived from wild type littermates. Our in silico model provides novel means to pharmacologically modulate select ubiquitously expressed ORs in the brain through high affinity ligand activation to prevent and eventually to treat concussion induced down regulation of ORs and subsequent cascade of tau pathology. J. Cell. Biochem. 117: 2241-2248, 2016. © 2016 Wiley Periodicals, Inc.

We previously showed that an HLA-DR variant containing arginine at position 74 of the DRβ1 chain (DRβ1-Arg74) is the specific HLA class II variant conferring risk for autoimmune thyroid diseases (AITD). We also identified 5 thyroglobulin (Tg) peptides that bound to DRβ1-Arg74. We hypothesized that blocking the binding of these peptides to DRβ1-Arg74 could block the continuous T-cell activation in thyroiditis needed to maintain the autoimmune response to the thyroid. The aim of the current study was to identify small molecules that can block T-cell activation by Tg peptides presented within DRβ1-Arg74 pockets. We screened a large and diverse library of compounds and identified one compound, cepharanthine that was able to block peptide binding to DRβ1-Arg74. We then showed that Tg.2098 is the dominant peptide when inducing experimental autoimmune thyroiditis (EAT) in NOD mice expressing human DRβ1-Arg74. Furthermore, cepharanthine blocked T-cell activation by thyroglobulin peptides, in particular Tg.2098 in mice that were induced with EAT. For the first time we identified a small molecule that can block Tg peptide binding and presentation to T-cells in autoimmune thyroiditis. If confirmed cepharanthine could potentially have a role in treating human AITD.

Over the past four decades, the number of people with Type 1 Diabetes (T1D) has increased by 4% per year, making it an important public health challenge. Currently, no curative therapy exists for T1D and the only available treatment is insulin replacement. HLA-DQ8 has been shown to present antigenic islet peptides driving the activation of CD4+ T-cells in T1D patients. Specifically, the insulin peptide InsB:9-23 activates self-reactive CD4+ T-cells, causing pancreatic beta cell destruction. The aim of the current study was to identify retro-inverso-D-amino acid based peptides (RI-D-peptides) that can suppress T-cell activation by blocking the presentation of InsB:9-23 peptide within HLA-DQ8 pocket. We identified a RI-D-peptide (RI-EXT) that inhibited InsB:9-23 binding to recombinant HLA-DQ8 molecule, as well as its binding to DQ8 expressed on human B-cells. RI-EXT prevented T-cell activation in a cellular antigen presentation assay containing human DQ8 cells loaded with InsB:9-23 peptide and murine T-cells expressing a human T-cell receptor specific for the InsB:9-23–DQ8 complex. Moreover, RI-EXT blocked T-cell activation by InsB:9-23 in a humanized DQ8 mice both ex vivo and in vivo, as shown by decreased production of IL-2 and IFN-γ and reduced lymphocyte proliferation. Interestingly, RI-EXT also blocked lymphocyte activation and proliferation by InsB:9-23 in PBMCs isolated from recent onset DQ8-T1D patients. In summary, we discovered a RI-D-peptide that blocks InsB:9-23 binding to HLA-DQ8 and its presentation to T-cells in T1D. These findings set the stage for using our approach as a novel therapy for patients with T1D and potentially other autoimmune diseases.