In this study, we identified a novel pyrazole-based derivative (P3C) that displayed potent cytotoxicity against 27 human cancer cell lines derived from different tissue origins with 50% cytotoxic concentrations (CC₅₀) in the low micromolar and nanomolar range, particularly in two triple-negative breast cancer (TNBC) cell lines (from 0.25 to 0.49 µM). In vitro assays revealed that P3C induces reactive oxygen species (ROS) accumulation leading to mitochondrial depolarization and caspase-3/7 and -8 activation, suggesting the participation of both the intrinsic and extrinsic apoptotic pathways. P3C caused microtubule disruption, phosphatidylserine externalization, PARP cleavage, DNA fragmentation, and cell cycle arrest on TNBC cells. In addition, P3C triggered dephosphorylation of CREB, p38, ERK, STAT3, and Fyn, and hyperphosphorylation of JNK and NF-kB in TNBC cells, indicating the inactivation of both p38MAPK/STAT3 and ERK1/2/CREB signaling pathways. In support of our in vitro assays, transcriptome analyses of two distinct TNBC cell lines (MDA-MB-231 and MDA-MB-468 cells) treated with P3C revealed 28 genes similarly affected by the treatment implicated in apoptosis, oxidative stress, protein kinase modulation, and microtubule stability.

G protein-coupled receptors (GPCRs) are the largest class of cell-surface receptor proteins with important functions in signal transduction and often serve as therapeutic drug targets. With the rapidly growing public data on three dimensional (3D) structures of GPCRs and GPCR-ligand interactions, computational prediction of GPCR ligand binding becomes a convincing option to high throughput screening and other experimental approaches during the beginning phases of ligand discovery. In this work, we set out to computationally uncover and understand the binding of a single ligand to GPCRs from several different families. Three-dimensional structural comparisons of the GPCRs that bind to the same ligand revealed local 3D structural similarities and often these regions overlap with locations of binding pockets. These pockets were found to be similar (based on backbone geometry and side-chain orientation using APoc), and they correlate positively with electrostatic properties of the pockets. Moreover, the more similar the pockets, the more likely a ligand binding to the pockets will interact with similar residues, have similar conformations, and produce similar binding affinities across the pockets. These findings can be exploited to improve protein function inference, drug repurposing and drug toxicity prediction, and accelerate the development of new drugs.

G protein-coupled receptors (GPCRs) constitute the largest group of membrane receptor proteins in eukaryotes. Due to their significant roles in various physiological processes such as vision, smell and inflammation, GPCRs are the targets of many prescription drugs. However, the functional and sequence diversity of GPCRs has kept their prediction and classification based on amino acid sequence data as a challenging bioinformatics problem. There are existing computational approaches, mainly using machine learning and statistical methods, to predict and classify GPCRs based on amino acid sequence and sequence derived features. In this paper, we describe a searchable MySQL database, named GPCR-PEnDB (GPCR Prediction Ensemble Database), of confirmed GPCRs and non-GPCRs. It was constructed with the goal of allowing users to conveniently access useful information of GPCRs in a wide range of organisms and to compile reliable training and testing datasets for different combinations of computational tools. This database currently contains 3129 confirmed GPCR and 3575 non-GPCR sequences collected from the UniProtKB/Swiss-Prot protein database, encompassing over 1200 species. The non-GPCR entries include transmembrane proteins for evaluating various prediction programs' abilities to distinguish GPCRs from other transmembrane proteins. Each protein is linked to information about its source organism, classification, sequence lengths and composition, and other derived sequence features. We present examples of using this database along with its graphical user interface, to query for GPCRs with specific sequence properties and to compare the accuracies of five tools for GPCR prediction. This initial version of GPCR-PEnDB will provide a framework for future extensions to include additional sequence and feature data to facilitate the design and assessment of software tools and experimental studies to help understand the functional roles of GPCRs.

