Ligand-based in silico target fishing can be used to identify the potential interacting target of bioactive ligands, which is useful for understanding the polypharmacology and safety profile of existing drugs. The underlying principle of the approach is that known bioactive ligands can be used as reference to predict the targets for a new compound.

Glioblastoma multiforme (GBM) is a lethal, therapy-resistant brain cancer consisting of numerous tumor cell subpopulations, including stem-like glioma-initiating cells (GICs), which contribute to tumor recurrence following initial response to therapy. Here, we identified miR-182 as a regulator of apoptosis, growth, and differentiation programs whose expression level is correlated with GBM patient survival. Repression of Bcl2-like12 (Bcl2L12), c-Met, and hypoxia-inducible factor 2α (HIF2A) is of central importance to miR-182 anti-tumor activity, as it results in enhanced therapy susceptibility, decreased GIC sphere size, expansion, and stemness in vitro. To evaluate the tumor-suppressive function of miR-182 in vivo, we synthesized miR-182-based spherical nucleic acids (182-SNAs); i.e., gold nanoparticles covalently functionalized with mature miR-182 duplexes. Intravenously administered 182-SNAs penetrated the blood-brain/blood-tumor barriers (BBB/BTB) in orthotopic GBM xenografts and selectively disseminated throughout extravascular glioma parenchyma, causing reduced tumor burden and increased animal survival. Our results indicate that harnessing the anti-tumor activities of miR-182 via safe and robust delivery of 182-SNAs represents a novel strategy for therapeutic intervention in GBM.

To assess telomerase as a cancer therapeutic target and determine adaptive mechanisms to telomerase inhibition, we modeled telomerase reactivation and subsequent extinction in T cell lymphomas arising in Atm(-/-) mice engineered with an inducible telomerase reverse transcriptase allele. Telomerase reactivation in the setting of telomere dysfunction enabled full malignant progression with alleviation of telomere dysfunction-induced checkpoints. These cancers possessed copy number alterations targeting key loci in human T cell lymphomagenesis. Upon telomerase extinction, tumor growth eventually slowed with reinstatement of telomere dysfunction-induced checkpoints, yet growth subsequently resumed as tumors acquired alternative lengthening of telomeres (ALT) and aberrant transcriptional networks centering on mitochondrial biology and oxidative defense. ALT+ tumors acquired amplification/overexpression of PGC-1β, a master regulator of mitochondrial biogenesis and function, and they showed marked sensitivity to PGC-1β or SOD2 knockdown. Genetic modeling of telomerase extinction reveals vulnerabilities that motivate coincidental inhibition of mitochondrial maintenance and oxidative defense mechanisms to enhance antitelomerase cancer therapy.

The fibroblast growth factor/fibroblast growth factor receptor (FGF/FGFR) signaling pathway plays crucial roles in cell proliferation, angiogenesis, migration, and survival. Aberration in FGFRs correlates with several malignancies and disorders. FGFRs have proved to be attractive targets for therapeutic intervention in cancer, and it is of high interest to find FGFR inhibitors with novel scaffolds. In this study, a combinatorial three-dimensional quantitative structure-activity relationship (3D-QSAR) model was developed based on previously reported FGFR1 inhibitors with diverse structural skeletons. This model was evaluated for its prediction performance on a diverse test set containing 232 FGFR inhibitors, and it yielded a SD value of 0.75 pIC50 units from measured inhibition affinities and a Pearson's correlation coefficient R2 of 0.53. This result suggests that the combinatorial 3D-QSAR model could be used to search for new FGFR1 hit structures and predict their potential activity. To further evaluate the performance of the model, a decoy set validation was used to measure the efficiency of the model by calculating EF (enrichment factor). Based on the combinatorial pharmacophore model, a virtual screening against SPECS database was performed. Nineteen novel active compounds were successfully identified, which provide new chemical starting points for further structural optimization of FGFR1 inhibitors.

Eukaryotic translation initiation factor eIF4AI, the founding member of DEAD-box helicases, undergoes ATP hydrolysis-coupled conformational changes to unwind mRNA secondary structures during translation initiation. However, the mechanism of its coupled enzymatic activities remains unclear. Here we report that a gating mechanism for Pi release controlled by the inter-domain linker of eIF4AI regulates the coupling between ATP hydrolysis and RNA unwinding. Molecular dynamic simulations and experimental results revealed that, through forming a hydrophobic core with the conserved SAT motif of the N-terminal domain and I357 from the C-terminal domain, the linker gated the release of Pi from the hydrolysis site, which avoided futile hydrolysis cycles of eIF4AI. Further mutagenesis studies suggested this linker also plays an auto-inhibitory role in the enzymatic activity of eIF4AI, which may be essential for its function during translation initiation. Overall, our results reveal a novel regulatory mechanism that controls eIF4AI-mediated mRNA unwinding and can guide further mechanistic studies on other DEAD-box helicases.

In this work, the relationship between cyclophilin A (CypA) and HCV prompted us to screen a series of small molecule CypA inhibitors which were previously reported by our group. Among them, compound 1, discovered as a non-immunosuppressive anti-HCV agent with an EC50 value of 0.67 μM in a virus assay, was selected for further study. Subsequent chemical modification by O-acylation led to a novel class of molecules, among which compound 25 demonstrated the most potent anti-HCV activity in the virus assay (EC50 = 0.19 μM), but low cytotoxicity and hERG cardiac toxicity. The following studies (a solution stability assay and a simple pharmacokinetic test together with a CypA enzyme inhibition assay) preliminarily indicated that 25 was a prodrug of 1. To the best of our knowledge, 25 is probably the most potent currently reported small molecule anti-HCV agent acting via the CypA inhibitory mechanism. Consequently, our study has provided a new potential small molecule for curing HCV infection.

The 2009 H1N1 influenza pandemic (pH1N1) led to record sales of neuraminidase (NA) inhibitors, which has contributed significantly to the recent increase in oseltamivir-resistant viruses. Therefore, development and careful evaluation of novel NA inhibitors is of great interest. Recently, a highly potent NA inhibitor, laninamivir, has been approved for use in Japan. Laninamivir is effective using a single inhaled dose via its octanoate prodrug (CS-8958) and has been demonstrated to be effective against oseltamivir-resistant NA in vitro. However, effectiveness of laninamivir octanoate prodrug against oseltamivir-resistant influenza infection in adults has not been demonstrated. NA is classified into 2 groups based upon phylogenetic analysis and it is becoming clear that each group has some distinct structural features. Recently, we found that pH1N1 N1 NA (p09N1) is an atypical group 1 NA with some group 2-like features in its active site (lack of a 150-cavity). Furthermore, it has been reported that certain oseltamivir-resistant substitutions in the NA active site are group 1 specific. In order to comprehensively evaluate the effectiveness of laninamivir, we utilized recombinant N5 (typical group 1), p09N1 (atypical group 1) and N2 from the 1957 pandemic H2N2 (p57N2) (typical group 2) to carry out in vitro inhibition assays. We found that laninamivir and its octanoate prodrug display group specific preferences to different influenza NAs and provide the structural basis of their specific action based upon their novel complex crystal structures. Our results indicate that laninamivir and zanamivir are more effective against group 1 NA with a 150-cavity than group 2 NA with no 150-cavity. Furthermore, we have found that the laninamivir octanoate prodrug has a unique binding mode in p09N1 that is different from that of group 2 p57N2, but with some similarities to NA-oseltamivir binding, which provides additional insight into group specific differences of oseltamivir binding and resistance.

Glucose transporter 4 (GLUT4) is a principal glucose transporter in response to insulin, and impaired translocation or decreased expression of GLUT4 is believed to be one of the major pathological features of type 2 diabetes mellitus (T2DM). Therefore, induction of GLUT4 translocation or/and expression is a promising strategy for anti-T2DM drug discovery. Here we report that the natural product (+)-Rutamarin (Rut) functions as an efficient dual inducer on both insulin-induced GLUT4 translocation and expression. Rut-treated 3T3-L1 adipocytes exhibit efficiently enhanced insulin-induced glucose uptake, while diet-induced obese (DIO) mice based assays further confirm the Rut-induced improvement of glucose homeostasis and insulin sensitivity in vivo. Subsequent investigation of Rut acting targets indicates that as a specific protein tyrosine phosphatase 1B (PTP1B) inhibitor Rut induces basal GLUT4 translocation to some extent and largely enhances insulin-induced GLUT4 translocation through PI3 kinase-AKT/PKB pathway, while as an agonist of retinoid X receptor α (RXRα), Rut potently increases GLUT4 expression. Furthermore, by using molecular modeling and crystallographic approaches, the possible binding modes of Rut to these two targets have been also determined at atomic levels. All our results have thus highlighted the potential of Rut as both a valuable lead compound for anti-T2DM drug discovery and a promising chemical probe for GLUT4 associated pathways exploration.