Filipino Americans show higher thyroid cancer recurrence rates compared to European Americans. Although they are likely to die of this malignancy, the molecular mechanism has not yet been determined. Recent studies demonstrated that small non-coding RNAs could be utilized to assess thyroid cancer prognosis in tumor models. The goal of this study is to determine whether microRNA (miRNA) signatures are differentially expressed in thyroid cancer in two different ethnic groups. We also determined whether these miRNA signatures are related to cancer staging. This is a retrospective study of archival samples from patients with thyroid cancer (both sexes) in the pathology division from the last ten years at Loma Linda University School of Medicine, California. Deidentified patient demographics were extracted from the patient chart. Discarded formalin-fixed paraffin-embedded tissues were collected post-surgeries. We determined the differential expressions of microRNA in archival samples from Filipino Americans compared to European Americans using the state-of-the-art technique, HiSeq4000. By ingenuity pathway analysis, we determined miRNA targets and the pathways that those targets are involved in. We validated their expressions by real-time quantitative PCR and correlated them with the clinicopathological status in a larger cohort of miRNA samples from both ethnicities. We identified the differentially upregulated/downregulated miRNA clusters in Filipino Americans compared to European Americans. Some of these miRNA clusters are known to target genes that are linked to cancer invasion and metastasis. In univariate analysis, ethnicity and tumor staging were significant factors predicting miR-4633-5p upregulation. When including these factors in a multivariate logistic regression model, ethnicity and tumor staging remained significant independent predictors of miRNA upregulation, whereas the interaction of ethnicity and tumor staging was not significant. In contrast, ethnicity remained an independent predictor of significantly downregulated miR-491-5p and let-7 family. We provide evidence that Filipino Americans showed differentially expressed tumor-tissue-derived microRNA clusters. The functional implications of these miRNAs are under investigation.

G protein-coupled receptors (GPCRs) are the largest class of cell-surface receptor proteins with important functions in signal transduction and often serve as therapeutic drug targets. With the rapidly growing public data on three dimensional (3D) structures of GPCRs and GPCR-ligand interactions, computational prediction of GPCR ligand binding becomes a convincing option to high throughput screening and other experimental approaches during the beginning phases of ligand discovery. Such predictions are cost-efficient and can be important aides for planning wet lab experiments to help elucidate signaling pathways and expedite drug discovery. There are existing computational tools for GPCR ligand binding prediction using various sequence and structural derived features. However, these methods have been typically tested on specific families of GPCRs and none has focused on features that characterize binding of a single ligand to multiple GPCR families. In this work, we set out to uncover and understand the binding of the same ligand to multiple GPCRs of different families. As expected only a few GPCRs share a conserved sequence motif. We observed local 3D structural similarities and local sequence similarities. We also observed that the GPCRs within the same family share similar binding pockets for the ligand. Molecular docking revealed that a ligand can bind to GPCRs of different families with different poses but similar conformations. Finally, the ligands bind to similar pockets with similar electrostatic and solvation properties and share similar residues. These findings can be exploited to improve protein function inference, drug repurposing and drug toxicity prediction, and accelerate the development of new drugs.

We devised a quantitative scoring function to assess the cumulative effects of somatic nonsynonymous single-nucleotide variants (SNVs) on protein-coding genes in patients with ovarian cancer (OvCa) and thyroid cancer (ThCa). The goal is to find novel candidate cancer-related genes for downstream bioinformatics analyses and wet-lab studies. With the Genomic Data Commons as primary data resource, SNV information was extracted from whole-exome sequencing data from patients with these cancers. A cumulative variant scoring function, Q(G), was developed to sum up the deleterious effects of the individual SNVs on gene G. While Q(G) can be computed using any popular functional effect analyzers such as FATHMM-XF, SIFT, PolyPhen, and CADD, we have also established an integrative scoring function iQ(G) that combines the deleterious assessments from different analyzers and demonstrated that iQ(G) is a more effective method for identifying likely cancer-related genes. Based on the iQ(G) rankings, the top three novel genes for OvCa are AHNAK2, UNC13A, and PCDHB4; and those for ThCa are PLEC, HECTD4, and CES1. Furthermore, the top 1% genes with highest iQ(G) scores for each cancer were submitted for KEGG pathway analysis. The results revealed that several genes of the CACNA1 family within the type II diabetes mellitus pathway are likely related to both OvCa and ThCa and suggested other molecular interactions that should be further studied in connection with OvCa prognosis and ThCa treatment.