The severe acute respiratory syndrome coronavirus (SARS-CoV) membrane protein is an abundant virion protein, and its interaction with the nucleocapsid protein is crucial for viral assembly and morphogenesis. Although the interacting region in the nucleocapsid protein was mapped to residues 168-208, the interacting region in the membrane protein and the interaction nature are still unclear. In this work, by using yeast two-hybrid and surface plasmon resonance techniques, the residues 197-221 of the membrane protein and the residues 351-422 of the nucleocapsid protein were determined to be involved in their interaction. Sequence analysis revealed that these two fragments are highly charged at neutral pH, suggesting that their interaction may be of ionic nature. Kinetic assays indicated that the endodomain (aa102-221) of the membrane protein interacts with the nucleocapsid protein with high affinity (K(D)=0.55+/-0.04 microM), however, this interaction could be weakened greatly by acidification, higher salt concentration (400 mM NaCl) and divalent cation (50 mM Ca2+), which suggests that electrostatic attraction might play an important role in this interaction. In addition, it is noted that two highly conserved amino acids (L218 and L219) in the membrane protein are not involved in this interaction. Here, we show that electrostatic interactions between the carboxyl termini of SARS-CoV membrane protein and nucleocapsid protein largely mediate the interaction of these two proteins. These results might facilitate therapeutic strategies aiming at the disruption of the association between SARS-CoV membrane and nucleocapsid proteins.

Severe acute respiratory syndrome coronavirus (SARS-CoV) is responsible for SARS infection. Nucleocapsid protein (NP) of SARS-CoV (SARS_NP) functions in enveloping the entire genomic RNA and interacts with viron structural proteins, thus playing important roles in the process of virus particle assembly and release. Protein-protein interaction analysis using bioinformatics tools indicated that SARS_NP may bind to human cyclophilin A (hCypA), and surface plasmon resonance (SPR) technology revealed this binding with the equilibrium dissociation constant ranging from 6 to 160nM. The probable binding sites of these two proteins were detected by modeling the three-dimensional structure of the SARS_NP-hCypA complex, from which the important interaction residue pairs between the proteins were deduced. Mutagenesis experiments were carried out for validating the binding model, whose correctness was assessed by the observed effects on the binding affinities between the proteins. The reliability of the binding sites derived by the molecular modeling was confirmed by the fact that the computationally predicted values of the relative free energies of the binding for SARS_NP (or hCypA) mutants to the wild-type hCypA (or SARS_NP) are in good agreement with the data determined by SPR. Such presently observed SARS_NP-hCypA interaction model might provide a new hint for facilitating the understanding of another possible SARS-CoV infection pathway against human cell.

Severe acute respiratory syndrome (SARS) coronavirus is a novel human coronavirus and is responsible for SARS infection. SARS coronavirus 3C-like proteinase (SARS 3CL(pro)) plays key roles in viral replication and transcription and is an attractive target for anti-SARS drug discovery. In this report, we quantitatively characterized the dimerization features of the full-length and N-terminal residues 1-7 deleted SARS 3CL(pro)s by using glutaraldehyde cross-linking SDS-PAGE, size-exclusion chromatography, and isothermal titration calorimeter techniques. Glutaraldehyde cross-linking SDS-PAGE and size-exclusion chromatography results show that, similar to the full-length SARS 3CL(pro), the N-terminal deleted SARS 3CL(pro) still remains a dimer/monomer mixture within a wide range of protein concentrations. Isothermal titration calorimeter determinations indicate that the equilibrium dissociation constant (K(d)) of the N-terminal deleted proteinase dimer (262 microm) is very similar to that of the full-length proteinase dimer (227 microm). Enzymatic activity assay using the fluorescence resonance energy transfer method reveals that N-terminal deletion results in almost complete loss of enzymatic activity for SARS 3CL(pro). Molecular dynamics and docking simulations demonstrate the N-terminal deleted proteinase dimer adopts a state different from that of the full-length proteinase dimer, which increases the angle between the two protomers and reduces the binding pocket that is not beneficial to the substrate binding. This conclusion is verified by the surface plasmon resonance biosensor determination, indicating that the model substrate cannot bind to the N-terminal deleted proteinase. These results suggest the N terminus is not indispensable for the proteinase dimerization but may fix the dimer at the active state and is therefore vital to enzymatic activity.

The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer beta-sheet composed of beta9 and beta10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O <--> C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.

SARS-CoV 3C-like protease (3CL(pro)) is an attractive target for anti-severe acute respiratory syndrome (SARS) drug discovery, and its dimerization has been extensively proved to be indispensable for enzymatic activity. However, the reason why the dissociated monomer is inactive still remains unclear due to the absence of the monomer structure. In this study, we showed that mutation of the dimer-interface residue Gly-11 to alanine entirely abolished the activity of SARS-CoV 3CL(pro). Subsequently, we determined the crystal structure of this mutant and discovered a complete crystallographic dimer dissociation of SARS-CoV 3CL(pro). The mutation might shorten the alpha-helix A' of domain I and cause a mis-oriented N-terminal finger that could not correctly squeeze into the pocket of another monomer during dimerization, thus destabilizing the dimer structure. Several structural features essential for catalysis and substrate recognition are severely impaired in the G11A monomer. Moreover, domain III rotates dramatically against the chymotrypsin fold compared with the dimer, from which we proposed a putative dimerization model for SARS-CoV 3CL(pro). As the first reported monomer structure for SARS-CoV 3CL(pro), the crystal structure of G11A mutant might provide insight into the dimerization mechanism of the protease and supply direct structural evidence for the incompetence of the dissociated monomer.

Growth hormone-releasing hormone (GHRH) regulates the secretion of growth hormone that virtually controls metabolism and growth of every tissue through its binding to the cognate receptor (GHRHR). Malfunction in GHRHR signaling is associated with abnormal growth, making GHRHR an attractive therapeutic target against dwarfism (e.g., isolated growth hormone deficiency, IGHD), gigantism, lipodystrophy and certain cancers. Here, we report the cryo-electron microscopy (cryo-EM) structure of the human GHRHR bound to its endogenous ligand and the stimulatory G protein at 2.6 Å. This high-resolution structure reveals a characteristic hormone recognition pattern of GHRH by GHRHR, where the α-helical GHRH forms an extensive and continuous network of interactions involving all the extracellular loops (ECLs), all the transmembrane (TM) helices except TM4, and the extracellular domain (ECD) of GHRHR, especially the N-terminus of GHRH that engages a broad set of specific interactions with the receptor. Mutagenesis and molecular dynamics (MD) simulations uncover detailed mechanisms by which IGHD-causing mutations lead to the impairment of GHRHR function. Our findings provide insights into the molecular basis of peptide recognition and receptor activation, thereby facilitating the development of structure-based drug discovery and precision medicine.

Elucidating the conformational heterogeneity of proteins is essential for understanding protein function and developing exogenous ligands. With the rapid development of experimental and computational methods, it is of great interest to integrate these approaches to illuminate the conformational landscapes of target proteins. SETD8 is a protein lysine methyltransferase (PKMT), which functions in vivo via the methylation of histone and nonhistone targets. Utilizing covalent inhibitors and depleting native ligands to trap hidden conformational states, we obtained diverse X-ray structures of SETD8. These structures were used to seed distributed atomistic molecular dynamics simulations that generated a total of six milliseconds of trajectory data. Markov state models, built via an automated machine learning approach and corroborated experimentally, reveal how slow conformational motions and conformational states are relevant to catalysis. These findings provide molecular insight on enzymatic catalysis and allosteric mechanisms of a PKMT via its detailed conformational landscape.

Leukotriene B4 receptor 1 (BLT1) plays crucial roles in the acute inflammatory responses and is a valuable target for anti-inflammation treatment, however, the mechanism by which leukotriene B4 (LTB4) activates receptor remains unclear. Here, we report the cryo-electron microscopy (cryo-EM) structure of the LTB4 -bound human BLT1 in complex with a Gi protein in an active conformation at resolution of 2.91 Å. In combination of molecule dynamics (MD) simulation, docking and site-directed mutagenesis, our structure reveals that a hydrogen-bond network of water molecules and key polar residues is the key molecular determinant for LTB4 binding. We also find that the displacement of residues M1013.36 and I2717.39 to the center of receptor, which unlock the ion lock of the lower part of pocket, is the key mechanism of receptor activation. In addition, we reveal a binding site of phosphatidylinositol (PI) and discover that the widely open ligand binding pocket may contribute the lack of specificity and efficacy for current BLT1-targeting drug design. Taken together, our structural analysis provides a scaffold for understanding BLT1 activation and a rational basis for designing anti-leukotriene drugs.

The interaction between Lymphocyte function-associated antigen 1 (LFA-1) and intercellular-adhesion molecule-1 (ICAM-1) plays important roles in the cell-mediated immune response and inflammation associated with dry eye disease. LFA-1/ICAM-1 antagonists can be used for the treatment of dry eye disease, such as Lifitegrast which has been approved by the FDA in 2016 as a new drug for the treatment of dry eye disease. In this study, we designed and synthesized some new structure compounds that are analogues to Lifitegrast, and their biological activities were evaluated by in vitro cell-based assay and also by in vivo mouse dry eye model. Our results demonstrated that one of these analogues of Lifitegrast (compound 1b) showed good LFA-1/ICAM-1 antagonist activity in in vitro assay; meanwhile, it also significantly reduced ocular surface epithelial cells damage, increased goblet cell density in dry eye mouse and highly improved the symptoms of dry eye mouse. Graphical abstract.

Identifying the potential compound-protein interactions (CPIs) plays an essential role in drug development. The computational approaches for CPI prediction can reduce time and costs of experimental methods and have benefited from the continuously improved graph representation learning. However, most of the network-based methods use heterogeneous graphs, which is challenging due to their complex structures and heterogeneous attributes. Therefore, in this work, we transformed the compound-protein heterogeneous graph to a homogeneous graph by integrating the ligand-based protein representations and overall similarity associations. We then proposed an Inductive Graph AggrEgator-based framework, named CPI-IGAE, for CPI prediction. CPI-IGAE learns the low-dimensional representations of compounds and proteins from the homogeneous graph in an end-to-end manner. The results show that CPI-IGAE performs better than some state-of-the-art methods. Further ablation study and visualization of embeddings reveal the advantages of the model architecture and its role in feature extraction, and some of the top ranked CPIs by CPI-IGAE have been validated by a review of recent literature. The data and source codes are available at https://github.com/wanxiaozhe/CPI-IGAE.

Glycolytic metabolism enzymes have been implicated in the immunometabolism field through changes in metabolic status. PGK1 is a catalytic enzyme in the glycolytic pathway. Here, we set up a high-throughput screen platform to identify PGK1 inhibitors. DC-PGKI is an ATP-competitive inhibitor of PGK1 with an affinity of K d = 99.08 nmol/L. DC-PGKI stabilizes PGK1 in vitro and in vivo, and suppresses both glycolytic activity and the kinase function of PGK1. In addition, DC-PGKI unveils that PGK1 regulates production of IL-1β and IL-6 in LPS-stimulated macrophages. Mechanistically, inhibition of PGK1 with DC-PGKI results in NRF2 (nuclear factor-erythroid factor 2-related factor 2, NFE2L2) accumulation, then NRF2 translocates to the nucleus and binds to the proximity region of Il-1β and Il-6 genes, and inhibits LPS-induced expression of these genes. DC-PGKI ameliorates colitis in the dextran sulfate sodium (DSS)-induced colitis mouse model. These data support PGK1 as a regulator of macrophages and suggest potential utility of PGK1 inhibitors in the treatment of inflammatory bowel disease.

Eukaryotic messenger RNA (mRNA) typically contains a methylated guanosine (m7G) cap, which mediates major steps of mRNA metabolism. Recently, some RNAs in both prokaryotic and eukaryotic organisms have been found to carry a non-canonical cap such as the NAD cap. Here we report that Arabidopsis DXO family protein AtDXO1, which was previously known to be a decapping enzyme for NAD-capped RNAs (NAD-RNA), is an essential component for m7G capping. AtDXO1 associates with and activates RNA guanosine-7 methyltransferase (AtRNMT1) to catalyze conversion of the guanosine cap to the m7G cap. AtRNMT1 is an essential gene. Partial loss-of-function mutations of AtRNMT1 and knockout mutation of AtDXO1 reduce m7G-capped mRNA but increase G-capped mRNAs, leading to similar pleiotropic phenotypes, whereas overexpression of AtRNMT1 partially restores the atdxo1 phenotypes. This work reveals an important mechanism in m7G capping in plants by which the NAD-RNA decapping enzyme AtDXO1 is required for efficient guanosine cap methylation.